Klara
Suchan
ab,
Justus
Just
ac,
Pascal
Becker
a,
Eva L.
Unger
bd and
Thomas
Unold
*a
aHelmholtz-Zentrum Berlin für Materialien und Energie GmbH, Structure and Dynamics of Energy Materials, Hahn-Meitner-Platz 1, 14109 Berlin, Germany. E-mail: unold@helmholtz-berlin.de
bChemical Physics and NanoLund, Lund University, Box 124, SE-22100, Lund, Sweden
cMAX IV Laboratory, Lund University, PO Box 118, SE-22100 Lund, Sweden
dHelmholtz-Zentrum Berlin für Materialien und Energie GmbH, Young Investigator Group Hybrid Materials Formation and Scaling, Kekuléstraße 5, 12489 Berlin, Germany
First published on 5th May 2020
The formation mechanism and the evolution of optoelectronic properties during annealing of chlorine-derived methylammonium lead iodide (MAPbI3−xClx) are investigated in detail combining in situ and ex situ optical and structural characterization. Using in situ optical reflectometry we are able to monitor the evolution of the MAPbI3−xClx phase as a function of time and processing temperature. The formation kinetics is fitted using an improved Johnson–Mehl–Avrami–Kolmogorov model and a delayed formation of MAPbI3−xClx is found when chlorine is present in the precursor. This is verified by X-ray diffraction and X-ray fluorescence measurements. From absolute photoluminescence measurements we determine the implied Voc during film formation, which exhibits a maximum at a specific time during the annealing process. In conjunction with ex situ time-resolved photoluminescence we deduce a decrease in the net doping density for increased annealing times, while the minority carrier lifetime stays constant. We thus demonstrate the potential of in situ optical spectroscopy to monitor and tailor the electronic properties of hybrid perovskites directly during film growth, which can be easily applied to different growth recipes and synthesis environments.
A deeper understanding of how precursors, intermediates and film preparation conditions determine the optoelectronic quality of the metal halide perovskite materials is crucial to: (1) ensure reproducibility of samples within individual laboratories (2) facilitate the adaptation of synthesis protocols and allow comparison between research laboratories and (3) translate device technology to larger areas.15–19 It has been found that the formation of MAPbI3 competes with decomposition into PbI2 such that there is only a relatively small time-window to achieve pure MAPbI3 without precursor residues or impurity phases such as PbI2.12,20 Obtaining materials with well-defined properties reproducibly is crucial to enable the technological utilization of metal-halide perovskites. Ideally, process monitoring and integrated process control would enable to adjust or terminate annealing when optimal material properties are achieved. Lab-based monitoring tools such as optical monitoring are very fast and thus may allow real-time control.21 Although real-time photoluminescence monitoring during perovskite growth has been reported previously22–24 this has so far not been performed in a realistic device fabrication environment e.g. in a glovebox and, even more importantly, so far neither quantitatively analyzed with respect to the evolution of structural and optoelectronic properties of the growing films.
Several studies have focused on understanding the basic mechanisms during the formation of perovskite films.12,25–30 Among the most investigated processes has been the formation of MAPbI3−xClx involving chloride-containing lead precursors such as PbCl2. This synthesis route results in thin-films with long minority charge carrier diffusion lengths and photoluminescence lifetimes as well as low series resistance and good morphology, comparable to materials grown from pure iodide precursors.31,32 Recently, a record Voc = 1.26 V has been demonstrated using such a synthesis route.13 In contrast to the PbI2-based route, there is no need for additional steps such as an anti-solvent drip. The unique combination of good optoelectronic properties and an easy synthesis protocol makes the process exceptionally well suited for up-scaling approaches. The role of chloride in the growth from chloride-containing precursors has been widely debated in literature. Chloride has been suggested to act as a dopant or trap state passivant to rationalize the longer charge carrier lifetimes and higher mobility.31,33 A photoluminescence microscopy study revealed a correlation between an enhanced Cl-content and the local photoluminescence lifetime in individual grains.34 While chlorine has been confirmed to be present in the precursor phase immediately after spin coating, the chlorine content decreases during annealing by loss of methylammoniumchloride (MACl). The chloride-content in annealed films was found to be only a few percent, varying throughout the depth of the film.30,35–38
The formation process itself and especially the crystalline intermediate phases containing chloride were found to be of significant importance, as they template the growth of large grain-size, oriented and high-quality material. Several studies have been carried out to elucidate the key intermediates and film formation kinetics.25,29,30,32,39–41 Stone et al. identified an intermediate phase as MA2PbClI3, consisting of a columnar structure of lead-iodide-chloride octahedra.28 While in previous studies mainly X-ray diffraction-based methods were applied to describe and monitor the growth process in hybrid perovskite thin films,12,20,28–30,40,41 there has been also a recent study in which solar cell properties were simultaneously monitored during growth, by use of a special device structure with interdigitated back contacts.42 Although this study gives information about the evolution of photovoltaic parameters during the growth of these specialized devices, the concept is not easily transferable to an in situ analysis of the optoelectronic properties of real planar devices.
In contrast, the approach followed in the present study focuses on lab-based in situ monitoring of UV/Vis reflectance and photoluminescence during the formation process in a real processing environment (glove box), which can be applied to any real growth process currently used for planar hybrid perovskite devices. In particular, by contact-less monitoring of the evolution of optoelectronic properties during MAPbI3−xClx film formation we are able to identify distinctively different phases during the film formation process and deduce the kinetics of precursor phase consumption, intermediate phase formation and consumption, as well as final product formation. Ex situ powder X-ray diffraction (XRD) as well as X-ray fluorescence (XRF) measurements on samples after defined annealing times are used to identify intermediate phases and quantify the relative chlorine content in the sample as a function of annealing time. Here we furthermore correlate the evolution of the band-gap and photoluminescence peak emission energy with the concentration of MAPbI3−xClx already formed, the photoluminescence quantum yield (PLQY) as well as photoluminescence lifetime to gain insight into the charge carrier dynamics as a function of annealing time and chlorine content. In particular, the PLQY allows us to directly estimate one of the most relevant photovoltaic parameters, the device open-circuit voltage, during the growth/annealing process. This demonstrates that by in situ optical monitoring of the MAPbI3−xClx formation, an optimal end-point of the annealing process can be defined for specific synthesis procedures, opening the way to further significant improvements of performance and reproducibility of halide perovskite based photovoltaics.
To investigate changes in the crystalline phases during the different processing stages, room temperature XRD measurements were performed on samples from process interruptions at defined annealing times of 0 min, 40 min, 45 min, 60 min and 70 min as shown in Fig. 2. Already during the drying stage (t ≤ 0 min), when strong PL emission at 1.7 eV is observed, the XRD pattern shows diffraction peaks from a crystalline phase, in agreement with previous reports.25,28 The XRD pattern exhibits diffraction peaks at 11°, 12°, 15.6°, 16°, 28.5° and 32.7°, which do not agree with those expected for the binary precursor or impurity phases MAI, MACl, PbI2 or PbCl2.
The presence of crystalline precursor phases during annealing and their subsequent decomposition, simultaneously to the formation of MAPbI3 has been discussed for different precursor salts.20,25,27,39 In the case of the PbCl2-derived synthesis route, diffraction peaks comparable with lattice distances in crystalline MAPbCl3 (ref. 25) have led to the conclusion that a pure or chloride-rich methylammonium lead halide perovskite could be a component of the crystalline precursor phase.25,32,43–45 The associated lattice plane distances of the diffraction peaks at 15.6° and 32.7° are indeed consistent with those of MAPbCl3. However, some of the diffraction peaks expected for MAPbCl3, such as the (011) and (022) peaks, are not observed in the powder diffractogram of our samples (see Fig. S1 and S2 in the ESI†). Instead, all diffraction peaks are consistent with a crystalline phase of MA2PbClI3 (ref. 28) which exhibits a similar near-range order but a lower local symmetry than MAPbCl3. In this proposed structure, chloride ions preferentially occupy the axial lattice sites of lead-iodide-chloride octahedra. Therefore, in the [002]-direction, the lattice constant is identical to that of MAPbCl3, resulting in the diffraction peak at 15.6°, which has sometimes been attributed to crystalline MAPbCl3.40,45
During annealing, weak diffraction peaks evolve, that can be assigned to the tetragonal phase of MAPbI3−xClx, while the diffraction peaks of the precursor phase(s) diminish in intensity accordingly. The final film only exhibits peaks corresponding to a pure MAPbI3−xClx phase. A direct comparison of the evolution of the peak at 15.6° attributed to the MA2PbClI3 (310)/(001) precursor and the MAPbI3−xClx (002)/(110) peaks is shown in Fig. 2b. It is noteworthy, that the decrease of the MA2PbClI3 peak occurs in parallel to the evolution of the MAPbI3−xClx peak but at no point in time a solid-solution of both phases exists. It can thus be concluded that the final film is not formed by ion-exchange from MA2PbClI3 but nucleates and grows while consuming the precursor phase.
Photoluminescence and optical reflection signals extracted from the in situ data sets at corresponding annealing times are shown in Fig. 3. The volume fraction of MAPbI3−xClx formed in the film can be estimated from the optical reflection data by employing a Beer–Lambert law based model,21 assuming that light is fully reflected at the silver back mirror and thus transmitted twice through the investigated layer of thickness d:
R(λ,t) = R0(t) + R1(t)exp[−α(λ,t)2dc(t)] | (1) |
No clear absorption edge can be discerned during the initial annealing stage between t = 0 min and t = 35 min. However, the luminescence spectrum of the dried film at t = 0 min, in Fig. 3a, exhibits a peak with a maximum around 1.73 eV and a significantly broadened high energy side.
At about t = 40 min an absorption edge can be recognized in Fig. 3b and the photoluminescence peak redshifts to about 1.6 eV, which agrees with literature values for band-to-band emission of MAPbI3−xClx at 100 °C.47 This is supported by the XRD patterns in Fig. 2a which show an increase of the diffraction peaks corresponding to MAPbI3−xClx between 40 min and 70 min. Thus, we propose that the variable c(t) obtained from fitting the reflection data can be associated with an increasing MAPbI3−xClx volume fraction in the film. Fig. 3 also shows the fits of the reflection spectra indicating that the spectra are well described by the model for annealing times > 35 min. Assuming that at t = 70 min the concentration c = 100% for the fully formed perovskite film with a film thickness of d = 300 nm (compare Fig. S3†), an absorption constant α0 = 1.6 × 105 cm−1 eV−0.5 is derived, which corresponds to an absorption coefficient of α = 2.5 × 104 cm−1 at 1.6 eV. This value is in good agreement with literature values for PbCl2-derived MAPbI3−xClx.48
In Fig. 4 the PL intensity, the PL peak position, the derived MAPbI3−xClx concentration c(t) and the total chlorine content of the film (from XRF) are shown as a function of annealing time. Note that for all ex situ measurements the errors of the individual measurements are very low, whereas a larger uncertainty arises from the sample to sample variation. A comparison of the ex situ measurements and in situ reproducibility for different synthesis runs can be found in the ESI, Fig. S7–S9.† It can be seen in Fig. 4a that although a significant PL yield is present during the whole process, there is a marked increase in PL yield after 40 min followed by a decrease which, will be further analyzed in the following. Correspondingly, the MAPbI3−xClx concentration c(t) obtained from optical reflection measurements (Fig. 4b) shows a sudden increase, agreeing with the increase of the MAPbI3−xClx phase determined from XRD measurements. Thus, the increase in PL yield can be explained by the formation of MAPbI3−xClx. The formation of MAPbI3−xClx does not occur continuously from the beginning of the annealing but only initiates after a certain delay time. This indicates that the formation does not follow a simple first order reaction but is related to a more complex transformation mechanism, as will be discussed further in Section 2.3.
From the evolution of the chlorine content determined from XRF measurements (Fig. 4d), it is noted that the onset of the MAPbI3−xClx formation coincides with a loss of about 50% of the initial chlorine content. This observation is consistent with previous studies, which showed that chlorine is released during formation of MAPbI3−xClx from mixed halide precursors.25,28 Our results indicate that no apparent transformation to MAPbI3−xClx occurs before an apparent “excess” of Cl has evaporated during the initial annealing phase. Please note that we did not observe this delay of MAPbI3 formation in a Cl-free processing route using only iodine-containing precursors (see Fig. S10 in the ESI†). This implies that the delay is caused by the presence of (excess) chlorine in the sample.
The PL peak position also shows an interesting behavior: at t = 0 min the maximum is at Epl ∼ 1.7 eV and thus cannot be attributed to pure MAPbI3, which has a bandgap of 1.6 eV at T = 100 °C.47 During the annealing, a red shift can be observed and the PL evolves continuously into the emission consistent with a pure MAPbI3 phase (or MAPbI3−xClx with small x). This continuous redshift is visible during the entire annealing time. This suggests a continuous evolution of the material luminescent at 1.7 eV at t = 0 into the final MAPbI3−xClx phase. In contrast, in XRD the diffraction peaks corresponding to the MAPbI3−xClx phase do not evolve by a continuous peak shift from the observed crystalline precursor phase. The observed luminescent precursor phase must thus be different from the crystalline precursor phase visible in XRD, but possibly arises from MAPbI3−xClx crystallites already present at the beginning of the annealing process, as will be discussed further in Section 2.3.
From the temperature-corrected PLQY, the quasi-Fermi level splitting μ (or implied open-circuit voltage iVoc) expected at room temperature can be estimated following the approach of Ross50
μ = μrad + kBTln(PLQY/ae) = qiVoc | (2) |
To estimate the minority carrier lifetime and possible changes in the charge carrier (doping) density n0 during annealing, time-resolved photoluminescence (TRPL) measurements were performed ex situ after different annealing times. To account for the concentration of MAPbI3−xClx formed after different times, the initial TRPL amplitude was corrected by the absorption of the incident laser beam. Measured TRPL transients are shown in the ESI in Fig. S13.† We find that the minority carrier lifetime stays nearly constant at τ ∼ 130 ns and even increases with longer times of annealing as seen in Fig. 5c. The externally observed photoluminescence yield is proportional to the radiative recombination coefficient krad, the photoexcited carrier density and the absorptivity at the emission wavelength YPL ∝ kradnpae, which for low injection conditions can be approximated by YPL ∝ τn0ae, where τ is the minority carrier lifetime and n0 is the net doping density.52–54 Although, without exact knowledge of additional constants, this does not allow the deduction of an absolute doping density, the relative change in doping density can be estimated from the PL-yield and carrier lifetime as shown in Fig. 5d. It can be seen that the doping density continuously decreases with processing time by about one order of magnitude. Thus, it can also be concluded that the observed decrease in the implied Voc for longer annealing time is related to a decrease in the doping density rather than to an increase in non-radiative recombination.
A comparison with the pure PbI2 derived synthesis route (ESI, Fig. S10†) shows that after the formation onset the photoluminescence yield shows a very similar behavior for both cases. Thus, it can be assumed that the decrease in doping density is not primarily caused by the decrease in the chlorine content. However, for both routes a ratio of 3:1 of MAI/PbI2 or MAI/PbCl2 was used. Thus MAPbI3/MAPbI3−xClx is grown from an initial excess of methylammonium in both cases. In the case of MAPbI3, this has been shown to result in highly p-doped MAPbI3 (ref. 55) which is primarily attributed to Pb-vacancies exhibiting a rather low formation enthalpy.56 Furthermore, it has been reported that annealing of pristine MAPbI3 reduces the hole density even leading to a conversion to n-type doping.55 Hence, the herein observed decrease in the effective doping density during the formation of MAPbI3−xClx can be assumed to be caused by a decrease of Pb-vacancies caused by the evaporation of MA during annealing.
Fig. 6 Kinetic data for the samples annealed at 80 °C, 85 °C, 90 °C, 95 °C and 100 °C. (a) showing the concentration of MAPbI3−xClx during formation together with the JMAK model fits with an assumed delay (solid lines) and without a delay (broken lines) (eqn (3)). (b) Arrhenius plot of the delay time tonset as determined by the global fit to the JMAK model to extract the activation energy of the delay process, resulting in EOa = 84 ± 7 kJ mol−1. As a comparison ln(kF0) + EFa/RT as obtained from the global fit of the formation kinetics for t > tonset is shown in blue. |
Kinetic parameters of the MAPbI3−xClx formation reaction can be extracted by applying an appropriate temperature-dependent formation model for the time evolution of the MAPbI3−xClx concentration c(t). The sigmoidal kinetics observed in Fig. 6a are characteristic for solid-state transformations and can be attributed to a nucleation and growth process. The widely used Johnson–Mehl–Avrami–Kolmogorov (JMAK) model57–59 describes such solid-state transformation reactions at constant temperature T and has been recently applied to the formation of perovskites from precursors.12 Model assumptions are that nucleation occurs homogenously distributed within the precursor phase, that the growth rate is constant in time and that growth occurs homogeneously in all directions.58 We apply the JMAK-model using a thermally activated rate equation that includes a delay time tonset to account for the observed delay in MAPbI3−xClx formation:
(3) |
The dataset shown in Fig. 6a is fitted for all temperatures simultaneously, varying k0, EFa and n globally for all processes. Alternatively, tonset is fixed to zero for the pure JMAK model while it assumes temperature-dependent finite values for the altered model given by eqn (3). From fitting the data in Fig. 6, it can be seen, that if the delay time is omitted (tonset = 0 min) the global fit insufficiently describes the measured data and results in a growth dimensionality of n = 5.5, which is physically difficult to rationalize. When introducing a delay time tonset(T), the experimental data can be fitted with higher accuracy for all temperatures, which can be also seen from the χ2-coefficient which decreases from 11 to 6. The resulting activation energy of the formation reaction is EFa = (94 ± 2) kJ mol−1 and the growth exponent is n = 2.1 ± 0.1. The growth exponent corresponds to the dimensionality of the growth, where n = 2 implies 2-dimensional growth of previously nucleated material. As a growth exponent n > 4 bears no physical meaning, this shows that the formation of MAPbI3−xClx can only be correctly described by including a temporal delay in the onset of MAPbI3−xClx formation (eqn (3)).
In addition, a second activation energy can be extracted from the delay times. Fig. 6b shows an Arrhenius plot of the delay times tonset yielding an activation energy of EOa = (84 ± 7) kJ mol−1. Although both activation energies, estimated from the temperature dependence of tonset(T) and from the temperature dependence of the subsequent formation kinetics, are very close to each other we believe that they may correspond to two distinct phenomena, as discussed further below.12,60
A comparison of the reaction kinetics determined in the present study with previously reported reaction kinetics observed by in situ GIWAXS shows some agreement and some differences (Fig. S14 in ESI†). While Chang,40 Unger25 and Stone28 observe a delay, this is not observed by Moore et al.12 The reported duration for complete conversion of the precursor to the final film varies from 30 min to 70 min at 100 °C and the reported activation energies vary from 69 kJ mol−160 to 86 ± 6 kJ mol−1.12 It has to be noted, that these activation energies were obtained without including a formation delay. A comparison of the previously presented data shows that even though the formation onset differs significantly in different studies, the reaction kinetics after the formation onset obtained in the present study agree with the data presented by Chang, Unger and Stone25,28,40 and hence seem to be more fundamental than the delay time (see Fig. S14†). Thus, we propose that variations observed for the formation onset strongly depend on the processing conditions. The annealing environment, N2 or dry air, varying relative humidity or ambient temperature is often not specified in detail but may have a crucial influence. This variance in process conditions and kinetics serves as a strong argument for the relevance of in situ process control during growth.
The strong variation of tonset visible in the literature data and in this study for different glovebox temperatures (see ESI, Fig. S7b†), together with the largely unaffected subsequent formation kinetics, indicate that the processes occurring during the delay and the formation phase itself may be distinctively different.
Considering all data from PL, UV-Vis, XRD and XRF we conclude, that MACl and MA2PbClI3 as well as MAPbI3−xClx – nuclei coexist within the film from the beginning. The fact that the MAPbI3−xClx cannot be observed by XRD in the beginning of the growth process can be explained by an initially small-domain size of the crystallites. We recall from Fig. 1 and 4 that a photoluminescence signal is observed from t = 0 min which starts at a peak energy of 1.7 eV and then continuously red-shifts to 1.6 eV. Although such a continuous red shift may be explained by a growth of the MAPbI3−xClx nuclei and thus decreasing of quantum confinement, it may also be due to a continuous substitution of chlorine with iodine within the MAPbI3−xClx seed crystals. Assuming a linear dependence of the band-gap of MAPbI3−xClx on the chlorine content x, the PL peak energy would indicate a Cl-concentration of 6% of the original seed crystals at t = 0 min, subsequently reducing to about 0.6% for the fully annealed MAPbI3−xClx film. This explanation is consistent with the Rietveld refinement of the powder diffractograms which show an increase of the unit cell volume approaching MAPbI3 literature values for t > 40 min (see Fig. S2 in the ESI†), indicating a decrease of the Cl content in MAPbI3−xClx from that point onwards.43,61,62
MACl was shown to evaporate during the process25 which correlates to a loss of chlorine within the film, as measured with XRF (Fig. 4). Once MACl is fully evaporated, further evaporation of chlorine leads to a decomposition of the MA2PbClI3 precursor. By this, the limiting elements (Pb/I) become available and thus a sudden formation onset of MAPbI3−xClx can be observed. Thus, we propose an activation energy of EOa = (84 ± 7) kJ mol−1 for the loss of Cl in form of MACl during the latent phase and an activation energy EFa = (94 ± 2) kJ mol−1 for the subsequent formation of MAPbI3−xClx and decomposition of MA2PbClI3.
Analyzing the MAI/PbCl2 based synthesis route it was found that a crystalline precursor phase develops which can be attributed to MA2PbClI3. Furthermore, the excess of chlorine is shown to delay the onset of MAPbI3−xClx formation. We propose that already from the beginning of the growth process MAPbI3−xClx nuclei exist, (with strong luminescence at 1.7 eV) from which the MAPbI3−xClx layer forms only after a delay time tonset. This is supported by a growth exponent of n = 2 for the formation of MAPbI3−xClx, indicating 2D growth with previous nucleation, and a PL emission which continuously approaches the final MAPbI3−xClx emission. By correlation with TRPL measurements we observed that the excited charge carrier lifetime stays nearly constant during annealing, however the effective doping density decreases significantly. This can be attributed to a loss of methylammoniumchloride during annealing. Thus, it is shown that through in situ monitoring the doping density could be fine-tuned by altering the annealing time. This work demonstrates the great potential of in situ optical spectroscopy in order to monitor and tailor the optoelectronic properties of hybrid perovskites directly during device processing. Furthermore, the possibility to analyze real growth processes in detail opens the way to an improved understanding and enhanced reproducibility of halide perovskite solution processing, which is a necessary prerequisite for up scaling approaches.
Lead(II) chloride was used as purchased from Sigma Aldrich. MAI and lead chloride were dissolved in anhydrous N,N-dimethylformamide (DMF) in a 3:1 ratio10 and stirred for 30 min at 20 °C in a nitrogen-filled glovebox. 150 μl of the solution were drop cast onto the TiO2 coated substrate and spin coated at 2000 rpm for 45 s. After spin coating, the samples were first left to dry for 30 min at 20 °C (ref. 63) and afterwards annealed on a hotplate at 100 °C until the perovskite was fully formed.10
During all room temperature measurements the temperature in the glovebox was kept constant at 20 ± 0.5 °C. The temperature of the samples was determined using a Pt100 resistor.
X-ray fluorescence was measured in a custom system (IFG) using a rhodium tube with a micro-focus capillary optics with a focus spot size of 100 μm. The fluorescence signal is detected with an energy dispersive detector with a resolution of 190 eV. For the measurements of the Cl-content a mean value and the standard deviation was calculated from 4 spots of each sample. The standard deviation of the Cl content lies between 0.1 and 0.2% and thus is smaller than the marker size in Fig. 4d.
Time-resolved photoluminescence was measured by excitation with a pulsed laser with 660 nm wavelength (pulse duration 100 ps) detection with a Si-Avalanche Photo Diode (time resolution 60 ps). Repetition rates of either 1 MHz or 250 kHz were used, depending on the transients measured. The excitation density, calculated from the time-averaged laser power, the repetition rate of the pulses, the spot size of the laser on the sample and absorption coefficient was Φ = 2.45 × 1015 photons per cm3 per pulse.
Footnote |
† Electronic supplementary information (ESI) available. See DOI: 10.1039/d0ta01237h |
This journal is © The Royal Society of Chemistry 2020 |