No Arabic abstract
Recent computational studies have reported evidence of a metastable liquid-liquid phase transition (LLPT) in molecular models of water under deeply supercooled conditions. A competing hypothesis suggests, however, that non-equilibrium artifacts associated with coarsening of the stable crystal phase have been mistaken for an LLPT in these models. Such artifacts are posited to arise due to a separation of time scales in which density fluctuations in the supercooled liquid relax orders of magnitude faster than those associated with bond-orientational order. Here, we use molecular simulation to investigate the relaxation of density and bond-orientational fluctuations in three molecular models of water (ST2, TIP5P and TIP4P/2005) in the vicinity of their reported LLPT. For each model, we find that density is the slowly relaxing variable under such conditions. We also observe similar behavior in the coarse-grained mW model of water. Our findings therefore challenge the key physical assumption underlying the competing hypothesis.e find that density relaxes significantly faster than bond-orientational order, as incorrectly predicted by this competing hypothesis.
We compare the spatial correlations of bond-breaking events and bond-orientational relaxation in a model two-dimensional liquid undergoing Newtonian dynamics. We find that the relaxation time of the bond-breaking correlation function is much longer than the relaxation time of the bond-orientational correlation function and self-intermediate scattering function. However, the relaxation time of the bond-orientational correlation function increases faster with decreasing temperature than the relaxation time of the bond-breaking correlation function and the self-intermediate scattering function. Moreover, the dynamic correlation length that characterizes the size of correlated bond-orientational relaxation grows faster with decreasing temperature than the dynamic correlation length that characterizes the size of correlated bond-breaking events. We also examine the ensemble-dependent and ensemble-independent dynamic susceptibilities for both bond-breaking correlations and bond-orientational correlations. We find that for both correlations, the ensemble-dependent and ensemble-independent susceptibilities exhibit a maximum at nearly the same time, and this maximum occurs at a time slightly shorter than the peak position of the dynamic correlation length.
We use the recently-proposed emph{compressible cell} Ising-like model [Phys. Rev. Lett. textbf{120}, 120603 (2018)] to estimate the ratio between thermal expansivity and specific heat (the Gruneisen parameter $Gamma$) in supercooled water. Near the critical pressure and temperature, $Gamma$ increases. The $Gamma$ value diverges near the pressure-induced finite-$T$ critical end-point [Phys. Rev. Lett. textbf{104}, 245701 (2010)] and quantum critical points [Phys. Rev. Lett. textbf{91}, 066404 (2003)], which indicates that two energy scales are governing the system. This enhanced behavior of $Gamma$ is caused by the coexistence of high- and low-density liquids [Science textbf{358}, 1543 (2017)]. Our findings support the proposed liquid-liquid critical point in supercooled water in the No-Mans Land regime, and indicates possible applications of this model to other systems.
Deeply supercooled water exhibits complex dynamics with large density fluctuations, ice coarsening and characteristic time scales extending from picoseconds to milliseconds. Here, we discuss implications of these time scales as they pertain to two-phase coexistence and to molecular simulations of supercooled water. Specifically, we argue that it is possible to discount liquid-liquid criticality because the time scales imply that correlation lengths for such behavior would be bounded by no more than a few nanometers. Similarly, it is possible to discount two-liquid coexistence because the time scales imply a bounded interfacial free energy that cannot grow in proportion to a macroscopic surface area. From time scales alone, therefore, we see that coexisting domains of differing density in supercooled water can be no more than nano-scale transient fluctuations.
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.
Water shows intriguing thermodynamic and dynamic anomalies in the supercooled liquid state. One possible explanation of the origin of these anomalies lies in the existence of a metastable liquid-liquid phase transition (LLPT) between two (high and low density) forms of water. While the anomalies are observed in experiments on bulk and confined water and by computer simulation studies of water-like models, the existence of a LLPT in water is still debated. Unambiguous experimental proof of the existence of a LLPT in bulk supercooled water is hampered by fast ice nucleation which is a precursor of the hypothesized LLPT. Moreover, the hypothesized LLPT, being metastable, in principle cannot exist in the thermodynamic limit (infinite size, infinite time). Therefore, computer simulations of water models are crucial for exploring the possibility of the metastable LLPT and the nature of the anomalies. In this work, we present new simulation results in the NVT ensemble for one of the most accurate classical molecular models of water, TIP4P/2005. To describe the computed properties and explore the possibility of a LLPT we have applied two-structure thermodynamics, viewing water as a non-ideal mixture of two interconvertible local structures (states). The results suggest the presence of a liquid-liquid critical point and a LLPT in this model for the simulated length and time scales. We have compared the behavior of TIP4P/2005 with other popular water-like models, namely mW and ST2, and with real water, all of which are well described by two-state thermodynamics. In view of the current debate involving different studies of TIP4P/2005, we discuss consequences of metastability and finite size in observing the liquid-liquid separation. We also address the relationship between the phenomenological order parameter of two-structure thermodynamics and the microscopic nature of the low-density structure.