Revealing the role of organic cations in hybrid halide perovskite CH3NH3PbI3

Abstract

The hybrid halide perovskite CH₃ NH₃ PbI₃ has enabled solar cells to reach an efficiency of about 20%, demonstrating a pace for improvements with no precedents in the solar energy arena. Despite such explosive progress, the microscopic origin behind the success of such material is still debated, with the role played by the organic cations in the light-harvesting process remaining unclear. Here van der Waals-corrected density functional theory calculations reveal that the orientation of the organic molecules plays a fundamental role in determining the material electronic properties. For instance, if CH₃ NH₃ orients along a (011)-like direction, the PbI₆ octahedral cage will distort and the bandgap will become indirect. Our results suggest that molecular rotations, with the consequent dynamical change of the band structure, might be at the origin of the slow carrier recombination and the superior conversion efficiency of CH₃ NH₃ PbI₃ .

ARTICLE

Received 1 May 2014 | Accepted 25 Mar 2015 | Published 27 Apr 2015

The hybrid halide perovskite CH3NH3PbI3 has enabled solar cells to reach an efciency of about 20%, demonstrating a pace for improvements with no precedents in the solar energy arena. Despite such explosive progress, the microscopic origin behind the success of such material is still debated, with the role played by the organic cations in the light-harvesting process remaining unclear. Here van der Waals-corrected density functional theory calculations reveal that the orientation of the organic molecules plays a fundamental role in determining the material electronic properties. For instance, if CH3NH3 orients along a(011)-like direction, the PbI6 octahedral cage will distort and the bandgap will become indirect. Our results suggest that molecular rotations, with the consequent dynamical change of the band structure, might be at the origin of the slow carrier recombination and the superior conversion efciency of CH3NH3PbI3.

DOI: 10.1038/ncomms8026 OPEN

Revealing the role of organic cationsin hybrid halide perovskite CH3NH3PbI3

Carlo Motta1,2, Fedwa El-Mellouhi2, Sabre Kais2,3, Nouar Tabet2, Fahhad Alharbi2 & Stefano Sanvito1

1 School of Physics, AMBER and CRANN Institute, Trinity College, Dublin Dublin 2, Ireland. 2 Qatar Environment and Energy Research Institute, PO Box 5825, Doha, Qatar. 3 Department of Chemistry and Physics, Purdue University, West Lafayette, Indiana 46323, USA. Correspondence and requests for materials should be addressed to C.M. (email: mailto:[email protected]

Web End [email protected] ) or to F.E.-M. (email: mailto:[email protected]

Web End [email protected] ).

NATURE COMMUNICATIONS | 6:7026 | DOI: 10.1038/ncomms8026 | http://www.nature.com/naturecommunications

Web End =www.nature.com/naturecommunications 1

& 2015 Macmillan Publishers Limited. All rights reserved.

ARTICLE NATURE COMMUNICATIONS | DOI: 10.1038/ncomms8026

Very few materials have taken a research eld by storm as the hybrid halide perovskites have done in the last 2 years in the solar cell community1. These are compounds with

the standard AMX3 perovskite structure, where X is the halide iodine, M is lead, while the remaining cation position, A, is taken by an organic molecule, in this case methylammonium (MA), CH3NH3. The success of such compounds is that high-efciency solar cells can be fabricated cheaply from the liquid phase2,3 and that the efciency can be tuned by controlling their structural order and composition4,5. Particularly intriguing is the fact that high efciencies are achieved even for planar cells6, indicating that the charge separation can occur in the hybrid perovskite absorber and that efcient charge diffusion takes place for both electrons and holes7.

Part of this behaviour can be related to the unusual situation concerning possible intrinsic defects. Recent electronic structure calculations8,9 have shown that the dominant simple or complex defects in CH3NH3PbI3 produce only shallow defect levels, while those able to act as deep traps have a high formation energy. Such observation suggests that, under ideal conditions, deep-level point defects cannot act as non-radiative recombination sitesa fact that can partially explain the relative long exciton lifetimes and diffusion lengths observed in these materials. At the same time, the absence of deep levels may explain also the large open-circuit voltage measured. However, despite all these hints, many fundamental questions remain open. In particular, it is still not clear what is the role played by the organic molecules in the efciency of these compounds. Theoretical investigations have now built a consensus that the MA group does not have any signicant contribution to the electronic structure around the band edges8,10,11 and suggest that the only role of the molecules is to stabilize electrostatically the perovskite structure.

Interestingly, not even the geometrical arrangement of the MA group is known in detail. At low temperature (below 150 K), the crystal adopts an orthorhombic (Pnma) structure, in which the PbI6 octahedra are strongly deformed assuming a rectangular basal plane. Such deformation restricts the rotational degrees of freedom of MA in the rhombus-shaped interstitial region, thus imposing a spatial ordering to CH3NH3. In this case, the organic cation is pinned and can only rotate along the CN axis. As the temperature is increased and the material progressively assumes a cubic structure (Pm-3m), by passing through a tetragonal one (I4/m), the CH3NH3 molecules become free to rotate between the octahedral cages. Above room temperature, such rotation is fast, to a point where both crystallographic analysis3,12 and NMR measurements13 have shown that the exact location of the MA groups cannot be determined.

Here we conduct density functional theory (DFT) calculations for CH3NH3PbI3 in the cubic phase, by taking into consideration the experimentally reported evidence of the fast rotation of CH3NH3. This means that we do not limit our analysis to

CH3NH3 oriented along the (100) or (111) direction, for which a higher symmetry is maintained, but we also explore cases where the symmetry is lowered. We nd that such symmetry lowering has profound consequences on the electronic structure, namely that the bandgap changes from direct to indirect. Crucially, such symmetry-lowering congurations represent local minima in the free energy surface of the crystal and they are stabilized by van der Waals (vdW) interactions. These are the key ingredients not only for obtaining accurate lattice parameters12,1417 (see Supplementary Fig. 1) but also for the internal geometry. Our calculations then return a picture of the CH3NH3PbI3 as a dynamical bandgap semiconductor, in which the exact position of the conduction band minimum depends on the particular spatial arrangement of the molecules. Importantly, our results are robust against bandgap corrections and spinorbit

interaction11,1821, and deliver an absorption spectrum in good agreement with experimental data near the absorbtion edge7,22.

ResultsGeometrical optimization. Let us start the discussion by presenting the relaxed crystal geometries. Initially, all the calculations have been performed at the level of the generalized gradient approximation (GGA)23 to the DFT exchange and correlation functional, including vdW interactions24 (see computational details section) and without explicitly considering the spin orbit coupling (SOC). Additional calculations beyond the GGA and including SOC will be discussed later in relation to the electronic structure.

In Fig. 1, the relaxed structures for the cubic phase of CH3NH3PbI3 are presented along the xy and xz planes for two different orientations of the CH3NH3 cation. Given the periodic boundary conditions, these structures effectively represent the specic case of an innite array of perfectly aligned molecular dipoles. The relaxation process is extremely sensitive to the initial conditions. When we start the geometrical optimization with the molecular cation oriented along the (111) direction, the structure relaxes maintaining the same orientation (Fig. 1a). In contrast, by starting from a (001)-oriented molecule, the relaxation may end up with CH3NH3 along (011) (Fig. 1b) or an equivalent direction.

However, it is important to point out that the cation energy landscape is very shallow. Indeed, both congurations represent local energy minima, with an energy difference of only about 20 meV in favour of (011). In both cases, the lattice parameter is calculated in excellent agreement with the experimental value of6.329 (ref. 25) and an error of o0.8% (see Supplementary Table 1 and Supplementary Fig. 1). If one does not include vdW forces in the relaxation, the discrepancy with the nominal value will become 42%, that is, the crystal structure is not properly described. Furthermore, the angle formed by the molecular axis and the (001) direction increases from 16 to 25 on inclusion of dispersion forces. To make a well-grounded comparison with previous calculations11,26, reporting no effect of MA on the structure, we started the relaxation from the X-ray diffraction solved structure with CH3NH3 along (100)27, whose band structure is presented in Supplementary Fig. 2.

By inspecting the relaxed structures, a prominent difference between the (111) and the other two geometries appears. In fact, in the case of (011), the rotation of the cation induces a deformation and a symmetry reduction of the inorganic PbI6 octahedra. The result of such PbI6 octahedra distortion is that the PbI bonds do not lie parallel to the crystal directions, but instead they slightly deviate by B6. Unlike the (111) case, the unit cell becomes pseudocubic with the lattice parameters along the y and z directions being 1% larger than that along x. Hence, it is clear that the dispersive forces are critical for the internal geometry optimization and consequently, as we will see next, for the electronic behaviour of hybrid perovskites.

Electronic properties. The ne differences in the crystal geometries have a qualitatively strong impact on the electronic structure. In Fig. 2, we show the electronic band structure along some high-symmetry points of the Brillouin zone calculated with and without including fully relativistic SOC. For both orientations, the energy gap computed without SOC lies in the1.51.7 eV range, close to the experimental value of about 1.55 eV (refs 22,28). However, as already shown5, this agreement is due to a fortuitous cancellation of errors, namely the GGA underestimates the bandgap, but this is counterbalanced by the lack of the spinorbit interaction. It is well known that the bottom of the conduction band mainly originates from the p

2 NATURE COMMUNICATIONS | 6:7026 | DOI: 10.1038/ncomms8026 | http://www.nature.com/naturecommunications

Web End =www.nature.com/naturecommunications

& 2015 Macmillan Publishers Limited. All rights reserved.

NATURE COMMUNICATIONS | DOI: 10.1038/ncomms8026 ARTICLE

a b

z

x

y

x

y

z

Figure 1 | Relaxed structures of the cubic phase of CH3NH3PbI3. Relaxed structures of the cubic phase of CH3NH3PbI3 for two different orientations of the cations, namely along (a) (111) and (b) (011). (b) The views along the xy and xz planes are shown in the upper and lower panels, respectively.

orbitals of Pb, while the top of the valence band is derived from the p orbitals of I. The highest occupied molecular orbital of MA is found deep below the valence band, B5 eV below the valence band maximum (VBM). Thus, one may argue that CH3NH3 does not play any role in the optical and electronic response of such materials, but that rather it does only contribute to their structural cohesion. However, a careful analysis of the density of states projected on the various atoms (see Fig. 3a and Supplementary Fig. 3) reveals that there is a small contribution B0.5 eV below the VBM attributable to the organic molecules.

This indicates that indeed there is an interaction between CH3NH3 and the inorganic PbI6 octahedra, in the form of hydrogen bonds between N and I.

A closer inspection to the band structures reveals that the orientation of CH3NH3 has a profound impact on the nature of the bandgap. In fact, while for the (111) case the gap is direct at the R point, the cationic rotation induces a transition to an indirect gap, and when MA is along (011) the conduction band mininum (CBM) shifts along the R-G line. Accordingly, the vertical energy difference between the CBM and the conduction band energy at R is 25 meV, as shown in the inset of Fig. 2b. The effect persists when SOC is turned on, as shown in Fig. 2c,d. In this case, a slight shift of the VBM towards G occurs as well, reducing the CB vertical energy difference approximately to half the non-SOC value. We stress here that the (011) orientation is not a special casethe indirect bandgap does appear for all the equivalent orientations, namely (011), (011), (101), (110), (101) and (011). As an example, in Supplementary Fig. 4, we show the (011) case, for which the CBM shifts in the direction opposite to that shown in Fig. 2, namely R-M.

Such a gap transformation is the most remarkable result of the interaction between the organic molecules and the inorganic framework, which is driven by the dispersive forces. As a result, CH3NH3 exerts strain on the PbI6 octahedra and changes the

nature of the bandgap. This is unusual in the perovskite world and it is a unique peculiarity of such hybrid systems. Typical inorganic perovskites such as SrTiO3 respond to very small uniaxial strains, as those induced by the antiferrodistortive phase transition, by splitting their threefold degenerate CBM29. To the best of our knowledge, very few semiconductors see their bandgap changing from direct to indirect under the effect of strain. An example is given by layered transition metal chalcogenides, like for instance, MoSe2 single layer, where the gap change is attributed to phonon softening30. Intriguingly, also for CH3NH3PbI3, we observe some softening of our calculated low-energy optical phonon modes, when the MA molecule rotates along the (011) directions. In particular, the calculated energy red shift is of the order of 10 cm 1 (see Supplementary Fig. 5).

We have carefully veried that the indirect bandgap is not an artifact of our computation method. First, we analysed the relaxation in a 2 2 2 (96 atoms) supercell to establish whether

further low-symmetry congurations are possible. The supercell converged geometries are very similar to the unit cell geometries repeated in space due to periodic boundary conditions. The calculations reveal no particular changes, although the size of the supercell would allow the PbI6 octahedra to rotate or tilt.

Consequently, the features of the bandgaps are preserved both at the qualitative as well as the quantitative level.

Second, we have tested different vdW functionals by comparing the geometries obtained with the TkatchenkoSchefer (TS) scheme24 to those calculated with the empirical Grimmes method31 and with a recently developed sophisticated many-body dispersion correction32. As shown in Supplementary Fig. 1, the unit cell volume compares well with the experimental value for all the vdW-corrected functionals, although that calculated with the TS scheme is closer. Thus, in general, the TS and Grimmes pairwise dispersion methods provide a sufciently good description of our system as found in a recent work on

NATURE COMMUNICATIONS | 6:7026 | DOI: 10.1038/ncomms8026 | http://www.nature.com/naturecommunications

Web End =www.nature.com/naturecommunications 3

& 2015 Macmillan Publishers Limited. All rights reserved.

ARTICLE NATURE COMMUNICATIONS | DOI: 10.1038/ncomms8026

molecular crystals17. In any case, the indirect bandgap is found for all the calculated geometries, including those obtained with the many-body dispersion correction. In contrast, when the dispersion interactions are not included, the gap remains direct since the structure does not present enough distortion. Then, we have performed calculations using the HSE06 (ref. 33) screened hybrid functional by starting from the relaxed structures. Also, in this case, the bandgap remains indirect (see Supplementary Fig. 3), although the overall energy distance between the CBM and the VBM increases, as expected.

Finally, we have performed additional calculations to establish whether or not the observed effect is the result of the interaction of the MA dipole with its periodic images. We have then veried that the position of the VBM and the CBM does not shift (i) on ipping the molecules dipole moment and (ii) on applying a rigid rotation of the molecule from the 111 to the 110 orientation (Supplementary Fig. 6). This conrms that it is the distortion of the inorganic PbI6 octahedra to govern the shape of the CBM and

VBM and not the dipole moment of the organic molecule.

As a last element of our analysis, we have investigated the relative stability of the different molecular orientations. In particular, we have performed nudged elastic band calculations and evaluated the barrier height for an MA rotation from (011) to (111), taking into account the possible changes in the cell geometry and the lattice vectors. We have found a barrier of

20 meV (1.929 kJ mol 1), which is less than kBT, and it is similar to values found recently in literature34. This conrms the concept that above room temperature the molecules continuously rotate. Given the strong dependence of the position in k-space of the CBM and the VBM on the cell geometry, the fast molecular rotation implies that the nature and shape of the bandgap are rapidly-varying functions of time. Recently, Mosconi et al.35 published CarParinello molecular dynamics simulations for a 2 2 2 CH3NH3PBI3 supercell at 319 K, showing uctuations

of 0.10.2 eV in the position of the VBM, and thus giving evidence that the size of the bandgap changes as function of the MA rotation. Here we go a step forward and demonstrate that for some molecular orientations the bandgap turns from direct to indirect.

Optical absorption. To seek further evidence of the presence of the indirect bandgap, we have calculated the perovskites optical absorption and compared it with available experimental data. The calculations have been performed over the electronic structure obtained with the GGA including vdW interaction in the independent particle approximation and without SOC. Although the fortuitous cancellation of errors returns a bandgap close to experiments, we concentrate here on the shape of the spectrum around the absorption edge. The calculations are performed for the two orientations of the CH3NH3 molecule, and the spectra are presented in Fig. 3c.

When MA is aligned along the (111) direction, the bandgap is direct and measures 1.423 eV 871 nm. The absorption prole

2

DOS (a.u.)PDOS (a.u.)

25 meV

I CH3NH3

Pb

Energy (eV)

1

0

(111)(011)

1 M R

R

M

1.5

4 2 0 2 4 6 Energy (eV)

1

0.20

0.15

0.10

0.05

0.00

700 750 800 850 900

2.0

Absorption (a.u.)

0.5

1.5

Energy (eV)

1.0

0

0.5

0.5

0.0

400 500 600 700 800 900Wavelength (nm)

1 M R

R

M

Figure 2 | Band structure of the fully relaxed CH3NH3PbI3 crystal. The bands are shown for molecule orientations along (a) (111) and (b) (011) direction. The insets show a magnication of the bands (which have been shifted in energy for convenience) around the bandgap and highlight the changes in the VBM and CBM caused by the rotation of CH3NH3. Note that for the (011) orientation the bandgap becomes indirect. The relativistic SOC band structures are shown in panels (c) and (d) for (111) and (011).

Figure 3 | Electronic and optical properties of the cubic phase CH3NH3PBI3. (a) Projected density of states (DOS) for the case of the (111)-oriented molecule on the Pb (cyan), I (black) species and the CH3NH3 (red). (b) DOS as a function of the chemical potential for the molecule oriented along the (111) and (011) directions are displayed. (c) The calculated absorption spectra are plotted for the different orientation of CH3NH3 clearly showing the impact the organic cation rotation. The inset in c is a magnication of the same quantity around the low-energy region.

4 NATURE COMMUNICATIONS | 6:7026 | DOI: 10.1038/ncomms8026 | http://www.nature.com/naturecommunications

Web End =www.nature.com/naturecommunications

& 2015 Macmillan Publishers Limited. All rights reserved.

NATURE COMMUNICATIONS | DOI: 10.1038/ncomms8026 ARTICLE

then increases monotonically from that value, indicating the possibility of direct transitions across the band edges. In contrast, when CH3NH3 is placed along (011) the bandgap is indirect, and in general slightly larger than that corresponding to the (111) orientation. In particular, for (011), we nd a direct gap of1.611 eV 769 nm and an indirect one of 1.629 eV 761 nm.

Here the shape of the absorption spectrum is quite different from the (111) case. We note a two-step absorption prole (see the inset of Fig. 3c), with a small absorption amplitude developing at the band edge followed by a signicantly more pronounced increase at around 750 nm. This reects the indirect nature of the bandgap and the fact that indirect transitions are not allowed, unless assisted by phonons. Certainly, the level of theory used here is not complete and one should use the many-body methods for deriving a more quantitative picture36; however, the inclusion of excitonic effects is not necessary to make our main point about the shape of the spectrum around the absorption edge. In fact, the differences between the direct and indirect bandgap are very evident from the inset of Fig. 3c, and should be compared with the typical experimental absorption spectrum measured at 300 K (ref. 4). This latter features a small amplitude around the bandgap followed by a rapid increase after about 50 nm, which seems to be a robust feature in the experimentally measured spectra, regardless of the morphology of the material7,22,37 and resistent to lm quality degradation38.

Whether or not such a near-edge feature can be attributed to the indirect bandgap is however difcult to establish with certainty. In fact, at 300 K, the MAs are free to rotate and sample the many quasi-equivalent (011) and (111) orientations. Thus, it is expected that no sharp indirect bandgap can be observed, but rather a smooth absorption spectra starting around the average bandgap value. Interestingly, we point out that our nding could help explaining the recently observed line broadening in photoluminescence experiments39, attributed to coupling with phonons. In fact, low-energy phonons are associated with cation rotations, which in turn would shift the CBM to different positions. As a result, the range of optical gap would change as a function of temperature, thus broadening the emission spectrum. Measurement of the optical absorption and thermal characteristics using very high-sensitivity techniques such as photothermal deection spectroscopy might reveal the existence of the indirect bandgap by ne lowering of temperature o300 K thus freezing the MA (011)-like orientations.

DiscussionAbsorbing light efciently and transporting the photocurrent with some gained voltage are the most important aspects in the design of a solar cell. The two processes can be characterized, respectively, by an empirical absorption length, La, inversely proportional to the absorption40,41, and by the carrier diffusion length, Ldiff. This is a measure of the mean distance travelled by carriers before recombining4244; it is directly proportional to the carrier lifetime, t, and it depends on several factors (mobility, carrier concentration and so on). In solar cells, it is essential to full the condition r Ldiff/La41 to ensure that the

photogenerated carriers can be extracted and collected. For practical reasons, rc1 is required for high power conversion efciency4042,45.

In general, direct bandgap semiconductors, such as GaAs, have strong absorption and thus short La, but this comes with fast carrier recombination and hence relatively small Ldiff. The

situation is opposite when the bandgap in indirect (for example, in Si). An ideal balance between La and Ldiff is present in some indirect bandgap semiconductors such as Cu3N (refs 41,42), where the direct transition starts a few kBT above the indirect gap,

yielding a short La. The small difference between the indirect gap and the lowest direct transition ensures that the absorption is not greatly affected by the indirect gap till the direct transitions start to dominate the absorption spectrum. In contrast, the lifetime of the photogenerated carriers, which relax to the band edges in the indirect gap semiconductors, is longer than that of the direct ones by orders of magnitudes43 ensuring a large Ldiff. This class of indirect semiconductors promises remarkable improvement in efciency via the large Ldiff/La ratio41,42.

From this short discussion and our results at hand, we can speculate on the possible origin of the high efciency in hybrid perovskites solar cells. Our calculations establish that, under certain conditions of molecule alignment, the bandgap of the high-temperature phase of CH3NH3PBI3 can be indirect and not direct. Note that an indirect bandgap in hybrid organic/inorganic perovskite may also arise as a result of Rashba splitting10,46. This feature originates from the strain exerted on the PbI6 octahedra by the organic molecules, when these points mostly towards one of the unit cell faces rather than along the diagonal, which constitute shallow minima in the total energy prole. The indirect bandgap is tens of meV smaller than the direct one and the position of the CBM with respect to the VBM depends on the specic orientation of the molecules. This analysis delivers a picture of CH3NH3PBI3 as a dynamical bandgap semiconductor, where the term dynamical refers to the fact that the position in k-space of the CBM depends on the molecular alignment, which in turn uctuates at room temperature. The calculated absorption spectrum near the band edge, characterized by a shallow amplitude at about 800 nm followed by a rapid enhancement within 50 nm, is very similar for (011) and equivalent orientations of the molecules, indicating that it is a feature robust to the molecular rotation and disorder.

The presence of an indirect bandgap, with only a small displacement in k-space with respect to the direct point, suggests that the rapid relaxation of the photoexcited electrons to the conduction band minimum is still possible. This process is incoherent, mediated by acoustic phonons and spontaneously breaks the excitons in free carriers. As a second effect of the collapse of the electrons to the indirect conduction band minimum is that the lifetime of the minority carriers is enhanced, so that radiative recombination is partially suppressed. Such suppression is not complete since the experiments reveal luminescence; however, it contributes to prolong the carrier lifetime. In fact, it is interesting to note that experiments report that the luminescence gets suppressed as the material crystallizes37 and that in any case it is sensitive to the sample morphology and quality47. Considering that the crystallization process leads to the geometries described here, such suppression of the luminescence on ordering supports our argument. To complete the picture, it is also worth reminding that the absence of strong scattering centres inhibits incoherent recombination, so that Ldiff remains long48.

ConclusionsDespite the tremendous advances in studying and understanding hybrid perovskites from both the theoretical and the experimental point of views, yet the reason behind the extended carrier lifetime in these materials remains a mystery. Our work suggests that the extended carrier lifetimes in CH3NH3PBI3 might be a consequence of the dynamical position of the conduction band minimum forming a basin around a high-symmetry point in the k-space. This proposed dynamical bandgap, induced by the molecular motion, prohibits carrier recombination and favours exciton separation, thus increasing the conversion efciency of the photvoltaic absorber CH3NH3PBI3. Crucially, all such

NATURE COMMUNICATIONS | 6:7026 | DOI: 10.1038/ncomms8026 | http://www.nature.com/naturecommunications

Web End =www.nature.com/naturecommunications 5

& 2015 Macmillan Publishers Limited. All rights reserved.

ARTICLE NATURE COMMUNICATIONS | DOI: 10.1038/ncomms8026

phenomenology appears to be initiated by the geometry assumed by the molecules following full geometry relaxation. A set of differently oriented CH3NH3 cation geometries, relaxed without any symmetry constrains using DFT including vdW forces, lead to geometries where the PbI6 octahedra distort. The band structure corresponding to such geometries present an indirect bandgap, which persists on bandgap corrections via screened hybrid functionals and SOC. Finally, we propose that a possible way for engineering the efciency of hybrid perovskites is to induce distortion to the PbI6 octahedra either by using strain or organic molecules strongly interacting with the inorganic cages.

Methods

DFT calculations. The electronic properties of the hybrid perovskites have been calculated with the all-electron FHI-aims code at the level of GGA in the Perdew BurkeErnzerhof23 parameterization. Long-range vdW interactions have been taken into account via the TS scheme24, which is also constructed over a GGA and a pairwise dispersive potential. The reciprocal space integration was performed over an 8 8 8 Monkhorst-Pack grid49, after performing convergence tests on

the total energy and the forces. A preconstructed high-accuracy all-electron basis set of numerical atomic orbitals was employed, as provided by the FHI-aims tight default option. This is enough to achieve well-converged forces as conrmed by comparing the results with those obtained with the tier2 basis set. The atomic zeroth-order regular approximation was applied to treat relativistic effects. Unconstrained structural optimization was performed with the BroydenFletcher GoldfarbShanno algorithm50 with with a tolerance of 10 3 eV 1.

Additional calculations, in particular, concerning the use of the hybrid functional HSE06 (ref. 33), were performed using the Gaussian suite51, with the periodic boundary condition52 code version. The Def2 (ref. 53) series of GAUSSIAN

basis sets were optimized following our own procedure, described in ref. 54, for the atomic species of interest in this work. Most numerical settings in Gaussian were left at the default values, for example, uncontrained geometry optimization settings, integral cutoffs and self-consistent convergence thresholds. The standard root mean squared force threshold in Gaussian for geometry optimizations is450 10 6 hartrees per Bohr. The default k-point mesh of 8 8 8 was used for

the 12 atoms unit cell of CH3NH3PbI3.

Finally calculations involving the spinorbit interaction have been carried out using the Vienna ab initio simulation package (VASP). A plane wave basis set energy cutoff of 520 eV was considered in the calculations and Brillouin zone integration was performed using 8 8 8 k-point mesh centred at the Gamma

point. Dispersion correction to the PerdewBurkeErnzerhof was done within the method of Tkatchenko and Schefer (DFT-TS)recently implemented in the 5.3.5 version of VASPused to account for vdW interactions. The projector augmented wave pseudopotentials were used for all elements 5d106s26p2 valence electron potential was used for the Pb atom. With the same set-up, VASP has been used to compute the band structures with fully relativistic SOC.

The linear optical absorption has been obtained with FHI-aims by computing the momentum matrix elements in a 5-eV energy window around the Fermi level and then by using the Fermi golden rule approximation. A dense k-point sampling of 15 15 15 has been employed to ensure converged spectra.

Phonon calculations. Phonon spectra were computed with the frozen-phonon approach for a fully relaxed 2 2 2 supercells using the Phonopy55 code in

combination with FHI-aims A displacement of 0.001 was applied to each atom in the three space directions to compute the dynamical matrix.The phonon DOS used the G point sampling, while the full phonon band structure was generated using the G, M, X and R high-symmetry points in the Brillouin zone.

References

1. Snaith, H. J. Perovskites: the emergence of a new era for low-cost, high-efciency solar cells. J. Phys. Chem. Lett. 4, 36233630 (2013).

2. Boix, P. P., Nonomura, K., Mathews, N. & Mhaisalkar, S. G. Current progress and future perspectives for organic/inorganic perovskite solar cells. Mater. Today 17, 1623 (2014).

3. Loi, M. A. & Hummelen, J. C. Hybrid solar cells: perovskites under the sun. Nat. Mater. 12, 10871089 (2013).

4. Deschler, F. et al. High photoluminescence efciency and optically pumped lasing in solution-processed mixed halide perovskite semiconductors. J. Phys. Chem. Lett. 5, 14211426 (2014).

5. Mosconi, E., Amat, A., Nazeeruddin, M. K., Gratzel, M. & De Angelis, F. First-principles modeling of mixed halide organometal perovskites for photovoltaic applications. J. Phys. Chem. C 117, 1390213913 (2013).

6. Liu, M., Johnston, M. B. & Snaith, H. J. Efcient planar heterojunction perovskite solar cells by vapour deposition. Nature 501, 395398 (2013).

7. Xing, G. et al. Long-range balanced electron- and hole-transport lengths in organic-inorganic CH3NH3PbI3. Science 342, 344347 (2013).

8. Yin, W.-J., Shi, T. & Yan, Y. Unusual defect physics in CH3NH3PbI3 perovskite solar cell absorber. Appl. Phys. Lett. 104, 063903 (2014).

9. 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. 5, 13121317 (2014).

10. Brivio, F., Butler, K. T., Walsh, A. & van Schilfgaarde, M. Relativistic quasiparticle self-consistent electronic structure of hybrid halide perovskite photovoltaic absorbers. Phys. Rev. B 89, 155204 (2014).

11. Umari, P., Mosconi, E. & De Angelis, F. Relativistic gw calculations on CH3NH3PbI3 and CH3NH3SnI3 perovskites for solar cell applications. Sci. Rep.

4, 4467 (2014).12. Baikie, T. et al. Synthesis and crystal chemistry of the hybrid perovskite (CH3NH3)PbI3 for solid-state sensitised solar cell applications. J. Mater. Chem. A 1, 56285641 (2013).

13. Knop, O., Wasylishen, R. E., White, M. A., Cameron, T. S. & Oort, M. J. M. V. Alkylammonium lead halides. Part 2. CH3NH3PbX3 (X Cl, Br, I)

perovskites: cuboctahedral halide cages with isotropic cation reorientation. Can.J. Chem. 68, 412422 (1990).14. Wang, Y. et al. Density functional theory analysis of structural and electronic properties of orthorhombic perovskite CH3NH3PbI3. Phys. Chem. Chem. Phys.

16, 14241429 (2014).15. Menndez-Proupin, E., Palacios, P., Wahnn, P. & Conesa, J. C. Self-consistent relativistic band structure of the CH3NH3PbI3 perovskite. Phys. Rev. B 90, 045207 (2014).

16. Egger, D. A. & Kronik, L. Role of dispersive interactions in determining structural properties of organic inorganic halide perovskites: insights from

rst-principles calculations. J. Phys. Chem. Lett. 5, 27282733 (2014).17. Kronik, L. & Tkatchenko, A. Understanding molecular crystals with dispersion-inclusive density functional theory: pairwise corrections and beyond. Acc. Chem. Res. 47, 32083216 (2014).

18. Even, J., Pedesseau, L., Jancu, J.-M. & Katan, C. Importance of spinorbit coupling in hybrid organic/inorganic perovskites for photovoltaic applications.J. Phys. Chem. Lett. 4, 29993005 (2013).19. Even, J., Pedesseau, L., Jancu, J.-M. & Katan, C. DFT and KP modelling of the phase transitions of lead and tin halide perovskites for photovoltaic cells. Phys. Status Solidi RRL 8, 3135 (2014).

20. Feng, J. & Xiao, B. Crystal structures, optical properties, and effective mass tensors of CH3NH3PbX3 (X I and Br) phases predicted from HSE06. J. Phys.

Chem. Lett. 5, 12781282 (2014).21. Giorgi, G., Fujisawa, J.-I., Segawa, H. & Yamashita, K. Cation role in structural and electronic properties of 3d organic-inorganic halide perovskites: a DFT analysis. J. Phys. Chem. C 118, 1217612183 (2014).

22. Lee, M. M., Teuscher, J., Miyasaka, T., Murakami, T. N. & Snaith, H. J. Efcient hybrid solar cells based on meso-superstructured organometal halide perovskites. Science 338, 643647 (2012).

23. Perdew, J. P., Burke, K. & Ernzerhof, M. Generalized gradient approximation made simple. Phys. Rev. Lett. 77, 38653868 (1996).

24. Tkatchenko, A. & Schefer, M. Accurate molecular van der waals interactions from ground-state electron density and free-atom reference data. Phys. Rev. Lett. 102, 073005 (2009).

25. Poglitsch, A. & Weber, D. Dynamic disorder in methylammoniumtrihalogenoplumbates (ii) observed by millimeter-wave spectroscopy. J. Chem. Phys. 87, 6373 (1987).

26. Brivio, F., Walker, A. B. & Walsh, A. Structural and electronic properties of hybrid perovskites for high-efciency thin-lm photovoltaics from rst-principles. APL Mater. 1, 042111 (2013).

27. 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. 52, 90199038 (2013).

28. Kojima, A., Teshima, K., Shirai, Y. & Miyasaka, T. Organometal halide perovskites as visible-light sensitizers for photovoltaic cells. J. Am. Chem. Soc. 131, 60506051 (2009).

29. El-Mellouhi, F., Brothers, E. N., Lucero, M. J., Bulik, I. W. & Scuseria, G. E. Structural phase transitions of the metal oxide perovskites SrTiO3, LaAlO3, and LaTiO3 studied with a screened hybrid functional. Phys. Rev. B 87, 035107 (2013).

30. Horzum, S. et al. Phonon softening and direct to indirect band gap crossover in strained single-layer MoSe2. Phys. Rev. B 87, 125415 (2013).

31. Grimme, S. Semiempirical GGA-type density functional constructed with a long-range dispersion correction. J. Comput. Chem. 27, 17871799 (2006).

32. Tkatchenko, A., DiStasio, R. A., Car, R. & Schefer, M. Accurate and efcient method for many-body van der Waals interactions. Phys. Rev. Lett. 108, 236402 (2012).

33. Krukau, A. V., Vydrov, O. A., Izmaylov, A. F. & Scuseria, G. E. Inuence of the exchange screening parameter on the performance of screened hybrid functionals. J. Chem. Phys. 125, 224106 (2006).

34. Frost, J. M. et al. Atomistic origins of high-performance in hybrid halide perovskite solar cells. Nano Lett. 14, 25842590 (2014).

6 NATURE COMMUNICATIONS | 6:7026 | DOI: 10.1038/ncomms8026 | http://www.nature.com/naturecommunications

Web End =www.nature.com/naturecommunications

& 2015 Macmillan Publishers Limited. All rights reserved.

NATURE COMMUNICATIONS | DOI: 10.1038/ncomms8026 ARTICLE

35. Mosconi, E., Quarti, C., Ivanovska, T., Ruani, G. & De Angelis, F. Structural and electronic properties of organo-halide lead perovskites: a combined IR-spectroscopy and ab initio molecular dynamics investigation. Phys. Chem. Chem. Phys. 16, 1613716144 (2014).

36. 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 118, 1156611572 (2014).

37. Burschka, J. et al. Sequential deposition as a route to high-performance perovskite-sensitized solar cells. Nature 499, 316319 (2013).

38. Niu, G. et al. Study on the stability of CH3NH3PbI3 lms and the effect of post-modication by aluminum oxide in all-solid-state hybrid solar cells. J. Mater. Chem. A 2, 705710 (2014).

39. Wehrenfennig, C., Liu, M., Snaith, H. J., Johnston, M. B. & Herz, L. M. Homogeneous emission line broadening in the organo lead halide perovskite CH3NH3PbI3XClX. J. Phys. Chem. Lett. 5, 13001306 (2014).

40. Loferski, J. J. Theoretical considerations governing the choice of the optimum semiconductor for photovoltaic solar energy conversion. J. Appl. Phys. 27, 777784 (1956).

41. Alharbi, F. et al. Abundant non-toxic materials for thin lm solar cells: Alternative to conventional materials. Renew. Energ. 36, 27532758 (2011).

42. Zakutayev, A. et al. Defect tolerant semiconductors for solar energy conversion.J. Phys. Chem. Lett. 5, 11171125 (2014).43. Fonash, S. Solar Cell Device Physics (Academic Press, 2010).44. Law, M., Solley, E., Liang, M. & Burk, D. Self-consistent model of minority-carrier lifetime, diffusion length, and mobility. IEEE Electron. Device Lett. 12, 401403 (1991).

45. Gratzel, M. Photoelectrochemical cells. Nature 414, 338344 (2001).46. Stroppa, A. et al. Tunable ferroelectric polarization and its interplay with spin orbit coupling in tin iodide perovskite. Nat. Commun. 5, 5900 (2014).

47. Choi, J. J., Yang, X., Norman, Z. M., Billinge, S. J. L. & Owen, J. S. Structure of methylammonium lead iodide within mesoporous titanium dioxide: active material in high-performance perovskite solar cells. Nano Lett. 14, 127133 (2013).

48. Stranks, S. D. et al. Electron-hole diffusion lengths exceeding 1 micrometer in an organometal trihalide perovskite absorber. Science 342, 341344 (2013).

49. Monkhorst, H. J. & Pack, J. D. Special points for brillouin-zone integrations. Phys. Rev. B 13, 51885192 (1976).

50. Press, W. H., Teukolsky, S. A., Vetterling, W. T. & Flannery, B. P. Numerical Recipes 3rd edn (Cambridge Univ. Press, 1997).

51. Frisch, M. J. et al. Gaussian 09 Revision {B}.01 (Gaussian Inc., 2009).52. Kudin, K. N. & Scuseria, G. E. Linear-scaling density-functional theory with gaussian orbitals and periodic boundary conditions: efcient evaluation of energy and forces via the fast multipole method. Phys. Rev. B 61, 1644016453 (2000).

53. Weigend, F. & Ahlrichs, R. Balanced basis sets of split valence, triple zeta valence and quadruple zeta valence quality for H to Rn: design and assessment of accuracy. Phys. Chem. Chem. Phys. 7, 32973305 (2005).

54. El-Mellouhi, F., Brothers, E. N., Lucero, M. J. & Scuseria, G. E. Modeling of the cubic and antiferrodistortive phases of SrTiO3 with screened hybrid density functional theory. Phys. Rev. B 84, 115122 (2011).

55. Togo, A., Oba, F. & Tanaka, I. First-principles calculations of the ferroelastic transition between rutile-type and CaCl2-type SiO2 at high pressures. Phys. Rev.

B 78, 134106 (2008).

Acknowledgements

This work is sponsored by the European Research Council, Quest project, (C.M. and S.S.) and by the Qatar Environment and Energy Research Institute (F.E.-M., F.A., N.T. and S.K.). Computational resources have been provided by the supercomputer facilities at the Trinity Center for High Performance Computing, at ICHEC (project tcphy038b) and the research computing at Texas A&M University at Qatar. We acknowledge fruitful discussions with Simon Gelinas, Marcelo Carignano, Mohamed Lamine Madjet and Felix Deschler.

Author contributions

C.M. and F.E.-M. contributed equally to the manuscript. C.M. and F.E.-M. performed the computational research, C.M., F.E.-M., S.K., N.T., F.A. and S.S. discussed the results and contributed to the preparation of the manuscript.

Additional information

Supplementary Information accompanies this paper at http://www.nature.com/naturecommunications

Web End =http://www.nature.com/ http://www.nature.com/naturecommunications

Web End =naturecommunications

Competing nancial interests: The authors declare no competing nancial interests.

Reprints and permission information is available online at http://npg.nature.com/reprintsandpermissions/

Web End =http://npg.nature.com/ http://npg.nature.com/reprintsandpermissions/

Web End =reprintsandpermissions/

How to cite this article: Motta, C. et al. Revealing the role of organic cations in hybrid halide perovskite CH3NH3PbI3. Nat. Commun. 6:7026 doi: 10.1038/ncomms8026 (2015).

This work is licensed under a Creative Commons Attribution 4.0 International License. The images or other third party material in this article are included in the articles Creative Commons license, unless indicated otherwise in the credit line; if the material is not included under the Creative Commons license, users will need to obtain permission from the license holder to reproduce the material. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/

Web End =http://creativecommons.org/licenses/by/4.0/

NATURE COMMUNICATIONS | 6:7026 | DOI: 10.1038/ncomms8026 | http://www.nature.com/naturecommunications

Web End =www.nature.com/naturecommunications 7

& 2015 Macmillan Publishers Limited. All rights reserved.