Low energy paths for octahedral tilting in inorganic halide perovskites

Instabilities relating to cooperative octahedral tilting is common in materials with perovskite structures, and in particular in the sub class of halide perovskites. In this work, the energetics of octahedral tilting in the inorganic metal halide perovskites CsPbI$_3$ and CsSnI$_3$ are investigated using first-principles density functional theory calculations. Several low energy paths between symmetry equivalent variants of the stable orthorhombic (\textit{Pnma}) perovskite variant are identified and investigated. The results are in favor of the presence of dynamic disorder in the octahedral tilting phase transitions of inorganic halide perovskites. In particular, one specific type of path, corresponding to an out-of-phase "tilt switch", is found to have significantly lower energy barrier than the others, which indicates the existence of a temperature range where the dynamic fluctuations of the octahedra follow only this type of path. This could produce a time averaged structure corresponding to the intermediate tetragonal (\textit{P4/mbm}) structure observed in experiments. Deficiencies of the commonly employed simple one-dimensional "double well" potentials in describing the dynamics of the octahedra are pointed out and discussed.Comment: Revised versio

Phase stability of Fe from first-principles: atomistic spin dynamics coupled with ab initio molecular dynamics simulations and thermodynamic integration

The calculation of free energies from first principles in materials is a formidable task which enables the prediction of phase stability with high accuracy; these calculations are complicated in magnetic materials by the interplay of electronic, magnetic, and vibrational degrees of freedom. In this work, we show the feasibility and accuracy of the calculation of phase stability in magnetic systems with ab initio methods and thermodynamic integration by sampling the magnetic and vibrational phase space with coupled atomistic spin dynamics-ab initio molecular dynamics (ASD-AIMD) simulations [Stockem et al., PRL 121, 125902 (2018)], where energies and interatomic forces are calculated with density functional theory (DFT). We employ the method to calculate the phase stability of Fe at ambient pressure from 800 K up to 1800 K. The Gibbs free energy difference between fcc and bcc Fe at zero pressure as a function of temperature is calculated carrying out thermodynamic integration over temperature on the energies at the DFT level from ASD-AIMD, using a reference free energy difference calculated in the paramagnetic state at temperatures much higher than the magnetic transition temperatures with thermodynamic integration over stress-strain variables with disordered local moment (DLM)-AIMD simulations. We show the importance of the magnetic ordering temperature of bcc Fe on the $\alpha$ to $\gamma$ structural transition temperature, whereas the $\gamma$ to $\delta$ transition is well reproduced independently of the exchange interactions. The Gibbs free energy difference between the two structures is within 5 meV/atom from the CALPHAD estimate, and both transition temperatures are reproduced within 150 K. The present work paves the way to free energy calculations in magnetic materials from first principles with accuracy in the order of 1 meV/atom

Ab initio Determination of Phase Stabilities of Dynamically Disordered Solids: rotational C2 disorder in Li2C2

The temperature-induced orthorhombic to cubic phase transition in Li2C2is a prototypical ex-ample of a solid to solid phase transformation between an ordered phase, which is well describedwithin the phonon theory, and a dynamically disordered phase with rotating molecules, for which the standard phonon theory is not applicable. The transformation in Li2C2 happens from a phase with directionally ordered C2 dimers to a structure, where they are dynamically disordered. We provide a description of this transition within the recently developed method (Klarbring et al.,Phys.Rev. Lett. 121, 225702 (2018)) employing ab initio molecular dynamics (AIMD) based stress-strain thermodynamic integration on a deformation path that connects the ordered and dynamically disordered phases. The free energy difference between the two phases is obtained. The entropy that stabilizes the dynamically disordered cubic phase is captured by the behavior of the stress on the deformation path

Structural Dynamics Descriptors for Metal Halide Perovskites

Metal halide perovskites have shown extraordinary performance in solar energy conversion technologies. They have been classified as "soft semiconductors" due to their flexible corner-sharing octahedral networks and polymorphous nature. Understanding the local and average structures continues to be challenging for both modelling and experiments. Here, we report the quantitative analysis of structural dynamics in time and space from molecular dynamics simulations of perovskite crystals. The compact descriptors provided cover a wide variety of structural properties, including octahedral tilting and distortion, local lattice parameters, molecular orientations, as well as their spatial correlation. To validate our methods, we have trained a machine learning force field (MLFF) for methylammonium lead bromide (CH$_3$NH$_3$PbBr$_3$) using an on-the-fly training approach with Gaussian process regression. The known stable phases are reproduced and we find an additional symmetry-breaking effect in the cubic and tetragonal phases close to the phase transition temperature. To test the implementation for large trajectories, we also apply it to 69,120 atom simulations for CsPbI$_3$ based on an MLFF developed using the atomic cluster expansion formalism. The structural dynamics descriptors and Python toolkit are general to perovskites and readily transferable to more complex compositions.Comment: 10 figure

Imperfections are not 0 K: free energy of point defects in crystals

Defects determine many important properties and applications of materials, ranging from doping in semiconductors, to conductivity in mixed ionic-electronic conductors used in batteries, to active sites in catalysts. The theoretical description of defect formation in crystals has evolved substantially over the past century. Advances in supercomputing hardware, and the integration of new computational techniques such as machine learning, provide an opportunity to model longer length and time-scales than previously possible. In this Tutorial Review, we cover the description of free energies for defect formation at finite temperatures, including configurational (structural, electronic, spin) and vibrational terms. We discuss challenges in accounting for metastable defect configurations, progress such as machine learning force fields and thermodynamic integration to directly access entropic contributions, and bottlenecks in going beyond the dilute limit of defect formation. Such developments are necessary to support a new era of accurate defect predictions in computational materials chemistry

Na–Ni–H phase formation at high pressures and high temperatures: hydrido complexes [NiH5]3– versus the perovskite NaNiH3

The Na-Ni-H system was investigated by in situ synchrotron diffraction studies of reaction mixtures NaH-Ni-H-2 at around 5, 10, and 12 GPa. The existence of ternary hydrogen-rich hydrides with compositions Na3NiH5 and NaNiH3, where Ni attains the oxidation state II, is demonstrated. Upon heating at similar to 5 GPa, face-centered cubic (fcc) Na3NiH5 forms above 430 degrees C. Upon cooling, it undergoes a rapid and reversible phase transition at 330 degrees C to an orthorhombic (Cmcm) form. Upon pressure release, Na3NiH5 further transforms into its recoverable Pnma form whose structure was elucidated from synchrotron powder diffraction data, aided by first-principles density functional theory (DFT) calculations. Na3NiH5 features previously unknown square pyramidal 18- electron complexes NiH53-. In the high temperature fcc form, metal atoms are arranged as in the Heusler structure, and ab initio molecular dynamics simulations suggest that the complexes are dynamically disordered. The Heusler-type metal partial structure is essentially maintained in the low temperature Cmcm form, in which NiH53- complexes are ordered. It is considerably rearranged in the low pressure Pnma form. Experiments at 10 GPa showed an initial formation of fcc Na3NiH5 followed by the addition of the perovskite hydride NaNiH3, in which Ni(II) attains an octahedral environment by H atoms. NaNiH3 is recoverable at ambient pressures and represents the sole product of 12 GPa experiments. DFT calculations show that the decomposition of Na3NiH5 = NaNiH3 + 2 NaH is enthalpically favored at all pressures, suggesting that Na3NiH5 is metastable and its formation is kinetically favored. Ni-H bonding in metallic NaNiH3 is considered covalent, as in electron precise Na3NiH5, but delocalized in the polyanion [NiH3](-).Funding Agencies|Swedish Research Council (VR)Swedish Research Council [2019-05551]; Swedish Government Strategic Research Area in Materials Science on Advanced Functional Materials at at Linkoping University (Faculty Grant SFO-Mat-LiU) [200900971]; Carl Tryggers Stiftelse (CTS) [16:198, 17:206]</p

Dynamic Local Structure in Caesium Lead Iodide: Spatial Correlation and Transient Domains

Metal halide perovskites are multifunctional semiconductors with tunable structures and properties. They are highly dynamic crystals with complex octahedral tilting patterns and strongly anharmonic atomic behaviour. In the higher temperature, higher symmetry phases of these materials, several complex structural features have been observed. The local structure can differ greatly from the average structure and there is evidence that dynamic two-dimensional structures of correlated octahedral motion form. An understanding of the underlying complex atomistic dynamics is, however, still lacking. In this work, the local structure of the inorganic perovskite CsPbI$_3$ is investigated using a new machine learning force field based on the atomic cluster expansion framework. Through analysis of the temporal and spatial correlation observed during large-scale simulations, we reveal that the low frequency motion of octahedral tilts implies a double-well effective potential landscape, even well into the cubic phase. Moreover, dynamic local regions of lower symmetry are present within both higher symmetry phases. These regions are planar and we report the length and timescales of the motion. Finally, we investigate and visualise the spatial arrangement of these features and their interactions, providing a comprehensive picture of local structure in the higher symmetry phases

A first-principles non-equilibrium molecular dynamicsstudy of oxygen diffusion in Sm-doped ceria

Solid oxide fuel cells are considered as one of the main alternatives for future sources of clean energy. To further improve their performance, theoretical methods able to describe the diffusion process in candidate electrolyte materials at finite temperatures are needed. The method of choice for simulating systems at finite temperature is molecular dynamics. However, if the forces are calculated directly from the Schrödinger equation (first-principles molecular dynamics) the computational expense is too high to allow long enough simulations to properly capture the diffusion process in most materials. This thesis introduces a method to deal with this problem using an external force field to speed up the diffusion process in the simulation. The method is applied to study the diffusion of oxygen ions in Sm-doped ceria, which has showed promise in its use as an electrolyte. Good agreement with experimental data is demonstrated, indicating high potential for future applications of the method

A First-Principles Study of Highly Anharmonic and Dynamically Disordered Solids

This thesis is a first-principles theoretical investigation of solid materials with high degrees of anharmonicity. These are systems where the dynamics of the constituent atoms is too complex to be well-described by a set of uncoupled harmonic oscillators. While theoretical studies of such systems pose a significant challenge, they are under increasing demand due to the prevalence of these sytems in next-generation technological applications. Indeed, very anharmonic systems are ubiquitous in envisioned materials for future solid-state batteries and fuel-cells, thermoelectrics and optoelectronics. In some of these cases, the anharmonicity is a “side-effect” that simply has to be dealt with in order to accurately model certain properties, while in other cases the anharmonicity is the origin of the high-performance of the material. There are two main parts to the thesis: The first is on materials with perovskite-related structures. This is a very large class of materials, members of which are typically highly anharmonic, not least in relation to a series of complex phase transformations between different structural modifications. In the thesis, I have studied a specific class of such phase-transformations that relate to tilting of the framework of octahedra that make up the structure. The oxide CaMnO3 and a set of inorganic halide perovskites were taken as model systems. The results shed some light on the experimentally observed differences between the local and average atomic structure in these systems. I have further studied Cs2AgBiBr6, a member of the so-called lead-free halide double perovskites. I rationalized its temperature induced phase transformation and found high degrees of anharmonicity and ultra-low thermal conductivity. Finally, I studied the influence of nuclear quantum effects, which are often ignored in computational modelling, on the structure and bonding in the hybrid organic-inorganic lead-halide perovskite, CH3NH3PbI3. The second part of the thesis deals with theoretical studies of the phase stability of dynamically disordered solids. These are solids which have some sort of time-averaged long-range order, characteristic of a crystalline solid, but where the anharmonicity is so strong that the basic concept of an equilibrium atomic position cannot be statically assigned to all atoms. Examples include certain solids with very fast ionic conduction, so called superionics, and solids with rotating molecular units. This absence of equilibrium atomic positions makes many standard computational techniques to evaluate phase-stability inapplicable. I outline a method to deal with this issue, which is based on a stress-strain thermodynamic integration on a deformation path from an ordered variant to the dynamically disordered phase itself. I apply the method to study the phase stability of the high-temperature δ-phase of Bi2O3, which is the fastest know solid oxide ion conductor, and to Li2C2, whose high temperature cubic phase contains rotating C2 dimers. The thesis exemplifies the need to go beyond many of the standard approximations used in first-principles computational materials science if accurate theoretical predictions are to be made. This is true, in particular, for many emerging material classes in the field of energy materials.I den konventionella teoretiska modellen för ett (kristallint) fast material antags varje atom kunna tillordnas en jämviktsposition. Rörelsen av atomerna runt dessa jämviktspositioner antags sedan ofta vara harmoniskt, d.v.s. hyfsat kunna beskrivs i termer av en samling (kvantmekaniska) fjädrar. Denna avhandling behandlar teori- och beräkningsstudier av material där ett eller båda av dessa antaganden inte är giltiga, så kallade anharmoniska material. En nogrann teoretisk behandling av sådana material är ofta ordentligt utmanande. I takt med en snabb teknologiska utveckling ställs allt mer specifika och stränga krav på de material som behövs för diverse applikationer. Inom flertalet områden dyker då denna typ av komplexa och anharmoniska material upp som potentiella kandidater. Till exempel som fastelektrolyter för batterier och bränsleceller eller som solcellsmaterial. Inom vissa applikationer är denna anharmonicitet en bieffekt som man helt enkelt måste ta hänsyn till för att kunna göra noggranna teoretiska förutsägelser om diverse materialegenskaper, i andra fall är anharmoniciteten själva ursprunget för materialets goda egenskaper. I den första delen av avhandlingen behandlar jag material som har, eller är nära relaterade till, den så kallade perovskitstrukturen. Detta är en väldigt stor klass av material, och strukturen dyker därför upp inom en mängd olika tillämpningar, inte minst i lovande solcellsmaterial. Dessa material är ofta mycket anharmoniska, vilket tar sig uttryck bland annat i en rad komplexa fastransformationer mellan olika typer av perovskitmodifikationer. I perovskitoxiden CaMnO3, samt i en samling halogenperovskiter, har jag har studerat en specifik typ av fastransformationer som tillkommer på grund av rotationer av de oktaedrar som utgör en del av strukturen. Jag har fortsatt studerat den väldigt kraftiga anharmoniciteten i den så kallade blyfria halogendubbelperovskiten Cs2AgBiBr6, och slutligen har jag studerat hur en kvantmekanisk behandling av atomkärnorna, något som oftast inte görs, påverkar materialegenskaper i CH3NH3PbI3, en så kallad hybrid organisk-inorganisk bly-halogenperovskit, som är ett extremt lovande solcellsmaterial. I den andra delen av avhandlingen studerar jag dynamiskt oordnade fasta material. I dessa material är atomernas rörelse för komplex för att varje atom ska kunna tilldellas en statisk jämviktsposition. Material i denna klass är, till exempel, lovande som fastelektrolyter i bränsleceller och batterier. Mer specifikt studerar jag en typ av fasövergång, från en ordnad fas till en fas med dynamisk oordning, som ofta sker när dessa material värms upp. Jag introducerar en beräkningsmetod för att utvärdera deras fasstabilitet. Metoden är baserad på en så kallad termodynamisk integration, utförd mellan en ordnad variant och den dynamiskt oordnade fasen själv. Metoden gör det möjligt att beräkna fastransformationstemperaturer i denna typ av material. Jag applicerar metoden på Bi2O3, som i sin δ-fas är det fasta material med högst känd syrejonledningsförmåga, samt på Li2C2 vars kubiska fas innehåller roterande C2 molekyler. Resultaten stämmer bra överens med kända experimentella mätningar