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.
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 behavi
our. 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 and sodium, 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 and Na 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 200 GPa (corresponding to 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.
Boron is a unique element, being the only element, all known polymorphs of which are superhard, and all of its crystal structures are distinct from any other element. The electron-deficient bonding in boron explains its remarkable sensitivity to even
small concentrations of impurity atoms and allows boron to form peculiar chemical compounds with very different elements. These complications made the study of boron a great challenge, creating also a unique and instructive chapter in the history of science. Strange though it may sound, the discovery of boron in 1808 was ambiguous, with pure boron polymorphs established only starting from the 1950s-1970s, and only in 2007 was the stable phase at ambient conditions determined. The history of boron research from its discovery to the latest discoveries pertaining to the phase diagram of this element, the structure and stability of beta-boron, and establishment of a new high-pressure polymorph, gamma-boron, is reviewed.
This Comment points out a number of errors in the recent paper by Zarechnaya, Dubrovinskaia, Dubrovinsky, et al. (Phys. Rev. Lett. 102, 185501 (2009)). Results and conclusions presented by Zarechnaya et al. (2009) are either incorrect or have been presented before.
