Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userCrystal structure prediction using ab initio evolutionary techniques: principles and applications
138
0
0.0
(
0
)
Ask ChatGPT about the research
No Arabic abstract
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.
rate research
Read More
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.
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.
We develop a theoretical and computational framework to study polarons in semiconductors and insulators from first principles. Our approach provides the formation energy, excitation energy, and wavefunction of both electron and hole polarons, and takes into account the coupling of the electron or hole to all phonons. An important feature of the present method is that it does not require supercell calculations, and relies exclusively on electron band structures, phonon dispersions, and electron-phonon matrix elements obtained from calculations in the crystal unit cell. Starting from the Kohn-Sham (KS) equations of density-functional theory, we formulate the polaron problem as a variational minimization, and we obtain a nonlinear eigenvalue problem in the basis of KS states and phonon eigenmodes. In our formalism the electronic component of the polaron is expressed as a coherent superposition of KS states, in close analogy with the solution of the Bethe-Salpeter equation for the calculation of excitons. We demonstrate the power of the methodology by studying polarons in LiF and Li2O2. We show that our method describes both small and large polarons, and seamlessly captures Frohlich-type polar electron-phonon coupling and non-Frohlich coupling to acoustic and optical phonons. To analyze in quantitative terms the electron-phonon coupling mechanisms leading to the formation of polarons, we introduce spectral decompositions similar to the Eliashberg spectral function. We validate our theory using both analytical results and direct calculations on large supercells. This study constitutes a first step toward complete ab initio many-body calculations of polarons in real materials.
First-principles calculation predict that olivine Li4MnFeCoNiP4O16 has ferrotoroidic characteristic and ferrimagnetic configuration with magnetic moment of 1.56 muB per formula unit. The ferrotoroidicity of this material makes it a potential candidate for magnetoelectric materials . Based on the orbital-resolved density of states for the transtion-metal ions in Li4MnFeCoNiP4O16, the spin configuration for Mn2+,Fe3+,Co2+, and Ni2+ is t2g3eg2, t2g3eg2,t2g1t2g3eg1eg2, and t2g2t2g3eg1eg2, respectively. Density functional theory plus U (DFT+U) shows a indirect band gap of 1.25 eV in this predicted material, which is not simply related to the electronic conductivity in terms of being used as cathode material in rechargeable Li-ion batteries.
We present the implementation of GAtor, a massively parallel, first principles genetic algorithm (GA) for molecular crystal structure prediction. GAtor is written in Python and currently interfaces with the FHI-aims code to perform local optimizations and energy evaluations using dispersion-inclusive density functional theory (DFT). GAtor offers a variety of fitness evaluation, selection, crossover, and mutation schemes. Breeding operators designed specifically for molecular crystals provide a balance between exploration and exploitation. Evolutionary niching is implemented in GAtor by using machine learning to cluster the dynamically updated population by structural similarity and then employing a cluster-based fitness function. Evolutionary niching promotes uniform sampling of the potential energy surface by evolving several sub-populations, which helps overcome initial pool biases and selection biases (genetic drift). The various settings offered by GAtor increase the likelihood of locating numerous low-energy minima, including those located in disconnected, hard to reach regions of the potential energy landscape. The best structures generated are re-relaxed and re-ranked using a hierarchy of increasingly accurate DFT functionals and dispersion methods. GAtor is applied to a chemically diverse set of four past blind test targets, characterized by different types of intermolecular interactions. The experimentally observed structures and other low-energy structures are found for all four targets. In particular, for Target II, 5-cyano-3-hydroxythiophene, the top ranked putative crystal structure is a $Z^prime$=2 structure with P$bar{1}$ symmetry and a scaffold packing motif, which has not been reported previously.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
