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
-
more » « less
Abstract Amorphous ices are usually classified as belonging to low-density or high-density amorphous ice (LDA and HDA) with densities
ρ L D A ≈ 0.94 g/cm3andρ H D A ≈ 1.15−1.17 g/cm3. However, a recent experiment crushing hexagonal ice (ball-milling) produced amedium -density amorphous ice (MDA,ρ M D A ≈ 1.06 g/cm3) adding complexity to our understanding of amorphous ice and the phase diagram of supercooled water. Motivated by the discovery of MDA, we perform computer simulations where amorphous ices are produced by isobaric cooling and isothermal compression/decompression. Our results show that, depending on the pressure employed, isobaric cooling can generate a continuum of amorphous ices with densities that expand in between those of LDA and HDA (briefly, intermediate amorphous ices, IA). In particular, the IA generated atP ≈ 125 MPa has a remarkably similar density and average structure as MDA, implying that MDA is not unique. Using the potential energy landscape formalism, we provide an intuitive qualitative understanding of the nature of LDA, HDA, and the IA generated at different pressures. In this view, LDA and HDA occupy specific and well-separated regions of the PEL; the IA prepared atP = 125 MPa is located in the intermediate region of the PEL that separates LDA and HDA. -
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
-
A series of experiments were conducted to investigate the flow field of a top-heated liquid fuel adjacent to an ice block. The experimental setup consisted of a borosilicate container containing an ice wall and a layer of n-heptane heated from above. Particle Image Velocimetry (PIV) and Background Oriented Schlieren (BOS) measurements were conducted on the liquid -phase. PIV measurements showed a surface flow toward the ice caused by surface -tension forces, which is driven by the horizontal temperature gradients on the liquid surface. A recirculation zone was observed under the free surface and near the ice. The combination of the two flow patterns caused lateral intrusion in the ice, instead of a uniform melting across ice surface. BOS measurements indicated presence of density gradients below the free surface of n-heptane and in regions near the ice block. These density gradients were created by local small-scale temperature gradients. The current experiments were conducted to explore the processes that influence the ice melting by immiscible liquid layers.more » « less
