Classical and quantum ordering of protons in cold solid hydrogen under megabar pressures

TL;DR

This study combines advanced theoretical methods to explore classical and quantum proton ordering in cold solid hydrogen under megabar pressures, revealing quantum effects in phase transitions and the importance of van der Waals forces.

Contribution

It provides new insights into the quantum and classical nature of phase transitions in high-pressure solid hydrogen, highlighting the role of quantum nuclear motion and van der Waals interactions.

Findings

Quantum effects strongly influence the I-II phase transition.

A P21c structure explains phase II vibrational spectra.

Van der Waals forces improve theoretical agreement with experiments.

Abstract

A combination of state-of-the-art theoretical methods has been used to obtain an atomic-level picture of classical and quantum ordering of protons in cold high-pressure solid hydro-gen. We focus mostly on phases II and III of hydrogen, exploring the effects of quantum nuclear motion on certain features of these phases (through a number of ab initio path-integral molecular dynamics (PIMD) simulations at particular points on the phase diagram). We also examine the importance of van der Waals forces in this system by performing calculations using the optB88-vdW density functional, which accounts for non-local correlations. Our calculations reveal that the transition between phases I and II is strongly quantum in nature, resulting from a competition between anisotropic inter-molecular interactions that restrict molecular rotation and thermal plus quantum fluctuations of the nuclear positions that facilitate it. The transition from phase II to III is more classical because quantum nuclear motion plays only a secondary role and the transition is determined primarily by the underlying potential energy surface. A structure of P21c symmetry with 24 atoms in the primitive unit cell is found to be stable when anharmonic quantum nuclear vibrational motion is included at finite temperatures using the PIMD method. This structure gives a good account of the infra-red (IR) and Raman vibron frequencies of phase II. We find ad-ditional support for C2/c structure as a strong candidate for phase III, since it remains transparent up to 300 GPa, even when quantum nuclear effects are included. Finally, we find that accounting for van der Waals forces improves the agreement between experiment and theory for the parts of the phase diagram considered, when compared to previous work which employed the widely-used Perdew-Burke-Ernzerhof (PBE) exchange-correlation functional.