No Arabic abstract
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 structure 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.
Based on the unbiased structure prediction, we showed that the stable form of NiSi compound under the pressure of 100 and 200 GPa is the Pmmn-structure. Furthermore, we discovered a new stable phase - the deformed tetragonal CsCl-type structure with a = 2.174 {AA} and c = 2.69 {AA} at 400 GPa. Specifically, the sequence of high-pressure phase transitions is the following: the Pmmn-structure - below 213 GPa, the tetragonal CsCl-type - in the range 213-522 GPa, and cubic CsCl - higher than 522 GPa. As the CsCl-type structure is considered as the model structure of FeSi compound at the conditions of the Earths core, this result implies restrictions on the Fe-Ni isomorphic miscibility in FeSi.
In the present work we have proposed the method that allows one to easily estimate hardness and bulk modulus of known or hypothetical solid phases from the data on Gibbs energy of atomization of the elements and corresponding covalent radii. It has been shown that hardness and bulk moduli of compounds strongly correlate with their thermodynamic and structural properties. The proposed method may be used for a large number of compounds with various types of chemical bonding and structures; moreover, the temperature dependence of hardness may be calculated, that has been performed for diamond and cubic boron nitride. The correctness of this approach has been shown for the recently synthesized superhard diamond-like BC5. It has been predicted that the hypothetical forms of B2O3, diamond-like boron, BCx and COx, which could be synthesized at high pressures and temperatures, should have extreme hardness.
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 structure 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 discuss the efficacy of evolutionary method for the purpose of structural analysis of amorphous solids. At present ab initio molecular dynamics (MD) based melt-quench technique is used and this deterministic approach has proven to be successful to study amorphous materials. We show that a stochastic approach motivated by Darwinian evolution can also be used to simulate amorphous structures. Applying this method, in conjunction with density functional theory (DFT) based electronic, ionic and cell relaxation, we re-investigate two well known amorphous semiconductors, namely silicon and indium gallium zinc oxide (IGZO). We find that characteristic structural parameters like average bond length and bond angle are within $sim$ 2% to those reported by ab initio MD calculations and experimental studies.
The layered crystal of EuSn$_2$As$_2$ has a Bi$_2$Te$_3$-type structure in rhombohedral ($Rbar{3}m$) symmetry and has been confirmed to be an intrinsic magnetic topological insulator at ambient conditions. Combining {it ab initio} calculations and emph{in-situ} x-ray diffraction measurements, we identify a new monoclinic EuSn$_2$As$_2$ structure in $C2/m$ symmetry above $sim$14 GPa. It has a three-dimensional network made up of honeycomb-like Sn sheets and zigzag As chains, transformed from the layered EuSn$_2$As$_2$ via a two-stage reconstruction mechanism with the connecting of Sn-Sn and As-As atoms successively between the buckled SnAs layers. Its dynamic structural stability has been verified by phonon mode analysis. Electrical resistance measurements reveal an insulator-metal-superconductor transition at low temperature around 5 and 15 GPa, respectively, according to the structural conversion, and the superconductivity with a textit{T}${rm {_C}}$ value of $sim 4$ K is observed up to 30.8 GPa. These results establish a high-pressure EuSn$_2$As$_2$ phase with intriguing structural and electronic properties and expand our understandings about the layered magnetic topological insulators.