{'fraction_non_alphanumeric': 0.035637549084981845, 'fraction_numerical': 0.025709416907460918, 'mean_word_length': 4.464777327935223, 'pattern_counts': {'":': 0, '<': 0, '<?xml version=': 0, '>': 0, 'https://': 0, 'lorem ipsum': 0, 'www.': 0, 'xml': 0}, 'pii_count': 0, 'substrings_counts': 0, 'word_list_counts': {'cursed_substrings.json': 0, 'profanity_word_list.json': 0, 'sexual_word_list.json': 0, 'zh_pornsignals.json': 0}}
Nonlinear Two-Photon Absorption in the Near-Infrared Band for Lead Bromide Perovskite Films Using an F-Scan Nonlinear Spectrometer

An F-scan nonlinear spectrometer is used to measure the two-photon absorption coefficient for CH3NH3PbBr3 perovskite films from 690 to 995 nm. This spectrometer uses an electrically focused tunable lens and a tunable femtosecond-pulse laser (Mai Tai-HP) with a resolution of 5 nm. Two-band models and saturation irradiance corrections are used to fit the experimental data. The nonlinear absorption in this wavelength range is of the order of cm/MW. We found that the best agreement between the experimental data, the reported values in the literature, and the theoretical model is obtained for an extended two-band model with irradiance saturation correction.


INTRODUCTION
Lead halide perovskites have emerged as desirable materials due to their particular optical and transport properties, 1,2 leading to the fabrication of high-performance solar cells and, more recently, light-emitting diodes (LEDs) 3,4 and batteries. 5 The typical structure of a three-dimensional (3D) perovskite is ABX 3 , being A an organic molecule such as methylammonium CH 3 NH 3 (MA), or formamidinium (FA), or an inorganic cation such as cesium (Cs). Lead (Pb) is commonly used in the B position, while I, Br, and Cl are used as X site anions. Combinations of these different atoms have been used to tune the bandgap and material stability. 6 Metal halide perovskites have received attention due to significant advances in the increment of efficiency in photovoltaic devices, 7 the defecttolerant characteristic produced onto the material structure, and remarkable optoelectronic properties such as high absorption coefficient and high photoluminescence in the linear regime.
On the other hand, nonlinear optical properties have been less studied, but reports show large third-order optical nonlinearities on thin-film halide perovskites; thus, there is enormous viability for applications in photonic devices. 8−12 In particular, Zhou et al. 8 listed the measured values of twophoton absorption (TPA) for MAPbBr 3 perovskites found in the literature, with values in the order of cm/GW for single crystals and in the order of cm/MW for films. The nonlinear optical measurements in these reports were made for a specific laser wavelength, yet it is desirable to measure these properties over a broader range of wavelengths to analyze some of the mechanisms that contribute to the nonlinear properties. In this work, we have been able to measure the nonlinear two-photon absorption (TPA) β of CH 3 NH 3 PbBr 3 thin-film perovskite in the near-infrared band, from 690 to 995 nm with a step of 5 nm using a recently proposed nonlinear absorption spectrometer (NAS) reported in ref 13, which is based on the TFscan technique. 14,15 The first part of Section 2 explains the technique used to fabricate the samples and present the physical and linear optical properties. In the second part of Section 2, we explain the working principle of the nonlinear spectrometer. In Section 3, first, we present the experimental results for the nonlinear two-photon absorption, and, second, we fit the data using two-band models. Finally, we analyze the results and give conclusions.

Fabrication and Characterization Methods of
Perovskite Films. For this work, films were prepared according to the method reported in ref 5. Methylammonium bromide (MABr) from dyesol was used as an organic cation. Dimethyl sulfoxide (DMSO, Sigma-Aldrich) and N,Ndimethylformamide (DMF, Alfa Aesar) were used as solvents, with lead bromide (PbBr 2 ) (Alfa Aesar) as the lead source. To obtain MAPbBr 3 , the precursor solution was deposited on top of 2.5 cm × 2.5 cm glass substrates by spin coating at 4000 rpm for 25 s. After 10 s of spin coating, 500 μL of diethyl ether was dropped to remove the DMF quickly. The films were then annealed at 65°C for 1 min plus 100°C for 10 min. Absorption measurements were taken in a Cary 100 Varian Xray diffractogram in a PANalytical diffractometer. Samples were scanned from 2θ = 5 to 60°in a Bragg-Brentano geometry, using Cu Kα (1.5408 Å) radiation with a step size of 0.04°and a speed of 5 degree/min. Figure 1a shows the absorption spectrum of the MAPbBr 3 film; the particular excitonic characteristic of the Br perovskite is observed near 521 nm. The absorption values near high wavelengths can be attributed to scattering due to the film's roughness. The Tauc plot from the absorption data is presented in the inset. By extrapolating the linear region of the absorption coefficient, the electronic bandgap energy is estimated to be 2.31 ± 0.01 eV. Figure 1b shows the characteristic signal of MAPbBr 3 perovskite in its cubic phase. A scanning electron microscopy (SEM) image of a MAPbBr 3 film fabricated with the same technique can be observed in Figure 1c in ref 5, showing a smooth morphology with faint grain boundaries.

Nonlinear Optical Spectrometry in the Near-Infrared Region.
A simplified scheme of an F-scan nonlinear spectrometer 13 is depicted in Figure 2. The laser pulse is modulated with an optical chopper and filtered with a computer-controlled neutral filter (NF). The light then impinges onto an electrically focused tunable lens (EFTL) that focalizes and defocalizes the beam on the thin-film sample of study. The sample is placed at a fixed position d s inside the range of the EFTL. All light transmitted through the sample, including scattered light, is collected and detected with a lens (L1) and an integrating sphere (D1). This light is then used to determine the two-photon absorption (TPA) coefficient.
For the experimental implementation, we used a Mai Tai tunable laser with a repetition rate of 80 MHz, an average pulse width of 75 fs, and laser emission from 690 to 995 nm, with a resolution of 5 nm. Respecting the use of high repetition lasers in measurements of TPA, several studies have addressed concerns about the correctness of the measurement, 16,17 particularly regarding the existence of cumulative thermal   Tai spectra-physics tunable pulsed laser, NF is the neutral density filter, OC is the optical chopper, EFTL is the electrically focused tunable lens, and M1 is a mirror. The inset shows the experimental parameters used in the theoretical calculation for TPA. effects due to linear absorption. 16,18−21 To address this problem, blanking has been suggested. 18,19,22,23 Blanking allows the material to transfer the heat before a new train of pulses arrives by blocking the laser beam periodically. Thus, because we are already using a 50% duty cycle optical chopper to modulate light for the lock-in amplifier, as stated by Nag et al., 23 to avoid thermal effects the chopper frequency has to be greater than 4D/2w 0 2 , where D is the thermal diffusion coefficient and w 0 is the beam waist at the sample. In our experiment, w 0 ≈ 0.016 mm, D = 0.145 mm 2 /s, 24 and the frequency used is approximately 1130 Hz. Possible effects regarding light polarization 19,25 were not addressed in this work. At the entrance surface of the sample, the average power was 20 mW and the irradiance varied as shown in Figure 3.
The EFTL is an OPTOTUNE-1030, controlled by an OPTOTUNE lens driver that gives a maximum current of 300 mA with a resolution of 0.1 mA. We used an integrating sphere to measure the transmitted laser light. The light collected by the integrating sphere is sent to a STANFORD RESEARCH 830 dual-channel lock-in amplifier, controlled through a general purpose interface bus (GPIB).
The nonlinear absorption spectrometer uses an electronically focused tunable lens (EFTL), with tunable focal distance f, as a dynamic system to generate, for each wavelength, the corresponding transmittance signal of an F-scan setup at two particular points: the point of maximum absorption (z = d s − f = 0), which corresponds to high-intensity (HI) irradiation, and the point of low-intensity (LI) irradiation (z = d s − f → ∞), the region where nonlinear effects are negligible. Using the conventional Z-scan open-aperture expression for the transmittance, 26 we can obtain an expression for the transmittance of light through the sample as a function of z = d s − f: Saturation of the TPA coefficient has been observed due to irradiance saturation. 27−29 To take into account this effect, we use eq 3 and the saturation value I sat = 0.57 GW/cm 2 reported by Liu et al. 27 ( ) where β 0 is the nonsaturated TPA coefficient.

Experimental Nonlinear Two-Photon Absorption.
The experimental TPA coefficients obtained with eq 2 correspond to the red dots in Figure 4. If we consider the correction proposed by Liu et al. 27 for a saturation irradiance of 0.57 GW/cm 2 , using eq 3, the corrected TPA coefficient corresponds to the black dots in Figure 4. As shown in Figure  4, the correction scales the value of TPA by an approximately constant factor. Using the common uncertainty 30 measured for CdS in a previous study, 13 for both curves, the estimated uncertainty is ±0.4 cm/MW; the maximum value of TPA is around two-thirds of the energy gap, as expected.

Nonlinear Two-Photon Absorption Model.
For MAPbI 3 , Wei et al. 31 report experimental TPA measurements without the characteristic maximum commonly observed in zinc blended semiconductors crystals at 0.7 times the bandgap energy of the material. Thus, they included relativistic effects induced by spin−orbit coupling to model this lack of a maximum. Nevertheless, by inspection of Figure 4, it is clear that for our MAPbBr 3 films, a maximum value of TPA exists around 0.7 times the bandgap energy. For this reason, we used two models based on the two-band model (TBM) presented in   where A is the scaling factor for β, E g is the bandgap energy in eV, ℏω is the photon energy in eV, E p is the Kane momentum energy (approximately 20 eV for direct bandgap materials), and C is the Sommerfeld factor and E b is the exciton binding energy in eV.
where n = 0, 1, or 2 correspond to a−a, a−f (or f−a), and f−f transitions, respectively. f n ≤ 1 is a numerical factor that depends on the angular averages of the interband matrix elements that are going to be used here as fitting parameters. Using the nonlinear least square method (python scipy.optimized.curve_fit), we have fitted the experimental data restricting the value of the bandgap energy to onestandard deviation (2.30−2.32 eV) and allowing the scaling factor to be as much as 40% below the accepted range (between 1940 and 5200). 36 If we consider the irradiance saturation of 0.57 GW/cm 2 and use an exciton binding energy of 20 meV, 37 the TBEM model will predict the blue dashdotted curve of Figure 6 for parameters A = 1880 cm/MW and E g = 2.30 eV. The scaling factor A has the expected numerical value but in cm/MW. To evaluate the quality of the fit, we calculated the reduced chi-squared statistic, finding a χ ν 2 = 1.6, with a probability of finding a better fit of more than 99%, corresponding to a very poor fitting. In contrast, the GTBM model predicts the red curve with reduced chi-square statistic of χ ν 2 = 0.6, corresponding to a probability of finding a better fit of only 2%; thus, it is an excellent fit. The fitting parameters are A = 2360 cm/MW, E g = 2.30 eV, f 1 = 0.13, f 2 = 0.77, and f 3 = 0.24, meaning that the main contribution is from the a−f and f−a transitions, while the a−a and f−f transitions contribute less. Again, the scaling factor A has the expected numerical value but in cm/MW. Now, when ignoring any possible irradiance saturation, the fitting result can be seen in Figure 7. Using an exciton binding energy of 20 meV, 37 the TBEM model will predict the blue dash-dotted curve for parameters A = 1182 cm/MW and E g = 2.31 eV. The scaling factor A is below the expected numerical value. The reduced chi-squared statistic is χ ν 2 = 1.7, with a probability of finding a better fit of more than 99%, which corresponds to a very poor fit. The GTBM model predicts the red curve with a reduced chi-square statistic of χ ν 2 = 1.0, corresponding to a probability of finding a better fit of only 62%; thus, it is a poor fit. For this case, the fitting parameters are A = 2100 cm/MW, E g = 2.30 eV, f 1 = 0.10, f 2 = 0.22, and f 3 = 0.93, meaning that the main contribution is from the f−f transition.
On behalf of the scattering of the data around the fitting curves, the smaller chi-square statistic, and the fact that an a− f(f−a) transition is more plausible than an f−f transition, we believe that the correct experimental values correspond to  We find that measuring the statistical distance between the values is 5.6 deviations; thus, they cannot be considered equal without a statistical doubt. Despite the disagreement with the value reported by Liu et al., it is important to remember that they did not report an uncertainty; if we suppose a 15% relative error, like in our data, the statistical distance reduces to only 2.5 deviations.

CONCLUSIONS
In conclusion, we have obtained an experimental dispersion curve for the two-photon absorption coefficient of MAPbBr 3 three-dimensional hybrid perovskite in a wavelength range of 690−995 nm, using a nonlinear absorption spectrometer based on an F-scan technique. We have found that the experimental TPA values corrected for irradiance saturation give the best possible fit, agreeing with the accepted fitting, and experimental parameters. We have found that the generalized two-band model of Garcia et al. 35