Charge Carrier Trapping at Surface Defects of Perovskite Solar Cell

Jan 27, 2017 - A facile method to evaluate the influence of trap densities on perovskite solar cell performance. Bingbing Chen , Hongwei Hu , Teddy Sa...
1 downloads 0 Views 624KB Size
Subscriber access provided by UB + Fachbibliothek Chemie | (FU-Bibliothekssystem)

Letter

Charge Carrier Trapping at Surface Defects of Perovskite Solar Cell Absorbers: A First-Principles Study Hiroki Uratani, and Koichi Yamashita J. Phys. Chem. Lett., Just Accepted Manuscript • DOI: 10.1021/acs.jpclett.7b00055 • Publication Date (Web): 27 Jan 2017 Downloaded from http://pubs.acs.org on January 28, 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 18

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

Charge Carrier Trapping at Surface Defects of Perovskite Solar Cell Absorbers: A First-Principles Study Hiroki Uratani*,†,‡ and Koichi Yamashita*,†,‡

†Department of Chemical System Engineering, School of Engineering, The University of Tokyo, 7-3-1, Hongo, Bunkyo-ku, 113-8656 Tokyo, Japan ‡CREST-JST, 7 Gobancho, Chiyoda-ku, 102-0076 Tokyo, Japan

AUTHOR INFORMATION Corresponding Author *[email protected] *[email protected] Notes The authors declare no competing financial interests.

ACS Paragon Plus Environment

1

The Journal of Physical Chemistry Letters

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

Page 2 of 18

ABSTRACT

The trapping of charge carriers at defects on surfaces or grain boundaries is detrimental for the performance of perovskite solar cells (PSCs). For example, it is the main limiting factor for carrier lifetime. Moreover, it causes hysteresis in the current-voltage curves, which is considered to be a serious issue for PSCs' operation. In this work, types of surface defects responsible for carrier trapping are clarified by a comprehensive first-principles investigation into surface defects of tetragonal CH3NH3PbI3 (MAPbI3). Considering defect formation energetics, it is proposed that a Pb-rich condition is preferred to an I-rich one, however, a moderate condition might possibly be the best choice. Our result paves the way for improving the performance of PSCs through rational strategy of suppressing carrier trapping at surface defects.

TOC GRAPHICS

ACS Paragon Plus Environment

2

Page 3 of 18

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

Organic-inorganic hybrid perovskites (OIHPs) such as CH3NH3PbI3 (MAPbI3) are attracting much interest for application as perovskite solar cells (PSCs). Since the report of the first PSC in 2009,1 their photoconversion efficiencies (PCEs) have increased from 3.8% to 22.1%2 to date. One of the detrimental factors in PSCs' performance is the trapping of charge carriers at defects on their surfaces or grain boundaries. First, the carrier trapping at surface defects causes nonradiative recombination loss and deteriorates the carrier lifetime. For example, Chen et al.3 passivated grain boundaries of MAPbI3 with PbI2, resulting in the improvement in the carrier lifetime and enhancement of PCE from 0.66% to 12.00%, because of the suppression of carrier trapping at defects on the grain boundaries. Noel et al.4 have prolonged carrier lifetime by almost seven times to achieve an improvement in PCE from 13.1% to 16.5% via surface passivation with Lewis bases, which donate electrons to fill midgap states that stem from surface defects, reducing carrier traps. Moreover, the carrier trapping on surfaces is one of main origins of the hysteresis in the current-voltage curves,5–7 which is a serious drawback of PSCs. Shao et al.6 eliminated the hysteresis by reducing electron traps using surface passivation by fullerene. Im et al.8 and Nie et al.9 also overcame the hysteresis through another approach: they fabricated grains with large size, which have reduced interfacial area associated with grain boundaries, i.e. reduced carrier trapping on surfaces. In conjunction with these successes, spectroscopic studies revealed that carrier trapping states actually exist on surfaces.10,11 However, to date, the types of surface defects that act as carrier traps are still not identified. Although theoretical works by Haruyama et al. revealed the absence of deep midgap states on surfaces without point defects,12,13 little is known about surfaces that have them. Defect physics of OIHPs has been substantially studied by first-principles,14–18 however, previous studies have focused on bulk point defects, not on surface defects.

ACS Paragon Plus Environment

3

The Journal of Physical Chemistry Letters

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

Page 4 of 18

The aim of this work is to clarify the types of surface defects responsible for carrier trapping, and to understand the relationship between surface defects formation and chemical conditions, in order to propose design rules for high-performance PSCs. We investigate the energy levels of defect states and the formation energy of various types of surface defects comprehensively, by first-principles calculations based on density-functional theory (DFT) for MAPbI3 systems containing surface defects, using slab models. We discuss the possible types of carrier-trapping surface defects for different conditions: I-rich, Pb-rich, and moderate. Finally, we show that the Pb-rich condition is more suitable than the I-rich condition and propose that the moderate condition might be the most appropriate. We focused on (001) surface, which is determined as one of the major types of MAPbI3 facet by first-principles calculations12,13 and powder X-ray diffraction measurements.19 2 × 2 surface models containing 384-552 atoms were constructed from the tetragonal MAPbI3 crystal structure. Three types of surface termination, that have been typically studied,12,13,20–22 were considered: "MAI", "flat", and "vacant" termination (Figure 1).

Figure 1. "MAI", "flat", and "vacant" type of termination. H, Pb, C, I, and N atoms are shown as pink, gray, brown, purple, and blue, respectively. The simulation cell is shown by dotted lines.

ACS Paragon Plus Environment

4

Page 5 of 18

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 outermost layers of "MAI" and "flat" termination are composed of MAI and PbI2, respectively. "Vacant" termination is constructed by eliminating eight PbI2 units from the "flat" termination. The cell parameters were fixed to the twice of experimental value (a = b = 17.602 Å)23 during calculations. A vacuum region of 15 Å or more was added to each model. Following the conventional approach,12,13,24 the directions of MA cations were alternated so that slabs had almost no net electric dipole moment. For each calculation of defect-containing systems, one side of the slab contained a defect. The types of surface defects considered are shown in Table 1. Here, VX (where X is an atom or a group of atoms) means that X is removed from the outermost layer, Xi means that excessive X is placed on the surface, and XY (where Y is an atom or a group of atoms) means that Y site of the outermost layer is substituted by X.

ACS Paragon Plus Environment

5

The Journal of Physical Chemistry Letters

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

Page 6 of 18

Table 1. Considered types of surface defects and their formation energy (in eV) under different conditions. termination

defect

I-rich

moderate

Pb-rich

MAI

VI

2.38

1.77

1.15

MAI

VMA

0.32

0.91

1.54

MAI

Pbi

4.69

3.45

2.24

MAI

PbMA

2.35

1.70

1.13

MAI

PbI

5.97

4.14

2.31

flat

Ii

-0.03

0.57

1.19

flat

VI

4.55

3.94

3.32

flat

Pbi

2.94

1.71

0.50

flat

VPb

-0.13

1.10

2.31

vacant

Ii

-0.12

0.48

1.10

vacant

VI

1.76

1.15

0.54

vacant

Pbi

2.71

1.48

0.26

vacant

VPb

-1.40

-0.17

1.04

The calculations were performed with Vienna Ab-initio Simulation Package25–28 using a planewave basis set with the projector augmented-wave (PAW) method.29 For H, Pb, C, I, and N atoms, 1, 14, 4, 7, and 5 valence electrons were explicitly treated, respectively. Structural relaxations were performed and then defect formation energy calculations were conducted at the relaxed structure, with 500 eV of a plane-wave cutoff. In structural relaxations, positions of all atoms were relaxed until the Hellman-Feynman forces became below 0.02 eV⁄Å. The Brillouin zone (BZ) was sampled at the Γ-point in the structural relaxations, and the 2 × 2 × 1 Γ-centered k-point mesh was adopted in the defect formation energy calculations. For structural relaxations and defect formation energy calculations, Perdew-Burke-Ernzerhof formalism of generalized

ACS Paragon Plus Environment

6

Page 7 of 18

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

gradient approximation (GGA-PBE)30 was adopted without the inclusion of spin-orbit coupling (SOC), as it is widely used for DFT studies of OIHPs including MAPbI3 and has provided reliable results for structure and thermodynamic properties.12–15,17,21,22,24,31–33 However, for defect level calculations, GGA-PBE without SOC fails because of valence band maxima (VBM)/conduction band minima (CBM) levels overestimation as a result of self-interaction error inherent in GGA and the neglect of SOC.16,18 For this reason, to calculate defect levels, we adopted Heyd-Scuseria-Ernzerhof (HSE) hybrid functional34 with 43% of Hartree-Fock exchange (denoted as  = 0.43) and the inclusion of SOC following the approach of Du,16,18 which provides reasonable positioning of VBM, CBM and defect levels for MAPbI3. In defect level calculations, we used the structure relaxed in the above-mentioned manner. Because the computational cost of the hybrid functional is much larger than that of GGA-PBE, the calculations were conducted using 400 eV of a plane-wave cutoff and BZ was sampled at the Γpoint. Note that VBM and CBM of MAPbI3 are located at the Γ-point.35 The energy levels of defect states were determined as the corresponding Kohn-Sham orbital levels. It should be noted that we did not determine the defect levels by the difference in total energy between the charged and the neutral system, which is the usual strategy for bulk defects,14,16–18 because, for slab models, the total energy of a charged system is meaningless because of the electrostatic repulsion between slabs. For visualization, VESTA software36 was used. The energy levels of defect states compared with the VBM and CBM of each system are shown in Figure 2.

ACS Paragon Plus Environment

7

The Journal of Physical Chemistry Letters

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

Page 8 of 18

Figure 2. Defect levels compared with VBM and CBM calculated by HSE ( = 0.43) with SOC. Energy levels are referenced to the VBM of each system. Defect levels are shown by red horizontal lines. VBM and CBM are shown by black horizontal lines. Note that defect states are described if and only if they are located between the VBM and CBM. These energy levels are referenced to the VBM of each system. The Kohn-Sham states corresponding to the VBM, CBM, and defect states were determined by their absolute energy levels, occupation numbers, and orbital shapes (for details, see the Supporting Information). The slight fluctuations in the VBM-CBM gap depending on the types of defects are presumably attributed to the finite-size error and the doping effect of defects. The defect formation energy was calculated taking into account the chemical potential of I, Pb, and MA denoted as  ,  , and  , respectively, following the approach introduced by Yin et al.14 (for details, see the

ACS Paragon Plus Environment

8

Page 9 of 18

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

Supporting Information). Assuming equilibrium among MAPbI3, MAI, and PbI2, the defect formation energy under the I-rich condition ( = 0 eV,  = −2.44 eV, and  = −3.15 eV), moderate condition ( = −0.61 eV,  = −1.21 eV, and  = −2.57 eV), and Pb-rich condition ( = −1.22 eV,  = 0 eV, and  = −1.93 eV) is shown in Table 1. On the "MAI" surface, only VMA has low formation energy, however, it has no deep midgap defect states (Figure 2). This means that no or very few carrier traps can form on the "MAI" surface regardless of the conditions, therefore "MAI" is the most suitable surface termination from the viewpoint of avoiding carrier trapping on surfaces. Remarkably, the Br counterpart of the “MAI” termination, i.e. “MABr” termination, was experimentally observed by Ohmann et al.37 on mechanically cleaved MAPbBr3 crystals. However, it should be noted that Lindblad et al. reported an understoichiometry in I and/or N at the annealed MAPbI3/TiO2 interfaces,38 suggesting that "MAI" is not the dominant surface termination in annealed samples. This observation can be supported by the fact that PbI2 species remains in grain boundaries after thermal annealing,3 because MAI has higher volatility than PbI2. From Table 1, it is suggested that different types of defects are dominant depending on the condition (I-rich, moderate, and Pbrich). Under the I-rich condition, MA vacancies (VMA) on “MAI” surface, excessive I atoms (Ii) on “flat” and “vacant” surfaces, and Pb vacancies (VPb) on “flat” and “vacant” surfaces are dominant. Among these, only Ii has a deep carrier-trapping state (Figure 2). Hence, only Ii can be the origin of carrier trapping on surfaces under the I-rich condition. Note that Ii have remarkably low formation energy (-0.03 eV on “flat” surface and -0.12 eV on “vacant” surface) under the I-rich condition (Table 1). The charge density distribution derived from the KohnSham orbitals corresponding to the defect states of Ii on "flat" and "vacant" surfaces are shown in

ACS Paragon Plus Environment

9

The Journal of Physical Chemistry Letters

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

Page 10 of 18

Figure 3a and b, respectively, revealing that these defect states are mainly composed of p orbitals of excessive I atoms, and Pb s orbital also contributes in the case of "vacant" surface.

Figure 3. Charge density distribution corresponding to defect states (left) and atomic structure of defects (right) of (a): Ii on "flat" surface, (b): Ii on "vacant" surface, (c): Pbi on "flat" surface, (d): Pbi on "vacant" surface, and (e): VI on "vacant" surface. Atoms in the outermost two layers are drawn. The charge density is illustrated by yellow isosurfaces. Blue isosurfaces indicate negative charge density in the core region, which is a common artifact of the PAW method. On the other hand, under the Pb-rich condition, excessive Pb atoms (Pbi) are dominant for "flat" and "vacant" surfaces, and I atom vacancies (VI) are also dominant for "vacant" surfaces. From

ACS Paragon Plus Environment

10

Page 11 of 18

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

The Journal of Physical Chemistry Letters

Figure 2, it is revealed that Pbi generates deep carrier trapping states near the middle point between VBM and CBM on both "flat" and "vacant" surfaces. Also, on "vacant" surfaces, VI creates carrier trapping states. However, it is important to point out that the formation energy of Pbi (on “flat” and “vacant” surfaces) and VI (on “vacant” surface) under the Pb-rich condition (+0.26 eV or more) is higher than that of Ii under the I-rich condition (-0.03 eV or less). The charge density distribution of the defect states of Pbi on "flat" and "vacant" surfaces, and VI on "vacant" surfaces are shown in Figure 3c, d, and e, respectively. For Pbi, charge localization in the bonding region between the excessive Pb and a Pb atom in the outermost layer can be observed, indicating that the carrier trapping states come from Pb-Pb bonding orbitals. Also, the defect state of VI on "vacant" surface is found to be mainly composed of p orbital of the undercoordinated Pb atom. From the above discussion, I-rich and Pb-rich conditions entail carriertrapping surface defects that have different origins: excessive I atoms under I-rich conditions, and excessive Pb atoms and I atom vacancies under Pb-rich conditions. Since the formation energy of carrier-trapping surface defects under the Pb-rich condition is higher than that under the I-rich condition, the Pb-rich condition is better than the I-rich condition. This can account for the fact that higher photovoltaic performance is realized by the two-step method39 than the onestep method1,40 because of the use of high PbI2 concentration (namely, excess Pb relative to I, for growth of MAPbI3) in the two-step method; in the two-step method, PbI2 is first introduced into mesoporous TiO2 by spin-coating and then the PbI2/TiO2 composite film is dipped into MAI solution, whereas in the one-step method, stoichiometric mixed solution of MAI and PbI2 is dropped onto TiO2. Under the moderate condition, only VPb on "vacant" surfaces, which does not act as carrier traps, has low formation energy, suggesting that no or very few carrier-trapping surface defects can form under the moderate condition. Hence, from the aspect of minimizing

ACS Paragon Plus Environment

11

The Journal of Physical Chemistry Letters

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

Page 12 of 18

carrier trapping at surface defects, the moderate condition might be the best choice. Notably, Hsu et al. reported that moderate concentration of excess PbI2 maximizes carrier lifetime in MAPbI3,41 suggesting that the moderate condition actually minimizes the charge carrier trapping at surface defects. In summary, we studied various types of surface defects on three types of terminations: "MAI", "flat", and "vacant", by first-principles calculations. Combining the calculated defect levels and the defect formation energy, our results can be summarized in three points as follows. (i) Under the I-rich condition, excessive I atoms on “flat” and “vacant” surfaces are responsible for the carrier trapping. On the other hand, under the Pb-rich condition, I atom vacancies on "vacant" surface and excessive Pb atoms on both "flat" and "vacant" surfaces act as carrier traps. (ii) The formation of carrier-trapping surface defects under the Pb-rich condition is thermodynamically more unfavorable than under the I-rich condition. (iii) Under the moderate condition, any surface defects that act as carrier traps have high formation energy, i.e. cannot easily form. From the above, to reduce carrier trapping on surfaces or grain boundaries so as to improve carrier lifetime and avoid hysteresis, the Pb-rich condition is better than the I-rich condition, consistent with the superiority of the two-step method39 over the one-step method1,40 in terms of photovoltaic performance. Also, the moderate condition can be potentially the most desirable. ACKNOWLEDGMENT Computational resources were provided by Reserch Center for Computational Science, National Institutes of Natural Sciences, Okazaki, Japan. ASSOCIATED CONTENT

ACS Paragon Plus Environment

12

Page 13 of 18

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

Supporting Information Available: Details of the determination of relevant Kohn-Sham states to obtain Figure 2, and details of the defect formation energy calculations (PDF)

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)

Zhou, Y.; Zhu, K. Perovskite Solar Cells Shine in the “Valley of the Sun.” ACS Energy

Lett. 2016, 1, 64–67. (3)

Chen, Q.; Zhou, H.; Song, T. Bin; Luo, S.; Hong, Z.; Duan, H. S.; Dou, L.; Liu, Y.;

Yang, Y. Controllable Self-Induced Passivation of Hybrid Lead Iodide Perovskites toward High Performance Solar Cells. Nano Lett. 2014, 14, 4158–4163. (4)

Noel, N. K.; Abate, A.; Stranks, S. D.; Parrott, E. S.; Burlakov, V. M.; Goriely, A.;

Snaith, H. J. Enhanced Photoluminescence and Solar Cell Performance via Lewis Base Passivation of Organic–Inorganic Lead Halide Perovskites. ACS Nano 2014, 8, 9815–9821. (5)

Snaith, H. J.; Abate, A.; Ball, J. M.; Eperon, G. E.; Leijtens, T.; Noel, N. K.; Stranks, S.

D.; Wang, J. T. W.; Wojciechowski, K.; Zhang, W. Anomalous Hysteresis in Perovskite Solar Cells. J. Phys. Chem. Lett. 2014, 5, 1511–1515. (6)

Shao, Y.; Xiao, Z.; Bi, C.; Yuan, Y.; Huang, J. Origin and Elimination of Photocurrent

Hysteresis by Fullerene Passivation in CH3NH3PbI3 Planar Heterojunction Solar Cells. Nat. Commun. 2014, 5, 1–7.

ACS Paragon Plus Environment

13

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

(7)

Page 14 of 18

Chen, B.; Yang, M.; Priya, S.; Zhu, K. Origin of J-V Hysteresis in Perovskite Solar Cells.

J. Phys. Chem. Lett. 2016, 7, 905–917. (8)

Im, J.-H.; Jang, I.-H.; Pellet, N.; Grätzel, M.; Park, N.-G. Growth of CH3NH3PbI3

Cuboids with Controlled Size for High-Efficiency Perovskite Solar Cells. Nat. Nanotechnol. 2014, 9, 927–932. (9)

Nie, W.; Tsai, H.; Asadpour, R.; Blancon, J.-C.; Neukirch, A. J.; Gupta, G.; Crochet, J. J.;

Chhowalla, M.; Tretiak, S.; Alam, M. A.; et al. High-Efficiency Solution-Processed Perovskite Solar Cells with Millimeter-Scale Grains. Science 2015, 347, 522–525. (10) Wu, X.; Trinh, M. T.; Niesner, D.; Zhu, H.; Norman, Z.; Owen, J. S.; Yaffe, O.; Kudisch, B. J.; Zhu, X. Y. Trap States in Lead Iodide Perovskites. J. Am. Chem. Soc. 2015, 137, 2089– 2096. (11) Kong, W.; Ding, T.; Bi, G.; Wu, H. Optical Characterizations of the Surface States in Hybrid Lead–halide Perovskites. Phys. Chem. Chem. Phys. 2016, 18, 12626–12632. (12) 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. (13) Haruyama, J.; Sodeyama, K.; Han, L.; Tateyama, Y. Surface Properties of CH3NH3PbI3 for Perovskite Solar Cells. Acc. Chem. Res. 2016, 49, 554–561. (14) Yin, W.-J.; Shi, T.; Yan, Y. Unusual Defect Physics in CH3NH3PbI3 Perovskite Solar Cell Absorber. Appl. Phys. Lett. 2014, 104, 063903/1-063903/4.

ACS Paragon Plus Environment

14

Page 15 of 18

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

(15) Kim, J.; Lee, S. H.; Lee, J. H.; Hong, K. H. The Role of Intrinsic Defects in Methylammonium Lead Iodide Perovskite. J. Phys. Chem. Lett. 2014, 5, 1312–1317. (16) Du, M. H. Efficient Carrier Transport in Halide Perovskites: Theoretical Perspectives. J. Mater. Chem. A 2014, 2, 9091. (17) Agiorgousis, M. L.; Sun, Y.; Zeng, H.; Zhang, S. Strong Covalency-Induced Recombination Centers in Perovskite Solar Cell Material CH3NH3PbI3. J. Am. Chem. Soc. 2014, 136, 14570–14575. (18) Du, M.-H. Density Functional Calculations of Native Defects in CH3NH3PbI3: Effects of Spin–Orbit Coupling and Self-Interaction Error. J. Phys. Chem. Lett. 2015, 6, 1461–1466. (19) Baikie, T.; Fang, Y.; Kadro, J. M.; Schreyer, M.; Wei, F.; Mhaisalkar, S. G.; Graetzel, M.; White, T. J. Synthesis and Crystal Chemistry of the Hybrid Perovskite (CH3NH3)PbI3 for Solid-State Sensitised Solar Cell Applications. J. Mater. Chem. A 2013, 1, 5628. (20) Torres, A.; Rego, L. G. C. Surface Effects and Adsorption of Methoxy Anchors on Hybrid Lead Iodide Perovskites: Insights for Spiro-MeOTAD Attachment. J. Phys. Chem. C 2014, 118, 26947–26954. (21) 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. (22) Zhang, L.; Liu, X.; Su, J.; Li, J. First Principles Study of Molecular Adsorption on Lead Iodide Perovskite Surface: A Case Study of Halogen Bond Passivation for Solar Cell Application. J. Phys. Chem. C 2016, 120, 23536−23541.

ACS Paragon Plus Environment

15

The Journal of Physical Chemistry Letters

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

Page 16 of 18

(23) Kawamura, Y.; Mashiyama, H.; Hasebe, K. Structural Study on Cubic–Tetragonal Transition of CH3NH3PbI3. J. Phys. Soc. Jpn. 2002, 71, 1694–1697. (24) Mosconi, E.; Amat, A.; Nazeeruddin, M. K.; Grätzel, M.; De Angelis, F. First Principles Modeling of Mixed Halde Organometal Perovskites for Photovoltaic Applications. J. Phys. Chem. C 2013, 117, 13902−13913. (25) Kresse, G.; Hafner, J. Ab Initio Molecular Dynamics for Liquid Metals. Phys. Rev. B 1993, 47, 558–561. (26) Kresse, G.; Hafner, J. Ab Initio Molecular-Dynamics Simulation of the Liquid-Metal– amorphous-Semiconductor Transition in Germanium. Phys. Rev. B 1994, 49, 14251–14269. (27) Kresse, G.; Furthmüller, J. Efficiency of Ab-Initio Total Energy Calculations for Metals and Semiconductors Using a Plane-Wave Basis Set. Comput. Mater. Sci. 1996, 6, 15–50. (28) 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. (29) Blöchl, P. E. Projector Augmented-Wave Method. Phys. Rev. B 1994, 50, 17953–17979. (30) Perdew, J. P.; Burke, K.; Ernzerhof, M. Generalized Gradient Approximation Made Simple. Phys. Rev. Lett. 1996, 77, 3865–3868. (31) Azpiroz, J. M.; Mosconi, E.; Bisquert, J.; De Angelis, F. Defect Migration in Methylammonium Lead Iodide and Its Role in Perovskite Solar Cell Operation. Energy Environ. Sci. 2015, 8, 2118–2127.

ACS Paragon Plus Environment

16

Page 17 of 18

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

The Journal of Physical Chemistry Letters

(32) Giorgi, G.; Yamashita, K. Zero-Dimensional Hybrid Organic–Inorganic Halide Perovskite Modeling: Insights from First Principles. J. Phys. Chem. Lett. 2016, 7, 888−899. (33) Mosconi, E.; Ronca, E.; De Angelis, F. First Principles Investigation of the TiO2/Organohalide Perovskites Interface: The Role of Interfacial Chlorine. J. Phys. Chem. Lett. 2014, 5, 2619–2625. (34) Heyd, J.; Scuseria, G. E.; Ernzerhof, M. Hybrid Functionals Based on a Screened Coulomb Potential. J. Chem. Phys. 2003, 118, 8207–8215. (35) Umari, P.; Mosconi, E.; De Angelis, F. Relativistic GW Calculations on CH3NH3PbI3 and CH3NH3SnI3 Perovskites for Solar Cell Applications. Sci. Rep. 2014, 4, 4467. (36) Momma, K.; Izumi, F. VESTA 3 for Three-Dimensional Visualization of Crystal, Volumetric and Morphology Data. J. Appl. Crystallogr. 2011, 44, 1272–1276. (37) Ohmann, R.; Ono, L. K.; Kim, H. S.; Lin, H.; Lee, M. V.; Li, Y.; Park, N. G.; Qi, Y. Real-Space Imaging of the Atomic Structure of Organic-Inorganic Perovskite. J. Am. Chem. Soc. 2015, 137, 16049–16054. (38) Lindblad, R.; Bi, D.; Park, B.; Oscarsson, J.; Gorgoi, M.; Siegbahn, H.; Odelius, M.; Johansson, E. M. J.; Rensmo, H. Electronic Structure of TiO2/CH3NH3PbI3 Perovskite Solar Cell Interfaces. J. Phys. Chem. Lett. 2014, 5, 648–653. (39) Burschka, J.; Pellet, N.; Moon, S.-J.; Humphry-Baker, R.; Gao, P.; Nazeeruddin, M. K.; Grätzel, M. Sequential Deposition as a Route to High-Performance Perovskite-Sensitized Solar Cells. Nature 2013, 499, 316–320.

ACS Paragon Plus Environment

17

The Journal of Physical Chemistry Letters

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

Page 18 of 18

(40) 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. (41) Hsu, H. Y.; Ji, L.; Du, M.; Zhao, J.; Yu, E. T.; Bard, A. J. Optimization of PbI2/MAPbI3 Perovskite Composites by Scanning Electrochemical Microscopy. J. Phys. Chem. C 2016, 120, 19890−19895.

ACS Paragon Plus Environment

18