BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

Artem R. Oganov

Artem R. Oganov appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record
cond-mat.mtrl-scicond-mat.supr-conphysics.comp-phphysics.chem-phphysics.geo-phastro-ph.EPcond-mat.othercond-mat.softcond-mat.str-elphysics.atm-clusphysics.optics

Catalog footprint

What is connected

71works
11topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

71 published item(s)

preprint2022arXiv

Effect of paramagnetic impurities on superconductivity in polyhydrides: $\textit{s}$-wave order parameter in Nd-doped LaH$_{10}$

Polyhydrides are a novel class of superconducting materials with extremely high critical parameters, which is very promising for applications. On the other hand, complete experimental study of the magnetic phase diagram for the best so far known superconductor, lanthanum decahydride LaH$_{10}$, encounters a serious complication because of the large upper critical magnetic field $\textit{H}$$_{C2}$(0), exceeding 120-160 T. Partial replacement of La atoms by magnetic Nd atoms results in a decrease of the upper critical field, which makes it attainable for existing pulse magnets. We found that addition of neodymium leads to significant suppression of superconductivity in LaH$_{10}$: each atomic % of Nd causes decrease in $\textit{T}$$_{C}$ by 10-11 K. Using strong pulsed magnetic fields up to 68 T, we constructed the magnetic phase diagram of the ternary (La,Nd)H$_{10}$ superhydride, which appears to be surprisingly linear with $\textit{H}$$_{C2}$ $\propto$ |$\textit{T}$ - $\textit{T}$$_C$|. The pronounced suppression of superconductivity in LaH$_{10}$ by magnetic Nd atoms and the robustness of $\textit{T}$$_C$ with respect to nonmagnetic impurities (e.g., Y, Al, C) under Anderson's theorem indicate the isotropic ($\textit{s}$-wave) character of conventional electron-phonon pairing in the synthesized superhydrides.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Electronegativity and chemical hardness of the elements under pressure

Abundant evidence has shown the emergence of exotic chemical phenomena under pressure, including the formation of unexpected compounds and strange crystal structures. In many cases, there is no convincing explanation for these phenomena and there are virtually no chemical rules or models capable of predicting or even rationalizing these phenomena. Here we calculate, as a function of pressure, two central chemical properties of atoms, electronegativity and chemical hardness, which can be seen as the first and second-order chemical potentials. Mulliken electronegativity, which equals minus the chemical potential of the electron relative to the vacuum, is appropriately modified - instead of taking the vacuum (impossible under high pressure), we take the homogeneous electron gas as reference. We find that for most elements, chemical hardness and electronegativity decrease with pressure, consistent with pressure-induced metallization. Furthermore, we discover that pressure-induced s-d orbital transfer makes Ni, Pd and Pt "pseudo-noble-gas" atoms with a closed d-shell configuration, and the elements preceding them (Fe and especially Co, Rh, Ir) electron acceptors, while the elements right after them (Cu, Ag, Zn, Cd, for example) become highly electropositive. We show the explicative and predictive power of our electronegativity and chemical hardness scales under pressure.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Fast general two- and three-body interatomic potential

We introduce a new class of machine learning interatomic potentials - fast General Two- and Three-body Potential (GTTP), which is as fast as conventional empirical potentials and require computational time that remains constant with increasing fitting flexibility. GTTP does not contain any assumptions about the functional form of two- and three-body interactions. These interactions can be modeled arbitrarily accurately, potentially by thousands of parameters not affecting resulting computational cost. Time complexity is O(1) per every considered pair or triple of atoms. The fitting procedure is reduced to simple linear regression on ab initio calculated energies and forces and leads to effective two- and three-body potential, reproducing quantum many-body interactions as accurately as possible. Our potential can be made continuously differentiable any number of times at the expense of increased computational time. We made a number of performance tests on one-, two- and three-component systems. The flexibility of the introduced approach makes the potential transferable in terms of size and type of atomic systems. We show that trained on randomly generated structures with just 8 atoms in the unit cell, it significantly outperforms common empirical interatomic potentials in the study of large systems, such as grain boundaries in polycrystalline materials.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

PNO: A Promising Deep-UV Nonlinear Optical Material with Extremely High Second Harmonic Generation Effect

In this work, the polar tetrahedron [PN$_2$O$_2$] was revealed as a new deep-ultraviolet (deep-UV) nonlinear optically active unit. Accordingly, a thermodynamically stable compound (PNO) consisting of the polar [PN$_2$O$_2$] units was predicted and suggested as a promising candidate of deep-UV nonlinear optical (NLO) material. Compared with other deep-UV materials known to date, PNO possesses the strongest second harmonic generation (SHG) coefficient (about 6 times that of KH$_2$PO$_4$ (KDP)). Moreover, its three-dimensional connectivity endows it with good mechanical and thermal properties. Therefore, PNO should be a new option for non-$π$-conjugated deep-UV NLO materials.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Ultrahigh-Pressure Magnesium Hydrosilicates as Reservoirs of Water in Early Earth

The origin of water on the Earth is a long-standing mystery, requiring a comprehensive search for hydrous compounds, stable at conditions of the deep Earth and made of Earth-abundant elements. Previous studies usually focused on the current range of pressure-temperature conditions in the Earth's mantle and ignored a possible difference in the past, such as the stage of the core-mantle separation. Here, using ab initio evolutionary structure prediction, we find that only two magnesium hydrosilicate phases are stable at megabar pressures, $α$-Mg$_2$SiO$_5$H$_2$ and $β$-Mg$_2$SiO$_5$H$_2$, stable at 262-338 GPa and >338 GPa,respectively (all these pressures now lie within the Earth's iron core). Both are superionic conductors with quasi-one-dimensional proton diffusion at relevant conditions. In the first 30 million years of Earth's history, before the Earth's core was formed, these must have existed in the Earth, hosting much of Earth's water. As dense iron alloys segregated to form the Earth's core, Mg$_2$SiO$_5$H$_2$ phases decomposed and released water. Thus, now-extinct Mg$_2$SiO$_5$H$_2$ phases have likely contributed in a major way to the evolution of our planet.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Crystal structure prediction at finite temperatures

Crystal structure prediction is a central problem of theoretical crystallography and materials science, which until mid-2000s was considered intractable. Several methods, based on either energy landscape exploration$^{1,2}$ or, more commonly, global optimization$^{3-8}$, largely solved this problem and enabled fully non-empirical computational materials discovery$^{9,10}$. A major shortcoming is that, to avoid expensive calculations of the entropy, crystal structure prediction was done at zero Kelvin and searched for the global minimum of the enthalpy, rather than free energy. As a consequence, high-temperature phases (especially those which are not quenchable to zero temperature) could be missed. Here we develop an accurate and affordable solution, enabling crystal structure prediction at finite temperatures. Structure relaxation and fully anharmonic free energy calculations are done by molecular dynamics with a force field (which can be anything from a parametric force field for simpler cases to a trained on-the-fly machine learning interatomic potential), the errors of which are corrected using thermodynamic perturbation theory to yield accurate ab initio results. We test the accuracy of this method on metals (probing the P-T phase diagram of Al and Fe), a refractory intermetallide (WB), and a significantly ionic ceramic compound (Earth-forming silicate MgSiO3 at pressures and temperatures of the Earth's lower mantle). We find that the hcp-phase of aluminum has a wider stability field than previously thought, and the temperature-induced transition $α$-$β$ in WB occurs at 2789 K. It is also found that iron has hcp structure at conditions of the Earth's inner core, and the much debated (and important for constraining Earth's thermal structure) Clapeyron slope of the post-perovskite phase transition in MgSiO3 is 5.88 MPa/K.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Sr-Doped Molecular Hydrogen: Synthesis and Properties of SrH$_{22}$

Recently, several research groups announced reaching the point of metallization of hydrogen above 400 GPa. Following the mainstream of extensive investigations of compressed polyhydrides, in this work we demonstrate that small (4 atom %) doping of molecular hydrogen by strontium leads to a dramatic reduction in the metallization pressure to about 200 GPa. Studying the high-pressure chemistry of the Sr-H system at 56-180 GPa, we observed the formation of several previously unknown compounds: C2/m-Sr$_3$H$_{13}$, pseudocubic SrH$_6$, SrH$_9$ with cubic F-43m Sr sublattice, and pseudotetragonal P1-SrH$_{22}$, the metal hydride with the highest hydrogen content discovered so far. Unlike Ca and Y, strontium forms molecular semiconducting polyhydrides, whereas calcium and yttrium polyhydrides are high-Tc superconductors with an atomic H sublattice. The latter phase, SrH$_{22}$ or Sr$_{0.04}$H$_{0.96}$, may be considered as a convenient model of the consistent bandgap closure and metallization of hydrogen. Using the impedance measurements in diamond anvil cells at 300-440 K, we estimated the direct bandgap of the Pm-3n-like compound P1-SrH$_6$ to be 0.44-0.51 eV at 150 GPa, and its metallization pressure to be 220 GPa. Together with the machine learning interatomic potentials, the impedance spectroscopy allowed us to estimate the diffusion coefficients of hydrogen D$_H$ = 1.0-2.8 E-10 m$^2$/s in SrH$_6$ and 1.2-2.1 E-9 m$^2$/s in P1-SrH$_{22}$ at 500-600 K.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Coevolutionary search for optimal materials in the space of all possible compounds

Over the past decade, evolutionary algorithms, data mining, and other methods showed great success in solving the main problem of theoretical crystallography: finding the stable structure for a given chemical composition. Here we develop a method that addresses the central problem of computational materials science: the prediction of material(s), among all possible combinations of all elements, that possess the best combination of target properties. This nonempirical method combines our new coevolutionary approach with the carefully restructured "Mendelevian" chemical space, energy filtering, and Pareto optimization to ensure that the predicted materials have optimal properties and a high chance to be synthesizable. The first calculations, presented here, illustrate the power of this approach. In particular, we find that diamond (and its polytypes, including lonsdaleite) are the hardest possible materials and that bcc-Fe has the highest zero-temperature magnetization among all possible compounds.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Collapse modes in SC and BCC arrangements of elastic beads

Collapse modes in compressed simple cubic (SC) and body-centered cubic (BCC) periodic arrangements of elastic frictionless beads were studied numerically using the discrete element method. Under pure hydrostatic compression, the SC arrangement tends to transform into a defective hexagonal close-packed or amorphous structure. The BCC assembly exhibits several modes of collapse, one of which, identified as cI16 structure, is consistent with the behavior of BCC metals Li and Na under high pressure. The presence of a deviatoric stress leads to the transformation of the BCC structure into face-centered cubic (FCC) one via the Bain path. The observed effects provide important insights on the origins of mechanical behavior of atomic systems, while the elastic spheres model used in our work can become a useful paradigm, expanding the capabilities of a hard sphere model widely used in many branches of science.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Helium Induced Nitrogen Salt at High Pressure

The energy landscape of helium-nitrogen mixtures is explored by ab initio evolutionary searches, which predicted several stable helium-nitrogen compounds in the pressure range from 25 to 100 GPa. In particular, the monoclinic structure of HeN$_{22}$ consists of neutral He atoms, partially ionic dimers N$_{2}$$^{δ-}$, and lantern-like cages N$_{20}$$^{δ+}$. The presence of helium not only greatly enhances structural diversity of nitrogen solids, but also tremendously lowers the formation pressure of nitrogen salt. The unique nitrogen framework of (HeN$_{20}$)$^{δ+}$N$_{2}$$^{δ-}$ may be quenchable to ambient pressure even after removing helium. The estimated energy density of N$_{20}$$^{δ+}$N$_{2}$$^{δ-}$ (10.44 kJ/g) is $\sim$2.4 times larger than that of trinitrotoluene (TNT), indicating a very promising high-energy-density material.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

High-Pressure Synthesis of Barium Superhydrides: Pseudocubic BaH12

Following the discovery of high-temperature superconductivity in the La-H system, where for the recently discovered fcc-LaH10 a record critical temperature Tc = 250 K was achieved [Drozdov et al., Nature, 569, 528 (2019) and Somayazulu et al., Phys. Rev. Lett. 122, 027001 (2019)], we studied the formation of new chemical compounds in the barium-hydrogen system at pressures up to 173 GPa. Using in situ generation of hydrogen from NH3BH3, we synthesized previously unknown superhydride BaH12 with a pseudocubic (fcc) Ba sublattice, which was observed in a wide range of pressures from 75 to 173 GPa in four independent experiments. DFT calculations indicate a close agreement between the theoretical and experimental equations of state. In addition to BaH12, we identified previously known P6/mmm BaH2 and possibly BaH10 and BaH6 as impurities in the samples. Ab initio calculations show that newly discovered semimetallic BaH12 contains H2, H3 molecular units and detached H12 chains. Barium dodecahydride is a unique molecular hydride with metallic conductivity which demonstrates a superconducting transition around 20 K at 140 GPa in agreement with calculations (19-32 K). The interpretation of the multiphase XRD data was possible thanks to the development of new Python scripts for postprocessing the results of evolutionary searches. These scripts help quickly identify the theoretical structures that explain the experimental data in the best way, among thousands of candidates.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

High-Pressure Synthesis of Magnetic Neodymium Polyhydrides

The current search for room-temperature superconductivity is inspired by the unique properties of the electron-phonon interaction in metal superhydrides. Encouraged by the recently found highest-$T_C$ superconductor fcc-$LaH_{10}$, here we discover several superhydrides of another lanthanide - neodymium. We identify three novel metallic Nd-H phases at pressure range from 85 to 135 GPa: $I4/mmm$-$NdH_4$, $C2/c$-$NdH_7$, $P6_3/mmc$-$NdH_9$, synthesized by laser-heating metal samples in NH3BH3 media for in situ generation of hydrogen. A lower trihydride $Fm\bar{3}m$-$NdH_3$ is found at pressures from 2 to 52 GPa. $I4/mmm$-$NdH_4$ and $C2/c$-$NdH_7$ are stable from 135 down to 85 GPa, and $P6_3/mmc$-$NdH_9$ from 110 to 130 GPa. Measurements of the electrical resistance of NdH9 demonstrate a possible superconducting transition at ~4.5 K in $P6_3/mmc$-$NdH_9$. Our theoretical calculations predict that all the neodymium hydrides have antiferromagnetic order at pressures below 150 GPa and represent one of the first discovered examples of strongly correlated superhydrides with large exchange spin-splitting in the electron band structure (> 450 meV). The critical N$é$el temperatures for new neodymium hydrides are estimated using the mean-field approximation as about 4 K ($NdH_4$), 251 K ($NdH_7$) and 136 K ($NdH_9$).

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Hydrocarbons under pressure: phase diagrams and surprising new compounds in the C-H system

Understanding the high-pressure behavior of C-H system is of great importance due to its key role in organic, bio-, petroleum and planetary chemistry. We have performed a systematic investigation of the pressure-composition phase diagram of the C-H system at pressures up to 400 GPa using evolutionary structure prediction coupled with ab initio calculations and discovered that only saturated hydrocarbons are thermodynamically stable. Several stable methane-hydrogen co crystals are predicted: 2CH4 * H2, earlier obtained experimentally, is predicted to have I4/m space group and 2-90 GPa stability range at 0 K, and two new thermodynamically stable compounds 2CH4 * 7H2 (P-3m1 space group) and CH4 * 9H2 (Cm space group), as potential energy storage materials. P21/c phase of methane is predicted to be stable at pressures < 8 GPa; bulk graphane (CH) was shown to be thermodynamically stable at 7-18 and 18-50 GPa and 0 K in the P-3m and Cmca phases, respectively; polyethylene is shown to have a narrow field of stability. We report the p-T-x phase diagram of the C-H system and p-T phase diagram of CH4.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Structure, Stability and Mechanical Properties of Boron-Rich Mo-B Phases: A Computational Study

Molybdenum borides were studied theoretically using first-principles calculations, empirical total energy model and global optimization techniques to determine stable crystal structures. Our calculations reveal the structures of known Mo-B phases, attaining close agreement with experiment. Following our developed lattice model, we describe in detail the crystal structure of boron-rich $MoB_x$ phases with 3<x<9 as the hexagonal $P6_3/mmc$-$MoB_3$ structure with Mo atoms partially replaced by triangular boron units. The most energetically stable arrangement of these $B_3$ units corresponds to their uniform distribution in the bulk of the crystal structure, which leads to the formation of a disordered nonstoichiometric phase, with ordering arising at compositions close to x=5 due to a strong repulsive interaction between neighboring $B_3$ units. The most energetically favorable structures of $MoB_x$ correspond to the compositions 4<x<5, with $MoB_5$ being the boron-richest stable phase. The estimated hardness of $MoB_5$ is 37-39 GPa, suggesting that the boron-rich phases are potentially superhard.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Superconductivity of $LaH_{10}$ and $LaH_{16}$ polyhydrides

We explore high-pressure phase stability and superconductivity of lanthanum hydrides $LaH_m$ (m=4-11,16). We predict stability of a hitherto unreported polyhydride $P6/mmm$-$LaH_{16}$ at pressures above 150 GPa; at 200 GPa its predicted superconducting $T_C$ is 156 K, critical field $μ_0$$H_C$(0) ~ 35 T and superconducting gap is up to 35 meV. We revisit superconductivity of the recently discovered $LaH_{10}$ and find its $T_C$ to be up to 259 K (170 GPa) from solving the Eliashberg equation and 271 K from solving the gap equation in SCDFT which also allowed us to compute the Coulomb pseudopotential $μ^*$ for $LaH_{10}$ and $LaH_{16}$. Presence of several polymorph modifications of LaH10 may explain the variety in the experimentally measured $T_C$ values for $LaH_{10}$ [1,2]

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Theoretical study of the pressure-induced structure, phase transition, mechanical and electronic properties in the V-N system

Stable compounds in the V-N system are systematically searched and four new high-pressure phases are found, including C2/m-V$_9$N, Pbam-V$_5$N$_2$, Pnma-V$_2$N and I4/mcm-VN$_2$. V$_2$N undergoes a phase transition from $\varepsilon$-Fe$_2$N-type V$_2$N (P$\bar{3}$1m) to $ζ$-Fe$_2$N-type V$_2$N (Pbcn) at 10 GPa and to Fe$_2$C-type V$_2$N (Pnnm) at 59 GPa, then to Pnma-V$_2$N at 96 GPa. Low-temperature tetragonal VN is theoretically proved to belong to space group P$\bar{4}$2m. The estimated Vickers hardnesses and fracture toughness of WC-type VN are around 37 GPa and 4.3-6.1 MPa m$^{1/2}$, respectively. Al$_2$Cu-type VN$_2$ (I4/mcm) with a Vickers hardness of 25-27 GPa and fracture toughness of 3.6-6.6 MPa m$^{1/2}$ also shows excellent mechanical properties. Elastic properties of WC-type mononitrides of transition metals from IVB group (Ti, Zr and Hf), VB group (V, Nb and Ta) and VIB (Cr, Mo and W) are calculated and compared. Both the bond strength and structural configuration determine the mechanical properties of a material.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Discovery of new boron-rich chalcogenides: orthorhombic B6X (X = S, Se)

New boron-rich sulfide B6S and selenide B6Se have been discovered from high pressure - high temperature synthesis combined with ab initio evolutionary crystal structure prediction, and studied by synchrotron X-ray diffraction and Raman spectroscopy at ambient conditions. As it follows from Rietveld refinement of powder X-ray diffraction data, both chalcogenides have orthorhombic symmetry and belongs to Pmna space group. All experimentally observed Raman bands have been attributed to the theoretically calculated phonon modes, and the mode assignment has been performed. Prediction of mechanical properties (hardness and elastic moduli) of new boron-rich chalcogenides have been made using ab initio routines, and both compounds were found to be members of a family of hard phases (Hv ~ 31 GPa).

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

High-Temperature Superconductivity in the Ti--H System at High Pressures

Search for stable high-pressure compounds in the Ti--H system reveals the existence of titanium hydrides with new stoichiometries, including Ibam-Ti$_2$H$_5$, I4/m-Ti$_5$H$_{13}$, I$\bar{4}$-Ti$_5$H$_{14}$, Fddd-TiH$_4$, Immm-Ti$_2$H$_{13}$, P$\bar{1}$-TiH$_{12}$, and C2/m-TiH$_{22}$. Our calculations predict I4/mmm $\rightarrow$ R$\bar{3}$m and I4/mmm $\rightarrow$ Cmma transitions in TiH and TiH$_2$, respectively. Phonons and the electron--phonon coupling of all searched titanium hydrides are analyzed at high pressure. It is found that Immm-Ti$_2$H$_{13}$ rather than the highest hydrogen content C2/m-TiH$_{22}$, exhibits the highest superconducting critical temperature T$_{c}$. The estimated T$_{c}$ of Immm-Ti$_2$H$_{13}$ and C2/m-TiH$_{22}$ are respectively 127.4--149.4 K ($μ^{*}$=0.1-0.15) at 350 GPa and 91.3--110.2 K at 250 GPa by numerically solving the Eliashberg equations. One of the effects of pressure on T$_{c}$ can be attributed to the softening and hardening of phonons with increasing pressure.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

On Distribution of Superconductivity in Metal Hydrides

Using the data on the superconducting critical temperature ($T_{C}$) for a number of metal hydrides, we found a rule that makes it possible to predict the maximum $T_{C}$ based only on the information about the electronic structure of metal atoms. Using this guiding principle, we explored the hydride systems for which no reliable information existed, predicted new higher hydrides in the K-H, Zr-H, Hf-H, Ti-H, Mg-H, Sr-H, Ba-H, Cs-H, and Rb-H systems at high pressures, and calculated their $T_{C}$. Results of the study of actinides and lanthanides show that they form highly symmetric superhydrides $XH_{7-9}$. However, actinide hydrides do not exhibit high-temperature superconductivity (except Th-H system) and might not be considered as promising materials for experimental studies, as well as all $d^m$-elements with m > 4, including metal hydrides of the noble elements. Designed neural network allowing the prediction of $T_{C}$ of various hydrides shows good accuracy and was used to estimate upper limit for $T_{C}$ of the materials with absence of the data. The developed rule, based on regular behavior of the maximum achievable critical temperature as a function of number of $d+f$ electrons, enables targeted predictions about the existence of new high-$T_{C}$ superconductors.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Superconducting Praseodymium Superhydrides

Superhydrides have complex hydrogenic sublattices and are important prototypes for studying metallic hydrogen and high-temperature superconductors. Encouraged by the results on LaH10, in consideration of the differences between La and Pr, Pr-H system is especially worth studying because of the magnetism and valence-band f-electrons in element Pr. Here we successfully synthesized praseodymium superhydrides (PrH9) in laser-heated diamond anvil cells. Synchrotron X-ray diffraction (XRD) analysis demonstrated the presence of previously predicted F43m-PrH9 and unexpected P63/mmc-PrH9 phases. Moreover, Fm3m-PrH3, P4/nmm-PrH(3-δ) and Fm3m-PrH(1+x) were found below 52 GPa. F43m-PrH9 and P63/mmc-PrH9 were stable above 100 GPa in experiment. Experimental studies of electrical resistance in the PrH9 sample showed the emergence of superconducting transition (Tc) below 9K and a dependent Tc on applied magnetic field. Theoretical calculations indicate that magnetic order and electron-phonon interaction coexist in a very close range of pressures in the PrH9 sample which may contribute to its low superconducting temperature Tc. Our results highlight the intimate connections among hydrogenic sublattices, density of states, magnetism and superconductivity in Pr-based superhydrides.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Diverse Chemistry of Stable Hydronitrogens, and Implications for Planetary and Materials Sciences

Nitrogen hydrides, including ammonia (NH3), hydrazine (N2H4), hydrazoic acid (HN3) and etc, are compounds of great fundamental and applied importance. Their high-pressure behavior is important because of their abundance in giant planets and because of the hopes of discoverying high-energy-density materials. Here, we have performed a systematic investi- gation on the structural stability of N-H system in a pressure range up to 800 GPa through evolutionary structure prediction simulations. Surprisingly, we found that high pressure stabilizes a series of previously unreported compounds with peculiar structural and electronic properties, such as the N4H, N3H, N2H and NH phases composed of nitrogen backbones, the N9H4 phase containing two dimensional metallic nitrogen planes and novel N8H, NH2, N3H7, NH4 and NH5 molecular phases. Another surprise is that NH3 becomes thermodynamically unstable above ~460 GPa. We found that high-pressure chemistry is much more diverse that hydrocarbon chemistry at normal conditions, leading to expectations that N-H-O and N-H-O-S systems under pressure are likely to possess richer chemistry than the known organic chemistry. This, in turn, opens a possibility of nitrogen-based life at high pressure. The predicted phase diagram of the N-H system also provides a reference for synthesis of high-energy-density materials.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Hydrogen sulfide at high pressure: change in stoichiometry

Hydrogen-sulfide (H2S) was studied by x-ray synchrotron diffraction (XRD) and Raman spectroscopy up to 150 GPa at 180-295 K and by quantum-mechanical variable-composition evolutionary simulations. The experiments show that H2S becomes unstable with respect to formation of new compounds with different structure and composition, including Cccm and a body-centered-cubic (bcc) like (R3m or Im-3m) H3S, the latter one predicted previously to show a record-high superconducting transition temperature, Tc of 203 K. These experiments provide experimental ground for understanding of this record high Tc. The experimental results are supported by theoretical structure searches that suggest the stability of new H3S, H4S3, H5S8, H3S5, and HS2 compounds at elevated pressures.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Prediction of a new ground state of superhard compound B6O at ambient conditions

Boron suboxide B6O, the hardest known oxide, has an R-3m crystal structure (α-B6O) that can be described as an oxygen-intercalated structure of α-boron, or, equivalently, as a cubic close packing of B12 icosahedra with two oxygen atoms occupying all octahedral voids in it. Here we show a new ground state of this compound at ambient conditions, Cmcm-B6O (\b{eta}-B6O), which in all quantum-mechanical treatments that we tested (GGA, LDA, and hybrid functional HSE06) comes out to be slightly but consistently more stable. Increasing pressure and temperature further stabilize it with respect to the known α-B6O structure. \b{eta}-B6O also has a slightly higher hardness and may be synthesized using different experimental protocols. We suggest that \b{eta}-B6O is present in mixture with α-B6O, and its presence accounts for previously unexplained bands in the experimental Raman spectrum.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Pressure-induced stabilization of carbonic acid and other compounds in the C-H-O system

The physicochemical behavior of elements and compounds is heavily altered by high pressure. The occurrence of pressure-induced reactions and phase transitions can be revealed by crystal structure prediction approaches. In this work, we explore the C-H-O phase diagram up to 400 GPa exploiting an evolutionary algorithm for crystal structure predictions along with ab initio calculations. Besides uncovering new stable polymorphs of high-pressure elements and known molecules, we predicted the formation of new compounds. A 2CH4:3H2 inclusion compound forms at low pressure and remains stable up to 215 GPa. Carbonic acid (H2CO3), highly unstable at ambient conditions, was predicted to form exothermically at mild pressure (about 1 GPa). As pressure rises, it polymerizes and, above 300 GPa, reacts with water to form orthocarbonic acid (H4CO4). This unexpected high-pressure chemistry is rationalized by analyzing charge density and electron localization function distributions, and implications for general chemistry and planetary science are also discussed.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Antiferromagnetic Stabilization in Ti8O12

Using the evolutionary algorithm USPEX and DFT+U calculations, we predicted a high-symmetry geometric structure of bare Ti8O12 cluster composed of 8 Ti atoms forming a cube, which O atoms are at midpoints of all of its edges, in excellent agreement with experimental results. Using Natural Bond Orbital analysis, Adaptive Natural Density Partitioning algorithm, electron localization function and partial charge plots, we find the origin of the particular stability of bare Ti8O12 cluster: unique chemical bonding where eight electrons of Ti atoms interacting with each other in antiferromagnetic fashion to lower the total energy of the system. The bare Ti8O12 is thus an unusual molecule stabilized by d-orbital antiferromagnetic coupling.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Boron phosphide under pressure: in situ study by Raman scattering and X-ray diffraction

Cubic boron phosphide BP has been studied in situ by X-ray diffraction and Raman scattering up to 55 GPa at 300 K in a diamond anvil cell. The bulk modulus of B0 = 174(2) GPa has been established, which is in excellent agreement with our ab initio calculations. The data on Raman shift as a function of pressure, combined with equation-of-state data, allowed us to estimate the Grüneisen parameters of the TO and LO modes of zinc-blende structure, γGTO = 1.16 and γGLO = 1.04, just like in the case of other AIIIBV diamond-like phases, for which γGTO > γGLO = 1. We also established that the pressure dependence of the effective electro-optical constant α is responsible for a strong change in relative intensities of the TO and LO modes from ITO/ILO ~0.25 at 0.1 MPa to ITO/ILO ~2.5 at 45 GPa, for which we also find excellent agreement between experiment and theory.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Effect of carbon vacancies on structural and mechanical properties of stable zirconium carbides: A first principles study

By using evolutionary algorithm USPEX, we have predicted a number of stable zirconium carbides. In addition to the well-known rocksalt-type stoichiometric ZrC (Fm-3m), present prediction also identifies five stable substoichiometric zirconium carbides adopting rocksalt-type structures with ordered carbon vacancies, Zr8C7 (P-1), Zr6C5 (C2/m), Zr5C4 (P-1), Zr3C2 (C2/m), and Zr2C (Fd-3m). The effects of carbon vacancies on structural and mechanical properties are investigated. We highlight that the distribution of carbon vacancies has significant influence on volume, Pugh's ratio, and hardness. We further propose that hardness can be enhanced by replacing carbon vacancies with suitable elements, in particular nitrogen and oxygen.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Effects of ferroelectric polarization on surface phase diagram: an evolutionary algorithm study of the BaTiO$_{3}$(001) surface

We have constructed the surface phase diagram of the BaTiO$_{3}$(001) surface by employing an evolutionary algorithm for surface structure prediction, where the ferroelectric polarization is included as a degree of freedom. Among over 1000 candidate structures explored, a surface reconstruction of (2$\times$1)-TiO is discovered to be thermodynamically stable and have the \emph{p2mm} plane group symmetry as observed experimentally. We find that the influence of ferroelectric polarization on the surface free energy can be either negligibly small or sizably large (over 1 eV per ($2 \times 1$) supercell), depending strongly on the surface structure and resulting in a significant distinction of surface phase diagram with varying ferroelectric polarization. It is therefore feasible to control the surface stability by applying an external electric field. Our results may have important implications in understanding the surface reconstruction of ferroelectric materials and tuning surface properties.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Emergence of novel hydrogen chlorides under high pressure

HCl, a 'textbook' example of a polar covalent molecule, is a well-known compound of hydrogen and chlorine. Inspired by the discovery of unexpected stable stoichiometries of sodium chlorides, we performed systematic searches for all stable compounds in the H-Cl system from ambient pressure to higher pressures up to 500 GPa using variable-composition ab initio evolutionary algorithm USPEX. We found several compounds that are stable under pressure, i.e. HCl, H$_2$Cl, H$_3$Cl, H$_5$Cl and H$_4$Cl$_7$, which display a rich variety of chemical bonding types. At ambient pressure, H$_2$, Cl$_2$ and HCl molecular crystals are formed by weak intermolecular van der Waals interactions and adjacent HCl molecules connect with each other to form asymmetric zigzag chains, which become symmetric under high pressure. In hydrogen-rich chlorides, H$_2$ and HCl react to form the thermodynamically stable H$_3$Cl crystalline compound in which molecular cyclic H$_3^+$ cations are stabilised by the Cl$^-$ sublattice. Increasing the amount of hydrogen leads to stable solid-state H$_5$Cl, in which H$_2$ formally combines with H$_3^+$ to form H$_5^+$ cations. Additionally, chlorine-based Kagomé layers are formed with intercalated zigzag HCl chains in chlorine-rich hydrides. These discoveries help to understand how varied bonding features can co-exist and evolve in one compound under extreme conditions.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Exotic stable calcium carbides: theory and experiment

It is well known that pressure causes profound changes in the properties of atoms and chemical bonding, leading to the formation of many unusual materials. Here we systematically explore all stable calcium carbides at pressures from ambient to 100 GPa using variable-composition evolutionary structure predictions. We find that Ca5C2, Ca2C, Ca3C2, CaC, Ca2C3, and CaC2 have stability fields on the phase diagram. Among these, Ca2C and Ca2C3 are successfully synthesized for the first time via high-pressure experiments with excellent structural correspondence to theoretical predictions. Of particular significance are the base-centered monoclinic phase (space group C2/m) of Ca2C, a quasi-two-dimensional metal with layers of negatively charged calcium atoms, and the primitive monoclinic phase (space group P21/c) of CaC with zigzag C4 groups. Interestingly, strong interstitial charge localization is found in the structure of R-3m-Ca5C2 with semimetallic behaviour.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Formation of As-As Interlayer Bonding in the collapsed tetragonal phase of NaFe$_2$As$_2$ under pressure

NaFe$_2$As$_2$ is investigated experimentally using powder x-ray diffraction and Raman spectroscopy at pressures up to 23 GPa at room temperature and using ab-initio calculations. The results reveal a pressure-induced structural modification at 4 GPa from the starting tetragonal to a collapsed tetragonal phase. We determined the changes in interatomic distances under pressure that allowed us to connect the structural changes and superconductivity. The transition is related to the formation of interlayer As-As bonds at the expense of weakening of Fe-As bonds in agreement with recent theoretical predictions.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Influences of stoichiometry deviations on geometry and electronic structures of TiO2 anatase (101)

We systematically investigated native defects and reconstructions of TiO2 anatase (101) using USPEX and VASP codes. For lightly reduced samples, the coverage of O vacancies at surface depends on its concentration in bulk, while Ti interstitials prefer to dissolve in bulk. Atomic configurations of reconstructed TiO2 anatase (101) structures were given for severe reduced samples at the first time. Increasing the reduction degree, the Fermi level can be modulated in the whole bandgap. For O-rich samples, due to the indirect-to-direct bandgap transition and bandgap narrowing, O interstitials are responsible for visible light photoactivity. In a word, stoichiometry deviations of TiO2 anatase (101) have broad implications for photocatalytic reactions.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Nitrogen oxides under pressure stability, ionization, polymerization, and superconductivity

Nitrogen oxides are textbook class of molecular compounds, with extensive industrial applications. Nitrogen and oxygen are also among the most abundant elements in the universe. We explore the N-O system at 0 K and up to 500 GPa though ab initio evolutionary simulations. Results show that two phase transformations of stable molecular NO2 exist at 7 and 64 GPa, and followed by decomposition of NO2 at 91 GPa. All of the NO+NO3- structures are found to be metastable at T=0 K, so experimentally reported ionic NO+NO3- is either metastable or stabilized by temperature. Upon increasing pressure, N2O5 transforms from P-1 to C2/c structure at 51 GPa. NO becomes thermodynamically stable at 198 GPa. This polymeric phase is superconducting (Tc = 2.0 K) and contains a -N-N- backbone.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Novel lithium-nitrogen compounds at ambient and high pressures

Using ab initio evolutionary simulations, we predict the existence of five novel stable Li-N compounds at pressures from 0 to 100 GPa (Li13N, Li5N, Li3N2, LiN2, and LiN5). Structures of these compounds contain of isolated N ions, N2 dimers, polyacetylene-like N chains and N5 rings, respectively. The structure of Li13N consists of Li atoms and Li12N icosahedra (with N atom in the center of the Li12 icosahedron) such icosahedra are notdescribed by Wade-Jemmis electron counting rules and are unique. Electronic structure of Li-N compounds is found to dramatically depend on composition and pressure, making this system ideal for studying metal-insulator transitions. For example, LiN3 undergoes a sequence of pressure-induced transitions: metal-insulator-metal-insulator. This work resolves the previous controversies of theory and experiment on Li2N2.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Novel phase of beryllium fluoride at high pressure

A previously unknown thermodynamically stable high-pressure phase of BeF$_{2}$ has been predicted using the evolutionary algorithm USPEX. This phase occurs in the pressure range 18-27 GPa. Its structure has $C2/c$ space group symmetry and contains 18 atoms in the primitive unit cell. Given the analogy between BeF$_{2}$ and SiO$_{2}$, silica phases have been investigated as well, but the new phase has not been observed to be thermodynamically stable for this system. However, it is found to be metastable and to have comparable energy to the known metastable phases of SiO$_{2}$, suggesting a possibility of its synthesis.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Phagraphene: A Low-energy Graphene Allotrope composed of 5-6-7 Carbon Rings with Distorted Dirac Cones

Using systematic evolutionary structure searching we propose a new carbon allotrope, phagraphene, standing for penta-hexa-hepta-graphene, because the structure is composed of 5-6-7 carbon rings. This two-dimensional (2D) carbon structure is lower in energy than most of the predicted 2D carbon allotropes due to its sp2-hybridization and density of atomic packing comparable to graphene. More interestingly, the electronic structure of phagraphene has distorted Dirac cones. The direction-dependent cones are further proved to be robust against external strain with tunable Fermi velocities.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Prediction of new thermodynamically stable aluminum oxides

Recently, it has been shown that under pressure, unexpected and counterintuitive chemical compounds become stable. Laser shock experiments (A. Rode, unpublished) on alumina (Al2O3) have shown non-equilibrium decomposition of alumina with the formation of free Al and a mysterious transparent phase. Inspired by these observations, with have explored the possibility of the formation of new chemical compounds in the system Al-O. Using the variable-composition structure prediction algorithm USPEX, in addition to the well-known Al2O3, we have found two extraordinary compounds Al4O7 and AlO2 to be thermodynamically stable in the pressure range 330-443 GPa and above 332 GPa, respectively. Both of these compounds at the same time contain oxide O2- and peroxide O22- ions, and both are insulating. Peroxo-groups are responsible for gap states, which significantly reduce the electronic band gap of both Al4O7 and AlO2.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Prediction of Novel Stable 2D-Silicon with Fivefold Coordination

Silicene, an analogue of graphene, was so far predicted to be the only two-dimensional silicon (2D-Si) with massless Dirac fermions. Here we predict a brand new 2D-Si Dirac semimetal, which we name siliconeet [silik'ni:t]. Unexpectedly, it has a much lower energy than silicene and robust direction-dependent Dirac cones with Fermi velocities comparable to those in graphene. Remarkably, its peculiar structure based on pentagonal rings and fivefold coordination plays a critical role in the novel electronic properties. Taking spin-orbit coupling into account, siliconeet can also be recognized as a 2D-topological insulator with a larger nontrivial band gap than silicene.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Prediction of novel stable compounds in the Mg-Si-O system under exoplanet pressures

The Mg-Si-O system is the major Earth and rocky planet-forming system. Here, through quantum variable-composition evolutionary structure explorations, we have discovered several unexpected stable binary and ternary compounds in the Mg-Si-O system. Besides the well-known SiO2 phases, we have found two extraordinary silicon oxides, SiO3 and SiO, which become stable at pressures above 0.51 TPa and 1.89 TPa, respectively. In the Mg-O system, we have found one new compound, MgO3, which becomes stable at 0.89 TPa. We find that not only the (MgO)x(SiO2)y compounds, but also two (MgO3)x(SiO3)y compounds, MgSi3O12 and MgSiO6, have stability fields above 2.41 TPa and 2.95 TPa, respectively. The highly oxidized MgSi3O12 can form in deep mantles of mega-Earths with masses above 20 M+ (M+:Earth's mass). Furthermore, the dissociation pathways of pPv-MgSiO3 are also clarified, and found to be different at low and high temperatures. The low-temperature pathway is MgSiO3 -> Mg2SiO4 + MgSi2O5 -> SiO2 + Mg2SiO4 -> MgO + SiO2, while the high-temperature pathway is MgSiO3 -> Mg2SiO4 + MgSi2O5 -> MgO + MgSi2O5 -> MgO + SiO2. Present results are relevant for models of the internal structure of giant exoplanets, and for understanding the high-pressure behavior of materials.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Pressure induced novel compounds in the Hf-O system from first-principles calculations

Using first-principles evolutionary simulations, we have systematically investigated phase stability in the Hf-O system at pressure up to 120 GPa. New compounds Hf5O2, Hf3O2, HfO and HfO3 are discovered to be thermodynamically stable at certain pressure ranges and a new stable high-pressure phase is found for Hf2O with space group Pnnm and anti-CaCl2-type structure. Both P62m-HfO and P4m2-Hf2O3 show semimetallic character. Pnnm-HfO3 shows interesting structure, simultaneously containing oxide O2- and peroxide [O-O]2- anions. Remarkably, it is P62m-HfO rather than OII-HfO2 that exhibits the highest mechanical characteristics among Hf-O compounds. Pnnm-Hf2O, Imm2-Hf5O2, P31m-Hf2O and P4m2-Hf2O3 phases also show superior mechanical properties, these phases can be quenched to ambient pressure and their properties can be exploited.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Stable magnesium peroxide at high pressure

Rocky planets are thought to comprise compounds of Mg and O as these are among the most abundant elements, but knowledge of their stable phases may be incomplete. MgO is known to be remarkably stable to very high pressure and chemically inert under reduced condition of the Earth's lower mantle. However, in 'icy' gas giants as well as in exoplanets oxygen may be a more abundant constituent (Ref. 1,2). Here, using synchrotron x-ray diffraction in laser-heated diamond anvil cells, we show that MgO and oxygen react at pressures above 94 GPa and T = 2150 K with the formation of the theoretically predicted I4/mcm MgO2 (Ref.3). Raman spectroscopy detects the presence of a peroxide ion (O22-) in the synthesized material as well as in the recovered specimen. Likewise, energy-dispersive x-ray spectroscopy confirms that the recovered sample has higher oxygen content than pure MgO. Our finding suggests that MgO2 may substitute MgO in rocky mantles and rocky planetary cores under highly oxidizing conditions.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Superconductivity of novel tin hydrides (Sn$_n$H$_m$) under pressure

With the motivation of discovering high-temperature superconductors, evolutionary algorithm is employed to search for all stable compounds in the Sn-H system. In addition to the traditional SnH$_4$, new hydrides SnH$_8$, SnH$_{12}$ and SnH$_{14}$ are found to be thermodynamically stable at high pressure. Dynamical stability and superconductivity of tin-hydrides are systematically investigated. I$\bar{4}$m2-SnH$_8$, C2/m-SnH$_{12}$ and C2/m-SnH$_{14}$ exhibit higher superconducting transition temperatures of 81, 93 and 97 K compared to the traditional compound SnH$_4$ with T$_c$ of 52 K at 200 GPa. An interesting bent H$_3^-$ in I$\bar{4}$m2-SnH$_8$ and novel liner H$_4^-$ in C2/m-SnH$_{12}$ are observed. All the new tin-hydrides remain metallic over their predicted range of stability. The intermediate-frequency wagging and bending vibrations have more contribution to electron-phonon coupling parameter than high-frequency stretching vibrations of H$_2$ and H$_3$.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

The Phase Diagram and Hardness of Carbon Nitrides

Novel superhard materials, especially those with superior thermal and chemical stability, are needed to replace diamond. Carbon nitrides (C-N), which are likely to possess these charac- teristics and have even been expected to be harder than diamond, are excellent candidates. Here we report three new superhard and thermodynamically stable carbon nitride phases. Based on a systematic evolutionary structure searches, we report a complete phase diagram of the C-N system at 0-300 GPa and analyze the hardest metastable structures. Surprising- ly, we find that at zero pressure, the earlier proposed graphitic-C3N4 structure (P-6m2) is dynamically unstable and the lowest energy form of structures based on s-triazine unit and s-heptazine unit have similar topology and belong to the same space group.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Two-Dimensional Magnetic Boron

We predict a two-dimensional (2D) antiferromagnetic (AFM) boron (designated as M-boron) by using ab initio evolutionary methodology. M-boron is entirely composed of B20 clusters in a hexagonal arrangement. Most strikingly, the highest valence band of M-boron is isolated, strongly localized, and quite flat, which induces spin polarization on each cap of the B20 cluster. This flat band originates from the unpaired electrons of the capping atoms, and is responsible for magnetism. M-boron is thermodynamically metastable and is the first cluster-based 2D magnetic material in the elemental boron system.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Backbone NxH Compounds at High Pressures

Optical and synchrotron x-ray diffraction diamond anvil cell experiments have been combined with first principles theoretical structure predictions to investigate mixed N2 and H2 up to 55 GPa. We found the formation of oligomeric NxH (x>1) compounds using mechano- and photochemistry at pressures above 47 and 10 GPa, respectively, and room temperature. These compounds can be recovered to ambient pressure at T<130 K, whereas at room temperature, they can be metastable down to 3.5 GPa. Our results suggest new pathways for synthesis of environmentally benign high energy-density materials and alternative planetary ice.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Exploration of stable compounds, crystal structures, and superconductivity in the Be-H system

Using first-principles variable-composition evolutionary methodology, we explored the high-pressure structures of beryllium hydrides between 0 and 400 GPa. We found that BeH$_2$ remains the only stable compound in this pressure range. The pressure-induced transformations are predicted as $Ibam$ $\rightarrow $ $P\bar{3}m1$ $\rightarrow $ $R\bar{3}m$ $ \rightarrow $ $Cmcm$ $ \rightarrow $ $P4/nmm$, which occur at 24, 139, 204 and 349 GPa, respectively. $P\bar{3}m1$ and $R\bar{3}m$ structures are layered polytypes based on close packings of H atoms with Be atoms filling octahedral voids in alternating layers. $Cmcm$ and $P4/nmm$ structures have 3D-networks of strong bonds, but also feature rectanular and squre, respectively, layers of H atoms with short H-H distances. $P\bar{3}m1$ and $R\bar{3}m$ are semiconductors while $Cmcm$ and $P4/nmm$ are metallic. We have explored superconductivity of both metallic phases, and found large electron-phonon coupling parameters of $ λ$=0.63 for $Cmcm$ (resulting in a $T_c$ of 32.1-44.1 K) at 250 GPa and $ λ$=0.65 for $P4/nmm$ ($T_c$ = 46.1-62.4 K) at 400 GPa. The dependence of $T_c$ on pressure indicates that $T_c$ initially increases to a maximum of 45.1 K for $Cmcm$ at 275 GPa and 97.0 K for $P4/nmm$ at 365 GPa, and then decreases with increasing pressure for both phases.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Phase stability, chemical bonding and mechanical properties of titanium nitrides: A first-principles study

We have performed first-principles evolutionary searches for all stable titanium nitrides and have found, in addition to the well-known rocksalt-type TiN, new ground states Ti$_3$N$_2$, Ti$_4$N$_3$, Ti$_6$N$_5$ at atmospheric pressure, and Ti$_2$N and TiN$_2$ at higher pressures. The latest nitrogen-rich structure presents encapsulated N$_2$ dumbbells with a N-N distance of 1.348 Å at 60 GPa and TiN$_2$ is predicted to be mechanically stable (quenchable). Our calculations of the mechanical properties (bulk modulus, shear modulus, Young's modulus, Poisson's ratio, and hardness) are in excellent agreement with the available experimental data and show that the hardness of titanium nitrides increases with increasing nitrogen content. The hardness of titanium nitrides is enhanced by strengthening directional covalent bonds and disappearance of Ti-Ti metallic bonds. Among the predicted compounds, TiN$_2$ has the highest hardness of 27.2 GPa.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Predicted Novel Hydrogen Hydrate Structures under Pressure from First Principles

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.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Semimetallic Two-Dimensional Boron Allotrope with Massless Dirac Fermions

It has been widely accepted that planar boron structures, composed of triangular and hexagonal motifs are the most stable two dimensional (2D) phases and likely precursors for boron nanostructures. Here we predict, based on ab initio evolutionary structure search, novel 2D boron structure with non-zero thickness, which is considerably, by 50 meV/atom lower in energy than the recently proposed alpha-sheet structure and its analogues. In particular, this phase is identified for the first time to have a distorted Dirac cone, after graphene and silicene the third elemental material with massless Dirac fermions. The buckling and coupling between the two sublattices not only enhance the energetic stability, but also are the key factors for the emergence of the distorted Dirac cone.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Stability of numerous novel potassium chlorides at high pressure

K-Cl is a simple system displaying all four main types of bonding, as it contains (i) metallic potassium, (ii) elemental chlorine made of covalently bonded Cl2 molecules held together by van der Waals forces, and (iii) an archetypal ionic compound KCl. Well-known octet rule (and its special case, the charge balance rule, assigning classical charges of +1 to K and -1 to Cl) predicts that no compounds other than KCl are possible. However, our quantum-mechanical variable-composition evolutionary simulations predict new stable compounds. The phase diagram of the K-Cl system turned out to be extremely complicated, featuring new thermodynamically compounds K3Cl, K2Cl, K3Cl2, K4Cl3, K5Cl4, K3Cl5, KCl3, KCl7. Of particular interest are 2D-metallic homologs Kn+1Cln, the presence of positively charged Cl atoms in KCl7, and the predicted stability of KCl3 already at the atmospheric pressure.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Synthesis of ultra-incompressible sp3-hybridized carbon nitride

Search of materials with C-N composition hold a great promise in creating materials which would rival diamond in hardness due to the very strong and relatively low-ionic C-N bond. Early experimental and theoretical works on C-N compounds were based on structural similarity with binary A3B4 structural- types; however, the synthesis of C3N4 remains elusive. Here we explored an unbiased synthesis from the elemental materials at high pressures and temperatures. Using in situ synchrotron X-ray diffraction and Raman spectroscopy we demonstrate synthesis of highly incompressible Pnnm CN compound with sp3 hybridized carbon is synthesized above 55 GPa and 7000 K. This result is supported by first principles evolutionary search, which finds that Pnnm CN is the most stable compound above 10.9 GPa. On pressure release below 6 GPa the synthesized CN compound amorphizes reattaining its 1:1 stoichiometry as confirmed by Energy-Dispersive X-ray Spectroscopy. This work underscores the importance of understanding of novel high-pressure chemistry rules and it opens a new route for synthesis of superhard materials.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Unexpected Reconstruction of the alpha-Boron (111) Surface

We report on a novel reconstruction of the alpha-boron (111) surface, discovered using an ab initio evolution structure search, and reveal that it has an unexpected neat structure and much lower surface energy than the recently proposed (111)-I_R,(a) surface. For this reconstruction, every single interstitial boron atom forms bridges with the unique polar-covalent bonds between neighboring B_12 icosahedra, which perfectly meet the electron counting rule and are responsible for the reconstruction-induced metal-semiconductor transition. The peculiar charge transfer between the interstitial atoms and the icosahedra plays an important role in stabilizing the surface.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Evolutionary Method for Predicting Surface Reconstructions with Variable Stoichiometry

We present a specially designed evolutionary algorithm for the prediction of surface reconstructions. This new technique allows one to automatically explore all the low-energy configurations with variable surface atoms and variable surface unit cells through the whole chemical potential range. The power of evolutionary search is demonstrated by the efficient identification of diamond 2*1 (100) and 2*1 (111) surfaces with a fixed number of surface atom and a fixed cell size. With further variation of surface unit cells, we study the reconstructions of the polar surface MgO (111). Experiment has detected an oxygen trimer (ozone) motif (Plass et al, 1998). We predict a new version of this motif which can be thermodynamically stable at extreme oxygen rich condition. Finally, we perform a variable stoichiometry search for a complex ternary system: semi-polar GaN (101bar1) with and without adsorbed oxygen. The search yields a non-intuitive reconstruction based on N3-trimers. These examples demonstrate that an automated scheme to explore the energy landscape of surfaces will improve our understanding of surface reconstructions. The method presented in this report can be generally applied to binary and multi-component systems.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

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 fitness function in conjunction with first-principles calculations and 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 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 of discovering novel high-k dielectrics with both fixed and variable compositions, and will speed up the process of materials discovery.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Photovoltaic effect in BiFeO3/TiO2 heterostructures tuned with epitaxial strain and an electric field

The photovoltaic effect in the BiFeO3/TiO2 heterostructures can be tuned by epitaxial strain and an electric field in the visible-light region which is manifested by the enhancement of absorption activity in the heterojunction under tensile strain and an electric field based on the first-principles calculations. It is suggested that there are coupling between photon, spin carrier, charge, orbital, and lattice in the interface of the bilayer film which makes the heterojunction an intriguing candidate towards fabricating the multifunctional photoelectric devices based on spintronics. The microscopic mechanism involved in the heterostruces is related deeply with the spin transfer and charge rearrangement between the Fe 3d and O 2p orbitals in the vicinity of the interface.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Prediction of stable hafnium carbides: their stoichiometries, mechanical properties, and electronic structure

Hafnium carbides are studied by a systematic search for possible stable stoichiometric compounds in the Hf-C system at ambient pressure using variable-composition ab initio evolutionary algorithm implemented in the USPEX code. In addition to well-known HfC, we predicted two additional compounds Hf3C2 and Hf6C5. The structure of Hf6C5 with space group C2/m contains 11 atoms in the primitive cell and this prediction revives the earlier proposal by A. I. Gusev. The stable structure of Hf3C2 also has space group C2/m, and is more energetically favorable than the Immm, P-3m1, P2 and C2221 structures put forward by A. I. Gusev. Dynamical and mechanical stability of the newly predicted structures have been verified by calculations of their phonons and elastic constants. The bulk and shear moduli of Hf3C2 are 195.8 GPa and 143.1 GPa, respectively, while for Hf6C5 they are 227.9 GPa and 187.2 GPa, respectively. Their mechanical properties are inferior to those of HfC due to the presence of structural vacancies. Chemical bonding, band structure, and Bader charge are presented and discussed.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Unexpected stable stoichiometries of sodium chlorides

Sodium chloride (NaCl), or rocksalt, is well characterized at ambient pressure. Due to the large electronegativity difference between Na and Cl atoms, it has highly ionic chemical bonding, with stoichiometry 1:1 dictated by charge balance, and B1-type crystal structure. Here, by combining theoretical predictions and diamond anvil cell experiments we show that new materials with different stoichiometries emerge at pressure as low as 20 GPa. Compounds such us Na3Cl, Na2Cl, Na3Cl2, NaCl3 and NaCl7 are theoretically stable and have unusual bonding and electronic properties. To test this prediction, at 55-80 GPa we synthesized cubic and orthorhombic NaCl3 at 55-70 GPa and 2D-metallic tetragonal Na3Cl. This proves that novel compounds, violating chemical intuition, can be thermodynamically stable even in simplest systems at non-ambient conditions.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Constrained evolutionary algorithm for structure prediction of molecular crystals: methodology and applications

Evolutionary crystal structure prediction proved to be a powerful approach for studying a wide range of materials. Here, we present a specifically designed algorithm for the prediction of the structure of complex crystals consisting of well-defined molecular units. The main feature of this new approach is that each unit is treated as a whole body, which drastically reduces the search space and improves the efficiency, but necessitates the introduction of new variation operators described here. To increase diversity of the population of structures, the initial population and part($\scriptsize{\sim}$20%) of the new generations are generated using space group symmetry combined with random cell parameters and random positions and orientations of molecular units. We illustrate the efficiency and reliability of this approach by number of tests (ice, ammonia, carbon dioxide, methane, benzene, glycine and butane-1,4-diammonium dibromide). This approach easily predicts the crystal structure of methane \emph{A} containing 21 methane molecules (105 atoms) per unit cell. We demonstrate that this new approach has also a high potential for the study of complex inorganic crystals on the examples of a complex hydrogen storage material Mg(BH$_4$)$_2$ and elemental boron.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Evolutionary Metadynamics: a Novel Method to Predict Crystal Structures

A novel method for crystal structure prediction, based on metadynamics and evolutionary algorithms, is presented here. This technique can be used to produce efficiently both the ground state and metastable states easily reachable from a reasonable initial structure. We use the cell shape as collective variable and evolutionary variation operators developed in the context of the USPEX method [Oganov, Glass, \textit{J. Chem. Phys.}, 2006, \textbf{124}, 244704; Lyakhov \textit{et al., Comp. Phys. Comm.}, 2010, \textbf{181}, 1623; Oganov \textit{et al., Acc. Chem. Res.}, 2011, \textbf{44}, 227] to equilibrate the system as a function of the collective variables. We illustrate how this approach helps one to find stable and metastable states for Al$_2$SiO$_5$, SiO$_2$, MgSiO$_3$, and carbon. Apart from predicting crystal structures, the new method can also provide insight into mechanisms of phase transitions.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Fe-C and Fe-H systems at pressures of the Earth's inner core

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 Birch's 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.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

First-Principles Determination of the Structure of Magnesium Borohydride

The energy landscape of Mg(BH4)2 under pressure is explored by ab initio evolutionary calculations. Two new tetragonal structures, with space groups P-4 and I41/acd are predicted to be lower in enthalpy, by 15.4 kJ/mol and 21.2 kJ/mol, respectively, than the earlier proposed P42nm phase. We have simulated X-ray diffraction (XRD) spectra, lattice dynamics, and equations of state (EOS) of these phases. The density, volume contraction, bulk modulus, and the simulated XRD patterns of I41/acd and P-4 structures are in excellent agreement with the experimental results.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Stability of xenon oxides at high pressures

Xenon, which is quite inert under ambient conditions, may become reactive under pressure. The possibility of formation of stable xenon oxides and silicates in the interior of the Earth could explain the atmospheric missing xenon paradox. Using the ab initio evolutionary algorithm, we predict the thermodynamical stabilization of Xe-O compounds at high pressures (XeO, XeO2 and XeO3 at pressures above 83, 102 and 114 GPa, respectively). Our calculations indicate large charge transfer in these oxides, suggesting that large electronegativity difference and pressure are the key factors favoring the formation of xenon compounds. Xenon compounds in the Earth's mantle, however, cannot directly explain the missing xenon paradox: xenon oxides are unstable in equilibrium with metallic iron in the Earth's lower mantle, while xenon silicates are predicted to spontaneously decompose at all mantle pressures (<136 GPa). This does not preclude Xe atoms from being retained in defects of mantle silicates and oxides.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Systematic search for low-enthalpy sp3 carbon using evolutionary metadynamics

We present a systematic search for low-energy metastable superhard carbon allotropes by using the recently developed evolutionary metadynamics technique. It is known that cold compression of graphite produces an allotrope at 15-20 GPa. Here we look for all low-enthalpy structures accessible from graphite. Starting from 2H- or 3R-graphite and applying the pressure of 20 GPa, a large variety of intermediate $sp^3$ carbon allotropes were observed in evolutionary metadynamics simulation. Our calculation not only found all the previous proposed candidates for `superhard graphite', but also predicted two allotropes (\emph{X}-carbon and \emph{Y}-carbon) showing novel 5+7 and 4+8 topologies. These superhard carbon allotropes can be classified into five families based on 6 (diamond/lonsdaleite), 5+7 (\emph{M/W}-carbon), 5+7 (\emph{X}-carbon), 4+8 (bct C$_4$), and 4+8 (\emph{Y}-carbon) topologies. This study shows evolutionary metadynamics is a powerful approach both to find the global minima and systematically search for low-energy metastable phases reachable from given starting materials.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Understanding the nature of "superhard graphite"

Numerous experiments showed that on cold compression graphite transforms into a new superhard and transparent allotrope. Several structures with different topologies have been proposed for this phase. While experimental data are consistent with these models, the only way to solve this puzzle is to find which structure is kinetically easiest to form. Using state-of-the-art molecular-dynamics transition path sampling simulations, we investigate kinetic pathways of the pressure-induced transformation of graphite to various superhard candidate structures. Unlike hitherto applied methods for elucidating nature of superhard graphite, transition path sampling realistically models nucleation events necessary for physically meaningful transformation kinetics. We demonstrate that nucleation mechanism and kinetics lead to $M$-carbon as the final product. $W$-carbon, initially competitor to $M$-carbon, is ruled out by phase growth. Bct-C$_4$ structure is not expected to be produced by cold compression due to less probable nucleation and higher barrier of formation.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Unexpected stable stoichiometries of sodium chlorides

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 Pauling's 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.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Unexpected Stoichiometries in Mg-O System under High Pressure

Using ab initio evolutionary simulations, we explore all the possible stoichiometries for Mg-O system at pressures 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 conditions at which MgO2 and Mg3O2 might exist in planetary conditions.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Comment on Zarechnaya et al. Pressure-induced isostructural phase transformation in gamma-B28. Phys. Rev. B 82, 184111 (2010)

Zarechnaya et al. claimed an isostructural transformation in gamma-B28 at about 40 GPa; below which the phase is more compressible (B0=227 GPa) and above which less compressible (B0=281 GPa) than in previous experiments or theory. Here we wish to point out some interesting questions related to these claims. Summarizing briefly our points, the suggestion of an isostructural transformation is inconsistent with ab initio calculations and experiment, we see no physical mechanism that could be responsible for such a transformation, and the evidence for it is self-contradictory.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Evolutionary search for novel superhard materials: Methodology and applications to forms of carbon and TiO2

We have developed a method for prediction of the hardest crystal structures in a given chemical system. It is based on the evolutionary algorithm USPEX (Universal Structure Prediction: Evolutionary Xtallography) and electronegativity-based hardness model that we have augmented with bond-valence model and graph theory. These extensions enable correct description of the hardness of layered, molecular, and low-symmetry crystal structures. Applying this method to C and TiO2, we have (i) obtained a number of low-energy carbon structures with hardness slightly lower than diamond and (ii) proved that TiO2 in any of its possible polymorphs cannot be the hardest oxide, its hardness being below 17 GPa.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

The high-pressure phase of boron, γ-B28: disputes and conclusions of 5 years after discovery

γ-B28 is a recently established high-pressure phase of boron. Its structure consists of icosahedral B12 clusters and B2 dumbbells in a NaCl-type arrangement (B2)δ+(B12)δ- and displays a significant charge transfer δ~0.5- 0.6. The discovery of this phase proved essential for the understanding and construction of the phase diagram of boron. γ-B28 was first experimentally obtained as a pure boron allotrope in early 2004 and its structure was discovered in 2006. This paper reviews recent results and in particular deals with the contentious issues related to the equation of state, hardness, putative isostructural phase transformation at ~40 GPa, and debates on the nature of chemical bonding in this phase. Our analysis confirms that (a) calculations based on density functional theory give an accurate description of its equation of state, (b) the reported isostructural phase transformation in γ-B28 is an artifact rather than a fact, (c) the best estimate of hardness of this phase is 50 GPa, (d) chemical bonding in this phase has a significant degree of ionicity. Apart from presenting an overview of previous results within a consistent view grounded in experiment, thermodynamics and quantum mechanics, we present new results on Bader charges in γ-B28 using different levels of quantum-mechanical theory (GGA, exact exchange, and HSE06 hybrid functional), and show that the earlier conclusion about significant degree of partial ionicity in this phase is very robust.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

How to quantify energy landscapes of solids

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 structure 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 Pauling's fifth rule.

0 citations0 reviewsNo code linkedOpen access