H2 Chemical Bond in a High-Pressure Crystalline Environment.

We show that the hydrogen in metal superhydride compounds can adopt two distinct states-atomic and molecular. At low pressures, the maximum number of atomic hydrogens is typically equal to the valency of the cation; additional hydrogens pair to form molecules with electronic states far below the Fermi energy causing low-symmetry structures with large unit cells. At high pressures, molecules become unstable, and all hydrogens become atomic. This study uses density functional theory, adopting BaH4 as a reference compound, which is compared with other stoichiometries and other cations. Increased temperature and zero-point motion also favor high-symmetry atomic states, and picosecond-timescale breaking and remaking of the bond permutations via intermediate H3- units.

High-temperature phase transitions in dense germanium.

Through a series of high-pressure x-ray diffraction experiments combined with in situ laser heating, we explore the pressure-temperature phase diagram of germanium (Ge) at pressures up to 110 GPa and temperatures exceeding 3000 K. In the pressure range of 64-90 GPa, we observe orthorhombic Ge-IV transforming above 1500 K to a previously unobserved high-temperature phase, which we denote as Ge-VIII. This high-temperature phase is characterized by a tetragonal crystal structure, space group I4/mmm. Density functional theory simulations confirm that Ge-IV becomes unstable at high temperatures and that Ge-VIII is highly competitive and dynamically stable at these conditions. The existence of Ge-VIII has profound implications for the pressure-temperature phase diagram, with melting conditions increasing to much higher temperatures than previous extrapolations would imply.

Synthesis of Weaire-Phelan Barium Polyhydride.

By combining pressures up to 50 GPa and temperatures of 1200 K, we synthesize the novel barium hydride, Ba8H46, stable down to 27 GPa. We use Raman spectroscopy, X-ray diffraction, and first-principles calculations to determine that this compound adopts a highly symmetric Pm3¯n structure with an unusual 534:1 hydrogen-to-barium ratio. This singular stoichiometry corresponds to the well-defined type-I clathrate geometry. This clathrate consists of a Weaire-Phelan hydrogen structure with the barium atoms forming a topologically close-packed phase. In particular, the structure is formed by H20 and H24 clathrate cages showing substantially weakened H-H interactions. Density functional theory (DFT) demonstrates that cubic Pm3¯n Ba8H46 requires dynamical effects to stabilize the H20 and H24 clathrate cages.

Unusually complex phase of dense nitrogen at extreme conditions.

Nitrogen exhibits an exceptional polymorphism under extreme conditions, making it unique amongst the elemental diatomics and a valuable testing system for experiment-theory comparison. Despite attracting considerable attention, the structures of many high-pressure nitrogen phases still require unambiguous determination. Here, we report the structure of the elusive high-pressure high-temperature polymorph ι-N2 at 56 GPa and ambient temperature, determined by single crystal X-ray diffraction, and investigate its properties using ab initio simulations. We find that ι-N2 is characterised by an extraordinarily large unit cell containing 48 N2 molecules. Geometry optimisation favours the experimentally determined structure and density functional theory calculations find ι-N2 to have the lowest enthalpy of the molecular nitrogen polymorphs that exist between 30 and 60 GPa. The results demonstrate that very complex structures, similar to those previously only observed in metallic elements, can become energetically favourable in molecular systems at extreme pressures and temperatures.

Helium-Iron Compounds at Terapascal Pressures.

We investigate the binary phase diagram of helium and iron using first-principles calculations. We find that helium, which is a noble gas and inert at ambient conditions, forms stable crystalline compounds with iron at terapascal pressures. A FeHe compound becomes stable above 4 TPa, and a FeHe_{2} compound above 12 TPa. Melting is investigated using molecular dynamics simulations, and a superionic phase with sublattice melting of the helium atoms is predicted. We discuss the implications of our predicted helium-iron phase diagram for interiors of giant (exo)planets and white dwarf stars.

Quantum and isotope effects in lithium metal.

The crystal structure of elements at zero pressure and temperature is the most fundamental information in condensed matter physics. For decades it has been believed that lithium, the simplest metallic element, has a complicated ground-state crystal structure. Using synchrotron x-ray diffraction in diamond anvil cells and multiscale simulations with density functional theory and molecular dynamics, we show that the previously accepted martensitic ground state is metastable. The actual ground state is face-centered cubic (fcc). We find that isotopes of lithium, under similar thermal paths, exhibit a considerable difference in martensitic transition temperature. Lithium exhibits nuclear quantum mechanical effects, serving as a metallic intermediate between helium, with its quantum effect-dominated structures, and the higher-mass elements. By disentangling the quantum kinetic complexities, we prove that fcc lithium is the ground state, and we synthesize it by decompression.

Stable all-nitrogen metallic salt at terapascal pressures.

The phase diagram and equation of state of dense nitrogen are of interest in understanding the fundamental physics and chemistry under extreme conditions, including planetary processes, and in discovering new materials. We predict several stable phases of nitrogen at multi-TPa pressures, including a P4/nbm structure consisting of partially charged N(2)(δ+) pairs and N(5)(δ-) tetrahedra, which is stable in the range 2.5-6.8 TPa. This is followed by a modulated layered structure between 6.8 and 12.6 TPa, which also exhibits a significant charge transfer. The P4/nbm metallic nitrogen salt and the modulated structure are stable at high pressures and temperatures, and they exhibit strongly ionic features and charge density distortions, which is unexpected in an element under such extreme conditions and could represent a new class of nitrogen materials. The P-T phase diagram of nitrogen at TPa pressures is investigated using quasiharmonic phonon calculations and ab initio molecular dynamics simulations.

Decomposition and Terapascal phases of water ice.

Computational searches for stable and metastable structures of water ice and other H:O compositions at TPa pressures have led us to predict that H(2)O decomposes into H(2)O(2) and a hydrogen-rich phase at pressures of a little over 5 TPa. The hydrogen-rich phase is stable over a wide range of hydrogen contents, and it might play a role in the erosion of the icy component of the cores of gas giants as H(2)O comes into contact with hydrogen. Metallization of H(2)O is predicted at a higher pressure of just over 6 TPa, and therefore H(2)O does not have a thermodynamically stable low-temperature metallic form. We have also found a new and rich mineralogy of complicated water ice phases that are more stable in the pressure range 0.8-2 TPa than any predicted previously.

Persistence and eventual demise of oxygen molecules at terapascal pressures.

Computational searches for structures of solid oxygen under high pressures in the multi-TPa range are carried out using density-functional-theory methods. We find that molecular oxygen persists to about 1.9 TPa at which it transforms into a semiconducting square-spiral-like polymeric structure (I4(1)/acd) with a band gap of ~3.0 eV. Solid oxygen forms a metallic zigzag chainlike structure (Cmcm) at about 3.0 TPa, but the chains in each layer gradually merge as the pressure is increased and a structure of Fmmm symmetry forms at about 9.3 TPa in which each atom has four nearest neighbors. The superconducting properties of molecular oxygen do not vary much with compression, although the structure becomes more symmetric. The electronic properties of oxygen have a complex evolution with pressure, swapping between insulating, semiconducting, and metallic.

Thermodynamically stable phases of carbon at multiterapascal pressures.

Phases of carbon are studied up to pressures of 1 petapascal (PPa) using first-principles density-functional-theory methods and a structure searching algorithm. Our extensive search over the potential energy surface supports the sequence of transitions diamond → BC8 → simple cubic under increasing pressure found in previous theoretical studies. At higher pressures we predict a soft-phonon driven transition to a simple hexagonal structure at 6.4 terapascals (TPa), and further transitions to the face centered cubic electride structure at 21 TPa, a double hexagonal close packed structure at 270 TPa, and the body centered cubic structure at 650 TPa.

Novel structures and superconductivity of silane under pressure.

Following the suggestion that hydrogen-rich compounds, and, in particular, silane (SiH4), might be high-T_{c} superconductors at moderate pressures, very recent experiments have confirmed that silane metallises and even becomes superconducting at high pressure. In this article, we present a structural characterization of compressed silane obtained with an ab initio evolutionary algorithm for crystal structure prediction. Besides the earlier molecular and chainlike structures of P2_{1}/c and I4_{1}/a symmetries, respectively, we propose two novel structures with space groups Fdd2 and Pbcn, to be stable at 25-55 and 220-250 GPa, respectively. According to our calculations, silane becomes metallic and superconducting at 220 GPa in the layered Pbcn structure, with a theoretical T_{c} of 16 K. Our calculations also show that the imaginary phonons of the recently proposed P6_{3} generate the Pbcn structure.

Superconducting high pressure phase of germane.

High-pressure structures of germane (GeH4) are explored through ab initio evolutionary methodology to reveal a metallic monoclinic structure of C2/c (4 molecules/cell). The C2/c structure consists of layerlike motifs containing novel "H2" units. Enthalpy calculations suggest a remarkably wide decomposition (Ge+H2) pressure range of 0-196 GPa, above which C2/c structure is stable. Perturbative linear-response calculations for C2/c GeH4 at 220 GPa predict a large electron-phonon coupling parameter lambda of 1.12 and the resulting superconducting critical temperature reaches 64 K.