research papers
Classification of perovskite structural types with dynamical octahedral tilting
aUniversity of Bern, Bern, Switzerland, and bLaboratory for Waste Management, Paul Scherrer Institute, Villigen-PSI, Switzerland
*Correspondence e-mail: donat.adams@geo.unibe.ch
Perovskites ABX3 with delocalized positions of the X atoms represent a distinct class of dynamically distorted structures with peculiar structural relations and physical properties. The delocalization originates from atoms crossing shallow barriers of the surface. Quantum mechanically, they can be treated similar to light atoms in diffusive states. Many of these perovskite structures are widely used functional materials thanks to their particular physical properties, such as superconductivity, ferroelectricity and photo-activity. A number of these properties are related to static or dynamic motion of octahedral units. Yet, a full understanding of the relationships between perovskite chemical bonding and physical properties is currently missing. Several studies indicate the existence of dynamic disorder generated by anharmonic motion of octahedral units, e.g. in halide perovskite structures. To simplify structural analysis of such systems we derive a set of space groups for simple perovskites ABX3 with dynamical octahedral tilting. The derived space groups extend the well established tables for static tiltings by Glazer [Acta Cryst. B (1972). 28, 3384–3392], Aleksandrov [Ferroelectrics (1976). 24, 801–805] and Howard & Stokes [Acta Cryst. B (1998). 54, 782–789]. Ubiquity of dynamical tilting is demonstrated by an analysis of the structural data for perovskites reported in recent scientific publications and the signature of dynamic tilting in the corresponding structures is discussed, which can be summarized as follows: (a) volume increase upon a lowering of temperature, (b) apparent distortion of octahedra (where Jahn–Teller distortions can be ruled out), (c) mismatch between observed instantaneous symmetry and average symmetry, (d) deviation of the experimental from the theoretically predicted structures for static tilting, (e) inconsistency of lattice parameters with those suggested by the theory of static tilts, and (f) large displacement parameters for atoms at the X and B sites. Finally, the possible influence of dynamic disorder on the physical properties of halide perovskites is discussed.
Keywords: perovskites; symmetry mode analysis; structure prediction; hybrid materials; decoupled anharmonic mode approximation; dynamical simulations; octahedral tilting; dynamic disorder.
1. Introduction
Interatomic interactions determine the symmetry and stability of et al., 2005). Strong interplay between the motion of atomic cores and electronic structure is particularly relevant when it comes to phonon–electron coupling, which is responsible for many remarkable properties of materials such as ferroelectricity [e.g. SrTiO3 (Choudhury et al., 2008)], superconductivity (Bednorz & Müller, 1986, 1988) and photoactivity (Iaru et al., 2017; Zhao et al., 2019). These electron–phonon interactions are particularly well documented in hybrid halide perovskites (Munson et al., 2019). Phonon–electron coupling is held responsible for the fact that hybrid halide perovskites materials have particularly interesting applications, including solar cells with particularly high efficiency.
in minerals and synthetic materials. Specific structural arrangement of the atoms and their interaction are linked to the electron distribution and thermal motion of atoms commonly described by vibrational ellipsoids (DemetriouA classical picture of a solid assumes that the atomic nuclei are localized at well defined minima of the T → 0 K, the probability density of the nucleons is not point-like but is extended over a finite range (Heisenberg uncertainty principle). Thermal excitations lead to the statistical occupation of higher vibrational energy levels in the corresponding discrete energy spectrum. Eventually the atomic probability density is smeared out over a broad spatial range, conventionally represented as vibrational ellipsoids (see Fig. 2 for how dynamical tilts could explain large vibrational ellipsoids in structures with dynamical tilting).
In many systems, such a minimum can be well described by a harmonic approximation. The quantum mechanical treatment of solids reveals that, even in the limit ofThe amplitude of atomic motion depends on the shape of the e.g. in the case of hopping of hydrogen atoms (Völkl & Alefeld, 1978). The associated wavefunction is delocalized over a large domain. The surface (PES) in such a system has several local minima separated by low activation barriers, allowing instantaneous localization of atoms in different positions (Churakov & Wunder, 2004). In this article, we will consider the systems in which the delocalized state of atoms is responsible for some outstanding macroscopic properties. The physical origin of these phenomena is related to the quantum mechanical description of the atomic motion and requires clear distinction of the PES and total energy of the system, as defined in Appendix A1. We will also further develop this notion and consider atoms as delocalized.
minima and the excitation probability. For light atoms, the spatial extension of the nuclei can be larger,In crystal structures, the atomic dynamic can be considered as a collective motion of several atoms, rather than dynamics of individual atoms. The resulting collective
are vibro-rotational modes, often addressed as the eigenmodes of the dynamical matrix. We find it particularly appealing to consider rotational modes, which in some perovskite structures are delocalized over a multi-well PES. This can be inherent to some perovskite phases within certain cell geometries.In such structures, atoms are confined by a PES with multiple local energy minima (see e.g. Figs. 1 or 2) and are separated by small energy barriers [40 meV, for CH3NH3PbI3 (Beecher et al., 2016), or 20–100 meV in CsPbI3 and CsSnI3 (Klarbring, 2019)]. If these barriers are low compared with the thermal energy, the transition between these minima is enhanced and can be excited by small perturbations (Adams & Passerone, 2016).
Furthermore, we find that dynamic octahedral rotations are ubiquitous in halide perovskites and have actually already been observed, for example in CsPbCl3 (Fujii et al., 1974), CH3NH3PbBr3 and CH3NH3PbCl3 (Chi et al., 2005; Swainson et al., 2003, 2015) – see also Beecher et al. (2016) and Marronnier et al. (2017). This kind of disorder can be characterized as dynamic, i.e. atoms instantaneously occupy specific positions but on average these positions are only partially occupied. Related to the observations in the cubic phase of CH3NH3PbI3 through inelastic X-ray scattering, the term `dynamic disorder' has been coined (Poglitsch & Weber, 1987; Egger et al., 2018) [see also the review of Whalley et al. (2017)].
Gao et al. (2021) found `structural dynamics resembling that of a liquid', describing a dynamic in a PES with multiple energy minima linked through shallow energy barriers, which can be crossed even at low temperature, either due to thermal excitations or due to quantum tunnelling.
In this article, we will refer to these octahedral rotations as `dynamic tilting'. The dynamic instabilities are involved in the series of phase transitions in CH3NH3PbI3: from the orthorhombic phase (Pbn21 space group), to tetragonal structure (I4cm space group) at ∼160 K, and finally to the pseudo-cubic tetragonal phase at ∼330 K (P4mm space group) (Onoda-Yamamuro et al., 1992).
As previous calculations show (Adams & Oganov, 2006), such dynamical tiltings can lead to particular phonon-band structures with a high (DOS) at low energies – see Fig. 2(a). These vibrational modes are relevant for the electron–phonon interaction, possibly explaining the unusually high susceptibility of hybrid halide perovskites to external factors such as electric fields, temperature, humidity or mechanical stress (Lin et al., 2021; Ugur et al., 2020). The observed properties can be related to the peculiar shape of the PES. Interestingly, it is precisely these instabilities and the resulting dynamic tilting (Zhu & Ertekin, 2019; Zhang et al., 2015) that can lead to strong electron-scattering rates, large spectral width, and in turn might be the premise of the high efficiency of perovskite solar cells (Wright et al., 2016; Herz, 2017).
We therefore systematically explore symmetry constraints on possible dynamic rotations in the octahedra in perovskite structures. Then we explore the most important experimentally accessible signatures of dynamic tilting. Our hypothesis is that many perovskites exhibit dynamic tilting, especially in the group of halide perovskites. Accordingly, we review the most recent publications reporting dynamic tilting, in order to substantiate our hypothesis. Finally, we discuss the possible implications of dynamic tilting for the physical properties of perovskites such as CH3NH3PbI3, considering the electronic structure, vibrational properties (especially phonon dispersion), the electron–phonon coupling and efficiency in photovoltaics.
2. Methods
To derive possible tilting systems, a 2 × 2 × 2 ABO3 perovskite containing 40 atoms was set up, in which the rigid octahedra were tilted in three spatial directions. The Glazer notation (Glazer, 1972) has been used to denote an octahedral tilting by angle a, where a+ is for two succeeding octahedra tilted in phase and a− is for succeeding octahedra tilted out of phase. These tilts lead to a contraction of the two lattice parameters, which are perpendicular to the tilting axis – see also Fig. 3. This reduction has been implemented using three-dimensional rotation matrices. They preserve the bonding distance within the octahedron without further (collective) constraints on the structure.
of the ideal cubicAs explained in the Introduction, recent studies have lead to the discovery of dynamic tilts (Marronnier et al., 2017; Adams & Passerone, 2016). They result from octahedra confined in the structure by a PES with multiple energy minima – see Fig. 2. In this study, we only consider systems with symmetric arrangements of minima of the PES relative to a mirror plane. In these structures the tilt angle oscillates between positive and negative amplitude, spending most of the time at |a| > 0, whereas the average position is still a = 0 – see Fig. 2. However, due to the non-zero average tilt amplitude, time averages of cell parameters perpendicular to the tilting axis appear shorted compared with the untilted cubic structure. We denote these tiltings by ad and consider them by leaving the atomic positions unchanged and contracting the cell parameters perpendicular to the tilting axis.
2.1. Classification of dynamic tilt systems
The possible combinations of dynamic tilts and their combinations with static tilts were applied and the resulting structures were classified using the crystallographic tool FINDSYM (Stokes & Hatch, 2005). It determines the primitive lattice vectors, the of the lattice and thus the unique of the For static tilts, Howard & Stokes (1998) showed that the 25 distinct isotropy subgroups of static tilts can be reduced to 15. They considered only `simple' tilt systems, in which the tilts around a particular axis have the same magnitude and either the same sign (the + pattern) or alternating sign (the − pattern). This results in the exclusion of tilt systems that accidently show the same tilt angle but different types of tilting (+, −) on different axes. From application of this concept to the dynamic tilts, i.e. excluding tilting systems showing the same tilting angle but different types of tilting along different axes, 19 different tilting systems have been derived. Their space groups are reported in Table 1.
3. Applications and results
3.1. Volume changes across phase transitions and negative thermal expansion
Positive
is common for many materials, a property which is maintained even across phase transitions and linked to extended thermal movement upon temperature increase or more precisely to lowering of the free energy upon lattice expansion. A negative is rather unusual and needs clarification.The explanation of the negative
comes naturally when dynamical tilting is considered. Note that any (dynamic or static) tilting leads to a decrease in the unit-cell volume – as long as the octahedra are considered as rigid. This is because within this approximation the volume is proportional towhere a, b and c are tilting angles. Thus, Φ has a maximum at a = b = c = 0, and therefore the cubic Pmm structure usually denoted as the tilting system a0a0a0 should have the largest volume, which however is not always observed. In the presence of dynamical tilting, the cubic structure is rather attributed to the adadad tilting system. Therefore, the sharp structural at 327.15 K in CH3NH3PbI3 from a tetragonal phase to the cubic phase (Jacobsson et al., 2015), which is linked to a volume drop, can be understood as an activation of additional dynamic tiltings leading to a volume decrease according to equation (1). This explanation is further supported by calculations that show large negative portions of the band structure around the R and M points of the corresponding to octahedral tilting (Brivio et al., 2015; Akbarzadeh et al., 2005), and observations revealing large thermal movement (Tyson et al., 2017).
Similar observations are available in other structures, e.g. in KNbO3. There the volume change is positive across both phase transitions observed between 300 and 750 K (Sakakura et al., 2011), with the cubic Pmm structure formed at high temperature after a series of transitions.
In hybrid halide perovskite semiconductors the presence of organic ions generally reduces the crystal symmetry compared with simple perovskites containing monoatomic ions. Therefore, the space groups listed in Table 1 cannot be expected to occur there. However, some hybrid halide (Lehmann et al., 2019; Mante et al., 2018), such as perovskites, show a (directional) negative or volume change across phase transitions, which can be viewed as a fingerprint of the dynamical tilting.
3.2. Agreement between lattice parameters and tilts system
The lattice parameters published e.g. by Brivio et al. (2015) for methylammonium lead iodide are in agreement with the underlying tilt system. For the tetragonal I4/mcm phase, c/a > 1, while for the orthogonal phase, the calculated c/a < 1 is consistent with their respective tilt system. In other materials, such as Pr0.5Sr0.5MnO3, which is reported in the Fmmm and I4/mcm symmetry, static tilts cannot explain the pseudo-cubic c/a ratio of 1.02 (tilt system a−b0b0). Static tilt would require rotation of more than 10°, whereas the observed rotation angle is only 3.8°. Consistent results can be obtained by considering dynamic tilts, which leads to modification of the lattice parameters without changing crystal coordinates and resulting in the Fmmm space group.
The Fmmm structure is also reported for other Pr/Sr ratios (Knížek et al., 2004) and for combined Pr—Sr—Ce doping (Heyraud et al., 2013). The Pr sites retain their statistically distributed fractional occupancies and no ionic ordering takes place. It is tempting to explain the discrepancy between tilt system and cell parameters by assuming distortion of octahedral sites. This, however, would undermine the success of the theory of static tilts based on rigid octahedra. Considering dynamical tilting, it is obvious that non-zero average tilt amplitude reduces cell parameters perpendicular to the tilting axis, while the instantaneous geometry of the octahedron remains rigid.
3.3. Instantaneous atomic positions and time-averaged structure
Egger et al. (2018) discussed the possibility for dynamic tilting in halide perovskites (`dynamical disorder' in their terminology). The arguments for the presence of dynamic disorder were gathered by analysing structural results obtained by complementary experimental methods. While X-ray diffraction yields an averaged structure with high symmetry, Raman spectroscopy shows a local structure with low symmetry. Similarly, Beecher et al. (2016) report `anharmonic modes […] with diffusive (order–disorder) dynamics persisting many tens to hundreds of Kelvin above the transition'. Still, as has been stated, these modes are unobservable by Bragg diffraction. Indeed, these simultaneous observations can be reconciled through the concept of dynamical tilting: we interpret these findings as snapshots of the dynamic tilting. It has been stated (Kassan-Ogly & Naish, 1986) that the multi-well nature of the atomic PES cannot only lead to structural phase transitions to new phases, which is the main subject of this article, but also to diffuse scattering: above the phase-transition temperature the atoms retain some characteristics of the interactions below the (e.g. coupling of the octahedra), which leads to correlated movement of the octahedra.
Gao et al. (2021) experimentally found `structural dynamics resembling that of a liquid' for CsMBr3 (M = Pb, Sn, Ge) based on diffuse inelastic that increases towards 0 cm−1. This can be generated by atomic dynamics on a PES with multiple energy minima linked through shallow energy barriers. Due to quantum tunnelling or thermal excitations, atoms can cross these barriers, even at low temperature. After tunnelling, there is no classical force to restore the original configuration and the corresponding vibrational excitation frequencies thus tend towards 0 cm−1.
The transition rate may vary, depending on the height of the energy barriers between adjacent f and i are the final state and the initial state, respectively, ρ is the density of the final states and 〈f|V|i〉 is the matrix element connecting the two states]. A high energy barrier for tilting can result in sluggish dynamics, at a time scale significantly larger than diffraction experiments, making visible the tilting of the octahedra (through the diffuse scattering) in crystals of otherwise higher symmetry. This high symmetry appears in the time-averaged X-ray patterns, or equivalently, in the averaged structure in molecular dynamics.
minima [this can be described in terms of Fermi's golden rule for the transition rate: , whereSimilar observations can be made for other structures, such as PbZrO3 and Zr-rich PbZr1−xTixO3, which are known to adopt a cubic structure above the Curie temperature of TC = 523 K (Zhang et al., 2015) [at room temperature, a centrosymmetric structure (space group Pbam) is observed, resulting from antiparallel displacements of the cations on the (110) planes and oxygen octahedral tilts of type a−a−c0 (Glazer et al., 1993)]. At a high temperature of T > 523 K, however, the diffraction pattern contains a considerable amount of diffuse scattering, which can be attributed to distortion modes at the M point in the i.e. correlated dynamical in-phase tiltings along the crystal main axis in agreement with molecular dynamic simulations (Zhang et al., 2015).
3.4. Possible space groups for perovskites with dynamical tilting
Often, the structure ), Aleksandrov (1976) or Howard & Stokes (1998), and therefore can be explained by static tilts.
for perovskites is based on the space groups that are listed in the tables of Glazer (1972As mentioned already, hybrid halide perovskite structural
should not be restricted to these structures, due to the presence of organic ions, which alter crystal symmetry compared with simple perovskites containing monoatomic ions. This has been widely accepted by the research community working with hybrid halide perovskite and helps to avoid wrong crystal symmetry assignment and misinterpretation of phase diagrams.It is described in the Methods that dynamical tilting can lead to space groups such as Fmmm, P4/mmm, Cmmm and Pmmm, which cannot be explained by static tilts. In the following sections (Sections 4.1–4.5), we review published perovskite structures with these space groups and show also that in perovskite CaTiO3 and cryolite Na3AlF6 the known data point towards the presence of dynamical tiltings in these systems.
4. Discussion
Dynamic tilting can explain some observed structural phase transitions and important physical properties within the perovskite class of structures. In many compounds, the composition forbids an atomic arrangement corresponding precisely to the perovskite ABX3 described above. This is often due to the chemical composition involving further chemical species or the ordering of the ions, frequently occurring in a rock-salt-like arrangement of the A and B cations or of the involved octahedra. Many of the resulting compounds are captured by the structure formula A2B′B′′X6 and are called double perovskites or layered perovskites. This class retains the stability and often also the rigidity of the octahedra (Hossain et al., 2018), while more chemical compositions are feasible than in the simple ABX3 composition. Layered perovskites are of great importance due to their strong and unusual magnetic interactions (Bristowe et al., 2015), superconductivity (Bednorz & Müller, 1986), and technical applications (Fan et al., 2015; Granados del Águila et al., 2020). We will therefore include some of these double perovskites in the discussion.
4.1. P4/mmm
Many structures are reported in the P4/mmm BaTiO3 (Buttner & Maslen, 1992) – the name of which is used for the whole of perovskites with P4/mmm symmetry – KCuF3 and KCrF3 (Edwards & Peacock, 1959), CeAlO3 (Tanaka et al., 1993), TlCuF3 (Rüdorff et al., 1963), SrFeO3 (Diodati et al., 2012), CsAuCl3 (Matsushita et al., 2007), and CeGaO3 (Shishido et al., 1997).
For BaTiO3 phase transitions from rhombohedral to orthorhombic, the tetragonal and cubic phases are known (Hayward & Salje, 2002). All phases except the rhombohedral phase show instabilities along several phonon modes, usually displayed as `negative frequencies', see e.g. Fig. 2 of Zhang et al. (2016) or Lebedev (2009). This indicates that the P4/mmm structure shows a negative curvature of the PES at zero tilting angle, e.g. for the longitudinal optic (LO) 1E mode, and the structure therefore might be stabilized by rather than the potential energy.
We would like to stress here that dynamic instabilities do not necessarily indicate a structural instability in the limit T → 0 K. In the framework of quantum mechanics, the ground state of the system at T → 0 K corresponds to the minimum of the total energy, which is the sum of V0 and of the nuclei (referred to as zero-point energy EZP) (see Appendix A1 for a detailed definition). Calculations (Adams & Passerone, 2016) suggest that EZP is particularly small in structures with delocalized atomic configurations represented by a multi-well PES. Therefore, sometimes the total energy (V0 + EZP) can have a minimum at a structure exhibiting a multi-well PES, rather than at a structure corresponding to the minimum of V0 (see also Fig. 4 and Appendix A2).
This important difference between the PES and the total energy could be overseen by atomistic simulations, where often only V0 is accounted for. Therefore, for the remaining structures reported (CeAlO3, TlCuF3, SrFeO3, CsAuCl3 and CeGaO3), the stabilization through the zero-point energy is plausible, as they remain stable in the P4/mmm phase down to low temperatures.
4.2. Cmmm
The Cmmm has been reported for (Li,La)TiO3-perovskite systems, which in turn give their name to the (Sanz et al., 2004). NaIO3 is also known in the Cmmm (Zachariasen, 1928), as well as Nd0.7TiO3 (Sefat et al., 2005).
4.3. Pmmm
Another Pmmm PbTiO3-perovskite (Cole & Espenschied, 1937). At ambient pressure it undergoes a to the cubic phase between 600 and 800 K (Zhu et al., 2011). The same is also found for NaNbO3 (Solovev et al., 1961), Mg0.5W0.5O3 (Zaslavskii & Bryzhina, 1963) and GdCoO3 (Ruggiero & Ferro, 1954). NaNbO3 particularly attracts our attention. It has been refined recently (Peel et al., 2012) for the high temperature phase at 773.15 K. Best refinements were achieved with space groups Pmmm (χ2 = 1.80) and Pnma (χ2 = 1.85).
type is4.4. Perovskite CaTiO3
In CaTiO3, a cascade of phase transitions from Pbnm to Cmcm (1380 K), I4/mcm (1500 K) and finally to the cubic Pmm phase (1580 K) is observed. For the first transition (Pbnm to Cmcm at 1380 K), small anomalies in the temperature dependence of the cell and structural parameters are observed (Kennedy et al., 1999).
However, other authors mention the large atomic displacement parameter of the cubic phase and state: `the high-temperature ). Quenching allows one to access the dynamic disorder of the cubic phase (Britvin et al., 2022). Furthermore, within the approximation of static tilts, the fading of tilting angles at the I4/mcm–Pmm phase boundary at 1500 K – Fig. 4 of Kennedy et al. (1999) – should lead to a significant increase of the volume in a temperature range of 20–50 K, which is not observed. Therefore, the cubic phase should be assigned to a dynamic tilting of oxygen octahedra (adadad) rather to an ideal perovskite structure (a0a0a0).
to cubic perovskite is triggered by the sudden increase of the mobility of the oxygen or at least of parts of it.' (Vogt & Schmahl, 19934.5. Cryolite Na3AlF6
The ; Spearing et al., 1994). The for the transition between the P21/n and the Immm (Anthony et al., 2005) is TC = 885 K (612°C). Due to excellent X-ray data, the atomic positions and the main axes of the vibrational ellipsoids are known (Yang et al., 1993) below and above TC. At low temperature, the system shows static tilts. The P21/n is generated by rigid octahedra and a tilt system with one in-phase tilting and two out-of-phase tiltings – see also Fig. 5. Due to the rock-salt ordering of the octahedra, the classification for simple perovskites (Lufaso & Woodward, 2001; Woodward, 1997a,b) – where a+b−c− corresponds to P21/m – cannot be applied. The experimental high-temperature orthorhombic Immm structure shows displacements of fluorine from the ideal cubic perovskite positions. By further inspection, they result from different bonding lengths in AlO6 and NaO6 octahedra. The pseudo-cubic lattice parameters, on the other hand, suggest a tilting by at least 2.5°
is well documented in cryolite (Steward & Rooksby, 1953Finally, the main axes of vibrational ellipsoids can increase by more than 200% over the
All this points toward a system with at least one dynamic tilting.4.6. Impact of dynamical tilting on the physical properties
The multi-well energy landscape and the amplification of thermal movement upon temperature increase lead to: (a) unusually large thermal anisotropic displacement factors for B and X sites, see e.g. Tyson et al. (2017) for CH3NH3PbI3; (b) negative (directional) temperature expansion, see e.g. Jacobsson et al. (2015) for the negative directional temperature expansion in CH3NH3PbI3; and (c) cell-volume decrease at the upon temperature increase, see e.g. Sakakura et al. (2011) discussing the volume contraction in Na0.5K0.5NbO3 at the phase transitions at 446 and 666 K.
Wright et al. (2016) investigated the electron–phonon coupling in hybrid lead halide perovskites in order to access the coupling transport properties, charge-carrier recombination and finally the charge-carrier mobility. It is surprising that the structural instabilities are overseen in these calculations. The harmonic approximation results in imaginary frequencies in these materials [see also Yang et al. (2017)], which is a signature of the negative curvature of the PES. However, the frequencies cannot be interpreted physically. Adams & Passerone (2016) have shown that dynamic instabilities can lead to high DOS in the vibrational dispersion relation close to ω = 0. Furthermore, the electron–phonon coupling coefficient depends on the atomic mean square deviation (MSD) (Antonius et al., 2015), which diverges for unstable modes in harmonic theory, and which is overestimated for stable modes in harmonic theory (Adams et al., 2020). It thus remains to be seen how the re-evaluation of the electron–phonon matrix element will change our picture of the electron–phonon in hybrid lead halide perovskites.
However, the importance of anharmonic vibrational excitations in the stabilization of the different phases of halide perovskites is known. They are accessed through Monte Carlo simulations [see e.g. Bechtel et al. (2019) for CsPbBr3] or Landau theory, which are computationally costly or semi-empirical, respectively [e.g. for CsPbI3 (Marronnier et al., 2017)]. However, the evaluation of the correct vibrational spectrum, e.g. in decoupled anharmonic mode approximation (DAMA) (Adams & Oganov, 2006; Adams et al., 2020), gives access to the free energy and thus to most physical and thermodynamic properties of the material. In some structures, the formation of dynamic instabilities seems to be fostered by pressure, e.g. in CsAuCl3 (Matsushita et al., 2007). In these structures, the evaluation of the free energy could clarify the role of pressure in the structure stabilization.
As mentioned in the Introduction, the multi-well character of the PES results in a high sensitivity of the ionic positions on perturbations. This is also reflected in the electron–phonon coupling constant
where ℏ is Planck's reduced constant, m is the effective mass of the charge carrier and ω is the frequency of LO phonons (Feynman, 1955). For degenerate eigenstates, the transition frequency vanishes and the coupling therefore diverges. Hence, the concept of polarons has to be extended.
5. Conclusions
5.1. Phase stabilization
Dynamic tilts can emerge upon temperature increase. Their mechanism can be explained based on the DAMA model presented by Adams & Passerone (2016). The double-well PES is the necessary condition for the onset of dynamic tilt (Fig. 2). The energies of the vibrational modes in the double-well potential are more dense than in a single well [Fig. 2(c), compare also the scale of the energy axis]. Most importantly, the ground state of the double well is almost degenerate, leading to a particularly high DOS at the lowest energy [this high density is calculated by Adams & Passerone (2016) for the vibrational DOS of cryolite]. Typically, the high temperature phase lies energetically higher than the low temperature phase, e.g. ΔV = 90 meV per CH3NH3PbCl3 unit has been reported for its cubic phase compared with its orthorhombic phase (Brivio et al., 2015). This energy difference is outweighed by the relevant thermodynamic potential at temperature, which is the free energy
with kB as the T as the temperature and εi as the energy of the vibrational states. A decreases whenever the vibrational energies εi decrease, e.g. through lattice expansion, or here when a leads to a high density of vibrational states at low energies in the spectrum, e.g. through degenerate low-energy phonon modes in the cubic phase of CsPbI3 (Marronnier et al., 2017).
The activation energy, i.e. the energy required to surmount the energy barrier between different energy minima, is generally comparable with the energy of thermal vibrations at the transition temperature [see Fig. 2(a)]. It is particularly low for structures with small tilting, which, according to Glazer (1972), appears mostly in structures where the size of the A cation matches that of the cavity. This leads to a small tilting and thus a low energy barrier with relatively low transition temperatures [e.g. 7.3 meV in CsPbI3 in the cubic phase, which transforms to the tetragonal β-phase at 533.15 K (Marronnier et al., 2017)]. If the A cation is small (e.g. MgSiO3), the tilting is larger, and the can potentially occur at higher temperatures. This type of transition is expected to be very common. In some systems, however, the lies above the melting point and therefore such transition is not observed.
5.2. Group–subgroup relations
In general, the appearance or disappearance of a dynamic tilting mode does not lead to the simple group/subgroup relation between resulting structures. This can be shown in the example of the
transition, where one b− transforms to the ad dynamic tilt. Considering group–subgroup relations, at least one intermediate is involved (either Pbcm or Pmmm) [data from Bilbao Crystallographic Server (Aroyo et al., 2006)]. In other structures, the onset of dynamic tilt does not lead to change, e.g. in the case of
The shortening of the crystal axis related to the onset of the dynamic tilt does not break any symmetry. As shown in these two examples, dynamic tilts do not correspond to irreducible representations of the i.e. the two end members relate in a simple or more complex way.
5.3. Crystallographic consequences of dynamical tilting
As shown in the few examples, the particular relationships in the structures at
can be indicative for transformations driven by dynamic tilts. Dynamic tilt should be considered if:(i) the observed instantaneous symmetry is in contradiction with the average symmetry;
(ii) the experimentally observed ), Woodward (1997a) or Aleksandrov (1976);
of a perovskite-type structure is not listed in the tables established by Glazer (1972(iii) the proportions of the lattice parameters are in contradiction with those suggested by the theory of static tilts, e.g. in the I4/mcm symmetry, which appears in the adadb− system as well as the a0a0b− tilting system, the first one corresponds strictly to a lattice ratio of c/a > 1 while the second one has no restrictions for this ratio;
(iv) the octahedra appear distorted. The octahedra represent a stable configuration of the BX6 chemical configuration. (Large) distortions are unlikely, except in the case of Jahn–Teller distortions.
The success of the theory of static tilts in numerous other systems (Aleksandrov, 1976; Aleksandrov & Bartolome, 1994; Glazer, 1972; Woodward, 1997a,b) underlines the correctness of these considerations.
5.4. Dynamical tilting and electron–phonon interaction
The interaction between phonons and charge carriers is important for the physical mechanisms controlling the mobility of the charge carriers. High mobility is desirable for applications that are based on the efficient separation of electrons and holes, e.g. in photovoltaics. It is significantly influenced by the interaction between photons and the charge carriers. Steele et al. (2019) confirm a strong electron–phonon coupling in CH3NH3PbCl3 and relate it to `strong anharmonicity and dynamic disorder'. Similar results are found for CsPbBr3 and CH3NH3PbBr3 (Sendner et al., 2016). In the work of Batignani et al. (2018), the interdependence between dynamic structural distortions, photo-carriers and photons in lead halide perovskites is also documented. The effect of this coupling on the mobility needs further investigation because at large coupling its temperature dependence can exhibit multiple extrema (Prodanović & Vukmirović, 2019). Such a coupling can have further unexpected effects on the such as the up-conversion, i.e. an increase of the luminescence frequency through from phonons (Granados del Águila et al., 2020).
Based on available observations, we can conclude that dynamical tilting can affect the following properties in perovskites, especially in lead halide perovskites:
(1) The electron–phonon coupling. The electron–phonon coupling coefficient depends on the atomic MSD, which is significantly increased in structures showing dynamical tilting compared with structures without structural instabilities.
(2) The charge-carrier mobility. Amongst others it is limited by the interaction of charge carriers with crystal vibrations (Herz, 2017). This interaction can be described in terms of polarons, which consists of the polarization of the ionic lattice by a mobile electron (Fröhlich, 1954; Feynman, 1955). The resulting field finally creates e.g. the Fröhlich interaction. As mentioned above, in a structure with dynamic tiltings the ionic eigenstates of a double well of the PES are degenerate, which results in high sensitivity of the ionic positions on perturbations, and presumably large electron–phonon coupling constant and polarizability.
(3) The spectral width of light-emitting semiconductor devices (Iaru et al., 2017) by electron–phonon interactions. Wright et al. (2016) underline that the `Fröhlich coupling to LO phonons is the predominant charge-carrier scattering mechanism in hybrid lead halide perovskites', leading to emission linewidth broadening – see also Herz (2017).
All these properties seem to be tightly related to the observed instability on the surface of the 3 (Kirschner et al., 2019) points to stabilization through the free energy [see also equation (3)], and thus supports our model.
and the resulting dynamical tilting. Therefore, the inherent structural sensibility might be the premise for the high efficiency of perovskite solar cells, which indicates that the challenge in the application of these materials is the stabilization of the structure against phase transitions. Temperature stabilization of CSPbIMany perovskite structures show interesting properties such as photoelectricity, but also magnetism, ferroelectricity or superconductivity. The dynamical tiltings result in almost degenerate vibrational modes (Adams & Passerone, 2016), which are worth further investigation. These modes could couple to the electronic states and thus modify the electron–electron interaction in the solid in an unexpected manner. For photoelectric materials, Marronnier et al. (2017) explicitly state that `the perovskite oscillations through the corresponding energy barrier could explain the underlying ferroelectricity and the dynamical Rashba effect predicted in halide perovskites for photovoltaics'. A number of physical properties of perovskite structures are yet to be explained. The space groups for dynamic tilting of cubic perovskites reported in this work will facilitate consideration of dynamical tilting in connection with physical properties of perovskites.
APPENDIX A
Structure stability
A1. The surface
In this article, we define the V0) as the sum of the of the electrons (Te), the due to the interaction of electrons (Uee), the Coulombic energy of electron–nucleus attraction (Vne) and the energy of interactions between nuclei (Vnn), i.e.
of an atomic configuration (For a given atomic configuration, this quantity is readily obtained by widely used ab initio codes such as VASP (Kresse & Furthmüller, 1996) or Quantum ESPRESSO (Giannozzi et al., 2009).
Often, the cohesive energy is estimated from the difference of V0 between different atomic configurations. These estimations correspond to T = 0 K and they most often neglect the of the nuclei, which at T = 0 K is called the zero-point energy EZP.
Minimization of V0 with respect to atomic arrangement and lattice parameters in the structure is used for prediction of structure stability in the limit of T → 0 K, assuming that the atoms behave as classical particles. However, in the framework of quantum mechanics, the stable atomic arrangement at T = 0 K corresponds to the minimum of the total energy, which is defined as the sum of V0 and the of the nuclei, referred to as the zero-point energy EZP (see Appendix A2).
A2. Quasi-harmonic approximation
The quasi-harmonic approximation is a widely used reference model for the description of thermodynamic and elastic properties of solids (Palumbo & Dal Corso, 2017; Togo & Tanaka, 2015). In this model, low-energy vibrational excitations are described through polarization vectors, which account for the contribution of individual atoms to a collective movement at a fixed lattice constant. In the harmonic approximation, the polarization vectors are chosen in such a way that the second-order interactions between modes vanish. Specifically, in perovskite structures, some polarization vectors correspond to rigid-body-like rotation of octahedra. The vibrational spectra of the system depend on the curvature of the PES along these polarization vectors.
The vibration spectra
where i is an integer and
are determined along the ) – see also Fig. 2(c). If , the frequency ω becomes imaginary and the harmonic approximation breaks down.
(Griffiths, 1995At T → 0 K, the stability of the system consisting of nuclei and electrons is defined by the minimum of the total energy of the system, which is different from the minimum of the V0. The difference between these entities is the zero-point energy, EZP = ℏω/2 with ω from equation (8). The value of EZP is thus determined by the shape of the PES in close vicinity to the equilibrium arrangement of atoms in the structure. In the harmonic approximation, EZP is related to the curvature of the PES. Independently of the approximation, the following rule holds: the wider the potential, the smaller the zero-point energy and the closer the total energy to the minimum of the EZP is particularly small in the multi-well PES. This can lead to the stabilization of delocalized configurations with partitioning of the atomic density between multiple minima of the PES (see also Fig. 4), even if the minimum of the V0 of the single-well potential is lower than the minimum of the multi-well potential.
A3. Statistical physics
The quasi-harmonic approximation is widely used to calculate the free energy of solids based on the energy spectrum of the phonons. For the sake of simplicity, the phonon dispersion (k-point sampling) can be accounted for by considering sufficiently large supercells. The energies of the spectra [equation (7)] can be used to obtain the
which in turn allows us to calculate the free energy in equation (3). Using the spectrum [equation (7)] and the Z, it is possible to calculate most of the thermal properties of a material. As temperature increases, the higher vibrational energy levels influence the material properties. The free energy is the thermodynamic potential, which is minimized by the atomic configuration at T > 0 K. Structures showing spectra with low phonon frequency energies minimize εi and thus also A. One implication of this dependence is where larger lattice parameters lead to weaker interactions of atoms and thus also affect the energy spectrum of atomic vibrations (red shift). Therefore, the increasing cell parameters minimize A when T increases (Guyot et al., 1996). On the other hand, it has been shown using the example of cryolite that anharmonic effects can lead to a quasi-degeneracy of numerous low-energy states, resulting in a decrease of free energy and structural stabilization. It is clear that this effect plays an important role in stabilization in various perovskites that show a negative curvature of the potential surface near the equilibrium structure.
Acknowledgements
The authors thank Rene Schliemann for help with geometry determination of the structure and Georgia Cametti for visualization. We are particularly grateful to Professor Dr Frank Nüesch for stimulating discussions and Professor Thomas Armbruster for careful reading of the manuscript.
Funding information
This research was supported by a grant from the Swiss National Supercomputing Centre (CSCS) under project ID s792.
References
Adams, D. J. & Oganov, A. R. (2006). Phys. Rev. B, 73, 184106. Web of Science CrossRef Google Scholar
Adams, D. J. & Passerone, D. (2016). J. Phys. Condens. Matter, 28, 305401. Google Scholar
Adams, D. J., Wang, L., Steinle-Neumann, G., Passerone, D. & Churakov, S. (2020). J. Phys. Condens. Matter 2, 666. Google Scholar
Akbarzadeh, A., Kornev, I., Malibert, C., Bellaiche, L. & Kiat, J.-M. (2005). Phys. Rev. B, 72, 205104. Web of Science CrossRef Google Scholar
Aleksandrov, K. (1976). Ferroelectrics, 14, 801–805. CrossRef CAS Web of Science Google Scholar
Aleksandrov, K. S. & Bartolome, J. (1994). J. Phys. Condens. Matter, 6, 8219–8235. CrossRef CAS Web of Science Google Scholar
Anthony, J. W., Bideaux, R. A., Bladh, K. W. & Nichols, M. C. (2005). Handbook of Mineralogy. Chantilly: Mineralogical Society of America. Google Scholar
Antonius, G., Poncé, S., Lantagne-Hurtubise, E., Auclair, G., Gonze, X. & Côté, M. (2015). Phys. Rev. B, 92, 085137. Web of Science CrossRef Google Scholar
Aroyo, M. I., Perez-Mato, J. M., Capillas, C., Kroumova, E., Ivantchev, S., Madariaga, G., Kirov, A. & Wondratschek, H. (2006). Z. Kristallogr. 221, 15–27. Web of Science CrossRef CAS Google Scholar
Batignani, G., Fumero, G., Kandada, A. R. S., Cerullo, G., Gandini, M., Ferrante, C., Petrozza, A. & Scopigno, T. (2018). Nat. Commun. 9, 1971. Web of Science CrossRef PubMed Google Scholar
Bechtel, J. S., Thomas, J. C. & Van der Ven, A. (2019). Phys. Rev. Mater. 3, 113605. Web of Science CrossRef Google Scholar
Bednorz, J. G. & Müller, K. A. (1986). Z. Phys. B Condens. Matter, 64, 189–193. CrossRef CAS Web of Science Google Scholar
Bednorz, J. G. & Müller, K. A. (1988). Rev. Mod. Phys. 60, 585–600. CrossRef CAS Web of Science Google Scholar
Beecher, A. N., Semonin, O. E., Skelton, J. M., Frost, J. M., Terban, M. W., Zhai, H., Alatas, A., Owen, J. S., Walsh, A. & Billinge, S. J. (2016). ACS Energy Lett. 1, 880–887. Web of Science CrossRef CAS Google Scholar
Bristowe, N., Varignon, J., Fontaine, D., Bousquet, E. & Ghosez, P. (2015). Nat. Commun. 6, 6677. Web of Science CrossRef PubMed Google Scholar
Britvin, S. N., Vlasenko, N. S., Aslandukov, A., Aslandukova, A., Dubrovinsky, L., Gorelova, L. A., Krzhizhanovskaya, M. G., Vereshchagin, O. S., Bocharov, V. N., Shelukhina, Y. S., Lozhkin, M. S., Zaitsev, A. N. & Nestola, F. (2022). Am. Mineral. 107, 1936–1945. Web of Science CrossRef ICSD Google Scholar
Brivio, F., Frost, J. M., Skelton, J. M., Jackson, A. J., Weber, O. J., Weller, M. T., Goñi, A. R., Leguy, A. M., Barnes, P. R. & Walsh, A. (2015). Phys. Rev. B, 92, 144308. Web of Science CrossRef Google Scholar
Buttner, R. H. & Maslen, E. N. (1992). Acta Cryst. B48, 764–769. CrossRef ICSD CAS Web of Science IUCr Journals Google Scholar
Chi, L., Swainson, I., Cranswick, L., Her, J.-H., Stephens, P. & Knop, O. (2005). J. Solid State Chem. 178, 1376–1385. Web of Science CSD CrossRef CAS Google Scholar
Choudhury, N., Walter, E. J., Kolesnikov, A. I. & Loong, C.-K. (2008). Phys. Rev. B, 77, 134111. Web of Science CrossRef Google Scholar
Churakov, S. & Wunder, B. (2004). Phys. Chem. Miner. 31, 131–141. Web of Science CrossRef CAS Google Scholar
Cole, S. S. & Espenschied, H. (1937). J. Phys. Chem. 41, 445–451. CrossRef ICSD CAS Google Scholar
Demetriou, D., Catlow, C., Chadwick, A., Mclntyre, G. & Abrahams, I. (2005). Solid State Ionics, 176, 1571–1575. Web of Science CrossRef ICSD CAS Google Scholar
Diodati, S., Nodari, L., Natile, M., Russo, U., Tondello, E., Lutterotti, L. & Gross, S. (2012). Dalton Trans. 41, 5517–5525. Web of Science CrossRef ICSD CAS PubMed Google Scholar
Edwards, A. J. & Peacock, R. D. (1959). J. Chem. Soc. pp. 4126–4127. Google Scholar
Egger, D. A., Bera, A., Cahen, D., Hodes, G., Kirchartz, T., Kronik, L., Lovrincic, R., Rappe, A. M., Reichman, D. R. & Yaffe, O. (2018). Adv. Mater. 30, 1800691. Web of Science CrossRef Google Scholar
Fan, Z., Xiao, J., Sun, K., Chen, L., Hu, Y., Ouyang, J., Ong, K. P., Zeng, K. & Wang, J. (2015). J. Phys. Chem. Lett. 6, 1155–1161. Web of Science CrossRef CAS PubMed Google Scholar
Feynman, R. P. (1955). Phys. Rev. 97, 660–665. CrossRef CAS Web of Science Google Scholar
Fröhlich, H. (1954). Adv. Phys. 3, 325–361. Google Scholar
Fujii, Y., Hoshino, S., Yamada, Y. & Shirane, G. (1974). Phys. Rev. B, 9, 4549–4559. CrossRef CAS Web of Science Google Scholar
Gao, L., Yadgarov, L., Sharma, R., Korobko, R., McCall, K. M., Fabini, D. H., Stoumpos, C. C., Kanatzidis, M. G., Rappe, A. M. & Yaffe, O. (2021). Mater. Adv. 2, 4610–4616. Web of Science CrossRef CAS Google Scholar
Giannozzi, P., Baroni, S., Bonini, N., Calandra, M., Car, R., Cavazzoni, C., Ceresoli, D., Chiarotti, G. L., Cococcioni, M., Dabo, I., Dal Corso, A., de Gironcoli, S., Fabris, S., Fratesi, G., Gebauer, R., Gerstmann, U., Gougoussis, C., Kokalj, A., Lazzeri, M., Martin-Samos, L., Marzari, N., Mauri, F., Mazzarello, R., Paolini, S., Pasquarello, A., Paulatto, L., Sbraccia, C., Scandolo, S., Sclauzero, G., Seitsonen, A. P., Smogunov, A., Umari, P. & Wentzcovitch, R. M. (2009). J. Phys. Condens. Matter, 21, 395502. Web of Science CrossRef PubMed Google Scholar
Glazer, A. M. (1972). Acta Cryst. B28, 3384–3392. CrossRef CAS IUCr Journals Web of Science Google Scholar
Glazer, A. M., Roleder, K. & Dec, J. (1993). Acta Cryst. B49, 846–852. CrossRef ICSD CAS Web of Science IUCr Journals Google Scholar
Granados del Águila, A., Do, T. T. H., Xing, J., Jee, W. J., Khurgin, J. B. & Xiong, Q. (2020). Nano Res. 13, 1962–1969. Google Scholar
Griffiths, D. J. (1995). Introduction to Quantum Mechanics. Upper Saddle River: Prentice-Hall. Google Scholar
Guyot, F., Wang, Y., Gillet, P. & Ricard, Y. (1996). Phys. Earth Planet. Inter. 98, 17–29. CrossRef CAS Web of Science Google Scholar
Hayward, S. & Salje, E. (2002). J. Phys. Condens. Matter, 14, L599–L604. Web of Science CrossRef CAS Google Scholar
Herz, L. M. (2017). ACS Energy Lett. 2, 1539–1548. Web of Science CrossRef CAS Google Scholar
Heyraud, S., Blanchard, P. E., Liu, S., Zhou, Q., Kennedy, B. J., Brand, H. E., Tadich, A. & Hester, J. R. (2013). J. Phys. Condens. Matter, 25, 335401. Web of Science CrossRef ICSD PubMed Google Scholar
Hossain, A., Bandyopadhyay, P. & Roy, S. (2018). J. Alloys Compd. 740, 414–427. Web of Science CrossRef CAS Google Scholar
Howard, C. J. & Stokes, H. T. (1998). Acta Cryst. B54, 782–789. Web of Science CrossRef CAS IUCr Journals Google Scholar
Iaru, C. M., Geuchies, J. J., Koenraad, P. M., Vanmaekelbergh, D. & Silov, A. Y. (2017). ACS Nano, 11, 11024–11030. Web of Science CrossRef CAS PubMed Google Scholar
Jacobsson, T. J., Schwan, L. J., Ottosson, M., Hagfeldt, A. & Edvinsson, T. (2015). Inorg. Chem. 54, 10678–10685. Web of Science CrossRef CAS PubMed Google Scholar
Kassan-Ogly, F. A. & Naish, V. E. (1986). Acta Cryst. B42, 325–335. CrossRef Web of Science IUCr Journals Google Scholar
Kennedy, B. J., Prodjosantoso, A. & Howard, C. J. (1999). J. Phys. Condens. Matter, 11, 6319–6327. Web of Science CrossRef ICSD CAS Google Scholar
Kirschner, M. S., Diroll, B. T., Guo, P., Harvey, S. M., Helweh, W., Flanders, N. C., Brumberg, A., Watkins, N. E., Leonard, A. A., Evans, A. M., Wasielewski, M. R., Dichtel, W. R., Zhang, X., Chen, L. X. & Schaller, R. D. (2019). Nat. Commun. 10, 504. Google Scholar
Klarbring, J. (2019). Phys. Rev. B, 99, 104105. Web of Science CrossRef Google Scholar
Knížek, K., Hejtmánek, J., Jirák, Z., Martin, C., Hervieu, M., Raveau, B., André, G. & Bourée, F. (2004). Chem. Mater. 16, 1104–1110. Google Scholar
Kresse, G. & Furthmüller, J. (1996). Phys. Rev. B, 54, 11169–11186. CrossRef CAS Web of Science Google Scholar
Lebedev, A. (2009). Phys. Solid State, 51, 2324–2333. Web of Science CrossRef CAS Google Scholar
Lehmann, F., Franz, A., Többens, D. M., Levcenco, S., Unold, T., Taubert, A. & Schorr, S. (2019). RSC Adv. 9, 11151–11159. Web of Science CrossRef CAS PubMed Google Scholar
Lin, Z., Zhang, Y., Gao, M., Steele, J. A., Louisia, S., Yu, S., Quan, L. N., Lin, C.-K., Limmer, D. T. & Yang, P. (2021). Matter, 4, 2392–2402. Web of Science CrossRef CAS Google Scholar
Lufaso, M. W. & Woodward, P. M. (2001). Acta Cryst. B57, 725–738. Web of Science CrossRef CAS IUCr Journals Google Scholar
Mante, P.-A., Stoumpos, C. C., Kanatzidis, M. G. & Yartsev, A. (2018). J. Phys. Chem. Lett. 9, 3161–3166. Web of Science CrossRef CAS PubMed Google Scholar
Marronnier, A., Lee, H., Geffroy, B., Even, J., Bonnassieux, Y. & Roma, G. (2017). J. Phys. Chem. Lett. 8, 2659–2665. Web of Science CrossRef CAS PubMed Google Scholar
Matsushita, N., Ahsbahs, H., Hafner, S. S. & Kojima, N. (2007). J. Solid State Chem. 180, 1353–1364. Web of Science CrossRef ICSD CAS Google Scholar
Munson, K. T., Swartzfager, J. R. & Asbury, J. B. (2019). ACS Energy Lett. 4, 1888–1897. Web of Science CrossRef CAS Google Scholar
Onoda-Yamamuro, N., Yamamuro, O., Matsuo, T. & Suga, H. (1992). J. Phys. Chem. Solids, 53, 277–281. CAS Google Scholar
Palumbo, M. & Dal Corso, A. (2017). J. Phys. Condens. Matter, 29, 395401. Web of Science CrossRef PubMed Google Scholar
Peel, M. D., Thompson, S. P., Daoud-Aladine, A., Ashbrook, S. E. & Lightfoot, P. (2012). Inorg. Chem. 51, 6876–6889. Web of Science CrossRef ICSD CAS PubMed Google Scholar
Poglitsch, A. & Weber, D. (1987). J. Chem. Phys. 87, 6373–6378. CrossRef CAS Web of Science Google Scholar
Prodanović, N. & Vukmirović, N. (2019). Phys. Rev. B, 99, 104304. Google Scholar
Rüdorff, W., Lincke, G. & Babel, D. (1963). Z. Anorg. Allg. Chem. 320, 150–170. Google Scholar
Ruggiero, A. & Ferro, R. (1954). Atti della Accademia Nazionale dei Lincei, Classe di Scienze Fisiche, Matematiche e Naturali, 8, 254–256. Google Scholar
Sakakura, T., Wang, J., Ishizawa, N., Inagaki, Y. & Kakimoto, K. (2011). IOP Conference Series: Materials Science and Engineering, Vol. 18, p. 022006. Bristol: IOP Publishing. Google Scholar
Sanz, J., Varez, A., Alonso, J. A. & Fernandez, M. T. (2004). J. Solid State Chem. 177, 1157–1164. Web of Science CrossRef ICSD CAS Google Scholar
Sefat, A. S., Amow, G., Wu, M.-Y., Botton, G. A. & Greedan, J. (2005). J. Solid State Chem. 178, 1008–1016. Web of Science CrossRef ICSD CAS Google Scholar
Sendner, M., Nayak, P. K., Egger, D. A., Beck, S., Müller, C., Epding, B., Kowalsky, W., Kronik, L., Snaith, H. J., Pucci, A. & Lovrinčić, R. (2016). Mater. Horiz. 3, 613–620. Web of Science CrossRef CAS Google Scholar
Shishido, T., Zheng, Y., Saito, A., Horiuchi, H., Kudou, K., Okada, S. & Fukuda, T. (1997). J. Alloys Compd. 260, 88–92. CrossRef ICSD CAS Web of Science Google Scholar
Solovev, S., Venevtsev, Y. N. & Zhanov, G. (1961). Sov. Phys. Crystallogr. 6, 171–175. Google Scholar
Spearing, D., Stebbins, J. & Farnan, I. (1994). Phys. Chem. Miner. 21, 373–386. CrossRef CAS Google Scholar
Steele, J. A., Puech, P., Monserrat, B., Wu, B., Yang, R. X., Kirchartz, T., Yuan, H., Fleury, G., Giovanni, D., Fron, E., Keshavarz, M., Debroye, E., Zhou, G., Sum, T. C., Walsh, A., Hofkens, J. & Roeffaers, M. B. J. (2019). ACS Energy Lett. 4, 2205–2212. Web of Science CrossRef CAS Google Scholar
Steward, E. G. & Rooksby, H. P. (1953). Acta Cryst. 6, 49–52. CrossRef IUCr Journals Web of Science Google Scholar
Stokes, H. T. & Hatch, D. M. (2005). J. Appl. Cryst. 38, 237–238. CrossRef CAS IUCr Journals Google Scholar
Swainson, I., Hammond, R., Soullière, C., Knop, O. & Massa, W. (2003). J. Solid State Chem. 176, 97–104. Web of Science CSD CrossRef CAS Google Scholar
Swainson, I., Stock, C., Parker, S., Van Eijck, L., Russina, M. & Taylor, J. (2015). Phys. Rev. B, 92, 100303. Web of Science CrossRef Google Scholar
Tanaka, M., Shishido, T., Horiuchi, H., Toyota, N., Shindo, D. & Fukuda, T. (1993). J. Alloys Compd. 192, 87–89. CrossRef ICSD CAS Web of Science Google Scholar
Togo, A. & Tanaka, I. (2015). Scr. Mater. 108, 1–5. Web of Science CrossRef CAS Google Scholar
Tyson, T., Gao, W., Chen, Y.-S., Ghose, S. & Yan, Y. (2017). Sci. Rep. 7, 9401. Web of Science CrossRef PubMed Google Scholar
Ugur, E., Alarousu, E., Khan, J. I., Vlk, A., Aydin, E., De Bastiani, M., Balawi, A. H., Gonzalez Lopez, S. P., Ledinský, M., De Wolf, S. & Laquai, F. (2020). Sol. RRL, 4, 2000382. Web of Science CrossRef Google Scholar
Vogt, T. & Schmahl, W. W. (1993). Europhys. Lett. 24, 281–285. CrossRef CAS Web of Science Google Scholar
Völkl, J. & Alefeld, G. (1978). Hydrogen in Metals I, pp. 321–348. Berlin: Springer. Google Scholar
Whalley, L. D., Frost, J. M., Jung, Y.-K. & Walsh, A. (2017). J. Chem. Phys. 146, 220901. Web of Science CrossRef PubMed Google Scholar
Woodward, P. M. (1997a). Acta Cryst. B53, 32–43. CrossRef CAS Web of Science IUCr Journals Google Scholar
Woodward, P. M. (1997b). Acta Cryst. B53, 44–66. CrossRef CAS Web of Science IUCr Journals Google Scholar
Wright, A. D., Verdi, C., Milot, R. L., Eperon, G. E., Pérez-Osorio, M. A., Snaith, H. J., Giustino, F., Johnston, M. B. & Herz, L. M. (2016). Nat. Commun. 7, 11755. Google Scholar
Yang, H., Ghose, S. & Hatch, D. (1993). Phys. Chem. Miner. 19, 528–544. CrossRef ICSD CAS Google Scholar
Yang, J., Wen, X., Xia, H., Sheng, R., Ma, Q., Kim, J., Tapping, P., Harada, T., Kee, T. W., Huang, F., Cheng, Y. B., Green, M., Ho-Baillie, A., Huang, S., Shrestha, S., Patterson, R. & Conibeer, G. (2017). Nat. Commun. 8, 14120. Google Scholar
Zachariasen, W. H. (1928). Matematisk Naturvidenskapelig Klasse, pp. 1–165. Google Scholar
Zaslavskii, A. & Bryzhina, M. (1963). Sov. Phys. Crystallogr. 7, 577–583. Google Scholar
Zhang, H.-Y., Zeng, Z.-Y., Zhao, Y.-Q., Lu, Q. & Cheng, Y. (2016). Z. Naturforsch. A, 71, 759–768. Web of Science CrossRef CAS Google Scholar
Zhang, N., Paściak, M., Glazer, A. M., Hlinka, J., Gutmann, M., Sparkes, H. A., Welberry, T. R., Majchrowski, A., Roleder, K., Xie, Y. & Ye, Z.-G. (2015). J. Appl. Cryst. 48, 1637–1644. Web of Science CrossRef CAS IUCr Journals Google Scholar
Zhao, D., Hu, H., Haselsberger, R., Marcus, R. A., Michel-Beyerle, M.-E., Lam, Y. M., Zhu, J.-X., La-O-Vorakiat, C., Beard, M. C. & Chia, E. E. (2019). ACS Nano, 13, 8826–8835. Web of Science CrossRef CAS PubMed Google Scholar
Zhu, J., Xu, H., Zhang, J., Jin, C., Wang, L. & Zhao, Y. (2011). J. Appl. Phys. 110, 084103. Web of Science CrossRef Google Scholar
Zhu, T. & Ertekin, E. (2019). Energy Environ. Sci. 12, 216–229. Web of Science CrossRef CAS Google Scholar
This is an open-access article distributed under the terms of the Creative Commons Attribution (CC-BY) Licence, which permits unrestricted use, distribution, and reproduction in any medium, provided the original authors and source are cited.