Manifestations of metastable criticality in the long-range structure of model water glasses

Liquid-liquid transition and ice crystallization in a machine-learned coarse-grained water model

Mounting experimental evidence supports the existence of a liquid-liquid transition (LLT) in high-pressure supercooled water. However, fast crystallization of supercooled water has impeded identification of the LLT line TLL(p) in experiments. While the most accurate all-atom (AA) water models display a LLT, their computational cost limits investigations of its interplay with ice formation. Coarse-grained (CG) models provide over 100-fold computational efficiency gain over AA models, enabling the study of water crystallization, but have not yet shown to have a LLT. Here, we demonstrate that the CG machine-learned water model Machine-Learned Bond-Order Potential (ML-BOP) has a LLT that ends in a critical point at pc = 170 ± 10 MPa and Tc = 181 ± 3 K. The TLL(p) of ML-BOP is almost identical to the one of TIP4P/2005, adding to the similarity in the equation of state of liquid water in both models. Cooling simulations reveal that ice crystallization is fastest at the LLT and its supercritical continuation of maximum heat capacity, supporting a mechanistic relationship between the structural transformation of water to a low-density liquid (LDL) and ice formation. We find no signature of liquid-liquid criticality in the ice crystallization temperatures. ML-BOP replicates the competition between formation of LDL and ice observed in ultrafast experiments of decompression of the high-density liquid (HDL) into the region of stability of LDL. The simulations reveal that crystallization occurs prior to the coarsening of the HDL and LDL domains, obscuring the distinction between the highly metastable first-order LLT and pronounced structural fluctuations along its supercritical continuation.

Free-energy landscape and spinodals for the liquid-liquid transition of the TIP4P/2005 and TIP4P/Ice models of water

Continued increases in computational power now make it possible to evaluate the free-energy landscape associated with the first-order liquid-liquid transition in realistic models of water for which an accurate estimate of the liquid-liquid critical point exists, and to explore its change with pressure near the coexistence line. We report the results of 50 μs-long NPT umbrella sampling simulations for two realistic models for water, TIP4P/2005 and TIP4P/ice, 3-9 K below their critical temperatures. The free energy profile at different pressures clearly shows the presence of two well-defined free energy basins and makes it possible to identify the liquid-liquid spinodal points, the limits of stability that define the (temperature dependent) pressure range within which two distinct free energy basins exist. The results show that for temperatures less than 10 K below the critical temperature, metastable states are possible across a very limited pressure interval, information that is relevant to the interpretation of experiments probing the metastable phase behavior of deeply supercooled water in the so-called no-man's land.

High-dimensional order parameters and neural network classifiers applied to amorphous ices

Amorphous ice phases are key constituents of water's complex structural landscape. This study investigates the polyamorphic nature of water, focusing on the complexities within low-density amorphous ice (LDA), high-density amorphous ice, and the recently discovered medium-density amorphous ice (MDA). We use rotationally invariant, high-dimensional order parameters to capture a wide spectrum of local symmetries for the characterization of local oxygen environments. We train a neural network to classify these local environments and investigate the distinctiveness of MDA within the structural landscape of amorphous ice. Our results highlight the difficulty in accurately differentiating MDA from LDA due to structural similarities. Beyond water, our methodology can be applied to investigate the structural properties and phases of disordered materials.

A continuum of amorphous ices between low-density and high-density amorphous ice

Amorphous ices are usually classified as belonging to low-density or high-density amorphous ice (LDA and HDA) with densities ρLDA ≈ 0.94 g/cm3 and ρHDA ≈ 1.15-1.17 g/cm3. However, a recent experiment crushing hexagonal ice (ball-milling) produced a medium-density amorphous ice (MDA, ρMDA ≈ 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 at P ≈ 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 at P = 125 MPa is located in the intermediate region of the PEL that separates LDA and HDA.

Complex phase diagram and supercritical matter

The supercritical region is often described as uniform with no definite transitions. The distinct behaviors of the matter therein, e.g., as liquidlike and gaslike, however, suggest "supercritical boundaries." Here we provide a mathematical description of these phenomena by revisiting the Yang-Lee theory and introducing a complex phase diagram, specifically a four-dimensional (4D) one with complex T and p. While the traditional 2D phase diagram with real temperature T and pressure p values (the physical plane) lacks Lee-Yang (LY) zeros beyond the critical point, preventing the occurrence of criticality, the off-plane zeros in this 4D scenario still induce critical anomalies in various physical properties. This relationship is evidenced by the correlation between the Widom line and LY edges in van der Waals, 2D Ising model, and water. The diverged supercritical boundaries manifest the high-dimensional feature of the phase diagram: e.g., when LY zeros of complex T or p are projected onto the physical plane, boundaries defined by isobaric heat capacity C_{p} or isothermal compression coefficient K_{T} emanates. These results demonstrate the incipient phase transition nature of the supercritical matter.

High-density liquid (HDL) adsorption at the supercooled water/vapor interface and its possible relation to the second surface tension inflection point

We investigate the properties of water along the liquid/vapor coexistence line in the supercooled regime down to the no-man's land. Extensive molecular dynamics simulations of the TIP4P/2005 liquid/vapor interface in the range 198-348 K allow us to locate the second surface tension inflection point with a high accuracy at 283 ± 5 K, close to the temperature of maximum density. This temperature also coincides with the appearance of a density anomaly at the interface known as the apophysis. We relate the emergence of the apophysis to the observation of high-density liquid (HDL) water adsorption in the proximity of the liquid/vapor interface.

Manifestations of the structural origin of supercooled water’s anomalies in the heterogeneous relaxation on the potential energy landscape

Liquid water is well-known for its intriguing thermodynamic anomalies in the supercooled state. The phenomenological two-state models—based on the assumption of the existence of two types of competing local states (or, structures) in liquid water—have been extremely successful in describing water’s thermodynamic anomalies. However, the precise structural features of these competing local states in liquid water still remain elusive. Here, we have employed a predefined structural order parameter-free approach to unambiguously identify two types of competing local states—entropically and energetically favored—with significantly different structural and energetic features in the TIP4P/2005 liquid water. This identification is based on the heterogeneous structural relaxation of the system in the potential energy landscape (PEL) during the steepest-descent energy minimization. This heterogeneous relaxation is characterized using order parameters inspired by the spin-glass transition in frustrated magnetic systems. We have further established a direct relationship between the population fluctuation of the two states and the anomalous behavior of the heat capacity in supercooled water. The composition-dependent spatial distribution of the entropically favored local states shows an interesting crossover from a spanning network-like single cluster to the spatially delocalized clusters in the close vicinity of the Widom line. Additionally, this study establishes a direct relationship between the topographic features of the PEL and the water’s thermodynamic anomalies in the supercooled state and provides alternate markers (in addition to the locus of maxima of thermodynamic response functions) for the Widom line in the phase plane.

Steady-like topology of the dynamical hydrogen bond network in supercooled water

We investigate the link between topology of the hydrogen bond network (HBN) and large-scale density fluctuations in water from ambient conditions to the glassy state. We observe a transition from a temperature-dependent topology at high temperatures, to a steady-like topology below the Widom temperature T_W ∼ 220 K signaling the fragile-to-strong crossover and the maximum in structural fluctuations. As a consequence of the steady topology, the network suppresses large-scale density fluctuations much more efficiently than at higher temperatures. Below T_W , the contribution of coordination defects of the kind A ₂ D ₁ (two acceptors and one donor) to the kinetics of the HBN becomes progressively more pronounced, suggesting that A ₂ D ₁ configurations may represent the main source of dynamical heterogeneities. Below the vitrification temperature, the freezing of rotational and translational degrees of freedom allow for an enhanced suppression of large-scale density fluctuations and the sample reaches the edges of nearly hyperuniformity. The formed network still hosts coordination defects, hence implying that nearly hyperuniformity goes beyond the classical continuous random network paradigm of tetrahedral networks and can emerge in scenarios much more complex than previously assumed. Our results unveil a hitherto undisclosed link between network topology and properties of water essential for better understanding water's rich and complex nature. Beyond implications for water, our findings pave the way to a better understanding of the physics of supercooled liquids and disordered hyperuniform networks at large.

Evidence of a liquid-liquid phase transition in H[Formula: see text]O and D[Formula: see text]O from path-integral molecular dynamics simulations

We perform path-integral molecular dynamics (PIMD), ring-polymer MD (RPMD), and classical MD simulations of H[Formula: see text]O and D[Formula: see text]O using the q-TIP4P/F water model over a wide range of temperatures and pressures. The density [Formula: see text], isothermal compressibility [Formula: see text], and self-diffusion coefficients D(T) of H[Formula: see text]O and D[Formula: see text]O are in excellent agreement with available experimental data; the isobaric heat capacity [Formula: see text] 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[Formula: see text]O and D[Formula: see text]O exhibit a liquid-liquid critical point (LLCP) at low temperatures and positive pressures. The data from PIMD/MD for H[Formula: see text]O and D[Formula: see text]O 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[Formula: see text]O, from PIMD simulations, is located at [Formula: see text] MPa, [Formula: see text] K, and [Formula: see text] g/cm[Formula: see text]. Isotope substitution effects are important; the LLCP location in q-TIP4P/F D[Formula: see text]O is estimated to be [Formula: see text] MPa, [Formula: see text] K, and [Formula: see text] g/cm[Formula: see text]. 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, [Formula: see text] MPa, [Formula: see text] K, and [Formula: see text] g/cm[Formula: see text]). 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 [Formula: see text] for D[Formula: see text]O and, particularly, H[Formula: see text]O suggest that improved water models are needed for the study of supercooled water.

Molecular Simulations of Aqueous Electrolytes: Role of Explicit Inclusion of Charge Transfer into Force Fields

We describe here simulations of aqueous salt solutions that are performed using an explicit charge transfer force field. The emphasis of the discussion is on the calculation of a dynamical property of the solutions: self-diffusion of water. While force fields that are based on pairwise additive potentials or on potentials with explicit inclusion of polarization or with scaled charges can provide at best a qualitative agreement with experiments, force fields with explicit inclusion of charge transfer can produce quantitative agreement with experiment for NaCl and KCl solutions. We argue that a force field with explicit charge transfer contains new physics absent in the previously used force fields described in recent reviews of molecular simulations of aqueous electrolytes.

Percolation transitions in compressed SiO2 glasses

Amorphous-amorphous transformations under pressure are generally explained by changes in the local structure from low- to higher-fold coordinated polyhedra^1-4. However, as the notion of scale invariance at the critical thresholds has not been addressed, it is still unclear whether these transformations behave similarly to true phase transitions in related crystals and liquids. Here we report ab initio-based calculations of compressed silica (SiO₂) glasses, showing that the structural changes from low- to high-density amorphous structures occur through a sequence of percolation transitions. When the pressure is increased to 82 GPa, a series of long-range ('infinite') percolating clusters composed of corner- or edge-shared tetrahedra, pentahedra and eventually octahedra emerge at critical pressures and replace the previous 'phase' of lower-fold coordinated polyhedra and lower connectivity. This mechanism provides a natural explanation for the well-known mechanical anomaly around 3 GPa, as well as the structural irreversibility beyond 10 GPa, among other features. Some of the amorphous structures that have been discovered mimic those of coesite IV and V crystals reported recently^5,6, highlighting the major role of SiO₅ pentahedron-based polyamorphs in the densification process of vitreous silica. Our results demonstrate that percolation theory provides a robust framework to understand the nature and pathway of amorphous-amorphous transformations and open a new avenue to predict unravelled amorphous solid states and related liquid phases^7,8.

The role of high-density and low-density amorphous ice on biomolecules at cryogenic temperatures: a case study with polyalanine

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).