Tobias
Haeger
a,
Ralf
Heiderhoff
*ab and
Thomas
Riedl
ab
aInstitute of Electronic Devices, University of Wuppertal, Wuppertal, Germany. E-mail: heiderho@uni-wuppertal.de
bWuppertal Center for Smart Materials & Systems, University of Wuppertal, Wuppertal, Germany
First published on 7th September 2020
Aside from photovoltaics, metal-halide perovskite semiconductors have also emerged as an attractive platform for LEDs and even lasers. For all of them, performance and operational stability are strongly influenced by thermally activated processes. As a result, studying the thermal properties of halide perovskites has become increasingly important. In addition to thermal conductivity measurements, thermal diffusivity and heat capacity studies are becoming more and more relevant for the dynamic operation of these devices. Here, we will review the current state of knowledge on the thermal properties of metal halide perovskites, including both experimental findings and theoretical insights. We will highlight the impact of the perovskite dimensionality (3D, 2D, 0D) on the thermal properties and how these properties change across the various phase transitions of these perovskites. Our review will also discuss thermal measurement methods, which can be favourably used to determine the thermal properties of both bulk and thin film samples. Advanced comprehensive tools and strategies for systematic thermal analyses and property optimizations are discussed.
The outstanding (opto-)electronic properties of these lead-halide perovskites have been re-discovered during the past decade. Most impressively, the power conversion efficiency of lab-scale perovskite solar cells (PSCs) has rapidly increased to a level of 25.2%,6 which rivals that of crystalline silicon. In addition, a wide range of other applications such as light emitting diodes (LEDs)7 and lasers8,9 have been explored. Moreover, this family of materials also holds great promise for thermoelectric generators.10 Here, the extremely low thermal conductivities of lead-halide perovskites11,12 may be advantageous, while for the vast majority of opto-electronic devices heat transfer may be an issue and thermal management requires particular attention. It is known that the crystalline structure of some representatives undergoes phase transitions T slightly above room temperature (RT); for example, CsPbCl3 is in the tetragonal phase at RT, whereas above 46.9 °C it becomes cubic.13 At early times high coefficients of thermal expansion were reported.14 Both the crystalline phase and thermal expansion influence electronic properties (for instance, carrier mobility, directionality) and carrier dynamics.15 More severely, substantial decomposition effects for MA-based perovskites have been found to occur already below 100 °C even in an inert atmosphere.16 Consequently, heating as a result of a lack of thermal management may turn out detrimental in various aspects.17–20 This holds particularly for (concentrator) solar cells and LEDs and the prospects of electrically or continuous wave operated perovskite lasers.21,22
Thus, in a wide range of these organic–inorganic and all-inorganic perovskites, the acquisition of extreme thermal properties has become more and more important, because both the lifetime and performance are influenced by temperature or temperature gradients.23,24 As detailed in the following, considerable research efforts have provided access to thermo-elastic and thermo-electric properties, as well as to intrinsic thermal conductivities, thermal diffusivities, and heat capacities.
It is important to note that thermal properties strongly depend on temperature themselves and phase transitions may cause significant and abrupt changes. By a similar token, the thermal properties of 0D and 2D materials may significantly differ from those of the bulk (3D).
In the following sections we will review the thermal properties of metal-halide perovskites and the experimental techniques that have been used to assess them. The experimental data will be complemented by the results of theoretical work that is indispensable to understand the particular thermal properties of this class of materials and to identify levers to overcome limitations.
(1) |
Fig. 1 An ultralow thermal conductivity can be obtained from the calculated phonon dispersion curves in the low-frequency domain of (a) tetragonal MAPbI3 at room temperature, while the thermal conductivity of (b) the cubic phase (>54 °C33) is slightly larger, because the phonon group velocities are given by the gradient.31 Copyright Wiley-VCH Verlag GmbH & Co. KGaA. Reproduced with permission. |
Only their absolute group velocities, mean-free paths, and lifetimes obtained from these relations in all directions provide access to the anisotropy of thermal transport. In addition, the phonon dispersion is essential for the description of heat dissipation in the temperature-dependent crystalline structure and the dimensionality given by its composition AxByXz in these perovskite semiconductor types.
The electron/phonon heat capacity can be calculated by taking the derivative of the total electron/phonon energy density with respect to the electron/phonon temperature: in other words, the amount of heat, ΔQ(), that must be added to a system in order to increase its temperature by ΔT(). A first-order phase transition in a crystalline structure is described by the Clausius–Clapeyron equation. This results in the divergence of the heat capacity csp(T) or the thermal expansion coefficient α(T), because the first derivative of the free energy is discontinuous. In contrast, a second-order (continuous) phase transition, which is described by the Ehrenfest equation, occurs with diverging correlation length as the crystalline phase transition temperature is approached. For material systems undergoing this phase transition due to continuous entropy changes such as magnetic and ferroelectric, the latent heat is zero, and no discontinuity in csp(T) and α(T) occurs.
Thus, the detection of the thermal properties λ, csp, cvol., a, and α depending on their temperature is challenging, because comprehensive thermal studies allow access to a variety of material properties and are indispensable for the heat management of semiconducting perovskite devices.
A very low thermal conductivity of 0.59 W (m K)−1 was found in the tetragonal phase at room temperature of MAPbI3, attributed to the low group velocity of acoustic phonons due to their low elastic stiffness and the strong anharmonicity, whereas a much higher thermal conductivity was predicted for the pseudo-cubic phase (1.80 W (m K)−1 at 57 °C).35 It was argued that the lower thermal conductivity in the tetragonal phase was due to the shorter lifetimes of the optical phonons and the smaller group velocity of acoustic phonons than those in the pseudo-cubic phase. In addition, the thermal conductivity of MAPbI3 was predicted to be below 1 W (m K)−1, and as low as 0.31 W (m K)−1 at room temperature.31 Such ultralow thermal conductivity, as a result of the enhanced phonon–phonon scattering from highly-overlapped phonon branches, hints to short phonon lifetimes <100 ps and consequently short phonon mean-free paths <10 nm. Compared with FAPbI3, a 3 times higher carrier-phonon relaxation rate is observed, while that in a cesium-based system is 10 times higher.36 In general, the thermal conductivities of lead-based perovskites are low, due to the low sound velocity of the PbX lattice.37 In this context it should be noted that the contribution of the average phonon group velocity to the thermal conductivity has been discussed elsewhere.30 Nevertheless, the degrees of freedom of the A-site cation also provide a route to tune λ.
Hence the role of the rotations of MA+ cations in thermal properties has been studied by molecular dynamics simulations in addition.38 The vibrational energy shifts as the temperature increases, where anharmonic couplings occur. Coupled translational and rotational motions of CH3NH3+ were found to interact with the Pb–I cages and afford couplings between the isolated lattice vibrations, which suppress the phonon transport in MAPbI3.
The anisotropy in the tensor of thermal conductivity for the tetragonal and orthorhombic phases was found to originate from the preferential orientations of MA+.31 It has been claimed that among all the atomistic interactions, electrostatic interactions dominate thermal conductivity in ionic MAPbI3 crystals. Hence, hybrid perovskites are sometimes also considered as “soft materials”.
Recent MDS show that the A-site component of two-dimensional metal halide perovskites induces a reduction of the stiffness and sound velocities along with giving rise to vibrational modes in the 5–15 THz range that are absent in the three-dimensional counterparts.32 This leads to ultralow thermal conductivities in the range of 0.10 to 0.19 W (m K)−1 depending on the molecular substructure of ammonium cations and owing to the weaker interactions in the layered structures.
Thermal transport investigations of perovskite alloys in particular, which have been successfully developed to design various layered structures,39,40 have been performed rarely, with the first reports published just recently.41–43 The electrical behavior, optical absorption, and heat conduction or convection are considered to give insight into heat dissipation across solar cells. These first inspections of heat accumulations and gradients in transparent front electrode and back electrode materials are promising for the higher thermal stabilities of these devices, since they also take the heat dissipation in air into account.
Note that dynamical heat transport analyses have not been carried out so far. Accordingly, their experimental thermal property studies for reliable dynamic device operation will be considered in addition in the following.
Perovskite | SC/PC | λ [W (m K)−1] | c s [J (kg K)−1] | c V [MJ (m3 K)−1] | a [mm2 s−1] | s [m s−1] | S [μV K−1] | Measurement technique | Measured [K] | Crystal phase | Perovskite manufacturing | Year | Ref. |
---|---|---|---|---|---|---|---|---|---|---|---|---|---|
MAPbI3 | SC | 0.5 | Steady-state method | RT | Tetragonal | ITC | 2014 | 12 | |||||
MAPbI3 | SC | 0.31 or <1 for all phases | Molecular dynamics simulation | RT | Orthorhombic, tetragonal, cubic | 2016 | 31 | ||||||
MAPbI3 | SC | 0.30–0.42 (i) | 920 ± 91 | Laser flash | 298–425 | Tetragonal–cubic | ITC/vapor diffusion | 2017 | 59 | ||||
MAPbI3 | SC | 0.34 ± 0.12 | SThM/3ω method | RT-313 | Tetragonal | ITC | 2017 | 61 | |||||
MAPbI3 | SC | 0.14 ± 0.02 | Micro-photoluminescence | RT | Tetragonal | Two-step CVD grown | 2017 | 220 | |||||
MAPbI3 | SC | 0.34 ± 0.08 | 311 calculated | 1.28221 | 0.28 calculated | 1390+43−66 | Frequency domain thermo-reflectance | RT | Tetragonal | Vapor diffusion | 2017 | 30 | |
MAPbI3 | SC | 0.59 | Molecular dynamics simulation | 160–330 | Tetragonal | 2016 | 35 | ||||||
MAPbI3 | SC | 1.8 | Molecular dynamics simulation | >330 | Pseudocubic | 2016 | 35 | ||||||
MAPbI3 | SC | 0.30 | Laser flash | RT | Tetragonal | ITC | 2018 | 63 | |||||
MAPbI3 | SC-NW | 0.22 | Microbridge – device | RT | Tetragonal | 222 | 2018 | 49 | |||||
MAPbI3 | PC | 0.6 | Longitudinal steady state technique | RT | Vacuum quenching | 2017 | 53 | ||||||
MAPbI3 | PC | 0.33 ± 0.12 | SThM/3ω method | RT-318 | Tetragonal | Spincoating/NIL | 2017 | 61 | |||||
MAPbI3 | PC | 0.3 | Steady-state method | RT | Tetragonal | ITC/glueing | 2014 | 12 | |||||
MAPbI3 | PC | 11.2 ± 0.8 | Time-domain thermal reflectance | RT | Tetragonal | 223 | 2016 | 62 | |||||
MAPbI3 | PC | 0.38 | 1350 ± 50 | Laser flash | 300–450 K | Tetragonal cubic | Mechanical alloying via high energy ball milling, followed by hot pressing | 2019 | 76 | ||||
MAPbI3 | PC | 0.36 ± 0.01 | Laser flash | 300 K | Tetragonal | — | 2019 | 77 | |||||
MAPbI3 | PC | 0.32 ± 0.03 | 6500/5500 | 3ω method | RT | Tetragonal | Two-step sequential vacuum | 2019 | 224 | ||||
Deposition | |||||||||||||
SVD | |||||||||||||
MAPbI3 | PC | 0.31 ± 0.03 for film thicknesses 65 ± 3 nm | 3ω method | RT | Coevaporation | 2019 | 225 | ||||||
MAPbI3 | PC | 0.44 ± 0.03 for film thicknesses 80 ± 5 nm | 3ω method | RT | Coevaporation | 2019 | 225 | ||||||
MAPbI3 | PC | 0.59 ± 0.04 for film thicknesses 100 ± 12 nm | 3ω method | RT | Coevaporation | 2019 | 225 | ||||||
MAPbI3 | PC | 820 (dark) | ITC | 2015 | 56 | ||||||||
540 (light) | |||||||||||||
MAPbI3 | PC | Average 0.4 | 1.200 ± 300 | TDTR | RT | Spincoating | 2016 | 60 | |||||
MAPbBr3 | SC | 0.44 ± 0.08 | SThM/3ω method | RT | ITC | 2017 | 61 | ||||||
MAPbBr3 | SC | 0.51 ± 0.12 | 378 calculated | 1.45221 | 0.35 calculated | 1717+79−1216 | Frequency domain thermo-reflectance | RT | Cubic | Vapor diffusion | 2017 | 30 | |
MAPbBr3 | SC | 0.37 | Laser flash | RT | Cubic | ITC | 2018 | 63 | |||||
MAPbBr3 | SC-NW | 0.32 | Microbridge – device | RT | Cubic | 222 | 2018 | 49 | |||||
MAPbBr3 | PC | 0.39 ± 0.05 | SThM/3ω method | RT | Spincoating/NIL | 2017 | 61 | ||||||
MAPbCl3 | SC | 0.73 ± 0.18 | 511 calculated | 1.62221 | 0.45 calculated | 2194+71−108 | Frequency domain thermo-reflectance | RT | Cubic | Vapor diffusion | 2017 | 30 | |
MAPbCl3 | SC | 0.50 ± 0.05 | SThM/3ω method | RT | ITC | 2017 | 61 | ||||||
MAPbCl3 | SC | 0.52 | Laser flash | RT | Cubic | ITC | 2018 | 63 | |||||
MAPbCl3 | PC | 0.50 ± 0.12 | SThM/3ω method | RT | Spincoating/NIL | 2017 | 61 | ||||||
CsPbI3 | SC-NW | 0.45 ± 0.05 (κl) | NW thermometry technique | RT | — | 2017 | 48 | ||||||
CsPbBr3 | PC | 0.43 ± 0.03 | 290 ± 90 calculated | 1.3 ± 0.4 | 0.3 ± 0.1 | SThM/3ω method | RT | Orthorhombic | Spincoating/NIL | 2019 | 50 | ||
CsPbBr3 | SC | 0.46 ± 0.12 | 220 calculated | 1.06 | 0.43 calculated | 1390+96−74 | Frequency domain thermo-reflectance | RT | Orthorhombic | Vapor diffusion | 2017 | 30 | |
CsPbBr3 | SC | 0.514 | 280 | 1.29 calculated | 0.399 | Pulse transient method | 37 °C | 2008 | 217 | ||||
CsPbBr3 | SC-NW | 0.36 | Microbridge – device | RT | Orthorhombic | 222 | 2018 | 49 | |||||
CsPbBr3 | SC-NW | 0.42 ± 0.04 (κl) | NW thermometry technique | RT | — | 2017 | 48 | ||||||
CsPbCl3 | SC | 0.530 | 272 | 1.12 calculated | 0.475 | Pulse transient method | 25 °C | 226 | 2008 | 217 | |||
CsPbCl3 | SC | 0.55 ± 0.03 | 272 ± 8 | 1.143 calculated | 0.49 ± 0.02 | Pulse transient method | 50 °C | 226 | 2000 | 11 | |||
CsPbCl3 | PC | 0.49 ± 0.04 | 220 ± 70 calculated | 0.9 ± 0.1 | 0.5 ± 0.1 | SThM/3ω method | RT | Monoclinic | Spincoating/NIL | 2020 | 51 | ||
FAPbBr3 | SC | 0.49 ± 0.12 | 444 calculated | 1.69 | 0.29 calculated | 1333+51−73 | Frequency domain thermo-reflectance | RT | Cubic | Vapor diffusion | 2017 | 30 | |
MASnI3 | PC | 0.09 | −720 | Steady state method | ITC | 2015 | 56 | ||||||
CsSnI3 | SC-NW | 0.38 ± 0.04 (κl) | NW thermometry technique | RT | — | 2017 | 48 | ||||||
CsSnI3 | SC-NW | 0.57 | 79 | NW thermometry technique | 300 K | — | 2017 | 48 | |||||
CsSnI3 | PC | 0.36 | 115 | 3ω method | RT | Spincoating | 2020 | 70 | |||||
(MA)3Bi2I9 | PC | 0.21 | 2600 ± 200 | Laser flash | 300–450 K | Mechanical alloying via high energy ball milling, followed by hot pressing | 2019 | 76 | |||||
(MA)3Bi2I9 | PC | 0.23 ± 0.02 | Laser flash | 300 K | — | 2019 | 77 | ||||||
Cs0.05FA0.95Pb(I0.83Br0.17)3 | PC | 0.26 ± 0.03 | 0.74 ± 0.11 calculated | 0.35 ± 0.05 | Thermo-reflectance | RT | Spincoating | 2020 | 85 | ||||
Low-dimensional perovskites | |||||||||||||
[Mn(C2H6OS)6]I4 | SC | 0.15 ± 0.01 | 580 ± 20 | 0.12 ± 0.01 | Laser flash differential calorimeter | RT | ITC | 2019 | 214 | ||||
CsPb2Br5 | PC | 0.33 ± 0.02 | 180 ± 60 calculated | 1.1 ± 0.3 | 0.3 ± 0.1 | SThM/3ω method | RT | Spincoating/NIL | 2019 | 50 | |||
CsPb2Cl5 | PC | 0.40 ± 0.02 | 0.7 ± 0.2 | 0.6 ± 0.1 | SThM/3ω method | RT | Spincoating/NIL | 2019 | 215 | ||||
Cs4PbCl6 | PC | 0.30 ± 0.03 | 150 ± 80 calculated | 0.5 ± 0.2 | 0.5 ± 0.1 | SThM/3ω method | RT | Spincoating/NIL | 2020 | 51 | |||
BA2PbI4 | PC | 0.10 ± 0.02 horizontal | 1.22 horizontal | Thermo-reflectance molecular dynamics simulation | RT | Spincoating | 2020 | 32 | |||||
0.18 ± 0.04 vertical | 1.18 vertical | ||||||||||||
BA2PbBr4 | PC | 0.11 ± 0.02 horizontal | 1.37 horizontal | Thermo-reflectance molecular dynamics simulation | RT | Spincoating | 2020 | 32 | |||||
0.19 ± 0.04 vertical | 1.37 vertical | ||||||||||||
PEAPbI4 | PC | 0.15 ± 0.03 | 1.19 | Thermo-reflectance molecular dynamics simulation | RT | Spincoating | 2020 | 32 | |||||
PEAPbBr4 | PC | 0.13 ± 0.03 | 1.22 | Thermo-reflectance molecular dynamics simulation | RT | Spincoating | 2020 | 32 | |||||
BA2MA3 | PC | 0.13 ± 0.03 vertical | Thermo-reflectance | RT | Spincoating | 2020 | 32 | ||||||
Pb4I13 | 0.11 ± 0.03 random |
The first detailed thermo-mechanical experimental investigations of both CsPbCl3 and CsPbBr3 perovskites were made in the 1970s in order to understand the nature of their structural phase transitions.44–47 But the corresponding specific heat csp and thermal diffusivity a values of the different phases as well at their transitions were obtained in transient measurements on single crystals much later (which will be discussed in more detail later in the Thermal properties across phase transitions section).11 In addition, thermal conductivity studies have been carried out for all-inorganic halide perovskite nanowires, CsPbI3 and CsPbBr3.48,49 The authors generally concluded that nanostructures or crystal grains with dimensions smaller than 100 nm significantly reduce the thermal transport because phonons have a mean-free path on the same order. Thus, somewhat higher thermal conductivities are reported for CsPbBr3 single crystals (0.46 W (m K)−1)30 and thin films containing large crystals (0.43 W (m K)−1)50 in comparison to nanowires (0.36 W (m K)−1).49 These thermal conductivities are somewhat lower than those recently found for CsPbCl3 thin films with large laterally extended crystallites (0.49 W (m K)−1).51 Note that, since the phonon mean-free-path is also reduced by defects, thermal conductivities are lower near grain boundaries in thin films, as shown in local thermal conductivity analyses with a high spatial resolution by an AFM based technique.50 This technique allows determining the local thermal diffusivity, aside from the local thermal conductivity; thus the local volumetric heat capacity can be calculated afterwards. Average a() of 0.3 mm2 s−1 and 0.5 mm2 s−1 as well as average cvol.() of 1.3 J cm−3 K−1 and 0.9 J cm−3 K−1 were determined for CsPbBr3 and CsPbCl3 at RT respectively.50,51
The findings of these studies suggested that heat dissipation in dynamically operated halide perovskite optoelectronic devices will be very inefficient. As such, the thermal management of light-emitting diodes and lasers based on halide perovskites is a challenge and requires particular attention (see strategies to enhance the thermal properties of transport layers in perovskite devices).
The specific heat capacities of MAPbX3 (X = I, Br, and Cl) perovskites were studied by Onoda-Yamamuro et al. in 1990 to be 19.0 J (K mol)−1, 8.2 J (K mol)−1, and 10.0 J (K mol)−1 at RT for X = I, Br, and Cl, respectively.52 Upon heating the MAPbI3 samples above room temperature, a tetragonal to pseudo-cubic polymorph transition (57 °C) was observed that led to a reduction of csp to 9.7 J (K mol)−1. Pisoni et al. performed the first thermal conductivity measurements of MAPbI3 for temperatures below RT.12 The analysis was performed by a steady state bar technique. The room temperature thermal conductivity λ was determined to be 0.5 W (m K)−1 and 0.3 W (m K)−1 for a single crystal and a polycrystalline sample, respectively. The same thermal conductivity was found for a polycrystalline sample when measuring the effect of the methylammonium ion on phonon scattering.53 The low thermal conductivity of the MAPbI3 single crystal was mainly attributed to the rotational motion of the MA cations54 and has attracted intense attention due to its possible thermoelectric applications.55 A high Seebeck coefficient S of 0.82 mV K−1,56 a decent charge carrier mobility in the range of 5–10 cm2 (V s)−1 for electrons and 1–5 cm2 (V s)−1 for holes,57 and a high carrier diffusion length >1 μm58 were reported. The thermal conductivity of a MAPbI3 perovskite was also evaluated using laser flash systems. Values of around 0.30–0.42 W (m K)−1 were observed for a single crystal in the temperature range of RT to 150 °C.59 Similar values were found for films using the time-domain thermo-reflectance (TDTR) method.60 In striking contrast to all other reports, an extremely high thermal conductivity of densely packed MAPbI3 films of 11.2 ± 0.8 W (m K)−1 at room temperature was claimed by Chen et al.62
On the other hand, consistent thermal conductivities were reported for single crystals (0.34 W (m K)−1) and densely packed thin perovskite films containing large crystals (0.33 W (m K)−1) using 3ω-scanning near-field thermal microscopy.61 In addition, the temperature-dependent thermal conductivity of MAPbI3 was first measured around the tetragonal–cubic phase transition temperature. It could be demonstrated that the thermal conductivity of a cleaved single crystal in the cubic structure increases to 1.1 W (m K)−1, in excellent agreement with the results of EMD simulations,31 which, however, could not be verified from other heat transport measurements.63 Furthermore, it could experimentally be shown that MAPbX3 (X = I, Br, and Cl) single crystals and thin films with the decreasing atomic number of the halide show a slightly increasing average thermal conductivity due to their significantly increasing elastic stiffness (see Fig. 2).30,61 Besides the above mentioned halide variations, studies have been carried out on FAPbBr3,30 on targeted thermo-electrics such as CsSnI3,48,70–72 and on estimated lead-free solid-state organic–inorganic halide perovskite solar cell materials such as MASnI373,74 and (MA)3Bi2I9,75,76 with which the influences of A- and B-cation substitutions on thermal transports were examined (see Fig. 3).
Fig. 2 Perovskites with the decreasing atomic number of the halide have slightly higher average thermal conductivities.30,61 |
The obtained λ values of 0.74–0.28 W (m K)−1 for CsSnI3 vary significantly for single crystals, nanowires, and polycrystalline films due to the enhanced phonon scattering by defects/vacancies and reduced crystal sizes. Polycrystalline MASnI3 seems to have a glassy behavior, which might be the reason for its ultralow thermal conductivity of 0.09 W (m K)−1.74 Lastly, it must be pointed out that the thermal properties of 0D and 2D perovskite materials significantly differ from those of the bulk (3D). 0D Bi based perovskites, for instance, have a lower thermal conductivity of 0.21 W (m K)−1 than Pb based systems, which is due to their soft phonon dispersion and weak bonds.77 Since these materials can also be useful for a wide range of applications,78 the Dimensionalities of metal-halide perovskites section will highlight their impact on thermal transport.
In summary, the mean values of thermal properties obtained from molecular dynamics calculations and experimental investigations of representative studies are listed in Table 2. This table is intended to provide an overview for the assessment of a suitable thermal management of halide perovskite devices, since these properties of these compounds are very different from those of other materials. For example, thermally insulating materials commonly have large heat capacities, whereas those with high thermal conductivities, such as metals, have low ones. Thus, afterwards, thermal measurement methods used will also be discussed briefly with regard to their applicability to determine these very low thermal quantities of both bulk and thin film samples. Advanced comprehensive tools for systematic thermal analyses are highlighted in this context.
Material | λ [W (m K)−1] | a [mm2 s−1] | C s [J (g K)−1] | C v [J (cm3 K)−1] | s [m s−1] |
---|---|---|---|---|---|
MAPbI3 | 0.35 ± 0.02 | 0.28 ± 0.07 | 0.31 ± 0.08 | 1.28 | 1295 ± 134.35 |
MAPbBr3 | 0.39 ± 0.07 | 0.35 ± 0.08 | 0.38 ± 0.09 | 1.45 | 1717+79−126 |
MAPbCl3 | 0.51 ± 0.11 | 0.45 ± 0.09 | 0.51 ± 0.13 | 1.62 | 2194+71−108 |
CsPbI3 | 0.45 ± 0.05 | — | — | — | — |
CsPbBr3 | 0.43 ± 0.06 | 0.40 ± 0.07 | 0.28 ± 0.04 | 1.29 ± 0.14 | 1361+74−96 |
CsPbCl3 | 0.53 ± 0.03 | 0.49 ± 0.01 | 0.27 ± 0.03 | 1.12 ± 0.13 | — |
Si | 15664 | 78 | 0.765 | 1.6 | 8500 |
GaAs | 4564 | 26 | 0.3366 | 1.7 | 5300 |
CZTS | 4.7267 | 2067 | 0.0567 | 0.24 | — |
CZTSe | 4.2667 | 2067 | 0.0467 | 0.21 | — |
Organic | 0.368 | 0.15 | 1.7 | 269 | — |
The common theoretical framework of these methods is based on the Fourier law:26,27
= −λ·gradT() | (2) |
The steady-state measurement technique can be described in the simplest way by the thermal characterization based on the analysis in space. The schematic of a thermal resistivity measurement setup in a standard four-point configuration is shown in Fig. 4. An electrical heater, where () is equal to the electrical power, is connected to a sample, represented by the black rectangle, using epoxy. The sample is connected to a calibrated reference material, which is glued to a heat sink for heat dissipation. The given temperature distribution is obtained using differential thermocouples fixed onto the sample and onto the reference surfaces to determine the thermal conductivity. This common technique was the very first one introduced to determine the thermal conductivity of materials and has been improved successively.80,81 Of course, it must be ensured that the thermal measurement itself does not lead to artifacts. For instance, the experiments must be carried out in the dark to avoid unwanted photo-induced effects. Special care must be taken in the case of polycrystalline perovskite structures, as the thermal transport is significantly influenced by the presence of many grain boundaries. Thus, large single crystal perovskites are initially preferred in order to explore the thermal properties in detail. Furthermore, solvents of glues used to mount heaters or thermocouples to the sample may cause degradation.
Fig. 4 Schematic of a thermal resistivity measurement setup in a standard four-point configuration. Reprinted with permission from ref. 12. Copyright (2014) American Chemical Society. |
Thus, a non-contact optical flash method was already established in 1961,82 which is widely used in commercial thermal conductivity measurement systems. Its application to the study of microelectronic materials by the detection of thermal waves through changes in their optical reflectance was primarily demonstrated in 1985.83 The initial thermal diffusivity measurements of thin films by means of an AC calorimetric method were performed in the same year.84 Similar to these analyses, a modulated thermo-reflectance microscopy technique (shown in Fig. 5) has recently been successfully applied for measurements of the thermal conductivity and the thermal diffusivity in mixed-cation perovskites.85
Fig. 5 Principle of a modulated photoreflectance experimental setup. Reproduced from ref. 85, with the permission of AIP Publishing. |
A modulated pump beam (532 nm) is focused onto a sample's surface to provide a modulated local heat source in a metallic film deposited on top of the perovskite film. The heat diffuses laterally on the surface (r direction) and perpendicularly to the sample's surface (z direction). A periodic modulation of the surface reflection coefficient ΔR as a result of the periodically induced surface temperature modulation is sensed by a CW probe beam (488 nm) and a lock-in amplifier.
Theoretical thermal models86 are used which describe the magnitude (ΔR) and the phase (ΔΦ) of the modulated photoreflectance signal measured as a function of the distance to the pump beam. To get access to the thermal conductivity and diffusivity of the perovskite film in the r direction and the z direction (depending on the crystalline orientation of the film) the different layer systems underneath must be analysed beforehand. Moreover, photo-induced degradation,87 slight changes in the optical properties of MAPbI3 perovskite across the tetragonal to cubic transition,88 and possible degradation under intense optical excitation89 must be taken into account. For the sake of completeness, it should be mentioned that, besides these modulated heating techniques, thermal diffusivity measurements using a pulsed laser heating source for transient picosecond time-domain thermal-reflectance (TDTR) were presented.90 Although the results of the first TDTR measurements of MAPbI3 thin films62 significantly differ from other results (thermal conductivity of MAPbI3 films of 11.2 W (m K)−1 at RT), this method is promising because a theoretical description for frequency modulated time delay photo-thermal and photo-acoustic wave spectroscopies is given.91 Note that the optical generation of excess carriers due to laser excitation may significantly affect the thermal conductivity measurements. Thus, a possible explanation for the higher thermal conductivity determined by this technique might be linked to the contribution of the optically generated charge carriers to the thermal transport. The previously only theoretically predicted anisotropies in perovskites correlated with direction-dependent variations of the phonon vibrational density of states31 could be recorded by a differential photo-acoustic method to measure the thermal conductivity of thin films.92 Details of heat flow analysis in layered structures using a picosecond optical pump-and-probe technique can be found elsewhere.93
In addition to the space and time domains, thermal properties have been measured in the frequency domain. Nowadays, frequency-domain thermo-reflectance (FDTR) is used besides TDTR to simultaneously measure in-plane and cross-plane thermal properties.94 This technique can determine the thermal conductivity and heat capacity of a sample at the same time and it combines the advantages of TDTR with the relative experimental simplicity of modulated photo-thermal methods. Using FDTR, the thermal conductivities of various single crystal lead-halide perovskites were found to be proportional to ca-cubic·vs/3, where ca-cubic is the volumetric heat capacity of acoustic phonons based on a cubic unit cell and vs is the speed of sound.30 An average acoustic phonon mean-free path a = 4.3 nm was determined. Thus, it could be shown that this technique is very promising for the determination of thermal properties.
Further advantages of these various techniques of thermal wave microscopy and their application for testing of micro- and nano-structured materials and devices are given in the review by Galovic et al.95 A review of other commonly used differential scanning calorimetry techniques is given by Gill et al.96 Note that, in addition to the drawbacks of the so far presented far-field thermal measurement methods, it is well known that thermal transport strongly depends on the dimensions of the materials and is of growing interest upon shrinking the film thickness.97–99 In this regard, innovative near-field thermal microscopy techniques open a non-destructive route to determine the thermal properties of perovskites beyond reach for classical thermal measurement techniques.
Good candidates are scanning probe microscopy based techniques with the opportunity to increase the spatial resolution of established analysis techniques on semiconducting perovskite materials and devices.100 The first atomic force microscopy based measurement technique for local thermal conductivity, called Scanning Thermal Microscopy (SThM), was invented in 1992.101 Shortly afterwards, the heat source and the temperature sensor were combined in one single element by means of resistive thermal probes for local heat transport analyses with high spatial resolution.102,103 This combined probe made SThM more and more important for a wide range of thermal property studies of crystalline solids, organic materials, and material compounds for (opto-)electronic applications because a range of thermal properties, such as heat transport, glass transition, melting, recrystallization, and thermal decomposition, can be obtained locally with high spatial resolution.104–109
The established macroscopic 3ω-technique110 by means of a modulated heating source was adapted to SThM for improved thermal transport measurements.111 The conventional 3ω technique, where the amplitude of the temperature oscillation of a line-shaped heat source mounted on a sample and powered using an AC thermal current (I0sin(ωt)) is measured logarithmically as a function of the angular frequency ω, has been used to determine the thermal conductivities of lead free perovskites.70–72 Meanwhile, the thermal conductivities of hybrid halide perovskite (MAPbX3; X = I, Br, and Cl) single crystals and thin films were first analysed by SThM without the need for extensive sample preparations.61 The thermal conductivity is determined from the slope of the measured amplitudes of the temperature oscillations probe of the thermal probe vs. ln(ω). In order to map local thermal conductivities with high resolution, the resistive thermal probe of the SThM is powered by a mixed signal of angular frequencies ω1 and ω2 (see Fig. 6).50ω1 and ω2 are chosen to be the lower frequency limit ωmin and the upper frequency limit ωmax of the thermal probe. The resulting probe are detected in the frequency domain at 3ωmin and 3ωmax by lock-in amplification for signal to noise improvement by using the output signals of a Wheatstone bridge independently of the ambient probe temperature. The thermal conductivity can be obtained from the difference of these signals using calibrated material references and is particularly suited to assess the thermal transport in thin films, as the interaction volume is ultra-small.112,113 The local thermal diffusivity and thus indirectly the local volumetric heat capacity can be determined simultaneously because the local heat flux into the perovskite sample is given by this experiment. Hence, the thermal transport properties of CsPbCl3 thin films across phase transitions have been analyzed by controlling the sample temperature.51 Since SThM is an AFM based microscopy technique, the thermal expansion coefficient can be measured simultaneously within the same range of temperatures in addition,61 expanding the capability of this instrument. Large negative and positive linear thermal expansion coefficients were found for the tetragonal and pseudo-cubic phases of MAPbI3 single crystals respectively.33,61 Since the resulting thermal stress and structural instability lead to the degradation of perovskites,114 the determination of the thermo-elastic quantities of perovskites depending on the temperature is focused on in the following.
Fig. 6 Schematic illustration of extended quantitative thermal conductivity measurements by SThM applying the so-called 3ω-technique in the frequency domain. Reprinted with permission from ref. 50. Copyright (2019) American Chemical Society. |
The melting temperature of CsPbI3 was determined in 1992 from the phase diagram obtained by differential thermal analysis (DTA) experiments to be higher than 480 °C at atmospheric pressure.115 A short time later, subsequent studies of crystal structures exhibited that CsPbBr3 melts at 567 °C and crystallizes from the melt at 514 °C.116 At the same time it could be demonstrated that melting points increase with decreasing atomic number, i.e. for CsPbCl3 up to 615 °C.117 All these compounds melt congruently as confirmed in these experiments. Hence, large CsPbBr3 crystals could be produced from the melt at high temperatures (600 °C) using the vertical Bridgman method.116 These single crystals have been proposed as promising candidates for X-ray and gamma-ray detection,118–120 since the first investigations on these direct semiconductors showed that most of their properties are comparable to, or even exceed, those of commercial state-of-the-art detector materials.
In this context it must be taken into account that the endothermic DTA peaks for the decomposition pathways of MAPbX3 compounds are reported to be low.121 However, although various decomposition reactions have been proposed:121–124
MAPbX3 → PbX2 + MA + HX | (3) |
MAPbX3 → PbX2 + MAX | (4) |
Fig. 7 Thermogravimetric analyses of MABr, MAPbBr3, PbBr2, CsPbBr3, and CsBr, showing the higher thermal stability of the inorganic perovskites compared to the hybrid organic−inorganic perovskites. Reprinted with permission from ref. 125. Copyright (2016) American Chemical Society. |
With respect to the evaporation or decomposition of CsPbBr3, CsBr is more stable than PbBr2, which is lost first. The small temperature shift between PbBr2 and CsPbBr3 can be explained by the slightly higher stability of the perovskite compound. For hybrid perovskites, it must be considered that all the methylammonium halides have lower melting points (for instance MAI: >190 °C123). Although they are mostly even reported to be over 200 °C121 for MAPbX3, it is still a challenge to have long-term stabilities at 85 °C (representing an elevated temperature on a roof on a hot summer day) required according to the International Standards (IEC 61646 climatic chamber tests) of solar cells.16 Thus the thermal activation energies (Ea) for the decomposition reactions have to be considered.126 From DTAs Eas of 68 ± 2 kJ mol−1, 60 ± 10 kJ mol−1, and 80 ± 20 kJ mol−1 have been reported for MAPbCl3, MAPbBr3, and MAPbI3, respectively.121 Somewhat larger Ea values have been determined by the first mass loss step observed in thermogravimetric analyses.127 Even much higher values have been found for lighter halides in MA and FA based compound calculations,128 summarized in Table 3.
Perovskite | E a [kJ mol−1] | α [×10−4 K−1] |
---|---|---|
MAPbI3 | 93 ± 8127 | α tet.(001): −5.1, αcub.(001): 4.661 |
83128 | α tet.(001): −1.06, αtet.(100): 1.32 | |
α cub.(100): 0.4833 | ||
α tet.(100): 0.58, αcub.(100): 0.3963 | ||
α tet.: 0.90135 | ||
MAPbBr3 | 98 ± 4127 | α cub.(100): 0.33136 |
365128 | α cub.(100): 0.3263 | |
α cub.: 0.67135 | ||
MAPbCl3 | 400128 | α cub.(100): 0.2963 |
CsPbI3 | 650 ± 90127 | α ort.: 0.39, αcub.: 0.39137 |
α ort.(001): 0.22 | ||
α ort.(110): 0.22 | ||
α tet.(001): −0.52 | ||
α cub.(100): 0.40138 | ||
CsPbBr3 | α ort.: 0.38, αtet.: 0.65, αcub.: 0.26139 | |
α ort.(101): 0.38140,141 | ||
CsPbCl3 | α mon.: 0.56, αort.: 0.50, αtet.: 1.27 | |
α cub.: 0.3014 | ||
α cub.: 0.22@T = 85 °C142 | ||
α mon.: 0.42, αort.: 0.38, αtet.: 0.78 | ||
α cub.: 0.29139 | ||
FAPbI3 | 115 ± 3127 | α tet.(100): 0.77 |
108128 | α tet.(010): 0.77 | |
α tet.(001): 0.49143 | ||
α ort.: 0.39135 | ||
FAPbBr3 | 133 ± 1127 | α cub.: 0.89135 |
332128 |
Note that these values do not reflect possible reliability issues related to other causes of degradation.129–131 Details of the chemical, structural, morphological, and optical stability and the impact of thermal stress on the charge transport layers of methylammonium lead halide based perovskites are available elsewhere.114 The main difference between the stability of hybrid and all-inorganic perovskite halides was found to be due to the different chemical properties of the A-site cations, because their entropic contribution to free Gibbs energy is of utmost importance to understand the intrinsic stability of these compounds.132 While a number of stability issues may be solved with proper encapsulation,133 it has to be noted that a critical intrinsic stability problem arises from the behavior under thermal mechanical stress. To make matters worse, it has long been known that perovskites have phase transitions at elevated temperatures up to 85 °C.13,134
The volume thermal expansion coefficients have been studied for CsPbI3 from the temperature dependency of the unit cell volume as (1/V0)·(dV/dT), where (dV/dT) is the change in the volume in the corresponding temperature interval and V0 is the volume at the reference temperature T0.137 The same volume thermal expansion coefficient of 1.18 × 10−4 K−1 has been determined for orthorhombic CsPbI3 (γ-phase) and cubic CsPbI3 (α-phase), with a large discontinuous volume change of about 6.9% at this transition. Assuming CsPbI3 to be isotropic in the first approximation, the volumetric expansion coefficient is three times the linear coefficient and would thus be slightly larger than those of CsPbBr3 and CsPbCl3. Note that it is still a challenge to stabilize the tetragonal phase (β-phase) of CsPbI3 perovskite at ambient temperature. Despite theoretical calculations predicting the β-phase to be more stable than the α-phase,144 experimentally the β-phase is difficult to obtain and to stabilize;145 thus thermal expansion measurements on β-CsPbI3 are extremely rare, to our knowledge.138,146
Large negative thermal expansion coefficients have been observed in the tetragonal phase for yellow and black CsPbI3 (see Fig. 8). This complicates the application in optoelectronic devices because only the CsPbI3 perovskite black phase is optically active. Due to their huge absolute values of thermal expansion coefficients, perovskite films fabricated by existing methods are strained, and the strain is caused by the different coefficients of thermal expansion of the perovskite and substrate during thermal annealing processes as well as due to the growth process itself. It should be noted that films have either compressive or tensile strains in the in-plane direction, which lead to tensile or compressive out-of-plane strains, respectively. Strain has been shown to accelerate the degradation of perovskite films under illumination, which can be explained by increased ion-migration in these films.147
Fig. 8 Structural phase kinetics of strained α-, β-, and γ-CsPbI3 perovskite thin films. From ref. 138. Reprinted with permission from AAAS. |
The strain might significantly influence the determination of thermal expansion coefficients. For example, although negative thermal expansion coefficients have been detected in the tetragonal phase for MAPbI3 perovskites, a five times larger value of αc-tet. = −5.1 × 10−4 K−1 has been found for single crystals61 compared to αc-tet. = −1.06 × 10−4 K−1 determined for thin films.33 Moreover, the thermal expansion coefficients in the cubic phase of MAPbI3 have been reported to be one order larger than those in single crystals. In this context, it must be noted that no discontinuity in thermal expansion has been detected at the phase transition of single crystals in the c-direction. Furthermore, it must be considered that the thermal expansion coefficients depend on the direction in the crystal. αa-tet. is reported to be positive and 1.32 × 10−4 K−1 for these thin films.33 The thermal expansion coefficients determined for the different crystal orientations are summarized in Table 2 in addition.
Note that negative thermal expansion coefficients are generally associated with materials that exhibit ferroelectric, ferromagnetic, or charge-transfer phase transitions that cause the material to contract on increasing the temperature due to either a redistribution or an ordering of ferroelectric/-magnetic domains.148 When discussing the ferroelectricity of perovskites, it must also be considered that as the film thickness decreases, stable domain patterns can change from a ferroelectric–ferroelastic domain to a ferroelectric one, which alters the domain size.149 A challenge for the detection of ferroelectric domains, for instance demonstrated by piezo-response force microscopy on MAPbI3 thin films,150 is that the domain size in thin films can be very small due to their rapidly decreasing thickness-dependence.
However, for the sake of completeness, it should also be mentioned that in comparison to films lower thermal expansion coefficients were determined for single crystalline MAPbI3 in the (100)-direction, αa-tet = 0.58 × 10−4 K−1 and αa-cub. = 0.39 × 10−4 K−1.63 A first order transition with a Δl/l of 2.64 × 10−4 was found in these directional studies. Nevertheless, it could also be demonstrated that the coefficient of thermal expansion of MAPbI3 is larger than those of MAPbBr3 and MAPbCl3 (0.32 × 10−4 K−1 and 0.29 × 10−4 K−1, respectively). In addition, the coefficients of thermal expansion of these hybrid perovskites are larger than those of related Cs–Pb-based all-inorganic perovskites. A larger thermal expansion coefficient was even detected for the lead-free double perovskite (MA)2AgBiBr6 compared to Cs2AgBiBr6 in their cubic phase (0.44 × 10−4 K−1 and 0.28 × 10−4 K−1 in the (111)-direction respectively).151,152 The higher thermal expansion coefficient was attributed to the weaker chemical bonds associated with the softer perovskite (MA)2AgBiBr6.63
Finally, it must be pointed out that for the black phase of γ-CsSnI3, the coefficient of the volumetric thermal expansion at RT was measured to be quite large (1.26 × 10−4 K−1).153 The coefficients of the linear thermal expansion were determined to be αa-ort.: 0.50 × 10−4 K−1, αb-ort.: 0.14 × 10−4 K−1, and αc-ort.: 7.71 × 10−4 K−1. The coefficient of volumetric thermal expansion in β-FAPbI3 (2.06 × 10−4 K−1) was reported to be two times larger than that in the α-phase.143 This ultra-large volumetric thermal expansion creates substantial challenges in device applications. Note that the volumetric thermal expansion is significantly reduced in the cubic α-phase, which corresponds to a negative contribution for large amplitude dynamic PbI6 octahedral tilts. In the α-phase, the linear thermal expansions are 0.77 × 10−4 K−1 in the ab-plane and 0.49 × 10−4 K−1 in the c-direction, suggesting again strong anisotropic mechanical properties. Two wide regions of large negative thermal expansion were found at low temperatures (−100 °C and −219 °C) for δ-FAPbI3.135 The linear thermal expansion for FAPbBr3 was determined to be 0.89 × 10−4 K−1 in the cubic phase at RT from these studies. The symmetry changes are linked again to the tilting of the octahedral (PbBr6) identified by complementary XRD analysis.
Besides that, the Young's modulus E and hardness H values of lead halide perovskite crystals are reported to be low.152,154 Thus, perovskite compounds can be considered as “soft materials”. The elastic moduli of single crystals vary in the range 10–20 GPa and a general trend of ECl > EBr> EI is observed at room temperature. The perovskite compounds can also be regarded as non-brittle, since H varies in the range 0.5–0.25 GPa with the trend HCl< HBr< HI. Metal-halide perovskite films can therefore facilitate applications in flexible devices that require plastic deformation. In this context, however, it should be kept in mind that their E and H strongly depend on the crystalline orientation. The implications of the thermo-mechanical behavior of metal-halide perovskites and interfaces in the reliability of solar cells, as well as opportunities for future research directions in this general area, are highlighted elsewhere.155
The first in-depth analysis of charge-carrier–phonon interactions in hybrid lead halide perovskites was carried out on FAPbI3, FAPbBr3, MAPbI3, and MAPbBr3.156 For FA-based perovskites it was already found that the Fröhlich interaction between charge carriers and LO phonons is the dominant contribution to the predominantly homogeneous broadening of the PL line width in these hybrid perovskites at room temperature (see Fig. 9). The Fröhlich interaction for bromide perovskites is larger than that for iodide perovskites, providing a link between composition and electron–phonon scattering that fundamentally limits the mobilities of the charge carriers. In addition to these spectral changes, the lattice strain causes non-radiative losses in halide perovskites,157 which has in addition been attributed to the locally induced crystal misorientation.158 The remaining nanoscale strain also reduces the stability in cesium-lead-bromide perovskites and shifts the local luminescence spectrum.159 In contradiction to these studies, it has been reported previously that light-induced lattice expansion leads to high-efficiency perovskite solar cells.160 While the authors of these studies claim to exclude heat-induced lattice expansion, it has been shown that the lattice expansion mechanism matches the heat-induced thermal expansion during illumination.161
Fig. 9 Temperature dependence of normalized steady-state PL on (a) FAPbI3, (b) FAPbBr3, (c) MAPbI3, and (d) MAPbBr3. From ref. 156. Reprinted with permission from Springer Nature. |
In summary, the information provided so far can be used to predict changes in the optical and charge transport properties of devices, since any structural change has a serious impact on the phononic and electronic structures.162 The properties of phonons are crucial for the understanding of not only charge transport in the form of electron–phonon interactions in polar semiconductors such as halide perovskites, but also for heat transport in functional materials. Phonon–phonon scattering plays a decisive role in determining acoustic phonon transport and thus the thermal conductivity of halide perovskites. Details about these coupling mechanisms have recently been summarized in the review article of Haque et al.25 In addition, in terms of the soft nature of halide perovskites, temperature and pressure are expected to have a strong influence on the crystal growth and the boundary fusion of these materials.
Wang et al. were able to show first that, similar to other “soft materials”, MAPbI3 can directly be patterned by thermal nanoimprint lithography (NIL) at low temperatures (100 °C), since electrostatic interactions dominate in ionic MAPbI3 crystals (see Fig. 10a).165 In these studies, nanoimprinted photodetectors were demonstrated with significantly improved performance compared to nanoimprinted devices. The performance enhancement was attributed to the NIL-induced improved crystallinity; particularly, the nanograting structure was favorable for better photon absorption and charge carrier transport. NIL has also been employed for the fabrication of meta-surfaces.166,167 By using a flat stamp in NIL, it could be shown that this imprinting process is able to substantially flatten an initially very rough polycrystalline hybrid perovskite layer and to transform it to a thin film consisting of laterally extended crystals on the order of tens of microns with a surface roughness below 1 nm (rms) (see Fig. 10b).61,168,169 A notable effect of this planar hot-pressing procedure was that these grains have similar physical properties to single crystals, e.g., the same thermal conductivity and PL-lifetime. Besides that, the films are well oriented and show enhanced stability when exposed to ambient air.169 The latter evidences earlier observations of enhanced stability of larger crystal grains.171 Recently, the planar NIL fabrication of continuous MAPbI3 quasi-single-crystal films composed of even 100 μm large-sized grains with deep boundary fusion has been reported.172 Here, a remarkable decrease of trap density (ntrap: 7.43 × 1011 cm−3) determines a long carrier lifetime (τ: 1.7 μs). Thus, planar hot-pressing has a lot of advantages in comparison to other procedures to produce large crystals.173–175 The ability to form nanostructures is particularly important for the preparation of resonator structures for metal-halide perovskite lasers.176 Hence, photonic nanostructures were patterned by thermal NIL directly into MAPbI3 perovskite films using nanopillar molds (illustrated in Fig. 10c).170 Apparently, surface defects which cause substantial recombination and luminescence quenching were eliminated by recrystallization upon thermal imprinting. The resulting 2D photonic crystal (2D-PC) lasers showed low lasing thresholds of 3.8 μJ cm−2. The first distributed feedback (DFB) lasers based on MAPbBr3 thin films with a linear second order photonic grating (periodicity: 300 nm) directly imprinted into the active layers were also presented.169 Very low lasing thresholds of 3.4 μJ cm−2 were achieved in these optically pumped DFB laser structures. Using a flat stamp to recrystallize the layer via NIL, it was shown that a significantly lowered threshold for the onset of amplified spontaneous emission (ASE) compared to pristine films was achieved due to reduced optical scattering. Note that low lasing threshold levels are essential on the way to the first electrically operated perovskite laser diode. A scalable lithographic scheme for thin organic metal-halide perovskite films has already been demonstrated and utilized for the realization of perovskite micro lasers.177
Fig. 10 (a) Schematic of a NIL process with a Si flat, nanopillar, and nanograting mold. Reprinted with permission from ref. 165. Copyright (2016) American Chemical Society. AFM images of (b) pristine and (c) planar hot-pressed MAPbBr3 layers. From ref. 169. Copyright Wiley-VCH Verlag GmbH & Co. KGaA. Reproduced with permission. (d) SEM image of the 2D-PC with dimensions of the triangular grating. Inset (bottom): Photograph of the sample under oblique while-light illumination. Inset (top): First Brillouin zone of the 2D-PC with labelled points of high symmetry. (e) Emission spectra upon optical pumping with increasing energy density (1.6–4.7 μJ cm−2). A scaled PL spectrum of a pristine MAPbI3 layer upon cw-excitation is shown for comparison. From ref. 170. Copyright Wiley-VCH Verlag GmbH & Co. KGaA. Reproduced with permission. |
It is noteworthy that all-inorganic metal halide perovskite films based on Cs can likewise be structured and flattened using NIL at temperatures <150 °C.50,51,178 The first CsPbBr3 thin-film DFB and vertical cavity surface emitting lasers (VCSELs) with an ultralow threshold (2.2 μJ cm−2, 300 ps pump pulses) operating at room temperature, which do not rely on the use of nanoparticles, were thus realized by NIL.178 In addition, the planar hot-pressing procedure enabled the creation of 0D, 2D, and 3D perovskite films depending on the stoichiometry of the pristine layer.50,51
In conclusion, these results indicate that NIL will be an excellent thermo-mechanically supported process platform for the patterning of metal-halide perovskites in the future.51 Thermal NIL can improve the crystal quality as well as create the required patterns. However, although there are initial approaches to describe the crystal quality improvement and the growth processes by NIL, a detailed simulation model of this hot-pressing method for perovskites has not been developed yet.
In addition, it must be noted that a strong temperature-dependence of stimulated emission has been observed in MAPbI3 thin films.179 A heating of up to ≈50 °C has been detected at RT in the first few hundred ns following pump turn-on, which leads to a dynamic increase in threshold and limits continuous-wave lasing in hybrid perovskite semiconductors. Furthermore, Joule heating prevents hybrid LEDs from reaching current densities greater than a few hundred A cm−2 without total emission quenching.180 Therefore, an improvement in heat dissipation and a reduction in the pulse width are inevitable for proper operation. The first effective thermal management strategies were proposed by L. Zhao et al., which enables bright and stable perovskite LEDs.181
Concepts for enhancing grain sizes, improving crystallinity, changing elastic properties and passivation of perovskite films by using organic additives were introduced.183,184 In these studies, the inclusion of various additives in MAPbI3 was examined and a significant improvement in mechanical robustness was observed. Besides, a dimensionality-controlled surface passivation for enhancing the performance and stability of perovskite solar cells has been pointed out recently.185 Note that in metal-halide perovskites the thermal properties of 0D and 2D materials significantly differ from those of the bulk (3D) and are thus discussed in the following section.
In order to combine charge transport and chemical protection in one, great efforts have also been made to produce functional layers with optimized optical, electrical and permeation barrier properties.186 Brinkmann et al. demonstrated that introducing a bi-layered electron-extraction interlayer consisting of aluminium-doped zinc oxide nano-particles and tin oxide, grown by low temperature atomic layer deposition (ALD), can significantly mitigate the decomposition of the perovskite in inverted perovskite solar cells even at elevated temperatures.187,188 It should be noted that in this example the ALD-grown tin oxide formed an outstandingly dense permeation barrier.
The reduction of stress and strain in films is also an important issue that should be addressed for the reliability of perovskite devices. Unfortunately, only a few studies have been carried out so far. Rolston et al. studied mechanical stresses that develop during the fabrication of perovskite films and a direct relationship between film stress and stability in both perovskite films and solar cells was determined.189 They developed simple strategies to reduce the stress and to improve the inherent moisture and thermal stability of perovskite films, e.g., by using low process temperatures, using polymeric substrates with thermal expansion coefficients that are closer to those of the perovskites, and pre-straining the substrate during perovskite processing. Otherwise highly textured morphologies are created by popular anti-solvent conversion methods due to in-plane compressive stress.190 The energy release leads to wrinkling during the intermediate phase of film formation, in which the substrate constrains the film from expanding. Ultra-smooth films can be obtained by slowing the rate of film formation since the extent of wrinkling is correlated with the rate of film conversion.
The regulation of strain in perovskite thin films by suitable selection of charge transport layers (CTLs) was also reported (the energy level alignment is shown in Fig. 11).191 Tensile strains in perovskite films were compensated for by transport layers with high thermal expansion coefficients and by elevating the processing temperature.
Fig. 11 Energy level alignment of various materials used as ETLs and HTLs in perovskite devices. From ref. 192. Reprinted with permission from Elsevier. |
Note that thermal effects can also impact other functional layers in perovskite devices, such as CTLs and electrodes. The extensive studies on their thermally enhanced chemical decompositions and strategies for their stabilization are reviewed in the article by Ava et al.114 For example, metal electrodes, which might contribute to thermal management, can be degraded by the corrosive decomposition products of the perovskite. These layers must therefore also be protected. This is even more problematic for semi-transparent solar cells, in which the semi-transparent electrodes are based on ultrathin metal films. Thus, a self-encapsulating thermostable and air-resilient electrically conductive permeation barrier was realized by an ultrathin Ag layer sandwiched between SnOx grown by low-temperature atomic layer deposition.193 These additional layers further complicate the design of components in terms of mechanical stability. The efforts to increase the energy of fracture of perovskite solar cells were compiled by Boyd et al.194
Although recent review articles discuss factors affecting the chemical and mechanical stability of perovskite solar cells as well as current issues and new strategies to push material and device stability,195,196 strategies for improved heat management in perovskite devices have not yet been engineered. It must be ensured that the other important layers within the devices can dissipate the heat because perovskites have a very low thermal conductivity. For example, in addition to a low coefficient of thermal expansion of 5 × 10−6 K−1, glass has a very low thermal conductivity of approximately 1 W (m K)−1 and a very low specific heat capacity in the range of 0.8 J (g K)−1, which prevents heat dissipation into this substrate.197 The commonly used transparent electrode material indium tin oxide (ITO) is degenerately doped (carrier density >1 × 1020 cm−3) and behaves like a metal.198 Nevertheless, the resulting in-plane heat transport in the ITO layer is strongly limited by its low λITO = 6 W (m K)−1 and a typically small film thickness of around 150 nm. In addition, it must also be noted that organic-based ETLs and HTLs, which are depicted in Fig. 11, have likewise ultra-low thermal conductivities (λPCBM = 0.031 W (m K)−1), and are highly anisotropic for the in-plane and out-of-plane directions (λPEDOT:PSS‖ = 1 W (m K)−1 and λPEDOT:PSS⊥ = 0.3 W (m K)−1, respectively).199–201 Thus, instead of organic-based charge transport layers, metal-oxide-based compounds, which have significantly higher thermal conductivities, could be used for improved thermal management. Polycrystalline Al-doped ZnO films, for instance, have more than one order higher thermal conductivities in the range of λZnO = 6 W (m K)−1 at RT and are therefore ideal for additional heat dissipation in ETLs.202 The same applies to TiO2 layers with λTiO2‖ = 12.6 W (m K)−1 and λTiO2⊥ = 8.8 W (m K)−1.203 Metal-oxides have also evolved as a powerful class of materials on the anode side.204 Transition metal oxides like MoO3, WO3, and V2O5 have successfully been introduced in organic and perovskite solar cells, which should enable enhanced thermal transports in HTLs.205 Charge transport layers based on graphene or silver nanowires are further potential candidates in heat management and have also been studied in recent years.206,207 Note that, in general, many strategies that have been proposed for the thermal management in organic devices can be used in the field of perovskites.68
In addition to the absolute values of thermal conductivities for the respective materials, it has to be considered that heat transport in the perovskites is due to phonons, while in the electrodes it is mostly associated with electrons. Therefore, it must be ensured that the heat transport at this interface is not limited by its thermal impedance. In order to guarantee an optimized heat flux, an enhancement of the thermal boundary conductance of such a system has recently been demonstrated by nm sized additional adhesion layers, either titanium or nickel.208 In this context it should likewise be mentioned that very thin layers can be considered as thermally transparent depending on their mean-free phonon and electron path lengths.112,113 This means that additional thin functional layers that do not impair the heat flow can be used in the design of metal-halide devices.
Although the effectiveness of heat spreading concepts is often directly related to the thermal conductivity of the materials used, alternative passive means of cooling should be considered, as well.209 One is to use CTLs with an elevated specific heat, which can absorb the heat generated during dynamic device operation. A different procedure is to effectively store heat dissipated from the electronic components, which is periodic in nature or suddenly temporary, using high latent heat materials (see the Thermal properties across phase transitions section). The use of materials with solid and liquid organic phase changes has previously been proposed with regard to the thermal management of mobile phones210 and the practical applications of liquid-phase-exfoliated graphene as filler material in phase change materials have been outlined.211 A review of the thermal properties of graphene, few-layer graphene and graphene nanoribbons has also been provided therein along with a discussion of practical applications of graphene in thermal management and energy storage.
An ultralow thermal conductivity of 0.15 W (m K)−1 was discovered in 0D hybrid lead-free [Mn(C2H6OS)6]I4 single crystals at room temperature; it was further decreased by introducing disorder into the system.214 Furthermore, the thermal conductivities of 2D hybrid perovskites could be varied in a range from 0.10 to 0.19 W (m K)−1 depending on the molecular substructure of ammonium cations (details of the thermal properties can be obtained from Table 1).32 The extremely low thermal conductivities can be attributed to the crystalline substructures of these perovskites. In a similar sense, substantially reduced thermal conductivities have been discussed for layered, 2D complex thermoelectric materials.216 This effect has been explained by increasingly glass-like properties, which hamper phonon propagation. Moreover, the low thermal conductivity found in 0D perovskites has earlier been associated with their soft phonon dispersion.77 It was found that the thermal conductivities of these perovskite films, due to the anisotropy of the crystalline structure, are significantly dictated by their relative orientation (see Fig. 12). The results presented so far suggest that organic based low-dimensional perovskites provide a unique opportunity to achieve ultralow thermal conductivity, which would be beneficial for thermo-electric applications.
Fig. 12 Schematic illustrations of vertically and randomly oriented BA2MA3Pb4I13 crystal orientations that lead to different thermal conductivities. Reprinted with permission from ref. 32. Copyright (2020) American Chemical Society. |
For all-inorganic perovskites, such as caesium lead chloride, it could be demonstrated for 3D phase CsPbCl3, 2D phase CsPb2Cl5, and 0D phase Cs4PbCl6 that the thermal conductivities as well as volumetric heat capacities become smaller upon reducing the dimensionality (see Table 4).51,215 Note that the thermal diffusivities, which are responsible for the dynamic heat transfer in devices, remain almost constant within the measurement uncertainties. Thus, the presented results offer a promising general route for designing 0D and 2D materials for future materials with exceptional thermal properties.
In another example, 2D phase CsPb2Br5 films show a low thermal conductivity of 0.33 W (m K)−1, a low thermal diffusivity of 0.3 mm2 s−1, and a small heat capacity of 0.18 J (g K)−1, which are drastically lower than those of their 3D counterparts (shown in Table 1).50
Fig. 13 Schematic representation of the phase transitions of hybrid and all-inorganic perovskites.44–47,135,218 |
Briefly summarized, first and second order transitions from a monoclinic to an orthorhombic to a tetragonal structure at 37 °C and 42 °C, respectively, and a first order transition to a cubic structure at 47 °C have been reported for CsPbCl3 single crystals.44,45 The measurement showed that CsPbCl3 is neither ferroelectric nor antiferroelectric in any of the phases.44,46 The phase transitions are basically associated with the condensation of rotational modes of PbCl6 octahedra around the three principal axes.47 A considerable acoustic attenuation was observed in these experiments in the orthorhombic phase. The specific heat csp and thermal diffusivity a values of the different phases as well at their transitions were obtained in transient measurements on single crystals.11 Similarly, two characteristic phase transitions have been found for CsPbBr3, the orthorhombic/tetragonal at 82 °C and the tetragonal/cubic at 132 °C for temperatures above RT.217 Significant differences were reported out of these macroscopic investigations for the heat transitions obtained between powder samples and bulk samples, which was assumed to be due to many defects in the powder particles. Furthermore, the pulse transient method used in these studies did not evidence the expected behavior of csp(T) when going from the monoclinic to the orthorhombic to the tetragonal phase of CsPbCl3.11 Specifically, at a first order transition, one would expect a divergence of csp(T) due to the latent heat involved. Note that it is surprising that all these transitions occur within a temperature range of ΔT ≈ 10 °C in the case of CsPbCl3 measurement. More exceptionally, for CsPbBr3 the two phase transitions (orthorhombic–tetragonal–cubic) are separated by a ΔT = 50 °C.217 However, even though clear changes in csp(T) could be observed, this behavior did not agree with theoretical predictions of an abrupt change in csp(T).
The variation of the thermal properties across the phase transitions in the range from RT to 65 °C was recently analyzed for CsPbCl3 thin films by SThM. A significant increase of the thermal conductivity was found for the cubic phase (T > 46 °C), while only slight variations in λ() were detectable for transitions from the monoclinic to the orthorhombic to the tetragonal phase. It was demonstrated that upon transition to the cubic phase, the relative increase of average cvol.() is significantly higher than that of λ(), because a() and cvol.() measurements are extremely sensitive to the amount of heat ΔQ() involved in the respective transition. Thus, the thermal diffusivity in the cubic phase is notably lower in comparison to that of the respective phase at room temperature.
Although phase transition studies using expansion measurements on all-inorganic perovskites were carried out at an early stage, tracking structural phase transitions in the corresponding hybrid lead-halide perovskites have only recently been carried out.135 On the other hand, the influence of phase transitions on the specific heat capacity of MAPbX3 (X = I, Br, and Cl) perovskites was already measured in 1990.52 A tetragonal to pseudo-cubic polymorph transition at 57 °C was observed for MAPbI3, while the orthorhombic to tetragonal phase change and the transitions of all other MAPbX3 perovskites occurred at substantially lower temperatures (−124 °C to −37 °C). In addition, anomalies in the anticipated thermo-physical properties at 57 °C were identified. Although all of these phase transitions were expected to be of the first order, the specific heat capacity characteristics in MAPbBr3 and MAPbI3 were found to rather indicate a second order transition.
The thermal conductivity characteristics of MAPbX3 perovskite single crystals were also evaluated above room temperature.59,61,219 Only slight variations of around 0.30–0.42 W (m K)−1 were observed in the temperature range of RT to 150 °C for MAPbI3.59 A similar behavior was found by using a laser flash technique for the indirect determination of thermal conductivity.219 The thermal diffusivity would thus decrease during the phase transition at 58.9 °C only due to latent heat. However, these experimental results contradict previous theoretical studies that predicted a significant increase of the thermal conductivity due to the higher group velocities of acoustic phonons in the cubic structure.31,35
As opposed to these indirect methods, it was demonstrated that the thermal conductivity of MAPbI3 single crystals in the cubic phase increased to 1.1 W (m K)−1 from 0.34 W (m K)−1 in the tetragonal phase by applying a direct approach to measure the thermal conductivity, i.e. the so called 3ω-technique.61 The measured thermal conductivities were in good agreement with the previously reported equilibrium molecular dynamics simulations.35 The phase transition at a somewhat lower temperature of 52 °C and the negative coefficient of thermal expansion found in the tetragonal phase of MAPbI3 were in agreement with those obtained from X-ray diffraction.33 The absolute values of the negative linear thermal expansion coefficients in the tetragonal structure were found to be quite similar to the positive linear thermal expansion coefficients in the cubic structure. In this context, however, it should be noted that, unlike single-crystal organic–inorganic hybrid perovskite expansion analyses,219 this transition showed no discontinuity in the Δl/l characteristics. The temperature-dependent thermal conductivity characteristics of the cleaved MAPbI3 single crystal were also reported to be continuous, which rather points to a second order transition, similar to the data of the specific heat capacity measurements.52
In view of the particular thermal properties of this class of materials, such as their notoriously low thermal conductivities, low specific heat capacities, and high thermal expansion coefficients, it was demonstrated that proper thermal management in devices will be key to overcoming some current limitations and may even unlock novel applications. Thus, advanced comprehensive tools and strategies for systematic thermal analyses and property optimizations were discussed. Furthermore, the impact of the perovskite dimensionality (3D, 2D, 0D) on the thermal properties and how these properties change across the various phase transitions of these perovskites were highlighted.
This journal is © The Royal Society of Chemistry 2020 |