An official website of the United States government
We perform path integral molecular dynamics (PIMD) simulations of a monatomic liquid that exhibits a liquid–liquid phase transition and liquid–liquid critical point. PIMD simulations are performed using different values of Planck’s constant h, allowing us to study the behavior of the liquid as nuclear quantum effects (NQE, i.e., atoms delocalization) are introduced, from the classical liquid ( h = 0) to increasingly quantum liquids ( h > 0). By combining the PIMD simulations with the ring-polymer molecular dynamics method, we also explore the dynamics of the classical and quantum liquids. We find that (i) the glass transition temperature of the low-density liquid (LDL) is anomalous, i.e., [Formula: see text] decreases upon compression. Instead, (ii) the glass transition temperature of the high-density liquid (HDL) is normal, i.e., [Formula: see text] increases upon compression. (iii) NQE shift both [Formula: see text] and [Formula: see text] toward lower temperatures, but NQE are more pronounced on HDL. We also study the glass behavior of the ring-polymer systems associated with the quantum liquids studied (via the path-integral formulation of statistical mechanics). There are two glass states in all the systems studied, low-density amorphous ice (LDA) and high-density amorphous ice (HDA), which are the glass counterparts of LDL and HDL. In all cases, the pressure-induced LDA–HDA transformation is sharp, reminiscent of a first-order phase transition. In the low-quantum regime, the LDA–HDA transformation is reversible, with identical LDA forms before compression and after decompression. However, in the high-quantum regime, the atoms become more delocalized in the final LDA than in the initial LDA, raising questions on the reversibility of the LDA–HDA transformation.
more »
« less
Potential energy landscape formalism for quantum molecular liquids
Abstract The potential energy landscape (PEL) formalism is a powerful tool within statistical mechanics to study the thermodynamic properties of classical low-temperature liquids and glasses. Recently, the PEL formalism has been extended to liquids/glasses that obey quantum mechanics, but applications have been limited to atomistic model liquids. In this work, we extend the PEL formalism to liquid/glassy water using path-integral molecular dynamics (PIMD) simulations, where nuclear quantum effects (NQE) are included. Our PIMD simulations, based on the q-TIP4P/F water model, show that the PEL of quantum water is both Gaussian and anharmonic. Importantly, the ring-polymers associated to the O/H atoms in the PIMD simulations, collapse at the local minima of the PEL (inherent structures, IS) for both liquid and glassy states. This allows us to calculate, analytically, the IS vibrational density of states (IS-VDOS) of the ring-polymer system using the IS-VDOS of classical water (obtained from classical MD simulations). The role of NQE on the structural properties of liquid/glassy water at various pressures are discussed in detail. Overall, our results demonstrate that the PEL formalism can effectively describe the behavior of molecular liquids at low temperatures and in the glass states, regardless of whether the liquid/glass obeys classical or quantum mechanics.
more »
« less
- PAR ID:
- 10558468
- Publisher / Repository:
- Nature Publishing Group
- Date Published:
- Journal Name:
- Communications Chemistry
- Volume:
- 7
- Issue:
- 1
- ISSN:
- 2399-3669
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
The potential energy landscape (PEL) formalism is a tool within statistical mechanics that has been used in the past to calculate the equation of states (EOS) of classical rigid model liquids at low temperatures, where computer simulations may be challenging. In this work, we use classical molecular dynamics (MD) simulations and the PEL formalism to calculate the EOS of the flexible q-TIP4P/F water model. This model exhibits a liquid–liquid critical point (LLCP) in the supercooled regime, at (Pc = 150 MPa, Tc = 190 K, and ρc = 1.04 g/cm3) [using the reaction field technique]. The PEL-EOS of q-TIP4P/F water and the corresponding location of the LLCP are in very good agreement with the MD simulations. We show that the PEL of q-TIP4P/F water is Gaussian, which allows us to calculate the configurational entropy of the system, Sconf. The Sconf of q-TIP4P/F water is surprisingly similar to that reported previously for rigid water models, suggesting that intramolecular flexibility does not necessarily add roughness to the PEL. We also show that the Adam–Gibbs relation, which relates the diffusion coefficient D with Sconf, holds for the flexible q-TIP4P/F water model. Overall, our results indicate that the PEL formalism can be used to study molecular systems that include molecular flexibility, the common case in standard force fields. This is not trivial since the introduction of large bending/stretching mode frequencies is problematic in classical statistical mechanics. For example, as shown previously, we find that such high frequencies lead to unphysical (negative) entropy for q-TIP4P/F water when using classical statistical mechanics (yet, the PEL formalism can be applied successfully).more » « less
-
We perform classical molecular dynamics (MD) and path-integral MD (PIMD) simulations of H2O and D2O using the q-TIP4P/F model over a wide range of temperatures and pressures to study the nuclear quantum effects (NQEs) on (i) the vitrification of liquid water upon isobaric cooling at different pressures and (ii) pressure-induced transformations at constant temperature between low-density amorphous and high-density amorphous ice (LDA and HDA) and hexagonal ice Ih and HDA. Upon isobaric cooling, classical and quantum H2O and D2O vitrify into a continuum of intermediate amorphous ices (IA), with densities in-between those of LDA and HDA (depending on pressure). Importantly, the density of the IA varies considerably if NQEs are included (similar conclusions hold for ice Ih at all pressures studied). While the structure of the IA is not very sensitive to NQE, the geometry of the hydrogen-bond (HB) is. NQE leads to longer and less linear HB in LDA, HDA, and ice Ih than found in the classical case. Interestingly, the delocalization of the H/D atoms is non-negligible and identical in LDA, HDA, and ice Ih at all pressures studied. Our isothermal compression/decompression MD/PIMD simulations show that classical and quantum H2O and D2O all exhibit LDA–HDA and ice Ih-HDA transformations, consistent with experiments. The inclusion of NQE leads to a softer HB-network, which lowers slightly the LDA/ice Ih-to-HDA transformation pressures. Interestingly, the HB in HDA is longer and less linear than in LDA, which is counterintuitive given that HDA is ≈25% denser than LDA. Overall, our results show that, while classical computer simulations provide the correct qualitative phenomenology of ice and glassy water, NQEs are necessary for a quantitative description.more » « less
-
The molecular origins of water’s anomalous properties have long been a subject of scientific inquiry. The liquid–liquid phase transition hypothesis, which posits the existence of distinct low-density and high-density liquid states separated by a first-order phase transition terminating at a critical point, has gained increasing experimental and computational support and offers a thermodynamically consistent framework for many of water’s anomalies. However, experimental challenges in avoiding crystallization near the postulated liquid–liquid critical point have focused attention to water’s canonical glassy states: low-density and high-density amorphous ice. Here, we use two Deep Potential machine-learning models, trained on the Strongly Constrained and Appropriately Normed density functional and the highly accurate Many-Body Polarizable potential, to conduct an investigation of water’s glassy phenomenology based on quantum mechanical calculations. Despite not being explicitly trained on amorphous ices, both models accurately capture the structure and transformation of the water glasses, including their interconversion along different thermodynamic paths. Isobaric quenching of liquid water at various pressures generates a continuum of intermediate amorphous ices and density fluctuations increase near the liquid–liquid critical pressure. The glass transition temperatures of the amorphous ices produced at different pressures exhibit two distinct branches, corresponding to low-density and high-density amorphous ice behaviors, consistent with experiment and the liquid–liquid transition hypothesis. Extrapolating transformation pressures from isothermal compressions to experimental compression rates brings our simulations into excellent agreement with data. Our findings demonstrate that machine-learning potentials trained on equilibrium phases can effectively model nonequilibrium glassy behavior and pave the way for studying long-timescale, out-of-equilibrium processes with quantum mechanical accuracy.more » « less
-
The discovery of ultrastable glasses raises novel challenges about glassy systems. Recent experiments studied the macroscopic devitrification of ultrastable glasses into liquids upon heating but lacked microscopic resolution. We use molecular dynamics simulations to analyze the kinetics of this transformation. In the most stable systems, devitrification occurs after a very large time, but the liquid emerges in two steps. At short times, we observe the rare nucleation and slow growth of isolated droplets containing a liquid maintained under pressure by the rigidity of the surrounding glass. At large times, pressure is released after the droplets coalesce into large domains, which accelerates devitrification. This two-step process produces pronounced deviations from the classical Avrami kinetics and explains the emergence of a giant lengthscale characterizing the devitrification of bulk ultrastable glasses. Our study elucidates the nonequilibrium kinetics of glasses following a large temperature jump, which differs from both equilibrium relaxation and aging dynamics, and will guide future experimental studies.more » « less
