Ice Ih anomalies: Thermal contraction, anomalous volume isotope effect, and pressure-induced amorphization
- Department of Chemistry, University of Illinois at Urbana-Champaign, 600 South Mathews Avenue, Urbana, Illinois 61801 (United States)
Ice Ih displays several anomalous thermodynamic properties such as thermal contraction at low temperatures, an anomalous volume isotope effect (VIE) rendering the volume of D{sub 2}O ice greater than that of H{sub 2}O ice, and a pressure-induced transition to the high-density amorphous (HDA) phase. Furthermore, the anomalous VIE increases with temperature, despite its quantum-mechanical origin. Here, embedded-fragment ab initio second-order many-body perturbation (MP2) theory in the quasiharmonic approximation (QHA) is applied to the Gibbs energy of an infinite, proton-disordered crystal of ice Ih at wide ranges of temperatures and pressures. The quantum effect of nuclei moving in anharmonic potentials is taken into account from first principles without any empirical or nonsystematic approximation to either the electronic or vibrational Hamiltonian. MP2 predicts quantitatively correctly the thermal contraction at low temperatures, which is confirmed to originate from the volume-contracting hydrogen-bond bending modes (acoustic phonons). It qualitatively reproduces (but underestimates) the thermal expansion at higher temperatures, caused by the volume-expanding hydrogen-bond stretching (and to a lesser extent librational) modes. The anomalous VIE is found to be the result of subtle cancellations among closely competing isotope effects on volume from all modes. Consequently, even ab initio MP2 with the aug-cc-pVDZ and aug-cc-pVTZ basis sets has difficulty reproducing this anomaly, yielding qualitatively varied predictions of the sign of the VIE depending on such computational details as the choice of the embedding field. However, the temperature growth of the anomalous VIE is reproduced robustly and is ascribed to the librational modes. These solid-state MP2 calculations, as well as MP2 Born–Oppenheimer molecular dynamics, find a volume collapse and a loss of symmetry and long-range order in ice Ih upon pressure loading of 2.35 GPa or higher. Concomitantly, rapid softening of acoustic phonons is observed starting around 2 GPa. They constitute a computational detection of a mechanical instability in ice Ih and the resulting pressure-induced amorphization to HDA.
- OSTI ID:
- 22657824
- Journal Information:
- Journal of Chemical Physics, Vol. 144, Issue 20; Other Information: (c) 2016 Author(s); Country of input: International Atomic Energy Agency (IAEA); ISSN 0021-9606
- Country of Publication:
- United States
- Language:
- English
Similar Records
Energy benchmarks for water clusters and ice structures from an embedded many-body expansion
Final Scientific/Technical Report
Related Subjects
ORGANIC
PHYSICAL AND ANALYTICAL CHEMISTRY
AMORPHOUS STATE
CONTRACTION
CRYSTAL GROWTH
EXPERIMENTAL DATA
HEAVY WATER
ICE
ISOTOPE EFFECTS
MOLECULAR DYNAMICS METHOD
PERTURBATION THEORY
PRESSURE RANGE GIGA PA
TEMPERATURE RANGE 0065-0273 K
TEMPERATURE RANGE 0400-1000 K
THERMAL EXPANSION
THERMODYNAMIC PROPERTIES