Finite-size scaling investigation of the liquid-liquid critical point in ST2 water and its stability with respect to crystallization

Crystallization and the liquid-liquid critical point in nonbonded modified-WAC models

For decades, it has been known that Liquid-Liquid Critical Points (LLCPs) can exist in one-component liquids, yet a comprehensive understanding of the conditions under which they arise remains elusive. To better comprehend the possible interplay between the LLCP and the crystalline phase, we conduct molecular dynamics simulations using the nonbonded family of modified-WAC (mWAC) models, which are known to exhibit a LLCP for certain parameter values. By comparing different versions of the mWAC model-those featuring a LLCP and those lacking one-we identify several key differences between the models relating to crystallization. Those models that do have a LLCP are found to have multiple stable crystalline phases, one of them being a solid-state ionic conductor similar to superionic ice. Moreover, we find that for models that do not have a LLCP, the liquid becomes a glass at a larger range of temperatures, possibly preventing the occurrence of a LLCP. Further studies are required to determine if these results are general or model-specific.

Unraveling thermodynamic anomalies of water: A molecular simulation approach to probe the two-state theory with atomistic and coarse-grained water models

Thermodynamic and dynamic anomalies of water play a crucial role in supporting life on our planet. The two-state theory attributes these anomalies to a dynamic equilibrium between locally favored tetrahedral structures (LFTSs) and disordered normal liquid structures. This theory provides a straightforward, phenomenological explanation for water's unique thermodynamic and dynamic characteristics. To validate this two-state feature, it is critical to unequivocally identify these structural motifs in a dynamically fluctuating disordered liquid. In this study, we employ a recently introduced structural parameter (θavg) that characterizes the local angular order within the first coordination shell to identify these LFTSs through molecular dynamics simulations. We employ both realistic water models with a liquid-liquid critical point (LLCP) and a coarse-grained water model without an LLCP to study water's anomalies in low-pressure regions below 2 kbar. The two-state theory consistently describes water's thermodynamic anomalies in these models, both with and without an LLCP. This suggests that the anomalies predominantly result from the two-state features rather than criticality, particularly within experimentally accessible temperature-pressure regions.

Supercooled solutions of sodium perchlorate in TIP4P/2005 water: The effect of martian solutes on thermodynamics and structure

We study the thermodynamic behavior of sodium perchlorate solutions in supercooled water through molecular dynamics numerical simulations. These solutions are of special interest because of the recent experimental results that led to hypothesize the presence of liquid water in perchlorate solutions beneath the Martian soil. We model water using the TIP4P/2005 potential. The results we obtain for solutions with concentrations 1.63 and 15.4 wt% are in agreement with those of a system undergoing a liquid-liquid phase transition where the liquid-liquid critical point shifts to slightly higher temperatures and lower pressures. The structure of the system is also analyzed, and we come to the conclusion that, even at the highest concentration considered, water retains its anomalous behavior.

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.

Thermodynamics and kinetics of crystallization in deeply supercooled Stillinger-Weber silicon

We study the kinetics of crystallization in deeply supercooled liquid silicon employing computer simulations and the Stillinger-Weber three-body potential. The free energy barriers to crystallization are computed using umbrella sampling Monte Carlo simulations and from unconstrained molecular dynamics simulations using a mean first passage time formulation. We focus on state points that have been described in earlier work [S. Sastry and C. A. Angell, Nat. Mater. 2, 739 (2003)] as straddling a liquid-liquid phase transition (LLPT) between two metastable liquid states. It was argued subsequently [Ricci et al., Mol. Phys. 117, 3254 (2019)] that the apparent transition is due to the loss of metastability of the liquid state with respect to the crystalline state. The presence of a barrier to crystallization for these state points is therefore of importance to ascertain, which we investigate, with due attention to ambiguities that may arise from the choice of order parameters. We find a well-defined free energy barrier to crystallization and demonstrate that both umbrella sampling and mean first passage time methods yield results that agree quantitatively. Our results thus provide strong evidence against the possibility that the liquids at state points close to the reported LLPT exhibit slow, spontaneous crystallization, but they do not address the existence of a LLPT (or lack thereof). We also compute the free energy barriers to crystallization at other state points over a broad range of temperatures and pressures and discuss the effect of changes in the microscopic structure of the metastable liquid on the free energy barrier heights.

A proposal for the structure of high- and low-density fluctuations in liquid water

Based on recent experimental data that can be interpreted as indicating the presence of specific structures in liquid water, we build and optimize two structural models which we compare with the available experimental data. To represent the proposed high-density liquid structures, we use a model consisting of chains of water molecules, and for low-density liquid, we investigate fused dodecahedra as templates for tetrahedral fluctuations. The computed infrared spectra of the models are in very good agreement with the extracted experimental spectra for the two components, while the extracted structures from molecular dynamics (MD) simulations give spectra that are intermediate between the experimentally derived spectra. Computed x-ray absorption and emission spectra as well as the O-O radial distribution functions of the proposed structures are not contradicted by experiment. The stability of the proposed dodecahedral template structures is investigated in MD simulations by seeding the starting structure, and remnants found to persist on an ∼30 ps time scale. We discuss the possible significance of such seeds in simulations and whether they can be viable candidates as templates for structural fluctuations below the compressibility minimum of liquid water.

Fluctuations near the liquid–liquid transition in a model of silica

Molecular dynamics simulations reveal anomalous small-angle scattering and liquid–liquid phase separation in an ionic model of silica.

The phase behavior study of human antibody solution using multi-scale modeling

Phase transformation in antibody solutions is of growing interest in both academia and the pharmaceutical industry. Recent experimental studies have shown that, as in near-spherical proteins, antibodies can undergo a liquid-liquid phase separation under conditions metastable with respect to crystallization. However, the phase diagram of the Y-shaped antibodies exhibits unique features that differ substantially from those of spherical proteins. Specifically, antibody solutions have an exceptionally low critical volume fraction (CVF) and a broader and more asymmetric liquid-liquid coexistence curve than those of spherical proteins. Using molecular dynamics simulation on a series of trimetric Y-shaped coarse-grained models, we investigate the phase behavior of antibody solutions and compare the results with the experimental phase diagram of human immunoglobulin G (IgG), one of the most common Y-shape typical of antibody molecules. With the fitted size of spheres, our simulation reproduces both the low CVF and the asymmetric shape of the experimental coexistence curve of IgG antibodies. The broadness of the coexistence curve can be attributed to the anisotropic nature of the inter-protein interaction. In addition, the repulsion between the inner parts of the spherical domains of IgG dramatically expands the coexistence region in the scaled phase diagram, while the hinge length has only a minor effect on the CVF and the overall shape of the coexistence curve. We thus propose a seven-site model with empirical parameters characterizing the exclusion volume and the hinge length of the IgG molecules, which provides a base for simulation studies of the phase behavior of IgG antibodies.

Failure and recovery in dynamical networks

Failure, damage spread and recovery crucially underlie many spatially embedded networked systems ranging from transportation structures to the human body. Here we study the interplay between spontaneous damage, induced failure and recovery in both embedded and non-embedded networks. In our model the network's components follow three realistic processes that capture these features: (i) spontaneous failure of a component independent of the neighborhood (internal failure), (ii) failure induced by failed neighboring nodes (external failure) and (iii) spontaneous recovery of a component. We identify a metastable domain in the global network phase diagram spanned by the model's control parameters where dramatic hysteresis effects and random switching between two coexisting states are observed. This dynamics depends on the characteristic link length of the embedded system. For the Euclidean lattice in particular, hysteresis and switching only occur in an extremely narrow region of the parameter space compared to random networks. We develop a unifying theory which links the dynamics of our model to contact processes. Our unifying framework may help to better understand controllability in spatially embedded and random networks where spontaneous recovery of components can mitigate spontaneous failure and damage spread in dynamical networks.

Evidence for a liquid-liquid critical point in supercooled water within the E3B3 model and a possible interpretation of the kink in the homogeneous nucleation line

Supercooled water exhibits many thermodynamic anomalies, and several scenarios have been proposed to interpret them, among which the liquid-liquid critical point (LLCP) hypothesis is the most commonly discussed. We investigated Widom lines and the LLCP of deeply supercooled water, by using molecular dynamics simulation with a newly reparameterized water model that explicitly includes three-body interactions. Seven isobars are studied from ambient pressure to 2.5 kbar, and Widom lines are identified by calculating maxima in the coefficient of thermal expansion and the isothermal compressibility (both with respect to temperature). From these data we estimate that the LLCP of the new water model is at 180 K and 2.1 kbar. The oxygen radial distribution function is calculated along the 2 kbar isobar. It shows a steep change in the height of its second peak between 180 and 185 K, which indicates a transition between the high-density liquid and low-density liquid phases and which is consistent with the ascribed location of the critical point. The good agreement of the height of the second peak of the radial distribution function between simulation and experiment at 1 bar, as a function of temperature, supports the validity of the model. The location of the LLCP within the model is close to the kink in the experimental homogeneous nucleation line. We use existing experimental data to argue that the experimental LLCP is at 168 K and 1.95 kbar and speculate how this LLCP and its Widom line might be responsible for the kink in the homogeneous nucleation line.

Water: A Tale of Two Liquids

Water is the most abundant liquid on earth and also the substance with the largest number of anomalies in its properties. It is a prerequisite for life and as such a most important subject of current research in chemical physics and physical chemistry. In spite of its simplicity as a liquid, it has an enormously rich phase diagram where different types of ices, amorphous phases, and anomalies disclose a path that points to unique thermodynamics of its supercooled liquid state that still hides many unraveled secrets. In this review we describe the behavior of water in the regime from ambient conditions to the deeply supercooled region. The review describes simulations and experiments on this anomalous liquid. Several scenarios have been proposed to explain the anomalous properties that become strongly enhanced in the supercooled region. Among those, the second critical-point scenario has been investigated extensively, and at present most experimental evidence point to this scenario. Starting from very low temperatures, a coexistence line between a high-density amorphous phase and a low-density amorphous phase would continue in a coexistence line between a high-density and a low-density liquid phase terminating in a liquid–liquid critical point, LLCP. On approaching this LLCP from the one-phase region, a crossover in thermodynamics and dynamics can be found. This is discussed based on a picture of a temperature-dependent balance between a high-density liquid and a low-density liquid favored by, respectively, entropy and enthalpy, leading to a consistent picture of the thermodynamics of bulk water. Ice nucleation is also discussed, since this is what severely impedes experimental investigation of the vicinity of the proposed LLCP. Experimental investigation of stretched water, i.e., water at negative pressure, gives access to a different regime of the complex water diagram. Different ways to inhibit crystallization through confinement and aqueous solutions are discussed through results from experiments and simulations using the most sophisticated and advanced techniques. These findings represent tiles of a global picture that still needs to be completed. Some of the possible experimental lines of research that are essential to complete this picture are explored.

Tunable Liquid-Liquid Critical Point in an Ionic Model of Silica

Recently, it was shown that the Woodcock-Angell-Cheeseman model for liquid silica [L. V. Woodcock, C. A. Angell, and P. Cheeseman, J. Chem. Phys. 65, 1565 (1976)] is remarkably close to having a liquid-liquid critical point (LLCP). We demonstrate that increasing the ion charge separates the global maxima of the response functions, while reducing the charge smoothly merges them into a LLCP, a phenomenon that might be experimentally observable with charged colloids. An analysis of the Si and O coordination numbers suggests that a sufficiently low Si/O coordination number ratio is needed to attain a LLCP.

Pressure-induced transformations in glassy water: A computer simulation study using the TIP4P/2005 model

We study the pressure-induced transformations between low-density amorphous (LDA) and high-density amorphous (HDA) ice by performing out-of-equilibrium molecular dynamics (MD) simulations. We employ the TIP4P/2005 water model and show that this model reproduces qualitatively the LDA-HDA transformations observed experimentally. Specifically, the TIP4P/2005 model reproduces remarkably well the (i) structure (OO, OH, and HH radial distribution functions) and (ii) densities of LDA and HDA at P = 0.1 MPa and T = 80 K, as well as (iii) the qualitative behavior of ρ(P) during compression-induced LDA-to-HDA and decompression-induced HDA-to-LDA transformations. At the rates explored, the HDA-to-LDA transformation is less pronounced than in experiments. By studying the LDA-HDA transformations for a broad range of compression/decompression temperatures, we construct a "P-T phase diagram" for glassy water that is consistent with experiments and remarkably similar to that reported previously for ST2 water. This phase diagram is not inconsistent with the possibility of TIP4P/2005 water exhibiting a liquid-liquid phase transition at low temperatures. A comparison with previous MD simulation studies of SPC/E and ST2 water as well as experiments indicates that, overall, the TIP4P/2005 model performs better than the SPC/E and ST2 models. The effects of cooling and compression rates as well as aging on our MD simulations results are also discussed. The MD results are qualitatively robust under variations of cooling/compression rates (accessible in simulations) and are not affected by aging the hyperquenched glass for at least 1 μs. A byproduct of this work is the calculation of TIP4P/2005 water's diffusion coefficient D(T) at P = 0.1 MPa. It is found that, for T ≥ 210 K, D(T) ≈ (T - T(MCT))(-γ) as predicted by mode coupling theory and in agreement with experiments. For TIP4P/2005 water, T(MCT) = 209 K and γ = 2.14, very close to the corresponding experimental values T(MCT) = 221 K and γ = 2.2.

Hydrophobicity: effect of density and order on water's rotational slowing down

Ab initio molecular dynamics (AIMD) simulations of over 4.5 ns were performed in the temperature range of T = 260-350 K with van der Waals corrections to investigate the relationship between local water density and tetrahedral order in bulk water and in the presence of a hydrophobe, tetramethylurea (TMU). We demonstrate that in bulk water, defects consisting of 5- and higher coordinated water are a major contributor to dynamics. Close to a hydrophobe, 3-coordinated defects take over. The co-existence of these defects gives rise to very different local densities. We propose that the slowing down of rotational motion close to a hydrophobe is induced by an interplay between density and order with the slowing down decreasing in the following order: (i) low-density ordered-water, (ii) normal-density ordered-water, (iii) high-density ordered-water and (iv) disordered-water. The proportions of these water environments vary with temperature. These local environments also support the idea of water's polymorphism, i.e., the existence of the high- and low-density states in supercooled water.

Correlation between thermodynamic anomalies and pathways of ice nucleation in supercooled water

The well-known classical nucleation theory (CNT) for the free energy barrier towards formation of a nucleus of critical size of the new stable phase within the parent metastable phase fails to take into account the influence of other metastable phases having density/order intermediate between the parent metastable phase and the final stable phase. This lacuna can be more serious than capillary approximation or spherical shape assumption made in CNT. This issue is particularly significant in ice nucleation because liquid water shows rich phase diagram consisting of two (high and low density) liquid phases in supercooled state. The explanations of thermodynamic and dynamic anomalies of supercooled water often invoke the possible influence of a liquid-liquid transition between two metastable liquid phases. To investigate both the role of thermodynamic anomalies and presence of distinct metastable liquid phases in supercooled water on ice nucleation, we employ density functional theoretical approach to find nucleation free energy barrier in different regions of phase diagram. The theory makes a number of striking predictions, such as a dramatic lowering of nucleation barrier due to presence of a metastable intermediate phase and crossover in the dependence of free energy barrier on temperature near liquid-liquid critical point. These predictions can be tested by computer simulations as well as by controlled experiments.

Relation between the two-body entropy and the relaxation time in supercooled water

The two-body excess entropy of supercooled water is calculated from the radial distribution functions obtained from computer simulation of the TIP4P model for different densities upon supercooling. This quantity is considered in connection with the relaxation time of the self intermediate scattering function. The relaxation time shows a mode coupling theory (MCT) behavior in the region of mild supercooling and a strong behavior in the deep supercooled region. We find here that the two-body entropy is connected to the relaxation time and shows a logarithmic behavior with an apparent asymptotic divergence at the mode coupling crossover temperature. There is also evidence of a change in behavior of the two-body entropy upon crossing from the fragile (hopping-free) state to the strong (hopping-dominated) state of supercooled water, and the relation that connects the two-body entropy and the relxation time in the MCT region no longer holds.

Supercooled water escaping from metastability

The return of supercooled water to a stable equilibrium condition is an irreversible process which, in large enough samples, takes place adiabatically. We investigated this phenomenon in water by fast imaging techniques. As water freezes, large energy and density fluctuations promote the spatial coexistence of solid and liquid phases at different temperatures. Upon synchronously monitoring the time evolution of the local temperature, we observed a sharp dynamic transition between a fast and a slow decay regime at about 266.6 K. We construe the observed phenomenon in terms of the temperature dependence of heat transfers from solid and liquid volumes already at their bulk coexistence temperature towards adjacent still supercooled liquid regions. These findings can be justified by observing that convective motions induced by thermal gradients in a supercooled liquid near coexistence are rapidly suppressed as the nucleated solid fraction overcomes, at low enough temperatures, a characteristic percolation threshold.

Two-state thermodynamics of the ST2 model for supercooled water

Thermodynamic properties of the ST2 model for supercooled liquid water exhibit anomalies similar to those observed in real water. A possible explanation of these anomalies is the existence of a metastable, liquid-liquid transition terminated by a critical point. This phenomenon, whose possible existence in real water is the subject of much current experimental work, has been unambiguously demonstrated for this particular model by most recent simulations. In this work, we reproduce the anomalies of two versions of the ST2 model with an equation of state describing water as a non-ideal "mixture" of two different types of local molecular order. We show that the liquid-liquid transition in the ST2 water is energy-driven. This is in contrast to another popular model, mW, in which non-ideality in mixing of two alternative local molecular orders is entropy-driven, and is not sufficiently strong to induce a liquid-liquid transition.

Critical behavior of a water monolayer under hydrophobic confinement

The properties of water can have a strong dependence on the confinement. Here, we consider a water monolayer nanoconfined between hydrophobic parallel walls under conditions that prevent its crystallization. We investigate, by simulations of a many-body coarse-grained water model, how the properties of the liquid are affected by the confinement. We show, by studying the response functions and the correlation length and by performing finite-size scaling of the appropriate order parameter, that at low temperature the monolayer undergoes a liquid-liquid phase transition ending in a critical point in the universality class of the two-dimensional (2D) Ising model. Surprisingly, by reducing the linear size L of the walls, keeping the walls separation h constant, we find a 2D-3D crossover for the universality class of the liquid-liquid critical point for L=h ≃ 50, i.e. for a monolayer thickness that is small compared to its extension. This result is drastically different from what is reported for simple liquids, where the crossover occurs for L=h ≃ 5, and is consistent with experimental results and atomistic simulations. We shed light on these findings showing that they are a consequence of the strong cooperativity and the low coordination number of the hydrogen bond network that characterizes water.