Autoionization in liquid water

Ab initio spectroscopy and ionic conductivity of water under Earth mantle conditions

The phase diagram of water at extreme conditions plays a critical role in Earth and planetary science, yet remains poorly understood. Here we report a first-principles investigation of the liquid at high temperature, between 11 GPa and 20 GPa-a region where numerous controversial results have been reported over the past three decades. Our results are consistent with the recent estimates of the water melting line below 1,000 K and show that on the 1,000-K isotherm the liquid is rapidly dissociating and recombining through a bimolecular mechanism. We found that short-lived ionic species act as charge carriers, giving rise to an ionic conductivity that at 11 GPa and 20 GPa is six and seven orders of magnitude larger, respectively, than at ambient conditions. Conductivity calculations were performed entirely from first principles, with no a priori assumptions on the nature of charge carriers. Despite frequent dissociative events, we observed that hydrogen bonding persists at high pressure, up to at least 20 GPa. Our computed Raman spectra, which are in excellent agreement with experiment, show no distinctive signatures of the hydronium and hydroxide ions present in our simulations. Instead, we found that infrared spectra are sensitive probes of molecular dissociation, exhibiting a broad band below the OH stretching mode ascribable to vibrations of complex ions.

Inverse Temperature Dependence of Nuclear Quantum Effects in DNA Base Pairs

Despite the inherently quantum mechanical nature of hydrogen bonding, it is unclear how nuclear quantum effects (NQEs) alter the strengths of hydrogen bonds. With this in mind, we use ab initio path integral molecular dynamics to determine the absolute contribution of NQEs to the binding in DNA base pair complexes, arguably the most important hydrogen-bonded systems of all. We find that depending on the temperature, NQEs can either strengthen or weaken the binding within the hydrogen-bonded complexes. As a somewhat counterintuitive consequence, NQEs can have a smaller impact on hydrogen bond strengths at cryogenic temperatures than at room temperature. We rationalize this in terms of a competition of NQEs between low-frequency and high-frequency vibrational modes. Extending this idea, we also propose a simple model to predict the temperature dependence of NQEs on hydrogen bond strengths in general.

Bonding in the superionic phase of water

The predicted superionic phase of water is investigated via ab initio molecular dynamics at densities of 2.0--3.0 g/cc (34-115 GPa) along the 2000 K isotherm. We find that extremely rapid (superionic) diffusion of protons occurs in a fluid phase at pressures between 34 and 58 GPa. A transition to a stable body-centered cubic O lattice with superionic proton conductivity is observed between 70 and 75 GPa, a much higher pressure than suggested in prior work. We find that all molecular species at pressures greater than 75 GPa are too short lived to be classified as bound states. Up to 95 GPa, we find a solid superionic phase characterized by covalent O-H bonding. Above 95 GPa, a transient network phase is found characterized by symmetric O-H hydrogen bonding with nearly 50% covalent character. In addition, we describe a metastable superionic phase with quenched O disorder.