Novel hydrogen hydrate structures under pressure.

Gas hydrates are systems of prime importance. In particular, hydrogen hydrates are potential materials of icy satellites and comets, and may be used for hydrogen storage. We explore the H2O-H2 system at pressures in the range 0-100 GPa with ab initio variable-composition evolutionary simulations. According to our calculation and previous experiments, the H2O-H2 system undergoes a series of transformations with pressure, and adopts the known open-network clathrate structures (sII, C0), dense "filled ice" structures (C1, C2) and two novel hydrate phases. One of these is based on the hexagonal ice framework and has the same H2O:H2 ratio (2:1) as the C0 phase at low pressures and similar enthalpy (we name this phase Ih-C0). The other newly predicted hydrate phase has a 1:2 H2O:H2 ratio and structure based on cubic ice. This phase (which we name C3) is predicted to be thermodynamically stable above 38 GPa when including van der Waals interactions and zero-point vibrational energy, and explains previously mysterious experimental X-ray diffraction and Raman measurements. This is the hydrogen-richest hydrate and this phase has a remarkable gravimetric density (18 wt.%) of easily extractable hydrogen.

Evolutionary search for new high-k dielectric materials: methodology and applications to hafnia-based oxides.

High-k dielectric materials are important as gate oxides in microelectronics and as potential dielectrics for capacitors. In order to enable computational discovery of novel high-k dielectric materials, we propose a fitness model (energy storage density) that includes the dielectric constant, bandgap, and intrinsic breakdown field. This model, used as a fitness function in conjunction with first-principles calculations and the global optimization evolutionary algorithm USPEX, efficiently leads to practically important results. We found a number of high-fitness structures of SiO2 and HfO2, some of which correspond to known phases and some of which are new. The results allow us to propose characteristics (genes) common to high-fitness structures--these are the coordination polyhedra and their degree of distortion. Our variable-composition searches in the HfO2-SiO2 system uncovered several high-fitness states. This hybrid algorithm opens up a new avenue for discovering novel high-k dielectrics with both fixed and variable compositions, and will speed up the process of materials discovery.

Unexpected stable stoichiometries of sodium chlorides.

Sodium chloride (NaCl), or rocksalt, is well characterized at ambient pressure. As a result of the large electronegativity difference between Na and Cl atoms, it has highly ionic chemical bonding (with 1:1 stoichiometry dictated by charge balance) and B1-type crystal structure. By combining theoretical predictions and diamond anvil cell experiments, we found that new materials with different stoichiometries emerge at high pressures. Compounds such as Na3Cl, Na2Cl, Na3Cl2, NaCl3, and NaCl7 are theoretically stable and have unusual bonding and electronic properties. To test this prediction, we synthesized cubic and orthorhombic NaCl3 and two-dimensional metallic tetragonal Na3Cl. These experiments establish that compounds violating chemical intuition can be thermodynamically stable even in simple systems at nonambient conditions.

Pressure-induced stabilization and insulator-superconductor transition of BH.

Diborane (B(2)H(6)), a high energy density material, was believed to be stable in a wide P, T interval. A systematic investigation of the B-H system using the ab initio variable-composition evolutionary simulations shows that boron monohydride (BH) is thermodynamically stable and can coexist with solid B, H(2), and B(2)H(6) in a wide pressure range above 50 GPa. B(2)H(6) becomes unstable and decomposes into the Ibam phase of BH and H(2) (C2/c) at 153 GPa. The semiconducting layered Ibam structure of BH at 168 GPa transforms into a metallic phase with space group P6/mmm and a 3D topology with strong B-B and B-H covalent bonds. The Ibam-P6/mmm transformation pathway suggests the possibility of obtaining the metastable Pbcm phase on cold decompression of the P6/mmm phase. The electron-phonon coupling calculations indicate that P6/mmm-BH is a phonon-mediated superconductor with a critical temperature of superconductivity (T(c)) of 14.1-21.4 K at 175 GPa.

Novel stable compounds in the Mg-O system under high pressure.

Using ab initio evolutionary simulations, we explore the entire range of possible stoichiometries for the Mg-O system at pressures of up to 850 GPa. In addition to MgO, our calculations find that two extraordinary compounds MgO2 and Mg3O2 become thermodynamically stable at 116 GPa and 500 GPa, respectively. Detailed chemical bonding analysis shows large charge transfer in all magnesium oxides. MgO2 contains peroxide ions [O-O](2-), while non-nuclear electron density maxima play the role of anions in the electride compound Mg3O2. The latter compound is calculated to have a much narrower band gap compared to MgO and MgO2. We discuss the conditions under which MgO2 and Mg3O2 could exist in planetary interiors.

Stability of xenon oxides at high pressures.

Xenon, which is quite inert under ambient conditions, may become reactive under pressure. The possibility of the formation of stable xenon oxides and silicates in the interior of the Earth could explain the atmospheric missing xenon paradox. Using an ab initio evolutionary algorithm, we predict the existence of thermodynamically stable Xe-O compounds at high pressures (XeO, XeO(2) and XeO(3) become stable at pressures above 83, 102 and 114 GPa, respectively). Our calculations indicate large charge transfer in these oxides, suggesting that large electronegativity difference and high pressure are the key factors favouring the formation of xenon compounds. However, xenon compounds cannot exist in the Earth's mantle: xenon oxides are unstable in equilibrium with the metallic iron occurring in the lower mantle, and xenon silicates are predicted to decompose spontaneously at all mantle pressures (<136 GPa). However, it is possible that xenon atoms may be retained at defects in mantle silicates and oxides.

How evolutionary crystal structure prediction works--and why.

Once the crystal structure of a chemical substance is known, many properties can be predicted reliably and routinely. Therefore if researchers could predict the crystal structure of a material before it is synthesized, they could significantly accelerate the discovery of new materials. In addition, the ability to predict crystal structures at arbitrary conditions of pressure and temperature is invaluable for the study of matter at extreme conditions, where experiments are difficult. Crystal structure prediction (CSP), the problem of finding the most stable arrangement of atoms given only the chemical composition, has long remained a major unsolved scientific problem. Two problems are entangled here: search, the efficient exploration of the multidimensional energy landscape, and ranking, the correct calculation of relative energies. For organic crystals, which contain a few molecules in the unit cell, search can be quite simple as long as a researcher does not need to include many possible isomers or conformations of the molecules; therefore ranking becomes the main challenge. For inorganic crystals, quantum mechanical methods often provide correct relative energies, making search the most critical problem. Recent developments provide useful practical methods for solving the search problem to a considerable extent. One can use simulated annealing, metadynamics, random sampling, basin hopping, minima hopping, and data mining. Genetic algorithms have been applied to crystals since 1995, but with limited success, which necessitated the development of a very different evolutionary algorithm. This Account reviews CSP using one of the major techniques, the hybrid evolutionary algorithm USPEX (Universal Structure Predictor: Evolutionary Xtallography). Using recent developments in the theory of energy landscapes, we unravel the reasons evolutionary techniques work for CSP and point out their limitations. We demonstrate that the energy landscapes of chemical systems have an overall shape and explore their intrinsic dimensionalities. Because of the inverse relationships between order and energy and between the dimensionality and diversity of an ensemble of crystal structures, the chances that a random search will find the ground state decrease exponentially with increasing system size. A well-designed evolutionary algorithm allows for much greater computational efficiency. We illustrate the power of evolutionary CSP through applications that examine matter at high pressure, where new, unexpected phenomena take place. Evolutionary CSP has allowed researchers to make unexpected discoveries such as a transparent phase of sodium, a partially ionic form of boron, complex superconducting forms of calcium, a novel superhard allotrope of carbon, polymeric modifications of nitrogen, and a new class of compounds, perhydrides. These methods have also led to the discovery of novel hydride superconductors including the "impossible" LiH(n) (n=2, 6, 8) compounds, and CaLi(2). We discuss extensions of the method to molecular crystals, systems of variable composition, and the targeted optimization of specific physical properties.

Exotic behavior and crystal structures of calcium under pressure.

Experimental studies established that calcium undergoes several counterintuitive transitions under pressure: fcc --> bcc --> simple cubic --> Ca-IV --> Ca-V, and becomes a good superconductor in the simple cubic and higher-pressure phases. Here, using ab initio evolutionary simulations, we explore the behavior of Ca under pressure and find a number of new phases. Our structural sequence differs from the traditional picture for Ca, but is similar to that for Sr. The beta-tin (I4(1)/amd) structure, rather than simple cubic, is predicted to be the theoretical ground state at 0 K and 33-71 GPa. This structure can be represented as a large distortion of the simple cubic structure, just as the higher-pressure phases stable between 71 and 134 GPa. The structure of Ca-V, stable above 134 GPa, is a complex host-guest structure. According to our calculations, the predicted phases are superconductors with Tc increasing under pressure and reaching approximately 20 K at 120 GPa, in good agreement with experiment.

High-pressure crystal structures and superconductivity of Stannane (SnH4).

There is great interest in the exploration of hydrogen-rich compounds upon strong compression where they can become superconductors. Stannane (SnH(4)) has been proposed to be a potential high-temperature superconductor under pressure, but its high-pressure crystal structures, fundamental for the understanding of superconductivity, remain unsolved. Using an ab initio evolutionary algorithm for crystal structure prediction, we propose the existence of two unique high-pressure metallic phases having space groups Ama2 and P6(3)/mmc, which both contain hexagonal layers of Sn atoms and semimolecular (perhydride) H(2) units. Enthalpy calculations reveal that the Ama2 and P6(3)/mmc structures are stable at 96-180 GPa and above 180 GPa, respectively, while below 96 GPa SnH(4) is unstable with respect to elemental decomposition. The application of the Allen-Dynes modified McMillan equation reveals high superconducting temperatures of 15-22 K for the Ama2 phase at 120 GPa and 52-62 K for the P6(3)/mmc phase at 200 GPa.

A little bit of lithium does a lot for hydrogen.

From detailed assessments of electronic structure, we find that a combination of significantly quantal elements, six of seven atoms being hydrogen, becomes a stable metal at a pressure approximately 1/4 of that required to metalize pure hydrogen itself. The system, LiH(6) (and other LiH(n)), may well have extensions beyond the constituent lithium. These hypothetical materials demonstrate that nontraditional stoichiometries can considerably expand the view of chemical combination under moderate pressure.

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.

Transparent dense sodium.

Under pressure, metals exhibit increasingly shorter interatomic distances. Intuitively, this response is expected to be accompanied by an increase in the widths of the valence and conduction bands and hence a more pronounced free-electron-like behaviour. But at the densities that can now be achieved experimentally, compression can be so substantial that core electrons overlap. This effect dramatically alters electronic properties from those typically associated with simple free-electron metals such as lithium (Li; refs 1-3) and sodium (Na; refs 4, 5), leading in turn to structurally complex phases and superconductivity with a high critical temperature. But the most intriguing prediction-that the seemingly simple metals Li (ref. 1) and Na (ref. 4) will transform under pressure into insulating states, owing to pairing of alkali atoms-has yet to be experimentally confirmed. Here we report experimental observations of a pressure-induced transformation of Na into an optically transparent phase at approximately 200 GPa (corresponding to approximately 5.0-fold compression). Experimental and computational data identify the new phase as a wide bandgap dielectric with a six-coordinated, highly distorted double-hexagonal close-packed structure. We attribute the emergence of this dense insulating state not to atom pairing, but to p-d hybridizations of valence electrons and their repulsion by core electrons into the lattice interstices. We expect that such insulating states may also form in other elements and compounds when compression is sufficiently strong that atomic cores start to overlap strongly.