Second-order many-body perturbation study of ice Ih

Ice Ih is arguably the most important molecular crystal in nature, yet our understanding of its structural and dynamical properties is still far from complete. We present embedded-fragment calculations of the structures and vibrational spectra of the three-dimensional, proton-disordered phase of ice Ih performed at the level of second-order many-body perturbation theory with a basis-set superposition error correction. Our calculations address previous controversies such as the one related to the O-H bond length as well as the existence of two types of hydrogen bonds with strengths differing by a factor of two. For the latter, our calculations suggest that the observed spectral features arise from the directionality or the anisotropy of collective hydrogen-bond stretching vibrations rather than the previously suggested vastly different force constants. We also report a capability to efficiently compute infrared and Raman intensities of a periodic solid. Our approach reproduces the infrared and Raman spectra, the variation of inelastic neutron scattering spectra with deuterium concentration, and the anomaly of heat capacities at low temperatures for ice Ih.