Rashba Band Splitting in Organohalide Lead Perovskites: Bulk and

May 3, 2017 - Francesco Ambrosio , Julia Wiktor , Filippo De Angelis , Alfredo Pasquarello. Energy & Environmental Science 2018 11 (1), 101-105 ...
1 downloads 0 Views 1MB Size
Subscriber access provided by HKU Libraries

Perspective

Rashba Band Splitting in Organohalide Lead Perovskites: Bulk and Surface Effects Edoardo Mosconi, Thibaud Etienne, and Filippo De Angelis J. Phys. Chem. Lett., Just Accepted Manuscript • Publication Date (Web): 03 May 2017 Downloaded from http://pubs.acs.org on May 4, 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 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 17

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

Rashba Band Splitting in Organohalide Lead Perovskites: Bulk and Surface Effects

Edoardo Mosconi,a,b Thibaud Etienne, a,c Filippo De Angelis a,b *

a

Computational Laboratory for Hybrid/Organic Photovoltaics (CLHYO), CNR-ISTM, via Elce di

Sotto, I-06123, Perugia, Italy. b

CompuNet, Istituto Italiano di Tecnologia, Via Morego 30, 16163 Genova, Italy.

c

Institut Charles Gerhardt - CNRS and Université de Montpellier, Place Eugène Bataillon - 34095

Montpellier, France.

E-mail: [email protected]

1 Environment ACS Paragon Plus

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

Abstract A Rashba/Dresselhaus band splitting has been recently measured in organohalide perovskites and invoked in various experiments as a possible cause for the reduced electron-hole recombination rates observed in this class of materials. In this perspective we discuss the interplay of electronic and nuclear degrees of freedom in defining such an effect in realistic methylammonium lead iodide (MAPbI3) models. We distinguish between bulk and surface effects and find that while a spatially local (in time and space) effect may be at work in the bulk, a “static” band-splitting effect is found at surfaces due to structural distortion. The proposed surface effect is consistent with the low surface recombination reported for MAPbI3 single crystals and might contribute to the success of organohalide perovskites.

2 Environment ACS Paragon Plus

Page 2 of 17

Page 3 of 17

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

Organohalide perovskites are dominating the landscape of emerging photovoltaic technologies.1-10 This class of materials exhibits a combination of physical properties which are hard to find even in well established inorganic semiconductors, including their optical11-14 and charge generation and transfer properties,9-10 the long carrier lifetimes15-17 and their sizable mobility.17-18 Having photovoltaic devices as a main applicative target, the most attractive property of organohalide lead perovskites is the low recombination rate of photogenerated carriers, which implies long diffusion length and allows efficient extraction of charge carriers even with modest mobility in polycrystalline films, thus high currents and voltages in optimized devices. On the one hand, it has been reported that methylammonium lead iodide (MAPbI3) and related compounds are generally defect-tolerant, thus scattering at defects or impurities does not severely limit the carrier lifetime and diffusion length, opposite to what found in typical semiconductors.

19-25,27-28

Keeping defect

tolerance in mind, it is then clear that the reasons behind the reduced carrier recombination rate in MAPbI326 are an inherent material property related to the combination of the constituent materials, i.e. lead, iodine (or different halide combinations) and the A-site cations. The presence of mobile dipolar organic cations has been related to the charge separation efficiency, by virtue of an exciton screening mechanism27 or by creating local screening domains,23, 28 although alloying of different cations, including the inorganic (spherical) caesium ion, has been shown to improve the solar cells efficiency.29 The presence of lead and to some extent of iodine imparts the material a giant spinorbit coupling (SOC),30 which dictates many materials properties, including band-gap variation in tin and lead perovskites,31 modulation of the band-gap with distortions from the cubic symmetry32 and the reduction of the carrier effective masses.33 SOC could be related to the observed spindependent charge recombination and spin-polarized carrier dynamics in MAPbI3,34-35 opening the way to spintronics based on organohalide perovskites.

3 Environment ACS Paragon Plus

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

A property intimately connected to SOC is the so-called Rashba/Dresselhaus effect, which has received considerable attention from a theoretical point of view.31-36 This effect16, 37-39 is the consequence of the breaking of inversion symmetry in the crystal in a direction orthogonal to a k-point sampling plane, and it is described by the so-called Bychkov-Rashba Hamiltonian:

where we find the Rashba primary correlation factor ߙ’, the electrostatic potential V, the linear momentum ‫( ∥ܘ‬defining the momentum space sampling, which is orthogonal to the gradient vector), the Pauli spin matrices σ and the gradient operator (સୄ ). Re-writing the Hamiltonian as:

the electric field projection can be included into the general Rashba interaction coefficient, hereafter ߙ.

The eigenvalues difference (ߝ ஷ ) at a particular k-point (∆k) sets the ߙ value:

±∆k points the respective position of the two ߝ± wells vertices for energy curves crossing at the k origin, and ߝ ஷ is defining the difference between the two energy curves at the vertex positions, see Scheme 1 for a definition of relevant parameters.

4 Environment ACS Paragon Plus

Page 4 of 17

Page 5 of 17

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

Scheme 1. Schematic representation of Rashba splitting for parabolic bands, along with relevant interaction parameters, see text for definitions. The valence and conduction bands (VB and CB) in the presence (absence) of Rashba splitting are the solid red and blue (dashed) lines, respectively. The shifted absorption and radiative recombination processes are also illustrated.

A single (parabolic) band will split into two bands by virtue of the Rashba effect, which are separated in k-space by ∆k and in energy by ε≠. The Rashba effect is intimately connected to the presence of heavy nuclei and to the presence of internal electric fields such as those due to the presence of polar (ferroelectric) structures. In MAPbI3 the valence (conduction) bands (VB and CB in Scheme 1) are mainly contributed by iodine (lead) and thus they undergo a different band 5 Environment ACS Paragon Plus

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

splitting, by virtue of the different SOC characterizing lead and iodine. This introduces a different splitting for the VB and CB, Scheme 1. The connection to recombination in a perovskite solar cell is then clear: as the photogenerated carriers thermalize at the band extrema which are located at different points in k-space, radiative recombination, which preserves k, will be quenched as in an indirect band-gap semiconductor. In MAPbI3 and in organohalide perovskites in general the VB and CB usually have different curvature, thus a different density of states, and spin helicity, which would further slow down recombination.16 While this band-splitting mechanism is an elegant framework that has gained considerable attention in the literature, it requires in principle a non-centrosymmetric crystal or a local-site symmetry breaking.40 Considering the I4/mcm space group of MAPbI3 it is likely that the material does not exhibit a macroscopic polarization at room temperature,41 though it could show an observable ferroelectric behaviour at lower temperature, i.e. just above the transition to the orthorhombic phase.42,43 This would in principle prevent the Rashba mechanism to be operative at room temperature, unless polar inclusions in polycrystalline samples are considered.44 Rappe and co-workers have however put forward the hypothesis that structural fluctuations, possibly accessible at room temperature, may give rise to local electric fields extending beyond the characteristic Rashba length scale (1-2 nm),16 thus not requiring the presence of macroscopic lack of inversion symmetry for this effect to be operative. We previously reported on the temporal and spatial scale of this “dynamical Rashba effect” in bulk MAPbI3,45 by means of Car-Parrinello Molecular Dynamics (CPMD),46-47 finding that even in a globally centrosymmetric structure, Figure 1a-b, the coupled inorganic/organic degrees of freedom evolution gives rise to a Rashba splitting which fluctuates on sub-ps time scale, Figure 1c-d. A similar effect was reported by Azarhoosh et al.48 The time average of this local Rashba effect is non zero on the investigated ∼2 ps time scale, longer than the expected time scale of electronic transitions, Figure 1c-d. This effect is however progressively quenched in globally centrosymmetric structures on increasing the probed

6 Environment ACS Paragon Plus

Page 6 of 17

Page 7 of 17

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

spatial scale from 1 to 4 formula units both in polar and apolar structures,45 compare the magnitude of the effect in Figure 1c and 1d. For the entire (initially apolar) 2x2x2 system a single band structure calculation carried out an a randomly selected structure from the dynamics delivers α values of ∼ 2 and ∼1 for the CB and VB, suggesting the persistence of the local Rashba effect up to this spatial scale.45 Notice that the data for the smaller sub-systems are purely indicative and are reasonably exaggerated since periodicity was enforced a posteriori for these sub-systems.

Figure 1. a-b) Structural models of MAPbI3viewed along the [001], a), and [110], b), directions. The native 2x2x2 supercell is enclosed by blue lines, while the reduced systems employed to sample different space scales are enclosed by green (4 formula units) and red (1 formula unit) lines. The notation 1-apolar/polar and 4-apolar/polar in c) and d) refers to globally polar and apolar

7 Environment ACS Paragon Plus

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

structures, respectively, made by 1 and 4 MAPbI3 units, in relation to the starting orientation of the methylammonium cations, whose dipoles are schematically represented by the arrows in panel b). The dashed lines in c) and d) represent the α time average, with values of 6.2-9.6 eV Å (c) and 2.72.4 eV Å (d), blue and red dashed lines, respectively. The φ dihedral angle discussed in Figure 2 is illustrated schematically.

From an experimental standpoint we are aware of a few studies which have directly or indirectly probed or invoked a Rashba effect to explain the observed phenomenology. Niesner et al. have directly observed a “giant Rashba splitting” in MAPbBr3 using angle-resolved photoelectron spectroscopy.49 They found a ring-shaped valence band and circular dichroism due to SOC in the low-temperature orthorhombic phase, which was maintained also in the high temperature cubic phase. A high α value of ∼10 eVÅ was reported for both phases, comparable to the highest values we found in our simulated dynamical Rashba effect for MAPbI3, see Figure 1c-d. Notably, these measurements are likely surface sensitive and probe a large fraction of the sample. Hutter et al. have reported time-resolved photo-conductance measurements showing that generation of free mobile charges is maximized for excitation energies just above what appears as an indirect band-gap, possibly related to the Rashba band splitting, Scheme 1.50 The employed microwave conductivity technique should probe photogenerated charge carriers in fairly small volumes, depending on the carrier mobility, though the signal is obtained from the entire sample. Wang et al. also suggested the presence of an indirect band-gap in MAPbI3 which they attributed to Rashba band splitting.51 Under hydrostatic pressure the Rashba splitting would be reduced due to a pressure-induced reduction in local electric field around the lead atom. Associated to this change, an increased carrier recombination and radiative efficiency was measured, consistent with the removal of the indirect band-gap. In this case, the optical measurements supposedly probe

8 Environment ACS Paragon Plus

Page 8 of 17

Page 9 of 17

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

the entire material volume, with some possible surface sensitivity due to the strong perovskite absorption coefficient limiting the light penetration depth. At the current stage of research, a question likely arises: Assuming a perfectly crystalline environment with no defects or other imperfections, are the experiments above probing an inhomogeneity in the sample due to local (time-dependent) distortions, or are they probing the sample surface, or a combination of the two? In other words, is the proposed Rashba effect an intrinsic bulk (but local) material property that can be measured or are the sample surfaces “interfering” with the possible observation of a bulk effect? Since the effect of surfaces on the Rashba effect in organohalide perovskites has not to our knowledge been reported, here we simulate two prototypical MAPbI3 (001) surfaces characterized by MAI- and PbI2-terminations, see Figure 2a and 2d, respectively, likely representing MAPbI3 surfaces when the material is grown in MAI-rich or PbI2-rich conditions.52 Our models are based on a 2x2x2 slab of the tetragonal unit cell with four layers of MA cations,53 with a starting apolar arrangement of the organic cations originated from the bulk structures of Figure 1. To simulate the investigated systems under thermal conditions we performed CPMD simulations on both slabs, setting the temperature at 330 ± 20 K.

9 Environment ACS Paragon Plus

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 2. a) - d). Optimized structures of the MAI-terminated and PbI2-terminated model slabs. b) e). Distribution of the angles formed between the methylammonium C-N axis and the perovskite ab plane during the CPMD simulation, see Figure 1 for the definition of the angle, for the bare MAIand PbI2-terminated MAPbI3 slabs. e) - f) SOC-DFT band structures calculated on the optimized structures of the MAI-terminated and PbI2-terminated slabs. The inset on the left side of panel c) represents the tetragonal Brillouin zone with labels of the high symmetry points. The inset on the right side of panel c) shows the calculated VB structure for the PbI2-terminated slab after cutting the capping PbI2 layers, i.e. for the system enclosed by the dashed lines in panel d) without further geometrical optimization.

10 Environment ACS Paragon Plus

Page 10 of 17

Page 11 of 17

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

An interesting parameter describing the impact of surfaces and of surface termination on the systems’ structure is the average arrangement of the organic cations, Figures 2b and 2e. As it can be noticed the organic cations in the outermost layers show a different orientation than those in the inner layers (bulk-like) in both MAI- and PbI2-terminated surfaces. This surface effect only marginally affects the inorganic sub-lattice in MAI-terminated slab, with almost perfectly preserved PbI6 octahedra, Figure 2a, while a pyramidalization of the outermost lead atoms is observed in the PbI2-terminated surface, Figure 2d. The band structure for the two slabs calculated by SOC-DFT on the optimized geometry obtained by the average dynamical structure reflects this geometrical difference, showing a larger Rashba splitting for the PbI2- than for the MAI-terminated slab, Figures 2f and 2c, respectively, despite the latter shows a stronger preferential orientation of the outermost MA cations. Notice the typical underestimate of the band-gap typical of SOC-DFT in both cases, with the PbI2-terminated slab showing a band-gap comparable to that of the bulk at the same level of theory.31 As in bulk MAPbI3, the largest Rashba splitting is observed along the Γ→X and Γ→M directions, see inset of Figure 2c for a representation of the tetragonal Brillouin zone. As expected, the bands are flat (and not split) along the Γ→Z direction, corresponding to the direction of surface truncation. A notable difference between the investigated slabs and the bulk is the comparable Rashba effect on the VB and CB in the former, while in the bulk the CB exhibits ca. a factor 2 larger splitting than the VB, due to the dominant lead contribution to the CB.45 This fairly large Rashba effect on the VB (∼2 eVÅ) calculated for the PbI2-terminated slab is smaller albeit of the same magnitude to that observed for MAPbBr349 and it clearly indicates that structural distortion due to surface termination may be responsible of a sizable fraction of the observed effect, especially when surface-sensitive techniques are employed. We cannot exclude that even higher effects could be observed for less ideal surfaces than the ones considered here, such as those predicted by Haruyama et al. characterized by partial loss of PbI2 units.54 We also investigated the effect of

11 Environment ACS Paragon Plus

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

removing the PbI2-capping layers from the PbI2-terminated slab, see right inset of Figure 2c, finding a splitting which is ca. a factor 2 larger than in the MAI-terminated case, having the same composition but different structure. This suggests that the PbI2-surface termination introduces structural distortions which are propagated to the inner layers. In summary, we have shown how surfaces will likely play a role in defining a bandsplitting effect in organohalide perovskites, possibly related to a Rashba/Dresselhaus effect. The bulk effect is likely a local (in space and time) effect which could play a role in diminishing electron-hole recombination and thus boost the solar cell efficiency. We believe, however, that this effect will be difficult to probe with conventional techniques, while it could be explored by microwave-conductivity experiments, which could locally probe the photogenerated charge carriers. On the other hand, the reported surface band-splitting is consistent with the low surface recombination reported MAPbI3 single crystals,55 ∼2–3 orders of magnitude lower than that in conventional (unpassivated) semiconductors employed in solar cells,55 and might be a key contribution to the success of organohalide perovskites.

12 Environment ACS Paragon Plus

Page 12 of 17

Page 13 of 17

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

Methods. For all calculations the experimental cell parameters reported by Poglitsch and Weber56 have been used. Geometry optimization and CPMD simulations have been carried out with the Quantum Espresso package57 along with the GGA-PBE58 functional. Electron-ion interactions were described by scalar relativistic ultrasoft pseudopotentials with electrons from N and C 2s, 2p; H 1s; I 5s, 5p; Pb 6s, 6p, 5d shells explicitly included in the calculations. Plane-wave basis set cutoffs for the smooth part of the wave functions and the augmented density were 25 and 200 Ry, respectively. CPMD simulations of the bulk have been performed with an integration time step of 5 au, for a total simulation time of ca. 12 ps after a few ps of thermal equilibration. The fictitious mass used for the electronic degrees of freedom is 500 au, and we set the atomic masses to their real values except for hydrogen atoms, for which the deuterium mass is used. CPMD simulations on perovskite slabs were performed with an integration time step of 10 au, fictitious mass for the electronic degrees of freedom 1000 au and atomic masses set to an identical values of 5 amu to enhance the dynamical sampling. Initial ion position randomization has been used to reach a temperature of ∼350 K, without further applying any thermostat. The electronic structure analysis has been carried out on the optimized geometries using the PWscf code with Plane-wave basis set cutoffs for the smooth part of the wave functions and the augmented density were 25 and 200 Ry. Spin-orbit coupling has been included using ultrasoft PBE-GGA pseudo potentials. 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) Im, J.-H.; Lee, C.-R.; Lee, J.-W.; Park, S.-W.; Park, N.-G. 6.5% efficient perovskite quantum-dot-sensitized solar cell. Nanoscale 2011, 3, 4088-4093. (3) 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. (4) Burschka, J.; Pellet, N.; Moon, S.-J.; Humphry-Baker, R.; Gao, P.; Nazeeruddin, M. K.; Gratzel, M. Sequential deposition as a route to high-performance perovskitesensitized solar cells. Nature 2013, 499, 316-319. (5) Heo, J. H.; Im, S. H.; Noh, J. H.; Mandal, T. N.; Lim, C.-S.; Chang, J. A.; Lee, Y. H.; Kim, H.-j.; Sarkar, A.; NazeeruddinMd, K., et al. Efficient inorganic-organic

13 Environment ACS Paragon Plus

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

hybrid heterojunction solar cells containing perovskite compound and polymeric hole conductors. Nat Photon 2013, 7, 486-491. (6) Liu, M.; Johnston, M. B.; Snaith, H. J. Efficient planar heterojunction perovskite solar cells by vapour deposition. Nature 2013, 501, 395-398. (7) Kim, H.-S.; Lee, C.-R.; Im, H.-H.; Lee, K.-B.; Moehl, T.; Marchioro, A.; Moon, S.-J.; Humpry-Baker, R.; Yum, J.-H.; Moser, J. E., et al. lead iodide perovskite sensitized all-solid-state submicron thin film mesoscopic solar cell with efficiency exceeding 9%. Sci. Rep. 2012, 2, 591. (8) 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. (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) 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. (11) Chiarella, F.; Zappettini, A.; Licci, F.; Borriello, I.; Cantele, G.; Ninno, D.; Cassinese, A.; Vaglio, R. Combined experimental and theoretical investigation of optical, structural, and electronic properties of CH3NH3SnX3 thin films (X=Cl, Br). Physical Review B 2008, 77, 045129. (12) Ogomi, Y.; Morita, A.; Tsukamoto, S.; Saitho, T.; Fujikawa, N.; Shen, Q.; Toyoda, T.; Yoshino, K.; Pandey, S. S.; Ma, T., et al. CH3NH3SnxPb(1–x)I3 Perovskite Solar Cells Covering up to 1060 nm. J. Phys. Chem. Lett. 2014, 5, 10041011. (13) Eperon, G. E.; Stranks, S. D.; Menelaou, C.; Johnston, M. B.; Herz, L. M.; Snaith, H. J. Formamidinium lead trihalide: a broadly tunable perovskite for efficient planar heterojunction solar cells. Energy Environ. Sci. 2014, 7, 982-988. (14) Stoumpos, C. C.; Malliakas, C. D.; Kanatzidis, M. G. Semiconducting Tin and Lead Iodide Perovskites with Organic Cations: Phase Transitions, High Mobilities, and Near-Infrared Photoluminescent Properties. Inorg. Chem. 2013, 52, 9019-9038. (15) Wehrenfennig, C.; Liu, M.; Snaith, H. J.; Johnston, M. B.; Herz, L. M. Chargecarrier dynamics in vapour-deposited films of the organolead halide perovskite CH3NH3PbI3-xClx. Energy & Environmental Science 2014, 7, 2269-2275. (16) Zheng, F.; Tan, L. Z.; Liu, S.; Rappe, A. M. Rashba Spin-Orbit Coupling Enhanced Carrier Lifetime in CH3NH3PbI3. Nano letters 2015, 15, 7794-7800. (17) Ponseca, C. S.; Savenije, T. J.; Abdellah, M.; Zheng, K.; Yartsev, A.; Pascher, T.; Harlang, T.; Chabera, P.; Pullerits, T.; Stepanov, A., 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. (18) Edri, E.; Kirmayer, S.; Mukhopadhyay, S.; Gartsman, K.; Hodes, G.; Cahen, D. Elucidating the charge carrier separation and working mechanism of CH3NH3PbI3−xClx perovskite solar cells. Nat. Commun. 2014, 5, 3461. (19) Kim, J.; Lee, S.-H.; Lee, J. H.; Hong, K.-H. The Role of Intrinsic Defects in Methylammonium Lead Iodide Perovskite. The Journal of Physical Chemistry Letters 2014, 5, 1312-1317. (20) Du, M. H. Efficient carrier transport in halide perovskites: theoretical perspectives. J. Mater. Chem. A 2014, 2, 9091-9098.

14 Environment ACS Paragon Plus

Page 14 of 17

Page 15 of 17

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) Yamada, Y.; Endo, M.; Wakamiya, A.; Kanemitsu, Y. Spontaneous Defect Annihilation in CH3NH3PbI3 Thin Films at Room Temperature Revealed by TimeResolved Photoluminescence Spectroscopy. J. Phys. Chem. Lett. 2015, 6, 482-486. (22) Buin, A.; Pietsch, P.; Xu, J.; Voznyy, O.; Ip, A. H.; Comin, R.; Sargent, E. H. Materials Processing Routes to Trap-Free Halide Perovskites. Nano letters 2014, 14, 6281-6286. (23) Ma, J.; Wang, L.-W. Nanoscale Charge Localization Induced by Random Orientations of Organic Molecules in Hybrid Perovskite CH3NH3PbI3. Nano letters 2015, 15, 248-253. (24) Frost, J. M.; Butler, K. T.; Brivio, F.; Hendon, C. H.; van Schilfgaarde, M.; Walsh, A. Atomistic Origins of High-Performance in Hybrid Halide Perovskite Solar Cells. Nano Lett. 2014, 14, 2584-2590. (25) Liu, S.; Zheng, F.; Koocher, N. Z.; Takenaka, H.; Wang, F.; Rappe, A. M. Ferroelectric Domain Wall Induced Band Gap Reduction and Charge Separation in Organometal Halide Perovskites. The Journal of Physical Chemistry Letters 2015, 6, 693-699. (26) 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, 1584-1589. (27) Even, J.; Pedesseau, L.; Katan, C. Analysis of Multivalley and Multibandgap Absorption and Enhancement of Free Carriers Related to Exciton Screening in Hybrid Perovskites. J. Phys. Chem. C 2014, 118, 11566-11572. (28) Quarti, C.; Mosconi, E.; De Angelis, F. Structural and electronic properties of organo-halide hybrid perovskites from ab initio molecular dynamics. Phys. Chem. Chem. Phys. 2015, 17, 9394-9409. (29) Saliba, M.; Matsui, T.; Seo, J.-Y.; Domanski, K.; Correa-Baena, J.-P.; Nazeeruddin, M. K.; Zakeeruddin, S. M.; Tress, W.; Abate, A.; Hagfeldt, A., et al. Cesium-containing triple cation perovskite solar cells: improved stability, reproducibility and high efficiency. Energy Environ. Sci. 2016, 9, 1989-1997. (30) 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. (31) Umari, P.; Mosconi, E.; De Angelis, F. Relativistic GW calculations on CH3NH3PbI3 and CH3NH3SnI3 Perovskites for Solar Cell Applications. Sci. Rep. 2014, 4, 4467. (32) Amat, A.; Mosconi, E.; Ronca, E.; Quarti, C.; Umari, P.; Nazeeruddin, M. K.; Grätzel, M.; De Angelis, F. Cation-Induced Band-Gap Tuning in Organohalide Perovskites: Interplay of Spin–Orbit Coupling and Octahedra Tilting. Nano Lett. 2014, 14, 3608-3616. (33) Giorgi, G.; Fujisawa, J.-I.; Segawa, H.; Yamashita, K. Small Photocarrier Effective Masses Featuring Ambipolar Transport in Methylammonium Lead Iodide Perovskite: A Density Functional Analysis. J. Phys. Chem. Lett. 2013, 4, 4213-4216. (34) Giovanni, D.; Ma, H.; Chua, J.; Grätzel, M.; Ramesh, R.; Mhaisalkar, S.; Mathews, N.; Sum, T. C. Highly Spin-Polarized Carrier Dynamics and Ultralarge Photoinduced Magnetization in CH3NH3PbI3 Perovskite Thin Films. Nano Lett. 2015, 15, 1553-1558. (35) Hsiao, Y.-C.; Wu, T.; Li, M.; Hu, B. Magneto-Optical Studies on SpinDependent Charge Recombination and Dissociation in Perovskite Solar Cells. Adv. Mater. 2015, 27, 2899-2906.

15 Environment ACS Paragon Plus

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

(36) Motta, C.; El-Mellouhi, F.; Kais, S.; Tabet, N.; Alharbi, F.; Sanvito, S. Revealing the role of organic cations in hybrid halide perovskite CH3NH3PbI3. Nat Commun 2015, 6. (37) Rashba, E. I. Properties of semiconductors with an extremum loop. 1. Cyclotron and combinational resonance in a magnetic field perpendicular to the plane of the loop. Soviet Physics-Solid State 1960, 2, 1109-1122. (38) Bychkov, Y. A.; Rashba, E. I. Oscillatory effects and the magnetic susceptibility of carriers in inversion layers. Journal of physics C: Solid state physics 1984, 17, 6039. (39) Kepenekian, M.; Robles, R.; Katan, C.; Sapori, D.; Pedesseau, L.; Even, J. Rashba and Dresselhaus Effects in Hybrid Organic–Inorganic Perovskites: From Basics to Devices. ACS nano 2015. (40) Zhang, X.; Liu, Q.; Luo, J.-W.; Freeman, A. J.; Zunger, A. Hidden spin polarization in inversion-symmetric bulk crystals. Nat Phys 2014, 10, 387-393. (41) Beilsten-Edmands, J.; Eperon, G. E.; Johnson, R. D.; Snaith, H. J.; Radaelli, P. G. Non-ferroelectric nature of the conductance hysteresis in CH3NH3PbI3 perovskite-based photovoltaic devices. Appl. Phys. Lett. 2015, 106, 173502. (42) Quarti, C.; Mosconi, E.; De Angelis, F. Interplay of Orientational Order and Electronic Structure in Methylammonium Lead Iodide: Implications for Solar Cell Operation. Chem. Mater. 2014, 26, 6557–6569. (43) Stroppa, A.; Quarti, C.; De Angelis, F.; Picozzi, S. Ferroelectric Polarization of CH3NH3PbI3: A Detailed Study Based on Density Functional Theory and Symmetry Mode Analysis. J. Phys. Chem. Lett. 2015, 6, 2223-2231. (44) Osherov, A.; Hutter, E. M.; Galkowski, K.; Brenes, R.; Maude, D. K.; Nicholas, R. J.; Plochocka, P.; Bulović, V.; Savenije, T. J.; Stranks, S. D. The Impact of Phase Retention on the Structural and Optoelectronic Properties of Metal Halide Perovskites. Adv. Mater. 2016, 28, 10757-10763. (45) Etienne, T.; Mosconi, E.; De Angelis, F. Dynamical Origin of the Rashba Effect in Organohalide Lead Perovskites: A Key to Suppressed Carrier Recombination in Perovskite Solar Cells? J. Phys. Chem. Lett. 2016, 7, 1638–1645. (46) Car, R.; Parrinello, M. Unified Approach for Molecular Dynamics and DensityFunctional Theory. Phys. Rev. Lett. 1985, 55, 2471-2474 (47) Giannozzi, P.; Angelis, F. D.; Car, R. First-principle molecular dynamics with ultrasoft pseudopotentials: Parallel implementation and application to extended bioinorganic systems. J. Chem. Phys. 2004, 120, 5903-5915. (48) Azarhoosh, P.; McKechnie, S.; Frost, J. M.; Walsh, A.; Schilfgaarde, M. v. Research Update: Relativistic origin of slow electron-hole recombination in hybrid halide perovskite solar cells. APL Mater. 2016, 4, 091501. (49) Niesner, D.; Wilhelm, M.; Levchuk, I.; Osvet, A.; Shrestha, S.; Batentschuk, M.; Brabec, C.; Fauster, T. Giant Rashba Splitting in CH3NH3PbBr3 Organic-Inorganic Perovskite. Physical Review Letters 2016, 117, 126401. (50) Hutter, E. M.; Gelvez-Rueda, M. C.; Osherov, A.; Bulovic, V.; Grozema, F. C.; Stranks, S. D.; Savenije, T. J. Direct-indirect character of the bandgap in methylammonium lead iodide perovskite. Nat Mater 2017, 16, 115-120. (51) Wang, T.; Daiber, B.; Frost, J. M.; Mann, S. A.; Garnett, E. C.; Walsh, A.; Ehrler, B. Indirect to direct bandgap transition in methylammonium lead halide perovskite. Energy Environ. Sci. 2017, 10, 509-515. (52) Mosconi, E.; Grancini, G.; Roldán-Carmona, C.; Gratia, P.; Zimmermann, I.; Nazeeruddin, M. K.; De Angelis, F. Enhanced TiO2/MAPbI3 Electronic Coupling by Interface Modification with PbI2. Chem. Mater. 2016, 28, 3612-3615.

16 Environment ACS Paragon Plus

Page 16 of 17

Page 17 of 17

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

(53) Mosconi, E.; Azpiroz, J. M.; De Angelis, F. Ab Initio Molecular Dynamics Simulations of Methylammonium Lead Iodide Perovskite Degradation by Water. Chem. Mater. 2015, 27, 4885-4892. (54) 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. (55) Yang, Y.; Yan, Y.; Yang, M.; Choi, S.; Zhu, K.; Luther, J. M.; Beard, M. C. Low surface recombination velocity in solution-grown CH3NH3PbBr3 perovskite single crystal. Nat Commun 2015, 6, 7961. (56) Poglitsch, A.; Weber, D. Dynamic disorder in methylammoniumtrihalogenoplumbates (II) observed by millimeter-wave spectroscopy. J. Chem. Phys. 1987, 87, 6373-6378. (57) Giannozzi, P.; Baroni, S.; Bonini, N.; Calandra, M.; Car, R.; Cavazzoni, C.; Ceresoli, D.; Chiarotti, G. L.; Cococcioni, M.; Dabo, I., et al. QUANTUM ESPRESSO: A Modular and Open-Source Software Project for Quantum Simulations of Materials. J. Phys.: Condens. Matter 2009, 21, 395502. (58) Perdew, J. P.; Burke, K.; Ernzerhof, M. Generalized Gradient Approximation Made Simple. Phys. Rev. Lett. 1996, 77, 3865-3868.

17 Environment ACS Paragon Plus