TiO2 Interface Facilitates Ultrafast

Feb 7, 2017 - Yaqing WeiZhaohui ZhouWei-Hai FangRun Long. The Journal ... Zhaosheng ZhangWei-Hai FangMarina V. TokinaRun LongOleg V. Prezhdo...
0 downloads 0 Views 1MB Size
Subscriber access provided by UB + Fachbibliothek Chemie | (FU-Bibliothekssystem)

Article 2

Strong Interaction at the Perovskite/TiO Interface Facilitates Ultrafast PhotoInduced Charge Separation: A Nonadiabatic Molecular Dynamics Study Run Long, Wei-Hai Fang, and Oleg V. Prezhdo J. Phys. Chem. C, Just Accepted Manuscript • DOI: 10.1021/acs.jpcc.6b12921 • Publication Date (Web): 07 Feb 2017 Downloaded from http://pubs.acs.org on February 8, 2017

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 free 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 accessible to all readers and 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 C 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

Strong Interaction at the Perovskite/TiO2 Interface Facilitates Ultrafast Photo-Induced Charge Separation: A Nonadiabatic Molecular Dynamics Study Run Long,1* Wei-Hai Fang,1 Oleg V. Prezhdo2† 1

College of Chemistry, Key Laboratory of Theoretical & Computational Photochemistry of Ministry of Education, Beijing Normal University, Beijing, 100875, P. R. China

2

Department of Chemistry, University of Southern California, Los Angeles, CA 90089, USA

ABSTRACT: Interfacial electron transfer (ET) plays a key role in the operation of solar cells based on TiO2 sensitized with organohalide perovskites, since it leads to separation of the photo-generated electrons and hole into different materials. The reported experimental ET times range by over three orders of magnitude, from sub-200 fs to over 300 ps. Using nonadiabatic molecular dynamics combined with ab initio time-domain density functional theory, we demonstrate that ET at a CH3NH3PbI3/TiO2 interface can be complete within 100 fs, indicating that the longer time scales reflect other processes, such as charge and exciton diffusion in perovskite bulk. The electron injection is fast because the interaction between the donor and acceptor species is strong at ambient conditions. Photoexcitation directly at the interface can create a charge separated state. Electrons generated farther away inject by a combination of adiabatic and nonadiabatic mechanisms. Thermally activated low frequency vibrational motions at the interface modulate the CH3NH3PbI3/TiO2 separation, creating opportunities for chemical bonding and generating channels for adiabatic ET. Higherfrequency modes create large nonadiabatic coupling. The interaction between perovskite and TiO2 is purely van der Waals at 0 K, whereas at ambient temperatures I-Ti covalent bonds

* †

Corresponding author, E-mail: [email protected] Corresponding author, E-mail: [email protected]

1 ACS Paragon Plus Environment

The Journal of Physical Chemistry

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

can form transiently at the interface. The covalent bonding is particularly important for photo-excitation of charge-separated states and adiabatic ET. The ET occur prior to nonradiative electronic energy losses that can lead to charge trapping and recombination. The ultrafast interfacial charge separation contributes to the high efficiencies of perovskitesensitized TiO2 solar cells. The reported simulations provide a detailed time-domain atomistic description of the interfacial ET and advance our understanding of carrier dynamics in perovskite solar cells.

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

1. INTRODUCTION Bringing

ideas

from

dye-sensitized

solar

cells,1

organohalide

perovskite

(CH3NH3PbI3)-sensitized TiO2 solar cells have been attracting intense attention in the past several years due to its rapidly increasing power conversion efficiency: starting at 3.8 % of the first perovskite solar cells,2 and reaching3-7 and exceeding 20%8-10 within five years. Low charge recombination rates, high charge-carrier mobilities, and long electron and hole diffusion lengths are believed to be responsible for the high performance of perovskite solar cells.11-28 In particular, the excited state lifetime in the CH3NH3PbI3 perovskite can last from several nanoseconds29 to several microseconds,30 and it further increases in the CH3NH3PbI3/TiO2 system.31-33 Slow electron-hole recombination reduces electronic energy dissipation into phonons and enhances photon-to-electron conversion efficiency. Realistic factors such as defects and environment, such as moisture, have a complicated influence on the electron-hole recombination. Local microstructures decrease excited electronic state lifetimes and deteriorate solar cells quality, while surface passivation can have a significant positive effect.34-35 Rational choice of dopants can control electron-hole recombination in pristine CH3NH3PbI336 and at the CH3NH3PbI3/TiO2 interface.37 Low humidity delays electron-hole recombination in CH3NH3PbI3 perovskites, while heavy humidity accelerates the recombination.38 The photo-induced interfacial electron transfer (ET) constitutes the initial step of conversion of solar energy into electrical or chemical energy, leading to generation of free charge carriers. Rapid charge separation is usually favorable, since separated charges cannot recombine easily. Further, if charges can separate before relaxing to their lowest energy states, they are less likely to get trapped at interfacial defects. Figure 1 shows the energy diagram of the CH3NH3PbI3/TiO2 interface. Photoexcitation of CH3NH3PbI3 promotes an electron from its ground state into its excited state that is located energetically within the TiO2 conduction 3 ACS Paragon Plus Environment

The Journal of Physical Chemistry

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

band (CB). The electron is subsequently injected into the TiO2 CB. Since the CB minimum (CBM) of perovskite CH3NH3PbI3 is higher than the CBM of TiO2, the electron can be efficiently injected into the TiO2 CB even if it relaxes rapidly down to the band edge of perovskite. This is because the excited electron requires several nanoseconds to recombine with the hole, needing to pass across a wide energy gap.31-33 The non-radiative electron-hole recombination is a high order process requiring excitation of many phonons.

Figure 1. Energy diagram for the photoinduced electron injection from CH3NH3PbI3 into a TiO2 surface. An absorbed photon promotes an electron from the valence band (VB) to the conduction band (CB) in CH3NH3PbI3, putting it in resonance with the TiO2 CB. Then, the excited electron is injected into the TiO2 CB.

Interestingly, experiments showed a variety of time scales of electron injection from CH3NH3PbI3 into TiO2, ranging from sub-200 fs to 300 ps, as measured by different techniques.39-43 Zhu et al. reported the injection times of 260-307 ps.41 Insertion of a layer of ultrathin graphene quantum dots between perovskite and TiO2 accelerates ET to 90-106 ps.41 Using time-resolved terahertz and microwave conductivity measurements, Ponseca and coauthors suggested that ET from CH3NH3PbI3 into TiO2 occurs within 2 ps.40 These ET time scales are several orders of magnitude longer than intraband carrier cooling.24 Sum et al. 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

pointed out that band bending due to strong interfacial dipoles may result in a barrier for ET between CH3NH3PbI3 and TiO2.44 The longer ET times can be attributed to electron diffusion inside perovskite, rather to the interfacial ET itself. More recently, Piatkowski and coauthors reported the photoinduced ultrafast ET from CH3NH3PbI3 into TiO2 to occur within 200 fs.39 Such rapid ET competes successfully with the electron-phonon energy relaxation inside CH3NH3PbI3,26 demonstrating that hot electrons can be efficiently extracted from the photosensitizer. The broad range of the experimental ET times reported for the CH3NH3PbI3-TiO2 system motivates theoretical investigations, which can separate the interfacial ET process from charge diffusion to and from the interface. By mimicking most directly the timeresolved experiments, ab initio time-domain simulations also establish the ET mechanisms, which have qualitatively different dependence on system properties, such as donor-acceptor coupling and local densities of states. The simulations identify the phonon modes that participate in the ET process. An ab initio atomistic description captures the realistic aspects of the CH3NH3PbI3/TiO2 interface, including the possibilities of chemical interactions, formation of defects, thermally induces disorder, charge trapping, etc. As demonstrated with the studies of TiO2 surfaces sensitized with molecule chromophores,45-51 water,52 semiconducting53-55 and metallic56 nanoparticles, and graphene,57 an atomistic description complements phenomenological models and provides many valuable and often critical insights into the photoinduced electron-vibrational dynamics.58 Presently, we report simulations of the photoexcited ET and related processes at the CH3NH3PbI3/TiO2 interface, using a combination of nonadiabatic molecular dynamics (NAMD) with real-time time-domain density functional theory (TDDFT).59-60 The simulations demonstrate that the photo-excited (PE) electron can be efficiently extracted by TiO2 within 100 fs, competing successfully with electron-phonon relaxation inside the 5 ACS Paragon Plus Environment

The Journal of Physical Chemistry

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 6 of 30

CH3NH3PbI3. The ET is ultrafast due to strong interaction between CH3NH3PbI3 and TiO2. The strong donor-acceptor coupling creates the possibility of optical excitation of chargeseparated states, if photoexcitation occurs near the interface. PE electrons generated farther from the interface inject with nearly equal probability by both adiabatic and nonadiabatic mechanisms. Several scenarios for individual ET events can be observed.

Thermally

activated low-frequency vibrational motions modulate chemical bonding at the interface and affect the donor-acceptor coupling strength. Higher frequency modes create large nonadiabatic coupling. The strong chemical-like donor-acceptor interaction at the CH3NH3PbI3/TiO2 interface favors ultrafast ET, which contributes to the observed high power conversion efficiency of solar cells based on perovskite sensitized TiO2.

2. METHODS The simulations of ET dynamics are performed using the mixed quantum-classical approach61 implementing NAMD59,

62-64

within TDDFT65-67 in the Kohn-Sham (KS)

representation.68 The faster and lighter electrons are treated quantum mechanically, while the slower and heavier nuclei are classical. The computational procedure is outlined below. A detailed description of this approach can be found in ref.59-60, 69-70

2.1 Time-Dependent Kohn-Sham Theory for Electron-Nuclear Dynamics The ET dynamics including the nonadiabatic effects71 are described by real-time TDDFT59-60 within the Kohn-Sham (KS)68 framework. The electron density, ,  , is expressed by the sum of the densities of the occupied time-dependent single-electron KS orbitals, , , Ne

ρ (r, t ) = ∑ | ϕ p (r, t ) |2

(1)

p =1

6 ACS Paragon Plus Environment

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

The evolution of the electron density is determined by the time-dependent variational principle, leading to a set of single-electron equations for the evolution of the KS orbitals:

ih

∂ϕ p (r, t ) ∂t

= H (r, R, t )ϕ p (r, t ); p = 1,2,..., Ne

(2)

The single-electron equations are coupled, because the Hamiltonian H depends on the electron density. The electron-vibrational interaction enters the Hamiltonian H through the external potential generated by the atomic motions. The time-dependent single-electron orbitals in eq 2 are expressed in the basis of the adiabatic KS orbitals, ϕ% p (r, R (t )) , which are calculated for the current atomic positions R. The adiabatic representation of the time-dependent KS orbital occupied by the PE electron is given in eq 3

ϕ PE (r , t ) = ∑ ck (t )ϕ%k (r; R (t ))

(3)

k

Inserting eq 3 into eq 2 gives the equation for evolution of the expansion coefficients.

ih

∂ c j (t ) = ∑ ck (t )(ε kδ jk + d jk ) ∂t k

(4)

where ε k is the energy of the adiabatic state k, and d jk is the nonadiabatic coupling between states k and j. The nonadiabatic coupling is created by atomic motions and reflects the strength of electron-vibrational interaction. It is calculated numerically as the overlap of orbitals j and k at sequential time steps

d jk = −ih < ϕ% j | ∇ R | ϕ%k > • ≈−

dR ∂ = −ih < ϕ% j | | ϕ%k > dt ∂t

ih (< ϕ% j (t ) | ϕ%k (t + ∆t ) > − < ϕ% j (t + ∆t ) | ϕ%k (t ) >) 2∆t

7 ACS Paragon Plus Environment

(5)

The Journal of Physical Chemistry

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 8 of 30

2.2 Electron Transfer Mechanisms: Adiabatic, Nonadiabatic and Direct by Photoexcitation The extent of ET from CH3NH3PbI3 into the TiO2 surface is computed by integrating the PE electron density over the region of the simulation cell occupied by CH3NH3PbI3, Figure 2.



CH3 NH3 PbI3

ρPE (r, t )dr = ∫

CH3 NH3 PbI3

= ∑ ck* (t )c j (t )∫

| ϕPE (r, t ) |2 dr

CH3 NH3 PbI3

k, j

ϕ%k* (r, R(t ))ϕ% j (r, R(t ))dr

(6)

The time-derivative of eq 6 gives the contributions of the adiabatic and nonadiabatic ET components: d ∫ CH 3 NH 3 PbI 3 ρ PE (r , t ) d r dt

=

* d ∫ CH 3 NH 3 PbI 3 φi*φ j d r   d ( ci c j ) * * φ φ + d r c c  i j i j ∫CH 3 NH 3 PbI3 dt  dt 

∑ k, j

The first term has fixed localizations of adiabatic states, ∫

CH 3 NH 3 PbI 3

φi*φ j dr , but changing

expansion coefficients, ci*c j . The second term has fixed adiabatic state occupations, but changing localizations. The first and second terms represent nonadiabatic ET and adiabatic ET, respectively. The adiabatic ET proceeds by a change in the localization of the adiabatic state occupied by the electron from CH3NH3PbI3 to TiO2, induced by atomic motions. To undergo a nonadiabatic transfer, the PE electron has to hop into a TiO2 state, causing a change in the state occupations. If the donor-acceptor coupling is very strong, photoexcitation can directly promote an electron from the donor state to an acceptor state, corresponding to the direct ET mechanism. Typically, nonadiabatic ET operates when the donor-acceptor coupling is weak , whereas adiabatic and direct ET require strong coupling. Adiabatic and nonadiabatic ET are typically described by different analytic expressions, in particular, the Marcus theory and the Fermi’s golden rule, respectively. Therefore, establishing the ET mechanism is of both fundamental and practical importance. A more 8 ACS Paragon Plus Environment

(7)

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

detailed discussion of the ET mechanisms can be found in our recent work49, 53, 56-57 and reviews.70, 72

2.3 Characterization of Active Phonon Modes: Autocorrelation Function and Fourier Transform Analysis of fluctuations of the electronic energy  caused by atomic motions allows one to characterize the phonon modes that couple to the electronic subsystem, promote charge transfer and are responsible for electron-phonon energy relaxation. The autocorrelation function (ACF), defined as (8)

 = 〈0〉 

indicates how many vibrational modes are coupled to the electronic degrees of freedom. Here, the brackets indicate canonical averaging. A rapidly decaying ACF suggests participation of many vibrations, while a periodically oscillating ACF shows that few vibrations influence the electronic subsystem. The energy ACF describes how the energy at a particular time depends on its value at earlier times. Generally, poorly correlated random motions make ACFs decay rapidly from 1 to 0. Well-correlated periodic vibrations result in ACFs that oscillate between 1 and −1. Often, the ACF is normalized (9)

  = 〈0〉 / 〈  0〉

by its initial value 0 = 〈  0〉 . The square root of this value gives the average energy fluctuation. Replacing the energy by localization in eqs. 8 and 9, gives ACF of the PE state localization, which is defined in eq. 6. Fourier transforms (FT) of the energy and localization ACFs characterize the phonon modes that couple to the electronic subsystem. The modes are identified by their frequencies. 9 ACS Paragon Plus Environment

The Journal of Physical Chemistry

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 electron-phonon coupling strength for the mode of a particular frequency is reflected in the magnitude of the FT signal at that frequency

2.4 Simulation Details The simulation cell, shown in Figure 2, contains a 150-atom anatase TiO2-(5 × 5) (001) surface interfaced with 108-atom CH3NH3PbI3 (3×3)(100) surface terminated with the organic cation CH3NH3. We also considered a CH3NH3PbI3/TiO2 interface terminated with PbI2, and the calculation showed it to be unstable.37 The lattice mismatch between the two slabs is only 0.01Å in both x and y direction. The system is separated from its image along the surface normal by a vacuum region of 20 Å. Additional calculations are performed with the simulation cell having a twice thicker perovskite layer, Figure S1 of Supporting Information (SI). The interfacial interaction between the perovskite and TiO2 subsystems includes the long-range van der Waals force described by the DFT-D2 method of Grimme.73 The geometry optimization, adiabatic MD, and nonadiabatic couplings are obtained using the Vienna ab initio simulation package (VASP).74 The electronic exchange and correlation energies are described by the Perdew-Burke-Ernzerhof (PBE) functional75 under generalized gradient approximation. The valence electron-ion interaction is treated by the projector augment wave method.76 It is important to note that spin-orbit interactions affect significantly the electronic structure of CH3NH3PbI3 due to presence of the heavy Pb atoms. It has been shown77-78 that DFT calculations using PBE and other generalized gradient approximation (GGA) functionals without spin-orbit interactions give a good agreement with the experimental band gap of the CH3NH3PbI3 perovskite,11,

24, 79

allowing one the use of these cost-efficient

functionals for exploring novel effects. 10 ACS Paragon Plus Environment

Page 10 of 30

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

After relaxing the geometry at 0 K, the CH3NH3PbI3/TiO2 system is heated up to 300 K using repeated velocity rescaling, corresponding to the temperature in the experiment.39 After that, a 3 ps adiabatic MD simulation is performed in the microcanonical ensemble with a 1 fs atomic time-step. The adiabatic state energies and nonadiabatic couplings are calculated for each step of the MD run. To simulate the photoinduced ET dynamics, 500 initial system geometries are selected randomly from the 3 ps adiabatic MD trajectory. The KS orbital corresponding to the PE state is chosen by considering optically active transitions near the perovskite band gap. In particular, excitations between pairs of KS orbitals from the valence and conduction are selected within the energy range from the perovskite band gap to 0.25 eV above it, and the excitation with the largest transition dipole moment is found. The energy fluctuation of the PE state is on the order of 100 meV. The inhomogeneous distribution of the photoexcited state energies due to thermal fluctuations is significantly smaller than the perovskite and the TiO2 band gaps; however, it has an influence on the position of the donor state with respect to the acceptor state due to frequent energy levels crossings. An electron is placed into the PE state orbital, and its evolution is tracked by solving the TDKS equation, eq 4, using the second-order differencing scheme and a 10-3 fs electronic time-step. This procedure is applied to study the ET dynamics within the manifold of CH3NH3PbI3 and TiO2 conduction band states.

3. RESULTS AND DISCUSSION 3.1 Geometric Structure Figure 2 presents the top (top panels) and side (bottom panels) views of the optimized geometry at 0 K (left panels) and a typical geometry form the MD simulation at 300 K (right panels). At 0 K, the TiO2 slab remains nearly intact, and all CH3NH3+ cations point in the 11 ACS Paragon Plus Environment

The Journal of Physical Chemistry

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 12 of 30

same direction, confirmed by both the side and top views. The distances between the I and Pb atoms and the cation orientations determine the globally stable geometry. The interaction between TiO2 and CH3NH3PbI3 is purely van der Waals. At ambient temperature, both CH3NH3PbI3 and TiO2 geometries change significantly. In particular, thermal fluctuations allow formation of several I-Ti covalent bonds at the interface. The chemical bonding enhances the donor-acceptor coupling and accelerates ET. The largest movement on the TiO2 surface is associated with surface bridging oxygen atoms approaching the CH3NH3PbI3 slab (side view). The motion of the bridging oxygens acts to reduce the donor-acceptor separation, and further enhances the donor-acceptor coupling. The organic groups do not directly participate in the interfacial bonding, rather, they affect the donor-acceptor interaction electrostatically. Generally, strong covalent bonds result in adiabatic ET,49,

53

while weak van der

Waals interaction leads to nonadiabatic ET.56-57 For instance, the adiabatic ET mechanism is dominant in a broad range of systems including TiO2 surface sensitized with the molecular chromophore49 and the semiconducting quantum dot53 due to covalent bonds. The nonadiabatic ET mechanism is dominant in the graphene/TiO257 and Au20/TiO256 systems, which predominantly interact by van der Waals forces. If the donor-acceptor interaction is sufficiently strong, the PE electron can be promoted directly from donor to acceptor as the photon is being absorbed. In the present CH3NH3PbI3/TiO2 system, the interfacial I-Ti covalent bonds and reduced donor-acceptor separation facilitate charge delocalization and are beneficial for both direct and adiabatic ET. The simulation cell with a thicker perovskite layer is shown in Figure S1 of SI. The strong perovskite/TiO2 interaction is seen at room temperature, similarly to that in the smaller simulation cell, Figure 2.

12 ACS Paragon Plus Environment

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

Figure 2. Top and side views of the simulation cell showing geometry of the CH3NH3PbI3 (100)-TiO2 anatase (001) interface at 0 K (left panels) and 300 K (right panels). Thermal fluctuations influence the system geometry and, hence, its electronic structure. A larger simulation cell is shown in Figure S1.

3.2 Electronic Structure The electronic interaction between the donor and acceptor subsystems mixes their orbital densities. Generally, the stronger is the interaction, the more significant is the mixing. Photoexcitation of an electron from a localized donor state into a state that exhibits a significant delocalization onto the acceptor results in partial ET by the direct photoexcitation mechanism. In the present system, both donor and acceptor energy levels are continuous bands, the PE state is delocalized between CH3NH3PbI3 and TiO2, and the direct mechanism contributes significantly to the overall ET event. Figure 3 presents the charge densities of the donor and acceptor levels. The electron donor state is delocalized significantly between CH3NH3PbI3 and TiO2 from, Figure 3a. 13 ACS Paragon Plus Environment

The Journal of Physical Chemistry

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

Within the CH3NH3PbI3 slab, the electron density is mostly localized on the surface Pb atoms and interfacial I atoms. The mixing is extremely strong at the covalent I-Ti bonds regions. In particular, about half of the PE electron is delocalized into TiO2 slab, indicating that the direct ET mechanism is responsible for nearly 50% of the overall injection event. The electron acceptor state is mostly localized within the TiO2 slab, and the tail of the wave function extends onto interfacial the I-Ti covalent bonding region, Figure 3b. Direct chemical-like bonding between CH3NH3PbI3 and TiO2 produces a strong donor-acceptor coupling, resulting in mixing of the donor and acceptor orbital densities.

Figure 3. Charge densities of (a) donor and (b) acceptor states of the CH3NH3PbI3-TiO2 system. Not only the photoexcited donor state is delocalized significantly into the TiO2 slab, but also the final acceptor state at the TiO2 conduction band edge is delocalized into the interfacial region. The latter fact is unexpected since the TiO2 conduction band edge resides in the CH3NH3PbI3 band gap. The wave function delocalization between the donor and acceptor materials leads to rapid electron transfer from CH3NH3PbI3 into TiO2. The projected density of states (PDOS) of the combined system, separated into the contribution from CH3NH3PbI3 (red line) and TiO2 (black line), is shown in the top panel of Figure 4. The perovskite band gap calculated for the optimized geometry of around 1.5 eV in agreement with the experimental value11, 24, 79 and previous DFT calculations.77-78 Doubling the thickness of the perovskite layer gives similar band gap and PDOS, compare Figure S2 of 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

SI with top panel of Figure 4. These results indicate that our simulation cell is able to capture the key electronic properties of the perovskite/TiO2 interface. The PDOS illustrates that that the CH3NH3PbI3 CB is well inside the TiO2 CB. The high density of TiO2 acceptor states facilitates efficient ET and favours the nonadiabatic ET mechanism. The bottom panel of Figure 4 shows the distribution of PE state energy and its localization. The two properties exhibit no correlation.

Figure 4. Top panel: Projected density of states (PDOS) of the CH3NH3PbI3 (100) / TiO2 anatase (001) interface, separated into contributions from CH3NH3PbI3 (red line) and TiO2 (black line). The Fermi level is set to zero. Bottom Panel: TiO2 PDOS in the excitation

15 ACS Paragon Plus Environment

The Journal of Physical Chemistry

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

energy range (black line) and localization of the photoexcited state on CH3NH3PbI3 (red squares). About 50% of the photoexcited state density delocalizes onto TiO2 on the average. The PE state localization varies little around 0.5 and remains almost constant in the excitation energy region (red dots), indicating that photoexcitation of CH3NH3PbI3 creates a state that is already about 50% delocalized onto the TiO2 acceptor. The local dip in the DOS correlates to the t2g–eg splitting of the Ti 3d orbitals that form the TiO2 CB.80 In the region of decreased DOS, the CH3NH3PbI3 interacts with fewer TiO2 CB states, and therefore, the interface properties become less averaged. Higher TiO2 DOS tends to average out the properties of the interface. The differences between the average behaviour and individual members of the canonical ensemble of CH3NH3PbI3/TiO2 systems are emphasized further in the ET dynamics section.

3.3 Electronic-Vibrational Interactions The atomic motions affect the energy and localization of the PE state, create an ensemble of initial conditions for the photoinduced ET, and to a large extent, determine whether the ET proceeds by the adiabatic or nonadiabatic mechanism. Nuclear fluctuations generate the nonadiabatic coupling, causing hops between adiabatic states: the nonadiabatic coupling is directly proportional to the nuclear velocity dR/dt. Nuclear motions also influence localization of adiabatic state densities: a shift of the density of the state occupied by the electron from donor to acceptor gives adiabatic ET. FTs of the evolving PE state energy and localization are shown in Figure 5. The normalized ACFs of the PE state energy and localization are presented in the insets. The FTs characterize the phonon modes that couple to electronic subsystem and promote the ET. The coupling is directly correlated to the second derivative of the energy along the nuclear 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

trajectory.81 Therefore, the vibrational modes that most strongly modulate the energy levels also create a large coupling.

Figure 5. Fourier transforms of phonon-induced fluctuation of the photoexcited state (a) energy and (b) localization. The insets show autocorrelation functions of these fluctuations. The FTs show that low-frequency vibrational modes couple to the electronic subsystems more strongly than high-frequency modes, Figure 5. Phonon modes with frequencies less than 200 cm-1 dominate both PE state energy and localization spectra. In addition to the largest amplitude at 100 cm-1, notable contributions from vibrations between 400 and 800 cm-1 are seen in the localization spectral density. The strongest peak at 100 cm-1 seen in both energy and localization spectra can be assigned to the I-Pb stretching mode at 94 cm-1.82 The peak can also be attributed to the librations of the organic cations at 119 cm-1.82 The organic group cations contribute to the ET by affecting the donor-acceptor interaction

17 ACS Paragon Plus Environment

The Journal of Physical Chemistry

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

electrostatically. Rotation of organic groups brings notable changes of the electrostatic field.83-84 The small peak at 200 cm-1 can be assigned to the Eg mode of TiO2 at 198 cm-1.85 A more accurate mode assignment requires a full normal mode analysis, which is quite expensive for systems of this size. Note that not only normal modes themselves but also electron-phonon coupling for each mode would need to be obtained. The phonon modes seen in the power spectra of the PE localization show significant differences from those of the PE state energy, in particular in the middle to high-frequency region. The localization is sensitive to changes in the wave function, while energy is computed by averaging over the entire wave function. Therefore, the localization shows many more frequencies in the power spectrum. Shown in Figure 5b, the peaks at 400 cm-1, 500 cm-1 and 615 cm-1 can be assigned to the 409 cm-1, 515 cm-1 and 633 cm-1 modes of anatase TiO2,85-86 respectively. The higher frequencies above 600 cm-1 are associated with fast motions of surface oxygen atoms. Because the electron donor state is shared by the CH3NH3PbI3 sensitizer and the TiO2 substrate, Figure 3a, its wave function is more sensitive to low-frequency phonons involving both subsystems. High-frequency CH and NH stretching modes of CH3NH3 molecules are not detected in the FTs, because these moieties do not contribute to the donor wave function, Figure 3a. As shown in the insets of Figure 5, the PE state energy ACF decays quite slowly. It reaches zero at approximately 100 fs and then rebounds to 30% of its initial magnitude at 300 fs. In comparison, the PE state localization ACF decays much faster and shows no significant rebounds. The difference is rationalized by the difference in the power spectra: the localization is sensitive to many more phonon modes than the energy.

3.4 Electron Transfer Dynamics 18 ACS Paragon Plus Environment

Page 18 of 30

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

The canonically averaged photo-induced ET dynamics from the CH3NH3PbI3 sensitizer into the TiO2 substrate is demonstrated in Figure 6a. The ET coordinate starts at 0.5 at t=0, indicating that absorption of a photon at the CH3NH3PbI3/TiO2 interface has a 50% probability of creating a charge transfer excited state. The remaining 50% of the ET events occur by both adiabatic and nonadiabatic mechanisms. The ET coordinate is defined by integrating the PE electron density over the region of the simulation cell occupied by the electron acceptor, i.e., TiO2, eq. 6. The relative amounts of adiabatic and nonadiabatic ET are calculated by separating the overall evolution of the ET coordinate into the contributions due to the changes in the localization and occupation, respectively, eq. 7. The schematic of Figure 6b illustrates the three ET mechanisms. A photon can promote an electron directly from donor to acceptor. This mechanism requires strong donoracceptor coupling.

Adiabatic ET also operates in the strong donor-acceptor interaction

regime. Nonadiabatic ET does not require strong donor-acceptor interaction and, thereby, occurs in a broader range of systems. The strong interaction between CH3NH3PbI3 and the TiO2 surface favors both direct and adiabatic ET, while the high density of TiO2 acceptor states facilitates nonadiabatic ET. All ET mechanisms lead to ultrafast electron injection. The total ET line is fitted to a bi-exponential function. The faster ET component arises almost entirely from the adiabatic mechanism. It can be assigned to a rapid shift of the electron density from the donor to the acceptor through the regions of strong chemical interaction, most notably I-Ti chemical bonds, Figures 2 and S1. The slower ET component reflects transfer of the electron that is localized in CH3NH3PbI3, corresponding to the experimental sub-200fs time scale.39 The CH3NH3PbI3 layer used here is very thin compared to experiment, and therefore, one can expect that the slower ET component obtained in the calculation dominates the experimental data.

19 ACS Paragon Plus Environment

The Journal of Physical Chemistry

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 6. (a) Average ET dynamics. The solid black, dashed blue, and dotted red lines represent the contributions from total, adiabatic and nonadiabatic ET, respectively. The open circles give a bi-exponential fit with the time scales shown. 50% of ET at the interface occurs by direct electron photoexcitation into TiO2. Photo-excited electrons generated inside perovskite transfer to TiO2 by an approximately equal combination of adiabatic and nonadiabatic mechanisms. (b) Schematic of the photoinduced ET mechanism. Adiabatic ET occurs by passing over a transition state barrier (curved red arrow). Non-adiabatic ET occurs via a hop between donor and acceptor states (downward blue arrow). Photoexcitation can promote the electron directly from the donor material to a state that is localized on the acceptor (upward green arrow).

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

Figure 7. Examples of individual ET events. The lines have the same meaning as in Figure 6a. (a) Adiabatic ET with a small nonadiabatic component. (b) Nonadiabatic ET with a small adiabatic contribution. (c) ET with equal adiabatic and nonadiabatic contributions. The adiabatic component exhibits a sequence of forward and backward transfer events due to several re-crossings of the transition state, Figure 6b. (d) An initial delay, associated with search for a transition state, is followed by rapid adiabatic ET. Figure 6 represents the average behavior of the system during the ET process. A more detailed picture is presented in Figure 7, which shows examples of individual ET events. The examples demonstrate interesting interplay between adiabatic, nonadiabatic, and direct ET mechanisms, and exhibit a broad range of possibilities for the ET dynamics. Figure 7a displays a rapid adiabatic ET without a nonadiabatic contribution. Figure 7b shows an example dominated by the nonadiabatic ET, which starts after a short initial delay. The example of Figure 7c demonstrates ET exhibiting a complex interplay of adiabatic and nonadiabatic contributions. The adiabatic component exhibits a sequence of forward and backward transfer events due to several transition state re-crossings. The oscillation in the 21 ACS Paragon Plus Environment

The Journal of Physical Chemistry

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 30

total ET process is driven by atomic motions that produces the re-crossings of the adiabatic ET barrier. Finally, Figure 7d demonstrates an initial delay, associated with search for a transition state, followed by rapid adiabatic ET. On the short time scale, the ET is driven entropically, by the larger density of TiO2 acceptor states within the relevant energy window, compared to the than CH3NH3PbI3 states state density, Figure 4. In the long time limit, prior to the electron-hole recombination,17, 37 the electron relaxes to the TiO2 CB minimum, which is lower than the CH3NH3PbI3 CB minimum. The energy difference between the CH3NH3PbI3 and TiO2 CB minima constitutes the main driving force for the photo-induced charge separation. The individual ET events shown in Figure 7 represent examples of electron injection that can be detected by single-molecule spectroscopies. Despite the diversity in the ET dynamics,

the

average

ET

is

straightforward.

Photoexcitation

of

the

CH3NH3PbI3/TiO2 interface has a high probability of creating a charge separated state. Electrons generated inside CH3NH3PbI3 transfer into TiO2 by the adiabatic mechanism if they are located in regions of chemical-like CH3NH3PbI3-TiO2 interaction facilitated by I-Ti bonds, or by the nonadiabatic mechanism if such interaction is absent.

4. CONCLUSIONS We have used real-time TDDFT combined with NAMD to simulate photo-induced ET at the CH3NH3PbI3/TiO2 interface. Showing good agreement with the experimental time scale,39 the reported study mimics the observed time-resolved ET process, establishes the mechanisms responsible for the electron injection, characterizes the atomic motions that participate in the movement of charge in this system, and generates a detailed understanding of the ET dynamics. 22 ACS Paragon Plus Environment

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

The simulation shows that the ET occurs on an sub-100 fs time scale in the presence of a strong interaction between CH3NH3PbI3 and TiO2. Photon absorption near the interface has a high probability of creating a charge-separated state. Electrons approaching the interface from CH3NH3PbI3 bulk can transfer to TiO2 by both adiabatic and nonadiabatic mechanisms. Thermal fluctuations of atoms of the donor and acceptor materials create opportunities for transient chemical bonding across the interface, notably, formation of I-Ti bonds. The bonding is not seen at 0K. The strong donor-acceptor coupling generated through such bonds makes charge transfer excitations optically active. It is also responsible for the adiabatic ET mechanism. In the absence of the strong donor-acceptor coupling, the ET is nonadiabatic. Because the interfacial chemical bonding is not seen at 0K, one expects primarily nonadiabatic ET at lower temperatures. The nonadiabatic mechanism also produces fast ET due to high density of acceptor states, and therefore, only a weak temperature dependence in the ET time can be expected. Adiabatic ET is often viewed as a thermally activated process.87 Transient formation of I-Ti bonds induced by thermal nuclear fluctuations can be regarded as part of the activation. The situation differs from ET at molecule/TiO2 interfaces,

45-51

where donor-acceptor chemical bonds are much more stable

and exist both at zero and ambient temperatures. The electron couples primarily to lower frequency vibrational motions including both I-Pb stretching mode and librations of the organic cations. These motions generate an inhomogeneous distribution of initial conditions for the photoinduced ET, modulate the donor-acceptor separation, and create nonadiabatic coupling. The higher frequency Ti-O vibrations available at the CH3NH3PbI3/TiO2 interface have a smaller influence on the ET dynamics. They shift the electronic density between the donor and acceptor species and facilitate adiabatic ET. The complex interplay of the electron transfer mechanisms, donoracceptor and electron-phonon interactions, and phonon dynamics creates a variety of 23 ACS Paragon Plus Environment

The Journal of Physical Chemistry

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

individual electron injection scenarios. All ET mechanisms seen in the simulations lead to ultrafast ET, indicating that the broad range of the reported experimental ET times is associated with charge and exciton diffusion in perovskite, rather than the interfacial dynamics. The ultrafast interfacial injection guarantees efficient photo-induced charge separation that creates favorable conditions for operation of perovskite solar cells.

Supporting Information. Geometric structure and projected density of states for the CH3NH3PbI3/TiO2 system with a twice thicker perovskite layer. These materials are available free of charge via the Internet at http://pubs.acs.org.

Acknowledgements R.L. is grateful to the National Science Foundation of China, Grant no. 21573022, the Recruitment Program of Global Youth Experts, and the Beijing Normal University Startup package. W.H.F. thanks the National Science Foundation of China, Grant nos. 21520102005 and 21421003. O.V.P. acknowledges support of the U.S. National Science Foundation, Grant no. CHE-1565704.

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

References 1. O'Regan, B.; Gratzel, M., A Low-Cost, High-Efficiency Solar Cell Based on Dye-Sensitized Colloidal TiO2 Films. Nature 1991, 353, 737-740. 2. Kojima, A.; Teshima, K.; Shirai, Y.; Miyasaka, T., Organometal Halide Perovskites as VisibleLight Sensitizers for Photovoltaic Cells. J. Am. Chem. Soc. 2009, 131, 6050-6051. 3. Tavakoli, M. M.; Tavakoli, R.; Hasanzadeh, S.; Mirfasih, M. H., Interface Engineering of Perovskite Solar Cell Using a Reduced-Graphene Scaffold. J. Phys. Chem. C 2016, 120, 19531-19536. 4. Sun, C.; Guo, Y.; Fang, B.; Yang, J.; Qin, B.; Duan, H.; Chen, Y.; Li, H.; Liu, H., Enhanced Photovoltaic Performance of Perovskite Solar Cells Using Polymer P(VDF-TrFE) as a Processed Additive. J. Phys. Chem. C 2016, 120, 12980-12988. 5. Son, D.-Y.; Bae, K.-H.; Kim, H.-S.; Park, N.-G., Effects of Seed Layer on Growth of Zno Nanorod and Performance of Perovskite Solar Cell. J. Phys. Chem. C 2015, 119, 10321-10328. 6. Paek, S.; Cho, N.; Choi, H.; Jeong, H.; Lim, J. S.; Hwang, J.-Y.; Lee, J. K.; Ko, J., Improved External Quantum Efficiency from Solution-Processed (CH3NH3)PbI3 Perovskite/PC71BM Planar Heterojunction for High Efficiency Hybrid Solar Cells. J. Phys. Chem. C 2014, 118, 25899-25905. 7. Wu, R.; Yang, B.; Zhang, C.; Huang, Y.; Cui, Y.; Liu, P.; Zhou, C.; Hao, Y.; Gao, Y.; Yang, J., Prominent Efficiency Enhancement in Perovskite Solar Cells Employing Silica-Coated Gold Nanorods. J. Phys. Chem. C 2016, 120, 6996-7004. 8. 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-480. 9. Guo, Z.; Manser, J. S.; Wan, Y.; Kamat, P. V.; Huang, L., Spatial and Temporal Imaging of Long-Range Charge Transport in Perovskite Thin Films by Ultrafast Microscopy. Nat. Commun. 2015, 6, 1-8. 10. Pazos-Outón, L. M., et al., Photon Recycling in Lead Iodide Perovskite Solar Cells. Science 2016, 351, 1430-1433. 11. Lee, M. M.; Teuscher, J.; Miyasaka, T.; Murakami, T. N.; Snaith, H. J., Efficient Hybrid Solar Cells Based on Meso-Superstructured Organometal Halide Perovskites. Science 2012, 338, 643-647. 12. Nie, W. Y., et al., High-Efficiency Solution-Processed Perovskite Solar Cells with MillimeterScale Grains. Science 2015, 347, 522-525. 13. Liao, W. Q., et al., Fabrication of Efficient Low-Bandgap Perovskite Solar Cells by Combining Formamidinium Tin Iodide with Methylammonium Lead Iodide. J. Am. Chem. Soc. 2016, 138, 1236012363. 14. Manser, J. S.; Christians, J. A.; Kamat, P. V., Intriguing Optoelectronic Properties of Metal Halide Perovskites. Chem. Rev. 2016, 116, 12956-13008. 15. Liu, M.; Johnston, M. B.; Snaith, H. J., Efficient Planar Heterojunction Perovskite Solar Cells by Vapour Deposition. Nature 2013, 501, 395-398. 16. 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. 17. Yang, Y.; Yang, M. J.; Li, Z.; Crisp, R.; Zhu, K.; Beard, M. C., Comparison of Recombination Dynamics in CH3NH3PbBr3and CH3NH3PbI3 Perovskite Films: Influence of Exciton Binding Energy. J. Phys. Chem. Lett. 2015, 6, 4688-4692. 18. Wehrenfennig, C.; Eperon, G. E.; Johnston, M. B.; Snaith, H. J.; Herz, L. M., High Charge Carrier Mobilities and Lifetimes in Organolead Trihalide Perovskites. Adv. Mater. 2014, 26, 15841589. 19. Etgar, L.; Gao, P.; Xue, Z.; Peng, Q.; Chandiran, A. K.; Liu, B.; Nazeeruddin, M. K.; Grätzel, M., Mesoscopic CH3NH3PbI3/TiO2 Heterojunction Solar Cells. J. Am. Chem. Soc. 2012, 134, 17396-17399. 20. Manser, J. S.; Saidaminov, M. I.; Christians, J. A.; Bakr, O. M.; Kamat, P. V., Making and Breaking of Lead Halide Perovskites. Acc. Chem. Res. 2016, 49, 330-338.

25 ACS Paragon Plus Environment

The Journal of Physical Chemistry

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

21. Xie, F. X.; Su, H.; Mao, J.; Wong, K. S.; Choy, W. C. H., Evolution of Diffusion Length and Trap State Induced by Chloride in Perovskite Solar Cell. J. Phys. Chem. C 2016, 120, 21248-21253. 22. Tsai, H. H., et al., High-Efficiency Two-Dimensional Ruddlesden-Popper Perovskite Solar Cells. Nature 2016, 536, 312-316. 23. Song, D., et al., Managing Carrier Lifetime and Doping Property of Lead Halide Perovskite by Postannealing Processes for Highly Efficient Perovskite Solar Cells. J. Phys. Chem. C 2015, 119, 2281222819. 24. 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. 25. Stewart, R. J.; Grieco, C.; Larsen, A. V.; Maier, J. J.; Asbury, J. B., Approaching Bulk Carrier Dynamics in Organo-Halide Perovskite Nanocrystalline Films by Surface Passivation. J J. Phys. Chem. Lett. 2016, 7, 1148-1153. 26. Yang, Y.; Ostrowski, D. P.; France, R. M.; Zhu, K.; van de Lagemaat, J.; Luther, J. M.; Beard, M. C., Observation of a Hot-Phonon Bottleneck in Lead-Iodide Perovskites. Nat. Photonics 2016, 10, 5359. 27. Dong, Q.; Fang, Y.; Shao, Y.; Mulligan, P.; Qiu, J.; Cao, L.; Huang, J., Electron-Hole Diffusion Lengths > 175 µm in Solution-Grown CH3NH3PbI3 Single Crystals. Science 2015, 347, 967-970. 28. Zhou, H.; Chen, Q.; Li, G.; Luo, S.; Song, T.-b.; Duan, H.-S.; Hong, Z.; You, J.; Liu, Y.; Yang, Y., Interface Engineering of Highly Efficient Perovskite Solar Cells. Science 2014, 345, 542-546. 29. Stranks, S. D.; Eperon, G. E.; Grancini, G.; Menelaou, C.; Alcocer, M. J. P.; Leijtens, T.; Herz, L. M.; Petrozza, A.; Snaith, H. J., Electron-Hole Diffusion Lengths Exceeding 1 Micrometer in an Organometal Trihalide Perovskite Absorber. Science 2013, 342, 341-344. 30. Shi, D., et al., Low Trap-State Density and Long Carrier Diffusion in Organolead Trihalide Perovskite Single Crystals. Science 2015, 347, 519-522. 31. Mei, A., et al., A Hole-Conductor–Free, Fully Printable Mesoscopic Perovskite Solar Cell with High Stability. Science 2014, 345, 295-298. 32. Qin, P.; Kast, H.; Nazeeruddin, M. K.; Zakeeruddin, S. M.; Mishra, A.; Bauerle, P.; Gratzel, M., Low Band Gap S,N-Heteroacene-Based Oligothiophenes as Hole-Transporting and Light Absorbing Materials for Efficient Perovskite-Based Solar Cells. Energ. Environ. Sci. 2014, 7, 2981-2985. 33. Schulz, P.; Edri, E.; Kirmayer, S.; Hodes, G.; Cahen, D.; Kahn, A., Interface Energetics in Organo-Metal Halide Perovskite-Based Photovoltaic Cells. Energ. Environ. Sci. 2014, 7, 1377-1381. 34. de Quilettes, D. W.; Vorpahl, S. M.; Stranks, S. D.; Nagaoka, H.; Eperon, G. E.; Ziffer, M. E.; Snaith, H. J.; Ginger, D. S., Impact of Microstructure on Local Carrier Lifetime in Perovskite Solar Cells. Science 2015, 348, 683-686. 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 Simulation. J. Am. Chem. Soc. 2016, 138, 3884-3890. 36. 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. 37. 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. 38. 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. 39. 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. 40. Ponseca, C. S., et al., Organometal Halide Perovskite Solar Cell Materials Rationalized: Ultrafast Charge Generation, High and Microsecond-Long Balanced Mobilities, and Slow Recombination. J. Am. Chem. Soc. 2014, 136, 5189-5192. 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

41. Zhu, Z., et al., Efficiency Enhancement of Perovskite Solar Cells through Fast Electron Extraction: The Role of Graphene Quantum Dots. J. Am. Chem. Soc. 2014, 136, 3760-3763. 42. Wang, L.; McCleese, C.; Kovalsky, A.; Zhao, Y.; Burda, C., Femtosecond Time-Resolved Transient Absorption Spectroscopy of CH3NH3PbI3 Perovskite Films: Evidence for Passivation Effect of PbI2. J. Am. Chem. Soc. 2014, 136, 12205-12208. 43. Marchioro, A.; Teuscher, J.; Friedrich, D.; Kunst, M.; van de Krol, R.; Moehl, T.; Gratzel, M.; Moser, J.-E., Unravelling the Mechanism of Photoinduced Charge Transfer Processes in Lead Iodide Perovskite Solar Cells. Nat Photon 2014, 8, 250-255. 44. Xing, G.; Wu, B.; Chen, S.; Chua, J.; Yantara, N.; Mhaisalkar, S.; Mathews, N.; Sum, T. C., Interfacial Electron Transfer Barrier at Compact TiO2/ CH3NH3PbI3 Heterojunction. Small 2015, 11, 3606-3613. 45. Stier, W.; Prezhdo, O. V., Thermal Effects in the Ultrafast Photoinduced Electron Transfer from a Molecular Donor Anchored to a Semiconductor Acceptor. Isr. J. Chem.2002, 42, 213-224. 46. Stier, W.; Prezhdo, O. V., Nonadiabatic Molecular Dynamics Simulation of Light-Induced, Electron Transfer from an Anchored Molecular Electron Donor to a Semiconductor Acceptor. J. Phys. Chem. B 2002, 106, 8047-8054. 47. Rego, L. G. C.; Batista, V. S., Quantum Dynamics Simulations of Interfacial Electron Transfer in Sensitized TiO2 Semiconductors. J. Am. Chem. Soc. 2003, 125, 7989-7997. 48. Abuabara, S. G.; Rego, L. G. C.; Batista, V. S., Influence of Thermal Fluctuations on Interfacial Electron Transfer in Functionalized TiO2 Semiconductors. J. Am. Chem. Soc.2005, 127, 18234-18242. 49. Duncan, W. R.; Stier, W. M.; Prezhdo, O. V., Ab Initio Nonadiabatic Molecular Dynamics of the Ultrafast Electron Injection across the Alizarin− TiO2 Interface. J. Am. Chem. Soc. 2005, 127, 7941-7951. 50. Duncan, W. R.; Prezhdo, O. V., Temperature Independence of the Photoinduced Electron Injection in Dye-Sensitized TiO2 Rationalized by Ab Initio Time-Domain Density Functional Theory. J. Am. Chem. Soc. 2008, 130, 9756-9762. 51. Monti, A.; Negre, C. F. A.; Batista, V. S.; Rego, L. G. C.; de Groot, H. J. M.; Buda, F., Crucial Role of Nuclear Dynamics for Electron Injection in a Dye-Semiconductor Complex. J. Phys. Chem. Lett. 2015, 6, 2393-2398. 52. Fischer, S. A.; Duncan, W. R.; Prezhdo, O. V., Ab Initio Nonadiabatic Molecular Dynamics of Wet-Electrons on the TiO2 Surface. J. Am. Chem. Soc. 2009, 131, 15483-15491. 53. 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. 54. Tafen, D. N.; Long, R.; Prezhdo, O. V., Dimensionality of Nanoscale TiO2 Determines the Mechanism of Photoinduced Electron Injection from a Cdse Nanoparticle. Nano Lett. 2014, 14, 17901796. 55. Tafen, D.; Prezhdo, O. V., Size and Temperature Dependence of Electron Transfer between CdSe Quantum Dots and a TiO2 Nanobelt. J. Phys. Chem. C 2015, 119, 5639-5647. 56. Long, R.; English, N. J.; Prezhdo, O. V., Minimizing Electron–Hole Recombination on TiO2 Sensitized with Pbse Quantum Dots: Time-Domain Ab Initio Analysis. J. Phys. Chem. Lett. 2014, 5, 2941-2946. 57. 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. 58. Wang, L. J.; Long, R.; Prezhdo, O. V., Time-Domain Ab Initio Modeling of Photoinduced Dynamics at Nanoscale Interfaces. In Annu. Rev. Phys. Chem., 2015, 66, 549-579. 59. Craig, C. F.; Duncan, W. R.; Prezhdo, O. V., Trajectory Surface Hopping in the TimeDependent Kohn-Sham Approach for Electron-Nuclear Dynamics. Phys. Rev. Lett. 2005, 95, 163001. 60. 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. 27 ACS Paragon Plus Environment

The Journal of Physical Chemistry

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

61. Akimov, A. V.; Prezhdo, O. V., Large-Scale Computations in Chemistry: A Bird's Eye View of a Vibrant Field. Chem. Rev. 2015, 115, 5797-5890. 62. Pal, S.; Trivedi, D. J.; Akimov, A. V.; Aradi, B.; Frauenheim, T.; Prezhdo, O. V., Nonadiabatic Molecular Dynamics for Thousand Atom Systems: A Tight-Binding Approach toward Pyxaid. J. Chem. Theory Comput. 2016, 12, 1436-1448. 63. Wang, L. J.; Akimov, A.; Prezhdo, O. V., Recent Progress in Surface Hopping: 2011-2015. J. Phys. Chem. Lett. 2016, 7, 2100-2112. 64. Nebgen, B.; Prezhdo, O. V., Fragment Molecular Orbital Nonadiabatic Molecular Dynamics for Condensed Phase Systems. J. Phys. Chem. A 2016, 120, 7205-7212. 65. Tretiak, S.; Igumenshchev, K.; Chernyak, V., Exciton Sizes of Conducting Polymers Predicted by Time-Dependent Density Functional Theory. Phys. Rev. B 2005, 71, 033201. 66. Baer, R.; Neuhauser, D., Real-Time Linear Response for Time-Dependent Density-Functional Theory. J. Chem. Phys. s 2004, 121, 9803-9807. 67. Chapman, C. T.; Liang, W.; Li, X., Ultrafast Coherent Electron–Hole Separation Dynamics in a Fullerene Derivative. J. Phys. Chem. Lett. 2011, 2, 1189-1192. 68. Kohn, W.; Sham, L. J., Self-Consistent Equations Including Exchange and Correlation Effects. Phys. Rev. 1965, 140, A1133-A1138. 69. 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. 70. 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. 71. Jasper, A. W.; Nangia, S.; Zhu, C.; Truhlar, D. G., Non-Born−Oppenheimer Molecular Dynamics. Acc. Chem. Res. 2 2006, 39, 101-108. 72. Sousa, C.; Tosoni, S.; Illas, F., Theoretical Approaches to Excited-State-Related Phenomena in Oxide Surfaces. Chem. Rev. 2013, 113, 4456-4495. 73. Grimme, S., Semiempirical Gga-Type Density Functional Constructed with a Long-Range Dispersion Correction. J. Comput. Chem. 2006, 27, 1787-1799. 74. Kresse, G.; Furthmüller, J., Efficient Iterative Schemes for Ab Initio Total-Energy Calculations Using a Plane-Wave Basis Set. Phys. Rev. B 1996, 54, 11169-11186. 75. Perdew, J. P.; Burke, K.; Ernzerhof, M., Generalized Gradient Approximation Made Simple. Phys. Rev. Lett. 1996, 77, 3865-3868. 76. Blöchl, P. E., Projector Augmented-Wave Method. Phys. Rev. B 1994, 50, 17953-17979. 77. Umari, P.; Mosconi, E.; De Angelis, F., Relativistic GW Calculations on CH3NH3PbI3 and CH3NH3SnI3 Perovskites for Solar Cell Applications. Sci. Rep. 2014, 4, 4467. 78. Nie, W., et al., Light-Activated Photocurrent Degradation and Self-Healing in Perovskite Solar Cells. Nat. Commun. 2016, 7, 11574. 79. Umebayashi, T.; Asai, K.; Kondo, T.; Nakao, A., Electronic Structures of Lead Iodide Based Low-Dimensional Crystals. Phys. Rev. B 2003, 67, 155405. 80. See, A. K.; Thayer, M.; Bartynski, R. A., Angle-Resolved Inverse-Photoemission Study of the Nearly Perfect TiO2(110) Surface. Phys. Rev. B 1993, 47, 13722-13729. 81. Prezhdo, O. V.; Kindt, J. T.; Tully, J. C., Perturbed Ground State Method for Electron Transfer. J. Chem. Phys. 1999, 111, 7818-7827. 82. 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. 83. Nie, W., et al., High-Efficiency Solution-Processed Perovskite Solar Cells with MillimeterScale Grains. Science 2015, 347, 522-525. 84. Gong, J., et al., Electron-Rotor Interaction in Organic-Inorganic Lead Iodide Perovskites Discovered by Isotope Effects. J. Phys. Chem. Lett. 2016, 7, 2879-2887.

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

85. Giarola, M.; Sanson, A.; Monti, F.; Mariotto, G.; Bettinelli, M.; Speghini, A.; Salviulo, G., Vibrational Dynamics of Anatase TiO2: Polarized Raman Spectroscopy and Ab Initio Calculations. Phys. Rev. B 2010, 81, 174305. 86. Zhang, W. F. H., Y. L.; Zhang, M. S.; Yin, Z.; Chen, Q., Raman Scattering Study on Anatase TiO2 Nanocrystals. J. Phys. D: Appl. Phys.2000, 33, 5. 87. Stier, W.; Duncan, W. R.; Prezhdo, O. V., Thermally Assisted Sub-10 Fs Electron Transfer in Dye-Sensitized Nanocrystalline TiO2 Solar Cells. Adv. Mater. 2004, 16, 240-244.

29 ACS Paragon Plus Environment

The Journal of Physical Chemistry

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

TOC Graphic

30 ACS Paragon Plus Environment

Page 30 of 30