Thank you for visiting nature.com. You are using a browser version with limited support for CSS. To obtain the best experience, we recommend you use a more up to date browser (or turn off compatibility mode in Internet Explorer). In the meantime, to ensure continued support, we are displaying the site without styles and JavaScript.

# Prediction of high thermoelectric performance in the low-dimensional metal halide Cs3Cu2I5

## Abstract

Metal halides have emerged as a new generation of semiconductors with applications ranging from solar cells to chemical sensors. We assess the thermoelectric potential of Cs3Cu2I5, which has a crystal structure formed of zero-dimensional [Cu2I5]3− anionic clusters that are separated by Cs+ counter cations. We find the compound exhibits the characteristics of a phonon-glass electron-crystal with a large imbalance in the conduction of heat and electrons predicted from first-principles transport theory. Strong anharmonic phonon–phonon scattering results in short-lived acoustic vibrations and an ultra-low lattice thermal conductivity (<0.1 W m−1 K−1). The dispersive conduction band leads to a high electron mobility (>10 cm2 V−1 s−1). For an n-type crystal at 600 K, a thermoelectric figure-of-merit ZT of 2.6 is found to be accessible, which for a cold source of 300 K corresponds to a thermodynamic heat-to-electricity conversion efficiency of 15%.

## Introduction

The thermoelectric effect enables the direct conversion of a temperature gradient into an electric voltage and vice versa1. Hence, there have been many studies to discover new materials that show high-energy conversion efficiency2,3,4. The performance of thermoelectric materials can conveniently be evaluated by the dimensionless figure of merit ZT, which is defined as

$$ZT=\frac{{S}^{2}\sigma T}{{\kappa }_{{\rm{elec}}}+{\kappa }_{{\rm{latt}}}}$$
(1)

where S is the Seebeck coefficient, σ is the electrical conductivity, T is the temperature, κelec is the electronic thermal conductivity, and κlatt is the lattice thermal conductivity. Achieving high Seebeck coefficient and electrical conductivity with low electronic and lattice thermal conductivity is therefore the formula for success. In practice, there are limits in the accessible property space for semiconductors, including the positive correlation between Seebeck coefficient, electrical conductivity, and electronic thermal conductivity5. Thus, low lattice thermal conductivity is one of promising design strategies for higher performing thermoelectric materials.

We have been inspired by observations of the thermoelectric effect in halide perovskites6,7,8, where a record ZT of 0.14 was recently achieved9. However, there is a strong desire to move away from heavy and toxic elements such as Pb. Cs3Cu2I5 is an alternative metal halide that has drawn attention for optoelectronics. Since the report of its application in blue light emitting diodes (LEDs)10, the optical properties of this material has been intensively investigated11,12,13. Cs3Cu2I5 has an unusual crystal structure where [Cu2I5]3− clusters are surrounded and isolated by Cs+ cations. Unlike perovskites, for example, this Pnma structure has no through-bond connectivity between neighboring anionic fragments. The atomic structure of Cs3Cu2I5 was recently linked to its optical properties including self-trapped excitons and a large Stokes shift in the photoluminescence11,13. Complex crystal structures have been associated with low thermal conductivity in the past14. Indeed, a value of 0.15 W m−1 K−1 was recently measured in [Mn(C2H6OS)6]I4, which shares a low-dimensional crystal structure15.

In this study, we assess the thermoelectric potential of Cs3Cu2I5. We calculate the thermal transport using anharmonic lattice dynamics, which highlights an ultra-low conductivity for the perfect crystal that is below the limit of standard models. We further investigate electron transport through an ab initio description of carrier scattering. The same low-dimensional crystal structure that gives rise to ultra-low thermal transport supports efficient electron transport with a room temperature electron mobility of 18.2 cm2 V−1 s−1. We predict that Cs3Cu2I5 exhibits the characteristics of an ideal thermoelectric semiconductor—a phonon-glass electron-crystal—with an accessible ZT of greater than 2.

## Results

### Basic bulk properties

The ground-state Pnma crystal structure of Cs3Cu2I5 is illustrated in Fig. 1a. Cs atoms reside at 4c and 8d Wyckoff positions, Cu atoms reside at 4c, and I atoms reside at three distinct 4c and 8d. An associated electron localization function (ELF) plot could be found in Supplementary Fig. 1. The atomic structure is similar to low-dimensional perovskite-inspired materials such as Cs4PbBr6 where metal halide anion clusters are surrounded and isolated by cation counter parts (effective 0D connectivity)16,17. However, due to distinct size and orbital configuration of Cu(I) 3d104s0, it does not form octahedral clusters like Pb(II) 5d106s218. The calculated equilibrium lattice constants of a0 = 10.06 Å b0 = 11.53 Å and c0 = 14.15 Å match well with values determined from X-ray diffraction ($${a}_{0}^{\exp }=10.19$$ Å $${b}_{0}^{\exp }=11.66$$ Å and $${c}_{0}^{\exp }=14.40$$Å)10. By fitting the energy-volume relationship to the third-order Birch-Murnaghan equation-of-state19, we obtained a bulk modulus of B0 = 10.6 GPa for Cs3Cu2I5, which is even softer than Cs4PbBr6 (B0 = 12.1 GPa)20. For comparison, the bulk modulus of Si is ~100 GPa.

The electronic band structure is shown in Fig. 1b. This has been calculated using density functional theory (DFT) with the hybrid HSE06 functional using Wannier interpolation21,22,23. The upper valence band, formed of Cu 3d and I 5p orbitals, is found to be flat in reciprocal space (i.e., short-range interactions in real space). At the same time, the lower conduction band is dispersive and formed of Cu 4s and I 5s (see Supplementary Fig. 2). The associated electron effective mass is 0.3 me, which is comparable to other wide band gap semiconductors such as In2O3 and ZnO. The behavior contrasts with other 0D metal halides that have a flat conduction band dispersion (e.g., Cs4PbBr6) where the lower conduction band is formed from Pb 6p)20. The dispersive nature and isotropic overlap of the Cu 4s orbitals contributes to the low electron effective mass in Cs3Cu2I5. For Cs3Cu2I5, the direct band gap of 3.82 eV at Γ point is obtained from HSE06 calculation, which matches well with measured values for the optical absorption onset12.

The phonon dispersion of Cs3Cu2I5 is shown in Fig. 1c. The lattice vibrations were calculated using the finite-displacement method and the PBEsol functional due to its good description of the crystal structure and phonon properties of semiconductors24,25,26. Since the primitive cell contains 40 atoms, there are 120 phonon bands where 3 acoustic branches are followed by a series of localized optic modes up to 7 THz. We find no imaginary frequencies, which confirms dynamic stability for this phase of Cs3Cu2I5. There is a dense block of optic modes that run from 0.5–2.5 THz involving the collection motion of group atoms, such as relative tilting of [Cu2I5]3− clusters. As a result, acoustic branches are strongly suppressed < 0.5 THz. The presence of this dense block of low-energy phonons just above the acoustic branch creates many energy and momentum allowed scattering pathways that significantly reduce the lifetime of the heat carrying acoustic modes. Such a high density of low-frequency modes also explains the small bulk modulus. Inspection of the eigenvectors at Γ point shows that acoustic modes are simple translation of the entire lattice sites, while the highest frequency optic mode involves Cu–Cu stretching within each cluster. Zoomed electronic and phonon band structures are shown in Supplementary Fig. 3.

### Ultra-low lattice thermal conductivity

Next we consider the lattice thermal conductivity (κlatt) by taking into account three-phonon scattering. This approach should provide an upper limit to the thermal transport as additional scattering events can be active in real crystals. Three-phonon scattering is subject to energy and momentum conservation rules for phonon creation and annihilation; however, the flat nature of the phonon dispersion for Cs3Cu2I5 opens up many allowed channels. Owing to structural anisotropy, the conductivity along three principal directions (κxx, κyy, and κzz) differs, as shown in Fig. 2a. The largest values are found along the b axis, which is parallel to the base of the trigonal-planar Cu atoms, and has larger group velocities. The isotropically-averaged conductivity (κiso) is 0.018 W m−1 K−1 at room temperature. We tested both the single mode relaxation time approximation (RTA) and direct solution of linearized phonon Boltzmann equation (LBTE)27, which are in good agreement. The LBTE value is ~10% lower across the relevant temperature range (see Supplementary Fig. 4).

The lattice thermal conductivity is among the lowest reported for bulk crystals. For instance, κiso of CH3NH3PbI328, Sn2S329, and Cs2BiAgBr630 was reported to be 0.05, 0.06, and 0.07 W m−1 K−1, respectively, at similar levels of theory. This values is also lower than the minimum lattice thermal conductivity ($${\kappa }_{\min }$$) expected from diffuson-mediated thermal transport31, which can be estimated from:

$${\kappa }_{\min }=0.76{k}_{{\rm{B}}}\frac{{v}_{{\rm{s}}}}{{V}_{{\rm{a}}}^{2/3}}$$
(2)

where kB is Boltzmann constant, vs is the speed of sound, and Va is the volume of unit cell per atom. The model predicts $${\kappa }_{\min }$$ of 0.1 W m−1 K−1 for Cs3Cu2I5. To understand the origin of such a low-lattice thermal conductivity, below the prediction of the analytical model, we analyzed the modal contributions to the net transport. Extremely short acoustic phonon lifetimes of sub 12 ps (Fig. 2b) result in short phonon mean-free paths of below 24 nm. Thus, we conclude that high-acoustic phonon-scattering rates are the main origin of the ultra-low lattice thermal conductivity in Cs3Cu2I5 despite their moderate group velocity and heat capacity. A full breakdown of the modal contributions is provided in Supplementary Fig. 5. The dense block of low-energy phonons have similarity to the states created by the CH3NH3 cation in CH3NH3PbI332. Although a similar computational approach has been shown to give good agreement with experimental measurements on the thermoelectric propoerties of metal halides such as CsSnI333, it should be noted that perturbation theory based on three-phonon scattering is pushed to its limits of applicability in such cases34,35. We propose that Cs3Cu2I5 is ideal for testing more sophisticated treatments of anharmonicity including self-consistent phonon theory36,37.

### Electronic transport properties

Given the wide electronic band gap of 3.82 eV, pristine Cs3Cu2I5 is an insulator as there is no significant thermal population of electrons and holes. However, high carrier concentrations (~1021 cm−3) can be achieved in such wide band gap semiconductors by doping (e.g. Al or Ga doped ZnO38,39 and Zr or W doped In2O340,41), which is an active research topic in transparent conductors42. Motivated by the small electron effective mass, we performed electron transport simulations under n-type doping conditions, in the range 1016–1020 cm−3. Our preliminary calculations suggest that Ba$${\,}_{{\rm{Cs}}}^{+}$$ has a shallow donor level, which could support an n-type doping density. The associated electrical transport properties are shown in Fig. 3. The electronic thermal conductivity remains lower than the lattice thermal conductivity for ne≤ 1018 cm−3, which means that thermal losses through phonons will be more critical for thermoelectric applications. Once the carrier concentration exceeds 1019 cm−3, the electronic component becomes dominant.

The electron mobility of Cs3Cu2I5 is found to be limited by optical phonon scattering, followed by ionized impurity and acoustic deformation potential scattering (see Supplementary Fig. 6a), which is consistent with the behavior of other polar semiconductors43,44. The electron relaxation times and associated materials properties are given in Table 1. The electron mobility at room temperature is 18.2 cm2 V−1 s−1 for an n-type doping concentration of 4 × 1018 cm−3. The simpler constant-relaxation time approximation (CRT), assuming the standard τ = 10−14 s, predicts an overestimated value of 60.1 cm2 V−1 s−1 (see Supplementary Fig. 6b). The electrical conductivity increases from around 3.5 to 2.3 × 104 Sm−1, while the absolute value of Seebeck coefficient decreases from around 750 to 60 μV K−1 in the given doping range at room temperature. The power factor (S2σ) is maximized around 1020 cm−3 at a value of 150 μW m−1 K−2.

### Thermoelectric performance prediction

By combining the results on thermal and electron transport, the predicted ZT of Cs3Cu2I5 as a function n-type doping concentration and temperature is shown in Fig. 4a. Here, the temperature range is set to 600 K as prediction for temperature region near the melting temperature of 663 K11 will not be reliable due to possibility of evaporation and mass loss, which is beyond limit of our simulation techniques. At 600 K maximum ZT is predicted to be 2.57 where a n-type doping concentration of 4 × 1018 cm−3 is assumed. Above this doping level, electronic thermal conductivity causes ZT to decrease. When the anisotropic nature of transport properties along the principal axes is considered, ZT can be maximized along c with a ZT of 3.35 due to a combination of lower lattice thermal conductivity and higher power factor. A full analysis of the anisotropy is found in Supplementary Fig. 7. The maximum ZT of Cs3Cu2I5 predicted is the highest for any metal halide compound. While the power factor is moderate, ZT is enhanced by the ultra-low lattice thermal conductivity. Although there have been attempts to make thermoelectric devices using metal halide perovksites8, the performance has been limited thus far. A value of ZT = 0.14 was achieved for CsSnI37, with a similar ZT = 0.13 for CH3NH3SnI36. On the other hand Pb-based halide perovskites have shown negligible thermoelectric response6,8. The predicted ZT of Cs3Cu2I5 is comparable to other emerging thermoelectric compounds. For instance, the maximum ZT of CsBi4Te6, Zn4Sb3, and, SnSe are 0.8 at 225 K45, 1.4 at 750 K46, 2.6 at 923 K47, respectively.

The efficiency of a thermoelectric device (η) is not determined by maximum ZT at a single temperature but depends on the average ZT ($$\overline{ZT}$$) over a wide temperature range following

$$\eta =\frac{{T}_{{\rm{H}}}-{T}_{{\rm{C}}}}{{T}_{{\rm{H}}}}\frac{\sqrt{1+\overline{ZT}}-1}{\sqrt{1+\overline{ZT}}+\frac{{T}_{{\rm{C}}}}{{T}_{{\rm{H}}}}}$$
(3)

where TH and TC are the temperature of hot side and cold side of the device, and $$\overline{ZT}=\frac{1}{({T}_{{\rm{H}}}-{T}_{{\rm{C}}})}\mathop{\int}\nolimits_{{T}_{{\rm{C}}}}^{{T}_{{\rm{H}}}}ZTdT$$. To further evaluate the performance of n-type Cs3Cu2I5-based thermoelectric devices, we plot efficiency over temperature gradient where TH is fixed as 600 K in Fig. 4b. The average $$\overline{ZT}$$ of 1.67 leads to a predicted efficiency of 14.9% over a 300–600 K temperature gradient, where ne is 4 × 1018 cm−3.

Beyond the performance of the bulk metal halide, existing crystal engineering approaches could also be applied to Cs3Cu2I5. Cation (e.g., Cs1−xRbx) and/or anion (e.g., I1−xBrx) substitution is one way to tune ZT, which has been widely adopted to modulate optoelectronic properties in metal halide perovskites11,48,49,50. Such ion substitution could further reduce the lattice thermal conductivity due to the additional alloy-scattering channels51,52. In practice, recent reports have shown improved values of ZT by partially replacing I in CsSnI3, and achieved ZT of 0.14 for CsSnCl3−xIx9 and 0.15 for CsSnBr3−xIx53.

## Discussion

We have demonstrated the potential of a metal halide based on a low-dimensional bonding network for heat-to-electricity conversion. Cs3Cu2I5 is predicted to support efficient electron transport with slow heat transport. An ultra-low lattice thermal conductivity of < 0.1 W m−1 K−1 is predicted, among the lowest reported for a crystalline solid. Owing to its effective thermal insulation, the maximum ZT of n-type Cs3Cu2I5 is predicted to be 2.57 at 600 K, which is over 10 times larger than the current record of 0.15 in CsSnBr3−xIx and close to the champion value of 2.6 achieved in SnSe. By alloying or nanostructuring, the ZT may be even further enhanced. Our work therefore provides an alternative avenue for developing high-performing thermoelectric materials. Those materials with low-dimensional crystal structure and dispersive band structure (e.g., band edges with s-s orbital overlap) could be attractive candidates. The main bottleneck to realizing such a high ZT will be achieving appropriate carrier concentrations. Further investigations are required to enhance the electrical conductivity of multinary metal halides beyond their natural limits.

## Methods

### Thermal transport

The thermal physics was simulated from lattice dynamics calculations considering the phonon eigenvectors, eigenvalues, and lifetimes due to three-phonon scattering. First, harmonic phonon analysis was performed using the Phonopy54 code with a finite-displacement step of 0.01 Å. A 2 × 2 × 2 supercell expansion of the relaxed unit cell, containing 320 atoms, was employed. In total 36 symmetry-reduced displacements were considered. Then, anharmonic calculations were performed using the Phono3py55 code with a finite-displacement step of 0.03 Å. For the third-order force constants, 7818 symmetry-reduced displacements were considered in the 40 atom unit cell. Γ-centered q-meshes of 24 × 24 × 18 and 8 × 8 × 6 were used to compute the phonon density of states and lattice thermal conductivity. Symmetrization of the force constants was applied where the maximum drift value of the third order force constants is 0.73 eV Å−1 before symmetrization. All norms of the atomic forces sets and a convergence test on our thermal transport calculations over q-meshes can be found in Supplementary Fig. 8.

The lifetimes are solved through a direct solution of the linearized Boltzmann transport equation (LBTE)27. Then the macroscopic lattice thermal conductivity κlatt is obtained by summing over phonon band indices (νλ) and wavevectors (q), while normalizing by the volume of the unit cell:

$${\kappa }_{{\rm{latt}}}=\frac{1}{N{V}_{0}}\sum _{\lambda }{\kappa }_{\lambda }=\frac{1}{N{V}_{0}}\sum _{\lambda }{C}_{\lambda }{{\bf{v}}}_{\lambda }\otimes {{\bf{v}}}_{\lambda }{\tau }_{\lambda }$$
(4)

where N is the number of unit cells in the crystal (equivalent to the number of wavevectors included in the Brillouin zone summation) and V0 is the volume of the crystallographic unit cell. κλ contains the product of the modal heat capacity (Cλ), group velocity (vλ), and phonon mean-free path (vλ × τλ, where τλ is the phonon lifetime). Graphical analyses of the modal contributions to lattice thermal conductivity and the atomic force sets used for the thermal conductivity calculations were performed using the Phono3py-Power-Tool code56.

### Electron transport

The transport of electrons was simulated by considering three scattering processes, with a characteristic scattering time following Matthiessen’s rule:

$$\frac{1}{{\tau }_{e}}=\frac{1}{{\tau }^{{\rm{ADP}}}}+\frac{1}{{\tau }^{{\rm{IMP}}}}+\frac{1}{{\tau }^{{\rm{POP}}}}$$
(5)

Each component of τe was calculated from Fermi’s golden rule:

$${\tilde{\tau }}_{n{\bf{k}}\to m{\bf{k}}+{\bf{q}}}^{-1}=\frac{2\pi }{\hslash }| {g}_{nm}({\bf{k}},{\bf{q}}){| }^{2}\delta ({\varepsilon }_{n{\bf{k}}}-{\varepsilon }_{m{\bf{k}}+{\bf{q}}})$$
(6)

where is the reduced Planck constant, ε is the electron energy, δ is the Dirac delta function and g is the coupling matrix element. Matrix elements due to the acoustic deformation potential (ADP), ionized impurity (IMP), and polar-optical phonon (POP) mechanisms are considered:

$${g}_{nm}^{{\rm{ADP}}}({\bf{k}},{\bf{q}})={\left[\frac{{k}_{{\rm{B}}}T{\alpha }_{v}^{2}}{{B}_{0}}\right]}^{\frac{1}{2}}\langle {\psi }_{m{\bf{k}}+{\bf{q}}}| {\psi }_{n{\bf{k}}}\rangle$$
(7)
$${g}_{nm}^{{\rm{IMP}}}({\bf{k}},{\bf{q}})={\left[\frac{{n}_{{\rm{ii}}}{Z}^{2}{e}^{2}}{{\epsilon }_{0}}\right]}^{\frac{1}{2}}\frac{\langle {\psi }_{m{\bf{k}}+{\bf{q}}}| {\psi }_{n{\bf{k}}}\rangle }{{\left|{\bf{q}}\right|}^{2}+{\beta }^{2}},$$
(8)
$${g}_{nm}^{{\rm{POP}}}({\bf{k}},{\bf{q}})={\left[\frac{\hslash {\omega }_{{\rm{po}}}}{2}\left(\frac{1}{{\epsilon }_{\infty }}-\frac{1}{{\epsilon }_{0}}\right)\right]}^{\frac{1}{2}}\frac{\langle {\psi }_{m{\bf{k}}+{\bf{q}}}| {\psi }_{n{\bf{k}}}\rangle }{\left|{\bf{q}}\right|}$$
(9)

The formalism is implemented in the ab initio scattering and transport (AMSET) package44. The required material parameters, including dielectric constants ϵ, phonon frequency ω, bulk modulus B0, deformation potential α, were all calculated from first-principles. The associated Seebeck coefficient, electrical conductivity, and electronic component of the thermal conductivity are calculated from the Onsager transport coefficients57. An interpolation factor of 10 was used during AMSET routine.

### Density functional theory

All of the underlying total energy and electronic structure calculations were performed based on Kohn–Sham density functional theory58, where periodic boundary conditions are considered to represent the extended crystal. The Vienna Ab Initio Simulation Package (VASP)59,60 was used with Projector augmented-wave (PAW)61,62 method where the valence states of Cs, Cu, and I are treated explicitly by 9(5s25p66s1), 17(2p63d104s1), and 7(5s25p5) electrons, respectively.

The Perdew–Burke–Ernzerhof exchange-correlation functional revised for solids (PBEsol)63 with the plane-wave kinetic energy cutoff of 700 eV was used to optimize crystal structure including the internal positions. A 3 × 3 × 2 Γ-centered k-mesh was adopted for structure optimization, where the convergence criteria were set to 10−6 eV and 10−3 eV Å−1 for total energy and atomic forces, respectively. This crystal structure was then used consistently for all transport calculations.

For thermal transport, the same DFT setup was employed and the Brillouin zone was sampled with 2 × 2 × 1 k-mesh for the 2 × 2 × 2 supercell. For electron transport, the required Kohn–Sham wavefunctions (ψ) were generated using the hybrid functional HSE0664,65. Here we used a dense 6 × 6 × 4 Γ-centered k-mesh for better sampling of the electronic Brillouin zone, but lowered the kinetic energy cutoff to 400 eV in order to reduce the computational load. A convergence test on our electronic transport calculations over k-meshes and interpolation factors can be found in Supplementary Fig. 9.

## Data availability

An online repository containing the optimized crystal structure, force constant sets, and raw AMSET input/output files have been made available at https://doi.org/10.5281/zenodo.4576211.

## Code availability

Post-processing scripts for AMSET package have been made available at https://doi.org/10.5281/zenodo.4576211.

## References

1. Goldsmid, H. J. Introduction to Thermoelectricity. (Springer-Verlag, Berlin Heidelberg, 2010).

2. Tan, G., Zhao, L.-D. & Kanatzidis, M. G. Rationally designing high-performance bulk thermoelectric materials. Chem. Rev. 116, 12123–12149 (2016).

3. Gorai, P., Stevanović, V. & Toberer, E. S. C. Nat. Rev. Mater. 2, 1–16 (2017).

4. Jain, A., Shin, Y. & Persson, K. A. Computational predictions of energy materials using density functional theory. Nat. Rev. Mater. 1, 1–13 (2016).

5. Snyder, G. J. in Energy Harvesting Technologies, 325–336 (Springer, 2009).

6. Mettan, X. et al. Tuning of the thermoelectric figure of merit of CH3NH3MI3 (M = Pb,Sn) photovoltaic perovskites. J. Phys. Chem. C. 119, 11506–11510 (2015).

7. Saini, S., Baranwal, A. K., Yabuki, T., Hayase, S. & Miyazaki, K. Growth of halide perovskites thin films for thermoelectric applications. MRS Adv. 4, 1719–1725 (2019).

8. Haque, M. A., Kee, S., Villalva, D. R., Ong, W. L. & Baran, D. Halide perovskites: thermal transport and prospects for thermoelectricity. Adv. Sci. 7, 1903389 (2020).

9. Liu, T. et al. Enhanced control of self-doping in halide perovskites for improved thermoelectric performance. Nat. Commun. 10, 5750 (2019).

10. Jun, T. et al. Lead-free highly efficient blue-emitting Cs3Cu2I5 with 0D electronic structure. Adv. Mater. 30, 1804547 (2018).

11. Roccanova, R. et al. Near-unity photoluminescence quantum yield in blue-emitting Cs3Cu2Br5−xIx (0 ≤ x ≤ 5). ACS Appl. Electron. Mater. 1, 269–274 (2019).

12. Chen, H. et al. Multiple self-trapped emissions in the lead-free halide Cs3Cu2I5. J. Phys. Chem. Lett. 11, 4326–4330 (2020).

13. Lian, L. et al. Photophysics in Cs3Cu2X5 (X = Cl, Br, or I): highly luminescent self-trapped excitons from local structure symmetrization. Chem. Mater. 32, 3462–3468 (2020).

14. Snyder, G. J. & Toberer, E. S. Complex thermoelectric materials. Nat. Mater. 7, 105–114 (2008).

15. Haque, M. A. et al. A 0D lead-free hybrid crystal with ultralow thermal conductivity. Adv. Funct. Mater. 29, 1809166 (2019).

16. Akkerman, Q. A. & Manna, L. What defines a halide perovskite? ACS Energy Lett. 5, 604–610 (2020).

17. Huang, Y.-T., Kavanagh, S. R., Scanlon, D. O., Walsh, A. & Hoye, R. L. Z. Perovskite-inspired materials for photovoltaics–from design to devices. Nanotechnology 32, 132004 (2021).

18. Gillespie, R. & Nyholm, R. Inorganic stereochemistry. Quart. Rev. Chem. Soc. 11, 339–380 (1957).

19. Birch, F. Finite elastic strain of cubic crystals. Phys. Rev. 71, 809–824 (1947).

20. Jung, Y.-K. et al. Intrinsic doping limit and defect-assisted luminescence in Cs4PbBr6. J. Mater. Chem. A 7, 20254–20261 (2019).

21. Marzari, N. & Vanderbilt, D. Maximally localized generalized Wannier functions for composite energy bands. Phys. Rev. B 56, 12847–12865 (1997).

22. Souza, I., Marzari, N. & Vanderbilt, D. Maximally localized Wannier functions for entangled energy bands. Phys. Rev. B 65, 1–13 (2002).

23. Pizzi, G. et al. Wannier90 as a community code: New features and applications. J. Phys.: Condens. Matter 32, 165902 (2020).

24. Wahl, R., Vogtenhuber, D. & Kresse, G. SrTiO3 and BaTiO3 revisited using the projector augmented wave method: Performance of hybrid and semilocal functionals. Phys. Rev. B 78, 104116 (2008).

25. Skelton, J. M. et al. Influence of the exchange-correlation functional on the quasi-harmonic lattice dynamics of II-VI semiconductors. J. Chem. Phys. 143, 064710 (2015).

26. Hua, X., Hao, S. & Wolverton, C. First-principles study of vibrational entropy effects on the PbTe-SrTe phase diagram. Phys. Rev. Mater. 2, 095402 (2018).

27. Chaput, L. Direct solution to the linearized phonon boltzmann equation. Phys. Rev. Lett. 110, 265506 (2013).

28. Whalley, L. D., Skelton, J. M., Frost, J. M. & Walsh, A. Phonon anharmonicity, lifetimes, and thermal transport in CH3NH3PbI3 from many-body perturbation theory. Phys. Rev. B 94, 220301 (2016).

29. Skelton, J. M. et al. Lattice dynamics of the tin sulphides SnS2, SnS and Sn2S3: Vibrational spectra and thermal transport. Phys. Chem. Chem. Phys. 19, 12452–12465 (2017).

30. Haque, E. & Hossain, M. A. Origin of ultra-low lattice thermal conductivity in Cs2BiAgX6 (X = Cl, Br) and its impact on thermoelectric performance. J. Alloy. Compd. 748, 63–72 (2018).

31. Agne, M. T., Hanus, R. & Snyder, G. J. Minimum thermal conductivity in the context of diffuson-mediated thermal transport. Energy Environ. Sci. 11, 609–616 (2018).

32. Liu, T. et al. Unusual thermal boundary resistance in halide perovskites: a way to tune ultralow thermal conductivity for thermoelectrics. ACS Appl. Mater. Interfaces 11, 47507–47515 (2019).

33. Lee, W. et al. Ultralow thermal conductivity in all-inorganic halide perovskites. Proc. Natl. Acad. Sci. 114, 8693–8697 (2017).

34. Mukhopadhyay, S. et al. Two-channel model for ultralow thermal conductivity of crystalline Tl3VSe4. Science 360, 1455–1458 (2018).

35. Xia, Y., Pal, K., He, J., Ozoliņš, V. & Wolverton, C. Particlelike phonon propagation dominates ultralow lattice thermal conductivity in crystalline Tl3VSe4. Phys. Rev. Lett. 124, 065901 (2020).

36. Hellman, O., Abrikosov, I. & Simak, S. Lattice dynamics of anharmonic solids from first principles. Phys. Rev. B 84, 180301 (2011).

37. Tadano, T. & Tsuneyuki, S. Self-consistent phonon calculations of lattice dynamical properties in cubic SrTiO3 with first-principles anharmonic force constants. Phys. Rev. B 92, 054301 (2015).

38. Ruske, F. et al. Improved electrical transport in Al-doped zinc oxide by thermal treatment. J. Appl. Phys. 107, 013708 (2010).

39. Kim, C. E. et al. Effect of carrier concentration on optical bandgap shift in ZnO: Ga thin films. Thin Solid Films 518, 6304–6307 (2010).

40. Aydin, E. et al. Zr-doped indium oxide (IZRO) transparent electrodes for perovskite-based tandem solar cells. Adv. Funct. Mater. 29, 1901741 (2019).

41. Liu, Y. et al. Solution processed W-doped In2O3 thin films with high carrier mobility. Ceram. Int. 46, 2173–2177 (2020).

42. Ellmer, K. Past achievements and future challenges in the development of optically transparent electrodes. Nat. Photon. 6, 809–817 (2012).

43. Frost, J. M. Calculating polaron mobility in halide perovskites. Phys. Rev. B 96, 195202 (2017).

44. Ganose, A. M. et al. Efficient calculation of carrier scattering rates from first principles. Preprint at http://arxiv.org/abs/2008.09734 (2020).

45. Chung, D.-Y. et al. CsBi4Te6: A high-performance thermoelectric material for low-temperature applications. Science 287, 1024–1027 (2000).

46. Lin, J. et al. Unexpected high-temperature stability of β-Zn4Sb3 opens the door to enhanced thermoelectric performance. J. Am. Chem. Soc. 136, 1497–1504 (2014).

47. Zhao, L.-D. et al. Ultralow thermal conductivity and high thermoelectric figure of merit in SnSe crystals. Nature 508, 373–377 (2014).

48. Jong, U.-G., Yu, C.-J., Ri, J.-S., Kim, N.-H. & Ri, G.-C. Influence of halide composition on the structural, electronic, and optical properties of mixed $${{\rm{CH}}}_{3}{{\rm{NH}}}_{3}{\rm{Pb}}{({{\rm{I}}}_{1-x}{{\rm{Br}}}_{x})}_{3}$$ perovskites calculated using the virtual crystal approximation method. Phys. Rev. B 94, 125139 (2016).

49. Jung, Y.-K., Lee, J.-H., Walsh, A. & Soon, A. Influence of Rb/Cs cation-exchange on inorganic Sn halide perovskites: from chemical structure to physical properties. Chem. Mater. 29, 3181–3188 (2017).

50. Benin, B. M. et al. Highly emissive self-trapped excitons in fully inorganic zero-dimensional tin halides. Angew. Chem. Int. Ed. 130, 11499–11503 (2018).

51. Daniels, L. M. et al. Phonon-glass electron-crystal behaviour by A site disorder in n-type thermoelectric oxides. Energy Environ. Sci. 10, 1917–1922 (2017).

52. Gurunathan, R., Hanus, R. & Snyder, G. J. Alloy scattering of phonons. Mater. Horiz. 7, 1452–1456 (2020).

53. Xie, H. et al. All-inorganic halide perovskites as potential thermoelectric materials: dynamic cation off-centering induces ultralow thermal conductivity. J. Am. Chem. Soc. 142, 9553–9563 (2020).

54. Togo, A. & Tanaka, I. First principles phonon calculations in materials science. Scr. Mater. 108, 1–5 (2015).

55. Togo, A., Chaput, L. & Tanaka, I. Distributions of phonon lifetimes in Brillouin zones. Phys. Rev. B 91, 094306 (2015).

56. Skelton, J. M. Phono3py-Power-Tools. https://github.com/skelton-group/Phono3py-Power-Tools (accessed 9 February 2021).

57. Madsen, G. K. H., Carrete, J. & Verstraete, M. J. BoltzTraP2, a program for interpolating band structures and calculating semi-classical transport coefficients. Comput. Phys. Commun. 231, 140–145 (2018).

58. Kohn, W. & Sham, L. J. Self-consistent equations including exchange and correlation effects. Phys. Rev. 140, A1133–A1138 (1965).

59. Kresse, G. & Furthmüller, J. Efficient iterative schemes for abinitio total-energy calculations using a plane-wave basis set. Phys. Rev. B 54, 11169–11186 (1996).

60. Kresse, G. & Furthmüller, J. Efficiency of ab-initio total energy calculations for metals and semiconductors using a plane-wave Basis Set. Comput. Mater. Sci. 6, 15–50 (1996).

61. Kresse, G. & Joubert, D. From ultrasoft pseudopotentials to the projector augmented-wave method. Phys. Rev. B 59, 1758–1775 (1999).

62. Blöchl, P. E. Projector augmented-wave method. Phys. Rev. B 50, 17953–17979 (1994).

63. Perdew, J. P. et al. Restoring the density-gradient expansion for exchange in solids and surfaces. Phys. Rev. Lett. 100, 136406 (2008).

64. Heyd, J., Scuseria, G. E. & Ernzerhof, M. Hybrid functionals based on a screened coulomb potential. J. Chem. Phys. 118, 8207–8215 (2003).

65. Krukau, A. V., Vydrov, O. A., Izmaylov, A. F. & Scuseria, G. E. Influence of the exchange screening parameter on the performance of screened hybrid functionals. J. Chem. Phys. 125, 224106 (2006).

## Acknowledgements

We thank J.M. Skelton and A.M. Ganose for useful discussions on phonon and electron transport. Via our membership of the UK’s HEC Materials Chemistry Consortium, which is funded by EPSRC (EP/L000202), this work used the ARCHER UK National Supercomputing Service (http://www.archer.ac.uk). This work was supported by a National Research Foundation of Korea (NRF) grant funded by the Korean government (MSIT) (No. 2018R1C1B6008728). This research was partially supported by the Graduate School of YONSEI University Research Scholarship Grants in 2019.

## Author information

Authors

### Contributions

Y.-K.J. performed the calculations, collected, and analyzed the data. I.T.H. and Y.C.K. initialized the project. A.W. led the project. All authors contributed to discussing the results and writing the manuscript.

### Corresponding authors

Correspondence to Yong Churl Kim or Aron Walsh.

## Ethics declarations

### Competing interests

The authors declare no competing interests.

Publisher’s note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

## Rights and permissions
{"ft_lang_label":"__label__en","ft_lang_prob":0.8341812,"math_prob":0.94563085,"size":32778,"snap":"2022-27-2022-33","text_gpt3_token_len":8690,"char_repetition_ratio":0.1384329,"word_repetition_ratio":0.013851561,"special_character_ratio":0.26136434,"punctuation_ratio":0.18557021,"nsfw_num_words":0,"has_unicode_error":false,"math_prob_llama3":0.97742087,"pos_list":[0,1,2,3,4],"im_url_duplicate_count":[null,null,null,1,null],"WARC_HEADER":"{\"WARC-Type\":\"response\",\"WARC-Date\":\"2022-08-16T06:14:02Z\",\"WARC-Record-ID\":\"<urn:uuid:824bf959-5408-4b58-aa22-c4dea310cfa3>\",\"Content-Length\":\"336014\",\"Content-Type\":\"application/http; msgtype=response\",\"WARC-Warcinfo-ID\":\"<urn:uuid:ebead8ed-6fc5-4d04-ba18-0a9b4a58e174>\",\"WARC-Concurrent-To\":\"<urn:uuid:eaccac2b-5e66-4e36-86e1-c5d783394463>\",\"WARC-IP-Address\":\"146.75.36.95\",\"WARC-Target-URI\":\"https://www.nature.com/articles/s41524-021-00521-9?error=cookies_not_supported&code=ccb7f751-215f-4b7c-b50a-bfd63e6f38f0\",\"WARC-Payload-Digest\":\"sha1:TLHYMUW67J27GYYE7XEGVJRRPMK2TCPY\",\"WARC-Block-Digest\":\"sha1:5AKIOY6KI2A7BHE5BFNBZJREFSQC7YQA\",\"WARC-Identified-Payload-Type\":\"text/html\",\"warc_filename\":\"/cc_download/warc_2022/CC-MAIN-2022-33/CC-MAIN-2022-33_segments_1659882572221.38_warc_CC-MAIN-20220816060335-20220816090335-00183.warc.gz\"}"}