Structural and topological changes across the liquid-liquid transition in water

Phase behavior of metastable water from large-scale simulations of a quantitatively accurate model near ambient conditions: The liquid-liquid critical point

The molecular mechanisms of water's unique anomalies are still debated upon. Experimental challenges have led to simulations suggesting a liquid-liquid (LL) phase transition, culminating in the supercooled region's LL critical point (LLCP). Computational expense, small system sizes, and the reliability of water models often limit these simulations. We adopt the CVF model, which is reliable, transferable, scalable, and efficient across a wide range of temperatures and pressures around ambient conditions. By leveraging the timescale separation between fast hydrogen bonds and slow molecular coordinates, the model allows a thorough exploration of the metastable phase diagram of liquid water. Using advanced numerical techniques to bypass dynamical slowing down, we perform finite-size scaling on larger systems than those used in previous analyses. Our study extrapolates thermodynamic behavior in the infinite-system limit, demonstrating the existence of the LLCP in the 3D Ising universality class in the low-temperature, low-pressure side of the line of temperatures of maximum density, specifically at TC = 186 ± 4 K and PC = 174 ± 14 MPa, at the end of a liquid-liquid phase separation stretching up to ∼200 MPa. These predictions align with recent experimental data and sophisticated models, highlighting that hydrogen bond cooperativity governs the LLCP and the origin of water anomalies. We also observe substantial cooperative fluctuations in the hydrogen bond network at scales larger than 10 nm, even at temperatures relevant to biopreservation. These findings have significant implications for nanotechnology and biophysics, providing new insights into water's behavior under varied conditions.

Hierarchy of topological transitions in a network liquid

The representation of complex systems as networks has become a critical tool across many fields of science. In the context of physical networks, such as biological neural networks, vascular networks, or network liquids where the nodes and edges occupy volume in three-dimensional space, the question of how they become densely packed is of special importance. Here, we investigate a model network liquid, which is known to densify via two successive liquid-liquid phase transitions (LLPTs). We elucidate the importance of rings-cyclic paths formed by bonded particles in the networks-and their spatial disposition in understanding the structural changes that underpin the increase in density across the LLPTs. Our analyses demonstrate that the densification of these networks is primarily driven by the formation of linked rings, and the LLPTs correspond to a hierarchy of topological transitions where rings form the fundamental building blocks. We envisage entanglement to emerge as a general mechanism for densification, with wide implications for the embedding of physical networks, especially in confined spaces.

Water Potential from Adaptive Force Matching for Ice and Liquid with Revised Dispersion Predicts Supercooled Liquid Anomalies in Good Agreement with Two Independent Experimental Fits

A revised version of the Water potential from Adaptive force matching for Ice and Liquid (WAIL) was developed by using the previous data set for fitting the WAIL model but with a dispersion term calculated using symmetry adapted perturbation theory (SAPT). The model has no adjustable parameters and relies solely on fitting first-principles information. The new model, named revised WAIL (rWAIL), shows improved predictions of most properties of water when compared to the previously published WAIL model. The rWAIL model also compares favorably to other first-principles-derived water models, such as MB-Pol, at only a fraction of the computational cost. The rWAIL model is used to study the properties of supercooled water. The model shows evidence of a liquid-liquid phase transition (LLPT) in the supercooled regimes with the liquid-liquid critical point (LLCP) at 203 K and 90 MPa. This estimate is in good agreement with a recent polynomial fit to the experimental density of water. Also, the fit to the surface tension of supercooled water based on the rWAIL model shows excellent agreement with the corresponding fit to the experimental data. Consistent with previously published molecular dynamics and experimental data, the surface tension of water exhibits exponential growth in the supercooled regime, which is likely a result of the emergence of a low-density liquid form of water. The simulation thus unites two separate experimental fits with one first-principles-based model, lending strong evidence of an LLPT in real water.

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.

Identification of local structures in water from supercooled to ambient conditions

Studies of water thermodynamics have long been tied to the identification of two distinct families of local structures, whose competition could explain the origin of the many thermodynamic anomalies and the hypothesized liquid-liquid critical point in water. Despite the many successes and insights gained, the structural indicators proposed throughout the years were not able to unequivocally identify these two families over a wide range of conditions. We show that a recently introduced indicator, Ψ, which exploits information on the hydrogen bond network connectivity, can reliably identify these two distinct local environments over a wide range of thermodynamic conditions (188-300 K and 0-13 kbar) and that close to the liquid-liquid critical point, the spatial correlations of density fluctuations are identical to those of the Ψ indicator. Our results strongly support the idea that water thermodynamic properties arise from the competition between two distinct and identifiable local environments.

A neural network potential with self-trained atomic fingerprints: A test with the mW water potential

We present a neural network (NN) potential based on a new set of atomic fingerprints built upon two- and three-body contributions that probe distances and local orientational order, respectively. Compared with the existing NN potentials, the atomic fingerprints depend on a small set of tunable parameters that are trained together with the NN weights. In addition to simplifying the selection of the atomic fingerprints, this strategy can also considerably increase the overall accuracy of the network representation. To tackle the simultaneous training of the atomic fingerprint parameters and NN weights, we adopt an annealing protocol that progressively cycles the learning rate, significantly improving the accuracy of the NN potential. We test the performance of the network potential against the mW model of water, which is a classical three-body potential that well captures the anomalies of the liquid phase. Trained on just three state points, the NN potential is able to reproduce the mW model in a very wide range of densities and temperatures, from negative pressures to several GPa, capturing the transition from an open random tetrahedral network to a dense interpenetrated network. The NN potential also reproduces very well properties for which it was not explicitly trained, such as dynamical properties and the structure of the stable crystalline phases of mW.

Correlated Fluctuations of Structural Indicators Close to the Liquid-Liquid Transition in Supercooled Water

Multiple numerical studies have unambiguously shown the existence of a liquid-liquid critical point in supercooled states for different numerical models of water, and various structural indicators have been put forward to describe the transformation associated with this phase transition. Here we analyze numerical simulations of near-critical supercooled water to compare the behavior of several of such indicators with critical density fluctuations. We show that close to the critical point most indicators are strongly correlated to density, and some of them even display identical distributions of fluctuations. These indicators probe the exact same free energy landscape, therefore providing a thermodynamic description of critical supercooled water which is identical to that provided by the density order parameter. This implies that close to the critical point, there is a tight coupling between many, only apparently distinct, structural degrees of freedom.

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.

Temperature- and pressure-dependence of the hydrogen bond network in plastic ice VII

We model, via classical molecular dynamics simulations, the plastic phase of ice VII across a wide range of the phase diagram of interest for planetary investigations. Although structural and dynamical properties of plastic ice VII are mostly independent on the thermodynamic conditions, the hydrogen bond network (HBN) acquires a diverse spectrum of topologies distinctly different from that of liquid water and of ice VII simulated at the same pressure. We observe that the HBN topology of plastic ice carries some degree of similarity with the crystal phase, stronger at thermodynamic conditions proximal to ice VII, and gradually lessening upon approaching the liquid state. Our results enrich our understanding of the properties of water at high pressure and high temperature, and may help in rationalizing the geology of

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 TW ∼ 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 TW , the contribution of coordination defects of the kind A 2 D 1 (two acceptors and one donor) to the kinetics of the HBN becomes progressively more pronounced, suggesting that A 2 D 1 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.

Liquid-Liquid Criticality in the WAIL Water Model

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 parametrized 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 [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.

Anomalies and Local Structure of Liquid Water from Boiling to the Supercooled Regime as Predicted by the Many-Body MB-pol Model

For the past 50 years, researchers have sought molecular models that can accurately reproduce water's microscopic structure and thermophysical properties across broad ranges of its complex phase diagram. Herein, molecular dynamics simulations with the many-body MB-pol model are performed to monitor the thermodynamic response functions and local structure of liquid water from the boiling point down to deeply supercooled temperatures at ambient pressure. The isothermal compressibility and isobaric heat capacity show maxima near 223 K, in excellent agreement with recent experiments, and the liquid density exhibits a minimum at ∼208 K. A local tetrahedral arrangement, where each water molecule accepts and donates two hydrogen bonds, is found to be the most probable hydrogen-bonding topology at all temperatures. This work suggests that MB-pol may provide predictive capability for studies of liquid water's physical properties across broad ranges of thermodynamic states, including the so-called water's "no man's land" which is difficult to probe experimentally.

The physics of empty liquids: from patchy particles to water

Empty liquids represent a wide class of materials whose constituents arrange in a random network through reversible bonds. Many key insights on the physical properties of empty liquids have originated almost independently from the study of colloidal patchy particles on one side, and a large body of theoretical and experimental research on water on the other side. Patchy particles represent a family of coarse-grained potentials that allows for a precise control of both the geometric and the energetic aspects of bonding, while water has arguably the most complex phase diagram of any pure substance, and a puzzling amorphous phase behavior. It was only recently that the exchange of ideas from both fields has made it possible to solve long-standing problems and shed new light on the behavior of empty liquids. Here we highlight the connections between patchy particles and water, focusing on the modelling principles that make an empty liquid behave like water, including the factors that control the appearance of thermodynamic and dynamic anomalies, the possibility of liquid-liquid phase transitions, and the crystallization of open crystalline structures.

Advances in the study of supercooled water

In this review, we report recent progress in the field of supercooled water. Due to its uniqueness, water presents numerous anomalies with respect to most simple liquids, showing polyamorphism both in the liquid and in the glassy state. We first describe the thermodynamic scenarios hypothesized for the supercooled region and in particular among them the liquid-liquid critical point scenario that has so far received more experimental evidence. We then review the most recent structural indicators, the two-state model picture of water, and the importance of cooperative effects related to the fact that water is a hydrogen-bonded network liquid. We show throughout the review that water's peculiar properties come into play also when water is in solution, confined, and close to biological molecules. Concerning dynamics, upon mild supercooling water behaves as a fragile glass former following the mode coupling theory, and it turns into a strong glass former upon further cooling. Connections between the slow dynamics and the thermodynamics are discussed. The translational relaxation times of density fluctuations show in fact the fragile-to-strong crossover connected to the thermodynamics arising from the existence of two liquids. When considering also rotations, additional crossovers come to play. Mobility-viscosity decoupling is also discussed in supercooled water and aqueous solutions. Finally, the polyamorphism of glassy water is considered through experimental and simulation results both in bulk and in salty aqueous solutions. Grains and grain boundaries are also discussed.

Structure of High-Pressure Supercooled and Glassy Water

We numerically investigate the structure of deep supercooled and glassy water under pressure, covering the range of densities corresponding to the experimentally produced high- and very-high-density amorphous phases. At T=188  K, a continuous increase in density is observed on varying pressure from 2.5 to 13 kbar, with no signs of first-order transitions. Exploiting a recently proposed approach to the analysis of the radial distribution function-based on topological properties of the hydrogen-bond network-we are able to identify well-defined local geometries that involve pairs of molecules separated by multiple hydrogen bonds, specific to the high- and very-high-density structures.