Hybrid Functional Study Rationalizes the Simple Cubic Phase of Calcium at High Pressures

Simple cubic (SC) phase has been long experimentally determined as the high-pressure phase III of elemental calcium (Ca) since 1984. However, recent density functional calculations within semi-local approximation showed that this SC phase is structurally unstable by exhibiting severely imaginary phonons, and is energetically unstable with respect to a theoretical body-centered tetragonal I41/amd structure over the pressure range of phase III. These calculations generated extensive debates on the validity of SC phase. Here we have re-examined the SC structure by performing more precise density functional calculations within hybrid functionals of Heyd-Scuseria-Erhzerhof (HSE) and PBE0. Our calculations were able to rationalize fundamentally the phase stability of SC structure over all other known phases by evidence of its actual energetic stability above 33 GPa and its intrinsically dynamical stability without showing any imaginary phonons in the entire pressure range studied. We further established that the long-thought theoretical I41/amd structure remains stable in a narrow pressure range before entering SC phase and is actually the structure of experimental Ca-III synthesized recently at low temperature 14 K as supported by the excellent agreement between our simulated X-ray diffraction patterns and the experimental data. Our results shed strong light on the crucial role played by the precise electron exchange energy in a proper description of the potential energy of Ca.

Room-Temperature Structures of Solid Hydrogen at High Pressures

By employing first-principles metadynamics simulations, we explore the 300 K structures of solid hydrogen over the pressure range 150-300 GPa. At 200 GPa, we find the ambient-pressure disordered hexagonal close-packed (hcp) phase transited into an insulating partially ordered hcp phase (po-hcp), a mixture of ordered graphene-like H2 layers and the other layers of weakly coupled, disordered H2 molecules. Within this phase, hydrogen remains in paired states with creation of shorter intra-molecular bonds, which are responsible for the very high experimental Raman peak above 4000 cm-1. At 275 GPa, our simulations predicted a transformation from po-hcp into the ordered molecular metallic Cmca phase (4 molecules/cell) that was previously proposed to be stable only above 400 GPa. Gibbs free energy calculations at 300 K confirmed the energetic stabilities of the po-hcp and metallic Cmca phases over all known structures at 220-242 GPa and >242 GPa, respectively. Our simulations highlighted the major role played by temperature in tuning the phase stabilities and provided theoretical support for claimed metallization of solid hydrogen below 300 GPa at 300 K.