High-pressure structures of solid hydrogen: Insights from ab initio molecular dynamics simulations

Understanding the structural behavior of solid hydrogen under high pressures is crucial for uncovering its unique properties and potential applications. In this study, starting from the phase I of solid hydrogen-free-rotator hcp structure, we conduct extensive ab initio molecular dynamics calculations to simulate the cooling, heating, and equilibrium processes within a pressure range of 80-260 GPa. Without relying on any structure previously predicted, we identify the high-pressure phase structures of solid hydrogen as P21/c for phase II, P6522 for phase III, and BG1BG2BG3 six-layer structure for phase IV, which are different from those proposed previously using the structure-search method. The reasonability of these structures are validated by Raman spectra and x-ray diffraction patterns by comparison with the experimental results. Our results actually show pronounced changes in the c/a ratio between phases I, III, and IV, which hold no brief for the experimental interpretation of an isostructural hcp transformations for phases I-III-IV.

Isotope Quantum Effects in the Metallization Transition in Liquid Hydrogen

Quantum effects in condensed matter normally only occur at low temperatures. Here we show a large quantum effect in high-pressure liquid hydrogen at thousands of Kelvins. We show that the metallization transition in hydrogen is subject to a very large isotope effect, occurring hundreds of degrees lower than the equivalent transition in deuterium. We examined this using path integral molecular dynamics simulations which identify a liquid-liquid transition involving atomization, metallization, and changes in viscosity, specific heat, and compressibility. The difference between H_{2} and D_{2} is a quantum mechanical effect that can be associated with the larger zero-point energy in H_{2} weakening the covalent bond. Our results mean that experimental results on deuterium must be corrected before they are relevant to understanding hydrogen at planetary conditions.

Understanding high pressure molecular hydrogen with a hierarchical machine-learned potential

The hydrogen phase diagram has several unusual features which are well reproduced by density functional calculations. Unfortunately, these calculations do not provide good physical insights into why those features occur. Here, we present a fast interatomic potential, which reproduces the molecular hydrogen phases: orientationally disordered Phase I; broken-symmetry Phase II and reentrant melt curve. The H₂ vibrational frequency drops at high pressure because of increased coupling between neighbouring molecules, not bond weakening. Liquid H₂ is denser than coexisting close-packed solid at high pressure because the favored molecular orientation switches from quadrupole-energy-minimizing to steric-repulsion-minimizing. The latter allows molecules to get closer together, without the atoms getting closer, but cannot be achieved within in a close-packed layer due to frustration. A similar effect causes negative thermal expansion. At high pressure, rotation is hindered in Phase I, such that it cannot be regarded as a molecular rotor phase.

New possible candidate structure for phase IV of solid hydrogen

It has been proved in experiments that there are at least five phases of solid hydrogen at high pressure, however, only the structure of phase I has been absolutely determined. We revisited the phase space of solid hydrogen in the pressure range of 200-500 GPa using the particle swarm optimization technique combined with first-principles simulations. A novel orthorhombic structure named Ama2 is proposed as a possible candidate structure for phase IV. The Ama2 structure is a 'mixed structure' with two different types of layers and is distinctly different from the previously reported Pc structure. Enthalpies and Gibbs free energies show that Ama2 and Pc are competitive in the pressure region of phase IV. Nevertheless, the Raman and infrared vibron frequencies of Ama2 calculated by using density functional perturbation theory based on first-principles lattice dynamics show a better agreement with the experimental measurements than those of the Pc structure. And the pressure dependence of these low-frequency Raman vibrons of Ama2 obtained from the first-principles molecular dynamics simulation shows a steeper slope, which resolves the long-standing issue of large discrepancies between the calculated Raman frequencies and the experimental ν ₁ [P. Loubeyre, F. Occelli and P. Dumas, Phys. Rev. B: Condens. Matter Mater. Phys., 2013, 87, 134101 and C. S. Zha, R. E. Cohen, H. K. Mao and R. J. Hemley, Proc. Natl. Acad. Sci. U.S.A., 2014, 111, 4792]. Structural and vibrational analyses show that the hydrogen molecules in the weakly bonded molecular layer of Ama2 form distorted hexagonal patterns, and their vibration can be used to explain the experimental ν ₁ vibron. It is found that the weakly bonded layer is almost the same as the layers in the C2/c structure. This confirms the experimental conclusion [P. Loubeyre, F. Occelli and P. Dumas, Phys. Rev. B: Condens. Matter Mater. Phys., 2013, 87, 134101] that the ordering of hydrogen molecules in the weakly bonded molecular layers of the 'mixed structure' for phase IV is similar to that in the layers of the C2/c structure.

Distinguishing the Structure of High-Pressure Hydrogen with Dielectric Constants

Identifying the structures of high-pressure hydrogen has been one of the central goals in high-pressure physics; however, it still presents a fundamental challenge because of the lack of an effective measure for distinguishing the structures. Herein, we address this issue by focusing on the potential candidates of phases II and III of high-pressure hydrogen. We find that the anisotropic dielectric constants of the different hydrogen solids and their responses to pressure behave differently depending on the atomic structures, corresponding to the different polarization responses of the structures to the external electric field. These findings are robust regardless of the quantum and thermal motion of hydrogen solids. Therefore, the anisotropic dielectric property can serve as a potential measure for probing the structures of high-pressure hydrogen as well as other high-pressure materials.