Search

We consider a sequence $\mathbf{T} = (\mathcal{T}_n : n \in \mathbb{N}^+)$ of trees $\mathcal{T}_n$ where, for some $\Delta \in \mathbb{N}^+$ every $\mathcal{T}_n$ has height at most $\Delta$ and as $n \to \infty$ the minimal number of children of a nonleaf tends to infinity. We can view every tree as a (first-order) $\tau$-structure where $\tau$ is a signature with one binary relation symbol. For a fixed (arbitrary) finite and relational signature $\sigma \supseteq \tau$ we consider the set $\mathbf{W}_n$ of expansions of $\mathcal{T}_n$ to $\sigma$ and a probability distribution $\mathbb{P}_n$ on $\mathbf{W}_n$ which is determined by a (parametrized/lifted) Probabilistic Graphical Model (PGM) $\mathbb{G}$ which can use the information given by $\mathcal{T}_n$. The kind of PGM that we consider uses formulas of a many-valued logic that we call $PLA^*$ with truth values in the unit interval $[0, 1]$. We also use $PLA^*$ to express queries, or events, on $\mathbf{W}_n$. With this setup we prove that, under some assumptions on $\mathbf{T}$, $\mathbb{G}$, and a (possibly quite complex) formula $\varphi(x_1, \ldots, x_k)$ of $PLA^*$, as $n \to \infty$, if $a_1, \ldots, a_k$ are vertices of the tree $\mathcal{T}_n$ then the value of $\varphi(a_1, \ldots, a_k)$ will, with high probability, be almost the same as the value of $\psi(a_1, \ldots, a_k)$, where $\psi(x_1, \ldots, x_k)$ is a ``simple'' formula the value of which can always be computed quickly (without reference to $n$), and $\psi$ itself can be found by using only the information that defines $\mathbf{T}$, $\mathbb{G}$ and $\varphi$. A corollary of this, subject to the same conditions, is a probabilistic convergence law for $PLA^*$-formulas., Comment: arXiv admin note: text overlap with arXiv:2401.04802

We present a pair of seismometers capable of measurement in all six axes of rigid motion. The vacuum-compatible devices implement compact interferometric displacement sensors to surpass the sensitivity of typical electrical readout schemes. Together with the capability to subtract the sensitivity-limiting coupling of ground tilt into horizontal motion, our seismometers can widen the sensing band towards mHz frequencies. This has notable applications across a range of fields requiring access to low-frequency signals, such as seismology and climate research. We particularly highlight their potential application in gravitational-wave observatories (LIGO) for observation of intermediate-mass black holes ($\sim 1000\,M_\odot$). The sensors are based on a near-monolithic fused-silica design consisting of a fused-silica mass and fibre, showing improved stability and robustness to tilt drifts, alignment, and control compared to all-metal or mixed metal-silica designs. We demonstrate tilt sensitivity that surpasses the best commercial alternatives in a significantly reduced footprint compared to our previous iterations of these sensors., Comment: 11 pages, 3 figures

We realise an intrinsic optically pumped magnetic gradiometer based on non-linear magneto-optical rotation. We show that our sensor can reach a gradiometric sensitivity of 18 $\text{fT}/\text{cm}/\sqrt{\text{Hz}}$ and can reject common mode homogeneous magnetic field noise with up to 30 dB attenuation. We demonstrate that our magnetic field gradiometer is sufficiently sensitive and resilient to be employed in biomagnetic applications. In particular, we are able to record the auditory evoked response of the human brain, and to perform real-time magnetocardiography in the presence of external magnetic field disturbances. Our gradiometer provides complementary capabilities in human biomagnetic sensing to optically pumped magnetometers, and opens new avenues in the detection of human biomagnetism., Comment: 7 Pages, 6 Figures

We consider finite relational signatures $\tau \subseteq \sigma$, a sequence of finite base $\tau$-structures $(\mathcal{B}_n : n \in \mathbb{N})$ the cardinalities of which tend to infinity and such that there is a fixed finite upper bound of the degree of (the Gaifman graph of) every $\mathcal{B}_n$. We let $\mathbf{W}_n$ be the set of all expansions to $\sigma$ of $\mathcal{B}_n$ we consider a probabilistic graphical model, a concept used in machine learning and artificial intelligence, to generate a probability distribution $\mathbb{P}_n$ on $\mathbf{W}_n$ for all $n$. We use a many-valued ``probability logic'' with truth values in the unit interval to express probabilities within probabilistic graphical models and to express queries on $\mathbf{W}_n$. This logic uses aggregation functions (e.g. the average) instead of quantifiers and it can express all queries (on finite structures) that can be expressed with first-order logic since the aggregation functions maximum and minimum can be used to express existential and universal quantifications, respectively. The main results concern asymptotic elimination of aggregation functions (the analogue of almost sure elimination of quantifiers for two-valued logics with quantifiers) and the asymptotic distribution of truth values of formulas (the analogue of logical convergence results for two-valued logics).

In this report we present the status of the Mainz group's lattice QCD calculation of the pion transition form factor, which describes the interaction of an on-shell pion with two off-shell photons. This form factor is the main ingredient in the calculation of the pion-pole contribution to hadronic light-by-light scattering in the muon $g-2$. We use the $N_f = 2 + 1$ CLS gauge ensembles, and we update our previous work by including a physical pion mass ensemble (E250). We compute the transition form factor in a moving frame as well as in the pion rest frame in order to have access to a wider range of photon virtualities. In addition to the quark-line connected correlator we also compute the quark-line disconnected diagrams that contribute to the form factor. At the final stage of the analysis, the result on E250 will be combined with the previous work published in 2019 to extrapolate the form factor to the continuum and to physical quark masses., Comment: Proceedings of the 40th International Symposium on Lattice Field Theory (Lattice 2023), July 31st - August 4th, 2023, Fermilab, Batavia, Illinois, USA

We present a comprehensive review of purely leptonic and semileptonic decays of $D^{0(+)}$, $D_s^{+}$, and charmed baryons (including $\Lambda_c^{+}$, $\Xi_c$ and $\Omega_c$). The precise studies of these decays help deepen our understanding and knowledge of quantum chromodynamics via measuring decay constants and form factors, and test the Standard Model through examining the unitarity of Cabibbo-Kobayashi-Maskawa matrix and lepton flavor universality. We give an overview of the theoretical and experimental tools before discussing the recent progress. The data sets collected by BESIII near the production thresholds of $D\bar{D}$, $D_s^{(*)+}D_s^{(*)-}$ and $\Lambda_c^{+}\bar{\Lambda}_c^{-}$ offer important opportunities for studies of charm physics., Comment: 24 pages, 10 figures

We calculate the form factors for the $B \to D^*\ell\nu_\ell$ decay in 2+1 flavor lattice QCD. For all quark flavors, we employ the M\"obius domain-wall action, which preserves chiral symmetry to a good precision. Our gauge ensembles are generated at three lattice cutoffs $a^{-1} \sim 2.5$, 3.6 and 4.5 GeV with pion masses as low as $M_\pi \sim 230$ MeV. The physical lattice size $L$ satisfies the condition $M_\pi L \geq 4$ to control finite volume effects (FVEs), while we simulate a smaller size at the smallest $M_\pi$ to directly examine FVEs. The bottom quark masses are chosen in a range from the physical charm quark mass to $0.7 a^{-1}$ to control discretization effects. We extrapolate the form factors to the continuum limit and physical quark masses based on heavy meson chiral perturbation theory at next-to-leading order. Then the recoil parameter dependence is parametrized using a model independent form leading to our estimate of the decay rate ratio between the tau ($\ell = \tau$) and light lepton ($\ell = e,\mu$) channels $R(D^*) = 0.252(22)$ in the Standard Model. A simultaneous fit with recent data from the Belle experiment yields $|V_{cb}| = 39.19(91)\times 10^{-3}$, which is consistent with previous exclusive determinations, and shows good consistency in the kinematical distribution of the differential decay rate between the lattice and experimental data., Comment: 45 pages, 19 figures; v2: version published in PRD

We consider a logic with truth values in the unit interval and which uses aggregation functions instead of quantifiers, and we describe a general approach to asymptotic elimination of aggregation functions and, indirectly, of asymptotic elimination of Mostowski style generalized quantifiers, since such can be expressed by using aggregation functions. The notion of ``local continuity'' of an aggregation function, which we make precise in two (related) ways, plays a central role in this approach.

A pixelated silicon detector, developed originally for particle physics experiments, was used for a beam profile measurement of a Co-60 irradiator in a water phantom. The beam profile was compared to a profile measured with a pinpoint ionization chamber. The differences in the pixel detector and pinpoint chamber relative profiles were within approximately 0.02, and after calculating correction factors with Monte Carlo simulations for the pixel detector, the differences were decreased to almost less than 0.005. The detector's capability to measure pulse-height was used to record an electron pulse-height spectrum in water in the Co-60 beam, and the results agreed well with simulations.

The isovector axial form factor of the nucleon plays a key role in interpreting data from long-baseline neutrino oscillation experiments. We present a lattice QCD calculation of this form factor, introducing a new method to directly extract its z-expansion from lattice correlators. Our final parameterization of the form factor, which extends up to spacelike virtualities of 0.7 GeV^2 with fully quantified uncertainties, agrees with previous lattice calculations but is significantly less steep than neutrino-deuterium scattering data suggests., Comment: Proceeding for the 39th International Symposium on Lattice Field Theory (Lattice2022), 8th-13th August 2022, Bonn, Germany

The isovector axial form factor of the nucleon plays a key role in interpreting data from long-baseline neutrino oscillation experiments. We perform a lattice-QCD based calculation of this form factor, introducing a new method to directly extract its $z$-expansion from lattice correlators. Our final parametrization of the form factor, which extends up to spacelike virtualities of $0.7\,{\rm GeV}^2$ with fully quantified uncertainties, agrees with previous lattice calculations but is significantly less steep than neutrino-deuterium scattering data suggests., Comment: 6 pages + 6 pages supplementary material; (3+6) figures; 2 tables in suppl.mat

Lattice quantum chromodynamics has proven to be an indispensable method to determine nonperturbative strong contributions to weak decay processes. In this white paper for the Snowmass community planning process we highlight achievements and future avenues of research for lattice calculations of weak $b$ and $c$ quark decays, and point out how these calculations will help to address the anomalies currently in the spotlight of the particle physics community. With future increases in computational resources and algorithmic improvements, percent level (and below) lattice determinations will play a central role in constraining the standard model or identifying new physics., Comment: contribution to Snowmass 2021; 19 pages; v2 corrected typo and added references

The $\beta$ decays of the ground state (gs) and isomeric state (m) of $^{96}$Y have been studied with the total absorption $\gamma$-ray spectroscopy technique at the Ion Guide Isotope Separator On-Line facility. The separation of the 8$^{+}$ isomeric state from the 0$^{-}$ ground state was achieved thanks to the purification capabilities of the JYFLTRAP double Penning trap system. The $\beta$-intensity distributions of both decays have been independently determined. In the analyses the de-excitation of the 1581.6 keV level in $^{96}$Zr, in which conversion electron emission competes with pair production, has been carefully considered and found to have significant impact on the $\beta$-detector efficiency, influencing the $\beta$-intensity distribution obtained. Our results for $^{96\text{gs}}$Y (0$^+$) confirm the large ground state to ground state $\beta$-intensity probability, although a slightly larger value than reported in previous studies was obtained, amounting to $96.6_{-2.1}^{+0.3}\%$ of the total $\beta$ intensity. Given that the decay of $^{96\text{gs}}$Y is the second most important contributor to the reactor antineutrino spectrum between 5 and 7 MeV, the impact of the present results on reactor antineutrino summation calculations has been evaluated. In the decay of $^{96\text{m}}$Y (8$^{+}$), previously undetected $\beta$ intensity in transitions to states above 6 MeV has been observed. This shows the importance of total absorption $\gamma$-ray spectroscopy measurements of $\beta$ decays with highly fragmented de-excitation patterns. $^{96\text{m}}$Y (8$^{+}$) is a major contributor to reactor decay heat in uranium-plutonium and thorium-uranium fuels around 10 s after fission pulses, and the newly measured average $\beta$ and $\gamma$ energies differ significantly from the previous values in evaluated databases (...)

We consider logics with truth values in the unit interval $[0,1]$. Such logics are used to define queries and to define probability distributions. In this context the notion of almost sure equivalence of formulas is generalized to the notion of asymptotic equivalence. We prove two new results about the asymptotic equivalence of formulas where each result has a convergence law as a corollary. These results as well as several older results can be formulated as results about the relative asymptotic expressivity of inference frameworks. An inference framework $\mathbf{F}$ is a class of pairs $(\mathbb{P}, L)$, where $\mathbb{P} = (\mathbb{P}_n : n = 1, 2, 3, \ldots)$, $\mathbb{P}_n$ are probability distributions on the set $\mathbf{W}_n$ of all $\sigma$-structures with domain $\{1, \ldots, n\}$ (where $\sigma$ is a first-order signature) and $L$ is a logic with truth values in the unit interval $[0, 1]$. An inference framework $\mathbf{F}'$ is asymptotically at least as expressive as an inference framework $\mathbf{F}$ if for every $(\mathbb{P}, L) \in \mathbf{F}$ there is $(\mathbb{P}', L') \in \mathbf{F}'$ such that $\mathbb{P}$ is asymptotically total variation equivalent to $\mathbb{P}'$ and for every $\varphi(\bar{x}) \in L$ there is $\varphi'(\bar{x}) \in L'$ such that $\varphi'(\bar{x})$ is asymptotically equivalent to $\varphi(\bar{x})$ with respect to $\mathbb{P}$. This relation is a preorder. If, in addition, $\mathbf{F}$ is at least as expressive as $\mathbf{F}'$ then we say that $\mathbf{F}$ and $\mathbf{F}'$ are asymptotically equally expressive. Our third contribution is to systematize the new results of this paper and several previous results in order to get a preorder on a number of inference systems that are of relevance in the context of machine learning and artificial intelligence.

The precision of lattice QCD calculations has been steadily improving for some time and is now approaching, or has surpassed, the 1% level for multiple quantities. At this level QED effects, i.e. the fact that quarks carry electric as well as color charge, come into play. In this report we will summarise results from the first lattice QCD+QED computations of the properties of ground-state charmonium and bottomonium mesons by the HPQCD Collaboration., Comment: Contribution to the proceedings of the 19th International Conference on Hadron Spectroscopy and Structure (HADRON 2021). Accepted for publication in Suplemento de la Revista Mexicana de F\'isica

Using a fully relativistic lattice fermion action, we compute the form factors of the semileptonic decay $B\to\pi\ell\nu$, which is required for the determination of the Cabibbo-Kobayashi-Maskawa matrix element $|V_{ub}|$. We employ the M\"{o}bius domain-wall fermion formalism for the generation of lattice ensembles with 2+1 sea quark flavours as well as for the valence heavy and light quarks. We compute the form factors at various values of the lattice spacing and multiple light and heavy quark masses, and extrapolate the results to the physical point. We combine our lattice results with the available experimental data to obtain $|V_{ub}| = (3.93\pm 0.41)\times 10^{-3}$., Comment: 34 pages, 20 figures, 12 tables. Plus supplemental material: 20 pages, 7 figures, 8 tables. Updated to reflect minor changes in published version, including additional references

Cadmium Telluride (CdTe) is a high-Z material with excellent photon radiation absorption properties, making it a promising material to include in radiation detection technologies. However, the brittleness of CdTe crystals as well as their varying concentration of defects necessitate a thorough quality assessment before the complex detector processing procedure. We present our quality assessment of CdTe as a detector material for multispectral medical imaging, a research which is conducted as part of the Consortium Project Multispectral Photon-counting for Medical Imaging and Beam characterization (MPMIB). The aim of the project is to develop novel CdTe detectors and obtain spectrum-per-pixel information that make the distinction between different radiation types and tissues possible. To evaluate the defect density inside the crystals -- which can deteriorate the detector performance -- we employ infrared microscopy (IRM). Posterior data analysis allows us to visualise the defect distributions as 3D defect maps. Additionally, we investigate front and backside differences of the material with current-voltage (IV) measurements to determine the preferred surface for the pixelisation of the crystal, and perform test measurements with the prototypes to provide feedback for further processing. We present the different parts of our quality assessment chain and will close with first experimental results obtained with one of our prototype PC detectors in a small tomographic setup., Comment: This is the Accepted Manuscript version of an article accepted for publication in Journal of Instrumentation. IOP Publishing Ltd is not responsible for any errors or omissions in this version of the manuscript or any version derived from it. The Version of Record is available online at https://doi.org/10.1088/1748-0221/17/01/C01070

We calculate the $B \to D^{(*)} \ell \nu$ form factors in 2+1 flavor relativistic lattice QCD by employing the M\"obius domain-wall action for all quark flavors. Our simulations are carried out at lattice cut-offs $a^{-1} \sim 2.5$, 3.6 and 4.5 GeV with the bottom quark masses up to 0.7 $a^{-1}$ to control discretization effects. We extrapolate the form factors to the continuum limit and physical quark masses, and discuss systematic uncertainties of the form factors., Comment: 9 pages, 4 figures, talk presented at the 38th International Symposium on Lattice Field Theory (LATTICE2021), 26th-30th July, 2021, Zoom/Gather@MIT, USA

We present the analysis of isovector axial vector nucleon form factors on a set of $N_f=2+1$ CLS ensembles with $\mathcal O(a)$-improved Wilson fermions and L\"uscher-Weisz gauge action. The set of ensembles covers a pion mass range of $130-353\,$MeV with lattice spacings between $0.05\,$fm and $0.09\,$fm. In particular, the set includes a $L/a=96$ ensemble at the physical pion mass. For the purpose of the form factor extraction, we employ both the summed operator insertion method (summation method) and explicit two-state fits in order to account for excited-state contributions to the nucleon correlation functions. To describe the $Q^{2}$-behavior of the form factors, we perform $z$-expansion fits. Finally, we present HBChPT-inspired chiral and continuum extrapolations of the axial charge and radius., Comment: 10 pages, 3 figures, 1 table. Talk presented at the 38th International Symposium on Lattice Field Theory (Lattice 2021), 26-30 July, 2021, Zoom/Gather@Massachusetts Institute of Technology

As tests of QCD in the bottomonium system, we give the most accurate results to date for the ground-state hyperfine splitting and the $\Upsilon$ leptonic width from full lattice QCD. These quantities are both accurately known from experiment, so can provide a good test of $b$ physics, but previous lattice results have been rather imprecise. We also test the impact on these quantities of the $b$ quark's electric charge. Our results are: $M_{\Upsilon}-M_{\eta_b} = $ 57.5(2.3)(1.0) MeV (where the second uncertainty comes from neglect of quark-line disconnected correlation functions) and decay constants, $f_{\eta_b} =$ 724(12) MeV and $f_{\Upsilon} =$ 677.2(9.7) MeV, giving $\Gamma(\Upsilon \rightarrow e^+e^-) =$ 1.292(37)(3) keV. We also give a new determination of the ratio of the masses for $b$ and $c$ quarks that is completely nonperturbative in lattice QCD and includes the calculation of QED effects for the first time. This gives a result for the $b$ quark mass of $\overline{m}_b(\overline{m}_b,n_f=5) =$ 4.202(21) GeV., Comment: 6 pages, The 38th International Symposium on Lattice Field Theory, LATTICE2021 26th-30th July, 2021 Zoom/Gather@Massachusetts Institute of Technology

Semileptonic $D \to{}K \ell \nu$ decays provide one angle of attack to get at the CKM matrix element $V_{cs}$, complementary to the study of leptonic $D_s$ decays. Here, HPQCD present the results of a recently published, improved determination of $V_{cs}$. We discuss a new, precise determination of $D\to K$ scalar and vector form factors from a lattice calculation on eight different $N_f=2+1+1$ MILC gluon field ensembles using the HISQ action, including three with physical light quark masses. When combined with experimental results, we are able to extract $|V_{cs}|=0.9663(80)$ to a sub percent level of precision for the first time. This is achieved using three different methods, which each combine our form factors with different sets of experimental results in different ways, with the results in very good agreement. Our primary method is to use $q^2$-binned data for the differential decay rate, but we also calculate $V_{cs}$ from the total branching fraction and from the value $|V_{cs}|f_+(0)$, which is also quoted by some experiments., Comment: 9 pages, 7 figures. Proceedings of the 38th International Symposium on Lattice Field Theory - LATTICE2021

In Statistical Relational Artificial Intelligence, a branch of AI and machine learning which combines the logical and statistical schools of AI, one uses the concept {\em para\-metrized probabilistic graphical model (PPGM)} to model (conditional) dependencies between random variables and to make probabilistic inferences about events on a space of "possible worlds". The set of possible worlds with underlying domain $D$ (a set of objects) can be represented by the set $\mathbf{W}_D$ of all first-order structures (for a suitable signature) with domain $D$. Using a formal logic we can describe events on $\mathbf{W}_D$. By combining a logic and a PPGM we can also define a probability distribution $\mathbb{P}_D$ on $\mathbf{W}_D$ and use it to compute the probability of an event. We consider a logic, denoted $PLA$, with truth values in the unit interval, which uses aggregation functions, such as arithmetic mean, geometric mean, maximum and minimum instead of quantifiers. However we face the problem of computational efficiency and this problem is an obstacle to the wider use of methods from Statistical Relational AI in practical applications. We address this problem by proving that the described probability will, under certain assumptions on the PPGM and the sentence $\varphi$, converge as the size of $D$ tends to infinity. The convergence result is obtained by showing that every formula $\varphi(x_1, \ldots, x_k)$ which contains only "admissible" aggregation functions (e.g. arithmetic and geometric mean, max and min) is asymptotically equivalent to a formula $\psi(x_1, \ldots, x_k)$ without aggregation functions.

In ultrasound tomography, the speed of sound inside an object is estimated based on acoustic measurements carried out by sensors surrounding the object. An accurate forward model is a prominent factor for high-quality image reconstruction, but it can make computations far too time-consuming in many applications. Using approximate forward models, it is possible to speed up the computations, but the quality of the reconstruction may have to be compromised. In this paper, a neural network -based approach is proposed, that can compensate for modeling errors caused by the approximate forward models. The approach is tested with various different imaging scenarios in a simulated two-dimensional domain. The results show that with fairly small training datasets, the proposed approach can be utilized to approximate the modelling errors, and to significantly improve the image reconstruction quality in ultrasound tomography, compared to commonly used inversion algorithms.

We provide a 0.8\%-accurate determination of $V_{cs}$ from combining experimental results for the differential rate of $D \rightarrow K$ semileptonic decays with precise form factors that we determine from lattice QCD. This is the first time that $V_{cs}$ has been determined with an accuracy that allows its difference from 1 to be seen. Our lattice QCD calculation uses the Highly Improved Staggered Quark (HISQ) action for all valence quarks on gluon field configurations generated by the MILC collaboration that include the effect of $u$, $d$, $s$ and $c$ HISQ quarks in the sea. We use eight gluon field ensembles with five values of the lattice spacing ranging from 0.15 fm to 0.045 fm and include results with physical $u/d$ quarks for the first time. Our calculated form factors cover the full $q^2$ range of the physical decay process and enable a Standard Model test of the shape of the differential decay rate as well as the determination of $V_{cs}$ from a correlated weighted average over $q^2$ bins. We obtain $|V_{cs}|= 0.9663(53)_{\text{latt}}(39)_{\text{exp}}(19)_{\eta_{EW}}(40)_{\text{EM}}$, where the uncertainties come from lattice QCD, experiment, short-distance electroweak and electromagnetic corrections, respectively. This last uncertainty, neglected for $D \rightarrow K \ell \nu$ hitherto, now needs attention if the uncertainty on $V_{cs}$ is to be reduced further. We also determine $V_{cs}$ values in good agreement using the measured total branching fraction and the rates extrapolated to $q^2=0$. Our form factors enable tests of lepton flavour universality violation. We find the ratio of branching fractions for $D^0 \rightarrow K^-$ with $\mu$ and $e$ in the final state to be $R_{\mu/e}=0.9779(2)_{\text{latt}}(50)_{\mathrm{EM}}$ in the Standard Model, with the uncertainty dominated by that from electromagnetic corrections., Comment: 27 pages, 28 figures. Small changes to text to improve readability. Updated to include the BES average for the D to K branching fraction from 2104.08081

We extend HPQCD's earlier $n_f=4$ lattice-QCD analysis of the ratio of $\overline{\mathrm{MSB}}$ masses of the $b$ and $c$ quark to include results from finer lattices (down to 0.03fm) and a new calculation of QED contributions to the mass ratio. We find that $\overline{m}_b(\mu)/\overline{m}_c(\mu)=4.586(12)$ at renormalization scale $\mu=3$\,GeV. This result is nonperturbative. Combining it with HPQCD's recent lattice QCD$+$QED determination of $\overline{m}_c(3\mathrm{GeV})$ gives a new value for the $b$-quark mass: $\overline{m}_b(3\mathrm{GeV}) = 4.513(26)$GeV. The $b$-mass corresponds to $\overline{m}_b(\overline{m}_b, n_f=5) = 4.202(21)$GeV. These results are the first based on simulations that include QED., Comment: 8 pages, 4 figures

We calculate the mass difference between the $\Upsilon$ and $\eta_b$ and the $\Upsilon$ leptonic width from lattice QCD using the Highly Improved Staggered Quark formalism for the $b$ quark and including $u$, $d$, $s$ and $c$ quarks in the sea. We have results for lattices with lattice spacing as low as 0.03 fm and multiple heavy quark masses, enabling us to map out the heavy quark mass dependence and determine values at the $b$ quark mass. Our results are: $M_{\Upsilon} -M_{\eta_b} = 57.5(2.3)(1.0) \,\mathrm{MeV}$ (where the second uncertainty comes from neglect of quark-line disconnected correlation functions) and decay constants, $f_{\eta_b}=724(12)$ MeV and $f_{\Upsilon} =677.2(9.7)$ MeV, giving $\Gamma(\Upsilon \rightarrow e^+e^-) = 1.292(37)(3) \,\mathrm{keV}$. The hyperfine splitting and leptonic width are both in good agreement with experiment, and provide the most accurate lattice QCD results to date for these quantities by some margin. At the same time results for the time moments of the vector-vector correlation function can be compared to values for the $b$ quark contribution to $\sigma(e^+e^- \rightarrow \mathrm{hadrons})$ determined from experiment. Moments 4--10 provide a 2\% test of QCD and yield a $b$ quark contribution to the anomalous magnetic moment of the muon of 0.300(15)$\times 10^{-10}$. Our results, covering a range of heavy quark masses, may also be useful to constrain QCD-like composite theories for beyond the Standard Model physics., Comment: 16 pages, 9 figures Version accepted by PRD: Appendix with tables of fit parameters that allow the dependence on heavy quark mass of the hyperfine splitting and decay constants to be reconstructed added

Heavy-tailed networks, which have degree distributions characterised by slower than exponentially bounded tails, are common in many different situations. Some interesting cases, where heavy tails are characterised by inverse powers $\lambda$ in the range $1<\lambda<2,$ arise for associative knowledge networks, and semantic and linguistic networks. In these cases, the differences between the networks are often delicate, calling for robust methods to characterise the differences. Here, we introduce a method for comparing networks using a density matrix based on q-generalised adjacency matrix kernels. It is shown that comparison of networks can then be performed using the q-generalised Kullback-Leibler divergence. In addition, the q-generalised divergence can be interpreted as a q-generalised free energy, which enables the thermodynamic-like macroscopic description of the heavy-tailed networks. The viability of the q-generalised adjacency kernels and the thermodynamic-like description in characterisation of complex networks is demonstrated using a simulated set of networks, which are modular and heavy-tailed with a degree distribution of inverse power law in the range $1<\lambda<2$.

The coloured dissolved organic matter (CDOM) concentration is the standard measure of humic substance in natural waters. CDOM measurements by remote sensing is calculated using the absorption coefficient (a) at a certain wavelength (e.g. 440nm). This paper presents a comparison of four machine learning methods for the retrieval of CDOM from remote sensing signals: regularized linear regression (RLR), random forest (RF), kernel ridge regression (KRR) and Gaussian process regression (GPR). Results are compared with the established polynomial regression algorithms. RLR is revealed as the simplest and most efficient method, followed closely by its nonlinear counterpart KRR., Comment: 7 pages, 4 figures

Game recommendation is an important application of recommender systems. Recommendations are made possible by data sets of historical player and game interactions, and sometimes the data sets include features that describe games or players. Collaborative filtering has been found to be the most accurate predictor of past interactions. However, it can only be applied to predict new interactions for those games and players where a significant number of past interactions are present. In other words, predictions for completely new games and players is not possible. In this paper, we use a survey data set of game likes to present content based interaction models that generalize into new games, new players, and both new games and players simultaneously. We find that the models outperform collaborative filtering in these tasks, which makes them useful for real world game recommendation. The content models also provide interpretations of why certain games are liked by certain players for game analytics purposes.

Objective: This work aims to reduce the acoustic noise level of transcranial magnetic stimulation (TMS) coils. TMS requires high currents (several thousand amperes) to be pulsed through the coil, which generates a loud acoustic impulse whose peak sound pressure level (SPL) can exceed 130 dB(Z). This sound poses a risk to hearing and elicits unwanted neural activation of auditory brain circuits. Methods: We propose a new double-containment coil with enhanced winding mounting (DCC), which utilizes acoustic impedance mismatch to contain and dissipate the impulsive sound within an air-tight outer casing. The coil winding is potted in a rigid block, which is mounted to the outer casing by its acoustic nodes that are subject to minimum vibration during the pulse. The rest of the winding block is isolated from the casing by an air gap, and sound is absorbed by foam within the casing. The casing thickness under the winding center is minimized to maximize the coil electric field output. Results: Compared to commercial figure-of-eight TMS coils, the DCC prototype has 10-33 dB(Z) lower SPL at matched stimulation strength, whilst providing 22% higher maximum stimulation strength than equally focal commercial coils. Conclusion: The DCC design greatly reduces the acoustic noise of TMS while increasing the achievable stimulation strength. Significance: The acoustic noise reduction from our coil design is comparable to that provided by typical hearing protection devices. This coil design approach can enhance hearing safety and reduce auditory co-activations in the brain and other detrimental effects of TMS sound., Comment: 8 pages, 5 figures

We review the present status of the Standard Model calculation of the anomalous magnetic moment of the muon. This is performed in a perturbative expansion in the fine-structure constant $\alpha$ and is broken down into pure QED, electroweak, and hadronic contributions. The pure QED contribution is by far the largest and has been evaluated up to and including $\mathcal{O}(\alpha^5)$ with negligible numerical uncertainty. The electroweak contribution is suppressed by $(m_\mu/M_W)^2$ and only shows up at the level of the seventh significant digit. It has been evaluated up to two loops and is known to better than one percent. Hadronic contributions are the most difficult to calculate and are responsible for almost all of the theoretical uncertainty. The leading hadronic contribution appears at $\mathcal{O}(\alpha^2)$ and is due to hadronic vacuum polarization, whereas at $\mathcal{O}(\alpha^3)$ the hadronic light-by-light scattering contribution appears. Given the low characteristic scale of this observable, these contributions have to be calculated with nonperturbative methods, in particular, dispersion relations and the lattice approach to QCD. The largest part of this review is dedicated to a detailed account of recent efforts to improve the calculation of these two contributions with either a data-driven, dispersive approach, or a first-principle, lattice-QCD approach. The final result reads $a_\mu^\text{SM}=116\,591\,810(43)\times 10^{-11}$ and is smaller than the Brookhaven measurement by 3.7$\sigma$. The experimental uncertainty will soon be reduced by up to a factor four by the new experiment currently running at Fermilab, and also by the future J-PARC experiment. This and the prospects to further reduce the theoretical uncertainty in the near future-which are also discussed here-make this quantity one of the most promising places to look for evidence of new physics., Comment: 196 pages, 103 figures, version published in Phys. Rept., bib files for the citation references are available from: https://muon-gm2-theory.illinois.edu

In $\beta$-decay studies the determination of the decay probability to the ground state of the daughter nucleus often suffers from large systematic errors. The difficulty of the measurement is related to the absence of associated delayed $\gamma$-ray emission. In this work we revisit the $4\pi\gamma-\beta$ method proposed by Greenwood and collaborators in the 1990s, which has the potential to overcome some of the experimental difficulties. Our interest is driven by the need to determine accurately the $\beta$-intensity distributions of fission products that contribute significantly to the reactor decay heat and to the antineutrinos emitted by reactors. A number of such decays have large ground state branches. The method is relevant for nuclear structure studies as well. Pertinent formulae are revised and extended to the special case of $\beta$-delayed neutron emitters, and the robustness of the method is demonstrated with synthetic data. We apply it to a number of measured decays that serve as test cases and discuss the features of the method. Finally, we obtain ground state feeding intensities with reduced uncertainty for four relevant decays that will allow future improvements in antineutrino spectrum and decay heat calculations using the summation method.

We have performed the first $n_f = 2+1+1$ lattice QCD computations of the properties (masses and decay constants) of ground-state charmonium mesons. Our calculation uses the HISQ action to generate quark-line connected two-point correlation functions on MILC gluon field configurations that include $u/d$ quark masses going down to the physical point, tuning the $c$ quark mass from $M_{J/\psi}$ and including the effect of the $c$ quark's electric charge through quenched QED. We obtain $M_{J/\psi}-M_{\eta_c}$ (connected) = 120.3(1.1) MeV and interpret the difference with experiment as the impact on $M_{\eta_c}$ of its decay to gluons, missing from the lattice calculation. This allows us to determine $\Delta M_{\eta_c}^{\mathrm{annihiln}}$ =+7.3(1.2) MeV, giving its value for the first time. Our result of $f_{J/\psi}=$ 0.4104(17) GeV, gives $\Gamma(J/\psi \rightarrow e^+e^-)$=5.637(49) keV, in agreement with, but now more accurate than experiment. At the same time we have improved the determination of the $c$ quark mass, including the impact of quenched QED to give $\overline{m}_c(3\,\mathrm{GeV})$ = 0.9841(51) GeV. We have also used the time-moments of the vector charmonium current-current correlators to improve the lattice QCD result for the $c$ quark HVP contribution to the anomalous magnetic moment of the muon. We obtain $a_{\mu}^c = 14.638(47) \times 10^{-10}$, which is 2.5$\sigma$ higher than the value derived using moments extracted from some sets of experimental data on $R(e^+e^- \rightarrow \mathrm{hadrons})$. This value for $a_{\mu}^c$ includes our determination of the effect of QED on this quantity, $\delta a_{\mu}^c = 0.0313(28) \times 10^{-10}$., Comment: Added extra discussion on QED setup, some new results to study the effects of strong isospin breaking in the sea (including new Fig. 1) and a fit stability plot for the hyperfine splitting (new Fig. 7). Version accepted for publication in PRD

We introduce a formal logical language, called conditional probability logic (CPL), which extends first-order logic and which can express probabilities, conditional probabilities and which can compare conditional probabilities. Intuitively speaking, although formal details are different, CPL can express the same kind of statements as some languages which have been considered in the artificial intelligence community. We also consider a way of making precise the notion of lifted Bayesian network, where this notion is a type of (lifted) probabilistic graphical model used in machine learning, data mining and artificial intelligence. A lifted Bayesian network (in the sense defined here) determines, in a natural way, a probability distribution on the set of all structures (in the sense of first-order logic) with a common finite domain $D$. Our main result is that for every "noncritical" CPL-formula $\varphi(\bar{x})$ there is a quantifier-free formula $\varphi^*(\bar{x})$ which is "almost surely" equivalent to $\varphi(\bar{x})$ as the cardinality of $D$ tends towards infinity. This is relevant for the problem of making probabilistic inferences on large domains $D$, because (a) the problem of evaluating, by "brute force", the probability of $\varphi(\bar{x})$ being true for some sequence $\bar{d}$ of elements from $D$ has, in general, (highly) exponential time complexity in the cardinality of $D$, and (b) the corresponding probability for the quantifier-free $\varphi^*(\bar{x})$ depends only on the lifted Bayesian network and not on $D$. The main result has two corollaries, one of which is a convergence law (and zero-one law) for noncritial CPL-formulas., Comment: This version of the article corrects an embarrassing typo in the journal version (in Theoretical Computer Science) of the article, as well as some other typos. The main correction is that '1 - e^{-cn}' in Theorem 3.15 in the journal version has been replaced with 'e^{-cn}' (which is the intended expression) in this version

We analyze the inclusive semileptonic $B \to X_c \ell\nu$ structure functions in 2+1-flavor lattice QCD. The M\"obius domain-wall fermion action is used for light, strange, charm and bottom quarks. The structure function receives contributions from various exclusive modes, including the dominant S-wave states $D^{(*)}_s$ as well as the P-wave states $D_s^{**}$. We can identify them in the lattice data, from which we put some constraints on the $B_s \to D_s^{**}\ell\nu$ form factors., Comment: 7 pages, 3 figures, talk presented at the 37th International Symposium on Lattice Field Theory (Lattice2019), 16-22 June 2019, Wuhan, China

We report on our calculation of the B \to D^(*) \ell \nu form factors in 2+1 flavor lattice QCD. The M\"obius domain-wall action is employed for light, strange, charm and bottom quarks. At lattice cutoffs 1/a \sim 2.4, 3.6 and 4.5 GeV, we simulate bottom quark masses up to 0.7/a to control discretization errors. The pion mass is as low as 230 MeV. We extrapolate the form factors to the continuum limit and physical quark masses, and make a comparison with recent phenomenological analyses., Comment: 7 pages, 4 figures, talk presented at the 37th International Symposium on Lattice Field Theory (Lattice2019), 16-22 June 2019, Wuhan, China

We report on a calculation of form factors for the semileptonic decay of $B$ meson to pion on $2+1$-flavour lattices with lattice spacings from 0.080 fm down to 0.044 fm. Using the M\"obius domain wall fermion action for both sea and valence quarks, we simulate pions with masses down to 225 MeV. By utilizing a range of heavy quark masses up to 2.44 times the mass of the charm quark we extrapolate to the physical $b$ quark mass. We discuss the dependence of the form factors on the pion mass, heavy quark mass, lattice spacing and the momentum-transfer. We extract the CKM matrix element $|V_{ub}|$ through a simultaneous fit with the $B \to \pi\ell\nu$ differential branching fractions provided by the Belle and BaBar collaborations after a chiral-continuum and physical $b$ quark extrapolations of our lattice data., Comment: 7 pages, Presented at the 37th International Symposium on Lattice Field Theory - Lattice2019, 16-22 June 2019, Wuhan, China

We present a lattice QCD determination of light quark masses with three sea-quark flavours ($N_f = 2+1$). Bare quark masses are known from PCAC relations in the framework of CLS lattice computations with a non-perturbatively improved Wilson-Clover action and a tree-level Symanzik improved gauge action. They are fully non-perturbatively improved, including the recently computed Symanzik counter-term $b_{\rm A} - b_{\rm P}$. The mass renormalisation at hadronic scales and the renormalisation group running over a wide range of scales are known non-perturbatively in the Schr\"odinger functional scheme. In the present paper we perform detailed extrapolations to the physical point, obtaining (for the four-flavour theory) $m_{u/d}(2{\rm GeV}) = 3.54(12)(9)$ MeV and $m_s(2{\rm GeV}) = 95.7(2.5)(2.4)$ MeV in the $\bar{MS}$ scheme. For the mass ratio we have $m_s/m_{u/d} = 27.0(1.0)(0.4)$. The RGI values in the three-flavour theory are $M_{u/d} = 4.70(15)(12)$ MeV and $M_s = 127.0(3.1)(3.2)$ MeV., Comment: 27 pages, 5 figures. Published version of the paper. The present version contains additions and improvements proposed by the referee; the results remain unchanged

The decays of the $\beta$-delayed neutron emitters $^{137}$I and $^{95}$Rb have been studied with the total absorption $\gamma$-ray spectroscopy technique. The purity of the beams provided by the JYFLTRAP Penning trap at the ion guide isotope separator on-line facility in Jyv\"askyl\"a allowed us to carry out a campaign of isotopically pure measurements with the decay total absorption $\gamma$-ray spectrometer, a segmented detector composed of eighteen NaI(Tl) modules. The contamination coming from the interaction of neutrons with the spectrometer has been carefully studied, and we have tested the use of time differences between prompt $\gamma$-rays and delayed neutron interactions to eliminate this source of contamination. Due to the sensitivity of our spectrometer, we have found a significant amount of $\beta$-intensity to states above the neutron separation energy that de-excite by $\gamma$-rays, comparable to the neutron emission probability. The competition between $\gamma$ de-excitation and neutron emission has been compared with Hauser-Feshbach calculations, and it can be understood as a nuclear structure effect. In addition, we have studied the impact of the $\beta$-intensity distributions determined in this work on reactor decay heat and reactor antineutrino spectrum summation calculations. The robustness of our results is demonstrated by a thorough study of uncertainties, and with the reproduction of the spectra of the individual modules and the module-multiplicity gated spectra. This work represents the state-of-the-art of our analysis methodology for segmented total absorption spectrometers.

We present a lattice QCD determination of the $B_s \to D_s \ell\nu$ scalar and vector form factors over the full physical range of momentum transfer. The result is derived from correlation functions computed using the Highly Improved Staggered Quark (HISQ) formalism, on the second generation MILC gluon ensembles accounting for up, down, strange and charm contributions from the sea. We calculate correlation functions for three lattice spacing values and an array of unphysically light $b$-quark masses, and extrapolate to the physical value. Using the HISQ formalism for all quarks means that the lattice current coupling to the $W$ can be renormalized non-perturbatively, giving a result free from perturbative matching errors for the first time. Our results are in agreement with, and more accurate than, previous determinations of these form factors. From the form factors we also determine the ratio of branching fractions that is sensitive to violation of lepton universality: $R(D_s) = \mathcal{B}(B_s\to D_s \tau \nu_{\tau})/\mathcal{B}(B_s\to D_s \ell \nu_{l})$, where $\ell$ is an electron or a muon. We find $R(D_s) = 0.2987(46)$, which is also more accurate than previous lattice QCD results. Combined with a future measurement of $R(D_s)$, this could supply a new test of the Standard Model. We also compare the dependence on heavy quark mass of our form factors to expectations from Heavy Quark Effective Theory., Comment: 22 pages, 17 figures

The $\beta$ intensity distributions of the decays of $^{100\text{gs},100\text{m}}$Nb and $^{102\text{gs},102\text{m}}$Nb have been determined using the Total Absorption $\gamma$-Ray Spectroscopy technique. The JYFLTRAP double Penning trap system was employed to disentangle the isomeric states involved, lying very close in energy, in a campaign of challenging measurements performed with the Decay Total Absorption $\gamma$-ray Spectrometer at the Ion Guide Isotope Separator On-Line facility in Jyv\"askyl\"a. The low-spin isomeric state of each niobium case was populated through the decay of the zirconium parent, that was treated as a contaminant. We have applied a method to extract this contamination, and additionally we have obtained $\beta$ intensity distributions for these zirconium decays. The $\beta$-strength distributions evaluated with these results were compared with calculations in quasiparticle random-phase approximation, suggesting a prolate configuration for the ground states of $^{100,102}$Zr. The footprint of the Pandemonium effect was found when comparing our results for the analyses of the niobium isotopes with previous decay data. The $\beta$-intensities of the decay of $^{102\text{m}}$Nb were obtained for the first time. A careful evaluation of the uncertainties was carried out, and the consistency of our results was validated taking advantage of the segmentation of our spectrometer. The final results were used as input in reactor summation calculations. A large impact on antineutrino spectrum calculations was already reported and here we detail the significant impact on decay heat calculations.

We present details of a lattice QCD calculation of the $B_s\to D_s^*$ axial form factor at zero recoil using the Highly Improved Staggered Quark (HISQ) formalism on the second generation MILC gluon ensembles that include up, down, strange and charm quarks in the sea. Using the HISQ action for all valence quarks means that the lattice axial vector current that couples to the $W$ can be renormalized fully non-perturbatively, giving a result free of the perturbative matching errors that previous lattice QCD calculations have had. We calculate correlation functions at three values of the lattice spacing, and multiple `$b$'-quark masses, for physical $c$ and $s$. The functional dependence on the $b$-quark mass can be determined and compared to Heavy Quark Effective Theory expectations, and a result for the form factor obtained at the physical value of the $b$-quark mass. We find $\mathcal{F}^{B_s\to D_s^*}(1) = h^s_{A_1}(1) = 0.9020(96)_{\text{stat}}(90)_{\text{sys}}$. This is in agreement with earlier lattice QCD results, which use NRQCD $b$ quarks, with a total uncertainty reduced by more than a factor of two. We discuss implications of this result for the $B\to D^*$ axial form factor at zero recoil and for determinations of $V_{cb}$., Comment: 18 pages, 9 figures

We present a nearly final analysis of the $u/d$ and $s$ quark masses, extracted using the PCAC quark masses reported in [PRD 95 (2017) 074504]. The data is based on the CLS $N_f = 2 + 1$ simulations with Wilson/Clover quarks and L\"uscher-Weisz gauge action, at four $\beta$ values (i.e. lattice spacings) and a range of quark masses. We use the ALPHA results of [EPJC 78 (2018) 387] for non-perturbative quark mass renormalisation and RG-running from hadronic to electroweak scales in the Schr\"odinger Functional scheme. Quark masses are quoted both in the $\overline{\rm MS}$ scheme and as RGI quantities., Comment: Talk presented at the 36th Annual International Symposium on Lattice Field Theory (Lattice 2018), 22-28 July 2018, Michigan, USA

Lattice QCD can provide a direct determination of meson electromagnetic form factors, making predictions for upcoming experiments at Jefferson Lab. The form factors are a reflection of the bound-state nature of the meson and so these calculations give information about how confinement by QCD affects meson internal structure. The region of high squared (space-like) momentum-transfer, $Q^2$, is of particular interest because perturbative QCD predictions take a simple form in that limit that depends on the meson decay constant. We previously showed in\cite{jonnaff} that, up to $Q^2$ of 6 $\mathrm{GeV}^2$, the form factor for a `pseudo-pion' made of strange quarks was significantly larger than the asymptotic perturbative QCD result and showed no sign of heading towards that value at higher $Q^2$. Here we give predictions for real mesons, the $K^+$ and $K^0$, in anticipation of JLAB results for the $K^+$ in the next few years. We also give results for a heavier meson, the $\eta_c$, up to $Q^2$ of 25 $\mathrm{GeV}^2$ for a comparison to perturbative QCD in a higher $Q^2$ regime., Comment: 6 pages, 5 figures, Proceedings of LATTICE2018

We present progress on an ongoing calculation of the $B_s\to D_s^{(*)} l \nu$ form factors calculated on the $n_f=2+1+1$ MILC ensembles and using the Highly Improved Staggered Quark action for all valence quarks. We perform the calculation at a range of $b$ quark masses (and lattice spacings) so that we can extrapolate to the physical $b$-quark mass., Comment: 6 pages, 3 figures, proceedings of the 36th Annual International Symposium on Lattice Field Theory - LATTICE2018

We report on a calculation of $B_c$ ground state and radial excitation energies, obtained from heavy-charm highly improved staggered quark (HISQ) correlators computed on MILC gauge ensembles, with lattice spacings down to $a=0.044$ fm. Using HISQ valence quarks on progressively finer lattices allows us to simulate up to the $b$-quark mass. In particular we focus on the $B_c(2S)$ energy, which we compare with O(\alpha_s)-improved non-relativistic QCD results computed on the same ensembles and recent experimental results from ATLAS., Comment: 7 pages, 3 figures. Presented at the 36th Annual International Symposium on Lattice Field Theory, July 22-28, 2018, East Lansing, MI, USA

We derive the Operator Product Expansion whose vacuum expectation value gives the time-moments of the pseudoscalar heavy-light current-current correlator up to and including terms in $\alpha_s^2$ multiplying $\langle\overline{\psi}\psi\rangle/M^3$ and terms in $\alpha_s$ multiplying $\langle \alpha_s G^2 \rangle/M^4$, where $M$ is the heavy-quark mass. Using lattice QCD results for heavy-strange correlators obtained for a variety of heavy quark masses on gluon field configurations including $u$, $d$ and $s$ quarks in the sea at three values of the lattice spacing, we are able to show that the contribution of the strange-quark condensate to the time-moments is very substantial. We use our lattice QCD time-moments and the OPE to determine a value for the condensate, fitting the 4th, 6th, 8th and 10th time-moments simultaneously. Our result, $\langle \overline{s}s \rangle^{\overline{\text{MS}}}(2 \text{GeV}) = -(296(11) \,\mathrm{MeV})^3$, agrees well with HPQCD's earlier, more direct, lattice QCD determination~\cite{McNeile:2012xh}. As well as confirming that the $s$ quark condensate is close in value to the light quark condensate, this demonstrates clearly the consistency of the Operator Product Expansion for fully nonperturbative calculations of matrix elements of short-distance operators in lattice QCD., Comment: 14 pages, 4 figures

We introduce a new method for calculating the ${\rm O}(\alpha^3)$ hadronic-vacuum-polarization contribution to the muon anomalous magnetic moment from ${ab-initio}$ lattice QCD. We first derive expressions suitable for computing the higher-order contributions either from the renormalized vacuum polarization function $\hat\Pi(q^2)$, or directly from the lattice vector-current correlator in Euclidean space. We then demonstrate the approach using previously-published results for the Taylor coefficients of $\hat\Pi(q^2)$ that were obtained on four-flavor QCD gauge-field configurations with physical light-quark masses. We obtain $10^{10} a_\mu^{\rm HVP,HO} = -9.3(1.3)$, in agreement with, but with a larger uncertainty than, determinations from $e^+e^- \to {\rm hadrons}$ data plus dispersion relations., Comment: Expanded and clarified discussion and revised Figure 4. Results unchanged. 11 pages, 5 tables, 5 figures. Version accepted to Physical Review D

DTAS is a segmented total absorption {\gamma}-ray spectrometer developed for the DESPEC experiment at FAIR. It is composed of up to eighteen NaI(Tl) crystals. In this work we study the performance of this detector with laboratory sources and also under real experimental conditions. We present a procedure to reconstruct offline the sum of the energy deposited in all the crystals of the spectrometer, which is complicated by the effect of NaI(Tl) light-yield non-proportionality. The use of a system to correct for time variations of the gain in individual detector modules, based on a light pulse generator, is demonstrated. We describe also an event-based method to evaluate the summing-pileup electronic distortion in segmented spectrometers. All of this allows a careful characterization of the detector with Monte Carlo simulations that is needed to calculate the response function for the analysis of total absorption {\gamma}-ray spectroscopy data. Special attention was paid to the interaction of neutrons with the spectrometer, since they are a source of contamination in studies of \b{eta}-delayed neutron emitting nuclei.

Background: $\beta$-delayed multiple neutron emission has been observed for some nuclei with A$\leq$100, with $^{100}$Rb being the heaviest $\beta$2n emitter measured to date. So far, only 25 P$_{2n}$ values have been determined for the $\sim$300 nuclei that may decay in this way. Accordingly, it is of interest to measure P$_{2n}$ values for the other possible multiple neutron emitters throughout the chart of the nuclides. It is of particular interest to make such measurement for nuclei with A$>$100 to test the predictions of theoretical models and simulation tools for the decays of heavy nuclei in the region of very neutron-rich nuclei. In addition, the decay properties of these nuclei are fundamental for the understanding of astrophysical nucleosynthesis processes such as the $r$-process, and safety inputs for nuclear reactors. Purpose: To determine for the first time the two neutron branching ratio, P$_{2n}$ value, for $^{136}$Sb through a direct neutron measurement, and to provide precise P$_{1n}$ values for $^{136}$Sb and $^{136}$Te. Method: Pure beams were provided by the JYFLTRAP at the IGISOL facility of the University of Jyv\"askyl\"a, Finland. The purified ions were implanted into a moving tape at the end of the beam line. The detection setup consisted of a plastic scintillator placed right behind the implantation point, and the BELEN detector, based on neutron counters embedded in a polyethylene matrix. The analysis was based on the study of the $\beta$- and neutron- growth-and-decay curves and the $\beta$-one-neutron and $\beta$-two-neutron time correlations. Results: The P$_{2n}$ value of $^{136}$Sb was found to be 0.14(3)\% and the measured P$_{1n}$ values for $^{136}$Sb and $^{136}$Te were found to be 32.2(15)\% and 1.47(6)\%, respectively. The measured P$_{2n}$ value is a factor 44 smaller than predicted by the FRDM+QRPA model used for $r$-process calculations., Comment: 11 pages, 10 figures