Robust Diffusive Proton Motions in Phase IV of Solid Hydrogen

Systematic first-principles molecular dynamics (MD) simulations with long simulation times (7-13 ps) for phase IV of solid hydrogen using different supercell sizes of 96, 288, 576, and 768 atoms established that the diffusive proton motions process in the graphene-like layer is an intrinsic property and independent of the simulation cell sizes. The present study highlights an often overlook issue in first-principles calculations that long time MD is essential to achieve ergodicity, which is mandatory for a proper description of dynamics of a system. The present results contradict a recent work [Phys. Rev. B 87, 174110 (2013)] in which the analysis was relied on short time slices (1-3 ps).

Hardness of FeB4: Density functional theory investigation

A recent experimental study reported the successful synthesis of an orthorhombic FeB4 with a high hardness of 62 GPa, which has reignited extensive interests on whether transition metal borides (TRBs) compounds will become superhard materials. However, it is contradicted with some theoretical studies suggesting transition metal boron compounds are unlikely to become superhard materials. Here, we examined structural and electronic properties of FeB4 using density functional theory. The electronic calculations show the good metallicity and covalent FeB bonding. Meanwhile, we extensively investigated stress strain relations of FeB4 under various tensile and shear loading directions. The calculated weakest tensile and shear stresses are 40 GPa and 25 GPa, respectively. Further simulations (e.g. electron localized function and bond length along the weakest loading direction) on FeB4 show the weak Fe-B bonding is responsible for this low hardness. Moreover, these results are consistent with the value of Vickers hardness (11.7 to 32.3 GPa) by employing different empirical hardness models and below the superhardness threshold of 40 GPa. Our current results suggest FeB4 is a hard material and unlikely to become superhard.

Anomalous melting behavior of solid hydrogen at high pressures

Hydrogen is the most abundant element in the universe, and its properties under conditions of high temperature and pressure are crucial to understand the interior of of large gaseous planets and other astrophysical bodies. At ultra high pressures solid hydrogen has been predicted to transform into a quantum fluid, because of its high zero point motion. Here we report first principles two phase coexistence and Z method determinations of the melting line of solid hydrogen in a pressure range spanning from 30 to 600 GPa. Our results suggest that the melting line of solid hydrogen, as derived from classical molecular dynamics simulations, reaches a minimum of 367 K at about 430 GPa, at higher pressures the melting line of the atomics Cs IV phase regain a positive slope. In view of the possible importance of quantum effects in hydrogen at such low temperatures, we also determined the melting temperature of the atomic CsIV phase at pressures of 400, 500, 600 GPa, employing Feynman path integral simulations. These result in a downward shift of the classical melting line by about 100 K, and hint at a possible secondary maximum in the melting line in the region between 500 and 600 GPa, testifying to the importance of quantum effects in this system. Combined, our results imply that the stability field of the zero temperature quantum liquid phase, if it exists at all, would only occur at higher pressures than previously thought.

Proton Transfer in Phase IV of Solid Hydrogen and Deuterium

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.

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.