Insights into the structure of liquid water from nuclear quantum effects on density and compressibility of ice polymorphs

TL;DR

This paper investigates nuclear quantum effects on the density and compressibility of ice and water, revealing anomalous isotope effects and their dependence on temperature and structure, which explain observed behaviors in liquid water.

Contribution

It provides a detailed analysis of nuclear quantum effects on various ice phases and liquid water, highlighting the transition from anomalous to normal isotope effects with temperature and density.

Findings

Heavy ice has a larger volume than light ice due to quantum effects.

Anomalous isotope effect observed in bulk modulus of hexagonal ice and clathrate hydrate.

Transition from anomalous to normal isotope effects correlates with increasing density.

Abstract

Nuclear quantum effects lead to an anomalous shift of the volume of hexagonal ice; heavy ice has a larger volume than light ice. This anomaly in ice increases with temperature and persists in liquid water up to the boiling point. We study nuclear quantum effects on the density and compressibility of several ice-like structures and crystalline ice phases. By calculating the anisotropic contributions to the stain tensor, we analyze how the compressibility changes along different directions in hexagonal ice, and find that hexagonal ice is softer along the x-y plane than the z-direction. Furthermore, by performing ab initio density functional theory calculations with a van der Waals functional and with the quasiharmonic approximation, we find an anomalous isotope effect in the bulk modulus of hexagonal ice: heavy ice has a smaller bulk modulus than light ice. In agreement with the experiments, we also obtain an anomalous isotope effect for clathrate hydrate structure I. For the rest of the ice polymorphs, the isotope effect is: i) anomalous for ice IX, Ih, Ic, clathrate, and low density liquid-like amorphous ice; ii) normal at T=0 K and becomes anomalous with increasing temperature for ice IX, II, high density liquid-like amorphous ices, and ice XV; iii) normal for ice VIII up to the melting point. There is a transition from an anomalous isotope effect to a normal isotope effect for both the volume and bulk modulus, as the density (compressibility) of the structures increases (decreases). This result can explain the anomalous isotope effect in liquid water: as the compressibility decreases from melting point to the compressibility minimum temperature, the difference between the volumes of the heavy and light water rapidly decreases, but the effect stays anomalous up to the boiling temperature as the hydrogen bond network is never completely broken by fully filling all the interstitial sites.