A family of magnetic halide double perovskites (HDPs) have recently attracted attention due to their potential to broaden application areas of halide double perovskites into, e.g., spintronics. Up to date the theoretical modeling of these systems have relied on primitive approximations to the density functional theory (DFT). In this paper, we study structural, electronic and magnetic properties of the Fe3+-containing HDPs Cs2AgFeCl6 and Cs2NaFeCl6 using a combination of more advanced DFT-based methods, including DFT + U, hybrid-DFT, and treatments of various magnetic states. We examine the effect of varying the effective Hubbard parameter, U-eff, in DFT + U and the mixing-parameter, alpha, in hybrid DFT on the electronic structure and structural properties. Our results reveal a set of localized Fe(d) states that are highly sensitive to these parameters. Cs2AgFeCl6 and Cs2NaFeCl6 are both antiferromagnets with Neel temperatures well below room temperature and are thus in their paramagnetic (PM) state at the external conditions relevant to most applications. Therefore, we have examined the effect of disordered magnetism on the electronic structure of these systems and find that while Cs2NaFeCl6 is largely unaffected, Cs2AgFeCl6 shows significant renormalization of its electronic band structure.

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.

Phonon-phonon and electron/exciton-phonon coupling play a vitally important role in thermal, electronic, as well as optical properties of metal halide perovskites. In this work, we evaluate phonon anharmonicity and coupling between electronic and vibrational excitations in novel double perovskite Cs2NaFeCl6 single crystals. By employing comprehensive Raman measurements combined with first-principles theoretical calculations, we identify four Raman-active vibrational modes. Polarization properties of these modes imply Fm (3) over barm symmetry of the lattice, indicative for on average an ordered distribution of Fe and Na atoms in the lattice. We further show that temperature dependence of the Raman modes, such as changes in the phonon line width and their energies, suggests high phonon anharmonicity, typical for double perovskite materials. Resonant multiphonon Raman scattering reveals the presence of high-lying band states that mediate strong electron-phonon coupling and give rise to intense nA(1g) overtones up to the fifth order. Strong electron-phonon coupling in Cs2NaFeCl6 is also concluded based on the Urbach tail analysis of the absorption coefficient and the calculated Frohlich coupling constant. Our results, therefore, suggest significant impacts of phonon-phonon and electron-phonon interactions on electronic properties of Cs2NaFeCl6, important for potential applications of this novel material.

Superionic conductors have great potential as solid-stateelectrolytes,but the physics of type-II superionic transitions remains elusive.In this study, we employed molecular dynamics simulations, using machinelearning force fields, to investigate the type-II superionic phasetransition in & alpha;-Li3N. We characterized Li3N above and below the superionic phase transition by calculatingthe heat capacity, Li+ ion self-diffusion coefficient,and Li defect concentrations as functions of temperature. Our findingsindicate that both the Li+ self-diffusion coefficient andLi vacancy concentration follow distinct Arrhenius relationships inthe normal and superionic regimes. The activation energies for self-diffusionand Li vacancy formation decrease by a similar proportion across thesuperionic phase transition. This result suggests that the superionictransition may be driven by a decrease in defect formation energeticsrather than changes in Li transport mechanism. This insight may haveimplications for other type-II superionic materials.

The calculation of free energies from first principles enables the prediction of phase stability of materials 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 simulations [Stockem et al., PRL 121, 125902 (2018)], where energies and interatomic forces are calculated with density functional theory. We employ the method to calculate the phase stability of Fe at ambient pressure from 800 up to 1800 K. The Gibbs free energy difference between fcc and bcc Fe at zero pressure is calculated with thermodynamic integration over temperature and over stress-strain variables and, for the best set of exchange interactions employed, the Gibbs free energy difference between the two structures is within 5 meV/atom from the CALPHAD estimate, corresponding to an error in transition temperature below 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.

Lead-free halide double perovskites (HDPs) have emerged as a new generation of thermochromic materials. However, further materials development and mechanistic understanding are required. Here, a highly stable HDP Cs2NaFeCl6 single crystal is synthesized, and its remarkable and fully reversible thermochromism with a wide color variation from light-yellow to black over a temperature range of 10 to 423 K is investigated. First-principles, density functional theory (DFT)-based calculations indicate that the thermochromism in Cs2NaFeCl6 is an effect of electron–phonon coupling. The temperature sensitivity of the bandgap in Cs2NaFeCl6 is up to 2.52 meVK−1 based on the Varshni equation, which is significantly higher than that of lead halide perovskites and many conventional group-IV, III–V semiconductors. Meanwhile, this material shows excellent environmental, thermal, and thermochromic cycle stability. This work provides valuable insights into HDPs' thermochromism and sheds new light on developing efficient thermochromic materials.

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 modeling 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 (CH3NH3PbBr3) 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 CsPbI3 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.

Lead-based halide perovskite crystals are shown to have strongly anharmonic structural dynamics. This behavior is important because it may be the origin of their exceptional photovoltaic properties. The double perovskite, Cs2AgBiBr6, has been recently studied as a lead-free alternative for optoelectronic applications. However, it does not exhibit the excellent photovoltaic activity of the lead-based halide perovskites. Therefore, to explore the correlation between the anharmonic structural dynamics and optoelectronic properties in lead-based halide perovskites, the structural dynamics of Cs2AgBiBr6 are investigated and are compared to its lead-based analog, CsPbBr3. Using temperature-dependent Raman measurements, it is found that both materials are indeed strongly anharmonic. Nonetheless, the expression of their anharmonic behavior is markedly different. Cs2AgBiBr6 has well-defined normal modes throughout the measured temperature range, while CsPbBr3 exhibits a complete breakdown of the normal-mode picture above 200 K. It is suggested that the breakdown of the normal-mode picture implies that the average crystal structure may not be a proper starting point to understand the electronic properties of the crystal. In addition to our main findings, an unreported phase of Cs2AgBiBr6 is also discovered below approximate to 37 K.

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 CaMnO₃, 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 Cs₂AgBiBr₆, och slutligen har jag studerat hur en kvantmekanisk behandling av atomkärnorna, något som oftast inte görs, påverkar materialegenskaper i CH₃NH₃PbI₃, 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å Bi₂O₃, som i sin δ-fas är det fasta material med högst känd syrejonledningsförmåga, samt på Li₂C₂ vars kubiska fas innehåller roterande C₂ molekyler. Resultaten stämmer bra överens med kända experimentella mätningar.

The lead-free halide double perovskite class of materials offers a promising venue for resolving issues related to toxicity of Pb and long-term stability of the lead-containing halide perovskites. We present a first-principles study of the lattice vibrations in Cs2AgBiBr6, the prototypical compound in this class and show that the lattice dynamics of Cs2AgBiBr6 is highly anharmonic, largely in regards to tilting of AgBr6 and BiBr6 octahedra. Using an energy- and temperature-dependent phonon spectral function, we then show how the experimentally observed cubic-to-tetragonal phase transformation is caused by the collapse of a soft phonon branch. We finally reveal that the softness and anharmonicity of Cs2AgBiBr6 yield an ultralow thermal conductivity, unexpected of high-symmetry cubic structures.

Environmentally friendly halide double perovskites with improved stability are regarded as a promising alternative to lead halide perovskites. The benchmark double perovskite, Cs2AgBiBr6, shows attractive optical and electronic features, making it promising for high-efficiency optoelectronic devices. However, the large band gap limits its further applications, especially for photovoltaics. Herein, we develop a novel crystal-engineering strategy to significantly decrease the band gap by approximately 0.26 eV, reaching the smallest reported band gap of 1.72 eV for Cs(2)AgBiBr(6)under ambient conditions. The band-gap narrowing is confirmed by both absorption and photoluminescence measurements. Our first-principles calculations indicate that enhanced Ag-Bi disorder has a large impact on the band structure and decreases the band gap, providing a possible explanation of the observed band-gap narrowing effect. This work provides new insights for achieving lead-free double perovskites with suitable band gaps for optoelectronic applications.

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](-).

Lead-free halide double perovskites (A(2)B(I)B(III)X(6)) with attractive optical and electronic features are considered to be a promising candidate to overcome the toxicity and stability issues of lead halide perovskites (APbX(3)). However, their poor absorption profiles limit device performance. Here the absorption band edge of Cs(2)AgBiBr(6)double perovskite to the near-infrared range is significantly broadened by developing doped double perovskites, Cs-2(Ag:Cu)BiBr6. The partial replacement of Ag ions by Cu ions in the crystal lattice is confirmed by the X-ray photoelectron spectroscopy (XPS) and solid-state nuclear magnetic resonance (ssNMR) measurements. Cu doping barely affects the bandgap of Cs2AgBiBr6; instead it introduces subbandgap states with strong absorption to the near-infrared range. More interestingly, the near-infrared absorption can generate band carriers upon excitation, as indicated by the photoconductivity measurement. This work sheds new light on the absorption modulation of halide double perovskites for future efficient optoelectronic devices.

CaMn1-xNbxO3 (x = 0, 0.5, 0.6, 0.7 and 0.10) thin films have been grown by a two-step sputtering/annealing method. First, rock-salt-structured (Ca,Mn1-x,Nb-x)O thin films were deposited on 11 & x304;00 sapphire using reactive RF magnetron co-sputtering from elemental targets of Ca, Mn and Nb. The CaMn1-xNbxO3 films were then obtained by thermally induced phase transformation from rock-salt-structured (Ca,Mn1-xNbx)O to orthorhombic during post-deposition annealing at 700 degrees C for 3 h in oxygen flow. The X-ray diffraction patterns of pure CaMnO3 showed mixed orientation, while Nb-containing films were epitaxially grown in [101] out of-plane-direction. Scanning transmission electron microscopy showed a Ruddlesden-Popper (R-P) secondary phase in the films, which results in reduction of the electrical and thermal conductivity of CaMn1-xNbxO3. The electrical resistivity and Seebeck coefficient of the pure CaMnO3 film were measured to 2.7 omega cm and -270 mu V K-1 at room temperature, respectively. The electrical resistivity and Seebeck coefficient were reduced by alloying with Nb and was measured to 0.09 omega cm and -145 mu V K-1 for x = 0.05. Yielding a power factor of 21.5 mu W K-2 m(-1) near room temperature, nearly eight times higher than for pure CaMnO3 (2.8 mu W K-2 m(-1)). The power factors for alloyed samples are low compared to other studies on phase-pure material. This is due to high electrical resistivity originating from the secondary R-P phase. The thermal conductivity of the CaMn1-xNbxO3 films is low for all samples and is the lowest for x = 0.07 and 0.10, determined to 1.6 W m(-1) K-1. The low thermal conductivity is attributed to grain boundary scattering and the secondary R-P phase.

The structural and vibrational properties of ammonium borohydride, NH4BH4, have been examined by first-principles density functional theory (DFT) calculations and inelastic neutron scattering (INS). The H disordered crystal structure of NH4BH4 is composed of the tetrahedral complex ions NH4+ and BH4-, which are arranged as in the fcc NaCl structure and linked by intermolecular dihydrogen bonding. Upon cooling, the INS spectra revealed a structural transition between 45 and 40 K. The reversible transition occurs upon heating between 46 and 49 K. In the low-temperature form reorientational dynamics are frozen. The libration modes for BH4- and NH4+ are near 300 and 200 cm(-1), respectively. Upon entering the fcc high-temperature form, NH4+ ions attain fast reorientational dynamics, as indicated in the disappearance of the NH4+ libration band, whereas BH4- ions become significantly mobile only at temperatures above 100 K. The vibrational behavior of BH4- ions in NH4BH4 compares well to the heavier alkali metal borohydrides, NaBH4-CsBH4. DFT calculations revealed a nondirectional nature of the dihydrogen bonding in NH4BH4 with only weak tendency for long-range order. Different rotational configurations of complex ions appear quasi-degenerate, which is reminiscent of glasses.

Instabilities relating to cooperative octahedral tilting is common in materials with perovskite structures, in particular in the subclass of halide perovskites. In this work the energetics of octahedral tilting in the inorganic metal halide perovskites CsPbI3 and CsSnI3 are investigated using first-principles density functional theory calculations. Several low-energy paths between symmetry equivalent variants of the stable orthorhombic (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 (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.

Recent theoretical investigations [A. B. Belonoshko et aL Nat. Geosci. 10, 312 (2017)] revealed the occurrence of the concerted migration of several atoms in bcc Fe at inner-core temperatures and pressures. Here, we combine first-principles and semiempirical atomistic simulations to show that a diffusion mechanism analogous to the one predicted for bcc iron at extreme conditions is also operative and of relevance for the high-temperature bcc phase of pure Ti at ambient pressure. The mechanism entails a rapid collective movement of numerous (from two to dozens) neighbors along tangled closed-loop paths in defect-free crystal regions. We argue that this phenomenon closely resembles the diffusion behavior of superionics and liquid metals. Furthermore, we suggest that concerted migration is the atomistic manifestation of vanishingly small co-mode phonon frequencies previously detected via neutron scattering and the mechanism underlying anomalously large and markedly non-Arrhenius self-diffusivities characteristic of bcc Ti.

Lithium carbide, Li2C2, is a fascinating material that combines strong covalent and weak ionic bonding resulting in a wide range of unusual properties. The mechanism of its phase transition from the ground-state orthorhombic (Immm) to the high-temperature cubic (Fm (3) over barm) crystal structure is not well understood and here we elucidate it with help of first-principles calculations. We show that stabilization of the cubic phase is a result of a temperature-induced disorientation of the C-C dumbbells and their further thermal rotations. Due to these rotations rather large deviatoric stress, which is associated with the dumbbell alignment along one of the crystallographic axes, averages out making the cubic structure mechanically stable. At high temperature we observe a type-II superionic transition to a state of high Li self-diffusion involving collective ionic motion mediated by the formation of Frenkel pairs.

Lead-free halide double perovskites with diverse electronic structures and optical responses, as well as superior material stability show great promise for a range of optoelectronic applications. However, their large bandgaps limit their applications in the visible light range such as solar cells. In this work, an efficient temperature-derived bandgap modulation, that is, an exotic fully reversible thermochromism in both single crystals and thin films of Cs2AgBiBr6 double perovskites is demonstrated. Along with the thermochromism, temperature-dependent changes in the bond lengths of Ag Symbol of the Klingon Empire Br (R-Ag Symbol of the Klingon Empire Br) and Bi Symbol of the Klingon Empire Br (R-Bi Symbol of the Klingon Empire Br) are observed. The first-principle molecular dynamics simulations reveal substantial anharmonic fluctuations of the R-Ag Symbol of the Klingon Empire Br and R-Bi Symbol of the Klingon Empire Br at high temperatures. The synergy of anharmonic fluctuations and associated electron-phonon coupling, and the peculiar spin-orbit coupling effect, is responsible for the thermochromism. In addition, the intrinsic bandgap of Cs2AgBiBr6 shows negligible changes after repeated heating/cooling cycles under ambient conditions, indicating excellent thermal and environmental stability. This work demonstrates a stable thermochromic lead-free double perovskite that has great potential in the applications of smart windows and temperature sensors. Moreover, the findings on the structure modulation-induced bandgap narrowing of Cs2AgBiBr6 provide new insights for the further development of optoelectronic devices based on the lead-free halide double perovskites.

Ab initio molecular dynamics (AIMD) in combination with the temperature dependent effective potential (TDEP) method has been used to go beyond the quasiharmonic approximation and study the lattice dynamics in ceria, CeO2, at finite temperature. The results indicate that the previously proposed connection between the B-1u phonon mode turning imaginary and the transition to the superionic phase in fluorite structured materials is an artifact of the failure of the quasiharmonic approximation in describing the lattice dynamics at elevated temperatures. We instead show that, in the TDEP picture, a phonon mode coupling to the E-u mode prevents the B-1u mode from becoming imaginary. We directly observe the superionic transition at high temperatures in our AIMD simulations and find that it is initiated by the formation of oxygen Frenkel pairs (FP). These FP are found to form in a collective process involving simultaneous motion of two oxygen ions.

The temperature-induced antiferrodistortive (AFD) structural phase transitions in CaMnO3, a typical perovskite oxide, are studied using first-principles density functional theory calculations. These transitions are caused by tilting of the MnO6 octahedra that are related to unstable phonon modes in the high-symmetry cubic perovskite phase. Transitions due to octahedral tilting in perovskites normally are believed to fit into the standard soft-mode picture of displacive phase transitions. We calculate phonon-dispersion relations and potential-energy landscapes as functions of the unstable phonon modes and argue based on the results that the phase transitions are better described as being of order-disorder type. This means that the cubic phase emerges as a dynamical average when the system hops between local minima on the potential-energy surface. We then perform ab initio molecular dynamics simulations and find explicit evidence of the order-disorder dynamics in the system. Our conclusions are expected to be valid for other perovskite oxides, and we finally suggest how to predict the nature (displacive or order-disorder) of the AFD phase transitions in any perovskite system.

Theoretical studies of phase stability in solid materials with dynamic disorder are challenging due to the failure of the standard picture of atoms vibrating around fixed equilibrium positions. Dynamically disordered solid materials show immense potential in applications. In particular, superionic conductors, where the disorder results in exceptionally high ionic conductivity, are very promising as solid state electrolytes in batteries and fuel cells. The biggest obstacle in living up to this potential is the limited stability of the dynamically disordered phases. Here, we outline a method to obtain the free energy of a dynamically disordered solid. It is based on a stress-strain thermodynamic integration on a deformation path between a mechanically stable ordered variant of the disordered phase, and the dynamically disordered phase itself. We show that the large entropy contribution associated with the dynamic disorder is captured in the behavior of the stress along the deformation path. We apply the method to Bi2O3, whose superionic delta phase is the fastest known solid oxide ion conductor. We accurately reproduce the experimental transition enthalpy and the critical temperature of the phase transition from the low temperature ground state a phase to the superionic d phase. The method can be used for a first-principles description of the phase stability of superionic conductors and other materials with dynamic disorder, when the disordered phase can be connected to a stable phase through a continuous deformation path.

The MgMnH system was investigated by in situ high pressure studies of reaction mixtures MgH2MnH2. The formation conditions of two complex hydrides with composition Mg3MnH7 were established. Previously known hexagonal Mg3MnH7 (h-Mg3MnH7) formed at pressures 1.52 GPa and temperatures between 480 and 500 degrees C, whereas an orthorhombic form (o-Mg3MnH7) was obtained at pressures above 5 GPa and temperatures above 600 degrees C. The crystal structures of the polymorphs feature octahedral [Mn(I)H-6](5) complexes and interstitial H-. Interstitial H- is located in trigonal bipyramidal and square pyramidal interstices formed by Mg2+ ions in h- and o-Mg3MnH7, respectively. The hexagonal form can be retained at ambient pressure, whereas the orthorhombic form upon decompression undergoes a distortion to monoclinic Mg3MnH7 (m-Mg3MnH7). The structure elucidation of o- and m-Mg3MnH7 was aided by first-principles density functional theory (DFT) calculations. Calculated enthalpy versus pressure relations predict m- and o-Mg3MnH7 to be more stable than h-Mg3MnH7 above 4.3 GPa. Phonon calculations revealed o-Mg3MnH7 to be dynamically unstable at pressures below 5 GPa, which explains its phase transition to m-Mg3MnH7 on decompression. The electronic structure of the quenchable polymorphs h- and m-Mg3MnH7 is very similar. The stable 18-electron complex [MnH6](5-) is mirrored in the occupied states, and calculated band gaps are around 1.5 eV. The study underlines the significance of in situ investigations for mapping reaction conditions and understanding phase relations for hydrogen-rich complex transition metal hydrides.

A first-principles nonequilibrium molecular dynamics (NEMD) study employing the color-diffusion algorithm has been conducted to obtain the bulk ionic conductivity and the diffusion constant of gadolinium-doped cerium oxide (GDC) in the 850-1150 K temperature range. Being a slow process, ionic diffusion in solids usually requires simulation times that are prohibitively long for ab initio equilibrium molecular dynamics. The use of the color-diffusion algorithm allowed us to substantially speed up the oxygen-ion diffusion. The key parameters of the method, such as field direction and strength as well as color-charge distribution, have been investigated and their optimized values for the considered system have been determined. The calculated ionic conductivity and diffusion constants are in good agreement with available experimental data.

Sm-doped ceria is a prospective electrolyte material for intermediate-temperature solid-oxide fuel cells (IT-SOFC). Equilibrium ab initio molecular dynamics (AIMD) studies of oxygen ion diffusion in this material are currently impractical due to the rareness of diffusive events on the accessible timescale. To overcome this issue we have performed ab initio non-equilibrium molecular dynamics calculations of Sm-doped ceria using the color diffusion algorithm. Applying an external force field we have been able to increase the frequency of diffusive events over the simulation time, while keeping the physical mechanism of diffusion intact. We have investigated the temperature dependence of the maximum strength of the applied external field that could be used while maintaining the response of the system in a linear regime. This allows one to obtain the diffusivity at zero field. The bulk ionic conductivity has been calculated and found to match the experimental data well. We have also compared the description of the diffusion process by our method to previous findings and show that the migration mechanism and site preference of oxygen vacancies with respect to the Sm dopants is well reproduced. (C) 2016 Elsevier B.V. All rights reserved.