Lewis Base Passivation of Hybrid Halide Perovskites Slows Electron

Feb 20, 2018 - Lihong Liu†, Wei-Hai Fang† , Run Long† , and Oleg V. Prezhdo‡. † Key Laboratory of Theoretical and Computational Photochemist...
0 downloads 0 Views 749KB Size
Subscriber access provided by UNIV OF DURHAM

Letter

Lewis Base Passivation of Hybrid Halide Perovskites Slows Electron-Hole Recombination: Time-Domain Ab Initio Analysis Lihong Liu, Weihai Fang, Run Long, and Oleg V. Prezhdo J. Phys. Chem. Lett., Just Accepted Manuscript • DOI: 10.1021/acs.jpclett.8b00177 • Publication Date (Web): 20 Feb 2018 Downloaded from http://pubs.acs.org on February 20, 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.

The Journal of Physical Chemistry Letters 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 26 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

Lewis Base Passivation of Hybrid Halide Perovskites Slows Electron-Hole Recombination: Time-Domain Ab Initio Analysis

Lihong Liu1, Wei-Hai Fang1, Run Long1*, and Oleg V. Prezhdo2 1

Key Laboratory of Theoretical and Computational Photochemistry, Ministry of

Education, College of Chemistry, Beijing Normal University, Beijing 100875, China 2

Department of Chemistry, University of Southern California, Los Angeles, California 90089, United States

Abstract: Non-radiative electron-hole recombination plays a key role in determining photon conversion efficiencies in solar cells. Experiments demonstrate significant reduction in the recombination rate upon passivation of methylammonium lead iodide perovskite with Lewis base molecules. Using non-adiabatic molecular dynamics combined with time-domain density functional theory, we find that the nonradiative charge recombination is decelerated by an order of magnitude upon adsorption of the molecules. Thiophene acts by the traditional passivation mechanism, forcing electron density away from the surface. In contrast, pyridine localizes the electron at the surface, while leaving it energetically near the conduction band edge. This is because pyridine creates a stronger coordinative bond with a lead atom of the perovskite, and has a lower energy unoccupied orbital, compared to thiophene, due to the more electronegative nitrogen atom, relative to thiophene’s sulfur. Both molecules reduce twofold the non-adiabatic coupling and electronic coherence time. A broad range of vibrational modes couple to the electronic subsystem, arising from inorganic and

*

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

organic components. The simulations reveal the atomistic mechanisms underlying the enhancement of the excited state lifetime achieved by the perovskite passivation, rationalize the experimental results, and advance our understanding of charge-phonon dynamics in perovskite solar cells.

TOC graphic

2 ACS Paragon Plus Environment

Page 2 of 26

Page 3 of 26 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 organic-inorganic halide perovskites, such as CH3NH3PbI3 (MAPbI3), have attracted great attention as a new class of optoelectronic semiconductors.1-8 The power conversion efficiency of the MAPbI3 perovskite solar cells has increased from 3.8%1 to 22.1%4, 6 within only six years of active research, due to excellent electronic and optical properties, including long carrier diffusion lengths9 and high light absorption coefficients.1 Material performance is limited by defects that provide charge scattering and recombination sites. Most often, defects are located at surfaces and crystal grain boundaries that contain under-coordinated atoms. Significant theoretical and experimental efforts are dedicated to studies on defect chemistry in hybrid perovskites.10-14 It has been discovered that the energy of many intrinsic defects, such as halide vacancies, is close to band edges.15 Charges trapped in such shallow states can be easily de-trapped and take a long time to recombine,16 providing a plausible rationalization for the remarkable performance of perovskite solar cells. Non-radiative recombination rates are still relatively high, as inferred from low photoluminescence efficiencies under small excitation densities.17 Perovskites suffer from low stability against humidity,18-19 which increases defect concentration and constitutes a major obstacle to large-scale applications. Replacing some iodine atoms in MAPbI3 with chlorines extends excited state lifetimes, in part due to surface passivation,20-21 triggering experimental22-23 and theoretical24-25 studies of chlorine influence on geometric, electronic and dynamical properties of MAPbI3, and allowing improvements in solar cell performance. Focusing on mixed iodine/chlorine perovskite MAPbI3-xClx, Noel et al. undertook 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

Page 4 of 26

further steps for improving material quality by passivating perovskite crystal surfaces with organic Lewis base molecules.26 By treating perovskite films with pyridine and thiophene, they observed a significantly inhibited nonradiative charge recombination, especially at low photoexcitation levels, leading to power conversion efficiencies of 15.3 and 16.5% for the thiophene and pyridine-treated devices, respectively.26 The authors postulated that the observed decrease in non-radiative decay occurred due to passivation of under-coordinated lead atoms on perovskite crystal surfaces via coordinate bonding to the nitrogen atom in pyridine or the sulfur atom in thiophene. Besides solar cells, the passivation technique can be used with other optoelectronic applications of perovskites, such as lasers and light-emitting diodes.26-27 Both solar cells

and

optoelectronic

devices

require

slow

non-radiative

electron-hole

recombination, because this process constitutes a major pathway for charge and energy losses. Further progress requires a thorough understanding of the mechanism by which Lewis base passivation reduces charge losses. Such understanding can be provided by atomistic ab initio studies, particularly in the time-domain, directly mimicking the time-resolved experiments. 26, 28 In this letter, we employ real-time time-dependent density functional theory (TDDFT) and non-adiabatic molecular dynamics (NAMD) in order to investigate the atomistic mechanisms underlying the observed reduction of non-radiative electron-hole recombination in MAPbI3 achieved by surface passivation with the Lewis base molecules.26 The simulations indicate that pyridine and thiophene achieve similar effects but through different mechanisms. Thiophene acts by removing 4 ACS Paragon Plus Environment

Page 5 of 26 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

electrons density from the surface. In contrast, pyridine attracts electron into an unoccupied orbital that appears near the bottom of the perovskite conduction band (CB) and that is localized around the strong coordinative bond formed by the pyridine nitrogen and surface lead atoms. It is essential that no trap states are created, and the electron can easily escape into the CB. Both Lewis bases reduce the NA electron-vibrational coupling and electronic coherence time by a factor of two, enhancing the excited state lifetime by an order of magnitude. A broad spectrum of vibrational modes, arising from the inorganic and organic parts, couple to the electronic subsystem. The state-of-the-art simulations highlight various factors affecting the non-radiative electron-hole recombination, and provide valuable chemical and physical insights that can be used to improve further the performance of perovskite solar cells. The simulations employ the decoherence induced surface hopping (DISH) approach29 implemented30 within real-time time-dependent density functional theory (TDDFT) in the Kohn-Sham framework.31-32 The DISH approach for NAMD incorporates loss of quantum coherence into the quantum-classical approximation by introducing transitions of classical trajectories between quantum states at decoherence events. The lighter and faster electrons are treated quantum mechanically, whereas the heavier and slower nuclei are described classically. Directly related to pure-dephasing of the optical response theory,20,

33-35

decoherence should be included in the

calculation, because it is significantly faster than the electron-hole recombination. We applied this approach to study excitation dynamics in diverse systems,33, 5 ACS Paragon Plus Environment

36-48

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

including perovskites containing dopants,49 grain boundaries35 and defects,50 in contact with TiO220, 51 and water,34 and exhibiting polarons52 and ordered phases.53 A detailed description of the theoretical approach can be found elsewhere.30, 54 Considering the previous experimental and theoretical studies of the MAPbI3 systems,55-56 we chose the Cl-doped PbI-terminated tetragonal MAPbI3 (001) surface to model electron-hole recombination with and without the adsorbed Lewis base molecules. The periodically repeated 72-atom (1×1) surface is composed of three MAPbI3 layers, with the bottom layer frozen in the bulk configuration. The geometry optimization, adiabatic MD, and NA coupling calculations are carried out using the Vienna ab initio simulation package (VASP). 57 The exchange-correlation interactions are treated using the Perdew-Burke-Ernzerhof (PBE) functional.58 The interaction between the ionic cores and the valence electrons are described by the projector-augmented wave approach.59 The van der Waals interactions are treated by the Grimme DFT-D2 method.60 Spin-orbit interaction is important in MAPbI3 due to in the presence of the heavy elements, Pb and I. However, it has been shown that semilocal PBE without spin-orbit interactions is able to provide good agreement with the experimental bandgaps61 and generate important insights into perovskite electronic structure,62-63 justifying the use of the less expensive methods for exploring novel effects. The errors due to lack of spin-orbit interaction and presence of electron self-interaction cancel in the simple PBE calculation. In order to obtain an accurate bandgap with spin-orbit interaction turned on, one needs to either use a hybrid functional or perform a GW calculation.64 6 ACS Paragon Plus Environment

Page 6 of 26

Page 7 of 26 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

It is impossible to perform quantum dynamics calculations on the current system with either hybrid functional or GW, because of large computational efforts. Therefore, we use simple PBE, which gives good agreement with experiment.61 This approach provided good results in our previous studies of hybrid organic-inorganic perovskites.20, 34-35, 51-52 The simulated systems are shown in Figure 1. The periodic images are separated along the surface normal by a vacuum region of 20 Å. The geometry optimization is performed using an 8×8×1 Monhorst-Pack k-point mesh,65 while the electronic structure calculations are carried out with a denser 10×10×1 mesh. After relaxing the system geometries at 0 K, velocity rescaling66-68 is employed to bring the temperature to 300 K. Then, a 6 ps adiabatic MD simulation is performed in the microcanonical ensemble with a 1 fs atomic time-step. The NA couplings are calculated at the Γ-point, since it corresponds to the fundamental bandgap of the three systems. The electron-hole recombination is simulated by the PYXAID code.30, 54 The initial 2000 fs from the 6 ps MD trajectories are used as initial geometries for the NAMD simulations.

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

Figure 1. Simulation cells showing optimized geometries of (a) the Cl-doped MAPbI3 (001) surface, and the same surface with adsorbed (b) pyridine and (c) thiophene molecules. Both molecules connect to the under-coordinated Pb atom; however, the coordinate (dative) bonding is stronger through the nitrogen atom of pyridine than the sulfur atom of thiophene. This difference has a notable influence on localization of the electron wavefunction.

Figure 1 shows the optimized geometries of the three systems under investigation. The Pb-I inorganic cages and the organic MA groups of the slabs maintain their bulk configurations. Adsorption of pyridine and thiophene creates a small perturbation to the topmost surface layer. The coordinate (dative) bonding between the N atom of pyridine and an under-coordinated Pb atom of the Cl-doped MAPbI3 (001) surface is quite strong. In contrast, the interaction of the S atom of thiophene with the Pb atom of the surface is weaker. The calculated binding energies are -0.892 eV and -0.541 eV for pyridine and thiophene, respectively. In combination with the larger electronegativity of N relative to S, the difference in the bonding strength is responsible for the different effects pyridine and thiophene have on the CB edge of the perovskite surface.

8 ACS Paragon Plus Environment

Page 8 of 26

Page 9 of 26

30

CH3NH3 CH 3NH3

(a)

I Cl Pb

20 10 0

PDOS (1/ev)

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

Bare Pyridine

(b)

20

0.3 0.2 0.1

10

0.0 1.0

0

1.5

2.0

Bare Thiophene

(c)

20

0.3 0.2 0.1

10

0.0 1.0

0

-2

1.5

0

2.0

2

Energy (eV) Figure 2. Partial density of states (PDOS) of (a) the bare Cl-doped MAPbI3 (001) surface split into contributions from CH3NH3, I, Cl, and Pb, (b) the same surface (black line) with adsorbed pyridine (red line), and (c) the same surface (black line) with adsorbed thiophene (blue line). The insets in (b) and (c) show magnified PDOS near the conduction band minimum (CBM). The Fermi level is set to zero. Molecular levels are far from the valence band maximum, but close to the CBM.

Figure 2a presents the projected density of states (PDOS) of the bare Cl-doped MAPbI3 (001) surface split into contributions from CH3NH3, I, Cl and Pb. The data demonstrate that the valence band (VB) and CB edges originate primarily from atomic orbitals of I and Pb, respectively. The energies of the orbitals of the organic 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

MA groups reside deep inside the bands. Therefore, these cations have no direct contribution to the NA coupling responsible for the electron-hole recombination. At the same time, they affect the band edge wavefunction through electrostatic interaction. The bandgap of the Cl-doped MAPbI3 (001) surface averaged canonically over the 6ps MD trajectory is 1.840 eV. This amount of energy is deposited into phonon modes during the non-radiative electron-hole recombination. Since phonon quanta available in the system are much smaller than the bandgap, even including C-H stretches of the organic molecules, the non-radiative electron-vibrational relaxation is a multi-phonon process. The PDOS of the Cl-doped MAPbI3 (001) surface with adsorbed pyridine and thiophene molecules, shown in parts (b) and (c) of Figure 2, are separated into the components arising from the surface and the molecules. The insets zoom onto the CB edges. Figure 2b demonstrates that pyridine has states near the CB edge. In comparison, thiophene has very little contribution near the edge. The difference arises due to the larger electronegativity of the N atom of pyridine, compared to the S atom of thiophene. The more electronegative atom lowers the energy of the electronic levels, bringing them closer to the CB edge. The stronger N-Pb coordinative bonding, compared to S-Pb, Figure 1, also contributes to hybridization of perovskite and pyridine orbitals near the CB edge. Neither molecule contributes to the VB edge. Both occupied and unoccupied molecular orbitals (MOs) of pyridine reside lower in pyridine than thiophene. The canonically averaged bandgaps at the Γ point obtained from the 6ps MD trajectory for the surface/pyridine and surface/thiophene systems are 10 ACS Paragon Plus Environment

Page 10 of 26

Page 11 of 26 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

1.788 eV and 1.818 eV, respectively, Table 1. They reflect thermal effects and are most relevant for comparison with experiment. The bandgap reduction is larger for pyridine, which is consistent with the stronger influence of pyridine on the CB edge. The small changes in the bandgap should have negligible influence on electron-hole recombination. In comparison, changes in localization of orbitals near the CB edge should have a more significant effect on the recombination.

Table 1 LUMO-HOMO bandgap, absolute value of NA coupling, pure-dephasing time T2*, and non-radiative electron-hole recombination time for the Cl-doped MAPbI3 (001) surface with and without the Lewis base molecules. The data are canonically averaged over of the 6 ps MD trajectories at the Γ point.

Bare

Bandgap (eV) 1.850

NA coupling (meV) 1.613

Dephasing (fs) 6.7

Recombination (ns) 1.27

Pyridine

1.788

0.765

3.5

9.22

Thiophene

1.818

0.770

3.3

16.6

The electron-hole recombination rate depends strongly on the magnitude of the NA coupling between the initial and final states. In turn, the strength of the coupling depends on the overlap of the corresponding wavefunctions, the sensitivity of the wavefunctions to atomic motions, and the speed of the motions. Since intraband electron and hole relaxation proceeds on a sub-picosecond timescale,9,

69-70

we

consider the NA coupling between the highest occupied MO (HOMO) and lowest unoccupied MO (LUMO), which support relaxed hole and electron. The adsorbed molecules have little effect on the orbitals near the VB edge, Figure S1 of the Supporting Information, because their occupied MOs reside deep inside the VB, Figure 2. The HOMO is created by I atomic orbitals and remains unchanged upon 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

adsorption of the Lewis bases. On the other hand, the unoccupied MOs of the adsorbed molecules are close to the CB edge, Figure 2, and the molecules affect significantly the perovskite LUMO, Figure 3.

Figure 3. LUMO charge density in (a) the bare Cl-doped MAPbI3 (001) surface, and the surface with adsorbed (b) pyridine and (c) thiophene. The LUMO of the bare surface is formed primarily by Pb atoms. Strong dative N-Pb bonding of pyridine to the under-coordinated Pb atom (Figure 1b) localizes the LUMO on pyridine and nearby Pb atoms, reducing the non-adiabatic coupling, accelerating decoherence, and hence, slowing down charge recombination. Thiophene acts as a traditional surface passivating agent by forcing the LUMO away from the surface, which is known to enhance charge recombination.20, 35

Both Lewis bases localize the LUMO, decreasing its overlap with the HOMO, and hence, reducing the NA coupling and the coherence time, Table 1. At the same time, the localization mechanism varies between pyridine and thiophene. Thiophene forces the LUMO away from the surface, Figure 3c. This is the expected behavior of surface passivating species: Unsaturated chemical bonds and defects encountered on surfaces create charge recombination sites, and surface passivation is designed to eliminate such sites and delocalize charges away from the surface. In contrast,

12 ACS Paragon Plus Environment

Page 12 of 26

Page 13 of 26 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

pyridine creates a localized surface state, Figure 3b, because the N atoms makes a strong bond with the perovskite, Figure 1b, and is more electronegative, contributing to the PDOS near the CB edge, Figure 2b. It is particularly important that the surface state created upon binding of pyridine to the perovskite is not a deep electron trap: Its energy is at the perovskite CB edge, and therefore, the electron can easily access the CB and travel to an electrode or a photochemical reaction site. At the same time, the localized LUMO reduces the NA coupling with the HOMO, slowing down the charge recombination. Figure 4 shows Fourier transforms (FTs) of the fluctuations of the HOMO-LUMO energy gaps for the three systems under investigation. Known as influence spectra, they characterize the phonon modes that couple to the electronic degrees of freedom during the non-radiative transition. The vibrational modes induce decoherence in the electronic subsystem and accommodate the excess energy lost during the non-radiative transition. A relatively broad range of phonon modes, arising from the both the I(Cl)-Pb inorganic cage and the organic MA cations, participate in the electron-hole recombination in all three systems. The frequencies extend from around 30 cm-1 to above 400 cm-1. The higher frequency vibrational motions become important in the presence of pyridine and thiophene, and especially pyridine that interacts strongly with the perovskite, Figure 1b. The peaks at 62 cm-1 and 100 cm-1 can be assigned to Pb-I bending and stretching modes, respectively.71 The peaks in the 200-400 cm-1 frequency range are attributed to the torsional motions of the MA cations.71 The LUMO and HOMO of the perovskites originate from the atomic 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

orbitals of the Pb and I atoms, Figures 2a and 3a. The energy levels of the MA cations reside deep inside the CB and VB, and therefore, molecular orbitals of the MA cations do not contribute directly to the electron and hole states near the band edge, Figure 2a. However, they affect the electron and hole wave functions near the band edge by electrostatic interactions, because the cations are charged and carry dipole moments. The Lewis base molecules introduce additional vibrational modes. The new signal at 300 cm-1 is particularly strong in the pyridine system, which bonds to the perovskite via the N-Pb coordinate bond, Figure 1b. Thiophene also extend the range of active vibrational motions to higher frequencies, although these peaks are less significant than for pyridine, Figure 4a, because the thiophene-perovskite bonding is not as strong.

Figure 4

(a) Fourier transforms (FT) of the unnormalized autocorrelation

functions for the HOMO-LUMO gap fluctuations in the Cl-doped MAPbI3(001) surface with and without the adsorbed pyridine and thiophene molecules. (b) Corresponding pure-dephasing functions. A broad range of phonons couple to the electronic subsystem, resulting in fast loss of quantum coherence. Decoherence accelerates in the presence of pyridine and thiophene molecules due to involvement of higher frequency modes. The inset in (b) shows the unnormalized autocorrelation functions of the three systems. The initial values represent the bandgap fluctuation 14 ACS Paragon Plus Environment

Page 14 of 26

Page 15 of 26 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

squared; the bigger the fluctuation, the faster the dephasing.72

The decoherence times were computed as the pure-dephasing times of the optical response theory,73 using the second-order cumulant approximation:20, 33-35, 52



′

 () =  − ħ   ′   ′′   ′′ 

(1)

Here,  () is the unnormalized autocorrelation function (ACF) of the

phonon-induced fluctuation of the energy gap  (), between electronic states i and

j, defined as  () = 〈 (  ) ( −   )〉 

(2)

The influence spectra shown in Figure 4a are obtained as Fourier transforms of the unnormalized autocorrelation ACF shown in the insert of Figure 4b. The pure-dephasing functions, Figure 4b, were fit by Gaussians, exp[-0.5(t/τ)2], and the times, τ, are summarized in Table 1. Phonon-induced loss of electronic coherence is fast, sub-10 fs in the present systems, because they contain a broad spectrum of vibrations stemming from both inorganic and organic parts. Many of these vibrations couple to the electronic subsystem, Figure 4a. The decoherence times are much shorter than the electron-hole recombination times reported in the experiment,26 requiring incorporation of decoherence into the NAMD simulation. Adsorption of the Lewis base molecules accelerates decoherence by a factor of two. Analysis of the second-order cumulant approximation for the optical response function leads one to conclude,72 the fast decoherence can arise either from a rapid decay or a large initial value of the energy gap ACF, shown in the insert in Figure 4b. 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

Since the three ACF decay on similar time scales, the differences in the pure-dephasing times are rationalized by the ACF initial values,  (0), which are equal to the canonically averaged squares of the energy gap fluctuation. The decoherence function is calculated by integrating and exponentiating the unnormalized ACF, Eq. (1), whose initial value,  (0), can be placed in front of the integral. Hence, a larger  (0)

leads to faster decoherence  ( ) . The gap

fluctuations have more significant amplitudes in the systems with the adsorbed molecules due to stronger electron-phonon coupling, as reflected by the taller peaks in the influence spectra, Figure 4a. Decoherence arises from elastic electron-phonon scattering, while inelastic scattering is represented in the NAMD simulation by the NA coupling. The NA matrix elements are smaller for the perovskite surfaces containing pyridine and thiophene than for the bare surface, Table 1, because the NA coupling depends on the wavefunction overlap, which decreases when the LUMOs become more localized in the presence of the molecules, Figure 3.

Figure 5. Electron-hole recombination in the Cl-doped MAPbI3(001) surface with and without the adsorbed Lewis base molecules. Both pyridine and thiophene slow down the non-radiative decay, due to smaller NA coupling and faster loss of coherence, Table 1. The coupling decreases because the electron wavefunctions 16 ACS Paragon Plus Environment

Page 16 of 26

Page 17 of 26 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

become more localized, Figure 3b,c.

Figure 5 demonstrates evolution of the excited state population in the three perovskite systems under consideration. The non-radiative decay times, summarized in Table 1, are obtained using the short-time, linear approximation to the exponential decay, f(t) = exp(-t/τ) ≈ 1-t/τ. The charge recombination takes a long time, over 1 ns in all three systems. This is beneficial for minimizing charge and energy losses, and achieving high photon-to-current conversion efficiencies of perovskite solar cells.4, 6, 9 The calculations demonstrate that the adsorbed Lewis bases reduce electron-hole recombination by an order of magnitude, rationalizing the experimental observation.26 The recombination is slowed down by the adsorbed molecules because the NA coupling and decoherence times are decreased twofold, Table 1. Generally, strong NA coupling and long-lived coherence lead to fast quantum transitions. Weaker coupling and shorter coherence retards the dynamics. The dynamics stops if the coherence time becomes infinitely small, giving rise to the quantum Zeno effect.74-75 Since population transfer requires a buildup of quantum coherence first, suppressing the coherence slows down population dynamics. The NA coupling is decreased and the decoherence is accelerated because the Lewis base molecules change localization of the electron wavefunctions, Figure 3. The dephasing is accelerated since wavefunctions localized in different parts of the system evolve independently and are not correlated. Pyridine and thiophene achieve this effect by different mechanisms. Pyridine forms a strong coordinative bond to the

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

perovskite surface, localizing the wavefunction around this bond. In comparison, thiophene acts by the traditional surface passivation mechanism, in which the charge density is removed from the surface. It is important that the pyridine energy level is still near the perovskite CB edge, Figure 2b.

Even though the pyridine DOS is small

compared to the perovskite DOS, the pyridine state mixes with the perovskite states, such that the LUMO of the combined system becomes delocalized between pyridine and two perovskite surface layers. Such surface state of electron decouples from the hole, leading to reduction of the electron-hole recombination rate. The calculated electron-hole recombination times are shorter than the experimental values,26 most likely because the experimental data account for charge diffusion inside perovskite. Due to computational limitations, the present work employs a small simulation cell, in which the electron and hole are kept artificially close to each other, enhancing the interaction and accelerating the recombination. Nevertheless, the experimental trend is reproduced very well, including both the significant slowing down of the recombination by the Lewis base molecules, and the small differences in the experimental timescales between pyridine and thiophene passivations, with the thiophene system showing a slightly slower recombination.26 The use of a hybrid functional would produce smaller NA coupling76 and better agreement with the experiment.26 However, in order to produce accurate bandgaps,61 hybrid functionals have to be combined with spin-orbit effects,62-63 both of which significantly increase the computational cost. In summary, using NAMD combined with ab initio real-time TDDFT, we 18 ACS Paragon Plus Environment

Page 18 of 26

Page 19 of 26 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

modeled the recent time-resolved experiments26 in order to elucidate the atomistic origin of the observed slowing of charge recombination in the MAPbI3 perovskite upon surface passivation with the Lewis base molecules. The non-radiative electron-hole recombination constitutes a major pathway for charge and energy losses in these materials, limiting the photon-to-electron conversion efficiency. The experimental data show that absorption of the Lewis bases extends the excited-state lifetime and increases the efficiency of the perovskite solar cells. Our simulations demonstrate that the Lewis base molecules change localization of the electron wavefunctions, thereby decreasing the NA coupling. Both pyridine and thiophene achieve the same effect, however, by different mechanisms. Thiophene acts according to the general expectation by forcing electron density away from the surface, which usually acts as a charge recombination site. In comparison, pyridine localizes electron near the surface, because it creates a strong coordinative bond with a surface Pb atom, and since its LUMO resides close to the perovskite CB edge, due to large electronegativity of the N atom. The localized state created by pyridine remains very close to the CB edge, avoiding a deep trap and allowing the charge to escape easily into the band. Both molecules decrease the NA electron-phonon coupling and coherence time by a factor of two, enhancing the computed excited state lifetime by an order of magnitude. A broad spectrum of vibrational modes, including those of both inorganic and organic components, couple to the electronic subsystem and induce rapid coherence loss. The insights obtained during the simulations establish the fundamental principles of charge recombination and reveal an unexpected 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

deceleration mechanism. The experimental trends are rationalized by considering chemical bonding, wavefunction localization, inelastic and elastic electron-vibrational scattering, and electronic coherence. The generated insights can be used to formulate practical guidelines for reduction of non-radiative charge recombination in perovskites via rational choice surface passivation agents. The study advances our understanding of the key factors influencing and controlling the performance of hybrid organic-inorganic perovskite solar cells.

Acknowledgements The first three authors acknowledge support of the National Science Foundation of China (Grant Nos. 21503017, 21573022, 21520102005, 21421003, and 21590801). R.L. is grateful to the Recruitment Program of Global Youth Experts of China, the Beijing Normal University Startup, and the Fundamental Research Funds for the Central Universities. O.V.P. acknowledges support of the U.S. Department of Energy (Grant No. DE-SC0014429).

Supporting Information Available: HOMO and HOMO-1 charge densities of the bare Cl-doped MAPbI3 (001) surface, and the surface with adsorbed pyridine and thiophene. The material is available free of charge via the Internet at http://pubs.acs.org.

20 ACS Paragon Plus Environment

Page 20 of 26

Page 21 of 26 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

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

Page 22 of 26

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.

Seo, G.; Seo, J.; Ryu, S.; Yin, W.; Ahn, T. K.; Seok, S. I., Enhancing the Performance of

Sensitized Solar Cells with PbS/CH3NH3PbI3 Core/Shell Quantum Dots. J. Phys. Chem. Lett. 2014, 5, 2015-2020. 3.

Wehrenfennig, C.; Liu, M.; Snaith, H. J.; Johnston, M. B.; Herz, L. M., Homogeneous Emission

Line Broadening in the Organo Lead Halide Perovskite CH3NH3PbI3–XClX. J. Phys. Chem. Lett. 2014, 5, 1300-1306. 4.

Jeon, N. J.; Noh, J. H.; Yang, W. S.; Kim, Y. C.; Ryu, S.; Seo, J.; Seok, S. I., Compositional

Engineering of Perovskite Materials for High-Performance Solar Cells. Nature 2015, 517, 476. 5.

Christians, J. A.; Manser, J. S.; Kamat, P. V., Multifaceted Excited State of CH3NH3PbI3. Charge

Separation, Recombination, and Trapping. J. Phys. Chem. Lett 2015, 6, 2086-2095. 6.

Pazos-Outón, L. M., et al., Photon Recycling in Lead Iodide Perovskite Solar Cells. Science 2016,

351, 1430-1433. 7.

Yamada, Y.; Yamada, T.; Shimazaki, A.; Wakamiya, A.; Kanemitsu, Y., Interfacial Charge-Carrier

Trapping

in

CH3NH3PbI3-Based

Heterolayered

Structures

Revealed

by

Time-Resolved

Photoluminescence Spectroscopy. J. Phys. Chem. Lett 2016, 7, 1972-1977. 8.

Guo, Z.; Wu, X. X.; Zhu, T.; Zhu, X. Y.; Huang, L. B., Electron-Phonon Scattering in Atomically

Thin 2D Perovskites. ACS Nano 2016, 10, 9992-9998. 9.

Xing, G.; Mathews, N.; Sun, S.; Lim, S. S.; Lam, Y. M.; Grätzel, M.; Mhaisalkar, S.; Sum, T. C.,

Long-Range Balanced Electron- and Hole-Transport Lengths in Organic-Inorganic CH3NH3PbI3. Science 2013, 342, 344-347. 10. Stewart, R. J.; Grieco, C.; Larsen, A. V.; Doucette, G. S.; Asbury, J. B., Molecular Origins of Defects in Organohalide Perovskites and Their Influence on Charge Carrier Dynamics. J. Phys. Chem. C 2016, 120, 12392-12402. 11. Peng, W. N.; Anand, B.; Liu, L. H.; Sampat, S.; Bearden, B. E.; Malko, A. V.; Chabal, Y. J., Influence of Growth Temperature on Bulk and Surface Defects in Hybrid Lead Halide Perovskite Films. Nanoscale 2016, 8, 1627-1634. 12. Wang, B. H.; Wong, K. Y.; Yang, S. F.; Chen, T., Crystallinity and Defect State Engineering in Organo-Lead Halide Perovskite for High-Efficiency Solar Cells. J. Mater. Chem. A 2016, 4, 3806-3812. 13. Tian, Y. X., et al., Enhanced Organo-Metal Halide Perovskite Photoluminescence from Nanosized Defect-Free Crystallites and Emitting Sites. J. Phys. Chem. Lett. 2015, 6, 4171-4177. 14. Cui, P.; Fu, P. F.; Wei, D.; Li, M. C.; Song, D. D.; Yue, X. P.; Li, Y. Y.; Zhang, Z. R.; Li, Y. F.; Mbengue, J. M., Reduced Surface Defects of Organometallic Perovskite by Thermal Annealing for Highly Efficient Perovskite Solar Cells. RSC Adv. 2015, 5, 75622-75629. 15. Kim, J.; Lee, S.-H.; Lee, J. H.; Hong, K.-H., The Role of Intrinsic Defects in Methylammonium Lead Iodide Perovskite. J. Phys. Chem. Lett. 2014, 5, 1312-1317. 16. Li, W.; Liu, J.; Bai, F.-Q.; Zhang, H.-X.; Prezhdo, O. V., Hole Trapping by Iodine Interstitial Defects Decreases Free Carrier Losses in Perovskite Solar Cells: A Time-Domain Ab Initio Study. ACS Energy Lett. 2017, 2, 1270-1278. 17. Kanemitsu, Y., Luminescence Spectroscopy of Lead-Halide Perovskites: Materials Properties and Application as Photovoltaic Devices. J. Mater. Chem. C 2017, 5, 3427-3437.

22 ACS Paragon Plus Environment

Page 23 of 26 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

18. Yang, J.; Siempelkamp, B. D.; Liu, D.; Kelly, T. L., Investigation of CH3NH3PbI3 Degradation Rates and Mechanisms in Controlled Humidity Environments Using in Situ Techniques. ACS Nano 2015, 9, 1955-1963. 19. Christians, J. A.; Miranda Herrera, P. A.; Kamat, P. V., Transformation of the Excited State and Photovoltaic Efficiency of CH3NH3PbI3 Perovskite Upon Controlled Exposure to Humidified Air. J. Am. Chem. Soc. 2015, 137, 1530-1538. 20. Long, R.; Prezhdo, O. V., Dopants Control Electron–Hole Recombination at Perovskite–TiO2 Interfaces: Ab Initio Time-Domain Study. ACS Nano 2015, 9, 11143-11155. 21. Liu, J.; Prezhdo, O. V., Chlorine Doping Reduces Electron–Hole Recombination in Lead Iodide Perovskites: Time-Domain Ab Initio Analysis. J. Phys. Chem. Lett. 2015, 6, 4463-4469. 22. Kim, H.-S.; Im, S. H.; Park, N.-G., Organolead Halide Perovskite: New Horizons in Solar Cell Research. J. Phys. Chem. C 2014, 118, 5615-5625. 23. Nie, W., et al., High-Efficiency Solution-Processed Perovskite Solar Cells with Millimeter-Scale Grains. Science 2015, 347, 522-525. 24. Yin, W.-J.; Yan, Y.; Wei, S.-H., Anomalous Alloy Properties in Mixed Halide Perovskites. J. Phys. Chem. Lett. 2014, 5, 3625-3631. 25. Zheng, F.; Takenaka, H.; Wang, F.; Koocher, N. Z.; Rappe, A. M., First-Principles Calculation of the Bulk Photovoltaic Effect in CH3NH3PbI3 and CH3NH3PbI3-XClX. J. Phys. Chem. Lett. 2015, 6, 31-37. 26. Noel, N. K.; Abate, A.; Stranks, S. D.; Parrott, E. S.; Burlakov, V. M.; Goriely, A.; Snaith, H. J., Enhanced Photoluminescence and Solar Cell Performance Via Lewis Base Passivation of Organic Inorganic Lead Halide Perovskites. ACS Nano 2014, 8, 9815-9821. 27. Shi, S. W.; Li, Y. F.; Li, X. Y.; Wang, H. Q., Advancements in All-Solid-State Hybrid Solar Cells Based on Organometal Halide Perovskites. Materials Horizons 2015, 2, 378-405. 28. Lin, Y., et al., π-Conjugated Lewis Base: Efficient Trap-Passivation and Charge-Extraction for Hybrid Perovskite Solar Cells. Adv. Mater. 2017, 29, 1604545. 29. Jaeger, H. M.; Fischer, S.; Prezhdo, O. V., Decoherence-Induced Surface Hopping. J. chem. phys. 2012, 137, 22A545. 30. 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. 31. Craig, C. F.; Duncan, W. R.; Prezhdo, O. V., Trajectory Surface Hopping in the Time-Dependent Kohn-Sham Approach for Electron-Nuclear Dynamics. Phys. Rev. Lett. 2005, 95, 163001. 32. Fischer, S. A.; Habenicht, B. F.; Madrid, A. B.; Duncan, W. R.; Prezhdo, O. V., Regarding the Validity of the Time-Dependent Kohn–Sham Approach for Electron-Nuclear Dynamics via Trajectory Surface Hopping. J. Chem. Phys. 2011, 134, 024102. 33. Long, R.; Fang, W.; Akimov, A. V., Nonradiative Electron–Hole Recombination Rate Is Greatly Reduced by Defects in Monolayer Black Phosphorus: Ab Initio Time Domain Study. J. Phys. Chem. Lett. 2016, 7, 653-659. 34. Long, R.; Fang, W.; 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. 35. 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

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

Simulation. J. Am. Chem. Soc. 2016, 138, 3884-3890. 36. Habenicht, B. F.; Craig, C. F.; Prezhdo, O. V., Time-Domain Ab Initio Simulation of Electron and Hole Relaxation Dynamics in a Single-Wall Semiconducting Carbon Nanotube. Phys. rev. lett. 2006, 96, 187401. 37. Chaban, V. V.; Prezhdo, V. V.; Prezhdo, O. V., Covalent Linking Greatly Enhances Photoinduced Electron Transfer in Fullerene-Quantum Dot Nanocomposites: Time-Domain Ab Initio Study. J. Phys. Chem. Lett. 2012, 4, 1-6. 38. Long, R.; English, N. J.; Prezhdo, O. V., Photo-Induced Charge Separation across the Graphene– TiO2 Interface Is Faster Than Energy Losses: A Time-Domain Ab Initio Analysis. J. Am. Chem. Soc. 2012, 134, 14238-14248. 39. Long, R.; English, N. J.; Prezhdo, O. V., Defects Are Needed for Fast Photo-Induced Electron Transfer from a Nanocrystal to a Molecule: Time-Domain Ab Initio Analysis. J. Am. Chem. Soc. 2013, 135, 18892-18900. 40. Akimov, A. V.; Neukirch, A. J.; Prezhdo, O. V., Theoretical Insights into Photoinduced Charge Transfer and Catalysis at Oxide Interfaces. Chem. Rev. 2013, 113, 4496-4565. 41. Habenicht, B. F.; Prezhdo, O. V., Nonradiative Quenching of Fluorescence in a Semiconducting Carbon Nanotube: A Time-Domain Ab Initio Study. Phys. Rev. Lett. 2008, 100, 197402. 42. Long, R.; Prezhdo, O. V., Ab Initio Nonadiabatic Molecular Dynamics of the Ultrafast Electron Injection from a PbSe Quantum Dot into the TiO2 Surface. J. Am. Chem. Soc. 2011, 133, 19240-19249. 43. Akimov, A. V.; Prezhdo, O. V., Nonadiabatic Dynamics of Charge Transfer and Singlet Fission at the Pentacene/C60 Interface. J. Am. Chem. Soc. 2014, 136, 1599-1608. 44. Long, R.; Prezhdo, O. V., Instantaneous Generation of Charge-Separated State on TiO2 Surface Sensitized with Plasmonic Nanoparticles. J. Am. Chem. Soc. 2014, 136, 4343-4354. 45. Hyeon-Deuk, K.; Prezhdo, O. V., Time-Domain Ab Initio Study of Auger and Phonon-Assisted Auger Processes in a Semiconductor Quantum Dot. Nano Lett. 2011, 11, 1845-1850. 46. Long, R.; Prezhdo, O. V., Asymmetry in the Electron and Hole Transfer at a Polymer–Carbon Nanotube Heterojunction. Nano Lett. 2014, 14, 3335-3341. 47. Long, R.; Prezhdo, O. V., Time-Domain Ab Initio Analysis of Excitation Dynamics in a Quantum Dot/Polymer Hybrid: Atomistic Description Rationalizes Experiment. Nano Lett. 2015, 15, 4274-4281. 48. Hyeon-Deuk, K.; Prezhdo, O. V., Multiple Exciton Generation and Recombination Dynamics in Small Si and CdSe Quantum Dots: An Ab Initio Time-Domain Study. ACS Nano 2012, 6, 1239-1250. 49. Liu, J.; Prezhdo, O. V., Chlorine Doping Reduces Electron-Hole Recombination in Lead Iodide Perovskites: Time-Domain Ab Initio Analysis. J. Phys. Chem. Lett. 2015, 6, 4463-4469. 50. Li, W.; Liu, J.; Bai, F. Q.; Zhang, H. X.; Prezhdo, O. V., Hole Trapping by Iodine Interstitial Defects Decreases Free Carrier Losses in Perovskite Solar Cells: A Time-Domain Ab Initio Study. ACS Energy Lett. 2017, 2, 1270-1278. 51. 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. 52. 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. 53. Jankowska, J.; Prezhdo, O. V., Ferroelectric Alignment of Organic Cations Inhibits Nonradiative Electron-Hole Recombination in Hybrid Perovskites: Ab Initio Nonadiabatic Molecular Dynamics. J.

24 ACS Paragon Plus Environment

Page 24 of 26

Page 25 of 26 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

Phys. Chem. Lett. 2017, 8, 812-818. 54. Akimov, A. V.; Prezhdo, O. V., The Pyxaid Program for Non-Adiabatic Molecular Dynamics in Condensed Matter Systems. J. Chem. Theory Comput. 2013, 9, 4959-4972. 55. Haruyama, J.; Sodeyama, K.; Han, L.; Tateyama, Y., Termination Dependence of Tetragonal CH3NH3PbI3 Surfaces for Perovskite Solar Cells. J. Phys. Chem. Lett. 2014, 5, 2903-2909. 56. Huang, J.; Tan, S.; Lund, P. D.; Zhou, H., Impact of H2O on Organic-Inorganic Hybrid Perovskite Solar Cells. Energy Environ. Sci. 2017, 10, 2284-2311. 57. Kresse, G.; Furthmuller, J., Efficient Iterative Schemes for Ab Initio Total-Energy Calculations using a Plane-Wave Basis Set, Phys. Rev. B 1996, 54, 11169. 58. Perdew, J. P.; Burke, K.; Ernzerhof, M., Generalized Gradient Approximation Made Simple. Phys. Rev. Lett. 1996, 77, 3865. 59. Blöchl, P. E., Projector Augmented-Wave Method. Phys. Rev. B 1994, 50, 17953. 60. Grimme, S., Semiempirical GGA-Type Density Functional Constructed with a Long-Range Dispersion Correction. J. Comput. Chem. 2006, 27, 1787-1799. 61. Chen, Q.; De Marco, N.; Yang, Y.; Song, T.-B.; Chen, C.-C.; Zhao, H.; Hong, Z.; Zhou, H.; Yang, Y., Under the Spotlight: The Organic–Inorganic Hybrid Halide Perovskite for Optoelectronic Applications. Nano Today 2015, 10, 355-396. 62. Even, J.; Pedesseau, L.; Jancu, J.-M.; Katan, C., Importance of Spin–Orbit Coupling in Hybrid Organic/Inorganic Perovskites for Photovoltaic Applications. J. Phys. Chem. Lett. 2013, 4, 2999-3005. 63. Brivio, F.; Butler, K. T.; Walsh, A.; van Schilfgaarde, M., Relativistic Quasiparticle Self-Consistent Electronic Structure of Hybrid Halide Perovskite Photovoltaic Absorbers. Phys. Rev. B 2014, 89, 155204. 64. Umari, P.; Mosconi, E.; De Angelis, F., Relativistic GW Calculations on CH3NH3PbI3 and CH3NH3SnI3 Perovskites for Solar Cell Applications. Sci. Rep-UK. 2014, 4, 4467. 65. Monkhorst, H. J.; Pack, J. D., Special Points for Brillouin-Zone Integrations. Phys. Rev. B 1976, 13, 5188. 66. Tully, J. C., Molecular Dynamics with Electronic Transitions. J. Chem. Phys. 1990, 93, 1061-1071. 67. Hammes-Schiffer, S.; Tully, J. C., Proton Transfer in Solution: Molecular Dynamics with Quantum Transitions. J. Chem. Phys. 1994, 101, 4657-4667. 68. Sifain, A. E.; Wang, L.; Prezhdo, O. V., Communication: Proper Treatment of Classically Forbidden Electronic Transitions Significantly Improves Detailed Balance in Surface Hopping. J. Chem. Phys. 2016, 144, 211102. 69. Madjet, M. E.-A.; Akimov, A. V.; El-Mellouhi, F.; Berdiyorov, G. R.; Ashhab, S.; Tabet, N.; Kais, S., Enhancing the Carrier Thermalization Time in Organometallic Perovskites by Halide Mixing. Phys. Chem. Chem. Phys. 2016, 18, 5219-5231. 70. Piatkowski, P.; Cohen, B.; Javier Ramos, F.; Di Nunzio, M.; Nazeeruddin, M. K.; Gratzel, M.; Ahmad, S.; Douhal, A., Direct Monitoring of Ultrafast Electron and Hole Dynamics in Perovskite Solar Cells. Phys. Chem. Chem. Phys. 2015, 17, 14674-14684. 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. Akimov, A. V.; Prezhdo, O. V., Persistent Electronic Coherence Despite Rapid Loss of Electron– Nuclear Correlation. J. Phys. Chem. Lett. 2013, 4, 3857-3864.

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

73. Mukamel, S., Principles of Nonlinear Optical Spectroscopyoxford University Press. New York 1995. 74. Prezhdo, O. V., Quantum Anti-Zeno Acceleration of a Chemical Reaction. Phys. Rev. Lett. 2000, 85, 4413-4417. 75. Kilina, S. V.; Neukirch, A. J.; Habenicht, B. F.; Kilin, D. S.; Prezhdo, O. V., Quantum Zeno Effect Rationalizes the Phonon Bottleneck in Semiconductor Quantum Dots. Phys. Rev. Lett. 2013, 110, 180404. 76. Lin, Y.; Akimov, A. V., Dependence of Nonadiabatic Couplings with Kohn–Sham Orbitals on the Choice of Density Functional: Pure vs Hybrid. The J. Phys. Chem. A 2016, 120, 9028-9041.

26 ACS Paragon Plus Environment

Page 26 of 26