Rationalizing the light-induced phase separation of mixed halide organic-inorganic perovskites.

Mixed halide hybrid perovskites, CH3NH3Pb(I1-x Br x )3, represent good candidates for low-cost, high efficiency photovoltaic, and light-emitting devices. Their band gaps can be tuned from 1.6 to 2.3 eV, by changing the halide anion identity. Unfortunately, mixed halide perovskites undergo phase separation under illumination. This leads to iodide- and bromide-rich domains along with corresponding changes to the material's optical/electrical response. Here, using combined spectroscopic measurements and theoretical modeling, we quantitatively rationalize all microscopic processes that occur during phase separation. Our model suggests that the driving force behind phase separation is the bandgap reduction of iodide-rich phases. It additionally explains observed non-linear intensity dependencies, as well as self-limited growth of iodide-rich domains. Most importantly, our model reveals that mixed halide perovskites can be stabilized against phase separation by deliberately engineering carrier diffusion lengths and injected carrier densities.Mixed halide hybrid perovskites possess tunable band gaps, however, under illumination they undergo phase separation. Using spectroscopic measurements and theoretical modelling, Draguta and Sharia et al. quantitatively rationalize the microscopic processes that occur during phase separation.

Introduction

Despite dramatic efficiency (η) increases of mixed organic–inorganic lead halide perovskite solar cells from 5 to 22%[1-4], η-values still lie below the Shockley limit[5, 6]. Ongoing efforts have therefore targeted bandgap tuning as a means of improving photovoltaic performance by better capturing the full solar spectrum[7]. Mixed halide perovskites (e.g., MAPb(I1−Br)3, where 0 < x < 1 and MA is the methylammonium cation (CH3NH3 +)) represent an appealing approach to realizing efficient bandgap tuning. Corresponding bandgaps are linear functions of mixing fraction x and span the range 1.57–2.27 eV[8]. Unfortunately, MAPb(I1−Br)3 exhibits unwanted instabilities under visible light irradiation, leading to I and Br segregation into separate iodide- and bromide-rich domains. While general aspects of photoinduced MAPb(I1−Br)3 phase separation have been established[9-16], unresolved issues remain. They include an explanation of dramatic spectral differences in the emission and absorption of phase-separated perovskites, the excitation intensity (I exc)-dependent growth of the iodide-rich domain emission, and non-linear intensity dependencies of the phase separation rate constant. Establishing a deeper understanding of light-induced MAPb(I1−Br)3 phase separation is thus key to engineering and controlling the stability of mixed halide perovskites for photovoltaic as well as other optoelectronic applications[17-19].

Here we rationalize the bulk of experimental observations regarding the light-induced phase segregation of MAPb(I1−Br)3. From the first principles density functional theory (DFT) calculations, we explain the driving forces leading to phase segregation. We then estimate the initial size of phase-separated domains and develop an analytical model that accounts for observed non-linear intensity-dependencies in kinetic rate constants. Most importantly, we suggest how light-induced phase segregation can be mitigated in MAPb(I1−Br)3 thin films.

Results

Spectroscopic observations

We first introduce spectroscopic observations that characterize the prolonged illumination of MAPb(I1−Br)3 thin films. In particular, Fig. 1a shows the evolution of a MAPb(I0.5Br0.5)3 thin film emission spectrum under continuous optical irradiation (λ exc = 405 nm, I exc = 20 mW cm−2). Evident is a decrease in the native MAPb(I0.5Br0.5)3 emission at λ mix = 652 nm accompanied by a corresponding rise of emission features at 725 and 527 nm. The former (latter) is associated with MAPbI3-like (MAPbBr3-like) photoluminescence. In what follows, the text therefore refers to the red (green) spectral feature as emission from iodide-rich (bromide-rich) domains. Furthermore, the emission growth kinetics of iodide- and bromide-rich domains are identical, as shown in Supplementary Fig. 1. Analogous behavior has been seen in films with other x-values (Supplementary Fig. 2).

Fig. 1

Spectroscopic observations of halide phase segregation. a Time evolution of MAPb(I0.5Br0.5)3 (MA = CH3NH3 +) emission spectra (PL, photoluminescence) under 405 nm continuous wave excitation (I exc = 20 mW cm−2). Times for selected spectra (from red to purple): 0.05, 1.41. 1.64, 1.69, 1.83, 1.93, 2.26, 2.40, 2.58, 2.68, 2.87, and 3.10 s. Inset: Emission spectra between 475 and 600 nm. b Corresponding time evolution of the MAPb(I0.5Br0.5)3 absorption spectra under 405 nm continuous wave (CW) excitation (I exc = 25 mW cm−2). Times for selected absorption spectra (from blue to green): 0, 1, 30 min. Inset: Absorption-based phase separation kinetics from absorption changes at 720 nm where the dashed red line represents an exponential fit to the data. c I iodide /I sat (I iodide – iodide-rich emission intensity, I sat – saturation emission intensity) under different I exc (I exc – excitation intensity). Dashed lines are fits using Eq. (4). The bottom solid red line shows I iodide /I sat when I exc = 40 μW cm−2. Excitation intensities for selected curves (from red to purple): 0.27, 0.79, 1.05, 1.46, 1.56, 2.27, 3.81, 5.04, 18.56, 57.28 mW cm−2. d I exc-dependent emission-based first order rate constant for phase separation. The dashed red line represents a fit to the data using Eq. (5)

These dramatic changes in emission are not evident in absorption. To illustrate, Fig. 1b shows the response of the MAPb(I0.5Br0.5)3 spectral absorption to irradiation. A slight decrease of the band-edge extinction is accompanied by growth of a tail at longer wavelengths.

Various attempts have been made to explain these photoinduced changes to the absorption/emission of MAPb(I1−Br)3 films[10, 11]. In brief, prior DFT calculations have suggested that photoirradiation causes lattice expansion, leading to MAPb(I1−Br)3 bandgap reduction[20]. MAPbI3 and MAPbBr3, however, are stable upon illumination. This suggests that a more universal explanation for the phenomenon exists. Photoinduced defect formation has also been invoked to explain subgap absorption/emission in MAPb(I0.4Br0.6)3 films[11]. However, the exact nature of these metastable states has never been explained.

Today, it is generally accepted that photoinduced changes to the absorption/emission of MAPb(I1−Br)3 films originates from light-induced halide phase segregation. This has been established by McGehee and co-workers using thin film X-ray diffraction measurements carried out during MAPb(I1−Br)3 illumination[10]. These measurements reveal the formation of two types of crystalline domains having lattice parameters consistent with those from iodide- (i.e., MAPb(I0.8Br0.2)3) and bromide-rich (i.e., MAPb(I0.3Br0.7)3) perovskite films with estimated conversion fractions ranging from ϕ = 20 to 23%[10]. Complementary absorption measurements suggest ϕ = 18%[12]. Temperature-dependent emission measurements have additionally yielded phase separation activation energies in good agreement with reported halide anion migration energies for inorganic perovskites (0.25–0.39 eV)[10], corroborating the above X-ray analysis.

Evidence for phase separation also comes from recent transient differential absorption measurements [12]. Namely, under external illumination a ground state bleach maximum associated with the absorption edge of MAPb(I0.5Br0.5)3 disappears and is replaced by two new bleach features, one at 530 nm and another at 720 nm. The former (latter) is associated with the absorption edge of bromide-rich (iodide-rich) perovskite films. These transient absorption results notably agree with the emission data in Fig. 1a.

Thus, partial phase separation of MAPb(I1−Br)3 thin films leads to near exclusive emission from iodide-rich domains. By contrast, corresponding absorption spectra remain nearly unchanged. Beyond this, little is known about why photoirradiation leads to MAPb(I1−Br)3 phase separation, let alone why complex phase separation kinetics are observed.

Figure 1c shows that MAPb(I1−Br)3 phase separation kinetics are I exc-dependent. Using emission data from an illuminated MAPb(I0.5Br0.5)3 thin film, Fig. 1c plots I iodide /I sat as a function of time; I iodide is the integrated intensity of the iodide-rich domain emission and I sat is an empirical saturation intensity. It is apparent that I iodide grows exponentially. Subsequent fitting of the data to yields an apparent, emission-based, first order rate constant for phase separation (k forward,em). In addition, Fig. 1c shows that below an excitation intensity threshold of I exc = 40 μW cm−2 no phase separation occurs.

Figure 1c thus makes three key points. First, k forward,em is excitation intensity-dependent. Namely, the rate of phase segregation increases with increasing I exc, as seen through a correspondingly faster exponential growth of I iodide /I sat. Corresponding phase separation rate constants, shown in Fig. 1d, range from k forward,em = 0.092–0.823 s−1 for I exc = 0.26–58.2 mW cm−2. Second, k forward,em saturates at k forward,em = 0.707 s−1 when I exc exceeds a critical excitation intensity of 20 mW cm−2. For comparison purposes, the associated absorption-based forward rate constant at I exc = 25 mW cm−2 is k forward,abs = 0.038 s−1 (inset, Fig. 1b). Finally, there exists an excitation intensity threshold below which no phase separation occurs.

Beyond this, there is a large asymmetry in MAPb(I1−Br)3 phase separation and recovery kinetics. In this regard, it is known that, under dark conditions, phase separation is reversible (confirmed by X-ray measurements[10]), leading to full recovery of the original MAPb(I1−Br)3 absorption and emission spectra over the course of 5 min [10-12]. The asymmetry is attributed to large activation energies for halide ion migration in the dark[10, 12, 15]. In the current study, we find an emission-based (absorption-based) recovery rate constant of k reverse,em = 6.8 × 10−3 s−1 (k reverse,abs = 2.3 × 10−3 s−1) by analyzing temporal kinetics of MAPb(I0.5Br0.5)3 emission (absorption) recovery (Supplementary Fig. 3). These values agree with existing literature values on the order of 10−3 s−1 [12] and reveal that k reverse is easily two orders of magnitude smaller than k forward.

A complete microscopic model that accounts for all these observations, specifically, the absorption/emission spectral asymmetry, the I exc-dependent forward rate constant, its saturation above a critical I exc-value, and the existence of a threshold intensity below which no phase separation occur, remains elusive. Towards addressing these open questions, Ginsberg and co-workers have recently proposed a microscopic model that explains some of these observations[15]. In it, photogenerated polarons localize in iodide-rich regions of MAPb(I1−Br)3 thin films, where pre-existing iodide-rich sites originate from stochastic composition fluctuations. Localized polarons induce local lattice strain, which, in turn, promotes further iodide migration to their immediate vicinity. This leads to iodide-rich domain growth. Self-limited growth arises from the finite size of the polaron deformation region and yields 8–10 nm iodide-rich domains. Subsequent emission from phase-separated films predominantly occurs from iodide-rich regions of the film due to their smaller bandgap relative to MAPb(I1−Br)3.

Despite successes of the polaron model in rationalizing phase segregation and self-limited domain size growth, many ancillary observations, including spectral asymmetries in absorption/emission response, as well as non-linear I exc-dependent kinetics, have not been addressed. Consequently, a more complete accounting of light-induced MAPb(I1−Br)3 phase segregation is needed to identify conditions that can increase the stability of MAPb(I1−Br)3 films under illumination.

Phase segregation model

In what follows, we describe a conceptual framework that rationalizes virtually all current experimental observations regarding MAPb(I1−Br)3 phase separation. The framework consists of a kinetic model that describes events following MAPb(I1−Br)3 photoexcitation. Implicit assumptions of the kinetic model are supported by first principles DFT calculations. The kinetic model is subsequently accompanied by analytical (probabilistic) estimates to describe the emission from phase-separated MAPb(I1−Br)3 films.

Figure 2a illustrates microscopic processes that underpin the kinetic model. First, light absorption creates an electron–hole pair. Along Path 1, these carriers recombine to produce native MAPb(I1−Br)3 emission (λ mix = 652 nm, x = 0.5). Along Path 2, photoexcitation provides the driving force to induce local halide anion rearrangement and microscopic MAPb(I1−Br)3 phase segregation into iodide- and bromide-rich domains. Although the exact final composition of these iodide-rich domains is not known, Hoke et al. have suggested the existence of a stable iodide-rich phase, following prolonged illumination, with a nominal composition of x~0.2[10]. Subsequent emission occurs from these iodide-rich domains (Path 2a, λ iodide) due to their favorable band offsets. In tandem, the phase separation reverses on significantly slower timescales, along Path 2b[12]. Finally, Path 3 considers the possibility of MAPb(I1−Br)3 photogenerated carriers diffusing to existing iodide-rich domains on length scales comparable to electron/hole diffusion lengths (ca. l e/h = 100 nm)[21]. Consequently, emission can occur from iodide-rich regions of phase-segregated films even when photoexcitation events themselves do not initiate phase segregation. Supplementary Note 1 provides a complete mathematical accounting of the model.

Fig. 2

Schematic of relevant kinetic processes during halide phase segregation. a low and b high excitation intensity light-induced phase separation of MAPb(I1−Br)3 (MA = CH3NH3 +). Microscopic rate constants (k 1−3) associated with the underlying kinetic model have been provided. Dark circles denote photogenerated electron–hole pairs. Empty circles denote electron–hole pairs which have induced phase separation. Filled red and white regions represent phase separated iodide-rich and corresponding bromide-rich domains

DFT calculations detailed at length in Supplementary Note 2 (see Supplementary Fig. 4) support major assumptions underlying the kinetic model. Specifically, calculations in four formula unit (FU) supercells show that the energy (ΔE GS(x)) to form mixed MAPb(I1−Br)3 from MAPbI3 and MAPbBr3 is small and slightly positive (on the order of 0.01 eV FU−1) up to x = 1/3 and modestly negative at larger x (Fig. 3a, blue circles). In this limit of nearly zero formation energies, mixing entropy controls phase stability. At ambient temperature, the ideal entropy of halide mixing substantially exceeds the mixing energy, so that, as illustrated in Fig. 3a (green squares), the mixed phase minimizes free energy (ΔF GS(x)) over a substantial composition range. Brivio et al.[9] report DFT-computed mixing energies of similar magnitude and nearly ideal mixing entropies over the entire composition range. Both results therefore imply that over the composition range considered, entropic considerations lead to MAPb(I1−Br)3 formation.

Fig. 3

Density functional theory modeling of halide phase segregation. a Estimated free energy of formation as a function of bromide fraction x. The blue line/circles are generalized gradient approximation (GGA)-computed 0 K ground state formation energies, including spin-orbit coupling. The green line/squares are free energies of mixing at 300 K, assuming ideal mixing on the halide sublattice. The red line/inverted triangles represent the free energy difference after single photon absorption. All energies are reported per two formula units. Blue and pink shaded areas indicate positive (phase separation favored) and negative (phase separation disfavored) free energy regions. b Approximate size of separated domains, n max, following single photon absorption

Next, we use DFT calculations to estimate the ground state valence band position of MAPb(I0.5Br0.5)3 relative to those of MAPbI3 and MAPbBr3. This comparison is shown in Supplementary Fig. 5, which reveals lower (higher) valence band edges within iodide-rich (bromide-rich) domains. Further, these calculations and available experiment[22, 23] indicate that mixed and iodide-rich conduction band edges are nearly isoenergetic. Altogether, these results imply that photogenerated holes will preferentially localize in iodide-rich regions of the film and will ultimately contribute to red-shifted, MAPbI3-like emission (Path 2, Fig. 2a).

This valence edge difference provides an energy driving force for demixing that counterbalances the above entropic preference for halide mixing. We can estimate this trade-off by expressing the phase separation energy between a photoexcited hole in MAPb(I1−Br)3 and a hole localized in an iodide-rich phase. Because of their near isoenergetic conduction band edges, experimentally-observed MAPb(I1−Br)3/MAPbI3 band gap differences[8] which increase with x, capture the demixing driving force. Consequently, the free energy change associated with the transformation of n formula units of photoexcited MAPb(I1−Br)3 into pure bromide and pure photoexcited iodide (ΔF *) can be written aswhere , F Br, and denote photoexcited MAPb(I1−Br)3, MAPbBr3, and MAPbI3 free energies and asterisks denote the presence of a photogenerated hole. ΔF is plotted in Fig. 3a (inverted triangles) for n = 2 at T = 300 K. Note that separation into iodide- and bromide-enriched phases yield correspondingly smaller, yet positive, ΔF GS values. In this regard, Fig. 3a shows ΔF to be close to zero in the composition region x < 0.1 with a predominantly iodide-rich composition of MAPb(I0.9Br0.1)3. Although the estimate is consistent with the results of McGehee et al.[10] who have previously found iodide-rich phases with x ~ 0.2 (i.e., MAPb(I0.8Br0.2)3), following prolonged photoexcitation, McGehee et al.[10], as well as Noh et al.[8] have also suggested that such iodide-enriched phases undergo a crystallographic phase transition from cubic to tetragonal. It is therefore possible that such a phase transition could stabilize an iodide-rich phase. Consequently, an exact quantitative accounting of light-induced phase separation, especially in the iodide-rich limit, would require both crystallographic forms to be considered explicitly within the current DFT and thermodynamic modeling[8, 10].

An approximate estimate for the initial domain size of phase-separated regions can be found from Eq. (1), which vanishes at n = n max, the maximum number of formula units capable of being phase separated by a single hole:

Figure 3b plots Eq. (2) and reveals that for compositions 0.083 < x < 0.5 and for temperatures 250 K < T < 400 K, n max ranges from 2 to 4. Exact n max values are sensitive to the actual composition of iodide- and bromide-rich phases and to the DFT energies. Consequently, they are difficult to predict precisely. Both Eqs. (1) and (2), however, point to the important role played by MAPb(I1−Br)3/MAPbI3 band gap differences in driving phase segregation under optical illumination.

At this point, returning to the kinetic model, we note that the existence of sizable MAPb(I1−Br)3 electron/hole diffusion lengths (l e/h = 100 nm)[21] means that photogenerated carriers created elsewhere in the film can readily access iodide-rich domains (Path 3, Fig. 2a). Consequently, radiative recombination in irradiated MAPb(I1−Br)3 films will be dominated by MAPbI3-like emission given (a) the existence of an intrinsic driving force for carriers to localize in iodide-rich domains, (b) the propensity of photoexcitation to induce MAPb(I1−Br)3 demixing, and (c) sizable l e/h values.

One of the more important outcomes of the kinetic model is an expression for the phase-segregated iodide-rich fraction, ϕ I(t), in the filmwhere at early times. The link between ϕ I(t) and the absorption-based forward rate constant for phase separation (k forward,abs) is made because the absorption measurement directly reports on the phase separation. Eq. (3) is derived in Supplementary Note 1. Using k forward,abs = 0.038 s−1 we find that phase separated fractions nominally start from a few percent and eventually grow to values on the order of 10%. The kinetic model, backed by DFT, thus rationalizes apparent emission/absorption asymmetries in Fig. 1a, b wherein partial phase segregation results in near exclusive iodide-rich emission despite little to no change of the absorption.

We now consider the origin of the I exc-dependent phase separation kinetics seen in Fig. 1c, d. This entails invoking a probabilistic consideration of the emission from iodide-rich domains with a focus on contributions from carrier diffusion (Path 3, Fig. 2a). In this regard, iodide-rich emission originating from diffusion becomes especially important at high excitation intensities, because carrier generation occurs in the presence of multiple segregated iodide-rich domains. A sizable probability therefore exists for carriers to find such a domain within a geometric volume defined by l e/h. This probability increases with increasing excitation intensity due to the slow timescale for phase separation recovery[12] and naturally introduces an I exc-dependence to I iodide. Figure 2b outlines the scenario.

Conceptually, we assume that iodide-rich domains are Poisson distributed within a given carrier diffusion volume, . We additionally assume that diffusing electrons and holes recombine with unit probability from an iodide-rich domain if one exists within V D. From this, we obtain the following expression for the emission intensity of iodide-rich domains

I sat is an empirical emission saturation intensity first seen in Fig. 1c, is the average volume of individual phase-separated domains and at short times. Equation (4) is derived in Supplementary Note 1. Of specific note, Eq. (4) immediately predicts that I iodide(t) grows exponentially. Furthermore, because l e/h is much greater than l (i.e., V D is much greater than ), I iodide saturates even for small values of ϕ I(t). The qualitative behavior seen in Fig. 1c is therefore immediately rationalized.

Next, we compare the phase-separated domain sizes, estimated using the relationship between ϕ I(t) and k (Eq. (4)), to experiment. In particular, from Eqs. (3) and (4) it can be shown that . Using l e/h = 100 nm, k forward,abs = 0.038 s−1 (Fig. 1b) and k forward,em = 0.7 s−1 (Fig. 1a), we estimate a long-timescale domain size of l = 37 nm. This order of magnitude length scale agrees well with the experimental observations of McGehee[10] and Ginsberg[15]. However, it is larger than initial sizes linked to DFT-estimated n max values (below 1 nm). The discrepancy is readily explained by invoking domain growth under prolonged illumination. Qualitatively, domain growth arises because phase separation recovery rates depend on specific surface area and hence increase with decreasing domain size. In support of this, ref.[12] reports a decrease of the recovery rate with increasing sample illumination time.

We now explain the non-linear I exc-dependence of k forward,em in Fig. 1d. Empirically, this behavior suggests that k forward,em is a non-linear function of photogenerated carrier density (ρ). On examining the link between k forward,em and ϕ I(t) in Eq. (4) and invoking Eq. (3) it is apparent that k forward,em takes the functional form where k sat is an empirical saturation rate constant and β is a fitting parameter. Physically, this implies that the probability of creating a phase-separated domain within a geometric volume, V D, remains the same irrespective of the number of photogenerated carriers already in that volume. Consequently, k forward,em can be derived independently by assuming that it is proportional to the probability that at least one carrier exists within V D. As described in Supplementary Note 3, the result is

Subsequent fitting of the data in Fig. 1d with Eq. (5) yields good agreement between model (dashed red line) and experiment.

Phase segregation versus photostability

Beyond qualitatively rationalizing the data, the most significant aspect of Eq. (5) is that it is predictive. This is because, Eq. (5) suggests that k forward,em can be tuned by varying either I exc or MAPb(I1−Br)3 carrier diffusion lengths. To illustrate, Eq. (5) can be used to estimate an intensity threshold below which no phase separation occurs. By finding I exc, where k forward,em is approximately an order of magnitude smaller than k reverse,em, a situation is established where phase separation becomes kinetically unfavorable and where any newly formed iodide- and bromide-rich domains immediately recover to MAPb(I1−Br)3. As shown in Supplementary Note 3, this occurs when I exc = 20 µW cm−2. Remarkably, this value is in excellent agreement with the experimentally-determined threshold value of I exc = 40 µW cm−2 in Fig. 1c.

Alternatively, Eq. (5) shows that lowering l e/h (i.e., V D) can yield k forward,em values comparable to k reverse,em. Under these conditions, phase separation in MAPb(I1−Br)3 films again becomes kinetically unfavorable. To directly test this latter prediction, we have fabricated CsPb(I0.5Br0.5)3 thin films along with their complementary nanocrystal-based films. Electron/hole diffusion lengths are shorter in nanocrystal-based films than in their thin film counterparts because grain boundaries between nanocrystals suppress charge transport[24]. Cesium lead halide perovskites have additionally been chosen because of the better size and composition control of corresponding nanocrystals[25]. Supplementary Note 4 details the synthesis of CsPbBr3 nanocrystals and thin films (see also Supplementary Figs. 6 and 7, Supplementary Table 1).

Figure 4a shows that CsPb(I0.5Br0.5)3 thin films exhibit similar PL spectral redshifts as their MAPb(I1−Br)3 counterparts under illumination (λ exc = 405 nm, I exc = 60 mW cm−2). An initial emission peak at λ mix = 637 nm decreases in intensity and is accompanied by prompt growth of a second feature at 687 nm. As before, this redshift arises due to halide phase separation[7].

Fig. 4

Halide phase segregation versus photostability. Time evolution of CsPb(I0.5Br0.5)3 a thin film and b nanocrystal-based film emission spectra (PL, photoluminescence) under 405 nm continuous wave (CW) excitation (I exc = 60 mW cm−2). c CsPb(I0.5Br0.5)3 thin film and nanocrystal-based film emission peak position during illumination. d Photostability phase diagram of MAPb(I1−Br)3 (MA = CH3NH3 +) obtained from Eq. (5)

By contrast, Fig. 4b reveals that no redshift occurs in CsPb(I0.5Br0.5)3 nanocrystal-based films under the same excitation conditions. If anything, a small 3 nm blueshift is observed, which does not arise from phase separation, as pure CsPbBr3 thin films respond identically under optical excitation. The blueshift thus likely relates to sample instabilities under ambient conditions. Figure 4c compares the emission peak position of a CsPb(I0.5Br0.5)3 thin film and a CsPb(I0.5Br0.5)3 nanocrystal-based film under illumination.

As an added control to verify that phase separation indeed arises from suppressed charge diffusion lengths and not from the limited nanocrystal size in nanocrystal-based films, we show in Supplementary Fig. 8 that at high excitation intensities (I exc = 500 W cm−2) phase separation in nanocrystal-based films can be turned on. We have additionally annealed nanocrystal-based films to sinter component nanoparticles together and have shown that sintered films behave in the same manner as their thin film MAPb(I0.5Br0.5)3 counterparts. At low excitation intensities (i.e., I exc = 2.5 W cm−2) no phase separation occurs. At high excitation intensities (i.e., I exc = 55 W cm−2), clear phase separation-induced redshifting of the emission becomes evident (Supplementary Fig. 9).

Discussion

We conclude that phase separation in mixed halide perovskites can be controlled by deliberately reducing either I exc or carrier diffusion lengths. Sensitivity to either parameter is predicted by Eq. (5) such that depending on application, a balance between the two can be found in order to improve the stability of mixed halide perovskites under illumination. Towards this, Fig. 4d provides a MAPb(I1-Br)3 phase segregation stability diagram constructed from Eq. (5). Assuming that stability occurs when k forward is less than k reverse, the condition k forward = k reverse sets the boundary between phase segregation and photostability. Employing k reverse = 6.8 × 10−3 s−1 (Supplementary Fig. 3) and k sat = 0.7 s−1 (Fig. 1d) yields complementary l e/h and I exc pairs which demark the MAPb(I1−Br)3 optical stability boundary (solid line, Fig. 4d). Green regions show parts of the l e/h/I exc phase diagram photostable under illumination. Associated red regions reveal conditions where phase segregation occurs.

We find that for MAPb(I1−Br)3 thin films under one sun illumination (I exc = 100 mW cm−2) l e/h should be smaller than 13 nm to suppress phase segregation. This insight could therefore aid the development of mesoporous hybrid perovskite solar cells where carrier diffusion lengths are already small and where power conversion efficiencies are currently η = 19%[26]. In addition, while reducing l e/h in corresponding planar architectures may not be optimal for improving η, the choice of cation (e.g., cesium, formamidinium) also affects k reverse. Hence, this represents an alternative way by which to enhance the photostability of mixed halide perovskite films in planar devices[7, 27].

In summary, we have established quantitative insights into the light-induced phase separation of MAPb(I1−Br)3 thin films. Bandgap reduction of iodide-rich domains is found to be the driving force that overcomes unfavorable formation energies to induce iodide and bromide segregation. A DFT-based thermodynamic model shows that entropy dominates formation free energies and provides estimates of the initial phase-separated domain size and composition. These predictions are sensitive to precise structural and entropic models, and quantitative predictions, especially in the iodide-rich limit, would demand a much broader survey of the perovskite crystallographic structure and composition spaces. In particular, it is possible that a crystallographic phase transition additionally stabilizes an iodide-rich phase. Initial phase-separated domains undergo self-limited growth to final sizes on the order of 30 nm. Despite partial phase segregation, the emission properties of MAPb(I1−Br)3 films are near exclusively dictated by nucleated iodide-rich domains. This stems from favorable energetics for hole localization and from the existence of large electron/hole diffusions lengths. Further modeling reveals that MAPb(I1−Br)3 films can be stabilized against phase segregation by reducing carrier diffusion lengths and/or excitation intensities with a phase stability diagram indicating complementary I exc and l e/h pairs yielding photostability. Insights from this study thus set the basis for future experiments to engineer and ultimately control the stability of mixed halide perovskites for photovoltaic as well as other optoelectronic applications.

Methods

Density functional theory

Plane-wave, supercell calculations were performed in the Vienna Ab initio Simulation Package (VASP)[28] within the generalized gradient approximation (GGA)[29] and the projector-augmented wave (PAW)[30] treatment of core states. Valence electrons were expanded to a 500 eV cut-off. Mixed halide structures were computed within orthorhombic supercells of approximate size 9.3 × 12.7 × 8.6 Å with the first Brillouin zone sampled using a 4 × 4 × 4 k-point mesh. Atomic positions and lattice parameters were relaxed simultaneously using the conjugate gradient method. Spin-orbit coupling and non-collinear magnetism were included in the total energy calculations. Further details can be found in the Supplementary Note 2.

Sample preparation

All sample processing was conducted in a N2 filled glove box. A 2 ml solution consisting of: CH3NH3Br (Dyesol) (0.15 M, 34 mg), CH3NH3I (Dyesol) (0.15 M, 48 mg), PbBr2 (98%, Alfa Aesar) (0.15 M, 110 mg), and PbI2 (99.9985%, Alfa Aesar) (0.15 M, 138 mg) in N,N-dimethylformamide (DMF, 99.8%, Sigma-Aldrich) was used as a precursor solution. A 0.2 μm PTFE filter was used to remove any particulates. Precursors were spincoated onto cleaned cover glass (Vetrini Coprioggetto) at 5500 rpm for 30 s with a 1000 rpm s−1 ramp rate. Polymethylmethacrylate (PMMA, Alfa Aesar, MW = 450–550 K) was subsequently deposited over the film by spincoating a 10 mg/ml solution in anhydrous chlorobenzene (2000 rpm for 30 s, 1000 rpm s−1 ramp rate). Resulting films were annealed on a hot plate at 100 °C for 5 min

Emission measurements

Emission spectra were collected with an electron multiplying charged coupled device (EMCCD) camera (Andor) connected to a spectrometer (Acton). A continuous wave 405 nm laser (Coherent, Obis) was used as the excitation source. For power-dependent emission intensity measurements, spectra were collected with variable EMCCD integration times, ranging from 0.1–0.4 s. All emission spectra were fit to two Gaussians (one for λ mix and the other for λ iodide) to estimate the integrated emission intensity of iodide-rich domains (I iodide). Integrated emission intensities of bromide-rich domains were obtained using both a longpass (425 nm, Chroma) and a shortpass (532 nm, Semrock) filter coupled to an avalanche photodiode (Perkin Elmer, SPCM-AQR-14).

For emission-based phase segregation recovery measurements, samples were first irradiated with a CW laser (λ exc = 520 nm, Coherent Obis) with an excitation intensity of I exc = 20 mW cm−2 for 1 min to induce phase segregation. Following this, the laser was turned off. It was subsequently turned on periodically for 50 ms at a significantly lower intensity (I exc = 50 μW cm−2) to acquire PL spectra during recovery under dark conditions. Emission spectra were acquired using a fiber-based spectrometer (Ocean Optics).

Absorption measurements

Complementary absorption measurements during phase separation were conducted by illuminating films with a 405 nm CW laser (I exc = 25 mW cm−2) and periodically acquiring absorption spectra over the course of 30 min using a commercial spectrometer (Varian).

Data availability

The data that support the findings of this study are available from the corresponding authors (M.K. or W.S.) upon request.