Abstract We perform path-integral molecular dynamics (PIMD), ring-polymer MD (RPMD), and classical MD simulations of H$$_2$$ O and D$$_2$$ O using the q-TIP4P/F water model over a wide range of temperatures and pressures. The density$$\rho (T)$$ , isothermal compressibility$$\kappa _T(T)$$ , and self-diffusion coefficientsD(T) of H$$_2$$ O and D$$_2$$ O are in excellent agreement with available experimental data; the isobaric heat capacity$$C_P(T)$$ obtained from PIMD and MD simulations agree qualitatively well with the experiments. Some of these thermodynamic properties exhibit anomalous maxima upon isobaric cooling, consistent with recent experiments and with the possibility that H$$_2$$ O and D$$_2$$ O exhibit a liquid-liquid critical point (LLCP) at low temperatures and positive pressures. The data from PIMD/MD for H$$_2$$ O and D$$_2$$ O can be fitted remarkably well using the Two-State-Equation-of-State (TSEOS). Using the TSEOS, we estimate that the LLCP for q-TIP4P/F H$$_2$$ O, from PIMD simulations, is located at$$P_c = 167 \pm 9$$ MPa,$$T_c = 159 \pm 6$$ K, and$$\rho _c = 1.02 \pm 0.01$$ g/cm$$^3$$ . Isotope substitution effects are important; the LLCP location in q-TIP4P/F D$$_2$$ O is estimated to be$$P_c = 176 \pm 4$$ MPa,$$T_c = 177 \pm 2$$ K, and$$\rho _c = 1.13 \pm 0.01$$ g/cm$$^3$$ . Interestingly, for the water model studied, differences in the LLCP location from PIMD and MD simulations suggest that nuclear quantum effects (i.e., atoms delocalization) play an important role in the thermodynamics of water around the LLCP (from the MD simulations of q-TIP4P/F water,$$P_c = 203 \pm 4$$ MPa,$$T_c = 175 \pm 2$$ K, and$$\rho _c = 1.03 \pm 0.01$$ g/cm$$^3$$ ). Overall, our results strongly support the LLPT scenario to explain water anomalous behavior, independently of the fundamental differences between classical MD and PIMD techniques. The reported values of$$T_c$$ for D$$_2$$ O and, particularly, H$$_2$$ O suggest that improved water models are needed for the study of supercooled water.
more »
« less
This content will become publicly available on August 1, 2026
Performance of an optical TPC Geant4 simulation with opticks GPU-accelerated photon propagation
Abstract We investigate the performance of , a 7.5 GPU-accelerated photon propagation tool compared with a single-threaded simulation. We compare the simulations using an improved model of the gaseous time projection chamber. Performance results suggest that improves simulation speeds by between$$58.47\pm {0.02}$$ and$$181.39\pm {0.28}$$ times relative to a CPU-only simulation and these results vary between different types of GPU and CPU. A detailed comparison shows that the number of detected photons, along with their times and wavelengths, are in good agreement between and .
more »
« less
- Award ID(s):
- 2004111
- PAR ID:
- 10645143
- Author(s) / Creator(s):
- ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; more »
- Publisher / Repository:
- Springer Nature
- Date Published:
- Journal Name:
- The European Physical Journal C
- Volume:
- 85
- Issue:
- 8
- ISSN:
- 1434-6052
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
Abstract We present the first unquenched lattice-QCD calculation of the form factors for the decay$$B\rightarrow D^*\ell \nu $$ at nonzero recoil. Our analysis includes 15 MILC ensembles with$$N_f=2+1$$ flavors of asqtad sea quarks, with a strange quark mass close to its physical mass. The lattice spacings range from$$a\approx 0.15$$ fm down to 0.045 fm, while the ratio between the light- and the strange-quark masses ranges from 0.05 to 0.4. The valencebandcquarks are treated using the Wilson-clover action with the Fermilab interpretation, whereas the light sector employs asqtad staggered fermions. We extrapolate our results to the physical point in the continuum limit using rooted staggered heavy-light meson chiral perturbation theory. Then we apply a model-independent parametrization to extend the form factors to the full kinematic range. With this parametrization we perform a joint lattice-QCD/experiment fit using several experimental datasets to determine the CKM matrix element$$|V_{cb}|$$ . We obtain$$\left| V_{cb}\right| = (38.40 \pm 0.68_{\text {th}} \pm 0.34_{\text {exp}} \pm 0.18_{\text {EM}})\times 10^{-3}$$ . The first error is theoretical, the second comes from experiment and the last one includes electromagnetic and electroweak uncertainties, with an overall$$\chi ^2\text {/dof} = 126/84$$ , which illustrates the tensions between the experimental data sets, and between theory and experiment. This result is in agreement with previous exclusive determinations, but the tension with the inclusive determination remains. Finally, we integrate the differential decay rate obtained solely from lattice data to predict$$R(D^*) = 0.265 \pm 0.013$$ , which confirms the current tension between theory and experiment.more » « less
-
A<sc>bstract</sc> We present a study of$$ {\Xi}_c^0\to {\Xi}^0{\pi}^0 $$ ,$$ {\Xi}_c^0\to {\Xi}^0\eta $$ , and$$ {\Xi}_c^0\to {\Xi}^0{\eta}^{\prime } $$ decays using the Belle and Belle II data samples, which have integrated luminosities of 980 fb−1and 426 fb−1, respectively. We measure the following relative branching fractions$$ {\displaystyle \begin{array}{c}\mathcal{B}\left({\Xi}_c^0\to {\Xi}^0{\pi}^0\right)/\mathcal{B}\left({\Xi}_c^0\to {\Xi}^{-}{\pi}^{+}\right)=0.48\pm 0.02\left(\textrm{stat}\right)\pm 0.03\left(\textrm{syst}\right),\\ {}\mathcal{B}\left({\Xi}_c^0\to {\Xi}^0\eta \right)/\mathcal{B}\left({\Xi}_c^0\to {\Xi}^{-}{\pi}^{+}\right)=0.11\pm 0.01\left(\textrm{stat}\right)\pm 0.01\left(\textrm{syst}\right),\\ {}\mathcal{B}\left({\Xi}_c^0\to {\Xi}^0{\eta}^{\prime}\right)/\mathcal{B}\left({\Xi}_c^0\to {\Xi}^{-}{\pi}^{+}\right)=0.08\pm 0.02\left(\textrm{stat}\right)\pm 0.01\left(\textrm{syst}\right)\end{array}} $$ for the first time, where the uncertainties are statistical (stat) and systematic (syst). By multiplying by the branching fraction of the normalization mode,$$ \mathcal{B}\left({\Xi}_c^0\to {\Xi}^{-}{\pi}^{+}\right) $$ , we obtain the following absolute branching fraction results$$ {\displaystyle \begin{array}{c}\mathcal{B}\left({\Xi}_c^0\to {\Xi}^0{\pi}^0\right)=\left(6.9\pm 0.3\left(\textrm{stat}\right)\pm 0.5\left(\textrm{syst}\right)\pm 1.3\left(\operatorname{norm}\right)\right)\times {10}^{-3},\\ {}\mathcal{B}\left({\Xi}_c^0\to {\Xi}^0\eta \right)=\left(1.6\pm 0.2\left(\textrm{stat}\right)\pm 0.2\left(\textrm{syst}\right)\pm 0.3\left(\operatorname{norm}\right)\right)\times {10}^{-3},\\ {}\mathcal{B}\left({\varXi}_c^0\to {\Xi}^0{\eta}^{\prime}\right)=\left(1.2\pm 0.3\left(\textrm{stat}\right)\pm 0.1\left(\textrm{syst}\right)\pm 0.2\left(\operatorname{norm}\right)\right)\times {10}^{-3},\end{array}} $$ where the third uncertainties are from$$ \mathcal{B}\left({\Xi}_c^0\to {\Xi}^{-}{\pi}^{+}\right) $$ . The asymmetry parameter for$$ {\Xi}_c^0\to {\Xi}^0{\pi}^0 $$ is measured to be$$ \alpha \left({\Xi}_c^0\to {\Xi}^0{\pi}^0\right)=-0.90\pm 0.15\left(\textrm{stat}\right)\pm 0.23\left(\textrm{syst}\right) $$ .more » « less
-
Abstract Radiogenic neutrons emitted by detector materials are one of the most challenging backgrounds for the direct search of dark matter in the form of weakly interacting massive particles (WIMPs). To mitigate this background, the XENONnT experiment is equipped with a novel gadolinium-doped water Cherenkov detector, which encloses the xenon dual-phase time projection chamber (TPC). The neutron veto (NV) can tag neutrons via their capture on gadolinium or hydrogen, which release$$\gamma $$ -rays that are subsequently detected as Cherenkov light. In this work, we present the first results of the XENONnT NV when operated with demineralized water only, before the insertion of gadolinium. Its efficiency for detecting neutrons is$$({82\pm 1}){\%}$$ , the highest neutron detection efficiency achieved in a water Cherenkov detector. This enables a high efficiency of$$({53\pm 3}){\%}$$ for the tagging of WIMP-like neutron signals, inside a tagging time window of$${250}~{\upmu }\hbox {s}$$ between TPC and NV, leading to a livetime loss of$${1.6}{\%}$$ during the first science run of XENONnT.more » « less
-
Abstract Ultra-pure NaI(Tl) crystals are the key element for a model-independent verification of the long standing DAMA result and a powerful means to search for the annual modulation signature of dark matter interactions. The SABRE collaboration has been developing cutting-edge techniques for the reduction of intrinsic backgrounds over several years. In this paper we report the first characterization of a 3.4 kg crystal, named NaI-33, performed in an underground passive shielding setup at LNGS. NaI-33 has a record low$$^{39}$$ K contamination of 4.3 ± 0.2 ppb as determined by mass spectrometry. We measured a light yield of 11.1 ± 0.2 photoelectrons/keV and an energy resolution of 13.2% (FWHM/E) at 59.5 keV. We evaluated the activities of$$^{226}$$ Ra and$$^{228}$$ Th inside the crystal to be$$5.9\pm 0.6~\upmu $$ Bq/kg and$$1.6\pm 0.3~\upmu $$ Bq/kg, respectively, which would indicate a contamination from$$^{238}$$ U and$$^{232}$$ Th at part-per-trillion level. We measured an activity of 0.51 ± 0.02 mBq/kg due to$$^{210}$$ Pb out of equilibrium and a$$\alpha $$ quenching factor of 0.63 ± 0.01 at 5304 keV. We illustrate the analyses techniques developed to reject electronic noise in the lower part of the energy spectrum. A cut-based strategy and a multivariate approach indicated a rate, attributed to the intrinsic radioactivity of the crystal, of$$\sim $$ 1 count/day/kg/keV in the [5–20] keV region.more » « less
An official website of the United States government
