Dual quantum locking: Dynamic coupling of hydrogen and water sublattices in hydrogen filled ice

Hydrogen hydrates (HH) are a unique class of materials composed of hydrogen molecules confined within crystalline water frameworks. Among their multiple phases, the filled ice structures, particularly the cubic C2 phase, exhibit exceptionally strong host-guest interactions due to ultra-short H2-H2O distances and a 1:1 stoichiometry leading to two interpenetrated identical diamond-like sublattices, one comprised of water molecules, the other of hydrogen molecules. At high pressures, nuclear quantum effects involving both hydrogen molecules and the water lattice become dominant, giving rise to a dual-lattice quantum system. In this work, we explore the sequence of pressure- and temperature-driven phase transitions in HH, focusing on the interplay between molecular rotation, orientational ordering, lattice symmetry breaking and hydrogen bond symmetrization. Using a combination of computational modeling based on classical and path-integral molecular dynamics, quantum embedding, and high pressure experiments, including Raman spectroscopy and synchrotron X-ray diffraction at low temperatures and high pressures, we identify signatures of quantum-induced ordering and structural transformations in the C2 phase. Our findings reveal that orientational ordering in HH occurs at much lower pressures than in solid hydrogen, by inducing structural changes in the water network and enhancing the coupling of water and hydrogen dynamics. This work provides new insights into the quantum behavior of hydrogen under extreme mechanochemical confinement and establishes hydrogen-filled ices as a promising platform for the design of hydrogen-rich quantum materials.

💡 Research Summary

This paper investigates the dual‑lattice quantum behavior of hydrogen‑filled ice (HH), focusing on the cubic C2 phase where water (H₂O) and hydrogen (H₂) each form interpenetrating diamond‑like sublattices with a 1:1 stoichiometry. The authors combine classical molecular dynamics, path‑integral molecular dynamics (PIMD), quantum‑embedding electronic structure calculations, Raman spectroscopy, and synchrotron X‑ray diffraction (XRD) to map the pressure‑temperature phase diagram and to elucidate the microscopic mechanisms behind the observed structural transformations.

At moderate pressures (<30 GPa) and high temperatures, the H₂ sublattice behaves as a “quantum plastic crystal”: H₂ molecules occupy well‑defined translational sites on a diamond lattice but retain nearly free rotational motion due to quantum fluctuations. The water sublattice remains proton‑disordered but structurally constrained by ice rules, preserving overall cubic (Fd‑3m) symmetry.

When temperature is lowered or pressure is increased beyond ~30 GPa, the rotational freedom of H₂ is progressively quenched. The authors introduce an orientation factor S (0 = complete disorder, 1 = perfect alignment) derived from quantum‑embedded calculations. S rises sharply in the pressure range 27–35 GPa, indicating a collective alignment of H₂ molecules along the crystallographic c‑axis. Two distinct orientational regimes are identified: a herringbone‑like arrangement at intermediate pressures (tilt angle decreasing with pressure) and a nematic phase at higher pressures where H₂ molecules point directly along c.

Concomitantly, the water lattice undergoes proton‑symmetrization (the O–H–O hydrogen bond becomes symmetric) at ~26 GPa, a pressure significantly lower than that required for symmetrization in pure ice VII. This symmetrization shortens O–O distances, intensifies the H₂–H₂O interaction, and creates an anisotropic crystal field that biases the external potential experienced by H₂. The authors solve the Schrödinger equation for a single H₂ molecule in this field, showing that the angular part of the potential evolves from multiple minima (favoring herringbone) to a single minimum at θ = 0,π (favoring nematic order).

Raman measurements reveal pressure‑induced broadening and splitting of the S₀(0) and S₀(1) rotational bands, with a characteristic 4 meV shift and loss of degeneracy that mirrors the loss of rotational symmetry. XRD data show a clear cubic‑to‑tetragonal transition (I4₁/amd) manifested by peak splitting; the tetragonal distortion quantified as c/√2a grows from ~0.5 % at 2.45 GPa, 150 K to ~4.6 % at 48 GPa, 300 K. The transition is sharper under compression than cooling, indicating that pressure more directly drives the coupling between H₂ orientation and lattice deformation.

Overall, the study demonstrates a feedback loop: proton symmetrization in the water framework modifies the crystal field, which locks H₂ rotations, and the resulting H₂ orientational order exerts anisotropic stress that further stabilizes the tetragonal distortion. This dual‑lattice quantum locking occurs at pressures far lower than those needed for similar ordering in pure solid hydrogen, highlighting the crucial role of host‑guest coupling. The work positions hydrogen‑filled ices as a versatile platform for exploring strongly coupled quantum lattices, with potential implications for high‑pressure superconductivity, quantum information materials, and the broader field of quantum confined molecular solids.