Low-Temperature Emission Dynamics of Methylammonium Lead Bromide Hybrid Perovskite Thin Films at the Sub-Micrometer Scale

We study the low-temperature (T = 4.7 K) emission dynamics of a thin film of methylammonium lead bromide (MAPbBr3), prepared via the anti-solvent method. Using intensity-dependent (over 5 decades) hyperspectral microscopy under quasi-resonant (532 nm) continuous wave excitation, we revealed spatial inhomogeneities in the thin film emission. This was drastically different at the band-edge (∼550 nm, sharp peaks) than in the emission tail (∼568 nm, continuum of emission). We are able to observe regions of the film at the micrometer scale where emission is dominated by excitons, in between regions of trap emission. Varying the density of absorbed photons by the MAPbBr3 thin films, two-color fluorescence lifetime imaging microscopy unraveled the emission dynamics: a fast, resolution-limited (∼200 ps) monoexponential tangled with a stretched exponential decay. We associate the first to the relaxation of excitons and the latter to trap emission dynamics. The obtained stretching exponents can be interpreted as the result of a two-dimensional electron diffusion process: Förster resonant transfer mechanism. Furthermore, the non-vanishing fast monoexponential component even in the tail of the MAPbBr3 emission indicates the subsistence of localized excitons. Finally, we estimate the density of traps in MAPbBr3 thin films prepared using the anti-solvent method at n∼1017 cm−3.


Introduction
Halide-based, hybrid organic-inorganic, perovskite materials have attracted the world's attention since 2012 through the demonstration of a low-cost solar cell approaching energy conversion records [1][2][3]. Since then, halide perovskite materials have successfully made their ways through multiple arrays of opto-electronics: solar cells [4], LEDs [5], lasers [6], photodetectors [7], solar concentrators [8], mechanical energy harvesting [9] or scintillation devices [10,11]. While methylammonium lead iodide (MAPI) has been arguably the most successful halide perovskite so far for application in solar cells [12,13], methylammonium lead bromide (MAPbBr 3 ) thin films or single crystals find applications in LEDs, scintillation and photodetector devices [10,14,15]. The dynamics and origin of their emission have thus been the center of several studies and are still highly debated subjects [16][17][18][19][20][21][22][23][24][25]. In order to improve device quality, it is also important to understand the spatial dependence of their emission, i.e., performing not only ensemble measurements, but also spatially resolved experiments, at the micrometer scale [26,27]. Indeed, while it was originally believed that the grain boundaries in polycrystalline thin films were causing no harm to device efficiency [28,29], it has been demonstrated otherwise [26,27,30]. Concerning the spatial characterization of the emission, much effort has been focused on the spatial inhomogeneity of the emission of MAPI [26,30,31], but information is more sparse for MAPbBr 3 thin films [27,32,33]. Usually, only two of the three following information sources are available: spatial, spectral and temporal. It is, however, efficient to identify the emitting entities with the spectral information and the corresponding emission dynamic mechanism with fast temporal information [16,17,20] and couple it to spatial information [26,30]. Concerning the emission dynamics of MAPbBr 3 thin films (contrary to MAPbBr 3 single crystals [16,[34][35][36]), most of the studies were performed at room temperature, where the thermal energy is large enough to ionize the shallow defects and observe free carrier relaxation [19]. There have been very few studies at low temperature, down to T = 10 K, where the shallow defect recombination dynamics are modified and free carrier relaxation is not expected [17,20,25,[37][38][39][40][41]. Low-temperature experiments enable one to probe thin film dynamics (and more specifically MAPbBr 3 in our case) without being hampered by the following: strong exciton-phonon interactions that take place at room temperature and create strong overlap between the emission peaks [25], organic cations' rotational motion [17], or free carrier emission [38,41]. It is thus a unique tool that allows one to focus solely on the excitons' and the traps' emissions. Furthermore, while most devices based on MAPbBr 3 (and especially LEDs) operate at room temperature, the differences between room-and low-temperature emissions [39] may lead to a better understanding of some physical aspects that govern the performances of LEDs [25]. Moreover, it is now well-known that MAPbBr 3 undergoes phase transitions as the temperature is lowered [17,38,42] and that working at cryogenic temperatures unlocks the possibility to probe the physics of the orthorombic crystalline phase, which is not accessible at room temperature. Finally, it is important to study the low-temperature response of hybrid perovskite thin films in the framework of new applications where they can serve as a passive medium for quantum optics applications [43]. Hyperspectral and fluorescence lifetime imaging microscopy (FLIM) performed at low temperature are interesting techniques to access altogether spectral, temporal and spatial information at the sub-micrometer scale and unravel the aforementioned properties. In this paper, we present a study which joins the four properties: spatial, spectral or temporal are observed at low temperature (T = 4.7 K) for the emission of MAPbBr 3 thin films prepared via the anti-solvent method. The first part is dedicated to continuous wave (CW) experiments performed with an excitation quasi-resonant with the MAPbBr 3 band-edge emission. Intensity-dependent hyperspectral imaging enables us to locate spatially and spectrally the emissions of excitons and defects. In the second part, two-color FLIM experiments are presented, which allow us to estimate the density of defects in the MAPbBr 3 thin film, as well as the main process of recombination through these traps and its evolution with the density of absorbed photons. Our work permits us to image and discriminate for the first time at T = 4 K the regions of emission arising from the excitons or the traps at the micrometer scale, and subsequently propose a model that explains the trap emission via a Förster resonant transfer of diffusing electrons recombining with static holes. We find that, contrary to another report, the excitons' dynamics are represented by a fast exponential and not a long-tail process [17]. These findings can have implications in improving the performances of devices based on light harvesting, light emission and future low-temperature applications of the hybrid perovskite thin films.

Methylamonium Lead Bromide Thin Film Preparation
The MAPbBr 3 thin film was prepared as described in Ref. [43]. Briefly, 112 mg of methylamonium bromide (MABr, 98%, Sigma Aldrich) powder and 367 mg of lead bromide (99.999%, perovskite grade, Sigma Aldrich) were dissolved in 2 mL of N-N-Dimethylformamide (anhydrous DMF 99.8%, Sigma Aldrich) leading to a 0.5 M solution of MAPbBr 3 precursor. Then, 100 µL of that solution was drop-cast on a microscope coverslip (BK7) at 5000 rpm for 30 s. Five seconds after starting the spin-coater, 500 µL of chlorobenzene (anhydrous 99.8%, Sigma Aldrich) was injected to enforce homogeneous crystallization of the MAPbBr 3 thin film. The result was a highly oriented, transparent, crystalline thin film (see Supplementary Materials Section S1) with a thickness of 90 nm measured with an atomic force microscope (AFM) [43].

Low-Temperature Confocal Microscopy
The setup is described in detail in the Supplementary Materials, Section S2. The sample was placed on three-axis steppers and scanners in a home-built confocal microscope coupled to a helium recycling cryostat. The system allows one to perform experiments at T = 4.7 K at the sample location. The excitation was provided by either a picosecond pulsed laser emitting at 445 nm, or a CW laser emitting at 532 nm. The excitation was sent to a half-waveplate and a polarizer, then entered a polarizing beam-splitter (PBS) which reflected the laser beam down to a quarter-waveplate and finally reached the microscope objective (which was at T = 4.7 K) as circularly polarized light. The microscope objective focused the beam onto the sample and collected the luminescence. The collected PL was redirected to the PBS and was coupled into a single mode fiber and sent to a 300 mm focal length spectrometer and a Handbury Brown and Twiss interferomter (HBT) via a 2 × 2 fiber coupler.

Hyperspectral Microscopy
For hyperspectral imaging measurements, the sample was excited with a 532 nm CW laser. The sample was raster-scanned through the laser beam by moving the xy lowtemperature scanner. The scanning range was 9 × 9 µm 2 at T = 4.7 K. Light collected from the microscope was sent to a 303 mm focal length monochromator and dispersed by a 600 line/mm grating. Spectral information was recorded by an emCCD camera cooled down at minus 80 • C. The system spectral and spatial resolution were 0.24 nm and 500 nm, respectively.

Two-Color Fluorescence Lifetime Imaging Microscopy
For FLIM measurements, the sample was excited with a picosecond pulsed laser at 445 nm. The sample was raster-scanned through the laser beam by moving the xy lowtemperature scanner. The scanning range was 9 × 9 µm 2 at T = 4.7 K. Light collected from the microscope was sent to an HBT interferometer consisting of two single-photon avalanche photodiode (SPAD) detectors (Micron Photon Devices, PDM) coupled by a 50:50 R/T beamsplitter. In front of the first SPAD (DET1), we have placed a longpass filter at 550 nm (Thorlabs, FELH0550) and a bandpass filter centered at 550 nm (Edmund optics), generating a 5 nm detection window ranging from 550 to 555 nm. In front of the second detector (DET2), a bandpass filter centered at 568 nm with 10 nm bandwidth was placed, creating the second detection window at 568 ± 5 nm. Both detectors recorded the PL signal from the microscope simultaneously, allowing for two-color FLIM experiments. The SPADs were connected to a multi-channel photon counting electronics with 80 ps time bin. The system timing and spatial resolution before reconvolution-based fitting were 240 ps and 700 nm, respectively. Up to four repetition rates were used: 5, 1, 0.5 and 0.1 MHz, with an integration time per pixel of 5 s for the first two, and 10 s for the last two.

Fitting the Fluorescence Lifetime Imaging Experiments
Each FLIM image contains 5000 decay curves, all of which were fitted using the "curve fit" algorithm from the SciPy library in Python. In order to keep the free parameters of the fit to a minimum, the photon background was estimated for every curve to the median value of the first 10 points in the decay curve. We used re-convolution with the IRF (see Supplementary Materials, Figure S4) of our system to obtain more sensible results for the fast component.

Ensemble Measurements
We present in Figure 1 ensemble characterizations of a MAPbBr 3 thin film during the sample cool-down from T = 300 to 4.7 K. Figure 1a presents the evolution of the normalized emission of the MAPbBr 3 thin film under an excitation of λ ext = 445 nm, as a function of time.
The emission spectrum has its maximum shifting from λ max = 532 nm at 300 K to 552 nm at 4.7 K. This is in accordance with previous results in the literature [17]. We note a discrete change in the emission wavelength at a temperature of T∼150 K, which might be the subtle sign of the MAPbBr 3 thin film transition from the tetragonal to orthorombic phase [17,38,42]. Figure 1b presents five examples of spectra recorded at five different temperatures: 300 K (violet), 180 K (red), 100 K (green), 20 K (blue) and 4.7 K (black). Intensities are presented in log scale. At high temperature, the spectrum is described by two peaks [16,17,43], while at low temperature there is the apparition of an emission tail at longer wavelength. At room temperature, the first high-energy peak at 532 nm is usually attributed to excitons [17,38,40], while the origin of the low-energy peak is still debated [17]. Figure 1c presents the evolution of the linewidth of the emission of the MAPbBr 3 thin film as a function of temperature. It varies from 35 meV at 300 K to 12 meV at 4.7 K. We note that we do not see a signature of the phase transition in the emission linewidth, while we see a hint of it on the emission wavelength.

Hyperspectral Microscopy
At T = 4.7 K, we have recorded a hyperspectral, 9 × 9 µm 2 image of Zone 1, at an excitation wavelength of 532 nm and an intensity of 100 nW/µm 2 , where each spatial pixel encompasses a spectrum ranging from 500 to 600 nm. Results obtained for different wavelength ranges are presented in Figure 2a-d. Figure 2a shows the signal recorded around 532 nm, i.e., the excitation laser wavelength. It is thus not an image of the MAPbBr 3 emission, but a confocal retro-reflectivity image of the sample surface. It gives us information on the topography of the sample. We observe that the thin film is locally structured with polycrystalline domains of ∼1 µm wide surrounded by the grain boundaries. This is in accordance with SEM images of our thin films [43]. Summing the spectra in every pixel for 540 nm ≤ λ ≤ 580 nm, we obtain the image in Figure 2b, which represents a confocal PL image of Zone 1. Except from a high-intensity spot in the bottom-left of the image, the PL signal from the thin film is rather homogeneous, varying by less than a factor of 2 over the whole image. We notice a moderate correspondence between the two images (topography and full-integrated PL) corresponding to a Peasron's r correlation coefficient between the two images of r = 0.41 (r takes any value between −1 and 1, r = −1 being a complete anticorrelation of the data, r = 0 the absence of correlation and r = 1 a complete correlation). We present in Figure 2c,d two images obtained by integrating the spectra in every pixel with a bandwidth of 0.8 nm at 545 nm ( Figure 2c) and 561 nm ( Figure 2d). The image in Figure 2c seems drastically different from the one with fully integrated spectra in Figure 2b. The corresponding Pearson's r between (c) and (b) reads at r = −0.01, i.e., no observable correlation. Figure 2c is structured with resolution-limited (∼500 nm), low-density emission spots, which have a strong signal-to-background dynamic of about 100. There are approximately 30 spots in the 81 µm 2 of Figure 2c, i.e., less than 0.5 spot per micron-squared. These emission spots appear to be randomly distributed over the Zone 1 surface. In contrast, Figure 2d (signal at 561 nm) exhibits a large collection of emission spots, leading to a continuous signal and a smaller dynamic in the image, with a max-to-min reading at about 6. There is a clear asymmetry in the emission signal, from low at the top-right, to high in the bottom-left. In this case, it more closely resembles Figure 2b, the fully integrated spectra image. Indeed, the Pearson's r between those two images ((b) and (d)) is r = 0.43, i.e., a moderate correlation. Finally, the Pearson's r between the high-( Figure 2c) and low-energy ( Figure 2d) seats is r = −0.12, denoting a weak anti-correlation between the two signals. Figure 2a-d thus leads us to conclude that the emission at low temperature of a MAPbBr 3 thin film is strongly spatially and wavelength dependent with an emission spectrum extremely richer than what has been measured before on ensemble measurements [17,40]. Figure 2e presents the sum of all spectra from all the pixels in Figure 2b. We obtain a continuous curve, with a strong asymmetric tail to the low-energy areaand a double-peak, whose maxima are located at 552 and 555 nm. Figure 2f presents spectra from single pixels located by colored-circle on the four images in Figure 2a-d. The colors of the spectra correspond to the pixels in the center of the similarly colored circles. Surprisingly, every pixel exhibits a different spectrum. Moreover, while the spectra are all quantitatively different, they all exhibit a structure with multi-peaks at high energy and a continuous tail at low energy. These narrow peaks (<1 nm) are located in the spectral region ∼545-560 nm. While in the region ∼545-550 nm the peaks seem to be isolated, those located in the region ∼550-560 nm are rising on top of what appears to be a continuous background.

Power Dependence
We present in Figure 3a an intensity dependence of the shape of the spectrum from the pixel in the center of the orange circle from Figure 2. The intensity varies over 5 orders of magnitude from 1 nW/µm 2 (black line) to 100 µW/µm 2 (red line). In between, we find the spectra recorded at 100 nW/µm 2 (blue line) and 10 µW/µm 2 (green line). The spectra are normalized to their respective maximum. The intensity of the low-energy part of the spectrum decreases relative to the band-edge emission. In order to characterize further this behavior, we have recorded hyperspectral images of Zone 1 for those six excitation intensities covering 5 decades, from to 100 µW/µm 2 down to 1 nW/µm 2 at an excitation wavelength of λ exc = 532 nm. We present in Figure 3b an example of the corresponding power dependence of the emission obtained on a single pixel from Zone 1 for two different emission wavelengths, namely, 549 nm (black dots) and 568 nm (blue triangles). The two curves are plotted in a log-log manner and their behavior can be well described by a linear interpolation, i.e., a power law in linear coordinates. We fit those two set of data with the following equation: where S, A, I and X are the PL signal, the amplitude, the excitation intensity and the power law exponent, respectively. We find X = 0.91 and 0.67 at 549 and 569 nm, respectively, indi-cating a quasi-linear increase in the PL of the MAPbBr 3 thin film at 549 nm (characteristic of exciton emissions [44,45]), but a strongly sub-linear behavior at 568 nm, suggesting a saturation effect (characteristic of a finite number of traps [17,45]). Extending this analysis to the entire Zone 1, we present in Figure 3c-f examples of the obtained hyperspectral images integrated over different spectral regions, with the same bandwidth of 0.8 nm.
The full hyperspectral data are presented in the Supplementary Materials in the form of a series of seven videos (one for every decade of powers and one for the final image of the spatial distribution of X), where one frame represents 0.08 nm: it highlights the strong spectral and spatial inhomogeneity of the MAPbBr 3 thin films, which cannot be fully presented in the manuscript. In Figure 3c,d, the spectra are integrated over 548.3 ± 0.4 nm for excitation intensities of 100 µW/µm 2 and 10 nW/µm 2 , respectively. Figure 3e,f was obtained by integrating the spectra over 568 ± 0.4 nm for the same intensities as (c) and (d), respectively. All contour plots are presented with intensity in log scale. A more complete figure is presented in the Supplementary Materials, Figure S3. More precisely, Figure S3d,f presents the imaging for spectra integrated over 561 ± 0.4 nm. The images in Figure 3c,d are highly comparable, pointing to a spatially homogeneous evolution of the signal at that emission wavelength. Indeed, while the PL signal is inhomogeneously distributed over the surface, its evolution as a function of intensity at the band-edge proves to be homogeneous. However, as the emission wavelength increases, the image changes from low to high intensities, indicating non-linearity in the evolution of its spatial inhomogeneity: as the power is increased, the image becomes more and more homogeneous, indicating that the emission from some regions in the image are experiencing a saturation that is not experienced for images at the band-edge. We present in Figure 3g,h and in Figure S3k, 100 × 100 pixel images representing the spatial distribution of X as a colorscale over Zone 1, for emission wavelengths 548 nm (Figure 3g), 561 nm ( Figure S3k) and 568 nm (Figure 3h). On a general level, we see that X decreases on average from high to low emission energy. We find the following average values: < X 548 > = 0.9 for the image at 548 nm, < X 561 > = 0.76 for the image at 561 nm and < X 568 > = 0.68. Looking at the details, we find that for the image at 548 nm ( Figure 3g) the signals from the inter-spaces between spots of high signal behave linearly with the excitation intensity (green color represents X∼1), while the spots themselves behave sub-linearly, with a red color, indicating X∼0.75. At an emission wavelength close to 560 nm ( Figure S3k), the power law exponent X reveals to be homogeneous for the whole image around X∼0.75. While we can conclude for Figure 3g that the green regions result from the emission of linear entities with excitation power, the equivalent straight conclusion about only sub-linear entities emitting in Figures S3k and 3h cannot easily be drawn. Indeed, the emission could be the superposition of linearly behaving entities and sub-linear ones, leading to X < 1. To support that assumption, one can still observe faint green regions over Figure S3k. Finally, at even larger emission wavelengths, 568 nm (Figure 3h) in the emission tail, we find that the region with the highest signal behaves similarly to the square root of the excitation intensity with X∼0.5 (blue regions), while the lower-signal regions behave like those at 561 nm and on the spots of high signal at the band-edge, i.e., with X∼0.75. From these power dependence experiments, we can conclude that part of the band-edge emission behaves linearly with the excitation intensity, as should the exciton. Whereas, as the wavelength increases, most of the emission in the image behaves sub-linearly, even similar to the square root of the intensity for the regions of Zone 1 at 568 nm. At the band-edge, we are able to locate spatially the micrometer-scale regions over which the excitons' emission dominated the signal. We conclude from the latter results that the emission at longer wavelength than the band-edge is a subtle mix between excitons and defects; however, these CW experiments do not permit us to unwrap the mixing, its evolution within the excitation intensity, nor to determine the emission mechanism at the tail of the emission. (e,f) The same intensities used in (c,d) but for emission spectra integrated over 568 ± 0.40 nm. (g,h): Image of 9 × 9 µm 2 , 100 × 100 pixel representing X as a colorscale obtained from the fitting procedure (following Equation (1) performed for every pixel of Zone 1 for 6 excitation intensities ranging over 5 decades from 100 µW/µm 2 down to 1 nW/µm 2 .

Fluorescence Lifetime Imaging Microscopy Under Pulsed, Non-Resonant Excitation
In order to understand in more detail this behavior, we have performed time-resolved measurements, which will give us access to the emission dynamics of excitons and traps. We present in the following section results from FLIM experiments performed on a 9 × 9 µm 2 , 50 × 50 pixel area (Zone 2). To perform these experiments, part of the signal is directed to a Handbury Brown and Twiss intererometer (HBT). In front of the first detector (DET1) the detection window ranges from 550 to 555 nm, i.e., close to the band-edge. In front of the second detector (DET2) the detection window is 568 ± 5 nm. The bandwidth of DET1 and DET2 are represented as dashed and dotted lines, respectively, in Figure 2e,f. Details of the setup are presented in the Supplementary Materials, Sections S2 and S4.

Exciton-and Trap-Related Emission Dynamics
We present in Figure 4a,b the PL intensity map of Zone 2 from DET1 ((a), 550-555 nm) and DET2 ((b), 568 ± 5 nm) obtained at a repetition rate of 500 kHz. Because of the larger bandwidth of the images in (a) and (b), compared to the images of Zone 1 (shown with 0.8 nm bandwidth), the visual correspondence between (a) and (b) appears to be larger, while taken at center wavelength 16 nm away. This is confirmed by the Pearson's r found to be r = 0.65, indicating a positive correlation between the two images. The correlation is, however, only partial (r < 1), highlighting the difference induced by the two spectral windows used to record the images. The signal over the surface is rather homogeneous since, for both images, the full-scale dynamic of the signal is close to 2. Figure 4c presents two examples of decay curves obtained for the same pixel in Zone 2, indicated by a dashed cross in both images in Figure 4a,b. Results are presented in log-log scale to capture visually the long-as well as the short-time dynamics. The red circles in Figure 4c correspond to the lifetime close to the band-edge around ∼550 nm (Figure 4a), and the blue circles the decay around 568 nm (Figure 4b). The curves are normalized to their respective maximum. Both curves exhibit a similar qualitative behavior: we can identify a fast and a slow process. However, the amplitude of the slow component appears to be largely diminished at wavelengths closer to the band-edge (500-550 nm) compared to wavelengths identified to be defect-related with the CW experiments (568 nm). From that correspondence, we assign the fast component to the exciton dynamics, and the stretched exponential to trap emission. While it is difficult to assign a precise behavior at first sight to the fast component, which is tangled with the instrument response function (IRF, see Supplementary Materials Figure S4), it is straightforward to conclude with the log-log plot that the slow process is neither mono-or bi-exponential but buries a complex, multi-exponential structure. Our different fitting attempts eliminated a strict power law to describe the behavior of the long-lived component, and we found that a stretched exponential was a better fit for both channels. Stretched exponential dynamics in hybrid perovskite have been observed in thin films of MAPbCl x I 3−x [26,46,47] and in MAPbBr 3 [17,45]. It is not surprising to observe complex decay dynamics in the emission tail (usually denoted an Urbach tail) since it is acknowledged that its emission is partially driven by intrinsic disorder and defects [19]. The association of those two processes, one fast mono-exponential in addition to a stretched exponential, was used to fit the data as follows: where y 0 , A 1 , A 2 , γ 1 , γ 2 and β are the background, the amplitudes of the fast and stretched components, the decay rate of the fast and stretched component and the stretching exponent. We obtained the results presented in Figure 4c where the fits are shown as black lines. For the two fits, we find γ 1 ∼ 4-5 ns −1 (corresponding to a decay time of the fast process equivalent to τ 1 = 1/γ 1 ∼ 200-250 ps). Similar fast decay times are observed in CsPbBr 3 perovskite quantum dots [48] at T = 4 K or in CsPbBr 3 single crystals [16]. Regarding the trap-related relaxation, we find β = 0.32 and 0.22 for 550 and 568 nm. Stretched exponential dynamics can be the mathematical description of a diffusion process [17,47,49] and do not describe typical donor-acceptor static relaxation [17,49], which resembles a stretched exponential if one is not careful enough to record the full dynamics over several orders of magnitude in the time and signal axes [17]. On the other hand, electrons diffusing over randomly distributed traps, and recombining with static holes on acceptors can give rise to a stretched exponential [49]. For example, the relaxation function probability φ(t) of an electron diffusing through a collection of randomly distributed traps via the Förster direct-transfer model, based on a R −6 interaction (where R is the distance to the trap), is written as [49]: where γ is the relaxation rate and β = D/6 for a process in D dimensions. Thus, we expect β to vary from 0.5 to 0.17 as the dimension of the diffusion goes from D = 3 to D = 1. For a 2D Förster direct transfer process, we expect β = 1/3. This is close to our experimental value found for defects emitting in the 550-555 nm region (β = 0.32), but significantly larger than what is observed around the 568 nm emission wavelength. The model described by Equation (3) takes into account transfer via parallel channels to possibly all the defects of the distribution [49]. Limiting the transfer possibilities to only the nearest neighbors, examining the D = 1 case (which can be solved analytically [49]), the relaxation function probability depends on β NN = 1/7∼0.14 [49] instead of β = 1/6∼0.17, i.e., a 20% decrease. We can speculate that the stretching exponent in the 2D case also decreases when only the nearest neighbors are taken into account. From these assumptions, we can conclude that a 2D Förster direct transfer diffusion process limited to the nearest neighbors could describe the dynamics of emission of the defects in the tail of the MAPbBr 3 emission. Whereas, the diffusion of the electrons in traps that exhibit an emission close to band-edge seem to have access to further-away traps, since the value of β matches that of the 2D case with no diffusion restriction. One feature which is not observed in other low-temperature reports [17], and only in a few at room temperature [45], is the remaining of the fast component associated to the exciton in the emission tail. This is a signature of the emission of localized excitons confined in minima of the local potential fluctuations, previously observed at room temperature in MAPbBr 3 [45]. We have repeated these experiments for different repetition rates, from 5 MHz to 100 kHz, keeping the average incident power constant. This enables us to vary the short-time density of photons absorbed per laser pulse by the MAPbBr 3 thin film, n i (with i = 5, 1, 0.5 and 0.1 being the repetition rates in MHz), without changing the number of photons absorbed by the thin film over a longer time average. It can be calculated as follows: where N i is the number of photons per laser pulse, i = 5, 1, 0.5 and 0.1 (after the laser repetition rate in MHz), < P > = 240 nW/µm 2 the average power onto the MAPbBr 3 thin film, E phot = 2.8 eV is the energy of a laser photon, f rep the laser repetition rate (5, 1, 0.5 or 0.1 MHz), α = 0.5 the fraction of light being absorbed by the MAPbBr 3 thin film at λ = 445 nm [43], and V abs = 10 13 cm 3 is the absorption volume of the incident laser light. This allowed us to have n i varying up to a factor of 50, by controlling the repetition rate from f rep = 5 MHz (n 5 = 5 × 10 15 cm −3 ) down to f rep = 100 kHz (n 0.1 = 2.5 × 10 17 cm −3 ). We present in Figure 5a (550-550 nm) and b (568 ± 5 nm) examples of decay curves obtained for different n i . We see that while close to the band-edge, as n i is increased, the stretched exponential component almost disappears, at 568 nm it decreases drastically by more than one order of magnitude but still exists, and the decays seem to have the same shape between 500 and 100 kHz but differ slightly from that of 5 MHz. The behavior close to the band-edge implies a saturation of the Förster diffusion mechanism at those wavelengths, which is the signature of defects filled by the laser excitation at n 0.1 . The second behavior at 568 nm points to a value of β that can vary. To characterize those two behaviors, we have performed FLIM experiments over Zone 2, at different value of n i . We first focus on the potential variation of the stretching exponent β in Figure 6. A pixel in Figure 6a-d represents the value of β as a colorscale, obtained from fitting the decay curves of every pixel with Equation (2), for n 5 , n 1 , n 0.5 and n 0.1 , respectively. From those figures, one can acknowledge the spatial dependence of β at the micrometer scale. We see that at 5 MHz, in Figure 6a, β has values larger than for the other smaller repetition rates, backing the qualitative assumption made for Figure 5. Its value varies from β = 0.2 up to β = 0.5 for a few points, but has the majority of the image close to β∼0.3, while for the other images, it seems to only vary with an amplitude of ±10%. Interestingly, the few blue regions (β∼0.2) in Figure 6a, which represent the minimum value for the images at 5 MHz, become red regions as the repetition frequency is decreased (i.e., n i is increased). However, those regions do not experience a large variation in the absolute value of β, which stays around 0.2-0.24 for the four values of n i . As the traps are filled, however, the green and red regions of Figure 6a tend to converge to that same range, and they are those which experience the largest change in absolute value of β. In the framework of the Förster direct transfer model, we can conclude that at 5 MHz, i.e., n 5 , some regions behave closer to a 3D Förster direct transfer diffusion (yellow to red pixels), while the rest (large majority of the image) are more compatible to the 2D model. At larger n i (smaller repetition rates, Figure 6b,c), the images are more compatible with a truncated 2D Förster direct transfer model, where only the nearest neighbors are available for diffusion. Thus, the diffusion mechanism depends on the relative population of the traps. This behavior is confirmed in Figure 6e, which presents the average value of β for the whole image at a particular n i , exhibiting that β = 0.32 at 5 MHz and falls to β = 0.23-0.24 for images recorded at 1, 0.5 and 0.1 MHz. It is worth noting that the value of β saturates at large n i .

Evolution of the Average Lifetime
We now focus on the second property observed in Figure 5: the evolution of the relative importance of the stretched exponential component compared to the fast component, as a function of n i , or the absorbed density of photons by the MAPbBr 3 thin films. First, we study the evolution of the average photon arrival time, < τ g > defined as follows [50]: where N f ast and N slow are the number of photons in the fast and slow components, Γ(1/β) is the Gamma function evaluated at 1/β, and the other parameters are obtained from the fits and introduced in Equation (2). We have computed < τ g > for the band-edge emission (550-555 nm) and the tail emission (568 ± 5 nm) over Zone 2 for n 5 , n 0.5 and n 0.1 . Results of that procedure are presented in Figure S5a-c (550-555 nm) and Figure S5d-f (568 ± 5 nm).
The average lifetime < τ g > of the emission closer to the band-edge shows little to no clear trend, but a slight acceleration at n 0.1 with < τ g > comprised ∼50-100 ps, while evolving in the range ∼100-200 ps for smaller n i . We saw, qualitatively, this behavior in Figure 5 where the stretched exponential component almost disappears from n 5 to n 0.1 for emission close to the band-edge. These results coupled to results from Figure S5a-c demonstrate that the weight of the fast component is dominating the decay at large n i . Moreover, we can see that < τ g > is spatially inhomogeneous and that the regions of longer lifetimes are the ones with smaller numbers of photons. In contrast, in the emission tail, Figure S5d-f exhibits a strong dependence of < τ g > with n i , going from dynamics ranging over < τ g >∼1-10 ns at n 5 , down to < τ g >∼0.2-0.5 ns for n 0.1 . The spatial variation of < τ g > at the emission tail of the MAPbBr 3 is anti-correlated to that of the band-edge: < τ g > is shorter in the region of lower signals. Finally, the dynamics of < τ g > over Figure S5f in the emission tail are close to that of Figure S5c near the band-edge, meaning that for n 0.1 = 2.5 × 10 17 cm −3 , the fast component dominates both decays. These results exhibit quantitatively the relative decrease in the stretched component with increasing n i and its spatial distribution.

Determination of the Density of Traps
We can go further in the analysis by focusing on the ratio We present in Figure 7 an image of 9 × 9 µm 2 , 50 × 50 pixel of Zone 2 where η is represented as a colorscale. . Figure 7d-f shows the corresponding values of η obtained for the tail if the emission (568 ± 5 nm). We can discriminate three cases for our interpretation: η 1 where the stretched dynamics (emission driven by electron diffusion) dominates, η 1 where the fast component (emission driven by excitons) and η ∼1 where both contributions are comparable. At n 5 close to the band-edge, Figure 7a, most of Zone 2 emission looks comparable, at least with the same order of magnitude contribution of both excitons and emission arising from the electron diffusion process. There is η ∼1-3 for most of the image, with a minority micron-sized region where η 1, confirming that the emission dynamics are spatially dependent at the micrometer scale. In the tail of the emission, Figure 7d, the whole image presents η 1, with a spatial modulation that shows η variation over one order of magnitude. In the red regions of Figure 7d, the emission due to the electron diffusion in the traps represents up to 99% of the emitted photons, i.e., virtually no contribution from excitons. As n i is increased, we observe a clear domination of exciton emissions at the band-edge, to finally obtain η ∼30-250 depending on the region in Zone 2, Figure 7c. We thus conclude that more than 99.5% of the emission is from excitons at the band-edge for some regions of Zone 2 at n 0.1 = 2.5 × 10 17 cm −3 , while its weight was comparable to the emission due to Förster transfer at n 5 = 5 × 10 15 cm −3 . The Förster direct transfer emission has thus been reduced by almost three orders of magnitude, relative to the exciton emission. This clearly points, coupled with the results obtained from CW experiments, to a saturation of the trap population close to the band-edge for n 0.1 . The trend is similar at 568 nm and n 0.1 , Figure 7e,f, where η∼10 in some regions, and excitons' emission represents up to 90% of the emitted photons, while it was around 1% for n 5 . We can thus be more specific combining the time-resolved and CW experiments and give an estimation of the population of traps in MAPbBr 3 thin films prepared via the anti-solvent method at n∼10 17 cm −3 , which is two orders of magnitude larger than for a single crystal of MAPbBr 3 [16].

Conclusions
To conclude, hyperspectral microscopy over 5 decades of excitation intensities performed on a MAPbBr 3 thin film at T = 4.7 K revealed the spatial inhomogeneity of the emission wavelength, with micrometer-scale regions dominated by excitons (band-edge emission) or defects (red emission tail). Two-color, time-resolved experiments revealed a complex dynamic composed by the exciton (fast mono-exponential process ∼100-200 ps) and trap emissions (stretched exponential). The stretched exponential corresponds to a resonant Förster transfer in two dimensions in a randomly distributed ensemble of traps, with the electron diffusing through the traps and recombining with immobile acceptors. Varying the density of absorbed photons n i by the MAPbBr 3 thin film, we find that the diffusion process shifts from a model where all the traps in the distribution are accessible for diffusion (low n i ) to a situation better described when only the nearest neighbors are available for the electron diffusion process (high n i ). We identified the spatial dependence of the stretching parameter β determining the dimensionality of the diffusion process. Tracking the evolution of the average emission lifetime < τ g > as a function of n i revealed its spatial inhomogeneities at the micrometer scale. Finally, we found that as n i was increased, the relative number of photons in the fast component dramatically increased to reach more than 99.5% of the emission close to the band-edge at n 0.1 = 2.5 × 10 17 cm −3 , indicating a saturation of the emission from traps. This enabled us to estimate the concentration of available diffusing traps to n∼10 17 cm −3 in MAPbBr 3 thin films prepared via the anti-solvent method.