skip to main content

Explore Research Products in the NSF-PAR

Title: Nuclear quantum effects on the dynamics and glass behavior of a monatomic liquid with two liquid states
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
Award ID(s):
2112550
NSF-PAR ID:
10352965
Author(s) / Creator(s):
; ;
Date Published:
Journal Name:
The Journal of Chemical Physics
Volume:
156
Issue:
20
ISSN:
0021-9606
Page Range / eLocation ID:
204502
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. ; ; ( , Physical Chemistry Chemical Physics)
    Experimental techniques, such as cryo-electron microscopy, require biological samples to be recovered at cryogenic temperatures ( T ≈ 100 K) with water being in an amorphous ice state. However, (bulk) water can exist in two amorphous ices at P < 1 GPa, low-density amorphous (LDA) ice at low pressures and high-density amorphous ice (HDA) at high pressures; HDA is ≈20–25% denser than LDA. While fast/plunge cooling at 1 bar brings the sample into LDA, high-pressure cooling (HPC), at sufficiently high pressure, produces HDA. HDA can also be produced by isothermal compression of LDA at cryogenic temperatures. Here, we perform classical molecular dynamics simulations to study the effects of LDA, HDA, and the LDA–HDA transformation on the structure and hydration of a small peptide, polyalanine. We follow thermodynamic paths corresponding to (i) fast/plunge cooling at 1 bar, (ii) HPC at P = 400 MPa, and (iii) compression/decompression cycles at T = 80 K. While process (i) produced LDA in the system, path (iii) produces HDA. Interestingly, the amorphous ice produced in process (ii) is an intermediate amorphous ice (IA) with properties that fall in-between those of LDA and HDA. Remarkably, the structural changes in polyalanine are negligible at all conditions studied (0–2000 MPa, 80–300 K) even when water changes among the low and high-density liquid states as well as the amorphous solids LDA, IA, and HDA. The similarities and differences in the hydration of polyalanine vitrified in LDA, IA, and HDA are described. Since the studied thermodynamic paths are suitable for the cryopreservation of biomolecules, we also study the structure and hydration of polyalanine along isobaric and isochoric heating paths, which can be followed experimentally for the recovery of cryopreserved samples. Upon heating, the structure of polyalanine remains practically unchanged. We conclude with a brief discussion of the practical advantages of (a) using HDA and IA as a cryoprotectant environment (as opposed to LDA), and (b) the use of isochoric heating as a recovery process (as opposed to isobaric heating). 
    more » « less
  2. ; ; ( , Scientific Reports)
    Abstract

    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
  3. ; ; ( , The Journal of Chemical Physics)
    As a liquid approaches the gas state, the properties of the potential energy landscape (PEL) sampled by the system become anomalous. Specifically, (i) the mechanically stable local minima of the PEL [inherent structures (IS)] can exhibit cavitation above the so-called Sastry volume, v S , before the liquid enters the gas phase. In addition, (ii) the pressure of the liquid at the sampled IS [i.e., the PEL equation of state, P IS ( v)] develops a spinodal-like minimum at v S . We perform molecular dynamics simulations of a monatomic water-like liquid and verify that points (i) and (ii) hold at high temperatures. However, at low temperatures, cavitation in the liquid and the corresponding IS occurs simultaneously and a Sastry volume cannot be defined. Remarkably, at intermediate/high temperatures, the IS of the liquid can exhibit crystallization, i.e., the liquid regularly visits the regions of the PEL that belong to the crystal state. The model liquid considered also exhibits a liquid–liquid phase transition (LLPT) between a low-density and a high-density liquid (LDL and HDL). By studying the behavior of P IS ( v) during the LLPT, we identify a Sastry volume for both LDL and HDL. The HDL Sastry volume marks the onset above which IS are heterogeneous (composed of LDL and HDL particles), analogous to points (i) and (ii) above. However, the relationship between the LDL Sastry volume and the onset of heterogeneous IS is less evident. We conclude by presenting a thermodynamic argument that can explain the behavior of the PEL equation of state P IS ( v) across both the liquid–gas phase transition and LLPT. 
    more » « less
  4. ; ; ( , The Journal of Chemical Physics)
    A fundamental and much-debated issue in glass science is the existence and nature of liquid–liquid transitions in glass-forming liquids. Here, we report the existence of a novel reentrant structural transition in a S-rich arsenic sulfide liquid of composition As 2.5 S 97.5 . The nature of this transition and its effect on viscosity are investigated in situ using a combination of differential scanning calorimetry and simultaneous Raman spectroscopic and rheometric measurements. The results indicate that, upon heating significantly above its glass transition temperature (261 K), the constituent [Formula: see text] sulfur chains in the structure of the supercooled liquid first undergo a [Formula: see text] chain-to-ring conversion near ∼383 K, which is exothermic in nature. Further heating above 393 K alters the equilibrium to shift in the opposite direction toward an endothermic ring-to-chain conversion characteristic of the well-known λ-transition in pure sulfur liquid. This behavior is attributed to the competing effects of enthalpy of mixing and conformational entropy of ring and chain elements in the liquid. The existence of reentrant structural transitions in glass-forming liquids could provide important insights into the thermodynamics of liquid–liquid transitions and may have important consequences for harnessing novel functionalities of derived glasses. 
    more » « less
  5. ; ; ; ; ( , Scientific Reports)
    Abstract

    The atomic structure of a germanium doped phosphorous selenide glass of composition Ge2.8P57.7Se39.5is determined as a function of pressure from ambient to 24 GPa using Monte-Carlo simulations constrained by high energy x-ray scattering data. The ambient pressure structure consists primarily of P4Se3molecules and planar edge shared phosphorus rings, reminiscent of those found in red phosphorous as well as a small fraction of locally clustered corner-sharing GeSe4tetrahedra. This low-density amorphous phase transforms into a high-density amorphous phase at ~6.3 GPa. The high-pressure phase is characterized by an extended network structure. The polyamorphic transformation between these two phases involves opening of the P3ring at the base of the P4Se3molecules and subsequent reaction with red phosphorus type moieties to produce a cross linked structure. The compression mechanism of the low-density phase involves increased molecular packing, whereas that of the high pressure phase involves an increase in the nearest-neighbor coordination number while the bond angle distributions broaden and shift to smaller angles. The entropy and volume changes associated with this polyamorphic transformation are positive and negative, respectively, and consequently the corresponding Clapeyron slope for this transition would be negative. This result has far reaching implications in our current understanding of the thermodynamics of polyamorphic transitions in glasses and glass-forming liquids.

     
    more » « less
  • Citation Formats
  • MLA
  • APA
  • Chicago
  • BibTeX
  • Save / Share this Record
  • Send to Email