Proton dynamics in ice VII at high pressures

Dimensionality dependence of the Kauzmann temperature: A case study using bulk and confined water

The Kauzmann temperature (T_K) of a supercooled liquid is defined as the temperature at which the liquid entropy becomes equal to that of the crystal. The excess entropy, the difference between liquid and crystal entropies, is routinely used as a measure of the configurational entropy, whose vanishing signals the thermodynamic glass transition. The existence of the thermodynamic glass transition is a widely studied subject, and of particular recent interest is the role of dimensionality in determining the presence of a glass transition at a finite temperature. The glass transition in water has been investigated intensely and is challenging as the experimental glass transition appears to occur at a temperature where the metastable liquid is strongly prone to crystallization and is not stable. To understand the dimensionality dependence of the Kauzmann temperature in water, we study computationally bulk water (three-dimensions), water confined in the slit pore of the graphene sheet (two-dimensions), and water confined in the pore of the carbon nanotube of chirality (11,11) having a diameter of 14.9 Å (one-dimension), which is the lowest diameter where amorphous water does not always crystallize into nanotube ice in the supercooled region. Using molecular dynamics simulations, we compute the entropy of water in bulk and under reduced dimensional nanoscale confinement to investigate the variation of the Kauzmann temperature with dimension. We obtain a value of T_K (133 K) for bulk water in good agreement with experiments [136 K (C. A. Angell, Science 319, 582-587 (2008) and K. Amann-Winkel et al., Proc. Natl. Acad. Sci. U. S. A. 110, 17720-17725 (2013)]. However, for confined water, in two-dimensions and one-dimension, we find that there is no finite temperature Kauzmann point (in other words, the Kauzmann temperature is 0 K). Analysis of the fluidicity factor, a measure of anharmonicity in the oscillation of normal modes, reveals that the Kauzmann temperature can also be computed from the difference in the fluidicity factor between amorphous and ice phases.

Hydrogen nuclear mean kinetic energy in water down the Mariana Trench: Competition of pressure and salinity

The Mariana Trench is one of the most famous and extreme environments on our planet. We report experimental values of the hydrogen nuclear mean kinetic energy in water samples at the same physical and chemical conditions than in the Challenger Deep within the Mariana Trench: a pressure of 1092 bars, a temperature of 1 °C, and a salinity of 35 g of salt per kg of water. Results were obtained by deep inelastic neutron scattering at the VESUVIO spectrometer at ISIS. We find that the effect of pressure is to increase the hydrogen nuclear mean kinetic energy with respect to ambient conditions, while ions in the solution have the opposite effect. These results confirm the recent state-of-the-art simulations of the nuclear hydrogen dynamics in water. The changes in the nuclear mean kinetic energy likely correspond to different isotopic fractionation values in the Challenger Deep compared to standard sea water.

Quantum behavior of water nano-confined in beryl

The proton mean kinetic energy, Ke(H), of water confined in nanocavities of beryl (Be₃Al₂Si₆O₁₈) at 5 K was obtained by simulating the partial vibrational density of states from density functional theory based first-principles calculations. The result, Ke(H) = 104.4 meV, is in remarkable agreement with the 5 K deep inelastic neutron scattering (DINS) measured value of 105 meV. This is in fact the first successful calculation that reproduces an anomalous DINS value regarding Ke(H) in nano-confined water. The calculation indicates that the vibrational states of the proton of the nano-confined water molecule distribute much differently than in ordinary H₂O phases, most probably due to coupling with lattice modes of the hosting beryl nano-cage. These findings may be viewed as a promising step towards the resolution of the DINS controversial measurements on other H₂O nano-confining systems, e.g., H₂O confined in single and double walled carbon nanotubes.

On the mean kinetic energy of the proton in strong hydrogen bonded systems

The mean atomic kinetic energies of the proton, Ke(H), and of the deuteron, Ke(D), were calculated in moderate and strongly hydrogen bonded (HB) systems, such as the ferro-electric crystals of the KDP type (XH2PO4, X = K, Cs, Rb, Tl), the DKDP (XD2PO4, X = K, Cs, Rb) type, and the X3H(SO4)2 superprotonic conductors (X = K, Rb). All calculations utilized the simulated partial phonon density of states, deduced from density functional theory based first-principle calculations and from empirical lattice dynamics simulations in which the Coulomb, short range, covalent, and van der Waals interactions were accounted for. The presently calculated Ke(H) values for the two systems were found to be in excellent agreement with published values obtained by deep inelastic neutron scattering measurements carried out using the VESUVIO instrument of the Rutherford Laboratory, UK. The Ke(H) values of the M3H(SO4)2 compounds, in which the hydrogen bonds are centro-symmetric, are much lower than those of the KDP type crystals, in direct consistency with the oxygen-oxygen distance ROO, being a measure of the HB strength.

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.