skip to main content

Title: Evidence of a liquid–liquid phase transition in H$$_2$$O and D$$_2$$O from path-integral molecular dynamics simulations

We perform path-integral molecular dynamics (PIMD), ring-polymer MD (RPMD), and classical MD simulations of H$$_2$$2O and D$$_2$$2O using the q-TIP4P/F water model over a wide range of temperatures and pressures. The density$$\rho (T)$$ρ(T), isothermal compressibility$$\kappa _T(T)$$κT(T), and self-diffusion coefficientsD(T) of H$$_2$$2O and D$$_2$$2O are in excellent agreement with available experimental data; the isobaric heat capacity$$C_P(T)$$CP(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$$2O and D$$_2$$2O exhibit a liquid-liquid critical point (LLCP) at low temperatures and positive pressures. The data from PIMD/MD for H$$_2$$2O and D$$_2$$2O 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$$2O, from PIMD simulations, is located at$$P_c = 167 \pm 9$$Pc=167±9 MPa,$$T_c = 159 \pm 6$$Tc=159±6 K, and$$\rho _c = 1.02 \pm 0.01$$ρc=1.02±0.01 g/cm$$^3$$3. Isotope substitution effects are important; the LLCP location in q-TIP4P/F D$$_2$$2O is estimated to be$$P_c = 176 \pm 4$$Pc=176±4 MPa,$$T_c = 177 \pm 2$$Tc=177±2 K, and$$\rho _c = 1.13 \pm 0.01$$ρc=1.13±0.01 g/cm$$^3$$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$$Pc=203±4 MPa,$$T_c = 175 \pm 2$$Tc=175±2 K, and$$\rho _c = 1.03 \pm 0.01$$ρc=1.03±0.01 g/cm$$^3$$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$$Tcfor D$$_2$$2O and, particularly, H$$_2$$2O suggest that improved water models are needed for the study of supercooled water.

more » « less
Award ID(s):
Author(s) / Creator(s):
; ;
Publisher / Repository:
Nature Publishing Group
Date Published:
Journal Name:
Scientific Reports
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract

    We report on a measurement of Spin Density Matrix Elements (SDMEs) in hard exclusive$$\rho ^0$$ρ0meson muoproduction at COMPASS using 160 GeV/cpolarised$$ \mu ^{+}$$μ+and$$ \mu ^{-}$$μ-beams impinging on a liquid hydrogen target. The measurement covers the kinematic range 5.0 GeV/$$c^2$$c2$$< W<$$<W<17.0 GeV/$$c^2$$c2, 1.0 (GeV/c)$$^2$$2$$< Q^2<$$<Q2<10.0 (GeV/c)$$^2$$2and 0.01 (GeV/c)$$^2$$2$$< p_{\textrm{T}}^2<$$<pT2<0.5 (GeV/c)$$^2$$2. Here,Wdenotes the mass of the final hadronic system,$$Q^2$$Q2the virtuality of the exchanged photon, and$$p_{\textrm{T}}$$pTthe transverse momentum of the$$\rho ^0$$ρ0meson with respect to the virtual-photon direction. The measured non-zero SDMEs for the transitions of transversely polarised virtual photons to longitudinally polarised vector mesons ($$\gamma ^*_T \rightarrow V^{ }_L$$γTVL) indicate a violation ofs-channel helicity conservation. Additionally, we observe a dominant contribution of natural-parity-exchange transitions and a very small contribution of unnatural-parity-exchange transitions, which is compatible with zero within experimental uncertainties. The results provide important input for modelling Generalised Parton Distributions (GPDs). In particular, they may allow one to evaluate in a model-dependent way the role of parton helicity-flip GPDs in exclusive$$\rho ^0$$ρ0production.

    more » « less
  2. Abstract

    Initially, vanadium dioxide seems to be an ideal first-order phase transition case study due to its deceptively simple structure and composition, but upon closer inspection there are nuances to the driving mechanism of the metal-insulator transition (MIT) that are still unexplained. In this study, a local structure analysis across a bulk powder tungsten-substitution series is utilized to tease out the nuances of this first-order phase transition. A comparison of the average structure to the local structure using synchrotron x-ray diffraction and total scattering pair-distribution function methods, respectively, is discussed as well as comparison to bright field transmission electron microscopy imaging through a similar temperature-series as the local structure characterization. Extended x-ray absorption fine structure fitting of thin film data across the substitution-series is also presented and compared to bulk. Machine learning technique, non-negative matrix factorization, is applied to analyze the total scattering data. The bulk MIT is probed through magnetic susceptibility as well as differential scanning calorimetry. The findings indicate the local transition temperature ($$T_c$$Tc) is less than the average$$T_c$$Tcsupporting the Peierls-Mott MIT mechanism, and demonstrate that in bulk powder and thin-films, increasing tungsten-substitution instigates local V-oxidation through the phase pathway VO$$_2\, \rightarrow$$2V$$_6$$6O$$_{13} \, \rightarrow$$13V$$_2$$2O$$_5$$5.

    more » « less
  3. Abstract

    We study the structure of the Liouville quantum gravity (LQG) surfaces that are cut out as one explores a conformal loop-ensemble$$\hbox {CLE}_{\kappa '}$$CLEκfor$$\kappa '$$κin (4, 8) that is drawn on an independent$$\gamma $$γ-LQG surface for$$\gamma ^2=16/\kappa '$$γ2=16/κ. The results are similar in flavor to the ones from our companion paper dealing with$$\hbox {CLE}_{\kappa }$$CLEκfor$$\kappa $$κin (8/3, 4), where the loops of the CLE are disjoint and simple. In particular, we encode the combined structure of the LQG surface and the$$\hbox {CLE}_{\kappa '}$$CLEκin terms of stable growth-fragmentation trees or their variants, which also appear in the asymptotic study of peeling processes on decorated planar maps. This has consequences for questions that do a priori not involve LQG surfaces: In our paper entitled “CLE Percolations” described the law of interfaces obtained when coloring the loops of a$$\hbox {CLE}_{\kappa '}$$CLEκindependently into two colors with respective probabilitiespand$$1-p$$1-p. This description was complete up to one missing parameter$$\rho $$ρ. The results of the present paper about CLE on LQG allow us to determine its value in terms ofpand$$\kappa '$$κ. It shows in particular that$$\hbox {CLE}_{\kappa '}$$CLEκand$$\hbox {CLE}_{16/\kappa '}$$CLE16/κare related via a continuum analog of the Edwards-Sokal coupling between$$\hbox {FK}_q$$FKqpercolation and theq-state Potts model (which makes sense even for non-integerqbetween 1 and 4) if and only if$$q=4\cos ^2(4\pi / \kappa ')$$q=4cos2(4π/κ). This provides further evidence for the long-standing belief that$$\hbox {CLE}_{\kappa '}$$CLEκand$$\hbox {CLE}_{16/\kappa '}$$CLE16/κrepresent the scaling limits of$$\hbox {FK}_q$$FKqpercolation and theq-Potts model whenqand$$\kappa '$$κare related in this way. Another consequence of the formula for$$\rho (p,\kappa ')$$ρ(p,κ)is the value of half-plane arm exponents for such divide-and-color models (a.k.a. fuzzy Potts models) that turn out to take a somewhat different form than the usual critical exponents for two-dimensional models.

    more » « less
  4. Abstract

    We provide moment bounds for expressions of the type$$(X^{(1)} \otimes \cdots \otimes X^{(d)})^T A (X^{(1)} \otimes \cdots \otimes X^{(d)})$$(X(1)X(d))TA(X(1)X(d))where$$\otimes $$denotes the Kronecker product and$$X^{(1)}, \ldots , X^{(d)}$$X(1),,X(d)are random vectors with independent, mean 0, variance 1, subgaussian entries. The bounds are tight up to constants depending ondfor the case of Gaussian random vectors. Our proof also provides a decoupling inequality for expressions of this type. Using these bounds, we obtain new, improved concentration inequalities for expressions of the form$$\Vert B (X^{(1)} \otimes \cdots \otimes X^{(d)})\Vert _2$$B(X(1)X(d))2.

    more » « less
  5. Abstract

    We present the first unquenched lattice-QCD calculation of the form factors for the decay$$B\rightarrow D^*\ell \nu $$BDνat nonzero recoil. Our analysis includes 15 MILC ensembles with$$N_f=2+1$$Nf=2+1flavors of asqtad sea quarks, with a strange quark mass close to its physical mass. The lattice spacings range from$$a\approx 0.15$$a0.15fm 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}|$$|Vcb|. 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}$$Vcb=(38.40±0.68th±0.34exp±0.18EM)×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$$χ2/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$$R(D)=0.265±0.013, which confirms the current tension between theory and experiment.

    more » « less