Negative Thermal Expansion in Cubic Ice: A Collective Quantum Effect of the hydrogen-bond network

Abstract

We report neutron powder diffraction measurements and path-integral molecular dynamics simulations of stacking-disorder-free cubic ice I$_c$, produced by topotactic degassing of C2 hydrogen hydrate. Across the cryogenic stability range, ice I$_c$ exhibits a density maximum near 70 K, closely matching that of hexagonal ice I$_h$ despite their different long-range stacking sequences. Negative thermal expansion in ice I is therefore not specific to hexagonal stacking, but arises from the shared open tetrahedral hydrogen-bond network. Simulations with the MB-pol potential quantitatively reproduce the experimental anomaly only when nuclear quantum effects are included. The density maximum coincides, within the temperature resolution, with maximal anisotropy of the proton quantum distribution. Neutron-derived displacement parameters independently reveal a strongly enhanced transverse proton displacement, while phonon calculations identify low-frequency transverse modes with the most negative Grüneisen parameters. Together, these results establish the negative thermal expansion of ice I as a collective quantum effect governed by nuclear statistics and the dynamics of the hydrogen-bond network.