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
This content will become publicly available on December 1, 2024
Liquid-liquid phase separation in supercooled water from ultrafast heating of low-density amorphous ice
Abstract Recent experiments continue to find evidence for a liquid-liquid phase transition (LLPT) in supercooled water, which would unify our understanding of the anomalous properties of liquid water and amorphous ice. These experiments are challenging because the proposed LLPT occurs under extreme metastable conditions where the liquid freezes to a crystal on a very short time scale. Here, we analyze models for the LLPT to show that coexistence of distinct high-density and low-density liquid phases may be observed by subjecting low-density amorphous (LDA) ice to ultrafast heating. We then describe experiments in which we heat LDA ice to near the predicted critical point of the LLPT by an ultrafast infrared laser pulse, following which we measure the structure factor using femtosecond x-ray laser pulses. Consistent with our predictions, we observe a LLPT occurring on a time scale < 100 ns and widely separated from ice formation, which begins at times >1 μs.
more »
« less
- Award ID(s):
- 2112550
- NSF-PAR ID:
- 10430913
- Date Published:
- Journal Name:
- Nature Communications
- Volume:
- 14
- Issue:
- 1
- ISSN:
- 2041-1723
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
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
-
A series of experiments were conducted to investigate the melting of ice adjacent to a water-immiscible liquid layer (n-dodecane) exposed to radiation from above. The experimental setup consisted of a borosilicate container containing an ice wall and a layer of n-dodecane heated from above. In addition to tracking the movement of the melt front, Particle Image Velocimetry (PIV) and Background Oriented Schlieren (BOS) measurements were conducted on the liquid-phase . Two distinct melting regimes were found to dominate the melting process. First was the uniform melting across the contact area with the immiscible liquid layer for low radiation levels (~1 kW/m 2 ). Second was the lateral intrusion regime, where a depression near free surface of the liquid forms in ice and grows laterally for radiation level greater than ~1 kW/m 2 . The ice surface remained flat and smooth in uniform melting regime, whereas in the lateral intrusion regime a series of rivulets were formed that carved valleys on the ice. PIV measurements showed a surface flow toward the ice for all heat flux levels caused by surface-tension forces. Increase of the heat flux levels caused a transition to multi-roll structure in the flow field. This multi-roll structure, which is accompanied by a recirculation zone near the ice, increased heat transfer coefficient near the surface of the liquid causing lateral intrusion regime. BOS measurements indicated presence of density gradients below the free surface of n-dodecane and in regions near ice that are caused by local small-scale temperature gradients. The current experiments were conducted to explore the melting dynamics and to shed light on the processes that influence the ice melting. Implications of such mechanisms in a real-life scenario, i.e. oil spill in ice-infested waters, needs to be explored further by using more liquids and improved accuracy with diagnostic techniques.more » « less
-
Abstract. Surface melting on the Antarctic Ice Sheet has been monitored by satellite microwave radiometry for over 40 years. Despite this long perspective, our understanding of the microwave emission from wet snow is still limited, preventing the full exploitation of these observations to study supraglacial hydrology. Using the Snow Microwave Radiative Transfer (SMRT) model, this study investigatesthe sensitivity of microwave brightness temperature to snow liquid water content at frequencies from 1.4 to 37 GHz. We first determine the snowpack properties for eight selected coastal sites byretrieving profiles of density, grain size and ice layers from microwave observations when the snowpack is dry during wintertime. Second, a series of brightness temperature simulations is run with added water. The results show that (i) a small quantity of liquid water (≈0.5 kg m−2) can be detected, but the actual quantity cannot be retrieved out of the full range of possible water quantities; (ii) the detection of a buried wet layer is possible up to a maximum depth of 1 to 6 m depending on the frequency (6–37 GHz) and on the snow properties (grain size, density) at each site; (iii) surface ponds and water-saturated areas may prevent melt detection, but the current coverage of these waterbodies in the large satellite field of view is presently too small in Antarctica to have noticeable effects; and (iv) at 1.4 GHz, while the simulations are less reliable, we found a weaker sensitivity to liquid water and the maximal depth of detection is relatively shallow (<10 m) compared to the typical radiation penetration depth in dry firn (≈1000 m) at this low frequency. These numerical results pave the way for the development of improved multi-frequency algorithms to detect melt intensity and the depth of liquid water below the surface in the Antarctic snowpack.more » « less
-
The transition between the gas-, supercritical-, and liquid-phase behavior is a fascinating topic, which still lacks molecular-level understanding. Recent ultrafast two-dimensional infrared spectroscopy experiments suggested that the vibrational spectroscopy of N 2 O embedded in xenon and SF 6 as solvents provides an avenue to characterize the transitions between different phases as the concentration (or density) of the solvent increases. The present work demonstrates that classical molecular dynamics (MD) simulations together with accurate interaction potentials allows us to (semi-)quantitatively describe the transition in rotational vibrational infrared spectra from the P-/R-branch line shape for the stretch vibrations of N 2 O at low solvent densities to the Q-branch-like line shapes at high densities. The results are interpreted within the classical theory of rigid-body rotation in more/less constraining environments at high/low solvent densities or based on phenomenological models for the orientational relaxation of rotational motion. It is concluded that classical MD simulations provide a powerful approach to characterize and interpret the ultrafast motion of solutes in low to high density solvents at a molecular level.more » « less
