skip to main content
US FlagAn official website of the United States government
dot gov icon
Official websites use .gov
A .gov website belongs to an official government organization in the United States.
https lock icon
Secure .gov websites use HTTPS
A lock ( lock ) or https:// means you've safely connected to the .gov website. Share sensitive information only on official, secure websites.


Title: Direct observation of liquid-liquid phase coexistence in deeply supercooled water using an accurate polarizable multipole model
Abstract Liquid water can be supercooled up to about 50~K below the melting point before undergoing homogeneous ice nucleation. Based on experimental thermodynamic observations and computer simulations it was hypothesized that below this temperature and at pressures of several kbar water undergoes a liquid-liquid phase transition (LLPT) and the transition line ends at a second critical point. However, challenges in experiments and simulations at such deep cooling leave doubts about the nature of the LLPT and the existence of the critical point.Here we use molecular dynamics simulations with a highly accurate and computationally efficient polarizable water model to establish the character of the LLPT and identify the location of the second critical point. Our microsecond-long simulations provide the first direct evidence of a well-defined moving interface between low-density and high-density water at conditions near the phase boundary. This is the ultimate proof of a first-order transition between two liquid phases with distinct free energy basins separated by a barrier, resolving a long-standing debate. These results provide new perspectives on supercooled water under pressure simulated with an accurate and realistic model suitable for studies of water in confined geological and biological environments.  more » « less
Award ID(s):
2053235
PAR ID:
10651804
Author(s) / Creator(s):
 ;  ;  
Publisher / Repository:
Research Square
Date Published:
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. 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
  2. The hypothesis that the anomalous behavior of liquid water is related to the existence of a second critical point in deeply supercooled states has long been the subject of intense debate. Recent, sophisticated experiments designed to observe the transformation between the two subcritical liquids on nano- and microsecond time scales, along with demanding numerical simulations based on classical (rigid) models parameterized to reproduce thermodynamic properties of water, have provided support to this hypothesis. A stronger numerical proof requires demonstrating that the critical point, which occurs at temperatures and pressures far from those at which the models were optimized, is robust with respect to model parameterization, specifically with respect to incorporating additional physical effects. Here, we show that a liquid–liquid critical point can be rigorously located also in the WAIL model of water [Pinnick et al., J. Chem. Phys. 137, 014510 (2012)], a model parameterized using ab initio calculations only. The model incorporates two features not present in many previously studied water models: It is both flexible and polarizable, properties which can significantly influence the phase behavior of water. The observation of the critical point in a model in which the water–water interaction is estimated using only quantum ab initio calculations provides strong support to the viewpoint according to which the existence of two distinct liquids is a robust feature in the free energy landscape of supercooled water. 
    more » « less
  3. Abstract Recent experiments and theoretical studies strongly indicate that water exhibits a liquid-liquid phase transition (LLPT) in the supercooled domain. An open question is how the LLPT of water can affect the properties of aqueous solutions. Here, we study the structural and thermodynamic properties of supercooled glycerol-water microdroplets at dilute conditions (χg = 3.2% glycerol mole fraction). The combination of rapid evaporative cooling with femtosecond X-ray scattering allows us to outrun crystallization and gain access to the deeply supercooled regime down toT = 229.3 K. We find that the density fluctuations of the glycerol-water solution or, equivalently, its isothermal compressibility,κT, increases upon cooling. This is confirmed by molecular dynamics simulations, which indicate that the presence of glycerol shifts the temperature of maximumκTfromT = 230 K in pure water down toT = 223 K in the solution. Our findings elucidate the interplay between the complex behavior of water, including its LLPT, and the properties of aqueous solutions at low temperatures, which can have practical consequences in cryogenic biological applications and cryopreservation techniques. 
    more » « less
  4. 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
  5. We investigate the microscopic pathway of spontaneous crystallization in the ST2 model of water under deeply supercooled conditions via unbiased classical molecular dynamics simulations. After quenching below the liquid–liquid critical point, the ST2 model spontaneously separates into low-density liquid (LDL) and high-density liquid phases, respectively. The LDL phase, which is characterized by lower molecular mobility and enhanced structural order, fosters the formation of a sub-critical ice nucleus that, after a stabilization time, develops into the critical nucleus and grows. Polymorphic selection coincides with the development of the sub-critical nucleus and favors the formation of cubic (Ic) over hexagonal (Ih) ice. We rationalize polymorphic selection in terms of geometric arguments based on differences in the symmetry of second neighbor shells of ice Ic and Ih, which are posited to favor formation of the former. The rapidly growing critical nucleus absorbs both Ic and Ih crystallites dispersed in the liquid phase, a crystal with stacking faults. Our results are consistent with, and expand upon, recent observations of non-classical nucleation pathways in several systems. 
    more » « less