Structure of a new dense amorphous ice

Generating interstitial water within the persisting tetrahedral H-bond network explains density increase upon compressing liquid water

Despite its ubiquitous nature, the atomic structure of water in its liquid state is still controversially debated. We use a combination of X-ray Raman scattering spectroscopy in conjunction with ab initio and path integral molecular dynamics simulations to study the local atomic and electronic structure of water under high pressure conditions. Systematically increasing fingerprints of non-hydrogen-bonded H[Formula: see text]O molecules in the first hydration shell are identified in the experimental and computational oxygen K-edge excitation spectra. This provides evidence for a compaction mechanism in terms of a continuous collapse of the second hydration shell with increasing pressure via generation of interstitial water within locally tetrahedral hydrogen-bonding environments.

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.

A continuum of amorphous ices between low-density and high-density amorphous ice

Amorphous ices are usually classified as belonging to low-density or high-density amorphous ice (LDA and HDA) with densities ρLDA ≈ 0.94 g/cm3 and ρHDA ≈ 1.15-1.17 g/cm3. However, a recent experiment crushing hexagonal ice (ball-milling) produced a medium-density amorphous ice (MDA, ρMDA ≈ 1.06 g/cm3) adding complexity to our understanding of amorphous ice and the phase diagram of supercooled water. Motivated by the discovery of MDA, we perform computer simulations where amorphous ices are produced by isobaric cooling and isothermal compression/decompression. Our results show that, depending on the pressure employed, isobaric cooling can generate a continuum of amorphous ices with densities that expand in between those of LDA and HDA (briefly, intermediate amorphous ices, IA). In particular, the IA generated at P ≈ 125 MPa has a remarkably similar density and average structure as MDA, implying that MDA is not unique. Using the potential energy landscape formalism, we provide an intuitive qualitative understanding of the nature of LDA, HDA, and the IA generated at different pressures. In this view, LDA and HDA occupy specific and well-separated regions of the PEL; the IA prepared at P = 125 MPa is located in the intermediate region of the PEL that separates LDA and HDA.

Kinetics and Mechanisms of Pressure-Induced Ice Amorphization and Polyamorphic Transitions in a Machine-Learned Coarse-Grained Water Model

Water glasses have attracted considerable attention due to their potential connection to a liquid-liquid transition in supercooled water. Here we use molecular simulations to investigate the formation and phase behavior of water glasses using the machine-learned bond-order parameter (ML-BOP) water model. We produce glasses through hyperquenching of water, pressure-induced amorphization (PIA) of ice, and pressure-induced polyamorphic transformations. We find that PIA of polycrystalline ice occurs at a lower pressure than that of monocrystalline ice and through a different mechanism. The temperature dependence of the amorphization pressure of polycrystalline ice for ML-BOP agrees with that in experiments. We also find that ML-BOP accurately reproduces the density, coordination number, and structural features of low-density (LDA), high-density (HDA), and very high-density (VHDA) amorphous water glasses. ML-BOP accurately reproduces the experimental radial distribution function of LDA but overpredicts the minimum between the first two shells in high-density glasses. We examine the kinetics and mechanism of the transformation between low-density and high-density glasses and find that the sharp nature of these transitions in ML-BOP is similar to that in experiments and all-atom water models with a liquid-liquid transition. Transitions between ML-BOP glasses occur through a spinodal-like mechanism, similar to ice crystallization from LDA. Both glass-to-glass and glass-to-ice transformations have Avrami-Kolmogorov kinetics with exponent n = 1.5 ± 0.2 in experiments and simulations. Importantly, ML-BOP reproduces the competition between crystallization and HDA→LDA transition above the glass transition temperature T_g, and separation of their time scales below T_g, observed also in experiments. These findings demonstrate the ability of ML-BOP to accurately reproduce water properties across various regimes, making it a promising model for addressing the competition between polyamorphic transitions and crystallization in water and solutions.

Free energy calculations and unbiased molecular dynamics targeting the liquid-liquid transition in water no man's land

The existence of a first-order phase transition between a low-density liquid (LDL) and a high-density liquid (HDL) form of supercooled water has been a central and highly debated issue of physics and chemistry for the last three decades. We present a computational study that allows us to determine the free-energy landscapes of supercooled water over a wide range of pressure and temperature conditions using the TIP4P/2005 force field. Our approach combines topology-based structural transformation coordinates, state-of-the-art free-energy calculation methods, and extensive unbiased molecular dynamics. All our diverse simulations cannot detect any barrier within the investigated timescales and system size, for a discontinuous transition between the LDL and HDL forms throughout the so-called "no man's land," until the onset of the solid, non-diffusive amorphous forms.

Evidence suggesting kinetic unfreezing of water mobility in two distinct processes in pressure-amorphized clathrate hydrates

Type II clathrate hydrates (CHs) with tetrahydrofuran (THF), cyclobutanone (CB) or 1,3-dioxolane (DXL) guest molecules collapse to an amorphous state near 1 GPa on pressurization below 140 K. On subsequent heating in the 0.2-0.7 GPa range, thermal conductivity and heat capacity results of the homogeneous amorphous solid show two glass transitions, first a thermally weak glass transition, GT1, near 130 K; thereafter a thermally strong glass transition, GT2, which implies a transformation to an ultraviscous liquid on heating. Here we compare the GTs of normal and deuterated samples and samples with different guest molecules. The results show that GT1 and GT2 are unaffected by deuteration of the THF guest and exchange of THF with CB or DXL, whereas the glass transition temperatures (T_gs) shift to higher temperatures on deuteration of water; T_g of GT2 increases by 2.5 K. These results imply that both GTs are associated with the water network. This is corroborated by the fact that GT2 is detected only in the state which is the amorphized CH's counterpart of expanded high density amorphous ice. The results suggest a rare transition sequence of an orientational glass transition followed by a glass to liquid transition, i.e., kinetic unfreezing of H₂O reorientational and translational mobility in two distinct processes.

Steady-like topology of the dynamical hydrogen bond network in supercooled water

We investigate the link between topology of the hydrogen bond network (HBN) and large-scale density fluctuations in water from ambient conditions to the glassy state. We observe a transition from a temperature-dependent topology at high temperatures, to a steady-like topology below the Widom temperature TW ∼ 220 K signaling the fragile-to-strong crossover and the maximum in structural fluctuations. As a consequence of the steady topology, the network suppresses large-scale density fluctuations much more efficiently than at higher temperatures. Below TW , the contribution of coordination defects of the kind A 2 D 1 (two acceptors and one donor) to the kinetics of the HBN becomes progressively more pronounced, suggesting that A 2 D 1 configurations may represent the main source of dynamical heterogeneities. Below the vitrification temperature, the freezing of rotational and translational degrees of freedom allow for an enhanced suppression of large-scale density fluctuations and the sample reaches the edges of nearly hyperuniformity. The formed network still hosts coordination defects, hence implying that nearly hyperuniformity goes beyond the classical continuous random network paradigm of tetrahedral networks and can emerge in scenarios much more complex than previously assumed. Our results unveil a hitherto undisclosed link between network topology and properties of water essential for better understanding water's rich and complex nature. Beyond implications for water, our findings pave the way to a better understanding of the physics of supercooled liquids and disordered hyperuniform networks at large.

Guiding epitaxial crystallization of amorphous solids at the nanoscale: interfaces, stress, and precrystalline order

The crystallization of amorphous solids impacts fields ranging from inorganic crystal growth to biophysics. Promoting or inhibiting nanoscale epitaxial crystallization and selecting its final products underpins applications in cryopreservation, semiconductor devices, oxide electronics, quantum electronics, structural and functional ceramics, and advanced glasses. As precursors for crystallization, amorphous solids are distinguished from liquids and gases by the comparatively long relaxation times for perturbations of the mechanical stress and for variations in composition or bonding. These factors allow experimentally controllable parameters to influence crystallization processes and to drive materials towards specific outcomes. For example, amorphous precursors can be employed to form crystalline phases, such as polymorphs of Al2O3, VO2, and other complex oxides, that are not readily accessible via crystallization from a liquid or through vapor-phase epitaxy. Crystallization of amorphous solids can further be guided to produce a desired polymorph, nanoscale shape, microstructure, and orientation of the resulting crystals. These effects can enable advances in applications in electronics, magnetic devices, optics, and catalysis. Directions for the future development of the chemical physics of crystallization from amorphous solids can be drawn from the impact of structurally complex and non-equilibrium atomic arrangements in liquids and the atomic-scale structure of liquid-solid interfaces.

Anomalous temperature dependence of the experimental x-ray structure factor of supercooled water

The structural changes of water upon deep supercooling were studied through wide-angle x-ray scattering at SwissFEL. The experimental setup had a momentum transfer range of 4.5 Å^-1, which covered the principal doublet of the x-ray structure factor of water. The oxygen-oxygen structure factor was obtained for temperatures down to 228.5 ± 0.6 K. Similar to previous studies, the second diffraction peak increased strongly in amplitude as the structural change accelerated toward a local tetrahedral structure upon deep supercooling. We also observed an anomalous trend for the second peak position of the oxygen-oxygen structure factor (q₂). We found that q₂ exhibits an unprecedented positive partial derivative with respect to temperature for temperatures below 236 K. Based on Fourier inversion of our experimental data combined with reference data, we propose that the anomalous q₂ shift originates from that a repeat spacing in the tetrahedral network, associated with all peaks in the oxygen-oxygen pair-correlation function, gives rise to a less dense local ordering that resembles that of low-density amorphous ice. The findings are consistent with that liquid water consists of a pentamer-based hydrogen-bonded network with low density upon deep supercooling.

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.

Superionicity, disorder, and bandgap closure in dense hydrogen chloride

Hydrogen bond networks play a crucial role in biomolecules and molecular materials such as ices. How these networks react to pressure directs their properties at extreme conditions. We have studied one of the simplest hydrogen bond formers, hydrogen chloride, from crystallization to metallization, covering a pressure range of more than 2.5 million atmospheres. Following hydrogen bond symmetrization, we identify a previously unknown phase by the appearance of new Raman modes and changes to x-ray diffraction patterns that contradict previous predictions. On further compression, a broad Raman band supersedes the well-defined excitations of phase V, despite retaining a crystalline chlorine substructure. We propose that this mode has its origin in proton (H+) mobility and disorder. Above 100 GPa, the optical bandgap closes linearly with extrapolated metallization at 240(10) GPa. Our findings suggest that proton dynamics can drive changes in these networks even at very high densities.

Structure of ice confined in silica nanopores

Observed anomalous thermodynamic properties of confined water such as deviations in the melting point and freezing point motivate the determination of the structure of confined water as a function of pore size and temperature. In this study, we investigate the dynamic evolution of the structure of confined ice in SBA-15 porous materials with pore diameters of 4 nm, 6 nm, and 8 nm at temperatures ranging from 183 K to 300 K using in operando Wide-Angle X-Ray Scattering (WAXS) measurements, X-Ray Partial Distribution Function (PDF) measurements, and classical Molecular Dynamics (MD) simulations. Formation of hexagonal ice structures is noted in all the three pore sizes. In silica nanopores with diameters of 4 nm, cubic ice formation is noted in addition to hexagonal ice. Longer lasting hydrogen bonds and longer residence times of the water molecules in the first coordination shell contribute to observed crystalline organization of ice in confinement. Self-diffusion coefficients of confined liquid water, predicted from classical MD simulations, are four orders of magnitude higher compared to ice formed in confinement. These experimental and simulation results provide comprehensive insights underlying the organization of confined water and ice in silica nanopores and the underlying physico-chemical interactions that contribute to the observed structures.

Crystal Structure of Symmetric Ice X in H2O-H2 and H2O-He under Pressure

Ice VII and ice X are the two most dominant phases, stable over a large pressure range between 2 and 150 GPa and made of fundamentally different chemical bonding. Yet, the two ice phases share a similar bcc-based crystal structure and lattice constants, resulting in a challenge to discern the crystal structure of ice VII and ice X. Here, we present well-resolved X-ray diffraction data of H₂O in quasi-hydrostatic H₂ and He pressure media, clearly resolving the two ice phases to 130 GPa and the dissociative nature of ice VII to X transition occurring at 20-50 GPa in H₂O-H₂ and 60-70 GPa in H₂O-He. The present diffraction data permits, for the first time, the accurate determination of the bulk moduli B₀ of 225 (or 228) GPa for ice X and 6.2 (or 4.5) GPa for ice VII, in H₂O-H₂ (or H₂O-He), which can provide new constraints for Giant planetary models.

Temperature-dependent kinetic pathways featuring distinctive thermal-activation mechanisms in structural evolution of ice VII

Ice amorphization, low- to high-density amorphous (LDA-HDA) transition, as well as (re)crystallization in ice, under compression have been studied extensively due to their fundamental importance in materials science and polyamorphism. However, the nature of the multiple-step "reverse" transformation from metastable high-pressure ice to the stable crystalline form under reduced pressure is not well understood. Here, we characterize the rate and temperature dependence of the structural evolution from ice VII to ice I recovered at low pressure (∼5 mTorr) using in situ time-resolved X-ray diffraction. Unlike previously reported ice VII (or ice VIII)→LDA→ice I transitions, we reveal three temperature-dependent successive transformations: conversion of ice VII into HDA, followed by HDA-to-LDA transition, and then crystallization of LDA into ice I. Significantly, the temperature-dependent characteristic times indicate distinctive thermal activation mechanisms above and below 110-115 K for both ice VIII-to-HDA and HDA-to-LDA transitions. Large-scale molecular-dynamics calculations show that the structural evolution from HDA to LDA is continuous and involves substantial movements of the water molecules at the nanoscale. The results provide a perspective on the interrelationship of polyamorphism and unravel its underpinning complexities in shaping ice-transition kinetic pathways.

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.

Is water one liquid or two?

The idea that water is a mixture of two distinct states is analyzed in some detail. It is shown that the known compressibility of water is in fact sufficiently small that for a volume of water of size 1 nm³, the density fluctuations are of order 4% of the average density. This is much smaller than the ≈25% density fluctuations that would be required for significant regions of high and low density water to occur on this volume scale. It is also pointed out that the density fluctuations in water are, if anything, smaller than those that occur in other common liquids which do not have the anomalous properties of water. It is shown that if the distribution of density fluctuations is unimodal, the system is in the one-phase region, and if bimodal, it is in the two-phase region. None of the liquid or amorphous phases of water explored in this work give any sign of being in the two-phase region. Existing neutron and X-ray scattering data on water in the amorphous phases, and in the stable liquid phases as a function pressure and temperature, are subject to a new set of empirical potential structure refinement simulations. These simulations are interrogated for their configurational entropy, using a spherical harmonic reconstruction of the full orientational pair correlation function. It is shown that the excess pair entropy derived from this function, plus the known perfect gas contributions, give a reasonable account of the total entropy of water, within the likely errors. This estimated entropy follows the expected declining trend with decreasing temperature. Evidence that higher density water will have higher entropy than lower density water emerges, in accordance with what is expected from the negative thermal expansion coefficient of water at low temperatures. However, this entropy increase is not large and goes through a maximum before declining at yet higher densities and pressures, in a manner reminiscent of what has been previously observed in the diffusion coefficient as a function of pressure. There is no evidence that ambient water can be regarded as patches of high density, high entropy and low density, low entropy liquid, as some have claimed, since high density water has a similar entropy to low density water. There is some evidence that the distinction between these two states will become more pronounced as the temperature is lowered. Extensive discussion of the use of order parameters to describe water structure is given, and it is pointed out that these indices generally cannot be used to infer two-state behavior.

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'.

New insights into decomposition characteristics of nanoscale methane hydrate below the ice point

In this paper, molecular dynamics simulation was used to study the decomposition process of nanoscale methane hydrate at 1 atm and 227 K. The results predict that methane hydrate decomposes into supercooled water (SCW) and methane gas and the resulting SCW turns into very high density amorphous ice (VHDA). The density of the VHDA is as high as 1.2-1.4 g cm^-3. The X-ray diffraction phase analysis showed that VHDA has a broad peak at 2θ of around 30°. The VHDA encapsulates the methane hydrate crystal so that the crystal can survive for a long time. The dissolved gas from the hydrate melt cannot escape out of the VHDA in a short time. The simulation results reveal new molecular insights into the decomposition behaviour of nanoscale methane hydrate below the ice point.

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.

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.

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.

Kinetically Controlled Two-Step Amorphization and Amorphous-Amorphous Transition in Ice

We report the results of in situ structural characterization of the amorphization of crystalline ice Ih under compression and the relaxation of high-density amorphous (HDA) ice under decompression at temperatures between 96 and 160 K by synchrotron x-ray diffraction. The results show that ice Ih transforms to an intermediate crystalline phase at 100 K prior to complete amorphization, which is supported by molecular dynamics calculations. The phase transition pathways show clear temperature dependence: direct amorphization without an intermediate phase is observed at 133 K, while at 145 K a direct Ih-to-IX transformation is observed; decompression of HDA shows a transition to low-density amorphous ice at 96 K and ∼1  Pa, to ice Ic at 135 K and to ice IX at 145 K. These observations show that the amorphization of compressed ice Ih and the recrystallization of decompressed HDA are strongly dependent on temperature and controlled by kinetic barriers. Pressure-induced amorphous ice is an intermediate state in the phase transition from the connected H-bond water network in low pressure ices to the independent and interpenetrating H-bond network of high-pressure ices.

Evidence of low-density and high-density liquid phases and isochore end point for water confined to carbon nanotube

Possible transition between two phases of supercooled liquid water, namely the low- and high-density liquid water, has been only predicted to occur below 230 K from molecular dynamics (MD) simulation. However, such a phase transition cannot be detected in the laboratory because of the so-called "no-man's land" under deeply supercooled condition, where only crystalline ices have been observed. Here, we show MD simulation evidence that, inside an isolated carbon nanotube (CNT) with a diameter of 1.25 nm, both low- and high-density liquid water states can be detected near ambient temperature and above ambient pressure. In the temperature-pressure phase diagram, the low- and high-density liquid water phases are separated by the hexagonal ice nanotube (hINT) phase, and the melting line terminates at the isochore end point near 292 K because of the retracting melting line from 292 to 278 K. Beyond the isochore end point (292 K), low- and high-density liquid becomes indistinguishable. When the pressure is increased from 10 to 600 MPa along the 280-K isotherm, we observe that water inside the 1.25-nm-diameter CNT can undergo low-density liquid to hINT to high-density liquid reentrant first-order transitions.

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.

Evolution of Hydrogen Dynamics in Amorphous Ice with Density

The single-particle dynamics of hydrogen atoms in several of the amorphous ices are reported using a combination of deep inelastic neutron scattering (DINS) and inelastic neutron scattering (INS). The mean kinetic energies of the hydrogen nuclei are found to increase with increasing density, indicating the weakening of hydrogen bonds as well as a trend toward steeper and more harmonic hydrogen vibrational potential energy surfaces. DINS shows much more pronounced changes in the O-H stretching component of the mean kinetic energy going from low- to high-density amorphous ices than indicated by INS and Raman spectroscopy. This highlights the power of the DINS technique to retrieve accurate ground-state kinetic energies beyond the harmonic approximation. In a novel approach, we use information from DINS and INS to determine the anharmonicity constants of the O-H stretching modes. Furthermore, our experimental kinetic energies will serve as important benchmark values for path-integral Monte Carlo simulations.

Structural characterization of eutectic aqueous NaCl solutions under variable temperature and pressure conditions

The structure of amorphous NaCl solutions produced by fast quenching is studied as a function of pressure, up to 4 GPa, by combined neutron diffraction experiments and classical molecular dynamics simulations. Similarly to LiCl solutions the system amorphizes at ambient pressure in a dense phase structurally similar to the e-HDA phase in pure water. The measurement of the static structure factor as a function of pressure allowed us to validate a new polarizable force field developed by Tazi et al., 2012, never tested under non-ambient conditions. We infer from simulations that the hydration shells of Na(+) cations form well defined octahedra composed of both H2O molecules and Cl(-) anions at low pressure. These octahedra are gradually broken by the seventh neighbour moving into the shell of first neighbours yielding an irregular geometry. In contrast to LiCl solutions and pure water, the system does not show a polyamorphic transition under pressure. This confirms that the existence of polyamorphism relies on the tetrahedral structure of water molecules, which is broken here.

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.

Connecting the Water Phase Diagram to the Metastable Domain: High-Pressure Studies in the Supercooled Regime

Pressure is extremely efficient to tune intermolecular interactions, allowing the study of the mechanisms regulating, at the molecular level, the structure and dynamics of condensed phases. Among the simplest molecules, water represents in many respects a mystery despite its primary role in ruling most of the biological, physical, and chemical processes occurring in nature. Here we report a careful characterization of the dynamic regime change associated with low-density and high-density forms of liquid water by measuring the line shape of the OD stretching mode of HOD in liquid water along different isotherms as a function of pressure. Remarkably, the high-pressure studies have been here extended down to 240 K, well inside the supercooled regime. Supported by molecular dynamics simulations, a correlation among amorphous and crystalline solids and the two different liquid water forms is attempted to provide a unified picture of the metastable and thermodynamic regimes of water.

Small-angle neutron scattering study of micropore collapse in amorphous solid water

Amorphous solid water (ASW) is shown to undergo a micropore collapse from cylindrical pores (3D) to lamellae (2D) at >120 K using small-angle neutron scattering.

Radical re-appraisal of water structure in hydrophilic confinement

The structure of water confined in MCM41 silica cylindrical pores is studied to determine whether confined water is simply a version of the bulk liquid which can be substantially supercooled without crystallisation. A combination of total neutron scattering from the porous silica, both wet and dry, and computer simulation using a realistic model of the scattering substrate is used. The water in the pore is divided into three regions: core, interfacial and overlap. The average local densities of water in these simulations are found to be about 20% lower than bulk water density, while the density in the core region is below, but closer to, the bulk density. There is a decrease in both local and core densities when the temperature is lowered from 298 K to 210 K. The radical proposal is made here that water in hydrophilic confinement is under significant tension, around -100 MPa, inside the pore.

Studies of cavitation and ice nucleation in 'doubly-metastable' water: time-lapse photography and neutron diffraction

High-speed photographic studies and neutron diffraction measurements have been made of water under tension in a Berthelot tube. Liquid water was cooled below the normal ice-nucleation temperature and was in a doubly-metastable state prior to a collapse of the liquid state. This transition was accompanied by an exothermic heat release corresponding with the rapid production of a solid phase nucleated by cavitation. Photographic techniques have been used to observe the phase transition over short time scales in which a solidification front is observed to propagate through the sample. Significantly, other images at a shorter time interval reveal the prior formation of cavitation bubbles at the beginning of the process. The ice-nucleation process is explained in terms of a mechanism involving hydrodynamically-induced changes in tension in supercooled water in the near vicinity of an expanding cavitation bubble. Previous explanations have attributed the nucleation of the solid phase to the production of high positive pressures. Corresponding results are presented which show the initial neutron diffraction pattern after ice-nucleation. The observed pattern does not exhibit the usual crystalline pattern of hexagonal ice [I(h)] that is formed under ambient conditions, but indicates the presence of other ice forms. The composite features can be attributed to a mixture of amorphous ice, ice-I(h)/I(c) and the high-pressure form, ice-III, and the diffraction pattern continues to evolve over a time period of about an hour.

Highly compressed two-dimensional form of water at ambient conditions

The structure of thin-film water on a BaF₂(111) surface under ambient conditions was studied using x-ray absorption spectroscopy from ambient to supercooled temperatures at relative humidity up to 95%. No hexagonal ice-like structure was observed in spite of the expected templating effect of the lattice-matched (111) surface. The oxygen K-edge x-ray absorption spectrum of liquid thin-film water on BaF₂ exhibits, at all temperatures, a strong resemblance to that of high-density phases for which the observed spectroscopic features correlate linearly with the density. Surprisingly, the highly compressed, high-density thin-film liquid water is found to be stable from ambient (300 K) to supercooled (259 K) temperatures, although a lower-density liquid would be expected at supercooled conditions. Molecular dynamics simulations indicate that the first layer water on BaF₂(111) is indeed in a unique local structure that resembles high-density water, with a strongly collapsed second coordination shell.