At ambient pressure, sodium, chlorine, and their only known compound NaCl, have well-understood crystal structures and chemical bonding. Sodium is a nearly-free-electron metal with the bcc structure. Chlorine is a molecular crystal, consisting of Cl2
molecules. Sodium chloride, due to the large electronegativity difference between Na and Cl atoms, has highly ionic chemical bonding, with stoichiometry 1:1 dictated by charge balance, and rocksalt (B1-type) crystal structure in accordance with Paulings rules. Up to now, Na-Cl was thought to be an ultimately simple textbook system. Here, we show that under pressure the stability of compounds in the Na-Cl system changes and new materials with different stoichiometries emerge at pressure as low as 25 GPa. In addition to NaCl, our theoretical calculations predict the stability of Na3Cl, Na2Cl, Na3Cl2, NaCl3 and NaCl7 compounds with unusual bonding and electronic properties. The bandgap is closed for the majority of these materials. Guided by these predictions, we have synthesized cubic NaCl3 at 55-60 GPa in the laser-heated diamond anvil cell at temperatures above 2000 K.
The solid inner core of the Earth is predominantly composed of iron alloyed with several percent Ni and some lighter elements, Si, S, O, H, and C being the prime candidates. There have been a growing number of papers investigating C and H as possible
light elements in the core, but the results are contradictory. Here, using ab initio simulations, we study the Fe-C and Fe-H systems at inner core pressures (330-364 GPa). Using the evolutionary structure prediction algorithm USPEX, we have determined the lowest-enthalpy structures of possible carbides (FeC, Fe2C, Fe3C, Fe4C, FeC2, FeC3, FeC4 and Fe7C3) and hydrides (Fe4H, Fe3H, Fe2H, FeH, FeH2, FeH3, FeH4) and have found that Fe2C (Pnma) is the most stable iron carbide at pressures of the inner core, while FeH, FeH3 and FeH4 are stable iron hydrides at these conditions. For Fe3C, the cementite structure (Pnma) and the Cmcm structure recently found by random sampling are less stable than the I-4 and C2/m structures found here. We found that FeH3 and FeH4 adopt chemically interesting thermodynamically stable structures, in both compounds containing trivalent iron. The density of the inner core can be matched with a reasonable concentration of carbon, 11-15 mol.percent (2.6-3.7 wt.percent) at relevant pressures and temperatures. This concentration matches that in CI carbonaceous chondrites and corresponds to the average atomic mass in the range 49.3-51.0, in close agreement with inferences from the Birchs law for the inner core. Similarly made estimates for the maximum hydrogen content are unrealistically high, 17-22 mol.percent (0.4-0.5 wt.percent), which corresponds to the average atomic mass in the range 43.8-46.5. We conclude that carbon is a better candidate light alloying element than hydrogen.
Recent studies showed that hardness, a complex property, can be calculated using very simple approaches or even analytical formulae. These form the basis for evaluating controversial experimental results (as we illustrate for TiO2-cotunnite) and enab
le a systematic search for novel hard materials, for instance, using global optimization algorithms (as we show on the example of SiO2 polymorphs).
Experimental studies established that calcium undergoes several counterintuitive transitions under pressure: fcc rightarrow bcc rightarrow simple cubic rightarrow Ca-IV rightarrow 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 (I41/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 ~20 K at 120 GPa, in good agreement with experiment.
Prediction of stable crystal structures at given pressure-temperature conditions, based only on the knowledge of the chemical composition, is a central problem of condensed matter physics. This extremely challenging problem is often termed crystal st
ructure prediction problem, and recently developed evolutionary algorithm USPEX (Universal Structure Predictor: Evolutionary Xtallography) made an important progress in solving it, enabling efficient and reliable prediction of structures with up to ~40 atoms in the unit cell using ab initio methods. Here we review this methodology, as well as recent progress in analyzing energy landscape of solids (which also helps to analyze results of USPEX runs). We show several recent applications - (1) prediction of new high-pressure phases of CaCO3, (2) search for the structure of the polymeric phase of CO2 (phase V), (3) high-pressure phases of oxygen, (4) exploration of possible stable compounds in the Xe-C system at high pressures, (5) exotic high-pressure phases of elements boron and sodium.
The temperature anomalies in the Earths mantle associated with thermal convection1 can be inferred from seismic tomography, provided that the elastic properties of mantle minerals are known as a function of temperature at mantle pressures. At present
, however, such information is difficult to obtain directly through laboratory experiments. We have therefore taken advantage of recent advances in computer technology, and have performed finite-temperature ab initio molecular dynamics simulations of the elastic properties of MgSiO3 perovskite, the major mineral of the lower mantle, at relevant thermodynamic conditions. When combined with the results from tomographic images of the mantle, our results indicate that the lower mantle is either significantly anelastic or compositionally heterogeneous on large scales. We found the temperature contrast between the coldest and hottest regions of the mantle, at a given depth, to be about 800K at 1000 km, 1500K at 2000 km, and possibly over 2000K at the core-mantle boundary.
The Earths lower mantle is believed to be composed mainly of (Mg,Fe)SiO3 perovskite, with lesser amounts of (Mg,Fe)O and CaSiO3). But it has not been possible to explain many unusual properties of the lowermost 150 km of the mantle (the D layer) with
this mineralogy. Here, using ab initio simulations and high-pressure experiments, we show that at pressures and temperatures of the D layer, MgSiO3 transforms from perovskite into a layered CaIrO3-type post-perovskite phase. The elastic properties of the post-perovskite phase and its stability field explain several observed puzzling properties of the D layer: its seismic anisotropy, the strongly undulating shear-wave discontinuity at its top and possibly the anticorrelation between shear and bulk sound velocities.
The post-perovskite phase of (Mg,Fe)SiO3 is believed to be the main mineral phase of the Earths lowermost mantle (the D layer). Its properties explain numerous geophysical observations associated with this layer - for example, the D discontinuity, it
s topography and seismic anisotropy within the layer. Here we use a novel simulation technique, first-principles metadynamics, to identify a family of low-energy polytypic stacking-fault structures intermediate between the perovskite and post-perovskite phases. Metadynamics trajectories identify plane sliding involving the formation of stacking faults as the most favourable pathway for the phase transition, and as a likely mechanism for plastic deformation of perovskite and postperovskite. In particular, the predicted slip planes are (010) for perovskite (consistent with experiment) and (110) for postperovskite (in contrast to the previously expected (010) slip planes). Dominant slip planes define the lattice preferred orientation and elastic anisotropy of the texture. The (110) slip planes in post-perovskite require a much smaller degree of lattice preferred orientation to explain geophysical observations of shear-wave anisotropy in the D layer.
We have developed an efficient and reliable methodology for crystal structure prediction, merging ab initio total-energy calculations and a specifically devised evolutionary algorithm. This method allows one to predict the most stable crystal structu
re and a number of low-energy metastable structures for a given compound at any P-T conditions without requiring any experimental input. Extremely high success rate has been observed in a few tens of tests done so far, including ionic, covalent, metallic, and molecular structures with up to 40 atoms in the unit cell. We have been able to resolve some important problems in high-pressure crystallography and report a number of new high-pressure crystal structures. Physical reasons for the success of this methodology are discussed.
We explore whether the topology of energy landscapes in chemical systems obeys any rules and what these rules are. To answer this and related questions we use several tools: (i)Reduced energy surface and its density of states, (ii) descriptor of stru
cture called fingerprint function, which can be represented as a one-dimensional function or a vector in abstract multidimensional space, (iii) definition of a distance between two structures enabling quantification of energy landscapes, (iv) definition of a degree of order of a structure, and (v) definitions of the quasi-entropy quantifying structural diversity. Our approach can be used for rationalizing large databases of crystal structures and for tuning computational algorithms for structure prediction. It enables quantitative and intuitive representations of energy landscapes and reappraisal of some of the traditional chemical notions and rules. Our analysis confirms the expectations that low-energy minima are clustered in compact regions of configuration space (funnels) and that chemical systems tend to have very few funnels, sometimes only one. This analysis can be applied to the physical properties of solids, opening new ways of discovering structure-property relations. We quantitatively demonstrate that crystals tend to adopt one of the few simplest structures consistent with their chemistry, providing a thermodynamic justification of Paulings fifth rule.
