Thermodynamic Analysis of the Two-Liquid Model for Anomalies of Water, HDL-LDL Fluctuations, and Liquid-Liquid Transition

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.

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.

On the existence of soliton-like collective modes in liquid water at the viscoelastic crossover

The problem of large-density variations in supercooled and ambient water has been widely discussed in the past years. Recent studies have indicated the possibility of nanometer-sized density variations on the subpicosecond and picosecond time scales. The nature of fluctuating density heterogeneities remains a highly debated issue. In the present work, we address the problem of possible association of such density variations with the dynamics of terahertz longitudinal acoustic-like modes in liquid water. Our study is based on the fact that the subpicosecond dynamics of liquid water are essentially governed by the structural relaxation. Using a mode coupling theory approach, we found that for typical values of parameters of liquid water, the dynamic mechanism coming from the combination of the structural relaxation process and the finiteness of the amplitude of terahertz longitudinal acoustic-like mode gives rise to a soliton-like collective mode on a temperature-dependent nanometer length scale. The characteristics of this mode are consistent with the estimates of the amplitudes and temperature-dependent correlation lengths of density fluctuations in liquid water obtained in experiments and simulations. Thus, the fully dynamic mechanism could contribute to the formation and dynamics of fluctuating density heterogeneities. The soliton-like collective excitations suggested by our analysis may be relevant to different phenomena connected with supercooled water and can be expected to be associated with some ultrafast biological processes.

Thermodynamics of supercooled and stretched water: Unifying two-structure description and liquid-vapor spinodal

We have applied a two-structure approach to the description of the thermodynamic properties of supercooled and stretched water, metastable toward vapor, ice, or both, by incorporating the stability limit of liquid with respect to vapor at negative pressures. In addition to the properties of water considered in previous studies, we include new data recently obtained in deeply supercooled and stretched regions. Our model reproduces the experimentally observed anomalies in metastable water up to 400 MPa and down to -140 MPa, and can provide a physically based extrapolation in regions where no measurements are available yet. Moreover, we are able to elucidate the thermodynamic nature of the alternative "states" of liquid water, namely, high-temperature denser water (state A) and "mother-of-ice" lighter water (state B). Based on the internal consistency of the described anomalies and new data on the isothermal compressibility, we exclude the critical-point-free scenario in which the first-order liquid-liquid transition line would continue into the stretched liquid state (doubly metastable) crossing the vapor-liquid spinodal. A "singularity-free" scenario remains an option for explaining supercooled water's anomalies within the framework of two-state thermodynamics; however, the extreme case of the singularity-free scenario, ideal mixing of A and B, seems improbable. We have also clarified the concept of fast interconversion of alternative states in supercooled water as a phenomenological representation of distribution of short-ranged local structures.

Pressure-Induced Densification of Ice Ih under Triaxial Mechanical Compression: Dissociation versus Retention of Crystallinity for Intermediate States in Atomistic and Coarse-Grained Water Models

Molecular-dynamics (MD) simulation of triaxially pressurized ice I_h up to 30 kbar at 240 K (with sudden mechanical pressurization from its ambient-pressure structure) has been carried out with both the single-particle mW and atomistic TIP4P-Ice water potentials on systems of up to ∼1 million molecules, for times of the order of 100 ns. It was found that the TIP4P-Ice systems adopted a high-density liquid state above ∼7 kbar, while densification of the mW systems retained essentially crystalline order, owing to a failure for the tetrahedral network to break down appreciably from its ice I_h lattice structure. Both are intermediate states adopted along the path toward respective thermodynamically stable states (and with pressure removal show reversion to I_h for mW and to supercooled liquid for TIP4P-Ice), similar to recent ice electro-freezing simulations in "No Man's Land". Densification kinetics showed faster mW-system adaptation.

Interstitial Voids and Resultant Density of Liquid Water: A First-Principles Molecular Dynamics Study

Many anomalous properties of water can be explained on the basis of the coexistence of more than one density states: high-density water (HDW) and low-density water (LDW). We investigated these two phases of water molecules through first-principles molecular dynamics simulations using density functional theory (DFT) in conjunction with various van der Waals-corrected exchange and correlation functionals. Different density regions were found to exist due to the difference in short-range and long-range forces present in DFT potentials. These density regions were identified and analyzed on the basis of the distribution of molecules and voids present. We defined a local structure index to distinguish and find the probability of occurrence of these different states. HDW and LDW arise due to the presence of "interstitial water" molecules in between the first and second coordination shells. The population of interstitial water molecules is found to affect the overall dynamics of the system as they change the hydrogen bond pattern.

Is High-Density Amorphous Ice Simply a "Derailed" State along the Ice I to Ice IV Pathway?

The structural nature of high-density amorphous ice (HDA), which forms through low-temperature pressure-induced amorphization of the "ordinary" ice I, is heavily debated. Clarifying this question is important for understanding not only the complex condensed states of H₂O but also in the wider context of pressure-induced amorphization processes, which are encountered across the entire materials spectrum. We first show that ammonium fluoride (NH₄F), which has a similar hydrogen-bonded network to ice I, also undergoes a pressure collapse upon compression at 77 K. However, the product material is not amorphous but NH₄F II, a high-pressure phase isostructural with ice IV. This collapse can be rationalized in terms of a highly effective mechanism. In the case of ice I, the orientational disorder of the water molecules leads to a deviation from this mechanism, and we therefore classify HDA as a "derailed" state along the ice I to ice IV pathway.

On the hydrogen-bond network and the non-Arrhenius transport properties of water

We study the structural and dynamic transformations of SPC/E water with temperature, through molecular dynamics (MD), and discuss the non-Arrhenius behavior of the transport properties and orientational dynamics, and the magnitude of the breakdown of the Stokes-Einstein (SE) and the Stokes-Einstein-Debye (SED) relations, in the light of these transformations. Our results show that deviations from Arrhenius behavior of the self-diffusion at low temperatures cannot be exclusively explained by the reduction of water defects (interstitial waters) and the increase of the local tetrahedrality, thus, suggesting the importance of the slowdown of collective rearrangements. Interestingly we find that at high temperatures (T  ⩾  340 K) water defects lead to a slight increase of the tetrahedrality and a decrease of the self-diffusion, opposite to water at low temperatures. The relative magnitude of the breakdown of the SE and the SED relations is found to be in accord with recent experiments (Dehaoui et al 2015 Proc. Natl Acad. Sci. USA 112 12020) resolving the discrepancy with previous MD results. Further, we show that SPC/E hydrogen-bond (HB) lifetimes deviate from Arrhenious behaviour at low temperatures in contrast with some previous MD studies. This deviation is nevertheless much smaller than that observed for the orientational dynamics and the transport properties of water, consistent with the relaxation times measured by several experimental methods. The HB acceptor exchange dynamics defined here by the acceptor switch and reform (librational dynamics) frequencies exhibit similar Arrhenius deviations, thus explaining to some extent the non-Arrhenius behavior of the transport properties and of the orientational dynamics of water. Our results also show that the fraction of HB switches through a bifurcated pathway follow a power law with the temperature decrease. Thus, at low temperatures HB acceptor switches are less frequent but occur on a faster time scale consistent with the temperature dependence of the ratio of the rotational relaxation times for the different Legendre polynomial ranks.

Growth rate of crystalline ice and the diffusivity of supercooled water from 126 to 262 K

Understanding deeply supercooled water is key to unraveling many of water's anomalous properties. However, developing this understanding has proven difficult due to rapid and uncontrolled crystallization. Using a pulsed-laser-heating technique, we measure the growth rate of crystalline ice, G(T), for 180 K < T < 262 K, that is, deep within water's "no man's land" in ultrahigh-vacuum conditions. Isothermal measurements of G(T) are also made for 126 K ≤ T ≤ 151 K. The self-diffusion of supercooled liquid water, D(T), is obtained from G(T) using the Wilson-Frenkel model of crystal growth. For T > 237 K and P ∼ 10-8 Pa, G(T) and D(T) have super-Arrhenius ("fragile") temperature dependences, but both cross over to Arrhenius ("strong") behavior with a large activation energy in no man's land. The fact that G(T) and D(T) are smoothly varying rules out the hypothesis that liquid water's properties have a singularity at or near 228 K at ambient pressures. However, the results are consistent with a previous prediction for D(T) that assumed no thermodynamic transitions occur in no man's land.

Sub-Tg features of glasses formed by cooling glycerol under pressure - Additional incompatibility of vibrational with configurational states in the depressurized, high density glass

The vibrational state of a glass is naturally incompatible with its configurational state, which makes the glass structurally unstable. When a glass is kept at constant temperature, both the vibrational and configurational states of a glass change with time until it becomes metastable (equilibrium) liquid and the two states become compatible. The process, known as structural relaxation, occurs at a progressively higher rate during heating, and the properties of a glass change accordingly. We add to this incompatibility by depressurizing a glass that had been formed by cooling a liquid under a high pressure, p, and then investigate the effects of the added incompatibility by studying thermal conductivity, κ, and the heat capacity per unit volume ρC_p of the depressurized glass. We use glycerol for the purpose and study first the changes in the features of κ and of ρC_p during glass formation on cooling under a set of different p. We then partially depressurize the glass and study the effect of the p-induced instability on the features of κ and ρC_p as the glass is isobarically heated to the liquid state. At a given low p, the glass configuration that was formed by cooling at high-p had a higher κ than the glass configuration that was formed by cooling at a low p. The difference is more when the glass is formed at a higher p and/or is depressurized to a lower p. On heating at a low p, its κ decreases before its glass-liquid transition range is reached. The effect is the opposite of the increase in κ observed on heating a glass at the same p under which it was formed. It is caused by thermally assisted loss of the added incompatibility of configurational and vibrational states of a high-p formed glass kept at low p. If a glass formed under a low-p is pressurized and then heated under high p, it would show the opposite effect, i.e., its κ would first increase to its high p value before its glass-to-liquid transition range.

Micro- and Nano-Liquid Phases Coexistent with Ice as Separation and Reaction Media

Ice has a variety of scientifically interesting features, some of which have not been reasonably interpreted despite substantial efforts by researchers. Most chemical studies of ice have focused on the elucidation of its physicochemical nature and its roles in the natural environment. Ice often contains impurities, such as salts, and in such cases, a liquid phase coexists with solid ice over a wide temperature range. This impure ice also acts as a cryoreactor, governing the circulation of chemical species of environmental importance. Reactions and phenomena occurring in this liquid phase show features different from those seen in normal bulk aqueous solutions. In the present account, we discuss the chemical characteristics of the liquid phase that develops in a frozen aqueous phase and show how novel analytical systems can be designed based on he features of the liquid phase which are predictable in some cases but unpredictable in others.

Diversity at the Water-Metal Interface: Metal, Water Thickness, and Confinement Effects

The structure and properties of water films in contact with metal surfaces are crucial to understand the chemical and electrochemical processes involved in energy-related technologies. The nature of thin water films on Pd, Pt, and Ru has been investigated by first-principles molecular dynamics to assess how the chemistry at the water-metal surface is responsible for the diversity in the behavior of the water layers closer to the metal. The characteristics of liquid water: the radial distribution functions, coordination, and fragment speciation appear only for unconfined water layers of a minimum of 1.4 nm thick. In addition, the water layer is denser in the region closest to the metal for Pd and Pt, where seven- and five-membered ring motifs appear. These patterns are identical to those identified by scanning tunneling microscopy for isolated water bilayers. On Ru densification at the interface is not observed, water dissociates, and protons and hydroxyl groups are locked at the surface. Therefore, the acid-base properties in the area close to the metal are not perturbed, in agreement with experiments, and the bulk water resembles an electric double layer. Confinement affects water making it closer to ice for both structural and dynamic properties, thus being responsible for the higher viscosity experimentally found at the nanoscale. All these contributions modify the solvation of reactants and products at the water-metal interface and will affect the catalytic and electrocatalytic properties of the surface.