Unravelling the Effects of A-Site Cations on Nonradiative Electron

Aug 10, 2018 - The electron–hole recombination is slow, ranging from subnanosecond to nanoseconds, because the NA coupling is small, less than 3 meV...
0 downloads 0 Views 909KB Size
Subscriber access provided by Kaohsiung Medical University

Energy Conversion and Storage; Plasmonics and Optoelectronics

Unravelling the Effects of A Site Cations on Nonradiative Electron-Hole Recombination in Lead Bromide Perovskites: Time-Domain Ab Initio Analysis Jinlu He, Wei-Hai Fang, and Run Long J. Phys. Chem. Lett., Just Accepted Manuscript • DOI: 10.1021/acs.jpclett.8b02115 • Publication Date (Web): 10 Aug 2018 Downloaded from http://pubs.acs.org on August 13, 2018

Just Accepted “Just Accepted” manuscripts have been peer-reviewed and accepted for publication. They are posted online prior to technical editing, formatting for publication and author proofing. The American Chemical Society provides “Just Accepted” as a service to the research community to expedite the dissemination of scientific material as soon as possible after acceptance. “Just Accepted” manuscripts appear in full in PDF format accompanied by an HTML abstract. “Just Accepted” manuscripts have been fully peer reviewed, but should not be considered the official version of record. They are citable by the Digital Object Identifier (DOI®). “Just Accepted” is an optional service offered to authors. Therefore, the “Just Accepted” Web site may not include all articles that will be published in the journal. After a manuscript is technically edited and formatted, it will be removed from the “Just Accepted” Web site and published as an ASAP article. Note that technical editing may introduce minor changes to the manuscript text and/or graphics which could affect content, and all legal disclaimers and ethical guidelines that apply to the journal pertain. ACS cannot be held responsible for errors or consequences arising from the use of information contained in these “Just Accepted” manuscripts.

is published by the American Chemical Society. 1155 Sixteenth Street N.W., Washington, DC 20036 Published by American Chemical Society. Copyright © American Chemical Society. However, no copyright claim is made to original U.S. Government works, or works produced by employees of any Commonwealth realm Crown government in the course of their duties.

Page 1 of 30 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry Letters

Unravelling the Effects of A Site Cations on Nonradiative Electron-Hole Recombination in Lead Bromide Perovskites: Time-Domain Ab Initio Analysis Jinlu He,1 Wei-Hai Fang1, Run Long1∗ 1

College of Chemistry, Key Laboratory of Theoretical & Computational

Photochemistry of Ministry of Education, Beijing Normal University, Beijing, 100875, P. R. China Abstract: Lead bromide perovskites APbBr3 (A = Cs, MA, FA) hold great promise in optoelectronics and photovoltaics. Because the bandgap of the three materials are similar, and also because the A site cation does not contribute to band edges, leading one to expect minor influence of A site cation on the excited-state lifetime of the perovskites. Experiments defy the expectation. By performing ab initio nonadiabatic (NA) molecular dynamics combined with time-domain density functional simulations, we demonstrate that the nonradiative electron-hole recombination times are on the order FAPbBr3 > MAPbBr3 > CsPbBr3, which are determined by the NA electron-phonon coupling because decoherence times are similar. The simulations show that the larger A site cation and the smaller NA due to larger A site cation suppressing the Pb-Br cages motions. The electron-hole recombination is slow, ranging from sub-nanosecond to nanoseconds, because the NA coupling is small, less than 3 meV, and because decoherence time is slow, less than 7 fs. Both the trend of recombination and the time scales show excellent agreement with experiments. The ∗

Corresponding author, E-mail: [email protected]

1

ACS Paragon Plus Environment

The Journal of Physical Chemistry Letters 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

time-domain atomistic simulations rationalize the experimental observations and advance our understanding of the cations on influencing perovskite excited-state lifetimes.

TOC only

2

ACS Paragon Plus Environment

Page 2 of 30

Page 3 of 30 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry Letters

Hybrid inorganic-organic perovskites (HIOPs) solar cells have attracted intense attention due to rapidly growing photon-to-current conversion efficiencies increasing from 3.8%1 to 22.7%2 within eight years. The significant improvement of power conversion efficiencies originates from the excellent electronic and optical properties, including

low-trap

density,3

long

charge

carrier

lifetimes4

and

high

photoluminescence (PL) yields.5 However, HOIPs suffer from degradation under humidity and photo-irradiation and those issues weaken their photovoltaic properties.6 This drawback stimulates synthesizing all inorganic perovskites, such as cesium metal halides perovskites, which show better light and moisture stability.7 Similar to HOIPs, the all inorganic perovskites contain many excellent properties such as high fluorescence quantum yields,8 and long carrier diffusion lengths.9 These advanced properties make all inorganic perovskites are excellent candidates for a number of practical applications, including light-emitting diodes,10-13 optically pumped lasers,14-16 photocatalytic water-splitting assemblies,17-19 and solar cells.20-22 While the power conversion efficiencies of inorganic perovskites are lower than organic– inorganic hybrid halide perovskites7 stemming from the short excited state lifetimes.23 It is interesting because the bandgap of the APbBr3 (A =Cs, CH3NH3 (MA), HC(NH2)2 (FA)) is close to each other (CsPbBr3: 2.36 eV;24 MAPbBr3: 2.34 eV;25 FAPbBr3: 2.26 eV26), and also because the A site cation

does not directly contribute

to band edges. As a result, one expect that they should have an insignificant influence on the excited state lifetime of the APbBr3. 3

ACS Paragon Plus Environment

The Journal of Physical Chemistry Letters 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Recently, Jana et al. have demonstrated using time-resolved photoluminescence spectroscopy correlated with scanning electron microscopy that the average fluorescence lifetime of lead bromine perovskites APbBr3 are significantly different, in which the organic cations carry a fundamental advantage over the inorganic cation for achieving the long-lived excited state lifetime.23 One may wonder what the factors lead to the different excited state lifetimes in APbBr3 perovskites, because the band edges are entirely composed by lead and bromine orbitals.27 In order to explore the underlying mechanism responsible for the disparity in the excited state lifetime, first-principles atomistic time-domain simulations are needed to simulate the nonradiative electron-hole recombination in APbBr3, generating the mechanistic understanding on the excited state quantum dynamics, and providing valuable guidelines for materials design and further improvement. Previous first-principles ab initio quantum-classical dynamics simulations showed that the rearrangements of organic cations affect the electron-hole recombination in MAPbBr3 because they change the bandgap, modify the mixing between electron and hole wave function, and alter the phonon-induced quantum coherence time.28 Wang and coauthors demonstrated, using linear scaling ab initio methods, that random orientation of organic cations causes a long-range potential fluctuations in HOIPs, which spatially locates electron and hole in different positons, extending the charge-carrier lifetime.29 Similarly, the photoinduced deformation of inorganic Pb-Br sublattice also extends the charge-carrier lifetime by several orders of 4

ACS Paragon Plus Environment

Page 4 of 30

Page 5 of 30 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry Letters

magnitude compared with the perfect crystal in all inorganic perovskite.30-31 Motivated by the recent experimental and theoretical works,23, 28-29, 32-34 we have performed ab initio time domain35 nonadiabatic molecular dynamics (NAMD)36 simulations in order to explore the atomistic mechanisms serving to the observed disparity in the nonradiative electron-hole recombination time scales among the three lead bromide perovskites. The calculations provide a mechanistic understanding of the influence of A site cation (Cs, MA, FA) on the excited state lifetimes and the perovskite solar cell performance. The obtained nonradiative electron−hole recombination times are ranging from sub-nanosecond to several nanoseconds, and both the trend in the recombination and the time scales show excellent agreement with available experimental data.32-34 Our calculations show that electron-hole recombination proceeding rapidly in CsPbBr3, followed by MAPbBr3 and FAPbBr3 can be attributed to the small magnitudes of NA coupling beucase decoherence times are similar. Despite cations (Cs, MA, FA) contribute neither to the highest occupied (HOMO) nor to the lowest unoccupied molecular orbital (LUMO), they influence the atomic motions of lead and bromine, which primarily constitute the LUMO and HOMO in lead bromine perovskites,27 corresponding to the initial and final states for electron-hole recombination in the three systems. The size of cation is on the order FA > MA > Cs, leading to suppression of vibrations of inorganic Pb-I sublattice on the same trend. As a result, the NA electron-phonon coupling has the inverse order, achieving the 5

ACS Paragon Plus Environment

The Journal of Physical Chemistry Letters 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

short-lived charge-carrier lifetime in CsPbBr3, and extending significantly in MAPbBr3 and FAPbBr3. The charge-carrier lifetimes are long, because the NA coupling is small, below 3 meV, and the time of quantum coherence is short, below 7 fs. The simulations advance our understanding of the influence of A site cation on the excited dynamics in lead bromine perovskites, and provide valuable insights into enhancing the performance of lead bromide perovskites solar cells. The NAMD simulations are performed by the decoherence-induced surface hopping (DISH) approach36 implemented within the framework of time-dependent Kohn-Sham density functional framework.37-38 The faster and lighter electrons are described quantum mechanically, while the slower and heavier nuclei are treated classically.38-40 DISH incorporates decoherence correction into quantum dynamics and captures the nature of nuclear trajectory branching resulting in surface hops.36 Decoherence correction is needed here, because the decoherence (pure-dephasing) time is significantly shorter than the electron−hole recombination time, which occur on sub-nanosecond to several nanoseconds.32-34 The decoherence time is calculated as the pure-dephasing time in the optical response theory.41 DISH36 is implemented40 within the classical path approximation, which significantly reduces the computational cost, allowing one to use a predetermined ground state trajectory to propagate the nonequilibrium electron-vibrational dynamics. The approach has been applied to study photoexcitation dynamics in a variety of systems, 28, 31, 42-47 including quantum dots,42 two-dimensional (2D) ruddlesden-popper perovskites,43 perovskite 6

ACS Paragon Plus Environment

Page 6 of 30

Page 7 of 30 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry Letters

passivated by water,44 containing grain boundaries,45 forming polarons28 and localized hole,31 and interfaced with TiO2.46-47 A detailed theoretical method can be found elsewhere.39-40 Geometry optimization, electronic structure, adiabatic MD, and NA coupling calculations were obtained by the Vienna ab initio simulation package (VASP)48. The generalized gradient approximation of Perdew-Burke-Ernzerhof (PBE)49 functional and projector-augmented wave (PAW)50 pseudopotentials were used. A 400 eV plane-wave basis energy cutoff was employed for geometry optimization and electronic structure calculations. The structure optimization for both lattice constants and atomic coordinates with Γ-centered 4 × 4 × 4 Monkhorst−Pack mesh stopped until the residual forces on each atom were below 0.01 eV/Å.51 A much denser 8 × 8 × 8 mesh was used for density of states calculations. The van der Waals interactions were described by the Grimme DFT-D3 method.52-53 To model electron-hole recombination, a 2 × 2 × 2 supercell was constructed using the optimized cubic phase unit cells of CsPbBr3, MAPbBr3, and FAPbBr3 in the present study, containing 40 atoms in the CsPbBr3 and 96 atoms in the MAPbBr3 and FAPbBr3 respectively. After geometry optimization at 0 K, the three systems were heated to 300 K by repeated velocity rescaling. Then, 6 ps adiabatic MD trajectories were generated at Γ point in the microcanonical ensemble with a 1 fs time step. To simulate the nonradiative electron−hole recombination, the first 1000 initial conditions were sampled from the adiabatic MD trajectories to perform NAMD simulations with the PYXAID code.39-40 7

ACS Paragon Plus Environment

The Journal of Physical Chemistry Letters 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The bandgaps of all three systems under investigation were scaled to the experimental value for each material for study of electron-hole recombination.24-26

Figure 1. Simulation cell showing the optimized geometry (top panel) of (a) CsPbBr3, (b) MAPbBr3 and (c) FAPbBr3 and three corresponding representative geometries during the MD run at 300 K (bottom panel). At room temperature, the distortion of the inorganic Pb-I sublattice is on the order CsPbBr3 > MAPbBr3> FAPbBr3 because larger cation suppresses Pb and I motions more significantly. Here, the size of cation is on the order FA > MA > Cs, and the deformation of the Pb-I framework has the same order.

Figure 1 presents the optimized geometries of the CsPbBr3, MAPbBr3, and FAPbBr3 at 0 K (top panels), and snapshots of the corresponding structures taken from the MD trajectory at 300 K (bottom panels). The optimized lattice constants of the CsPbBr3, MAPbBr3, and FAPbBr3 systems are 5.900, 5.920, and 5.990 Å, agreeing well with experimental data of 5.870,54 5.919,55 5.99456 for each material 8

ACS Paragon Plus Environment

Page 8 of 30

Page 9 of 30 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry Letters

respectively. The calculated averaged bond length of the Pb-Br in the corresponding material is 2.950, 2.974, and 3.002 Å respectively, also showing an excellent agreement with experiments, in which they are 2.960 Å,57 2.970 Å57 and 3.000 Å58 for CsPbBr3, MAPbBr3, and FAPbBr3 systems. At room temperature, the bond length increases to 3.087, 3.088, and 3.098 Å respectively. The increase in the Pb-Br bond length of the MAPbBr3 and FAPbBr3 is smaller than that of in the CsPbBr3 after heating at 300 K, because the larger A site organic cations suppress the deformation of inorganic sublattice, in particular for the FAPbBr3. The change in the Pb-Br bond lengths reflects the strong electron-vibrational coupling existing in the CsPbBr3, followed by MAPbI3 and FAPbBr3. This point can be further demonstrated by the lattice free volume, calculated by the unit cell volume subtracting the constituent ions volumes.59 In general, a larger free volume provides more space for atomic motions, enhancing electron-phonon coupling. The volume of CsPbBr3, MAPbBr3, and FAPbBr3 of the optimized unit cell is calculated to be 205.38, 207.47, 214.92 Å3. The ionic radii of Cs+, MA+, FA+, Pb2+, and Br- are 1.88, 2.17, 2.53, 1.16, and 1.96 Å, respectively.60-61 Each APbBr3 unit cell contains one A (Cs+, MA+, FA+), one Pb2+, and three Br-, adding the spheres volumes in the unit cell gives the total ion volumes of 128.99 Å3 for CsPbBr3, 143.96 Å3 for MAPbBr3, and 168.99 Å3 for FAPbBr3. Therefore, the free volumes for the bromine and lead in the CsPbBr3, MAPbBr3, and FAPbBr3 are 76.39, 63.51, 45.93 Å3, respectively. The ratio between the free volumes and the total volumes denotes the specific free volumes,62 which are 37.19%, 30.61%, 9

ACS Paragon Plus Environment

The Journal of Physical Chemistry Letters 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 10 of 30

and 21.37% for the perovskites containing Cs, MA, and FA, respectively. The smaller specific free volume indicates that the amplitude of ion motions is smaller and the electron-phonon coupling is weaker. Hence, it is expected that the electron-phonon coupling decreases in the sequence CsPbBr3 > MAPbBr3 > FAPbBr3. In addition to the reported average bond lengths and specific free volumes, we calculate the canonically averaged standard deviation of the position of each atoms

 ,  =  − 〈 〉 , where represents the location of atom  at time t along the 6 ps MD trajectories. A smaller value of standard deviation reflects smaller atomic fluctuation and weaker electron-phonon coupling. The standard deviations in the positions of Pb and Br atoms are calculated for the three systems, because the charge densities of HOMO and LUMO are primarily consisted by Pb and Br orbitals. Shown in Table 1, the standard deviations in the positions of Br/Pb atoms decrease in the following order: CsPbBr3 > MAPbBr3 > FAPbBr3. The decrease in atomic fluctuations reflects electron-phonon couplings decreasing in the same order. Thus, the analysis of bond length, specific free volume, and atomic fluctuation agrees with each other and confirms that electron-phonon coupling is strongest in the CsPbBr3, with reduction in the MAPbBr3 and FAPbBr3. Table 1. Standard Deviations in the Positions of Pb and Br Atoms in CsPbBr3, MAPbBr3 and FAPbBr3.

Pb

CsPbBr3

MAPbBr3

FAPbBr3

0.343

0.332

0.258

10

ACS Paragon Plus Environment

Page 11 of 30 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry Letters

Br

0.603

0.440

0.412

Figure 2 shows the projected density of states (PDOS) of the CsPbBr3, MAPbBr3 and FAPbBr3 systems, calculated at 0 K using the optimized geometries. The PDOS illustrates that the HOMO is primarily composed by the Br atoms, with secondary contributions of Pb atoms, and the LUMO arises primarily from Pb atoms, with small contributions of Br atoms, agreeing with previous works.

27, 63

The two key orbitals

constitute the initial and final states for electron-hole recombination in the three systems under investigation. Although the cations Cs+, MA+, and FA+ do not contribute to the band edges, they do affect the local inorganic Pb-Br structure via electrostatic interactions. The calculated direct bandgaps of CsPbBr3, MAPbBr3 and FAPbBr3 are 1.60 eV, 1.94 eV and 1.91 eV at the Γ point respectively, showing a good agreement with the bandgap for each material obtained using the same PBE functional.64-65 Due to well-known DFT problem,66 the obtained values are smaller than the experimental bandgap for CsPbBr3 (2.36 eV),24 MAPbBr3 (2.34 eV),25 and FAPbBr3 (2.26 eV).26 To directly mimic the experiments,23 we scaled the calculated bandgaps to their experimental values for simulation of electron-hole recombination in these materials. The influence of bandgap on the recombination time is small, because the recombination rate depends linearly on the bandgap according to gap law,67 and also because the difference in bandgaps and decoherence times (Table 2) of the three materials is small. Consequently, the electron-hole recombination times are 11

ACS Paragon Plus Environment

The Journal of Physical Chemistry Letters 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

determined by the NA electron-phonon coupling.

Figure 2. Projected density of states (PDOS) of (a) CsPbBr3, (b) MAPbBr3, and (c) FAPbBr3 systems, split into contributions from the Pb, Br, and A cation. In all three systems, cation A contributes neither to the LUMO nor to the HOMO.

The NA matrix element relies significantly on both the HOMO/LUMO mixing −iħϕ ϕ  and nuclear velocity d/dt . A larger mixing of the corresponding charge densities and faster nuclear dynamics typically generate stronger NA electron-phonon coupling. Shown in Figure 3 a-c, the HOMO is localized primarily on Br atoms, while the LUMO is localized mainly on Pb atoms. The wave function mixing is similar in all of the three perovskites. The difference in NA electron-phonon coupling stems from the nuclear velocity, Table 1. Because FA is larger than MA and Cs, and hence, the strength of the NA electron-phonon coupling decreases from CsPbBr3, MAPbBr3, and to FAPbBr3.

12

ACS Paragon Plus Environment

Page 12 of 30

Page 13 of 30 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry Letters

Figure 3. Charge densities of the HOMO and LUMO. (a) CsPbBr3, (b) MAPbBr3 and (c) FAPbBr3 systems.

In order to characterize the phonons that couple to the electronic transition, we computed Fourier transforms (FT) of autocorrelation functions (ACF) of the HOMO−LUMO bandgap fluctuations in the three systems. Figure 4 shows that low-frequency vibrations in the 30−300 cm−1 region govern the nonradiative relaxation in all the three systems. Figure 4a shows in the CsPbBr3 system that the major peak at 67 cm-1 and the secondary peak at 133 cm−1 can be attributed to the vibrations of the inorganic Pb-Br octahedron and motion of Cs+ cations, respectively.68 Figure 4b demonstrates in the MAPbBr3 that the dominant peak at 33 cm-1 can be assigned to the Pb-Br inorganic framework distortion.69 The peaks at 133 cm-1 and 233 cm-1 can be assigned to the lurching of MA+ cations70 and torsional mode of the MA cations.71 For Figure 4c illustrates in the FAPbBr3 that the primary peak at 33 cm-1 can be also assigned to the Pb-Br sublattice distortion.69 The peak at 121 cm-1 can be ascribed to the vibration of 13

ACS Paragon Plus Environment

The Journal of Physical Chemistry Letters 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

FA+ cations.72 The dominant vibrational mode creates the largest NA electron-phonon coupling and induces rapid loss of quantum coherence. Compared the spectral density of CsPbBr3 with the MAPbBr3 and FAPbBr3, the main peak in the CsPbBr3 has higher frequency because the size of Cs+ cations are smaller the MA+ and FA+ cations, which allow inorganic Pb-Br framework move more in space in CsPbBr3 than that in MAPbBr3 and FAPbBr3 and generate the strongest NA electron-phonon coupling and shortest decoherence time. The range of active vibrations in MAPbBr3 is wider than those in FAPbBr3, leading to stronger NA electron-phonon coupling and shorter decoherence time. The phonon modes originating from A site cation affect the interaction between themselves and the inorganic lead and bromide atoms via electric filed they created, affecting the electron-phonon coupling and quantum decoherence time. Overall, the NA electron-coupling decreases in perovskite containing Cs, MA, and FA, and the docoherence time decreases in the reverse order, Table 2. One should be noted that decoherence time of the three systems are similar and have insignificant influence on the difference of the charge recombination times.

14

ACS Paragon Plus Environment

Page 14 of 30

Page 15 of 30 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry Letters

Figure 4. Spectral densities obtained by Fourier transforms of the autocorrelation functions, inset of Figure 5, for the fluctuations of the HOMO−LUMO energy gap in (a) CsPbBr3, (b) MAPbBr3, and (c) FAPbBr3.

The pure-dephasing functions computed using the second-order cumulant approximation in optical response theory41 are shown in Figure 5. The pure-dephasing times, τ, were obtained by fitting the pure-dephasing functions with a Gaussian, !" −0.5&/' . The pure-dephasing times decrease from FAPbBr3, to the MAPbBr3, to the CsPbBr3, Table 2, corresponding to 6.64, 6.19 and 5.73 fs, respectively. The moderate drop of dephasing time from FAPbBr3 to CsPbBr3 can be attributed to the different phonon modes discussed above. The un-ACF, calculated from the fluctuations of the HOMO-LUMO energy gap, inset of Figure 5 further gives the origin of pure-dephasing times differing each other. Generally, the larger initial value of the un-ACF and the shorter dephasing time, because under cumulative approximation, the initial value of the un-ACF equals to the variance of the energy gap, reflecting the strength of electron-phonon coupling. As a result, the pure 15

ACS Paragon Plus Environment

The Journal of Physical Chemistry Letters 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

dephasing times of three systems decrease as the sequence: FAPbBr3 > MAPbBr3 > CsPbBr3. It is necessary to incorporate decoherence corrections into NAMD simulations because the 5-7 fs pure-dephasing times are extremely shorter than the sub-nanosecond to several nanoseconds electron−hole recombination times reported experimentally.32-34

Figure 5. Pure-dephasing functions for the LUMO−HOMO transition in CsPbBr3, MAPbBr3, and FAPbBr3. The timescales obtained by Gaussian fits are summarized in Table 1. The inset shows the un-ACF, whose initial values are equal to the bandgap fluctuation squared. In general, The bigger initial value results in the faster dephasing.41

Figure 6 shows the evolution of the population of the first excited state in the three perovskite systems. The nonradiative electron-hole recombination corresponds to a single-molecule process and decay exponentially. Calculation of the whole charge recombination dynamics up to tens of nanoseconds demands heavily computational 16

ACS Paragon Plus Environment

Page 16 of 30

Page 17 of 30 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry Letters

cost, and thus we ran 20 ps NAMD simulations, in which the populations decay linearly due to the limited range, and extrapolated the recombination times based on the short dynamics. Shown in Table 2, the times, ', are obtained with the short-time linear approximation to the exponential decay, (& = !"−&/' = 1 − −&/', which are 0.88, 1.28, and 3.11 ns for CsPbBr3, MAPbBr3 and FAPbBr3 respectively, agreeing well with the data reported experimentally.32-34 Generally, pure DFT functional such as PBE used here, underestimates the band gaps of all the lead perovskites, and thus enhances the NA electron-phonon coupling because small band gap favors electron and hole wave functions mixing, leading to a faster electron-hole recombination. As discussed above, the difference in electron-hole recombination times for the three perovskites are primarily stemming from the NA electron-phonon coupling, because the decoherence times of three system are similar. The NA coupling is in the order: CsPbBr3 > MAPbBr3 > FAPbBr3, the recombination as result of same trend. The research demonstrates that A site cation plays a key role in determining the excited state lifetime of lead bromide perovskites through mainly influencing the NA electron-phonon coupling, suggesting that rational choice of A site cation can significantly reduce charge and energy losses and increase the conversion efficiency of perovskite solar cells.

17

ACS Paragon Plus Environment

The Journal of Physical Chemistry Letters 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 18 of 30

Figure 6. Electron-hole recombination dynamics in CsPbBr3, MAPbBr3 and FAPbBr3.

Table 2 Experimental Bandgap, Averaged Absolute NA Coupling, Pure-Dephasing Time, and Nonradiative Electron-Hole Recombination Time for CsPbBr3, MAPbBr3 and FAPbBr3. Bandgap

NA coupling

Dephasing

Recombination

(eV)

(meV)

(fs)

(ns)

CsPbBr3

2.36

2.28

5.73

0.88

MAPbBr3

2.34

1.62

6.19

1.27

FAPbBr3

2.26

1.01

6.64

3.10

In summary, we have investigated the nonradiative electron−hole recombination in CsPbBr3, MAPbBr3 and FAPbBr3 perovskites by NAMD combined with TDDFT, in order to resolve the microscopic mechanism serving to the difference in excited state lifetimes of the three materials. The simulations support the experimental 18

ACS Paragon Plus Environment

Page 19 of 30 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry Letters

observations that perovskite containing larger A site cation favors longer excited state lifetime. The ionic radii are in the sequence: FA > MA > Cs, and as a result of the excited state lifetimes follow the same trend. This is surprising because the bandgaps and dephasing times in the three materials are similar. The simulations demonstrate that A site cations steer the atomic fluctuations of inorganic Pb-Br sublattice and affect the NA coupling although A site cations do not directly contribute to the band edge states. The larger the A site cation and the smaller atomic motions in Pb and Br, which reduces NA electron-phonon coupling and leads the recombination to proceed slow in perovskite containing FA, moderate in MA, and fast in Cs. The reported simulations provide a detailed description of the complex quantum dynamics in the lead bromine perovskites, generating important insights and suggesting design principles for improve perovskite solar cells performance via a rational choice of cations.

Acknowledgements This work was supported by the National Science Foundation of China, Grant Nos. 21573022, 51861135101, 21520102005, and 21421003. R. L. is grateful to the Fundamental Research Funds for the Central Universities, the Recruitment Program of Global Youth Experts of China, and the Beijing Normal University Startup.

19

ACS Paragon Plus Environment

The Journal of Physical Chemistry Letters 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

References (1) Kojima, A.; Teshima, K.; Shirai, Y.; Miyasaka, T. Organometal Halide Perovskites as Visible-Light Sensitizers for Photovoltaic Cells. J. Am. Chem. Soc. 2009, 131, 6050−6051. (2) Fu, F.; Pisoni, S.; Weiss, T. P.; Feurer, T.; Wäckerlin, A.; Fuchs, P.; Nishiwaki, S.; Zortea, L.; Tiwari, A. N.; Buecheler, S. Compositionally Graded Absorber for Efficient and Stable Near‐Infrared‐Transparent Perovskite Solar Cells. Adv. Sci. 2018, 5, 1700675. (3) Shi, D.; Adinolfi, V.; Comin, R.; Yuan, M. J.; Alarousu, E.; Buin, A.; Chen, Y.; Hoogland, S.; Rothenberger, A.; Katsiev, K.; et al. Low Trap-State Density and Long Carrier Diffusion in Organolead Trihalide Perovskite Single Crystals. Science 2015,

347, 519−522. (4) Bi, Y.; Hutter, E. M.; Fang, Y.; Dong, Q.; Huang, J.; Savenije, T. J. Charge Carrier Lifetimes Exceeding 15 Ms in Methylammonium Lead Iodide Single Crystals. J. Phys.

Chem. Lett. 2016, 7, 923−928. (5) Deschler, F.; Price, M.; Pathak, S.; Klintberg, L. E.; Jarausch, D.-D.; Higler, R.; Hüttner, S.; Leijtens, T.; Stranks, S. D.; Snaith, H. J.; et al. High Photoluminescence Efficiency and Optically Pumped Lasing in Solution-Processed Mixed Halide Perovskite Semiconductors. J. Phys. Chem. Lett. 2014, 5, 1421−1426. (6) Niu, G. D.; Guo, X. D.; Wang, L. D. Review of Recent Progress in Chemical Stability of Perovskite Solar Cells. J. Mater. Chem. A, 2014, 3, 8970-8980. 20

ACS Paragon Plus Environment

Page 20 of 30

Page 21 of 30 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry Letters

(7) Liang, J.; Zhao, P. Y.; Wang, C. X.; Wang, Y. R.; Hu, Y.; Zhu, G. Y.; Ma, L. B.; Liu, J.; Jin, Z. CsPb0.9Sn0.1IBr2 Based All-Inorganic Perovskite Solar Cells with Exceptional Efficiency and Stability. J. Am. Chem. Soc. 2017, 139, 14009-14012. (8) Swarnkar, A.; Chulliyil, R.; Ravi, V. K.; Irfanullah, M.; Chowdhury, A.; Nag, A. Colloidal CsPbBr3 Perovskite Nanocrystals: Luminescence beyond Traditional Quantum Dots. Angew. Chem. Int. Ed. 2016, 54, 15424-15428. (9) Yettapu, G. R.; Talukdar, D.; Sarkar, S.; Swarnkar, A.; Nag, A.; Ghosh, P.; Mandal, P. THz Conductivity within Colloidal CsPbBr3 Perovskite Nanocrystals: Remarkably High Carrier Mobilities and Large Diffusion Lengths. Nano Lett. 2016, 16, 4838-4848. (10) Wong, A. B.; Lai, M. L.; Eaton, S. W.; Yu, Y.; Lin, E.; Dou, L. T.; Fu, A.; Yang, P. D. Growth and Anion Exchange Conversion of CH3NH3PbX3 Nanorod Arrays for Light-Emitting Diodes. Nano Lett. 2015, 15, 5519-5524. (11) Leyden, M. R.; Meng, L. Q.; Jiang, Y.; Ono, L. K.; Qiu, L. B.; Juarez-Perez, E. J.; Qin, C. J.; Adachi, C.; Qi, Y. B. Methylammonium Lead Bromide Perovskite Light Emitting Diodes by Chemical Vapor Deposition. J. Phys. Chem. Lett. 2017, 8, 3193-3198. (12) Song, L.; Guo, X. Y.; Hu, Y. S.; Lv, Y.; Lin, J.; Liu, Z. Q.; Fan, Y.; Liu, X. Y. Efficient Inorganic Perovskite Light-Emitting Diodes with Polyethylene Glycol Passivated Ultrathin CsPbBr3 Films. J. Phys. Chem. Lett. 2017, 8, 4148-4154. (13) Yao, J. S.; Ge, J.; Han, B. N.; Wang, K. H.; Yao, H. B.; Yu, H. L.; Li, J. H.; Zhu, 21

ACS Paragon Plus Environment

The Journal of Physical Chemistry Letters 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

B. S.; Song, J. Z.; Chen, C.; et al. Ce3+-Doping to Modulate Photoluminescence Kinetics for Efficient CsPbBr3 Nanocrystals Based Light-Emitting Diodes. J. Am.

Chem. Soc. 2018, 140, 3626-3634. (14) Zhu, H. M.; Fu, Y. P.; Meng, F.; Wu, X. X.; Gong, Z. Z.; Ding, Q.; Gustafsson, M. V.; Trinh, M. T.; Jin, S.; Zhu, X. Y. Lead Halide Perovskite Nanowire Lasers with Low Lasing Thresholds and High Quality Factors. Nat. Mater. 2015, 14, 636-642. (15) Fu, Y. P.; Zhu, H. M.; Schrader, A. W.; Liang, D.; Ding, Q.; Joshi, P.; Hwang, L.; Zhu, X. Y.; Jin, S. Nanowire Lasers of Formamidinium Lead Halide Perovskites and Their Stabilized Alloys with Improved Stability. Nano Lett. 2016, 16, 1000-1008. (16) Xu, Y. P.; Chen, Q.; Zhang, C. F.; Wang, R.; Wu, H.; Zhang, X. Y.; Xing, G. H.; Yu, W. W.; Wang, X. Y.; Zhang, Y.; et al. Two-PhotonPumped Perovskite Semiconductor Nanocrystal Lasers. J. Am. Chem. Soc. 2016, 138, 3761-3768. (17) Gurudayal; Sabba, D.; Kumar, M. H.; Wong, L. H.; Barber, J.; Grätzel, M.; Mathews, N. Perovskite–Hematite Tandem Cells for Efficient Overall Solar Driven Water Splitting. Nano Lett. 2015, 15, 3833-3839. (18) Chen, Y.-S.; Manser, J. S.; Kamat, P. V. All Solution-Processed Lead Halide Perovskite-BiVO4 Tandem Assembly for Photolytic Solar Fuels Production. J. Am.

Chem. Soc. 2015, 137, 974-981. (19) Kim, J. H.; Jo, Y.; Kim, J. H.; Jang, J. W.; Kang, H. J.; Lee, Y. H.; Kim, D. S.; Jun, Y.; Lee, J. S. Wireless Solar Water Splitting Device with Robust Cobalt-Catalyzed, Dual-Doped BiVO4 Photoanode and Perovskite Solar Cell in 22

ACS Paragon Plus Environment

Page 22 of 30

Page 23 of 30 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry Letters

Tandem: A Dual Absorber Artificial Leaf. ACS Nano 2015, 9, 11820-11829. (20) Hoffman, J. B.; Zaiats, G.; Wappes, I.; Kamat, P. V. CsPbBr3 Solar Cells: Controlled Film Growth through Layer-by-Layer Quantum Dot Deposition. Chem.

Mater. 2017, 29, 9767–9774. (21) Panigrahi, S.; Jana, S.; Calmeiro, T.; Nunes, D.; Martins, R.; Fortunato, E. Imaging the Anomalous Charge Distribution inside CsPbBr3 Perovskite Quantum Dots Sensitized Solar Cells. ACS Nano 2017, 11, 10214−10221. (22) Kulbak, M.; Gupta, S.; Kedem, N.; Levine, I.; Bendikov, T.; Hodes, G.; Cahen, D. Cesium Enhances Long-Term Stability of Lead Bromide Perovskite-Based Solar Cells.

J. Phys. Chem. Lett. 2016, 7, 167-172. (23) Jana, A.; Mittal, M.; Singla, A.; Sapra, S. Solvent-Free, Mechanochemical Syntheses of Bulk Trihalide Perovskites and Their Nanoparticles. Chem. Commun. 2017, 53, 3046-3049. (24) Hoffman, J. B.; Schleper, A. L.; Kamat, P. V. Transformation of Sintered CsPbBr3 Nanocrystals to Cubic CsPbI3 and Gradient CsPbBrxI3−x through Halide Exchange. J. Am. Chem. Soc. 2016, 138, 8603−8611. (25) Hills-Kimball, K.; Nagaoka, Y.; Cao, C.; Chaykovsky, E.; Chen, O. Synthesis of Formamidinium Lead Halide Perovskite Nanocrystals through Solid–Liquid–Solid Cation Exchange. J. Mater. Chem. C, 2017, 5, 5680--5684. (26) Arora, N.; Orlandi, S.; Dar, M. I.; Aghazada, S.; Jacopin, G. n.; Cavazzini, M.; Mosconi, E.; Gratia, P.; De Angelis, F.; Pozzi, G.; et al. High Open-Circuit Voltage: 23

ACS Paragon Plus Environment

The Journal of Physical Chemistry Letters 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Fabrication of Formamidinium Lead Bromide Perovskite Solar Cells Using Fluorene– Dithiophene Derivatives as Hole-Transporting Materials. ACS Energy Lett. 2016, 1, 107-112. (27) Protesescu, L.; Yakunin, S.; Bodnarchuk, M. I.; Krieg, F.; Caputo, R.; Hendon, C. H.; Yang, R. X.; Walsh, A.; Kovalenko, M. V. Nanocrystals of Cesium Lead Halide Perovskites (CsPbX3, X = Cl, Br, and I): Novel Optoelectronic Materials Showing Bright Emission with Wide Color Gamut. Nano Lett. 2015, 15, 3692−3696. (28) Zhang, Z. S.; Long, R.; Tokina, M. V.; Prezhdo, O. V. Interplay between Localized and Free Charge Carriers Can Explain Hot Fluorescence in the CH3NH3PbBr3 Perovskite: Time-Domain Ab Initio Analysis. J. Am. Chem. Soc. 2017,

139, 17327−17333. (29) Ma, J.; Wang, L.-W. Nanoscale Charge Localization Induced by Random Orientations of Organic Molecules in Hybrid Perovskite CH3NH3PbI3. Nano Lett. 2015, 15, 248−253. (30) Santomauro, F. G.; Grilj, J.; Mewes, L.; Nedelcu, G.; Yakunin, S.; Rossi, T.; Capano, G.; Al Haddad, A.; Budarz, J.; Kinschel, D.; et al. Localized Holes and Delocalized Electrons in Photoexcited Inorganic Perovskites: Watching Each Atomic Actor by Picosecond X-Ray Absorption Spectroscopy. Struct. Dyn. 2017, 4, 044002. (31) He, J. L.; Guo, M.; Long, R. Photoinduced Localized Hole Delays Nonradiative Electron-Hole Recombination in Cesium-Lead Halide Perovskites: A Time-Domain Ab Initio Analysis. J. Phys. Chem. Lett. 2018, 9, 3021-3028. 24

ACS Paragon Plus Environment

Page 24 of 30

Page 25 of 30 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry Letters

(32) Huang, H.; Li, J.; Yi, Y.; Wang, J. H.; Kang, Y. H.; Chu, P. K.; Ong, H. C.; Yu, X.-F. In Situ Growth of All-Inorganic Perovskite Nanocrystals on Black Phosphorus Nanosheets. Chem. Commun. 2018, 54, 2365-2368. (33) Zhang, Z.-Y.; Wang, H.-Y.; Zhang, Y.-X.; Hao, Y.-W.; Sun, C.; Zhang, Y.; Gao, B.-R.; Chen, Q.-D.; Sun, H.-B. The Role of Trap-assisted Recombination in Luminescent Properties of Organometal Halide CH3NH3PbBr3 Perovskite Films and Quantum Dots.

Sci. Rep. 2016, 6, 27286. (34) Das, J.; Subbiah, A. S.; Mahuli, N.; Singh, R.; Sarkar, S. K. One-Step Solution-Processed

Formamidinium

Lead

Tribromide

Formation

for

Better

Reproducible Planar Perovskite Solar Cells. Energy Technol. 2017, 5, 1807–1813. (35) Runge, E.; Gross, E. K. U. Density-Functional Theory for Time-Dependent Systems. Phys. Rev. Lett. 1984, 52, 997-1000. (36) Jaeger, H. M.; Fischer, S.; Prezhdo, O. V. Decoherence-Induced Surface Hopping.

J. Chem. Phys. 2012, 137, 22A545. (37) Kohn, W.; Sham, L. J. Self-Consistent Equations Including Exchange and Correlation Effects. Phys. Rev. 1965, 140, A1133−A1138. (38) Craig, C. F.; Duncan, W. R.; Prezhdo, O. V. Trajectory Surface Hopping in the Time-Dependent Kohn-Sham Approach for ElectronNuclear Dynamics. Phys. Rev.

Lett. 2005, 95, 163001. (39) Akimov, A. V.; Prezhdo, O. V. The Pyxaid Program for NonAdiabatic Molecular 25

ACS Paragon Plus Environment

The Journal of Physical Chemistry Letters 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Dynamics in Condensed Matter Systems. J. Chem. Theory Comput. 2013, 9, 4959−4972. (40) Akimov, A. V.; Prezhdo, O. V. Advanced Capabilities of the Pyxaid Program: Integration Schemes, Decoherence Effects, Multiexcitonic States, and Field-Matter Interaction. J. Chem. Theory Comput. 2014, 10, 789−804. (41) Mukamel, S.: Principles of Nonlinear Optical Spectroscopy; Oxford University Press: New York, 1995. (42) He, J. L.; Vasenko, A. S.; Long, R.; Prezhdo, O. V. Halide Composition Controls Electron−Hole Recombination in Cesium−Lead Halide Perovskite Quantum Dots: A Time Domain Ab Initio Study. J. Phys. Chem. Lett. 2018, 9, 1872−1879. (43) Zhang, Z. S.; Fang, W. H.; Tokina, M. V.; Long, R.; Prezhdo, O. V. Rapid Decoherence Suppresses Charge Recombination in Multi-Layer 2D Halide Perovskites: Time-Domain Ab Initio Analysis. Nano Lett. 2018, 18, 2459-2466. (44) Long, R.; Fang, W. H.; Prezhdo, O. V. Moderate Humidity Delays Electron-Hole Recombination in Hybrid Organic-Inorganic Perovskites: Time-Domain Ab Initio Simulations Rationalize Experiments. J. Phys. Chem. Lett. 2016, 7, 3215−3222. (45) Long, R.; Liu, J.; Prezhdo, O. V. Unravelling the Effects of Grain Boundary and Chemical Doping on Electron-Hole Recombination in CH3NH3PbI3 Perovskite by Time-Domain Atomistic Simulation. J. Am. Chem. Soc. 2016, 138, 3884−3890. (46) Long, R.; Prezhdo, O. V. Dopants Control Electron-Hole Recombination at Perovskite-TiO2 Interfaces: Ab Initio Time-Domain Study. ACS Nano 2015, 9, 26

ACS Paragon Plus Environment

Page 26 of 30

Page 27 of 30 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry Letters

11143-11155. (47) Long, R.; Fang, W.-H.; Prezhdo, O. V. Strong Interaction at the Perovskite/TiO2 Interface Facilitates Ultrafast Photoinduced Charge Separation: A Nonadiabatic Molecular Dynamics Study. J. Phys. Chem. C 2017, 121, 3797−3806. (48) Kresse, G.; Furthmuller, J. Efficient Iterative Schemes for Ab Initio Total-Energy Calculations using a Plane-Wave Basis Set. Phys. Rev. B: Condens. Matter Mater.

Phys. 1996, 54, 11169−1186. (49) Perdew, J. P.; Burke, K.; Ernzerhof, M. Generalized Gradient Approximation Made Simple. Phys. Rev. Lett. 1996, 77, 3865−3868. (50) Blöchl, P. E. Projector Augmented-Wave Method. Phys. Rev. B: Condens. Matter

Mater. Phys. 1994, 50, 17953−17979. (51) Monkhorst, H. J.; Pack, J. D. Special Points for Brillouin-Zone Integrations. Phys.

Rev. B 1976, 13, 5188−5192. (52) Grimme, S.; Antony, J.; Ehrlich, S.; Krieg, H. A Consistent and Accurate Ab Initio Parametrization of Density Functional Dispersion Correction (Dft-D) for the 94 Elements H-Pu. J. Chem. Phys. 2010, 132, 154104. (53) Stefan Grimme, S. E.; Goerigk, L. Effect of the Damping Function in Dispersion Corrected Density Functional Theory. J. Comput. Chem. 2011, 32, 1456−1465. (54) Rodova, M.; Broz ́ ek, J.; Kníz ek, K.; Nitsch, K. Phase Transitions in Ternary Caesium Lead Bromide. J. Therm. Anal. Cal.Orim. 2003, 71, 667−673. (55) Elbaz, G. A.; Straus, D. B.; Semonin, O. E.; Hull, T. D.; Paley, D. W.; Kim, P.; 27

ACS Paragon Plus Environment

The Journal of Physical Chemistry Letters 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Owen, J. S.; Kagan, C. R.; Roy, X. Unbalanced Hole and Electron Diffusion in Lead Bromide Perovskites. Nano Lett. 2017, 17, 1727−1732. (56) Zhumekenov, A. A.; Saidaminov, M. I.; Haque, M. A.; Alarousu, E.; Sarmah, S. P.; Murali, B.; Dursun, I.; Miao, X. H.; Abdelhady, A. L.; Wu, T. Formamidinium Lead Halide Perovskite Crystals with Unprecedented Long Carrier Dynamics and Diffusion Length. ACS Energy Lett. 2016, 1, 32-37. (57) Knop, O.; Wasylishen, R. E.; White, M. A.; Cameron, S.; Michiel, J. M.; Oort, V. Alkylammonium Lead Halides. Part 2. CH3NH3PbX3(X = C1, Br, I) Perovskites: Cuboctahedral Halide Cages with Isotropic Cation Reorientation. Can. J. Chem. 1990,

68, 412-422. (58) Sun, S. J.; Isikgor, F. H.; Deng, Z. Y.; Wei, F. X.; Kieslich, G.; Bristowe, P. D.; Ouyang, J.; Cheetham, A. K. Factors Influencing the Mechanical Properties of Formamidinium Lead Halides and Related Hybrid Perovskites. ChemSusChem 2017,

10, 3740 – 3745. (59) Sammells, A. F.; Cook, R. L.; White, J. H.; Osborne, J. J.; MacDuff, R. C. Rational Selection of Advanced Solid Electrolytes for Intermediate Temperature Fuel Cells. Solid State Ionics 1992, 52, 111−123. (60) Nedelcu, G.; Protesescu, L.; Yakunin, S.; Bodnarchuk, M. I.; Grotevent, M. J.; Kovalenko, M. V. Fast Anion-Exchange in Highly Luminescent Nanocrystals of Cesium Lead Halide Perovskites (CsPbX3, X= Cl, Br, I). Nano Lett. 2015, 15, 5635−5640. 28

ACS Paragon Plus Environment

Page 28 of 30

Page 29 of 30 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry Letters

(61) Kieslich, G.; Sun, S.; Cheetham, A. K. An Extended Tolerance Factor Approach for Organic-Inorganic Perovskites. Chem. Sci. 2015, 6, 3430-3433. (62) Hayashi, H.; Inaba, H.; Matsuyama, M.; Lan, N. G.; Dokiya, M.; Tagawa, H. Structural Consideration on the Ionic Conductivity of Perovskite-Type Oxides. Solid

State Ionics 1999, 122, 1−15. (63) Jacobsson, T. J.; Correabaena, J. P.; Pazoki, M.; Saliba, M.; Schenk, K.; Grätzel, M.; Hagfeldt, A. Exploration of the Compositional Space for Mixed Lead Halogen Perovskites for High Efficiency Solar Cells. Energy Environ. Sci. 2016, 9, 1706-1724. (64) Motta, C.; El-Mellouhi, F.; Sanvito, S. Exploring the Cation Dynamics in Lead-Bromide Hybrid Perovskites. Phys. Rev. B 2016, 93,235412. (65) Yuan, Y.; Xu, R.; Xu, H. T.; Hong, F.; Xu, F.; Wang, L. J. Nature of the Band Gap of Halide Perovskites ABX3 (A = CH3NH3, Cs; B = Sn, Pb; X = Cl, Br, I): First-Principles Calculations. Chin. Phys. B 2015, 24, 116302−116303. (66) Yakovkin, I. N.; Dowben, P. A. The Problem of the Band Gap in Lda Calculations. Surf. Rev. Lett. 2007, 14, 481-487. (67) Englman, R.; Jortner, J. The Energy Gap Law for Non-Radiative Decay in Large Molecules. J. Lumin. 1970, 1-2, 134-142. (68) Cha, J. H.; Han, J. H.; Yin, W. P.; Park, C.; Park, Y.; Ahn, T. K.; Cho, J. H.; Jung, D. Y. Photoresponse of CsPbBr3 and Cs4PbBr6 Perovskite Single Crystals. J. Phys.

Chem. Lett. 2017, 8, 565−570. (69) Yaffe, O.; Guo, Y.; Tan, L. Z.; Egger, D. A.; Hull, T.; Stoumpos, C. C.; Zheng, F.; 29

ACS Paragon Plus Environment

The Journal of Physical Chemistry Letters 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Heinz, T. F.; Kronik, L.; Kanatzidis, M. G.; et al. Local Polar Fluctuations in Lead Halide Perovskite Crystals. Phys. Rev. Lett. 2017, 118, 13600. (70) Leguy, A. M.; Goñi, A. R.; Frost, J. M.; Skelton, J.; Brivio, F.; Rodríguez-Martínez, X.; Weber, O. J.; Pallipurath, A.; Alonso, M. I.; Campoy-Quiles, M.; et al. Dynamic Disorder, Phonon Lifetimes, and the Assignment of Modes to the Vibrational Spectra of Methylammonium Lead Halide Perovskites. Phys. Chem.

Chem. Phys. 2016, 18, 27051-27066. (71) Quarti, C.; Grancini, G.; Mosconi, E.; Bruno, P.; Ball, J. M.; Lee, M. M.; Snaith, H. J.; Petrozza, A.; Angelis, F. D. The Raman Spectrum of the CH3NH3PbI3 Hybrid Perovskite: Interplay of Theory and Experiment. J. Phys. Chem. Lett. 2014, 5, 279−284. (72) Han, Q. F.; Bae, S. H.; Sun, P. Y.; Hsieh, Y. T.; Yang, Y.; Rim, Y. S.; Zhao, H. X.; Chen, Q.; Shi, W. Z.; Li, G.; et al. Single Crystal Formamidinium Lead Iodide (FAPbI3): Insight into the Structural, Optical, and Electrical Properties. Adv. Mater. 2016, 28, 2253–2258.

30

ACS Paragon Plus Environment

Page 30 of 30