Helium and methane are major components of giant icy planets and are abundant in the universe. However, helium is the most inert element in the periodic table and methane is one of the most hydrophobic molecules, thus whether they can react with each other is of fundamental importance. Here, our crystal structure searches and first-principles calculations predict that a HeCH compound is stable over a wide range of pressures from 55 to 155 GPa and a HeCH compound becomes stable around 105 GPa.
View Article and Find Full Text PDFWe 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.
View Article and Find Full Text PDFA 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.
View Article and Find Full Text PDFIce is one of the most extensively studied condensed matter systems. Yet, both experimentally and theoretically several new phases have been discovered over the last years. Here we report a large-scale density-functional-theory study of the configuration space of water ice.
View Article and Find Full Text PDFWe 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.
View Article and Find Full Text PDFRoom-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.
View Article and Find Full Text PDFActa Crystallogr B Struct Sci Cryst Eng Mater
August 2016
The sixth blind test of organic crystal structure prediction (CSP) methods has been held, with five target systems: a small nearly rigid molecule, a polymorphic former drug candidate, a chloride salt hydrate, a co-crystal and a bulky flexible molecule. This blind test has seen substantial growth in the number of participants, with the broad range of prediction methods giving a unique insight into the state of the art in the field. Significant progress has been seen in treating flexible molecules, usage of hierarchical approaches to ranking structures, the application of density-functional approximations, and the establishment of new workflows and `best practices' for performing CSP calculations.
View Article and Find Full Text PDFSurface 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.
View Article and Find Full Text PDFThe only known compound of sodium and hydrogen is archetypal ionic NaH. Application of high pressure is known to promote states with higher atomic coordination, but extensive searches for polyhydrides with unusual stoichiometry have had only limited success in spite of several theoretical predictions. Here we report the first observation of the formation of polyhydrides of Na (NaH3 and NaH7) above 40 GPa and 2,000 K.
View Article and Find Full Text PDFThe noble gases are the most inert group of the periodic table, but their reactivity increases with pressure. Diamond-anvil-cell experiments and ab initio modelling have been used to investigate a possible direct reaction between xenon and oxygen at high pressures. We have now synthesized two oxides below 100 GPa (Xe2O5 under oxygen-rich conditions, and Xe3O2 under oxygen-poor conditions), which shows that xenon is more reactive under pressure than predicted previously.
View Article and Find Full Text PDFThe 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).
View Article and Find Full Text PDFElectron-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.
View Article and Find Full Text PDFProc Natl Acad Sci U S A
June 2015
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.
View Article and Find Full Text PDFWe 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.
View Article and Find Full Text PDFAlkali 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.
View Article and Find Full Text PDFStructures 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.
View Article and Find Full Text PDFPhys Rev E Stat Nonlin Soft Matter Phys
February 2014
A compression algorithm is introduced for multideterminant wave functions which can greatly reduce the number of determinants that need to be evaluated in quantum Monte Carlo calculations. We have devised an algorithm with three levels of compression, the least costly of which yields excellent results in polynomial time. We demonstrate the usefulness of the compression algorithm for evaluating multideterminant wave functions in quantum Monte Carlo calculations, whose computational cost is reduced by factors of between about 2 and over 25 for the examples studied.
View Article and Find Full Text PDFWe 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.
View Article and Find Full Text PDFA 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.
View Article and Find Full Text PDFThe 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.
View Article and Find Full Text PDFThe melting temperature of solid hydrogen drops with pressure above ~65 GPa, suggesting that a liquid state might exist at low temperatures. It has also been suggested that this low-temperature liquid state might be non-molecular and metallic, although evidence for such behaviour is lacking. Here we report results for hydrogen at high pressures using ab initio methods, which include a description of the quantum motion of the protons.
View Article and Find Full Text PDFComputational 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.
View Article and Find Full Text PDFSymmetric 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.
View Article and Find Full Text PDF