We study the effect of temperature up to 1000K on the structure of dense molecular para-hydrogen and ortho-deuterium, using the path-integral Monte Carlo method. We find a structural phase transition from orientationally disordered hexagonal close packed (hcp) to an orthorhombic structure of Cmca symmetry before melting. The transition is basically induced by thermal fluctuations, but quantum fluctuations of protons (deuterons) are important in determining the transition temperature through effectively hardening the intermolecular interaction. We estimate the phase line between hcp and Cmca phases as well as the melting line of the Cmca solid.
The recent discovery of phase IV of solid hydrogen and deuterium consisting of two alternate layers of graphenelike three-molecule rings and unbound H2 molecules have generated great interests. However, vibrational nature of phase IV remains poorly understood. Here, we report a peculiar proton transfer and a simultaneous rotation of three molecule rings in graphenelike layers predicted by ab initio variable cell molecular dynamics simulations for phase IV of solid hydrogen and deuterium at pressure ranges of from 250 to 350 GPa and temperature range of from 300 to 500 K. This proton transfer is intimately related to the particular elongation of molecules in graphenelike layers, and it becomes more pronounced with increasing pressure at the course of larger elongation of molecules. As the consequence of proton transfer, hydrogen molecules in graphenelike layers are short lived and hydrogen vibration is strongly anharmonic. Our findings provide direct explanations on the observed abrupt increase of Raman width at the formation of phase IV and its large increase with pressure.
Quantum nuclear zero-point motions in solid H$_2$ and D$_2$ under pressure are investigated at 80 K up to 160 GPa by first-principles path-integral molecular dynamics calculations. Molecular orientations are well-defined in phase II of D$_2$, while solid H$_2$ exhibits large and very asymmetric angular quantum fluctuations in this phase, with possible rotation in the (bc) plane, making it difficult to associate a well-identified single classical structure. The mechanism for the transition to phase III is also described. Existing structural data support this microscopic interpretation.
Being the simplest element with just one electron and proton the electronic structure of the Hydrogen atom is known exactly. However, this does not hold for the complex interplay between them in a solid and in particular not at high pressure that is known to alter the crystal as well as the electronic structure. Back in 1935 Wigner and Huntington predicted that at very high pressure solid molecular hydrogen would dissociate and form an atomic solid that is metallic. In spite of intense research efforts the experimental realization, as well as the theoretical determination of the crystal structure has remained elusive. Here we present a computational study showing that the distorted hexagonal P6$_3$/m structure is the most likely candidate for Phase III of solid hydrogen. We find that the pairing structure is very persistent and insulating over the whole pressure range, which suggests that metallization due to dissociation may precede eventual bandgap closure. Due to the fact that this not only resolve one of major disagreement between theory and experiment, but also excludes the conjectured existence of phonon-driven superconductivity in solid molecular hydrogen, our results involve a complete revision of the zero-temperature phase diagram of Phase III.
Optical properties of compressed fluid hydrogen in the region where dissociation and metallization is observed are computed by ab-initio methods and compared to recent experimental results. We confirm that above 3000 K both processes are continuous while below 1500K the first order phase transition is accompanied by a discontinuity of the DC conductivity and the thermal conductivity, while both the reflectivity and absorption coefficient vary rapidly but continuously. Our results support the recent analysis of NIF experiments (P. Celliers et al, Science 361, 677-682 (2018)) which assigned the inception of metallization to pressures where the reflectivity is about 0.3. Our results also support the conclusion that the temperature plateau seen in laser-heated DAC experiments at temperatures higher than 1500 K corresponds to the onset of of optical absorption, not to the phase transition.
We present an accurate study of the static-nucleus electronic energy band gap of solid molecular hydrogen at high pressure. The excitonic and quasiparticle gaps of the $C2/c$, $Pc$, $Pbcn$, and $P6_3/m$ structures at pressures of 250, 300, and 350~GPa are calculated using the fixed-node diffusion quantum Monte Carlo (DMC) method. The difference between the mean-field and many-body band gaps at the same density is found to be almost independent of system size and can therefore be applied as a scissor correction to the mean-field gap of an infinite system to obtain an estimate of the many-body gap in the thermodynamic limit. By comparing our static-nucleus DMC energy gaps with available experimental results, we demonstrate the important role played by nuclear quantum effects in the electronic structure of solid hydrogen. Our DMC results suggest that the metallization of high-pressure solid hydrogen occurs via a structural phase transition rather than band gap closure.