PNNL

Abstract

Ice Ih is arguably the most important molecular crystal in nature, yet our understanding of its structural and dynamical properties is still incomplete. To explain the origin of two peaks in the hydrogen-bond-stretching region of the inelastic neutron scattering (INS) spectra, the existence of two types of hydrogen bonds with strengths differing by a factor of two was previously hypothesized. We present first-principles calculations based on diagrammatic many-body perturbation theory of the structures and vibrational spectra of ice Ih, which suggest that the observed spectral features arise from the directionality or anisotropy of the hydrogen-bond stretching vibrations rather than their vastly different force constants, disproving the previous hypothesis. Our calculations also reproduce the infrared and Raman spectra, the variation of INS spectra with deuterium concentration, and the anomaly of heat capacities at low temperatures, together rendering our calculations a paradigm for "crystals from first principles" as envisioned by Maddox.