Manipulating crystallization dynamics through chelating molecules for bright perovskite emitters

Molecular additives are widely utilized to minimize non-radiative recombination in metal halide perovskite emitters due to their passivation effects from chemical bonds with ionic defects. However, a general and puzzling observation that can hardly be rationalized by passivation alone is that most of the molecular additives enabling high-efficiency perovskite light-emitting diodes (PeLEDs) are chelating (multidentate) molecules, while their respective monodentate counterparts receive limited attention. Here, we reveal the largely ignored yet critical role of the chelate effect on governing crystallization dynamics of perovskite emitters and mitigating trap-mediated non-radiative losses. Specifically, we discover that the chelate effect enhances lead-additive coordination affinity, enabling the formation of thermodynamically stable intermediate phases and inhibiting halide coordination-driven perovskite nucleation. The retarded perovskite nucleation and crystal growth are key to high crystal quality and thus efficient electroluminescence. Our work elucidates the full effects of molecular additives on PeLEDs by uncovering the chelate effect as an important feature within perovskite crystallization. As such, we open new prospects for the rationalized screening of highly effective molecular additives.

I ntroducing molecular additives into perovskite precursors has become one of the most effective and prevailing strategies to improve the performance of metal halide perovskite lightemitting diodes (PeLEDs), and recently has boosted the external quantum efficiency (EQE) to high values above~20% [1][2][3][4][5][6][7] . The performance enhancement results from suppressed nonradiative recombination and/or improved light-out coupling efficiency due to morphological variation 6,8,9 . The former is generally believed to be associated with defect passivation-surface dangling bonds in perovskites are said to be healed by additional coordination or ionic bonding with the additives 1,10,11 , leading to the elimination of trap states. Accordingly, a wide range of passivating molecular additives has been investigated for PeLEDs, especially those containing Lewis base moieties such as amino-and carboxylfunctionalized molecules 8,[12][13][14][15] . In addition, the critical factors underlying passivation effectiveness have been well identified, including hardness/softness of functional moieties 16 , steric hindrance determined by molecular configurations 17,18 , as well as hydrogen bonding between passivating functional groups and organic cations 7,19 .
Despite these advances in understanding molecular passivation, a widely observed phenomenon yet to be rationalized by passivation alone is that almost all the effective additives for perovskite emitters are chelating molecules that contain more than one electron-rich functional moiety 6,7,20 . Conversely, monofunctionalized counterparts have received far less attention despite exhibiting similar molecular configurations and identical passivating groups. This contrast implies additional factors are at play (beyond the influence of passivation) that critically affect radiative recombination in the resulting halide perovskites.
Inspirations can be drawn from recent work on perovskite photovoltaics and point toward the possibility that molecular additives can modify crystallization dynamics 21,22 , resulting in discrepancies in defect densities and hence different radiative recombination processes. For instance, modulating crystallization dynamics with Lewis basic solvents, such as dimethyl sulfoxide (DMSO) and 1-methyl-2-pyrrolidinone (NMP), have been intensively investigated [21][22][23] . However, the chemical composition of perovskite emitters is remarkably different from that of perovskite photovoltaics due to the required enhancement of quantum and/or dielectric confinement 4,[24][25][26] . Consequently, the crystallization control strategies developed in perovskite photovoltaics cannot be expected to be straightforwardly translated to PeLEDs. A thorough understanding of how the molecular additives influence perovskite crystallization in PeLEDs is thus required to better rationalize where the optoelectronic improvements originate.
Here, we reveal that the effects of prevailing molecular additives on manipulating the crystallization process of perovskite emitters are largely underestimated. By comparing two highly effective chelating molecules with their respective monofunctionalized counterparts, we find that chelation enhanced lead-additive binding affinity enables the formation of thermodynamically stable intermediates, leading to effective inhibition of perovskite nucleation and retarded crystal growth. The retarded crystallization process results in highly emissive perovskite thin films and high-efficiency near-infrared (NIR) PeLEDs with EQEs over 20%. We further generalize our conclusions by investigating a wide range of additives with different coordination numbers and functional moieties. Our results fill the knowledge gap between chelating molecular additives and bright perovskite emitters, through establishing the correlation between leadadditive coordination affinity and perovskite crystallization dynamics. As such, we provide a universal guideline for tuning the crystallization dynamics for high-performance PeLEDs by leveraging the chelate effect.

Results
PeLED device performance and film characterization. We compare the effects of two popular chelating additives (CAs) with their respective mono-functionalized additives (MFAs) on formamidinium lead tri-iodide (FAPbI 3 ) emitters. The CAs we use are 5-aminovaleric acid (5AVA) and 3,6,9,12-tetraoxatetradecane-1,14-diamine (NH 2 -PEG 4 -NH 2 ), which consist of commonly used amino-and/or carboxyl-groups; their respective mono-functionalized counterparts are n-butylamine (BA) and propionic acid (PA) for the former, and 2-(2-methoxyethoxy)ethanamine (m-PEG 2 -NH 2 ) for the latter (Fig. 1a). Here, m-PEG 2 -NH 2 is regarded as an MFA despite the electrondonating nature of the O atoms. This is based on the fact that the effects of oxygen atoms on crystallization are not as prominent as amino groups, as verified by comparing a counterpart with a longer polyethylene glycol (PEG) chain (2,5,8,11-tetraoxatridecan-13-amine, m-PEG 4 -NH 2 ) (Supplementary Fig. 1 and Supplementary Note 1). All the additives were introduced into perovskite precursors at a mole ratio of PbI 2 :formamidinium iodide (FAI):additive = 1:2:x (x = 0-0.8). The use of excess organic components aims to passivate the vacancy defects, as widely demonstrated in PeLEDs 27,28 . As shown in the X-ray diffraction (XRD) patterns ( Supplementary Fig. 2), the addition of additives does not change the three-dimensional structure of FAPbI 3 at the optimal feed ratio for device fabrication. Besides, all the perovskite films show nano-island morphological features, although the addition of CA leads to smaller crystal sizes and denser surface coverage compared to the MFA-based counterparts ( Supplementary Fig. 3). Our further investigations indicate that the morphological differences between MFA-and CA-films cannot be straightforwardly linked to the discrepancies in the device performance (Supplementary Fig. 4 and Supplementary Note 2).
We fabricate PeLEDs with additives and that without (named as the "control") by subsequential deposition of zinc oxide (ZnO; 10 nm)/polyethyleneimine ethoxylated (PEIE)/FAPbI 3 /poly(9,9-dioctylfluorene-co-N-(4-(3-methylpropyl)) diphenylamine (TFB)/ MoO 3 (7 nm)/Au (80 nm) on indium tin oxide (ITO) substrates (see Fig. 1b and a cross-sectional image in Supplementary Fig. 5a). Devices based on the combined use of BA and PA at the mole ratio of 1:1 (named as the "BA-PA") are also studied for comparison. The energy levels of the perovskite films are determined by ultraviolet photoelectron spectroscopy (UPS) spectra and a flat-band energy level diagram of devices (with control, m-PEG 2 -NH 2 and NH 2 -PEG 4 -NH 2 devices as the examples) are shown in Supplementary Fig. 5b-d and further discussed in Supplementary Note 3.
As the additives reduce the charge carrier injection to some extent due to their insulating nature ( Supplementary Fig. 6), we carefully control their loading contents (see the details in Supplementary Table 1) to maximize the peak EQEs and simultaneously ensure a radiance as large as possible. We show a summary of average EQE values and detailed performance parameters with the optimized feed ratio in Fig. 1c and Supplementary Table 2, respectively, from which we find that CA addition results in considerable enhancement in terms of both radiance and EQE values compared to control and the respective MFA counterparts. Specifically, control devices show an average peak EQE of 5.3 ± 0.7%; 5AVA-devices give an average peak EQE of 15.5 ± 0.8% compared to 7.1 ± 0.7%, 6.5 ± 0.6%, and 7.6 ± 0.5% of BA, PA, and BA-PA mixture, respectively; NH 2 -PEG 4 -NH 2 devices have an average peak EQE of 19.2 ± 0.6% compared to 11.0 ± 0.5% of m-PEG 2 -NH 2 . The NH 2 -PEG 4 -NH 2 -and 5AVA-devices exhibit average radiance values of 360 ± 27 W sr −1 m −2 and 173 ± 20 W sr −1 m −2 , which contrast sharply with the respective MFA devices (Supplementary Table 2). In addition, the champion NH 2 -PEG 4 -NH 2 -and 5AVA-devices reach high EQE values of 20.2% and 17.2%, respectively (see device characteristics in Supplementary Figs. 7 and 8). Intriguingly, BA-PA mixed precursors and double amounts of m-PEG 2 -NH 2 do not give rise to comparable performance as the respective CA counterpart, suggesting the presence of additional underlying factors determining the device performance in addition to the amounts of Lewis base moieties. Moreover, CA devices show a prolonged operational lifetime with respect to the MFA counterparts measured under a constant current density of 20 mA cm −2 ( Supplementary Fig. 9).
Given that m-PEG 2 -NH 2 shares identical moieties with NH 2 -PEG 4 -NH 2 , we choose them as the main candidates for the rest of our investigations. We explore the photophysical properties of the perovskite films to further evaluate the radiative recombination and provide a rationale for the superior performance of CA devices. As shown in Fig. 1d, the NH 2 -PEG 4 -NH 2 films display a much less pronounced dependence of photoluminescence quantum efficiencies (PLQEs) on excitation fluences with respect to its MFA counterparts and control films, giving the highest peak PLQE value up to 50%. More efficient radiative recombination in NH 2 -PEG 4 -NH 2 films is also evident in PL mapping measurements 29,30 , from which we observe very homogeneous emission in space for the NH 2 -PEG 4 -NH 2 film with time but distinct PL blinking in the control and m-PEG 2 -NH 2 cases (Supplementary Fig. 10 and Supplementary Movies 1-3).
To further understand the discrepancy in the PL properties, we proceed to investigate the electronic states of the perovskite films using transient absorption (TA) measurements ( Fig. 1e-g). In all the cases, we observe two photobleaching (PB) bands situated near 780 nm (PB1) and 490 nm (PB2), suggesting no significant alteration in their photophysical behavior with additives addition. The appearance of PB2 has been previously observed in bulk CH 3 NH 3 PbI 3 perovskites, but the assignment remains unclear 31,32 . One previous report proposes that PB1 and PB2 arise from the transitions between two valence bands (VBs) to a common conduction band minimum (CBM) 31 . Notably, the control film shows a distinct band (>800 nm) in different delayed time scales, suggesting a relatively large number of sub-bandgap trap states 33 . The band tail absorption is reduced to some extent in m-PEG 2 -NH 2 films but remarkably weakened in NH 2 -PEG 4 -NH 2 incorporated films. Following the decreases in the band tail absorption, the PB1 is slightly blue-shifted compared to control and m-PEG 2 -NH 2 cases, which can be assigned to the decoupling between PB1 and the less pronounced sub-bandgap trap states. Combining the results of PLQEs, PL mapping, and TA measurements, we conclude that trap-assisted nonradiative recombination is largely suppressed in NH 2 -PEG 4 -NH 2 films, resulting in superior device performance. Understanding the effects of additives. The reduced trapassisted nonradiative recombination could be attributed to efficient defect passivation, improved crystal quality, or a collective effect of them. By Fourier transform-infrared spectroscopy (FT-IR) measurements, we demonstrate that both m-PEG 2 -NH 2 and NH 2 -PEG 4 -NH 2 are retained in the perovskite films, as confirmed by the appearance of stretching vibration of ether bonds (C-O-C) at around 1110 cm −1 ( Supplementary Fig. 11). The presence of additives provides opportunities for them to interact with defects chemically, contributing to the mitigated nonradiative recombination. It is further confirmed by the density of functional theory (DFT) calculations (Supplementary Figs. 12 and 13). Notably, our theoretical calculations suggest little differences in passivation effectiveness between these two additives, as further discussed in Supplementary Fig. 13 and Supplementary Note 4. This is also consistent with the fact that the passivation effects are usually determined by the type of functional groups and their chemical environment, both of which are identical in m-PEG 2 -NH 2 and NH 2 -PEG 4 -NH 2 , as confirmed by calculation of the electronic property of the molecules ( Supplementary  Fig. 14). These results exclude the possibility that the significant difference in radiative recombination and device performance are the results of different passivation effectiveness between these two molecules.
With little difference in passivation effectiveness expected to occur between the use of CA and MFA additives, we focus our investigations on the effects of different molecules on the crystallization process. The first question that arises is whether these additives have different interaction strengths with the perovskite precursors; the interaction strength determines the intermediate phases and thus plays the key role in governing the crystallization process 21,34,35 . Amines can form hydrogen bonds with charged ammonium cations (FA + ) 7 , which are characterized by 1 H NMR. The number of Lewis base moieties is constant in all the samples, that is, at a mole ratio of 1:0.2 for FAI: m-PEG 2 -NH 2 mixtures and 1:0.1 for FAI:NH 2 -PEG 4 -NH 2 ones. The resonance signal of protonated ammonium peaks at 9.3 ppm in the pure FAI deuterated-DMF solution (Fig. 2a). In both cases with additive addition, the resonance signals of ammonium shift to 8.8 ppm, implying similar interaction strengths between both additives and FA + . These results are in line with the DFT calculations that the chemical environments of the amino groups and oxygen atoms in m-PEG 2 -NH 2 and NH 2 -PEG 4 -NH 2 are very close. In contrast, for the samples containing PbI 2 , the resonance of ammonium becomes less susceptible to NH 2 -PEG 4 -NH 2 addition compared to m-PEG 2 -NH 2 , as adding the latter results in a much larger shift toward high field (Fig. 2b). These results imply that NH 2 -PEG 4 -NH 2 has a stronger tendency than m-PEG 2 -NH 2 to work with Pb 2+ and thus leads to a weaker interaction with FA + in the Pb 2+ -contained precursor solutions.
In this regard, we proceed to directly investigate the interaction between the additives and Pb 2+ . Adding NH 2 -PEG 4 -NH 2 into PbI 2 solution leads to the formation of white precipitation ( Supplementary Fig. 15). This observation is consistent with the fact that amines are well-known building blocks for constructing metal amide complexes. Moreover, the yielded products will gradually diminish with FAI addition, indicating a competitive ligand exchange between the electron-rich moieties of additives and iodide anions from FAI.
We next quantified the interaction strength between Lewis base ligands and Pb 2+ on the basis that iodide ions compete with Lewis base moieties to coordinate with Pb 2+ ions 36 . We perform Benesi−Hildebrand (BH) analyses and determine the formation constant (K f ) of iodoplumbate complex PbI 3 − in the solution of interest 37 , by monitoring the Ultraviolet-visible (UV-vis) absorption evolution with increasing FAI content in the PbI 2 :additive parental solutions. We double the amounts of m-PEG 2 -NH 2 with regards to NH 2 -PEG 4 -NH 2 to keep Lewis base moieties constant in the solutions. By following the absorption spectra evolution with increasing FAI content, we observe gradually enhanced absorption intensity of PbI 3 − in both cases ( Supplementary  Fig. 16). Notably, PbI 2 :NH 2 -PEG 4 -NH 2 parental solution requires much higher amounts of FAI to achieve comparable absorbance to m-PEG 2 -NH 2 ones. Figure 2c shows the BH plots, from which our calculations yield a larger K f value (1/slope) of 31.1 M −1 for m-PEG 2 -NH 2 solution compared to 3.1 M −1 for that with NH 2 -PEG 4 -NH 2 . These results confirm that the iodoplumbate complex is more difficult to form with the presence of NH 2 -PEG 4 -NH 2 , indicating a higher affinity between NH 2 -PEG 4 -NH 2 and Pb 2+ . We attribute this to the chelate effect that the chelating ligands (e.g. NH 2 -PEG 4 -NH 2 ) have a stronger affinity to the metal ions than a collection of their nonchelating or weak chelating counterparts (e.g. m-PEG 2 -NH 2 ) despite identical functional groups, leading to kinetically more stable lead-additive coordination in the solutions.
Another important effect of the Pb 2+ -CA coordination with chelation enhanced affinity is the formation of thermodynamically stable lead complexes in the solid states. We show the UV-vis absorption spectra of PbI 2 : additive and neat PbI 2 films in Fig. 2d and their XRD patterns in Supplementary Fig. 17. All these films were annealed at 100°C for 10 min before measurements. Compared to the neat PbI 2 film, the PbI 2 : NH 2 -PEG 4 -NH 2 film turns to colorless ( Supplementary Fig. 18), where the absorption onset shifts to a shorter wavelength of around~450 nm (Fig. 2d). It also displays an amorphous feature with no discernible diffraction peaks, in sharp contrast to the PbI 2 neat film which shows a typical crystalline (110) diffraction peak at 2θ = 12.6° (Supplementary Fig. 17). These results suggest that the addition of NH 2 -PEG 4 -NH 2 reconstructs Pb 2+ coordination to form complexes that are stable at high temperatures for perovskite crystallization. In contrast, the PbI 2 :m-PEG 2 -NH 2 film shows no difference compared to the neat PbI 2 film regarding both absorption and XRD patterns, indicating that m-PEG 2 -NH 2 only has labile coordination to lead cations and that the formed complexes are readily decomposed at high temperatures needed for annealing perovskite films. We observe similar results in 5AVA and its MFA counterparts (BA, PA, and BA-PA combination). In particular, 5AVA:PbI 2 films show identical blue-shifted absorption spectra and amorphous features as the NH 2 -PEG 4 -NH 2 cases, while BA, PA, and BA-PA additions give no difference compared to neat PbI 2 films ( Supplementary  Fig. 19).
The formation of kinetically and thermodynamically more stable complexes, both attributed to the strong affinity of CAs to Pb 2+ , collectively affect the crystallization process of perovskites. Our characterizations show that the additives which exhibit a stronger affinity to lead cations are more effective in suppressing perovskite nucleation. The XRD patterns for the control and m-PEG 2 -NH 2 (PbI 2 :FAI:m-PEG 2 -NH 2 = 1:2:0.2) precursor films before annealing display the diffraction peaks from α-phase FAPbI 3 ( Supplementary Fig. 20a, b). We observe broad PL emission bands ( Supplementary Fig. 21) peaking at shorter wavelengths compared to the typical 3D FAPbI 3 in the control and m-PEG 2 -NH 2 -based precursor films, as well as two distinct excitonic absorption features peaking at~450 and~495 nm (Fig. 3a). Although the assignment of~450 nm excitonic absorption is not clear, the other results clearly suggest the formation of α-FAPbI 3 nanocrystals with size smaller than the Bohr radius. This can be directly observed by the transmittance  (Fig. 3b, c), where (100) and (200) planes of FAPbI 3 crystals are discernable in control-and MFA-precursor samples. It is worthy to mention that the perovskite nanocrystals form during the film formation process rather than in the precursor solutions, as no "Tyndall effect" is visible in the precursor solutions ( Supplementary Fig. 22). In comparison, no signature of crystalline perovskites can be found from XRD patterns (Supplementary Fig. 20c) and TEM images (Fig. 3d) in the optimal NH 2 -PEG 4 -NH 2 precursor films without annealing (PbI 2 :FAI:NH 2 -PEG 4 -NH 2 = 1:2:0.2). Consistently, we do not observe any discernible excitonic absorption features from UV-vis absorption spectra (Fig. 3a). Due to the absence of perovskites in the precursor films, no PL can be detected. All these results indicate that the perovskite nucleation is completely suppressed at room temperature with the addition of NH 2 -PEG 4 -NH 2 by forming kinetically and thermodynamically stable intermediate phases. Notably, although these results indicate that chelating ligands are effective in suppressing perovskite nucleation, we notice that suppression of nucleation at room temperature does not necessarily indicate a good film quality, as simply increasing the amount of nonchelating additive can also suppress perovskite nucleation yet without benefiting the film quality ( Supplementary Fig. 23). We further investigate the in situ crystallization process with different additives by monitoring the grazing-incidence wideangle X-ray scattering (GIWAXS) pattern evolution during thermal annealing. Since the conversion from the intermediate phases to α-FAPbI 3 perovskites is very fast at a high temperature (100°C) for device fabrication, we slow down the process by continuously increasing the temperature from 88 to 100°C in 60 s. We show the evolution of integrated diffraction intensities on annealing time in Fig. 3e, f. The m-PEG 2 -NH 2 film (0.4 eq to lead content) initially shows two distinct scattering rings with q value of~0.98 and~0.83 Å −1 , which correspond to α-FAPbI 3 and nonperovskite δ-FAPbI 3 , respectively. The signatures of α-FAPbI 3 and δ-FAPbI 3 can be detected in the NH 2 -PEG 4 -NH 2 sample at the initial stage as well, but with much weaker diffraction intensities, consistent with the XRD results that the perovskite nucleation can hardly occur at a low temperature ( Supplementary  Fig. 20c). In addition, five more weak diffraction signals at the q values of 0.40, 0.65, 0.74, 0.86, and 0.90 Å are visible, implying the coexistence of multiple crystalline intermediate phases during crystallization. Following the increase of annealing time, the diffraction signals from δ-FAPbI 3 and intermediate phases gradually diminish until they eventually disappear, accompanied by the increases in α-FAPbI 3 diffraction intensity. Notably, a much longer time is required for the saturation of (110) plane diffraction peak in the NH 2 -PEG 4 -NH 2 film compared to m-PEG 2 -NH 2 . This is in line with the proposed competitive ligand exchange process that iodide anions need to replace Lewis base ligands to coordinate with lead cations, and thus higher energy is required for the formation of PbI 6 4− octahedra and following perovskite nucleation 38 . In addition, we suggest that the ligand exchange process continually occurs during the crystallization process, leading to constrained crystal growth and hence small grain sizes ( Supplementary Fig. 3). As such, we conclude that the chelating additives serve as more efficient crystallization inhibitors by leveraging their strong affinity to lead cations.
Investigations of other Lewis base additives. In order to complete the whole picture and generalize our findings, we investigate a wide range of Lewis base additives. We start with DMSO, the most intensively used Lewis base additive in perovskite photovoltaics, where a colorless DMSO-PbI 2 -FAI intermediate phase is expected to form 21 . It is believed that the formation of the intermediate phase can retard perovskite crystallization and hence enable high film qualities. In contrast, we observe the appearance of red DMSO-precursor films in our case. The excitonic absorption features in the UV-vis absorption spectrum ( Supplementary  Fig. 24a) are identical to the films without any additives, indicating the formation of perovskite nanocrystals in the precursor films rather than DMSO-PbI 2 -FAI intermediate phases. The difference between the DMSO-precursor film in photovoltaics and that in LEDs in our case can be ascribed to the employment of excess FAI and relevant changes in coordination equilibrium in the films. In other words, as a result of the weak DMSO-PbI 2 binding affinity, the DMSO-PbI 2 -FAI intermediate phase is not stable in a haliderich environment that is typically used in PeLEDs. Consistently, we observe no positive effect on the device performance of PeLEDs with the DMSO addition ( Supplementary Fig. 24).
We proceed to investigate other Lewis base additives with different donating functional moieties and coordination numbers, including acrylates in Group A, carboxylic acids in Group B and aromatic amines in Group C (Fig. 4a). To estimate the dependence of lead-additive affinity on the coordination number of the additives in each group, we measure the UV-vis absorption spectra of the precursor films and keep the total amounts of functional moieties constant in each group by modifying the ratio of additives. With increasing coordination numbers, the absorption signatures from high-order iodoplumbate complex and/or crystalline perovskites gradually decrease (Fig. 4b), indicating enhanced lead-additive binding affinity. Due to the weak interaction between acrylates and lead cations, only those with high coordination numbers (e.g., tetra-and penta-acrylates) can influence the constituents of iodoplumbate complexes within the precursor films. We show the representative device characteristics in Supplementary Fig. 25. The averaged peak EQE values and all the detailed performance parameters are summarized in Fig. 4c and Supplementary Table 2, respectively. As expected, the device performance improves with increasing coordination numbers in each group. In addition, the dependence of device performance on the coordination numbers of additives is also visible in aromatic acids (Supplementary Fig. 26). All these results further confirm that increasing the lead-additive affinity by leveraging the chelate effect is a universal and effective strategy for preparing highly emissive perovskite thin films.
Having concluded the effectiveness and general applicability of chelation enhanced lead-additive affinity on improving film quality and device performance, the remaining question is if there is an upper limit to the affinity strength before film quality no longer benefits. We thus start with the bidentate chelating additive, hexamethylenediamine (HMDA), from which we enhance the affinity strength by moving to a multidentate chelating ligand, that is, triethylenetetramine (Trien) (see molecular structures in Supplementary Fig. 27). Trien is wellknown for its superstrong affinity to metal ions and hence has been widely used to construct environmentally stable metal complexes 39 . Adding Trien into the perovskite precursor solutions leads to the immediate formation of white precipitation ( Supplementary Fig. 28), indicating that the lead-Trien binding is too strong to be replaced by iodides. This result also indicates that an additive with too strong affinity may bring about impurities into the films. In contrast, HMDA addition has no negative impact on perovskite formation, and its binding affinity to lead cations is strong enough to suppress perovskite nucleation at room temperature by leveraging the chelate effect (Supplementary Fig. 29). The optimal HMDA-based devices give an average EQE of~11% (Supplementary Fig. 30). These results indicate that an optimal range of coordination strength should be targeted, with a superstrong lead-additive affinity negatively impacting the perovskite formation.

Discussion
Based on our collected findings, we present an overview of the correlation between lead-additive coordination affinity and crystallization dynamics in Fig. 5. In all scenarios, the additives compete with halide anions to coordinate with lead cations, resulting in an enhanced activation energy barrier (ΔE a ) for perovskite nucleation due to the additional interactions involved in the formation of PbI 6 4− octahedra. However, the halide-rich environment typically used for fabricating PeLEDs significantly hinders the coordination of additives and the relevant formation of intermediate phases. This becomes dominant for the case of using additives with weak coordination affinities, such as DMSO and the above-mentioned MFAs, resulting in little difference in the energy barrier for perovskite crystallization compared to those without using any additives. Consequently, perovskite nucleation and crystallization occur rapidly at the initial stage during annealing or even happen at room temperature, leading to a high density of defect sites and severe nonradiative recombination. In Fig. 5 Schematic illustration of perovskite crystallization process with different lead-additive binding affinity. The additives with weak coordination affinities (e.g., MFAs) cause rapid perovskite nucleation at the initial stage of thermal annealing or even at room temperature. And those with strong affinity (e.g., CAs) bring about thermodynamically stable intermediates, resulting in highly emissive perovskites. The pink and yellow boxes illustrate the competitive ligand exchange between halides and additives with different binding affinity, as well as the respective evolution of thermal energy barrier (solid line) of perovskite crystallization compared to that without using additive (dash line). Here, E a denotes the activation energy for crystallization, and ΔE a is the enhanced energy caused by additive coordination. comparison, the strong lead-CA binding affinity enables effective ligand exchange, which takes over halide coordination. It thus gives rise to an enlarged energy barrier for perovskite nucleation by the formation of more thermodynamically and kinetically stable intermediate phases, slowing down the release of lead cations for perovskite crystallization during the annealing process. We also note that a multidentate additive with superstrong affinity may suppress the formation of perovskite crystals or bring about impurities that deteriorate the device performance. As such, for a specific precursor with a particular halide stoichiometry, tuning the lead-additive coordination affinity will be very useful for achieving optimal crystal quality. From a material design point of view, besides tuning coordination number of chelating molecular discussed here, other factors influencing the stability of metal complexes such as hard-soft nature of functional groups, macrocyclic effect, resonance effect, steric hindrance are all potential protocols for tuning the affinity and relevant crystallization dynamics of perovskites.
In summary, we have elucidated how molecular additives improve electroluminescence of perovskite emitters and developed a complete crystallization model which goes beyond defect passivation, by establishing the correlation between lead-additive coordination affinity and perovskite crystallization dynamics. We conclude that the additives with a strong binding affinity that leads to kinetically and thermodynamically stable intermediate phases are more effective in inhibiting perovskite nucleation and slowing down crystallization dynamics, which enables high crystal quality and thus efficient electroluminescence. This also indicates a different requirement in lead-additive binding affinity between perovskite emitters and photovoltaic absorbers due to the discrepancies in precursor stoichiometry. Our results rationalize the general and puzzling observation that chelating additives are much more widely used in state-of-the-art PeLEDs, and provide a universal guideline for tuning the crystallization dynamics for high-performance PeLEDs by leveraging the chelate effect.
Preparation of perovskite solution. All the additives were introduced into perovskite precursors with a mole ratio of PbI 2 :formamidinium iodide (FAI):additives = 1:2:x (x = 0-0.8) in DMF to prepare the perovskite precursor solutions. The optimized feed ratio of additives for each type of device has been summarized in Supplementary Table 1.
Device fabrication. Patterned ITO glass substrates were cleaned with deionized water and ethanol step by step, followed by UV-ozone treatment for 20 min before use. ZnO nanoparticles were spin-coated on the cleaned ITO substrates with 3000 r.p.m. for 30 s in ambient. Then, the ZnO-coated substrates were moved into an N 2 -filled glove box to deposit PEIE (1.1 mg mL −1 in isopropanol) at 3000 r.p.m. for 30 s, followed by annealing at 100°C for 10 min. After cooling down, perovskite precursors were spin-coated at 3000 r.p.m. for 30 s and then annealed on a preheated hotplate at 100°C for 10 min. TFB (12 mg mL −1 in chlorobenzene) was spin-coated on top of the perovskite films at 3000 r.p.m. for 30 s. After that, 7 nm MoO 3 and 80 nm Au were deposited as the electrode in a thermal evaporator, the deposition rates are~0.2 and~0.5 Å s −1 , respectively. For the electron-only devices, phenyl-C61-butyric acid methyl ester (PCBM, 30 mg mL −1 in chlorobenzene) was spin-coated on the top of the perovskite films at 3000 r.m.p for 30 s, followed by annealing at 70°C for 10 min, followed by deposition of 100 nm Al. The active area of our PeLEDs is 7.25 mm 2 .
Device characterization. The PeLEDs were measured at room temperature in an N 2 -filled glove box. A Keithley 2400 source meter incorporated with a QE Pro spectrometer (Ocean Optics) was used to measure the electrical output of PeLEDs. The applied voltage started from 0 V and increased with a step of 0.05 V, lasting for 300 ms at each voltage step for stabilization and measurements. The operational stability of the PeLEDs was measured under a constant current density of 20 mA cm −2 in the glove box.
Grazing incident wide-angle X-ray scattering (GIWAXS). GIWAXS data were obtained from NCD-SWEET beamline, ALBA Synchrotron, Spain. The energy of the X-ray beam was set to 12.95 keV using a Si (1 1 1) channel-cut monochromator and further collimated with an array of Be lenses. The incidence angle was 1°and the diffraction patterns were collected using a Rayonix ® LX255-HS area detector, which consists of a pixel array of 5760 × 1920 (V × H) with a pixel size of 88.54 × 88.54 μm 2 for the pixel binning employed of 2 × 2. The scattering vector q was calibrated using Cr 2 O 3 as standard, obtained using a sample to detector distance of 220 mm. All the samples were measured using a calibrated Linkam THMS600 heating stage adapted for grazing-incidence experiments under N 2 atmosphere. The scan rate was one frame per two seconds.
Transient absorption (TA). The femtosecond transient absorption spectra were collected using a pump-probe configuration. A femtosecond oscillator (Mai Tai, Spectra Physics) is used as a seed laser and injected into a regenerative amplifier (Spitfire XP Pro, Spectra Physics) to generate a pulsed laser (800 nm, 80 fs pulse duration, 1 kHz repetition frequency). The second-harmonic generation (400 nm) from a BBO crystal was used as the pump. And the white light continuum (WLC) as the probe was produced by focusing the 800 nm fs pulse on a thin CaF 2 plate. Polarization between the pump and probe was set to the magic angle (54.7°).
During the experiments, gates on the pump intensity and probe spectra were set to avoid the influence of laser fluctuation on the measurements.
Thin-film characterization. X-ray diffraction (XRD) patterns were conducted by an X-ray diffractometer (Panalytical X'Pert Pro) with CuKα radiation. The X-ray was generated on a Copper target, and the wavelength of the X-ray was 1.54 Å. Morphology images of the perovskite films were measured by an FEI (Quanta 200 FEG) scanning electron microscopy under a voltage of 5 kV. UV-vis absorbance spectra were collected from a PerkinElmer Lambda 900. All the precursor films for UV-vis measurements were encapsulated by glass slides and epoxy. Steady-state PL spectra were obtained with a 450-nm excitation laser and an Andor spectrometer (Shamrock sr-303i-B, coupled to a Newton EMCCD Si array detector). The fluence-dependent PLQEs were measured by a typical three-step technique with a combination of a 445-nm continuous-wave laser, spectrometer, and an integrating sphere.
The PL mapping measurements were conducted from a home-built wide-field epifluorescence microscope based on Olympus IX71 with a ×40 objective lens (Olympus LUCPlanFL, NA = 0.6). A 405-nm diode laser (CW) was used for excitation. The power density was approximately 0.25 W cm −2 . The PL intensity as a function of time was recorded by a charge-coupled device camera (ProEM 512B, Princeton Instruments) with 100 ms of exposure time.
Fourier transform-infrared (FT-IR) spectra of perovskite films on ITO/ZnO/ PEIE were conducted by a PIKE MIRacle Vertex 70 Spectrometer (Bruker) using a DLaTGS detector at room temperature. The spectrometer is equipped with a PIKE MIRacle TM attenuated total reflectance (ATR) accessory as the sampling accessory. The measuring system was continuously purged with N 2 before and during the measurements. The spectra were acquired at 2 cm −1 resolution with 30 scans over a wavenumber range between 4000 and 600 cm −1 .
Transmittance electron microscopy (TEM) of precursor films was performed by a TECNAI G2 F20 transmission electron microscope with an accelerating voltage of 200 kV and a Gatan SC200 CCD camera. Ultraviolet photoelectron spectroscopy (UPS) of perovskite films were recorded by a multifunctional photoelectron spectrometer (PHI 5000 Versaprobe II) with a He I ultraviolet radiation source (21.2 eV).
Nuclear magnetic resonance (NMR). High-resolution 1 H NMR data were recorded on a Bruker 600 MHz Avance III spectrometer at Larmor frequencies ν ( 1 H) = 600.274 MHz. 1 H NMR measurements single pulse experiments with 90°p ulse of 18 μs, relaxation delay 10 s, 64 scans were applied. Dried D 7 -DMF (99.5%) was used as the solvent and was also used to calibrate isotropic chemical shifts. The temperature was kept constant at 298 ± 0.2 K with a BVT 3000 temperature unit during measurements. All samples were dissolved in dried DMF at room temperature in an inert (Ar) atmosphere. The concentration for Pb 2+ and FAI is 0.12 and 0.24 M, respectively.
UV-vis absorbance spectra for solution and Benesi−Hildebrand (BH) analyses. UV-vis absorbance spectra for solutions were collected by a PerkinElmer Lambda 35 at room temperature. Anhydrous DMF was used as the solvent. The concentration of Pb 2+ for all the solutions is fixed at 0.25 mM. The mole ratio of PbI 2 :m-PEG 2 -NH 2 and PbI 2 :NH 2 -PEG 4 -NH 2 parent solution is 1:10 and 1:5, respectively. The formation equilibrium constant (K f ) of PbI 3 − was deduced by the Benesi−Hildebrand equation stated below: Here, A 0 is PbI 3 − absorbance of the parent solution without FAI addition; A is the collected absorbance with FAI; A max is the saturated absorbance with FAI addition; K f is deduced from the slope of Benesi−Hildebrand plot.
Density-functional theory (DFT) calculations. We performed the densityfunctional theory (DFT) calculations by using the projector-augmented wave (PAW) method as implemented in Vienna Ab initio Simulation Package (VASP) 41,42 . The generalized gradient approximation (GGA) with Perdew-Burke-Ernzerhof (PBE) method was used to optimize the crystal structures of cubic-phase FAPbI 3 bulk, slabs, and interfaces. The wave functions were expanded on a plane-wave basis set up to a kinetic energy cutoff of 450 eV. The van der Waals functional (vdw) was also included in the structural optimizations and electronic calculations using the zero-damping DFT-D3 method of Grimme. Monkhorst-Pack-type K-meshes of 6 × 6 × 6 for the bulk FAPbI 3 , 4 × 4 × 1 for the FAPbI 3 (110) slabs with different terminations, and 2 × 2 × 1 for the super-slabs and interfaces with m-PEG 2 -NH 2 or m-PEG 4 -NH 2 modification. All the structures were optimized until the forces on every single atom were less than 0.01 eV/Å. The binding energies were defined as E B = E interface -E FAPbI3 -E ligand , where E interface is the total energy of the interface after m-PEG 2 -NH 2 or m-PEG 4 -NH 2 passivation, and E FAPbI3 and E ligand represent the total energies of the individual parts. The molecular graphics viewer VESTA was used to plot the crystal structures. The geometrical and electronic properties of additives were calculated by B3LYP (Becke Three Parameters Hybrid Functional with Lee-Yang-Perdew correlation functions) functional and 6-31G* basis set with Gaussian 09 program package. The electrostatic potential surfaces were visualized and drawn using Multiwfn and VMD software.

Data availability
The data that support the plots within this paper and other findings of this study are available from the corresponding author upon reasonable request.