Quantum hydrogen-bond symmetrization in the superconducting hydrogen sulfide system.

The quantum nature of the proton can crucially affect the structural and physical properties of hydrogen compounds. For example, in the high-pressure phases of H2O, quantum proton fluctuations lead to symmetrization of the hydrogen bond and reduce the boundary between asymmetric and symmetric structures in the phase diagram by 30 gigapascals (ref. 3). Here we show that an analogous quantum symmetrization occurs in the recently discovered sulfur hydride superconductor with a superconducting transition temperature Tc of 203 kelvin at 155 gigapascals--the highest Tc reported for any superconductor so far. Superconductivity occurs via the formation of a compound with chemical formula H3S (sulfur trihydride) with sulfur atoms arranged on a body-centred cubic lattice. If the hydrogen atoms are treated as classical particles, then for pressures greater than about 175 gigapascals they are predicted to sit exactly halfway between two sulfur atoms in a structure with Im3m symmetry. At lower pressures, the hydrogen atoms move to an off-centre position, forming a short H-S covalent bond and a longer H···S hydrogen bond in a structure with R3m symmetry. X-ray diffraction experiments confirm the H3S stoichiometry and the sulfur lattice sites, but were unable to discriminate between the two phases. Ab initio density-functional-theory calculations show that quantum nuclear motion lowers the symmetrization pressure by 72 gigapascals for H3S and by 60 gigapascals for D3S. Consequently, we predict that the Im3m phase dominates the pressure range within which the high Tc was measured. The observed pressure dependence of Tc is accurately reproduced in our calculations for the phase, but not for the R3m phase. Therefore, the quantum nature of the proton fundamentally changes the superconducting phase diagram of H3S.

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.

Evidence from Quantum Monte Carlo Simulations of Large-Gap Superfluidity and BCS-BEC Crossover in Double Electron-Hole Layers.

We report quantum Monte Carlo evidence of the existence of large gap superfluidity in electron-hole double layers over wide density ranges. The superfluid parameters evolve from normal state to BEC with decreasing density, with the BCS state restricted to a tiny range of densities due to the strong screening of Coulomb interactions, which causes the gap to rapidly become large near the onset of superfluidity. The superfluid properties exhibit similarities to ultracold fermions and iron-based superconductors, suggesting an underlying universal behavior of BCS-BEC crossovers in pairing systems.

Structure and Metallicity of Phase V of Hydrogen.

A new phase V of hydrogen was recently claimed in experiments above 325 GPa and 300 K. Because of the extremely small sample size at such record pressures the measurements were limited to Raman spectroscopy. The experimental data on increase of pressure show decreasing Raman activity and darkening of the sample, which suggests band gap closure and impending molecular dissociation, but no definite conclusions could be reached. Furthermore, the available data are insufficient to determine the structure of phase V, which remains unknown. Introducing saddle-point ab initio random structure searching, we find several new structural candidates of hydrogen which could describe the observed properties of phase V. We investigate hydrogen metallization in the proposed candidate structures, and demonstrate that smaller band gaps are associated with longer bond lengths. We conclude that phase V is a stepping stone towards metallization.

Prediction of 10-fold coordinated TiO2 and SiO2 structures at multimegabar pressures.

We predict by first-principles methods a phase transition in TiO2 at 6.5 Mbar from the Fe2P-type polymorph to a ten-coordinated structure with space group I4/mmm. This is the first report, to our knowledge, of the pressure-induced phase transition to the I4/mmm structure among all dioxide compounds. The I4/mmm structure was found to be up to 3.3% denser across all pressures investigated. Significant differences were found in the electronic properties of the two structures, and the metallization of TiO2 was calculated to occur concomitantly with the phase transition to I4/mmm. The implications of our findings were extended to SiO2, and an analogous Fe2P-type to I4/mmm transition was found to occur at 10 TPa. This is consistent with the lower-pressure phase transitions of TiO2, which are well-established models for the phase transitions in other AX2 compounds, including SiO2. As in TiO2, the transition to I4/mmm corresponds to the metallization of SiO2. This transformation is in the pressure range reached in the interiors of recently discovered extrasolar planets and calls for a reformulation of the equations of state used to model them.

Hydrogen Clathrate Structures in Rare Earth Hydrides at High Pressures: Possible Route to Room-Temperature Superconductivity.

Room-temperature superconductivity has been a long-held dream and an area of intensive research. Recent experimental findings of superconductivity at 200 K in highly compressed hydrogen (H) sulfides have demonstrated the potential for achieving room-temperature superconductivity in compressed H-rich materials. We report first-principles structure searches for stable H-rich clathrate structures in rare earth hydrides at high pressures. The peculiarity of these structures lies in the emergence of unusual H cages with stoichiometries H_{24}, H_{29}, and H_{32}, in which H atoms are weakly covalently bonded to one another, with rare earth atoms occupying the centers of the cages. We have found that high-temperature superconductivity is closely associated with H clathrate structures, with large H-derived electronic densities of states at the Fermi level and strong electron-phonon coupling related to the stretching and rocking motions of H atoms within the cages. Strikingly, a yttrium (Y) H_{32} clathrate structure of stoichiometry YH_{10} is predicted to be a potential room-temperature superconductor with an estimated T_{c} of up to 303 K at 400 GPa, as derived by direct solution of the Eliashberg equation.

Vibrational effects on surface energies and band gaps in hexagonal and cubic ice.

Surface energies of hexagonal and cubic water ice are calculated using first-principles quantum mechanical methods, including an accurate description of anharmonic nuclear vibrations. We consider two proton-orderings of the hexagonal and cubic ice basal surfaces and three proton-orderings of hexagonal ice prism surfaces, finding that vibrations reduce the surface energies by more than 10%. We compare our vibrational densities of states to recent sum frequency generation absorption measurements and identify surface proton-orderings of experimental ice samples and the origins of characteristic absorption peaks. We also calculate zero point quantum vibrational corrections to the surface electronic band gaps, which range from -1.2 eV for the cubic ice basal surface up to -1.4 eV for the hexagonal ice prism surface. The vibrational corrections to the surface band gaps are up to 12% smaller than for bulk ice.

Metallic icosahedron phase of sodium at terapascal pressures.

Alkali metals exhibit unexpected structures and electronic behavior at high pressures. Compression of metallic sodium (Na) to 200 GPa leads to the stability of a wide-band-gap insulator with the double hexagonal hP4 structure. Post-hP4 structures remain unexplored, but they are important for addressing the question of the pressure at which Na reverts to a metal. Here, we report the reentrant metallicity of Na at the very high pressure of 15.5 terapascal (TPa), predicted using first-principles structure searching simulations. Na is therefore insulating over the large pressure range of 0.2-15.5 TPa. Unusually, Na adopts an oP8 structure at pressures of 117-125 GPa and the same oP8 structure at 1.75-15.5 TPa. The metallization of Na occurs on the formation of a stable and striking body-centered cubic cI24 electride structure consisting of Na_{12} icosahedra, each housing at its center about one electron that is not associated with any Na ions.

High-pressure hydrogen sulfide from first principles: a strongly anharmonic phonon-mediated superconductor.

We use first-principles calculations to study structural, vibrational, and superconducting properties of H_{2}S at pressures P≥200 GPa. The inclusion of zero-point energy leads to two different possible dissociations of H2S, namely 3H2S→2H3S+S and 5H2S→3H3S+HS2, where both H3S and HS2 are metallic. For H3S, we perform nonperturbative calculations of anharmonic effects within the self-consistent harmonic approximation and show that the harmonic approximation strongly overestimates the electron-phonon interaction (λ≈2.64 at 200 GPa) and Tc. Anharmonicity hardens H─S bond-stretching modes and softens H─S bond-bending modes. As a result, the electron-phonon coupling is suppressed by 30% (λ≈1.84 at 200 GPa). Moreover, while at the harmonic level Tc decreases with increasing pressure, the inclusion of anharmonicity leads to a Tc that is almost independent of pressure. High-pressure hydrogen sulfide is a strongly anharmonic superconductor.

Calcium peroxide from ambient to high pressures.

Structures of calcium peroxide (CaO2) are investigated in the pressure range 0-200 GPa using the ab initio random structure searching (AIRSS) method and density functional theory (DFT) calculations. At 0 GPa, there are several CaO2 structures very close in enthalpy, with the ground-state structure dependent on the choice of exchange-correlation functional. Further stable structures for CaO2 with C2/c, I4/mcm and P21/c symmetries emerge at pressures below 40 GPa. These phases are thermodynamically stable against decomposition into CaO and O2. The stability of CaO2 with respect to decomposition increases with pressure, with peak stability occurring at the CaO B1-B2 phase transition at 65 GPa. Phonon calculations using the quasiharmonic approximation show that CaO2 is a stable oxide of calcium at mantle temperatures and pressures, highlighting a possible role for CaO2 in planetary geochemistry. We sketch the phase diagram for CaO2, and find at least five new stable phases in the pressure-temperature ranges 0 ≤ P ≤ 60 GPa, 0 ≤ T ≤ 600 K, including two new candidates for the zero-pressure ground state structure.

Vibrational renormalisation of the electronic band gap in hexagonal and cubic ice.

Electron-phonon coupling in hexagonal and cubic water ice is studied using first-principles quantum mechanical methods. We consider 29 distinct hexagonal and cubic ice proton-orderings with up to 192 molecules in the simulation cell to account for proton-disorder. We find quantum zero-point vibrational corrections to the minimum electronic band gaps ranging from -1.5 to -1.7 eV, which leads to improved agreement between calculated and experimental band gaps. Anharmonic nuclear vibrations play a negligible role in determining the gaps. Deuterated ice has a smaller band-gap correction at zero-temperature of -1.2 to -1.4 eV. Vibrations reduce the differences between the electronic band gaps of different proton-orderings from around 0.17 eV to less than 0.05 eV, so that the electronic band gaps of hexagonal and cubic ice are almost independent of the proton-ordering when quantum nuclear vibrations are taken into account. The comparatively small reduction in the band gap over the temperature range 0 - 240 K of around 0.1 eV does not depend on the proton ordering, or whether the ice is protiated or deuterated, or hexagonal, or cubic. We explain this in terms of the atomistic origin of the strong electron-phonon coupling in ice.

Temperature effects in first-principles solid state calculations of the chemical shielding tensor made simple.

We study the effects of atomic vibrations on the solid-state chemical shielding tensor using first principles density functional theory calculations. At the harmonic level, we use a Monte Carlo method and a perturbative expansion. The Monte Carlo method is accurate but computationally expensive, while the perturbative method is computationally more efficient, but approximate. We find excellent agreement between the two methods for both the isotropic shift and the shielding anisotropy. The effects of zero-point quantum mechanical nuclear motion are important up to relatively high temperatures: at 500 K they still represent about half of the overall vibrational contribution. We also investigate the effects of anharmonic vibrations, finding that their contribution to the zero-point correction to the chemical shielding tensor is small. We exemplify these ideas using magnesium oxide and the molecular crystals L-alanine and ß-aspartyl-L-alanine. We therefore propose as the method of choice to incorporate the effects of temperature in solid state chemical shielding tensor calculations using the perturbative expansion within the harmonic approximation. This approach is accurate and requires a computational effort that is about an order of magnitude smaller than that of dynamical or Monte Carlo approaches, so these effects might be routinely accounted for.

Excitations and benchmark ensemble density functional theory for two electrons.

A new method for extracting ensemble Kohn-Sham potentials from accurate excited state densities is applied to a variety of two-electron systems, exploring the behavior of exact ensemble density functional theory. The issue of separating the Hartree energy and the choice of degenerate eigenstates is explored. A new approximation, spin eigenstate Hartree-exchange, is derived. Exact conditions that are proven include the signs of the correlation energy components and the asymptotic behavior of the potential for small weights of the excited states. Many energy components are given as a function of the weights for two electrons in a one-dimensional flat box, in a box with a large barrier to create charge transfer excitations, in a three-dimensional harmonic well (Hooke's atom), and for the He atom singlet-triplet ensemble, singlet-triplet-singlet ensemble, and triplet bi-ensemble.

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.

Excitons and biexcitons in symmetric electron-hole bilayers.

Symmetric electron-hole bilayer systems have been studied at zero temperature using the diffusion quantum Monte Carlo method. A flexible trial wave function is used that can describe fluid, excitonic, and biexcitonic phases. We calculate condensate fractions and pair correlation functions for a large number of densities r(s) and layer separations d. At small d we find a one-component fluid phase, an excitonic fluid phase, and a biexcitonic fluid phase, and the transitions among them appear to be continuous. At d=0, excitons appear to survive down to about r(s)=0.5 a.u., and biexcitons form at r(s)>2.5 a.u.

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.

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.

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.

High pressure ionic and molecular crystals of ammonia monohydrate within density functional theory.

A combination of first-principles density functional theory calculations and a search over structures is used to predict the stability of a proton-transfer modification of ammonia monohydrate with space group P4∕nmm. The phase diagram is calculated with the Perdew-Burke-Ernzerhof (PBE) density functional, and the effects of a semi-empirical dispersion correction, zero point motion, and finite temperature are investigated. Comparison with MP2 and coupled cluster calculations shows that the PBE functional over-stabilizes proton transfer phases because too much electronic charge moves with the proton. This over-binding is partially corrected by using the PBE0 hybrid exchange-correlation functional, which increases the enthalpy of P4∕nmm by about 0.6 eV per formula unit relative to phase I of ammonia monohydrate and shifts the transition to the proton transfer phase from the PBE pressure of 2.8 GPa to about 10 GPa. This is consistent with experiment as proton transfer phases have not been observed at pressures up to ∼9 GPa, while higher pressures have not yet been explored experimentally.