The structural validity of various thermodynamical models of supercooled water.

The thermodynamic response functions of water exhibit an anomalous increase upon cooling that becomes strongly amplified in the deeply supercooled regime due to structural fluctuations between disordered and tetrahedral local structures. Here, we compare structural data from recent x-ray laser scattering measurements of water at 1 bar and temperatures down to 227 K with structural properties computed for several different water models using molecular dynamics simulations. Based on this comparison, we critically evaluate four different thermodynamic scenarios that have been invoked to explain the unusual behavior of water. The critical point-free model predicts small variations in the tetrahedrality with decreasing temperature, followed by a stepwise change at the liquid-liquid transition around 228 K at ambient pressure. This scenario is not consistent with the experimental data that instead show a smooth and accelerated variation in structure from 320 to 227 K. Both the singularity-free model and ice coarsening hypothesis give trends that indirectly indicate an increase in tetrahedral structure with temperature that is too weak to be consistent with experiment. A model that includes an apparent divergent point (ADP) at high positive pressure, however, predicts structural development consistent with our experimental measurements. The terminology ADP, instead of the commonly used liquid-liquid critical point, is more general in that it focuses on the growing fluctuations, whether or not they result in true criticality. Extrapolating this model beyond the experimental data, we estimate that an ADP in real water may lie around 1500 ± 250 bars and 190 ± 6 K.

Ultrafast X-ray probing of water structure below the homogeneous ice nucleation temperature.

Water has a number of anomalous physical properties, and some of these become drastically enhanced on supercooling below the freezing point. Particular interest has focused on thermodynamic response functions that can be described using a normal component and an anomalous component that seems to diverge at about 228 kelvin (refs 1-3). This has prompted debate about conflicting theories that aim to explain many of the anomalous thermodynamic properties of water. One popular theory attributes the divergence to a phase transition between two forms of liquid water occurring in the 'no man's land' that lies below the homogeneous ice nucleation temperature (TH) at approximately 232 kelvin and above about 160 kelvin, and where rapid ice crystallization has prevented any measurements of the bulk liquid phase. In fact, the reliable determination of the structure of liquid water typically requires temperatures above about 250 kelvin. Water crystallization has been inhibited by using nanoconfinement, nanodroplets and association with biomolecules to give liquid samples at temperatures below TH, but such measurements rely on nanoscopic volumes of water where the interaction with the confining surfaces makes the relevance to bulk water unclear. Here we demonstrate that femtosecond X-ray laser pulses can be used to probe the structure of liquid water in micrometre-sized droplets that have been evaporatively cooled below TH. We find experimental evidence for the existence of metastable bulk liquid water down to temperatures of 227(-1)(+2) kelvin in the previously largely unexplored no man's land. We observe a continuous and accelerating increase in structural ordering on supercooling to approximately 229 kelvin, where the number of droplets containing ice crystals increases rapidly. But a few droplets remain liquid for about a millisecond even at this temperature. The hope now is that these observations and our detailed structural data will help identify those theories that best describe and explain the behaviour of water.

Communication: ab initio simulations of hydrogen-bonded ferroelectrics: collective tunneling and the origin of geometrical isotope effects.

Ab initio simulations that account for nuclear quantum effects have been used to examine the order-disorder transition in squaric acid, a prototypical H-bonded antiferroelectric crystal. Our simulations reproduce the >100 K difference in transition temperature observed upon deuteration as well as the strong geometrical isotope effect observed on intermolecular separations within the crystal. We find that collective transfer of protons along the H-bonding chains - facilitated by quantum mechanical tunneling - is critical to the order-disorder transition and the geometrical isotope effect. This sheds light on the origin of isotope effects and the importance of tunneling in squaric acid which likely extends to other H-bonded ferroelectrics.

A transferable H2O interaction potential based on a single center multipole expansion: SCME.

A transferable potential energy function for describing the interaction between water molecules is presented. The electrostatic interaction is described rigorously using a multipole expansion. Only one expansion center is used per molecule to avoid the introduction of monopoles. This single center approach turns out to converge and give close agreement with ab initio calculations when carried out up to and including the hexadecapole. Both dipole and quadrupole polarizability are included. All parameters in the electrostatic interaction as well as the dispersion interaction are taken from ab initio calculations or experimental measurements of a single water molecule. The repulsive part of the interaction is parametrized to fit ab initio calculations of small water clusters and experimental measurements of ice Ih. The parametrized potential function was then used to simulate liquid water and the results agree well with experiment, even better than simulations using some of the point charge potentials fitted to liquid water. The evaluation of the new interaction potential for condensed phases is fast because point charges are not present and the interaction can, to a good approximation, be truncated at a finite range.

Wide-angle X-ray diffraction and molecular dynamics study of medium-range order in ambient and hot water.

We have developed wide-angle X-ray diffraction measurements with high energy-resolution and accuracy to study water structure at three different temperatures (7, 25 and 66 °C) under normal pressure. Using a spherically curved Ge crystal an energy resolution better than 15 eV has been achieved which eliminates influence from Compton scattering. The high quality of the data allows for a reliable Fourier transform of the experimental data resolving shell structure out to ~12 Å, i.e. 5 hydration shells. Large-scale molecular dynamics (MD) simulations using the TIP4P/2005 force-field reproduce excellently the experimental shell-structure in the range 4-12 Å although less agreement is seen for the first peak in the intermolecular pair-correlation function (PCF). The Shiratani-Sasai Local Structure Index [J. Chem. Phys. 104, 7671 (1996)] identifies a tetrahedral minority giving the intermediate-range oscillations in the O-O PCF and a disordered majority providing a more featureless background in this range. The current study supports the proposal that the structure of liquid water, even at high temperatures, can be described in terms of a two-state fluctuation model involving local structures related to the high-density and low-density forms of liquid water postulated in the liquid-liquid phase transition hypothesis.

Spatially inhomogeneous bimodal inherent structure of simulated liquid water.

In the supercooled regime at elevated pressure two forms of liquid water, high-density (HDL) and low-density (LDL), have been proposed to be separated by a coexistence line ending at a critical point, but a connection to water at ambient conditions has been lacking. Here we perform large-scale molecular dynamics simulations and demonstrate that the underlying potential energy surface gives a strictly bimodal characterization of the molecules at all temperatures and pressures, including the biologically and technologically important ambient regime, as spatially inhomogeneous either LDL- or HDL-like with a 3 : 1 predominance for HDL under ambient conditions. The Widom line in the supercooled regime, where maximal structural fluctuations take place, coincides with a 1 : 1 distribution. Although our results are based on molecular dynamics force-field simulations the close agreement with recent analyses of experimental X-ray spectroscopy and scattering data indicates a unified description also of real liquid water covering supercooled to ambient conditions.

Enhanced small-angle scattering connected to the Widom line in simulations of supercooled water.

We present extensive simulations on the TIP4P∕2005 water model showing significantly enhanced small-angle scattering (SAS) in the supercooled regime. The SAS is related to the presence of a Widom line (T(W)) characterized by maxima in thermodynamic response functions and Ornstein-Zernike correlation length. Recent experimental small-angle x-ray scattering data [Huang et al., J. Chem. Phys. 133, 134504 (2010)] are excellently reproduced, albeit with an increasing temperature offset at lower temperatures. Assuming the same origin of the SAS in experiment and model this suggests the existence of a Widom line also in real supercooled water. Simulations performed at 1000 bar show an increased abruptness of a crossover from dominating high-density (HDL) to dominating low-density (LDL) liquid and strongly enhanced SAS associated with crossing T(W), consistent with a recent determination of the critical pressure of TIP4P∕2005 at 1350 bar. Furthermore, good agreement with experimental isothermal compressibilities at 1000, 1500, and 2000 bar shows that the high pressure supercooled thermodynamic behavior of water is well described by TIP4P∕2005. Analysis of the tetrahedrality parameter Q reveals that the HDL-LDL structural transition is very sharp at 1000 bar, and that structural fluctuations become strongly coupled to density fluctuations upon approaching T(W). Furthermore, the tetrahedrality distribution becomes bimodal at ambient temperatures, an observation that possibly provides a link between HDL-LDL fluctuations and the structural bimodality in liquid water indicated by x-ray spectroscopic techniques. Computed x-ray absorption spectra are indeed found to show sensitivity to the tetrahedrality parameter.

Increasing correlation length in bulk supercooled H2O, D2O, and NaCl solution determined from small angle x-ray scattering.

Using small angle x-ray scattering, we find that the correlation length of bulk liquid water shows a steep increase as temperature decreases at subzero temperatures (supercooling) and that it can, similar to the thermodynamic response functions, be fitted to a power law. This indicates that the anomalous properties of water are attributable to fluctuations between low- and high-density regions with rapidly growing average size upon supercooling. The substitution of H(2)O with D(2)O, as well as the addition of NaCl salt, leads to substantial changes of the power law behavior of the correlation length. Our results are consistent with the proposed existence of a liquid-liquid critical point in the deeply supercooled region but do not exclude a singularity-free model.

The inhomogeneous structure of water at ambient conditions.

Small-angle X-ray scattering (SAXS) is used to demonstrate the presence of density fluctuations in ambient water on a physical length-scale of approximately 1 nm; this is retained with decreasing temperature while the magnitude is enhanced. In contrast, the magnitude of fluctuations in a normal liquid, such as CCl(4), exhibits no enhancement with decreasing temperature, as is also the case for water from molecular dynamics simulations under ambient conditions. Based on X-ray emission spectroscopy and X-ray Raman scattering data we propose that the density difference contrast in SAXS is due to fluctuations between tetrahedral-like and hydrogen-bond distorted structures related to, respectively, low and high density water. We combine our experimental observations to propose a model of water as a temperature-dependent, fluctuating equilibrium between the two types of local structures driven by incommensurate requirements for minimizing enthalpy (strong near-tetrahedral hydrogen-bonds) and maximizing entropy (nondirectional H-bonds and disorder). The present results provide experimental evidence that the extreme differences anticipated in the hydrogen-bonding environment in the deeply supercooled regime surprisingly remain in bulk water even at conditions ranging from ambient up to close to the boiling point.