The Raman Spectrum of the CH3NH3PbI3 Hybrid Perovskite: Interplay

The broad, unstructured 200–400 cm–1 features are assigned to the torsional mode of the methylammonium cations, which we propose as a marker of th...
24 downloads 4 Views 3MB Size
Letter pubs.acs.org/JPCL

The Raman Spectrum of the CH3NH3PbI3 Hybrid Perovskite: Interplay of Theory and Experiment Claudio Quarti,*,† Giulia Grancini,‡ Edoardo Mosconi,† Paola Bruno,‡ James M. Ball,§ Michael M. Lee,§ Henry J. Snaith,§ Annamaria Petrozza,*,‡ and Filippo De Angelis*,† †

Computational Laboratory for Hybrid/Organic Photovoltaics (CLHYO), CNR-ISTM, I-06123, Perugia, Italy Center for Nano Science and Technology@PoliMi, Istituto Italiano di Tecnologia, Via Giovanni Pascoli 70/3, 20133, Milano, Italy § Department of Physics, University of Oxford, Clarendon Laboratory, Parks Road, Oxford, OX1 3PU, United Kingdom ‡

S Supporting Information *

ABSTRACT: We report the low-frequency resonant Raman spectrum of methylammonium lead-iodide, a prototypical perovskite for solar cells applications, on mesoporous Al2O3. The measured spectrum assignment is assisted by DFT simulations of the Raman spectra of suitable periodic and model systems. The bands at 62 and 94 cm−1 are assigned respectively to the bending and to the stretching of the Pb−I bonds, and are thus diagnostic modes of the inorganic cage. We also assign the librations of the organic cations at 119 and 154 cm−1. The broad, unstructured 200−400 cm−1 features are assigned to the torsional mode of the methylammonium cations, which we propose as a marker of the orientational disorder of the material. Our study provides the basis to interpret the Raman spectra of organohalide perovskites, which may allow one to further understand the properties of this important class of materials in relation to their full exploitation in solar cells. SECTION: Energy Conversion and Storage; Energy and Charge Transport

H

related organohalide MAPbCl3 perovskite shows various bands between 40 and 321 cm−1, associated to the librations of the MA organic cations and a band at 483 cm−1, assigned to the MA torsional mode.27 The assignment of the MA torsional mode represents a very important issue in the characterization of hybrid organic−inorganic systems, with reported frequencies in the 281−487 cm−1 range.27,33−35 Such large spread of MA torsional frequencies is quite surprising. To the best of our knowledge, no experimental data are available for such mode in solution or in the gas phase, but a few theoretical calculations have set it at 267−324 cm−1.36,37 Hybrid perovskites are intrinsically complex materials, where the presence of various types of interactions and structural disorder15 may play an important role in the materials properties. In this work, the assignment of the Raman spectrum of MAPbI3 is supported by first-principles DFT calculations38,39 carried out on extended periodic systems and, for interpretative purposes, on simplified models. While this type of analyses have been applied to perovskites, especially in relation to their ferroelectric properties,40−47 to the best of our knowledge this is the first report on the Raman spectrum of the prototypical MAPbI3 perovskite, along with the first DFT simulation of the vibrational spectra of hybrid perovskites.

ybrid lead−halide perovskites are revolutionizing the photovoltaic landscape, having been claimed as “the next big thing in photovoltaics”.1 From their first application in 2009 by Kojima et al. as solar cells sensitizers,2 photovoltaic devices based on these materials showed a fast and continuous increase in their performance,2−14 with recent top efficiency exceeding 15%. Methylammonium lead-iodide, hereafter MAPbI3, represents a prototype system for photovoltaic applications.5,6,11,13 In spite of the wide interest in MAPbI3, the characterization of its hybrid structure and the understanding of its chemical-physical properties are still in an early stage. If, on the one hand, the study of the structure and of the properties of this material have been carried out by X-ray diffraction and UV−vis spectroscopy,15−18 along with electronic structure calculations,18−21 a comprehensive vibrational spectroscopy study of MAPbI3 is still missing. Here we report a combined experimental and theoretical Raman vibrational analysis of MAPbI3 in the low-frequency region, where the bands associated to the vibrations of the interacting inorganic/organic constituents are expected. To the best of our knowledge, the only vibrational investigation on MAPbI3 was carried out by IR spectroscopy in the 700−3500 cm−1 region,22 where only the internal modes of the MA organic cations are present. Several vibrational spectroscopy studies of inorganic23−26 and hybrid (layered) perovskites27−32 have been reported, which may provide useful guidelines for the MAPbI3 spectral assignment. The Raman spectrum of the © 2013 American Chemical Society

Received: November 28, 2013 Accepted: December 26, 2013 Published: December 26, 2013 279

dx.doi.org/10.1021/jz402589q | J. Phys. Chem. Lett. 2014, 5, 279−284

The Journal of Physical Chemistry Letters

Letter

Figure 1. Left: Experimental resonant (a) and nonresonant calculated (tet-1, tet-2, ortho) Raman spectra of MAPbI3 (b−d). The tet-1, tet-2, and ortho optimized structures are shown on the right panel parallel to the ab plane. The periodic structures and Raman spectra were calculated by the Quantum Espresso program package.38

features, thus allowing us to establish some qualitatively interesting structure/property correlations. The experimental band falling at 62 cm−1, as well as that at 154 cm−1, find a good correspondence in the three theoretical spectra. The band at 94 cm−1 is well accounted in tet-2 and in ortho; in tet-1, many vibrational normal modes are found in the 80−100 cm−1 region, although the convoluted spectrum does not show a well-defined maximum. The band at 119 cm−1 instead is quite puzzling. This band does not find any parallel in tet-2 nor in ortho, while it can be possibly associated with the band at 141 cm−1 found in tet-1. Since the three models differ mainly in the orientation of the organic cations, the fact that the region around the experimental band at 119 cm−1 is differently described in the three simulated structures suggests that this is probably associated with the MA libration modes. The assignment of the Raman spectrum is reported in detail on the basis of the results obtained for the tet-1 structure. The assignments made for the tet-1 structure also hold for tet-2 and ortho, unless differently noted. The experimental band falling at 62 cm−1 can be assigned to the theoretical band predicted around 53 cm−1. The present band is mainly due to four normal modes falling at 50, 51, 56, and 60 cm−1, and the visual inspection of the eigenvectors associated to such modes (Figure 2a) demonstrates that they consist mainly in the bending of the I−Pb−I bonds and in the consequent libration of the cations due to the deformation of the inorganic cage. The assignment of the band at 62 cm−1 to the inorganic structure is confirmed by previous studies on related perovskites.28−31 Thus, the band measured for MAPbI3 at 62 cm−1 is a clear marker of the inorganic component. The slight underestimation of this vibrational frequency, which we calculated at 53 cm−1, is possibly due to the neglect of spin−orbit interactions,19,20 which are expected to lead to frequency increases.48,49 Inspection of the normal modes calculated in the region

The Raman spectrum measured under resonant conditions (excitation wavelength 532 nm) for a MAPbI3 film deposited on mesoporous Al2O3 is reported in Figure 1a (see Supporting Information for experimental and sample preparation details). The spectrum is quite similar to the spectra of related hybrid perovskites,28,30 and is characterized by a broad Raman signal in the 50−170 cm−1 region, with two well-defined maxima at 62 and 119 cm−1 and two less intense maxima at 94 and 154 cm−1. A broad and not resolved signal between 200 and 340 cm−1 is also found, plus an additional weak feature at 390 cm−1. In Figure 1b−d the experimental resonant Raman spectrum is compared to the nonresonant Raman spectra calculated in the harmonic approximation for two tetragonal structures (tet-1 and tet-2)19 and for one orthorhombic structure (ortho) composed by four MAPbI3 units (see models in the Supporting Information), representative of the room temperature and lowtemperature dominant phases.15,16 The two tetragonal structures differ for the initial orientation of the organic cations: tet-2 shows an ordered head-to-tail arrangement, while tet-1 shows a more disordered structure,19 inducing a different structural distortion, see Figure 1b−d. The lattice parameters for the investigated structures calculated at the LDA level used for the calculation of Raman intensities (see computational details in the Supporting Information) are in good agreement (within 1−2%) with experimental data15 (see Supporting Information). The tet-1 structure is 0.1 eV more stable than tet2, in agreement with our previous results.19 We stress here that while calculated vibrational frequencies can be directly compared to experimental values, the calculated Raman intensities are obtained in a nonresonant regime, thus in principle a different intensity distribution compared to resonant experimental data could be expected. Nevertheless, the calculated spectra reproduce quite well the experimental 280

dx.doi.org/10.1021/jz402589q | J. Phys. Chem. Lett. 2014, 5, 279−284

The Journal of Physical Chemistry Letters

Letter

normal mode, neither Raman active or inactive, in the range 200−280 cm−1, though our results could be sensitive to both DFT functional and anharmonic effects, as we will explain below. Furthermore, since our calculations are performed in the nonresonant regime, we cannot exclude that second-order (two phonons) effects may contribute to the observed intensity.24−26 The torsional vibrational frequencies of the organic cations, Figure 2d, are the only modes predicted in the 200−400 cm−1 range. In particular, LDA places these modes at 295−315, 288− 301, and 391−399 cm−1 in tet-1, tet-2, and ortho, respectively. Notably, computing the structure and vibrational frequencies by PBE, a down-shift of the MA torsional modes to 271−309 cm−1 is predicted for tet-1. The large torsional frequencies difference predicted for the two tetragonal structures compared to those calculated for the orthorhombic structure may seem surprising at first glance, but it is actually not totally unexpected considering the large distribution of torsional frequencies reported in the literature for MA-based compounds.27,33−35 Furthermore, comparing the theoretical spectra of Figure 1, we observe that the Raman intensity predicted for the torsional modes increases with the structural disorder, in the order tet-1 > tet-2 > ortho. It is thus reasonable to wonder whether the experimental bands at 200−390 cm−1 can be associated to the torsional mode of the MA cations. To check this possibility and to explain the large differences in torsional frequencies found in the different considered structures, we have carried out further analyses on simplified models obtained by “decomposition” of the optimized crystal structures: (i) the single isolated MA cations; (ii) the four MA cations; and (iii) the four MA cations, plus the three nearest iodide anions, all at the crystal structure geometries (see Figure 3 and Table 1). The calculated vibrational frequency of the torsional mode of the isolated MA cation at its optimized structure is predicted at 309 cm−1 by LDA (278 cm−1 when optimized by the same method used for periodic simulations in a large supercell) in good agreement with previous results.36,37 This mode is predicted to be Raman inactive, consistent with the C3v symmetry of the isolated MA cation. The torsional frequency down-shifts to 302 cm−1 when calculated by PBE, paralleling the frequency shift observed for the perovskite models. Interestingly, anharmonic effects50 lead to further frequency down-shifts for the MA torsional mode, delivering a calculated frequency of 277 (298) cm−1 by PBE (LDA). Based on this data, it is not unreasonable to speculate that a combination of DFT functional and anharmonic effects could down-shift the calculated torsional modes for the tetragonal structures to ∼250 cm−1, as in the experimental Raman spectrum. The torsional mode frequency calculated in the harmonic approximation by LDA for the isolated MA cations at their crystal structures geometries are strongly blue-shifted, by more than 150 cm−1. In this case, the Raman intensities slightly differ from zero because of the partial loss of C3v symmetry. The observed shift is caused by the distortion of the molecule from the gas phase geometry, as suggested by a detailed structural analysis, see Table S1 and Figure S4, Supporting Information. From Table 1, it is interesting to notice that when the interactions between the MA cations and the three closest iodides are considered, the frequencies of the torsional modes shift to lower values. This frequency down-shift amounts to ∼100 cm−1 for tet-1 and tet-2, and to ∼60 cm−1 for the ortho structure and is originated by hydrogen bonding between the ammonium hydrogen and the iodide atoms. Alongside, the

Figure 2. Normal modes computed for the tet-1 structure associated to the most intense Raman bands of MAPbI3. (a) I−Pb−I bending mode at 60 cm−1 (Raman intensity ∼791 Å4/amu); (b) Pb−I stretching mode at 88 cm−1 (∼21.82 Å4/amu); (c) MA libration mode at 141 cm−1 (Raman intensity 789 Å4/amu); (d) MA torsional mode at 295 cm−1 (∼2042 Å4/amu).

around 94 cm−1 reveals a major component from librations of the organic cations and from Pb−I stretchings, Figure 2b. In particular, normal modes mainly corresponding to Pb−I stretching are predicted at 72, 73, 87, and 88 cm−1 and they are characterized by small or medium Raman intensity (15%, or less, of the intensity of the band at 60 cm−1). Thus, the 94 cm−1 band is reasonably associated with both the Pb−I stretching and to libration modes of the cations, and it mainly brings information on the inorganic component of the material. This is consistent with the assignment of bands in the 110−121 cm−1 range to the Pb−Cl stretching mode25,28 and with the assignment of a 100 cm−1 band of a putative KPbI3 perovskite to Pb−I stretching.26 As pointed out previously, the fact that the Raman spectra calculated on the different models are quite different in the 100−150 cm−1 region suggests that the band at 119 cm−1 is probably associated to the motion of the organic cations. The visual inspection of the eigenvectors associated to the modes calculated for tet-1 (Figure 2c) confirms this point, showing that the modes predicted above 90 cm−1 consist mainly in the librations of the organic cations. The band falling at 119 cm−1 can be tentatively associated to the band calculated at 141 cm−1. A similar result holds for the weak band measured at 154 cm−1, that can be safely associated with the band calculated at 156 cm−1, corresponding to the libration of the organic cations.27 The assignment of the broad and unresolved band ranging from 200 to 350 cm−1 and the additional weak feature at ∼390 cm−1 represent a very subtle point in the interpretation of the Raman spectrum of MAPbI3. In fact, the theoretical spectra calculated for the three investigated structures do not show any 281

dx.doi.org/10.1021/jz402589q | J. Phys. Chem. Lett. 2014, 5, 279−284

The Journal of Physical Chemistry Letters

Letter

Raman intensity increases, going from ∼10 −3 to ∼20 Å4/amu for tet-1 and tet-2 and to ∼1 Å4/amu for the ortho structure. Notably, we find different hydrogen bonding patterns for the investigated structures (see Figure 3 and Table S2, Supporting Information). Tet-1 and tet-2 show two short (2.49−2.67 Å) and one long (2.82 Å) hydrogen bonds; the ortho structure shows instead three almost equivalent hydrogen bonds (∼2.5 Å), which preserve a C3v-like environment, leading to a decrease of the Raman intensity of the torsional modes. Our analysis shows that the frequency of the MA torsional modes in the crystal structures may be the result of two opposite effects, associated to the interactions between the cation and the inorganic cage. On the one hand, the deformation of the MA molecule shifts the torsional mode toward higher frequencies; on the other, the formation of specific hydrogen bonding interactions, typical of these compounds20,51 shifts this mode toward lower frequencies. Because of the known orientational disorder of the organic cations in MAPbI3,15,21 the different MA molecules are subject to different balances between molecular deformation and hydrogen bonding interactions and thus their torsional frequency is spread over a large range of values. It is also reasonable to expect that the torsional mode could be strongly influenced by all those parameters that affect the orientational order of the cations, like temperature and morphology. In particular, we expect that increasing the ordering of the cations, by, e.g., decreasing the temperature, a shift of the torsional mode to higher frequencies and an associated decrease in Raman intensity could be observed. In summary, we have reported the resonant Raman spectrum of the widely used MAPbI3 perovskite measured between 60 and 450 cm−1, assigning its main vibrational features with the help of electronic structure calculations. Important markers of the inorganic component have been found at 60 and at 94 cm−1, which are mainly associated to the I−Pb−I bending and Pb−I stretching. The region between 100 and 200 cm−1 is instead associated to librational motions of the MA cations. On the basis of electronic structure calculations, we have elucidated how the interactions between the organic cations and the inorganic counterpart may affect the vibrational frequency and the Raman intensity of the MA torsional mode in MAPbI3. The broad and unresolved band at 200−340 cm−1 is consequently assigned to the torsional mode of the MA cations, and we propose this mode as a possible marker of the orientational order of the organic cations in the material, and thus of the whole crystal. We believe the present study to provide the required assignments for the interpretation of the Raman spectra of organohalide perovskites, which may be extremely useful to understand in detail the properties of this class of materials in relation to their full exploitation in solar cells.

Figure 3. Reduced models used to investigate the torsional frequency of the MA cation in MAPbI3. Left the four isolated MA cations at their geometry in the crystal structures. Middle: the four isolated MA cations plus the three nearest neighbor iodide anions at the crystal structures geometries. Right: the periodic crystal models.

Table 1. Vibrational Frequencies (ω) and Raman Intensities Computed on the Isolated Cations with the Same Geometry They Have in the Crystals (Isolated MA), Including the Three Closest Iodides to the NH3+ Group (isolated MA +3 iodide) and on the Periodic Crystal Structuresa isolated MA ω (cm−1)

Raman int. (Å4/amu)

A A′ B B′

465 465 473 473

0.0050 0.0050 0.0011 0.0011

A A′ B B′

467 467 469 469

0.0003 0.0003 0.0003 0.0003

A A′ B B′

495 495 490 490

0.0003 0.0003 0.0044 0.0044

isolated MA + 3 iodide ω (cm−1)

Raman int. (Å4/amu)

Tet-1 344 344 327 327 Tet-2 367 367 355 355 Ortho 437 437 429 429

periodic structures ω (cm−1)

Raman int. (Å4/amu)

28.1350 28.1341 19.2620 19.2618

295 299 312 315

2042.4988 1542.5836 305.6091 128.0141

18.8574 18.9417 18.4804 18.4844

288 296 297 301

814.8004 1436.1303 15.4563 527.7518

1.5272 1.5271 1.5509 1.5509

391 392 398 399

3.3586 42.5448 3.2853 49.2765



ASSOCIATED CONTENT

S Supporting Information *

Sample preparation, details on the Raman spectrum measurement and theoretical and computational details. Additional data and structural analyses. This material is available free of charge via the Internet at http://pubs.acs.org.

a

The AB−A′B′ labels, defined in Figure 3, for the periodic structures refer to normal modes which are maximally localized on one of the MA molecules. Calculations on the isolated models are carried out by Gaussian 09.39



AUTHOR INFORMATION

Corresponding Authors

*E-mail: [email protected]. *E-mail: [email protected]. 282

dx.doi.org/10.1021/jz402589q | J. Phys. Chem. Lett. 2014, 5, 279−284

The Journal of Physical Chemistry Letters

Letter

*E-mail: fi[email protected].

Sensitised Solar Cell Applications. J. Mater. Chem. A 2013, 1, 5628− 5641. (17) 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. (18) Colella, S.; Mosconi, E.; Fedeli, P.; Listorti, A.; Gazza, F.; Orlandi, F.; Ferro, P.; Besagni, T.; Rizzo, A.; Calestani, G.; et al. MAPbI3−xClx Mixed Halide Perovskite for Hybrid Solar Cells: The Role of Chloride as Dopant on the Transport and Structural Properties. Chem. Mater. 2013, 25, 4613−4618. (19) Mosconi, E.; Amat, A.; Nazeeruddin, Md. K.; Grätzel, M.; De Angelis, F. First-Principles Modeling of Mixed Halide Organometal Perovskites for Photovoltaic Applications. J. Phys. Chem. C 2013, 117, 13902−13913. (20) 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. (21) Brivio, F.; Walker, A. B.; Walsh, A. Structural and Electronic Properties of Hybrid Perovskites for High-Efficiency Thin-Film Photovoltaics from First-Principles. Appl. Phys. Lett. 2013, 1, 042111. (22) Onoda-Yamamuro, N.; Matsuo, T.; Suga, H. Calorimetric and IR Spectroscopic Studies of Phase Transitions in Methylammonium Trihalogenoplumbates (II). J. Phys. Chem. Solids 1990, 51, 1383−1395. (23) Khilji, M. Y.; Sherman, W. F.; Wilkinson, G. R. Raman Study of Three Polytypes of PbI2. J. Raman. Spectrosc. 1982, 13, 127−133. (24) Baltog, I.; Lefrant, S.; Mihut, L.; Mondescu, R. Resonant Raman Scattering near the Band Gap in 4H-PbI2 Crystals. Phys. Status Solidi B 1993, 176, 247−254. (25) Calistru, D. M.; Mihut, L.; Lefrant, S.; Baltog, I. Identification of the Symmetry of Phonon Modes in CsPbCl3 in Phase IV by Raman and Resonance-Raman Scattering. J. Appl. Phys. 1997, 82, 5391−5395. (26) Baibarac, M.; Preda, N.; Mihut, L.; Baltog, I.; Lefrant, S.; Mevellec, J. Y. On the Optical Properties of Micro- and Nanometric Size PbI2 Particles. J. Phys.: Condens. Matter 2004, 16, 2345−2356. (27) Maalej, a.; Abid, Y.; Kellel, A.; Daoud, A.; Lautié, A.; Romain, F. Phase Transitions and Crystal Dynamics in the Cubic Perovskite CH3NH3PbCl3. Solid State Commun. 1997, 5, 279−284. (28) Dammak, T.; Fourati, N.; Boughzala, H.; Mlayah, A.; Abid, Y. Xray Diffraction, Vibrational and Photoluminescence Studies of the SelfOrganized Quantum Well Crystal H3N(CH2)6NH3PbBr4. J. Lumin. 2007, 127, 404−408. (29) Elleuch, S.; Abid, Y.; Mlayah, A.; Boughzala, H. Vibrational and Optical Properties of a One-Dimensional Organic−Inorganic Crystal [C6H14N]PbI3. J. Raman Spectrosc. 2008, 39, 786−792. (30) Dammak, T.; Elleuch, S.; Bougzhala, H.; Mlayah, A.; Chtourou, R.; Abid, Y. Synthesis, Vibrational and Optical Properties of a New Three-Layered Organic−Inorganic Perovskite (C4H9NH3)4Pb3I4Br6. J. Lumin. 2009, 129, 893−897. (31) Samet, A.; Ben Ahmed, A.; Mlayah, A.; Boughzala, H.; Hlil, E. K.; Abid, Y. Optical Properties and Ab Initio Study on the Hybrid Organic−Inorganic Material [(CH3)2NH2]3[BiI6]. J. Mol. Struct. 2011, 977, 72−77. (32) Kessentini, A.; Belhouchet, M.; Suñol, J. J.; Abid, Y.; Mhiri, T. Synthesis, Crystal Structure, Vibrational Spectra, Optical Properties and Theoretical Investigation of Bis(2-aminobenzimidazolium) Tetraiodocadmate. J. Mol. Struct. 2013, 1039, 207−213. (33) Waldron, R. D. The Infrared Spectra of Three Solid Phases of Methyl Ammonium Chloride. J. Chem. Phys. 1953, 21, 734−741. (34) Théorêt, A.; Sandorfy, C. The Infrared Spectra of Solid Methylammonium Halide - II. Spectrochim. Acta 1967, 23A, 519−542. (35) Jiang, Z. T.; James, B. D.; Liesegang, J.; Tan, K. L.; Gopalakrishnan, R.; Novak, I. Investigation of the Ferroelectric Phase Transition in (CH3NH3)HgCl3. J. Phys. Chem. Solids 1995, 56, 277−283. (36) DeFrees, D. J.; McLean, A. D. Molecular Orbital Predictions of the Vibrational Frequencies of Some Molecular Ions. J. Chem. Phys. 1985, 82, 333−341.

Notes

The authors declare no competing financial interest.



ACKNOWLEDGMENTS We thank FP7-NMP-2013 project 604032 “MESO” for financial support.



REFERENCES

(1) Bisquert, J. The Swift Surge of Perovskite Photovoltaics. J. Phys. Chem. Lett. 2013, 4, 2597−2598. Snaith, H. J. Perovskites: The Emergence of a New Era for Low-Cost, High-Efficiency Solar Cells. J. Phys. Chem. Lett. 2013, 4, 3623−3630. (2) 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. (3) Im, J.-H.; Lee, C.-R.; Lee; Lee, J.-W.; Park, S.-W.; Park, N.-G. 6.5% Efficient Perovskite Quantum-Dot-Sensitized Solar Cell. Nanoscale 2011, 3, 4088−4093. (4) Chung, I.; Lee, B.; He, J.; Chang, R. P. H.; Kanatzidis, M.; et al. All-Solid-State Dye-Sensitized Solar Cells With High Efficiency. Nature 2012, 485, 486−489. (5) 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, 463−467. (6) Etgar, L.; Gao, P.; Xue, Z.; Peng, Q.; Chandiran, A. K.; Liu, B.; Nazeeruddin, Md. K.; Grätzel, M. Mesoscopic CH3NH3PbI3/TiO2 Heterojunction Solar Cells. J. Am. Chem. Soc. 2012, 134, 17396− 17399. (7) Crossland, E. J. W.; Noel, N.; Sivaram, V.; Leijtens, T.; Alexander-Webber, J. A.; Snaith, H. J. Mesoporous TiO2 Single Crystals Delivering Enhanced Mobility And Optoelectronic Device Performance. Nature 2013, 495, 215−219. (8) Abrusci, A.; Stranks, S. D.; Docampo, P.; Yip, H.-L; Jen, A. K.-J; Snaith, H. J. High-Performance Perovskite-Polymer Hybrid Solar Cells via Electronic Coupling with Fullerene Monolayers. Nano Lett. 2013, 13, 3124−3128. (9) M. Liu, M.; Johnston, M. B.; Snaith, H. J. Efficient Planar Heterojunction Perovskite Solar Cells By Vapour Deposition. Nature 2013, 501, 395−398. (10) J. Burschka, J.; Pellet, N.; Moon, S.-J.; Humphry-Baker, R.; Peng, G.; Nazeeruddin, Md. K.; Grätzel, M. Sequential Deposition as a Route to High-Performance Perovskite-Sensitized Solar Cells. Nature 2013, 499, 316−319. (11) 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.; Nazeeruddin, Md. K.; et al. Efficient Inorganic−Organic Hybrid Heterojunction Solar Cells Containing Perovskite Compound and Polymeric Hole Conductors. Nat. Photon 2013, 7, 486−491. (12) 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, 234, 341− 344. (13) 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. (14) Edri, E.; Kirmayer, S.; Cahen, D.; Hodes, G. High Open-Circuit Voltage Solar Cells Based on Organic−Inorganic Lead Bromide Perovskite. J. Phys. Chem. Lett. 2013, 4, 897−902. (15) Poglitsch, A.; Weber, D. Dynamic Disorder in Methylammoniumtrihalogenoplumbates (II) Observed by Millimeterwave Spectroscopy. J. Chem. Phys. 1987, 87, 6373−6378. (16) Baikie, T.; Fang, Y.; Kadro, J. M; Schreyer, M.; Wei, F.; Mhaisalkar, S. G.; Grätzel, M.; White, T. J. Synthesis and Crystal Chemistry of the Hybrid Perovskite (CH3NH3)PbI3 for Solid-State 283

dx.doi.org/10.1021/jz402589q | J. Phys. Chem. Lett. 2014, 5, 279−284

The Journal of Physical Chemistry Letters

Letter

(37) Zeroka, D.; Jensen, J. O. Infrared Spectra of Some Isotopomers of Methylamine and the Methylammonium Ion: A Theoretical Study. J. Mol. Struct.: THEOCHEM 1998, 425, 181−192. (38) Giannozzi, P.; Baroni, S.; Bonini, N.; Calandra, M.; Car, R.; Cavazzoni, C.; Ceresoli, D.; Chiarotti, G. L.; Coccioni, 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. (39) Frisch, M. J.; Trucks, G. W.; Schlegel, H. B.; Scuseria, G. E.; Robb, M. A.; Cheeseman, J. R.; Scalmani, G.; Barone, V.; Mennucci, B.; Petersson, G. A., et al. Gaussian 09, revision C.01; Gaussian, Inc.: Wallingford, CT, 2009. (40) Zhong, W.; King-Smith, R. D.; Vanderbilt, D. Giant LO−TO Splittings in Perovskite Ferroelectrics. Phys. Rev. B 1994, 72, 3618− 3621. (41) Waghmare, U. V.; Rabe, K. M. Ab Initio Statistical Mechanics of the Ferroelectric Phase Transition in PbTiO3. Phys. Rev. B 1997, 55, 6161−6173. (42) Corà, F.; Catlow, C. R. A. QM Investigations on PerovskiteStructured Transition Metal Oxides: Bulk, Surfaces and Interfaces. Faraday Discuss. 1999, 114, 421−442. (43) Ghosez, Ph.; Cockayne, E.; Waghmare, U. V.; Rabe, K. M. Lattice Dynamics of BaTiO3, PbTiO3, and PbZrO3: A Comparative First-Principles Study. Phys. Rev. B 1999, 60, 836−843. (44) Piskunov, S.; Heifets, E.; Eglitis, R. I.; Borstel, G. Bulk Properties and Electronic Structure of SrTiO3, BaTiO3, PbTiO3 Perovskites: An Ab Initio HF/DFT Study. Comput. Mater. Sci. 2004, 29, 165−178. (45) Wahl, R.; Vogtenhuber, D.; Kresse, G. SrTiO3 and BaTiO3 Revisited Using the Projector Augmented Wave Method: Performance of Hybrid and Semilocal Functionals. Phys. Rev. B 2008, 78, 104116. (46) Evarestov, R. A. Hybrid Density Functional Theory LCAO Calculations on Phonons in Ba(Ti,Zr, Hf,)O3. Phys. Rev. B 2011, 83, 014105. (47) Moreira, E.; Henriques, J. M.; Alzevedo, D. L.; Caetano, E. W. S.; Freire, V. N.; Albuquerque, E. L. Structural, Optoelectronic, Infrared and Raman Spectra of Orthorhombic SrSnO3 from DFT Calculations. J. Solid State Chem. 2011, 184, 921−928. (48) Pyykkö, P. Relativistic Effects in Structural Chemistry. Chem. Rev. 1988, 88, 563−594. (49) Berces, A. Harmonic Vibrational Frequencies and Force Constants of M(CO)5CX (M) Cr, Mo, W; X) O, S, Se). The Performance of Density Functional Theory and the Influence of Relativistic Effects. J. Phys. Chem. 1996, 100, 16538−16544. (50) Barone, V. Vibrational Zero-Point Energies and Thermodynamic Functions beyond the Harmonic Approximation. J. Chem. Phys. 2004, 120, 3059−3065. (51) Mitzi, D. B. Templating and Structural Engineering in Organic− Inorganic Perovskites. J. Chem. Soc., Dalton Trans. 2001, 1−12.

284

dx.doi.org/10.1021/jz402589q | J. Phys. Chem. Lett. 2014, 5, 279−284