Water polyamorphism: Reversibility and (dis)continuity

Unveiling a common phase transition pathway of high-density amorphous ices through time-resolved x-ray scattering

Here, we investigate the hypothesis that despite the existence of at least two high-density amorphous ices, only one high-density liquid state exists in water. We prepared a very-high-density amorphous ice (VHDA) sample and rapidly increased its temperature to around 205 ± 10 K using laser-induced isochoric heating. This temperature falls within the so-called "no-man's land" well above the glass-liquid transition, wherein the IR laser pulse creates a metastable liquid state. Subsequently, this high-density liquid (HDL) state of water decompresses over time, and we examined the time-dependent structural changes using short x-ray pulses from a free electron laser. We observed a liquid-liquid transition to low-density liquid water (LDL) over time scales ranging from 20 ns to 3 μs, consistent with previous experimental results using expanded high-density amorphous ice (eHDA) as the initial state. In addition, the resulting LDL derived both from VHDA and eHDA displays similar density and degree of inhomogeneity. Our observation supports the idea that regardless of the initial annealing states of the high-density amorphous ices, the same HDL and final LDL states are reached at temperatures around 205 K.

Phase Diagram of Aqueous Solutions of LiCl: a Study of Concentration Effects on the Anomalies of Water

We perform molecular dynamics simulations in order to study thermodynamics and the structure of supercooled aqueous solutions of lithium chloride (LiCl) at concentrations c = 0.678 and 2.034 mol/kg. We model the solvent using the TIP4P/2005 potential and the ions using the Madrid-2019 force field, a force field particularly suited for studying this solution. We find that, for c = 0.678 mol/kg, the behavior of the equation of state, studied in the P-T plane, indicates the presence of a liquid-liquid phase transition, similar to what was previously found for bulk water. We estimate the position of the liquid-liquid critical point to be at T_c ≈ 174 K, P_c ≈ 1775 bar, and ρ_c ≈ 1.065 g/cm³. When the concentration is tripled to c = 2.034 mol/kg, no critical point is observed, indicating its possible disappearance at this concentration. We also study the water-water and water-ions structure in the two solutions, and we find that at the concentrations examined the effect of ions on the water-water structure is not strong, and all the features found in bulk water are preserved. We also calculate the hydration number of the Li and Cl ions, and in line with experiments, we find the value of 4 for Li⁺ and between 5.5 and 6 for Cl^-, confirming the good performances of the Madrid-2019 force field.

Liquid-liquid phase separation in supercooled water from ultrafast heating of low-density amorphous ice

Recent experiments continue to find evidence for a liquid-liquid phase transition (LLPT) in supercooled water, which would unify our understanding of the anomalous properties of liquid water and amorphous ice. These experiments are challenging because the proposed LLPT occurs under extreme metastable conditions where the liquid freezes to a crystal on a very short time scale. Here, we analyze models for the LLPT to show that coexistence of distinct high-density and low-density liquid phases may be observed by subjecting low-density amorphous (LDA) ice to ultrafast heating. We then describe experiments in which we heat LDA ice to near the predicted critical point of the LLPT by an ultrafast infrared laser pulse, following which we measure the structure factor using femtosecond x-ray laser pulses. Consistent with our predictions, we observe a LLPT occurring on a time scale < 100 ns and widely separated from ice formation, which begins at times >1 μs.

Molecular dynamics analysis of elastic properties and new phase formation during amorphous ices transformations

Unlike conventional first-order phase transitions, the kinetics of amorphous-amorphous transitions has been much less studied. The ultrasonic experiments on the transformations between low-density and high-density amorphous ice induced by pressure or heating provided the pressure and temperature dependencies of elastic moduli. In this article, we make an attempt to build a microscopic picture of these experimentally studied transformations using the molecular dynamics method with the TIP4P/Ice water model. We study carefully the dependence of the results of elastic constants calculations on the deformation rates. The system size effects are considered as well. The comparison with the experimental data enriches our understanding of the transitions observed. Our modeling gives new information about the formation mechanisms of new phase clusters during the transition between low-density and high-density amorphous ices. We analyse the applicability of the term “nucleation” for these processes.

Oxygen NMR of high-density and low-density amorphous ice

Using oxygen-17 as a nuclear probe, spin relaxometry was applied to study the high-density and low-density states of amorphous ice, covering temperatures below and somewhat above their glass transitions. These findings are put in perspective with results from deuteron nuclear magnetic resonance and with calculations based on dielectrically detected correlation times. This comparison reveals the presence of a wide distribution of correlation times. Furthermore, oxygen-17 central-transition echo spectra were recorded for wide ranges of temperature and pulse spacing. The spectra cannot be described by a single set of quadrupolar parameters, suggesting a distribution of H–O–H opening angles that is broader for high-density than for low-density amorphous ice. Simulations of the pulse separation dependent spin-echo spectra for various scenarios demonstrate that a small-step frequency diffusion process, assigned to the presence of homonuclear oxygen–oxygen interactions, determines the shape evolution of the pulse-separation-dependent spectra.

Direct observation of reversible liquid-liquid transition in a trehalose aqueous solution

Recent studies on liquid water suggest that the two liquid waters exist in the supercooled temperature region and that their existence relates to the anomalous behavior of low-temperature liquid water such as the maximum density at 4 °C. However, the experimental investigation of two liquid waters is difficult because of the rapid crystallization. In this study, a reversible liquid–liquid transition in a trehalose aqueous solution by the change in pressure was observed directly. This result suggests strongly that two liquid waters exist in the aqueous solution. This study has implications for wide fields related to liquid water, such as solution chemistry, cryobiology, meteorology, and food engineering.

Water forms two glassy waters, low-density and high-density amorphs, which undergo a reversible polyamorphic transition with the change in pressure. The two glassy waters transform into the different liquids, low-density liquid (LDL) and high-density liquid (HDL), at high temperatures. It is predicted that the two liquid waters also undergo a liquid–liquid transition (LLT). However, the reversible LLT, particularly the LDL-to-HDL transition, has not been observed directly due to rapid crystallization. Here, I prepared a glassy dilute trehalose aqueous solution (0.020 molar fraction) without segregation and measured the isothermal volume change at 0.01 to 1.00 GPa below 160 K. The polyamorphic transition and the glass-to-liquid transition for the high-density and low-density solutions were examined, and the liquid region where both LDL and HDL existed was determined. The results show that the reversible polyamorphic transition induced by the pressure change above 140 K is the LLT. That is, the transition from LDL to HDL is observed. Moreover, the pressure hysteresis of LLT suggests strongly that the LLT has a first-order nature. The direct observation of the reversible LLT in the trehalose aqueous solution has implications for understanding not only the liquid–liquid critical point hypothesis of pure water but also the relation between aqueous solution and water polyamorphism.

The impact of temperature and unwanted impurities on slow compression of ice

For slowly compressed hexagonal ice pressure-induced amorphisation to high-density amorphous ice (HDA) takes place below and at 130 K, but polymorphic transformation to ice IX takes place at 140–170 K. Stable ice II only forms above 170 K. Ice IX impurities trigger ice IX growth even at 120 K. HDA and ice IX are equally long-lived, where both can be regarded as metastable phases.

Slow compression of ice at 100 K usually results in pressure-amorphisation, unless there are ice IX seeds and temperature gradients.

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.

Quantum polyamorphism in compressed distinguishable helium-4

We demonstrate that two amorphous solid states can exist in ⁴He consisting of distinguishable Boltzmann atoms under compressed conditions. The isothermal compression of normal or supercritical fluid ⁴He was conducted at 3-25 K using the isobaric-isothermal path integral centroid molecular dynamics simulation. The compression of fluid first produced the low-dispersion amorphous (LDA) state possessing modest extension of atomic necklaces. Further isothermal compression up to the order of 10 kbar to 1 Mbar or an isobaric cooling of LDA induced the transition to the high-dispersion amorphous (HDA) state. The HDA was characterized by long quantum wavelengths of atoms extended over several Angstroms and the promotion of atomic residual diffusion. They were related to the quantum tunneling of atoms bestriding the potential saddle points in this glass. The change in pressure or temperature induced the LDA-HDA transition reversibly with hysteresis, while it resembled the coil-globule transition of classical polymers. The HDA had lower kinetic and higher Gibbs free energies than the LDA at close temperature. The HDA was absent at T ≥ 13 K, while the LDA-HDA transition pressure significantly decreased with lowering temperature. The LDA and HDA correspond to the trapped and tunneling regimes proposed by Markland et al. [J. Chem. Phys. 136, 074511 (2012)], respectively. The same reentrant behavior as they found was observed for the expansion factor of the quantum wavelength as well as for atomic diffusivity.

Glass polyamorphism in gallium: Two amorphous solid states and their transformation on the potential energy landscape

Using the potential energy landscape (PEL) formalism and molecular dynamics simulations, we investigate a phase transformation between two amorphous solid states of gallium, namely, a low-density amorphous solid (LDA) and a high-density amorphous solid (HDA), and compare with its equilibrium counterpart, the liquid-liquid phase transition (LLPT). It is found that on the PEL, the signatures of the out-of-equilibrium LDA-HDA transition are reminiscent of those of the equilibrium LLPT in terms of pressure, inherent structure pressure, inherent structure energy, and shape function, indicating that the LDA-HDA transformation is a first-order-like transition. However, differences are also found between the out-of-equilibrium phase transition and the equilibrium one, for example, the path from LDA to HDA on the PEL cannot be accessed by the path from LDL to HDL. Our results also suggest that the signatures of the out-of-equilibrium transition in gallium are rather general features of systems with an accessible LLPT-not only systems with pairwise interactions but also those with many-body interactions. This finding is of crucial importance for obtaining a deeper understanding of the nature of transitions in the polyamorphic family.

Absence of the liquid-liquid phase transition in aqueous ionic liquids

The anomalies of supercooled water may be explained by an underlying liquid-liquid phase transition (LLPT) between high- and low-density states. Recently, its observation at 185 K was inferred using solutions containing aqueous ionic liquids at a solute mole fraction of x=0.156 [Woutersen et al., Science 359, 1127 (2018)10.1126/science.aao7049]. We employ x-ray diffraction, calorimetry, and dilatometry on these hydrazinium trifluoroacetate solutions at x=0.00-0.40 to show that the transition at 185 K is not related to a genuine LLPT of water. Continuous densification upon compression, continuous changes of halo position, and absence of thermal signatures for a high- to low-density transition rule out the possibility of an LLPT for x≥0.13. The data show that employing sophisticated solutions adds a layer of complexity that hampers extrapolation of the LLPT concept from one- to two-component systems. The possibility of an LLPT can only be probed for pure water or sufficiently dilute aqueous solutions.

Glass polymorphism and liquid-liquid phase transition in aqueous solutions: experiments and computer simulations

One of the most intriguing anomalies of water is its ability to exist as distinct amorphous ice forms (glass polymorphism or polyamorphism). This resonates well with the possible first-order liquid-liquid phase transition (LLPT) in the supercooled state, where ice is the stable phase. In this Perspective, we review experiments and computer simulations that search for LLPT and polyamorphism in aqueous solutions containing salts and alcohols. Most studies on ionic solutes are devoted to NaCl and LiCl; studies on alcohols have mainly focused on glycerol. Less attention has been paid to protein solutions and hydrophobic solutes, even though they reveal promising avenues. While all solutions show polyamorphism and an LLPT only in dilute, sub-eutectic mixtures, there are differences regarding the nature of the transition. Isocompositional transitions for varying mole fractions are observed in alcohol but not in ionic solutions. This is because water can surround alcohol molecules either in a low- or high-density configuration whereas for ionic solutes, the water ion hydration shell is forced into high-density structures. Consequently, the polyamorphic transition and the LLPT are prevented near the ions, but take place in patches of water within the solutions. We highlight discrepancies and different interpretations within the experimental community as well as the key challenges that need consideration when comparing experiments and simulations. We point out where reinterpretation of past studies helps to draw a unified, consistent picture. In addition to the literature review, we provide original experimental results. A list of eleven open questions that need further consideration is identified.

Structural differences between unannealed and expanded high-density amorphous ice based on isotope substitution neutron diffraction

We here report isotope substitution neutron diffraction experiments on two variants of high-density amorphous ice (HDA): its unannealed form prepared via pressure-induced amorphization of hexagonal ice at 77 K, and its expanded form prepared via decompression of very-high density amorphous ice at 140 K. The latter is about 17 K more stable thermally, so that it can be heated beyond its glass-to-liquid transition to the ultraviscous liquid form at ambient pressure. The structural origin for this large thermal difference and the possibility to reach the deeply supercooled liquid state has not yet been understood. Here we reveal that the origin for this difference is found in the intermediate range structure, beyond about 3.6 Å. The hydration shell markedly differs at about 6 Å. The local order, by contrast, including the first as well as the interstitial space between first and second shell is very similar for both. 'eHDA' that is decompressed to 0.20 GPa instead of 0.07 GPa is here revealed to be rather far away from well-relaxed eHDA. Instead it turns out to be roughly halfway between VHDA and eHDA - stressing the importance for decompressing VHDA to at least 0.10 GPa to make an eHDA sample of good quality.

Glass polymorphism in TIP4P/2005 water: A description based on the potential energy landscape formalism

The potential energy landscape (PEL) formalism is a statistical mechanical approach to describe supercooled liquids and glasses. Here, we use the PEL formalism to study the pressure-induced transformations between low-density amorphous ice (LDA) and high-density amorphous ice (HDA) using computer simulations of the TIP4P/2005 molecular model of water. We find that the properties of the PEL sampled by the system during the LDA-HDA transformation exhibit anomalous behavior. In particular, at conditions where the change in density during the LDA-HDA transformation is approximately discontinuous, reminiscent of a first-order phase transition, we find that (i) the inherent structure (IS) energy, e_IS(V), is a concave function of the volume and (ii) the IS pressure, P_IS(V), exhibits a van der Waals-like loop. In addition, the curvature of the PEL at the IS is anomalous, a nonmonotonic function of V. In agreement with previous studies, our work suggests that conditions (i) and (ii) are necessary (but not sufficient) signatures of the PEL for the LDA-HDA transformation to be reminiscent of a first-order phase transition. We also find that one can identify two different regions of the PEL, one associated with LDA and another with HDA. Our computer simulations are performed using a wide range of compression/decompression and cooling rates. In particular, our slowest cooling rate (0.01 K/ns) is within the experimental rates employed in hyperquenching experiments to produce LDA. Interestingly, the LDA-HDA transformation pressure that we obtain at T = 80 K and at different rates extrapolates remarkably well to the corresponding experimental pressure.

X-ray studies of the transformation from high- to low-density amorphous water

Here we report about the structural evolution during the conversion from high-density amorphous ices at ambient pressure to the low-density state. Using high-energy X-ray diffraction, we have monitored the transformation by following in reciprocal space the structure factor S_OO( Q) and derived in real space the pair distribution function g_OO( r). Heating equilibrated high-density amorphous ice (eHDA) at a fast rate (4 K min^-1), the transition to the low-density form occurs very rapidly, while domains of both high- and low-density coexist. On the other hand, the transition in the case of unannealed HDA (uHDA) and very-high-density amorphous ice is more complex and of continuous nature. The direct comparison of eHDA and uHDA indicates that the molecular structure of uHDA contains a larger amount of tetrahedral motives. The different crystallization behaviour of the derived low-density amorphous states is interpreted as emanating from increased tetrahedral coordination present in uHDA. This article is part of the theme issue 'The physics and chemistry of ice: scaffolding across scales, from the viability of life to the formation of planets'.

Link between molecular mobility and order parameter during liquid-liquid transition of a molecular liquid

Liquid-liquid transition (LLT) is the transformation of one liquid to another via first-order phase transition. For example, LLT in a molecular liquid, triphenyl phosphite, is macroscopically the transformation from liquid I in a supercooled state to liquid II in a glassy state. Reflecting the transformation from the liquid to glassy state, the LLT is accompanied by considerable slowing down of overall molecular dynamics, but little is known about how this proceeds at a molecular level coupled with the evolution of the order parameter. We report such information by performing time-resolved simultaneous measurements of dielectric spectroscopy and phase contrast microscopy/Raman spectroscopy by using a dielectric cell with transparent electrodes. We find that the temporal change in molecular mobility crucially depends on whether LLT is nucleation growth type occurring in the metastable state or SD type occurring in the unstable state. Furthermore, our results suggest that the molecular mobility is controlled by the local order parameter: more specifically, the local activation energy of molecular rotation is controlled by the local fraction of locally favored structures formed in the liquid. Our study sheds light on the temporal change in the molecular dynamics during LLT and its link to the order parameter evolution.

Thermodynamically driven assemblies and liquid-liquid phase separations in biology

The sustenance of life depends on the high degree of organization that prevails through different levels of living organisms, from subcellular structures such as biomolecular complexes and organelles to tissues and organs. The physical origin of such organization is not fully understood, and even though it is clear that cells and organisms cannot maintain their integrity without consuming energy, there is growing evidence that individual assembly processes can be thermodynamically driven and occur spontaneously due to changes in thermodynamic variables such as intermolecular interactions and concentration. Understanding the phase separation in vivo requires a multidisciplinary approach, integrating the theory and physics of phase separation with experimental and computational techniques. This paper aims at providing a brief overview of the physics of phase separation and its biological implications, with a particular focus on the assembly of membraneless organelles. We discuss the underlying physical principles of phase separation from its thermodynamics to its kinetics. We also overview the wide range of methods utilized for experimental verification and characterization of phase separation of membraneless organelles, as well as the utility of molecular simulations rooted in thermodynamics and statistical physics in understanding the governing principles of thermodynamically driven biological self-assembly processes.

X-ray Scattering and O-O Pair-Distribution Functions of Amorphous Ices

The structure factor and oxygen–oxygen pair-distribution functions of amorphous ices at liquid nitrogen temperature (T = 77 K) have been derived from wide-angle X-ray scattering (WAXS) up to interatomic distances of r = 23 Å, where local structure differences between the amorphous ices can be seen for the entire range. The distances to the first coordination shell for low-, high-, and very-high-density amorphous ice (LDA, HDA, VHDA) were determined to be 2.75, 2.78, and 2.80 Å, respectively, with high accuracy due to measurements up to a large momentum transfer of 23 Å^–1. Similarities in pair-distribution functions between LDA and supercooled water at 254.1 K, HDA and liquid water at 365.9 K, and VHDA and high-pressure liquid water were found up to around 8 Å, but beyond that at longer distances, the similarities were lost. In addition, the structure of the high-density amorphous ices was compared to high-pressure crystalline ices IV, IX , and XII, and conclusions were drawn about the local ordering.

Heating- and pressure-induced transformations in amorphous and hexagonal ice: A computer simulation study using the TIP4P/2005 model

We characterize the phase behavior of glassy water by performing extensive out-of-equilibrium molecular dynamics simulations using the TIP4P/2005 water model. Specifically, we study (i) the pressure-induced transformations between low-density (LDA) and high-density amorphous ice (HDA), (ii) the pressure-induced amorphization (PIA) of hexagonal ice (I_h), (iii) the heating-induced LDA-to-HDA transformation at high pressures, (iv) the heating-induced HDA-to-LDA transformation at low and negative pressures, (v) the glass transition temperatures of LDA and HDA as a function of pressure, and (vi) the limit of stability of LDA upon isobaric heating and isothermal decompression (at negative pressures). These transformations are studied systematically, over a wide range of temperatures and pressures, allowing us to construct a P-T phase diagram for glassy TIP4P/2005 water. Our results are in qualitative agreement with experimental observations and with the P-T phase diagram obtained for glassy ST2 water that exhibits a liquid-liquid phase transition and critical point. We also discuss the mechanism for PIA of ice I_h and show that this is a two-step process where first, the hydrogen-bond network (HBN) is distorted and then the HBN abruptly collapses. Remarkably, the collapse of the HB in ice I_h occurs when the average molecular orientations order, a measure of the tetrahedrality of the HBN, is of the same order as in LDA, suggesting a common mechanism for the LDA-to-HDA and I_h-to-HDA transformations.

Relaxation dynamics and transformation kinetics of deeply supercooled water: Temperature, pressure, doping, and proton/deuteron isotope effects

Above its glass transition, the equilibrated high-density amorphous ice (HDA) transforms to the low-density pendant (LDA). The temperature dependence of the transformation is monitored at ambient pressure using dielectric spectroscopy and at elevated pressures using dilatometry. It is found that near the glass transition temperature of deuterated samples, the transformation kinetics is 300 times slower than the structural relaxation, while for protonated samples, the time scale separation is at least 30 000 and insensitive to doping. The kinetics of the HDA to LDA transformation lacks a proton/deuteron isotope effect, revealing that this process is dominated by the restructuring of the oxygen network. The x-ray diffraction experiments performed on samples at intermediate transition stages reflect a linear combination of the LDA and HDA patterns implying a macroscopic phase separation, instead of a local intermixing of the two amorphous states.

Thermally induced amorphous to amorphous transition in hot-compressed silica glass

In situ Raman and Brillouin light scattering techniques were used to study thermally induced high-density amorphous (HDA) to low-density amorphous (LDA) transition in silica glass densified in hot compression (up to 8 GPa at 1100 °C). Hot-compressed silica samples are shown to retain structural and mechanical stability through 600 °C or greater, with reduced sensitivity in elastic response to temperature as compared with pristine silica glass. Given sufficient thermal energy to overcome the energy barrier, the compacted structure of the HDA silica reverts back to the LDA state. The onset temperature for the HDA to LDA transition depends on the degree of densification during hot compression, commencing at lower temperatures for samples with higher density, but all finishing within a temperature range of 250-300 °C. Our studies show that the HDA to LDA transition at high temperatures in hot-compressed samples is different from the gradual changes starting from room temperature in cold-compressed silica glass, indicating greater structural homogeneity achieved by hot compression. Furthermore, the structure and properties of hot-compressed silica glass change continuously during the thermally induced HDA to LDA transition, in contrast to the abrupt and first-order-like polyamorphic transitions in amorphous ice. Different HDA to LDA transition mechanisms in amorphous silica and amorphous ice are explained by their different energy landscapes.

On the crystallisation temperature of very high-density amorphous ice

The influence of the protocol of preparation on the crystallisation temperature TX of very high-density amorphous ice (VHDA) was studied by varying the annealing pressure (1.1, 1.6 and 1.9 GPa) and temperature (160, 167 and 175 K, respectively). TX increases by up to 4 K in the pressure range of 0.7 to 1.8 GPa for samples annealed at 1.9 GPa compared to samples annealed at 1.1 GPa. Concomitantly, secondary crystallisation channels are suppressed, indicating the absence of structural inhomogeneities. For VHDA prepared at 1.1 GPa and 1.6 GPa our results indicate such inhomogeneities, which we regard to be incompletely amorphized, distorted nanodomains of hexagonal ice that cannot be detected through X-ray diffraction experiments. VHDA prepared at high pressures and temperatures thus represents the amorphous state of water at >0.7 GPa least affected by nanocrystals that has been described so far. We expect the TX obtained for the samples prepared in this manner to be close to the ultimate limit, i.e., we do not consider it possible to raise the low-temperature border to the no-man's land notably further by changing the preparation protocol. An additional, considerable increase in this border will only be possible by working at much shorter time-scales, e.g., by employing fast heating experiments.

Calorimetric study of water's two glass transitions in the presence of LiCl

A DSC study of dilute glassy LiCl aqueous solutions in the water-dominated regime provides direct evidence of a glass-to-liquid transition in expanded high density amorphous (eHDA)-type solutions. Similarly, low density amorphous ice (LDA) exhibits a glass transition prior to crystallization to ice Ic. Both glass transition temperatures are independent of the salt concentration, whereas the magnitude of the heat capacity increase differs. By contrast to pure water, the glass transition endpoint for LDA can be accessed in LiCl aqueous solutions above 0.01 mole fraction. Furthermore, we also reveal the endpoint for HDA's glass transition, solving the question on the width of both glass transitions. This suggests that both equilibrated HDL and LDL can be accessed in dilute LiCl solutions, supporting the liquid-liquid transition scenario to understand water's anomalies.

High-density amorphous ice: nucleation of nanosized low-density amorphous ice

The pressure dependence of the crystallization temperature of different forms of expanded high-density amorphous ice (eHDA) was scrutinized. Crystallization at pressures 0.05-0.30 GPa was followed using volumetry and powder x-ray diffraction. eHDA samples were prepared via isothermal decompression of very high-density amorphous ice at 140 K to different end pressures between 0.07-0.30 GPa (eHDA^0.07-0.3). At 0.05-0.17 GPa the crystallization line T _x (p) of all eHDA variants is the same. At pressures  >0.17 GPa, all eHDA samples decompressed to pressures  <0.20 GPa exhibit significantly lower T _x values than eHDA^0.2 and eHDA^0.3. We rationalize our findings with the presence of nanoscaled low-density amorphous ice (LDA) seeds that nucleate in eHDA when it is decompressed to pressures  <0.20 GPa at 140 K. Below ~0.17 GPa, these nanosized LDA domains are latent within the HDA matrix, exhibiting no effect on T _x of eHDA^<0.2. Upon heating at pressures  ⩾0.17 GPa, these nanosized LDA nuclei transform to ice IX nuclei. They are favored sites for crystallization and, hence, lower T _x . By comparing crystallization experiments of bulk LDA with the ones involving nanosized LDA we are able to estimate the Laplace pressure and radius of ~0.3-0.8 nm for the nanodomains of LDA. The nucleation of LDA in eHDA revealed here is evidence for the first-order-like nature of the HDA  →  LDA transition, supporting water's liquid-liquid transition scenarios.

Supercooled and glassy water: Metastable liquid(s), amorphous solid(s), and a no-man's land

We review the recent research on supercooled and glassy water, focusing on the possible origins of its complex behavior. We stress the central role played by the strong directionality of the water-water interaction and by the competition between local energy, local entropy, and local density. In this context we discuss the phenomenon of polyamorphism (i.e., the existence of more than one disordered solid state), emphasizing both the role of the preparation protocols and the transformation between the different disordered ices. Finally, we present the ongoing debate on the possibility of linking polyamorphism with a liquid-liquid transition that could take place in the no-man's land, the temperature-pressure window in which homogeneous nucleation prevents the investigation of water in its metastable liquid form.

Slow Dynamics and Structure of Supercooled Water in Confinement

We review our simulation results on properties of supercooled confined water. We consider two situations: water confined in a hydrophilic pore that mimics an MCM-41 environment and water at interface with a protein. The behavior upon cooling of the α relaxation of water in both environments is well interpreted in terms of the Mode Coupling Theory of glassy dynamics. Moreover, we find a crossover from a fragile to a strong regime. We relate this crossover to the crossing of the Widom line emanating from the liquid-liquid critical point, and in confinement we connect this crossover also to a crossover of the two body excess entropy of water upon cooling. Hydration water exhibits a second, distinctly slower relaxation caused by its dynamical coupling with the protein. The crossover upon cooling of this long relaxation is related to the protein dynamics.

Effect of water polyamorphism on the molecular vibrations of glycerol in its glassy aqueous solutions

A glassy dilute glycerol-water solution undergoes a mutual polyamorphic transition relating to the transition between high- and low-density amorphous ices of solvent water. The polyamorphic transition behavior depends on the glycerol concentration, indicating that the glycerol affects the water polyamorphism. Here, we used the glassy dilute glycerol-water solution of the solute molar fraction of 0.07 and examined the effect of the polyamorphic change in solvent water on the molecular vibrations of glycerol via Raman spectroscopy. It is found that the molecular vibration of glycerol in high-density liquid like solvent water is different from that in the low-density liquid like solvent water and that the change in the molecular vibration of glycerol is synchronized with the polyamorphic transition of solvent water. The dynamical change of the solute molecule relates to the polyamorphic state of solvent water. This result suggests that the polyamorphic fluctuation of water structure emanated from the presumed liquid-liquid critical point plays an important role for the function of aqueous solution under an ambient condition such as the conformational stability of solute, the functional expression of solute, and so on.

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.

Glass polymorphism in glycerol-water mixtures: I. A computer simulation study

We perform out-of-equilibrium molecular dynamics (MD) simulations of water-glycerol mixtures in the glass state. Specifically, we study the transformations between low-density (LDA) and high-density amorphous (HDA) forms of these mixtures induced by compression/decompression at constant temperature. Our MD simulations reproduce qualitatively the density changes observed in experiments. Specifically, the LDA-HDA transformation becomes (i) smoother and (ii) the hysteresis in a compression/decompression cycle decreases as T and/or glycerol content increase. This is surprising given the fast compression/decompression rates (relative to experiments) accessible in MD simulations. We study mixtures with glycerol molar concentration χ(g) = 0-13% and find that, for the present mixture models and rates, the LDA-HDA transformation is detectable up to χ(g) ≈ 5%. As the concentration increases, the density of the starting glass (i.e., LDA at approximately χ(g) ≤ 5%) rapidly increases while, instead, the density of HDA remains practically constant. Accordingly, the LDA state and hence glass polymorphism become inaccessible for glassy mixtures with approximately χ(g) > 5%. We present an analysis of the molecular-level changes underlying the LDA-HDA transformation. As observed in pure glassy water, during the LDA-to-HDA transformation, water molecules within the mixture approach each other, moving from the second to the first hydration shell and filling the first interstitial shell of water molecules. Interestingly, similar changes also occur around glycerol OH groups. It follows that glycerol OH groups contribute to the density increase during the LDA-HDA transformation. An analysis of the hydrogen bond (HB)-network of the mixtures shows that the LDA-HDA transformation is accompanied by minor changes in the number of HBs of water and glycerol. Instead, large changes in glycerol and water coordination numbers occur. We also perform a detailed analysis of the effects that the glycerol force field (FF) has on our results. By comparing MD simulations using two different glycerol models, we find that glycerol conformations indeed depend on the FF employed. Yet, the thermodynamic and microscopic mechanisms accompanying the LDA-HDA transformation and hence, our main results, do not. This work is accompanied by an experimental report where we study the glass polymorphism in glycerol-water mixtures prepared by isobaric cooling at 1 bar.

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.

The glass transition in high-density amorphous ice

There has been a long controversy regarding the glass transition in low-density amorphous ice (LDA). The central question is whether or not it transforms to an ultraviscous liquid state above 136 K at ambient pressure prior to crystallization. Currently, the most widespread interpretation of the experimental findings is in terms of a transformation to a superstrong liquid above 136 K. In the last decade some work has also been devoted to the study of the glass transition in high-density amorphous ice (HDA) which is in the focus of the present review. At ambient pressure HDA is metastable against both ice I and LDA, whereas at > 0.2 GPa HDA is no longer metastable against LDA, but merely against high-pressure forms of crystalline ice. The first experimental observation interpreted as the glass transition of HDA was made using in situ methods by Mishima, who reported a glass transition temperature Tg of 160 K at 0.40 GPa. Soon thereafter Andersson and Inaba reported a much lower glass transition temperature of 122 K at 1.0 GPa. Based on the pressure dependence of HDA's Tg measured in Innsbruck, we suggest that they were in fact probing the distinct glass transition of very high-density amorphous ice (VHDA). Very recently the glass transition in HDA was also observed at ambient pressure at 116 K. That is, LDA and HDA show two distinct glass transitions, clearly separated by about 20 K at ambient pressure. In summary, this suggests that three glass transition lines can be defined in the p-T plane for LDA, HDA, and VHDA.