Abstract
Methylammonium lead iodide perovskite (MAPbI_{3}) exhibits long charge carrier lifetimes that are linked to its high efficiency in solar cells. Yet, the mechanisms governing these unusual carrier dynamics are not completely understood. A leading hypothesis—disproved in this work—is that a large, static bulk Rashba effect slows down carrier recombination. Here, using second harmonic generation rotational anisotropy measurements on MAPbI_{3} crystals, we demonstrate that the bulk structure of tetragonal MAPbI_{3} is centrosymmetric with I4/mcm space group. Our calculations show that a significant Rashba splitting in the bandstructure requires a noncentrosymmetric lead iodide framework, and that incorrect structural relaxations are responsible for the previously predicted large Rashba effect. The small Rashba splitting allows us to compute effective masses in excellent agreement with experiment. Our findings rule out the presence of a large static Rashba effect in bulk MAPbI_{3}, and our measurements find no evidence of dynamic Rashba effects.
Introduction
Intense research efforts have enabled a dramatic increase in the power conversion efficiency of methylammonium lead iodide (MAPbI_{3})^{1} solar cells^{2,3,4,5}. These rapid advances relied on the many desirable properties of MAPbI_{3}, including its easy deposition from solution^{2,6} or vapor^{7}, a nearideal band gap^{8,9} that can be tuned with composition^{10,11}, and strong sunlight absorption in submicrometre thick films^{12}. Despite the rise in device performance, the origin of the exceptional transport properties of MAPbI_{3} remains the subject of debate. Long carrier diffusion lengths of up to 1 μm have been recognized early on as a key factor for device performance^{13}, and diffusion lengths as high as tens of microns in single crystals (and more recently in polycrystalline thin films)^{14,15,16} have been reported. The associated long carrier lifetimes of up to tens of microseconds are also particularly surprising given that MAPbI_{3} films are solution processed and polycrystalline^{14,17}.
Understanding carrier dynamics in MAPbI_{3} requires a close examination of its atomic and electronic structure, which are unusually complex compared to conventional solar cell materials^{18}. Similar to other perovskites, crystals of MAPbI_{3} exhibit multiple phases: cubic with a space group of \(Pm\bar 3m\) above 327 K^{19}, tetragonal with an elongated c axis and alternating inplane rotations of the PbI_{6} octahedra between 162 and 327 K^{19}, and orthorhombic with a space group of Pnma below 162 K^{20}. The rotations of the methylammonium (MA) ions^{19,21} and their interactions with the lead iodide (PbI) framework further complicate the structure. The space group of the room temperature tetragonal phase—the most relevant for practical applications—has been the subject of extensive debate, with some works reporting a centrosymmetric space group of I4/mcm (point group 4/mmm)^{19,20} and others a noncentrosymmetric space group of I4cm (point group 4mm)^{22,23}.
The presence of heavy elements such as lead and iodine introduces a strong spin–orbit coupling (SOC) in the electronic structure. When inversion symmetry is absent, as is the case in the bulk of noncentrosymmetric crystals and at surfaces or interfaces, the electrons feel an effective magnetic field due to the SOC. This interaction, known as the Rashba effect^{24,25}, removes the electron spin degeneracy and splits each of the valence and conduction band edges, resulting in an indirect band gap. Recent work proposed that a bulk Rashba effect—i.e., not due to the surface or interfaces—may be the cause of long carrier lifetimes and diffusion lengths in MAPbI_{3}^{26,27,28}. This idea has received surprising attention in spite of the fact that, with the exception of a few reports on the tetragonal phase^{22,23}, most diffraction experiments attribute centrosymmetric space groups to all three crystal phases^{19,20}.
There are two main mechanisms that can break inversion symmetry in the bulk of MAPbI_{3}, one static and one dynamic. The PbI framework can statically distort from a set of ideal cornersharing octahedra, forming a noncentrosymmetric phase as has been recently suggested^{23}, and resulting in a static bulk Rashba effect. Alternatively, inversion symmetry can be broken dynamically by the MA ion rotations in the absence of a static distortion of the PbI framework. Since the point group of the MA ions is C_{3v}, the relatively free MA ion rotation at room temperature^{19} implies that the highly symmetric point groups assigned to MAPbI_{3} may actually be approximations of a more complex dynamic structure^{19,20}. Dynamic structural fluctuations can also occur due to phonon modes or the interaction of the MA ions with the PbI framework^{28,29,30}. All these scenarios can induce a socalled dynamic bulk Rashba effect with a characteristic time scale shorter than those of the MA ion rotation (~10 ps)^{21} or the vibration of the highestfrequency PbI framework optical mode (~250 fs)^{31,32}. The strong interplay between the SOC, the PbI framework and the MA ion rotation intertwines the structural, electronic and vibrational degrees of freedom in MAPbI_{3}.
In this work, we set out to resolve the debate on the tetragonal space group and the related bulk Rashba effects by combining second harmonic generation rotational anisotropy (SHGRA) experiments, which have recently been established as a sensitive probe of point group symmetry^{33}, and density functional theory (DFT) calculations. We find that the SHGRA data fits excellently to the centrosymmetric I4/mcm space group, and perform several sets of DFT calculations that establish the absence of a large bulk Rashba effect in tetragonal MAPbI_{3}. When the PbI framework is constrained to the observed I4/mcm symmetry and only the MA ions are relaxed, we find a very small Rashba splitting due to the random MA ion orientation, making the band gap nearly direct and the band edge dispersions nearly parabolic. However, when two widely used approaches to relax tetragonal and cubic cells in DFT are employed, noncentrosymmetric crystal structures with distorted PbI framework and atomic positions inconsistent with experiment are formed. The large Rashba effect found in these structures is shown to be an artifact of DFT. The SHGRA measurements do not find evidence of a dynamic Rashba effect, and the computational results show that the MA ion rotations can generate only a small dynamic Rashba effect. Our results unravel complex structure–property relationships in tetragonal MAPbI_{3}, providing evidence of its centrosymmetric nature, direct band gap, and negligible bulk Rashba effect.
Results
Experimental determination of the point group of MAPbI_{3}
While it is difficult to distinguish between the centrosymmetric I4/mcm and noncentrosymmetric I4cm space groups with diffraction based techniques, the presence or absence of inversion symmetry has dramatic consequences on SHG. The bulk electric dipole (ED)induced SHG is only allowed in noncentrosymmetric crystals^{34}; in centrosymmetric crystals, the EDinduced SHG is forbidden in the bulk due to inversion symmetry, but still allowed at the surface where inversion symmetry is necessarily broken. Therefore it is necessary to determine whether an SHG signal has surface or bulk origin using a technique like SHGRA, which was not done in ref. ^{23}, before drawing conclusions about whether or not the bulk crystal structure breaks inversion symmetry. We perform SHGRA measurements on single crystals of MAPbI_{3} in the tetragonal phase at room temperature (see Methods) to determine its point group symmetry, and compare the SHGRA data with both the bulk ED response expected from a I4cm crystal and the surface ED response from a I4/mcm crystal, which has surface point group mm2. The space groups and associated point groups are summarized in Table 1.
Figure 1 shows the SHGRA results for different combinations of incident (in) and outgoing (out) polarizations, both of which can be either parallel (P) or perpendicular (S) to the light scattering plane, and compares the experimental data to simulated SHGRA patterns calculated based on the structure of the ED induced SHG susceptibility tensor for the mm2 and 4mm point groups^{33,35}.
For all polarization channels, the data fit excellently the angular dependence of the SHG intensity for the surface mm2 point group of a crystal with space group I4/mcm. On the other hand, the expected response from a crystal with space group I4cm is qualitatively incompatible with the measurements. In particular, the measured number and orientation of the lobes in the P_{in}−P_{out}, P_{in}−S_{out}, and S_{in}−P_{out} channels and the absence of any SHG intensity in the S_{in}−S_{out} channel cannot be reconciled with a 4mm point group. This result is consistent with several diffraction experiments^{19,20}, and directly disproves the recently hypothesized noncentrosymmetric space group I4cm^{22,23}. On this basis, given the centrosymmetric nature of MAPbI_{3}, the static bulk Rashba effect is expected to be negligible.
The resolution of our SHGRA measurements allows us to additionally comment on dynamic Rashba effects. If there were dynamic inversion symmetrybreaking due to the MA ions^{36} or the PbI framework phonons^{31,32}, then one would expect the SHG contribution from a 4mm point group to fluctuate between zero and some finite value, resulting in a nonzero timeaveraged signal. Given that a complete SHGRA scan takes several minutes, some admixture of the mm2 and 4mm contributions shown in Fig. 1 would be expected, with the latter likely being dominant since bulk ED signals are usually far stronger than surface signals. Given that we do not detect any 4mm contribution above our noise floor, we find no evidence of dynamic inversion symmetry breaking. With regards to the possibility of spatial domains, the SHGRA patterns shown in Fig. 1 were reproducible upon changing the beam position over the entire surface of the sample. Therefore we find no evidence of 4mm domains on the length scale of our resolution of ~10 microns.
Ab initio structural relaxation of MAPbI_{3}
The presence of bulk inversion symmetry in tetragonal MAPbI_{3} is crucial to correctly predict its properties using firstprinciples calculations. The computational determination of crystal structures is a difficult problem that is rich with subtleties. Starting from a guessed initial structure, an ideal computational approach would be able to relax the atomic positions of tetragonal MAPbI_{3} to our measured space group symmetry. DFT calculations have been used extensively to study the structure of MAPbI_{3}. However, since DFT is inherently a ground state (zerotemperature) method, it cannot be applied straightforwardly to obtain the structure of hightemperature phases (such as the tetragonal and cubic phases of MAPbI_{3}), which are stabilized by vibrational entropy. For example, the socalled vibrational soft modes in perovskite materials^{37}, which possess frequencies that vary continuously with temperature and exhibit peculiar trends at phase transitions, cannot be accurately described in DFT phonon calculations. This fact is a clear indication that probing the potential energy surface at zero temperature, as is done in DFT structural relaxations, is inadequate to compute the structures of the tetragonal and cubic MAPbI_{3} phases. While ab initio molecular dynamics based on DFT can capture thermal effects, the accurate firstprinciples methods available to compute bandstructures, absorption spectra and more recently carrier transport are applied on a snapshot of the structure, which clearly needs to be representative of the space group symmetry of the material. Despite these complexities and pitfalls, DFT structural relaxations of tetragonal and cubic MAPbI_{3} have been used extensively as a starting point to compute their bandstructures^{38,39,40,41}.
To study the impact of the structure on bandstructure and Rashba effect in tetragonal and cubic MAPbI_{3}, and understand the importance of using structure snapshots with the correct space group symmetry measured here, we carry out DFT structural relaxations using a variety of protocols, and then compute the bandstructure and bulk Rashba effect. We first relax initial tetragonal and cubic MAPbI_{3} structures without constraining the symmetry and lattice parameters during the relaxation, as is commonplace in the literature^{38,39,40,41}. These calculations employ \(\sqrt 2 \times \sqrt 2 \times 2\) cells with four formula units; each of the four MA ions is randomly aligned in one of the equivalent \(\left\langle {111} \right\rangle\) directions^{19}, and the localdensity approximation (LDA) of DFT is employed (see Methods).
Figure 2a, b shows the lattice parameters as a function of relaxation step for such cubic and tetragonal structures, respectively. In both phases, the three lattice parameters initially shrink by 1–2% compared to the initial (experimental) value, a wellknown trend for the LDA. As the relaxation continues, the three lattice parameters become distinct and both the initially tetragonal and cubic structures become orthorhombic, with three different lattice constants and nearly orthogonal lattice vectors. We interpret this result as a consequence of the lack of vibrational entropy in the calculation, and as an attempt by DFT to capture the zerotemperature orthorhombic phase of MAPbI_{3}. In these fully relaxed structures (see Fig. 2c, d) the Pb atoms move from the center of the octahedra, distorting the PbI framework and resulting in alternating long and short Pb−I bonds. This offcentering of the Pb atoms breaks inversion symmetry and generates a large structural distortion. While in previous work such orthorhombic structures have been described as pseudocubic^{26,41} and employed to investigate the electronic structure of tetragonal and cubic MAPbI_{3}, it is clear that the DFT ground state is entirely different from that of these hightemperature phases, which are stabilized by vibrational entropy. In a second set of DFT relaxations, the lattice parameters are fixed but the atoms are allowed to relax, as is also commonplace in the literature^{26,27,39}. The resulting structures (see Supplementary Fig. 1) are similar to the fully relaxed case, combining Pb offcentering, a distorted PbI framework, and extensive inversion symmetry breaking. Note that if the low temperature orthorhombic structure is fully relaxed, the octahedra remain undistorted ^{42}. We verify this result by fully relaxing an experimentally measured orthorhombic structure^{43}, and find that the structure remains centrosymmetric upon relaxation (see Supplementary Data 3 and the related bandstructure in Supplementary Fig. 2).
Bandstructure of cubic and tetragonal MAPbI_{3}
Next, we adopt an alternative procedure to obtain hightemperature MAPbI_{3} structures where the I4/mcm tetragonal and \(Pm\bar 3m\) cubic symmetries consistent with experiment are enforced by fixing the PbI framework and allowing only the MA ions to relax (see Methods). The resulting structures, referred to in the following as MArelaxed, are shown in Fig. 2e, f (data files of example MArelaxed structures are provided in Supplementary Data 1 and 2). The electronic bandstructures of the fully relaxed and MArelaxed structures are compared in Fig. 3 for both singleformulaunit cells and 4formulaunit cells with each of the 4 MA ions aligned randomly in one of the \(\left\langle {111} \right\rangle\) equivalent directions (see Methods). Note that the PbI framework is noncentrosymmetric in the fully relaxed structures and centrosymmetric in the MArelaxed structures, which are however overall noncentrosymmetric due to the random MA orientations. The bandstructures of the fully relaxed tetragonal and cubic structures exhibit a large Rashba splitting of the band edges. By contrast, in the MArelaxed structures the splitting is considerably reduced in all cases, and the residual Rashba effect is due to the parallel alignment (in the singleformulaunit cell in Fig. 3a, b) and random orientation (in the supercell in Fig. 3c) of the MA ions in our simulation cells. We expect that in the real crystal the net electric field from many randomly oriented polar MA ions, and thus the bulk Rashba effect, would be on average even smaller than that in our calculations on cells with only a handful of MA ions.
Given that the LDA plus SOC significantly underestimates the band gap, we recompute the bandstructure for the singleformulaunit MAPbI_{3} cubic structure using the HeydScuseriaErnzerhof (HSE) hybrid functional (see Methods)^{44}, which can improve the accuracy of the computed bandstructure. Figure 3b shows that the qualitative behavior of the Rashba splitting in the fully relaxed and MArelaxed structures is the same in the LDA and HSE calculations, a clear sign that the Rashba splitting is a consequence of the structural distortions of the PbI framework, regardless of the approximation employed to describe electron interactions. Employing centrosymmetric snapshots of the structure is clearly important since the magnitude of the Rashba splitting is fictitiously large in the fully relaxed structures with noncentrosymmetric PbI framework.
In the bandstructures in Fig. 3, the bulk Rashba effect removes the spin degeneracy of both the valence and conduction bands, leading to the formation of distinct parabolic band minima shifted from a highsymmetry point of the Brillouin zone (R point in the singleformulaunit cell and Γ point in the 4formulaunit supercells). This trend can be captured by a Rashba Hamiltonian (H_{R}) acting on each kpoint with an effective SOC interaction^{24}:
where σ are the Pauli spin matrices, and \(\widehat {\bf{n}}\) is a unit vector along the electric field induced by the breaking of inversion symmetry. The Rashba parameter α, which measures the magnitude of the splitting, depends on the crystallographic direction and on the magnitude of the macroscopic polarization due to inversion symmetry breaking^{45}. The resulting lowenergy bandstructure along a given direction consists of two split parabolic bands with energy \(E_ \pm ({k}) = \frac{{\hbar ^2{k}^2}}{{2m^ \ast }} \pm \alpha {k}\) and, for the conduction band, with minima at k_{ m } = α × m^{*}/ħ^{2} (where m^{*} is the effective mass), from which the Rashba parameter can be extracted using \(\alpha = 2{\kern 1pt} \left {E_  \left( {{k}_m} \right)} \right{\mathrm{/}}{k}_m\). Analogous expressions hold for the valence band and its Rashba parameter. A system that is only slightly perturbed from its centrosymmetric structure—such as our MArelaxed structures—exhibits a small value of α, and accordingly a small band splitting^{24}.
Quantifying the Rashba effect
To quantify the bulk Rashba splitting in MAPbI_{3}, we first compute the bandstructures along high symmetry directions for 10 tetragonal and 10 cubic 4formulaunit MArelaxed supercells. In each structure, the 4 MA ions are initially randomized in one of the 8 equivalent \(\left\langle {111} \right\rangle\) directions and then relaxed, thus allowing us to extract a distribution of Rashba parameters and effective masses over multiple sets of MA ion orientations^{19}. In each of these MArelaxed structures, the PbI framework is fixed to the experimental symmetry and only the MA ions are allowed to relax. The average and standard deviation of the Rashba parameters for these structures are given in Table 2, and an example bandstructure is given in Fig. 3c.
Our computed Rashba parameters are of order (0.1–0.2) ± 0.1 eV Å in the cubic structures and (0.2–0.4) ± 0.2 eV Å in the tetragonal structures. These values are an order of magnitude smaller than those of fully relaxed structures with large distortions of the PbI framework. For example, ref. ^{41} reported α ≈ 1 eV Å in fully relaxed cubic structures, and ref. ^{25} a value of α = 4 eV Å in a distorted cubic structure with P4mm space group. While refs. ^{26,40} did not compute the Rashba coefficient of their fully relaxed structures, a large Rashba splitting can be seen in their bandstructures. We conclude that once the experimental crystal structure is maintained, the bulk Rashba effect is much smaller than previously predicted and only due to the random MA ion orientations, leading to a nearly direct band gap for tetragonal and cubic MAPbI_{3} and to an expected negligible effect of SOC on charge transport. In fact, since carriers undergo diffusive transport with scattering events spaced apart by tens of unit cells, they feel a space and time averaged structure, and are insensitive to the small Rashba effect induced by the local MA ion orientations. Lastly, note that the bandstructure of orthorhombic MAPbI_{3} similarly exhibits a negligible Rashba effect (see Supplementary Fig. 2).
Discussion
The presence of a surface Rashba effect is clearly not ruled out by our work, given that inversion symmetry is inherently broken at surfaces and interfaces; its role in devices may be significant^{46}. Since inversion symmetry is also broken at grain boundaries, Rashba effects in polycrystalline films may also differ greatly from our single crystal results. The influence of Rashba effects due to surfaces, interfaces and grain boundaries depends on the penetration depth of the surface electric field that couples to the electrons through the SOC. The Debye length (L_{D}) provides an estimate of the length scale over which such surface Rashba effects are significant:
where \(\epsilon _0\) and \(\epsilon _\infty\) are the dielectric constant of the vacuum and of the material, respectively, k_{B} the Boltzmann constant, T is temperature, e the electron charge and n the carrier concentration. Using the reported room temperature high frequency dielectric constant of 23.3^{47}, and an intrinsic charge carrier density of 1 × 10^{9} cm^{−3} ^{22}, we find a Debye length value of ~180 μm in the dark, and ~180 nm under solar illumination assuming the carrier concentration reaches ~10^{15} cm^{−3} ^{48}. Surface Rashba effects could thus be particularly relevant in solar cells, where perovskite films are typically several hundred nanometers thick^{49,50}, and in polycrystalline thin films, where typical grain sizes in high performing devices are of order 100 nm–1 μm^{49,50}. Our results imply that the source of any Rashba effects that may have been observed experimentally should be sought at such surfaces and interfaces rather than in the bulk of MAPbI_{3}.
In fact, it is worth reconciling our results with recent reports of giant Rashba splitting in MAPbBr_{3} single crystals using angleresolved photoelectron spectroscopy (ARPES)^{51} and Rashba splitting in CsPbBr_{3} nanocrystals^{52}. ARPES is well known to probe electronic states within 1−10 nm of the surface, where the surface Rashba effect dominates due to the long Debye screening length, as noted above. We thus argue that ref. ^{51} accurately measures the surface Rashba effect, and the large Rashba splitting observed is not indicative of the bulk bandstructure of the material. Similarly, since the nanocrystals measured in ref. ^{52} are ~5 nm in diameter, their electronic and optical properties are expected to be dominated by surface Rashba effects.
We additionally note that previous work extracted effective masses through parabolic fits of bandstructures with strong Rashba splitting^{26,40}. However, in the presence of a strong Rashba effect transport cannot be described by a free carrier approximation due to the presence of a spin texture^{53}, and extracting effective masses is an illdefined procedure due to the nonparabolic dispersions. Since our MArelaxed structures exhibit a very small Rashba splitting, we can obtain effective masses for a range of MA ion orientations. The electron and hole effective masses of the 10 MArelaxed tetragonal and cubic structures discussed above are shown in Fig. 4 The computed effective masses are in excellent agreement with recent experiments^{54}. They change appreciably for different MA ion random orientations in spite of the fact that the PbI framework is identical in all the structures, thus justifying the spread in computed effective mass values reported in the literature. While the cubic phase exhibits electron masses lighter than hole masses (by ~50%) as is common in semiconductors, we find an unconventional result in the tetragonal phase, where the electron and hole masses are nearly identical in each of the random MA orientations. These results appear to support a hypothesis advanced recently of balanced transport between the electrons and holes in tetragonal MAPbI_{3}^{55,56}, though accurate calculations of the mobility are necessary to further ascertain this aspect.
Finally, it is interesting to compare the carrier dynamics in MAPbI_{3} with perhaps the bestperforming conventional solar cell semiconductor, GaAs. The effective masses found here for tetragonal MAPbI_{3} are comparable with those of electrons and light holes in GaAs^{57}—the electron and light hole masses in GaAs are slightly smaller than those in MAPbI_{3}, but the heavy hole masses in GaAs are considerably larger^{57}. Commonly reported values for the nonradiative (τ_{nr} = 100 ns) and radiative (τ_{rad} = 2 μs) carrier lifetimes in MAPbI_{3}^{48}, combined with the effective masses, give estimated diffusion lengths in the range of 1–10 μm, similar to values reported experimentally^{13,14,15,16,17} and comparable with those of moderately doped GaAs^{58}. Finally, both GaAs and tetragonal MAPbI_{3} luminesce strongly and with a high quantum yield in surface passivated samples, clearly suggesting that the band gap of MAPbI_{3} is direct^{14}. The parallel between tetragonal MAPbI_{3} and GaAs, two materials defining a gold standard for solar cells, is intriguing and worthy of further investigation.
In summary, our measurements establish that tetragonal MAPbI_{3} at room temperature is centrosymmetric with space group I4/mcm. Our calculations show that the bulk Rashba effect is very small in MAPbI_{3} structures with symmetry consistent with experiment, and the band gap nearly direct. Greatly overlooked structural artifacts due to DFT relaxation result in a fictitiously large bulk Rashba effect due to structural distortions of the PbI framework, both in tetragonal and cubic MAPbI_{3}. The DFT fully relaxed structures are inconsistent with the centrosymmetric point groups of these phases, calling for a critical revision of existing bandstructure calculations. While the Rashba effect^{26,27} has been invoked in order to explain the favorable transport properties of tetragonal MAPbI_{3}, we show that SOC effects are negligible, and argue that the favorable transport in MAPbI_{3} is mainly the result of its intrinsic small effective masses and long carrier lifetimes, similar to GaAs.
Methods
DFT calculations
We carry out DFT calculations within the local density approximation (LDA)^{59} using the QUANTUM ESPRESSO code^{60}. The calculations employ fully relativistic ultrasoft pseudopotentials^{61} that include the SOC, a kinetic energy cutoff of 60 Ry for the planewave basis set and a uniform 6 × 6 × 6 kpoint grid. Calculations are performed on both the cubic and tetragonal phases of MAPbI_{3}; both singleformulaunit structures and \(\sqrt 2 \times \sqrt 2 \times 2\) supercells with 4 formula units are employed. In the supercells, the initial structures are prepared using experimental lattice parameters, and each of the 4 methylammonium (MA) ions is oriented randomly in one of the 8 equivalent \(\left\langle {111} \right\rangle\) directions^{19}. The structures are relaxed with and without symmetry constraints within the LDA until the force on each atom is less than 0.1 eV/Å. For calculations where only the MA ions are relaxed and the PbI framework is fixed, the lattice parameters are chosen as the ones minimizing the DFT forces for the given symmetry. In these calculations, the PbI framework is fixed and the MA ions are allowed to relax until the average force on each atom less than 0.1 eV/Å. Since LDA plus SOC is known to drastically underestimate the band gap, for bandstructure calculations employed within the LDA an empirical scissor shift is applied in Figs. 3 and 4 to match the experimental band gap. The effective masses of the electrons and holes are extracted using fine grids near the conduction band minimum and valence band maximum, respectively. The DFT calculations using the HSE hybrid exchangecorrelation functional are performed using the Vienna Ab initio Simulation Program (VASP) code^{62} and projector augmented wave (PAW)^{63} pseudopotentials. A kinetic energy cutoff of 400 eV and uniform 6 × 6 × 6 kpoint grids are employed, together with the SOC, as implemented in the VASP code.
Crystal growth
MAPbI_{3} single crystals were grown by the inverse temperature crystallization method as previously reported^{64}. Briefly, 1 M MAPbI_{3} solution was first prepared by dissolving equimolar MAI and PbI_{2} in gammabutyrolactone (GBL) at 70 °C. The temperature of the asprepared solution was then stepwise increased from 70 to 130 °C with a rate of 5 °C/h to crystallize MAPbI_{3} by utilizing the inverse solubility of MAPbI_{3} in GBL. Asgrown MAPbI_{3} crystals were collected from the hot solution with filter papers to remove remaining solution on the surface. The MAPbI_{3} crystals were further polished using a Unipo1210 polisher with SiC papers (grit size of 600–200) to produce smooth surfaces for the SHG study. The photograph of a sample is given in the Supplementary Fig. 3.
SHGRA measurements
SHGRA measurements are performed on the polished (100) surface of MAPbI_{3} at room temperature using a fast rotating scattering plane based technique (see ref. ^{33}). The incident light is derived from a regeneratively amplified Ti:sapphire laser system operating at a 100 kHz repetition rate with a center wavelength of 800 nm and a pulse duration of 100 fs. The incident fluence is maintained at below 1.5 mJ/cm^{2} and samples are never heated above the structural phase transition temperature of 327 K to avoid sample damage.
Expressions for SHGRA patterns
The mathematical expressions for the electricdipole induced SHGRA patterns for a surface point group mm2 are listed below for different polarization geometries:
where ϕ and θ are the scattering plane angle and angle of laser incidence, respectively. Using the experimental value of θ = 10°, a best fit of Eqs. 3, 4, and 5 to the data in Fig. 1 yields the following parameters:
The mathematical expressions for the electricdipole induced SHGRA patterns for a bulk point group 4mm along the [100] direction are listed below for different polarization geometries:
Data availability
The authors declare that the data supporting the findings of this study are available within the paper and its Supplementary Information files.
References
 1.
Kojima, A., Teshima, K., Shirai, Y. & Miyasaka, T. Organometal halide perovskites as visiblelight sensitizers for photovoltaic cells. J. Am. Chem. Soc. 131, 6050–6051 (2009).
 2.
Kim, H.S. et al. Lead iodide perovskite sensitized allsolidstate submicron thin film mesoscopic solar cell with efficiency exceeding 9%. Sci. Rep. 2, 591 (2012).
 3.
McGehee, M. D. Perovskite solar cells: continuing to soar. Nat. Mater. 13, 845–846 (2014).
 4.
Yang, W. S. et al. Iodide management in formamidiniumleadhalide–based perovskite layers for efficient solar cells. Science 356, 1376–1379 (2017).
 5.
National Renewable Energy Lab. Research Cell Record Efficiency Chart (2018). URL https://www.nrel.gov/pv/assets/images/efficiencychart.png. Accessed 19 March 2018.
 6.
Burschka, J. et al. Sequential deposition as a route to highperformance perovskitesensitized solar cells. Nature 499, 316–320 (2013).
 7.
Liu, M., Johnston, M. B. & Snaith, H. J. Efficient planar heterojunction perovskite solar cells by vapour deposition. Nature 501, 395–398 (2013).
 8.
Quarti, C. et al. Structural and optical properties of methylammonium lead iodide across the tetragonal to cubic phase transition: implications for perovskite solar cells. Energy Environ. Sci. 16, 155–163 (2016).
 9.
Rühle, S. Tabulated values of the ShockleyQueisser limit for single junction solar cells. Sol. Energy 130, 139–147 (2016).
 10.
Eperon, G. E. et al. Formamidinium lead trihalide: a broadly tunable perovskite for efficient planar heterojunction solar cells. Energy Environ. Sci. 7, 982–988 (2014).
 11.
Eperon, G. E. et al. Perovskiteperovskite tandem photovoltaics with optimized band gaps. Science 354, 861–865 (2016).
 12.
Ball, J. M., Lee, M. M., Hey, A. & Snaith, H. J. Lowtemperature processed mesosuperstructured to thinfilm perovskite solar cells. Energy Environ. Sci. 6, 1739 (2013).
 13.
Stranks, S. D. et al. Electronhole diffusion lengths exceeding 1 micrometer in an organometal trihalide perovskite absorber. Science 342, 341–344 (2013).
 14.
Brenes, R. et al. Metal halide perovskite polycrystalline films exhibiting properties of single crystals. Joule 1, 155–167 (2017).
 15.
Alarousu, E. et al. Ultralong radiative states in hybrid perovskite crystals: compositions for submillimeter diffusion lengths. J. Phys. Chem. Lett. 8, 4386–4390 (2017).
 16.
Shi, D. et al. Low trapstate density and long carrier diffusion in organolead trihalide perovskite single crystals. Science 347, 519–522 (2015).
 17.
DeQuilettes, D. W. et al. Photoluminescence lifetimes exceeding 8 μs and quantum yields exceeding 30% in hybrid perovskite thin films by ligand passivation. ACS Energy Lett. 1, 438–444 (2016).
 18.
Bernardi, M. & Grossman, J. C. Computer calculations across time and length scales in photovoltaic solar cells. Energy Environ. Sci. 9, 2197–2218 (2016).
 19.
Kawamura, Y., Mashiyama, H. & Hasebe, K. Structural study on cubictetragonal transition of CH_{3}NH_{3}PbI_{3}. J. Phys. Soc. Jpn 71, 1694–1697 (2002).
 20.
Weller, M. T., Weber, O. J., Henry, P. F., Di Pumpo, A. M. & Hansen, T. C. Complete structure and cation orientation in the perovskite photovoltaic methylammonium lead iodide between 100 and 352 K. Chem. Commun. 51, 4180–4183 (2015).
 21.
Leguy, A. M. et al. The dynamics of methylammonium ions in hybrid organic–inorganic perovskite solar cells. Nat. Commun. 6, 7124 (2015).
 22.
Stoumpos, C. C., Malliakas, C. D. & Kanatzidis, M. G. Semiconducting tin and lead iodide perovskites with organic cations: phase transitions, high mobilities, and nearinfrared photoluminescent properties. Inorg. Chem. 52, 9019–9038 (2013).
 23.
Rakita, Y. et al. Tetragonal CH_{3}NH_{3}PbI_{3} is ferroelectric. Proc. Natl Acad Sci. USA 114, 5504–5512 (2017).
 24.
Bychkov, Y. A. & Rashba, E. Properties of a 2D electron gas with spectral degeneracies lifted. JETP Lett. 39, 78–81 (1984).
 25.
Kepenekian, M. et al. Rashba and Dresselhaus effects in hybrid organicinorganic perovskites: from basics to devices. ACS Nano 9, 15557–11567 (2015).
 26.
Azarhoosh, P., McKechnie, S., Frost, J. M., Walsh, A. & van Schilfgaarde, M. Relativistic origin of slow electronhole recombination in hybrid halide perovskite solar cells. APL Mater. 4, 091501 (2016).
 27.
Zheng, F. et al. Rashba spin–orbit coupling enhanced carrier lifetime in CH_{3}NH_{3}PbI_{3}. Nano Lett. 15, 7794–7800 (2015).
 28.
Etienne, T., Mosconi, E. & De Angelis, F. Dynamical origin of the Rashba effect in organohalide lead perovskites: a key to suppressed carrier recombination in perovskite solar cells? J. Phys. Chem. Lett. 7, 1638–1645 (2016).
 29.
Motta, C. et al. Revealing the role of organic cations in hybrid halide perovskite CH_{3}NH_{3}PbI_{3}. Nat. Commun. 6, 7026 (2015).
 30.
Wu, X. et al. Lightinduced picosecond rotational disordering of the inorganic sublattice in hybrid perovskites. Sci. Adv. 3, e160238 (2017).
 31.
Leguy, A. M. et al. Dynamic disorder, phonon lifetimes, and the assignment of modes to the vibrational spectra of methylammonium lead halide perovskites. Phys. Chem. Chem. Phys. 18, 27051–27066 (2016).
 32.
Sendner, M. et al. Optical phonons in methylammonium lead halide perovskites and implications for charge transport. Mater. Horiz. 3, 613–620 (2016).
 33.
Harter, J., Zhao, Z., Yan, J.Q., Mandrus, D. & Hsieh, D. A paritybreaking electronic nematic phase transition in the spinorbit coupled metal Cd_{2}Re_{2}O_{7}. Science 356, 295–299 (2017).
 34.
GuyotSionnest, P. & Shen, Y. R. Bulk contribution in surface secondharmonic generation. Phys. Rev. B 38, 7985–7989 (1988).
 35.
Boyd, R. W. Nonlinear Optics (Academic Press, San Diego, CA, USA, 2003).
 36.
Yaffe, O. et al. Local polar fluctuations in lead halide perovskite crystals. Phys. Rev. Lett. 118, 136001 (2017).
 37.
Glazer, A. The classification of tilted octahedra in perovskites. Acta Cryst. B 28, 3384–3392 (1972).
 38.
Agiorgousis, M. L., Sun, Y.Y., Zeng, H. & Zhang, S. Strong covalencyinduced recombination centers in perovskite solar cell Material CH_{3}NH_{3}PbI_{3}. J. Am. Chem. Soc. 136, 14570–14575 (2014).
 39.
Buin, A. et al. Materials processing routes to trapfree halide perovskites. Nano Lett. 14, 6281–6286 (2014).
 40.
Umari, P., Mosconi, E. & De Angelis, F. Relativistic GW calculations on CH_{3}NH_{3}PbI_{3} and CH_{3}NH_{3}SnI_{3} perovskites for solar cell applications. Sci. Rep. 4, 4467 (2014).
 41.
He, Y. & Galli, G. Perovskites for solar thermoelectric applications: a first principle study of CH_{3}NH_{3}AI_{3} (A=Pb and Sn). Chem. Mater. 26, 5394–5400 (2014).
 42.
Filip, M. R., Verdi, C. & Giustino, F. GW band structures and carrier effective masses of CH_{3}NH_{3}PbI_{3} and hypothetical perovskites of the type APbI_{3}: A=NH_{4}, PH_{4}, AsH_{4}, and SbH_{4}. J. Phys. Chem. C 119, 25209–25219 (2015).
 43.
Whitfield, P. S. et al. Structures, phase transitions and tricritical behavior of the hybrid perovskite methyl ammonium lead iodide. Sci. Rep. 6, 35685 (2016).
 44.
Heyd, J., Scuseria, G. E. & Ernzerhof, M. Hybrid functionals based on a screened Coulomb potential. J. Chem. Phys. 118, 8207–8215 (2003).
 45.
Leppert, L., ReyesLillo, S. E. & Neaton, J. B. Electric field and straininduced Rashba effect in hybrid halide perovskites. J. Phys. Chem. Lett. 7, 3683–3689 (2016).
 46.
Mosconi, E., Etienne, T. & De Angelis, F. Rashba band splitting in organohalide lead perovskites: bulk and surface effects. J. Phys. Chem. Lett. 8, 2247–2252 (2017).
 47.
OnodaYamamuro, N., Matsuo, T. & Suga, H. Dielectric study of CH_{3}NH_{3}PbX_{3} (X=Cl, Br, I). J. Phys. Chem. Solids 53, 935–939 (1992).
 48.
Johnston, M. B. & Herz, L. M. Hybrid perovskites for photovoltaics: chargecarrier recombination, diffusion, and radiative efficiencies. Acc. Chem. Res. 49, 146–154 (2016).
 49.
Saliba, M. et al. Incorporation of rubidium cations into perovskite solar cells improves photovoltaic performance. Science 354, 206–209 (2016).
 50.
Zhang, W. et al. Enhanced optoelectronic quality of perovskite thin films with hypophosphorous acid for planar heterojunction solar cells. Nat. Commun. 6, 10030 (2015).
 51.
Niesner, D. et al. Giant Rashba splitting in CH_{3}NH_{3}PbBr_{3} organicinorganic perovskite. Phys. Rev. Lett. 117, 126401 (2016).
 52.
Isarov, M. et al. Rashba effect in a single colloidal CsPbBr_{3} perovskite nanocrystal detected by magnetooptical measurements. Nano Lett. 17, 5020–5026 (2017).
 53.
Molenkamp, L. W., Schmidt, G. & Bauer, G. E. W. Rashba Hamiltonian and electron transport. Phys. Rev. B 64, 121202 (2001).
 54.
Miyata, A. et al. Direct measurement of the exciton binding energy and effective masses for charge carriers in organicinorganic trihalide perovskites. Nat. Phys. 11, 582–587 (2015).
 55.
Xing, G. et al. Longrange balanced electron and holetransport lengths in organicinorganic CH_{3}NH_{3}PbI_{3}. Science 342, 344–347 (2013).
 56.
Giorgi, G., Fujisawa, J.I., Segawa, H. & Yamashita, K. Small photocarrier effective masses featuring ambipolar transport in methylammonium lead iodide perovskite: a density functional analysis. J. Phys. Chem. Lett. 4, 4213–4216 (2013).
 57.
Yu, P. Y. & Cardona, M. Fundamentals of Semiconductors: Physics and Materials Properties (Springer, New York, USA, 2010).
 58.
Aukerman, L. W., Millea, M. F. & McColl, M. Diffusion lengths of electrons and holes in GaAs. J. Appl. Phys. 38, 685–690 (1967).
 59.
Perdew, J. P. & Zunger, A. Selfinteraction correction to densityfunctional approximations for manyelectron systems. Phys. Rev. B 23, 5048–5079 (1981).
 60.
Giannozzi, P. et al. QUANTUM ESPRESSO: a modular and opensource software project for quantum simulations of materials. J. Phys. Condens. Matter 21, 395502 (2009).
 61.
Vanderbilt, D. Soft selfconsistent pseudopotentials in a generalized eigenvalue formalism. Phys. Rev. B 41, 7892–7895 (1990).
 62.
Kresse, G. & Furthmüller, J. Efficient iterative schemes for ab initio totalenergy calculations using a planewave basis set. Phys. Rev. B 54, 11169–11186 (1996).
 63.
Blöchl, P. E. Projector augmentedwave method. Phys. Rev. B 50, 17953–17979 (1994).
 64.
Saidaminov, M. I. et al. Highquality bulk hybrid perovskite single crystals within minutes by inverse temperature crystallization. Nat. Commun. 6, 7586 (2015).
 65.
Setyawan, W. & Curtarolo, S. Highthroughput electronic band structure calculations: challenges and tools. Comput. Mater. Sci. 49, 299–312 (2010).
Acknowledgements
M.B. acknowledges partial support from startup funds and from the Space Solar Program Initiative at the California Institute of Technology. K.F. thanks the California Institute of Technology for support through the SURF fellowship program, and Peter Foley and Linn Leppert for fruitful discussions. SHGRA measurements were supported by the U. S. Department of Energy under grant DESC0010533. D.H. also acknowledges funding for instrumentation from the Institute for Quantum Information and Matter, an NSF Physics Frontiers Center (PHY1125565) with support of the Gordon and Betty Moore Foundation through grant GBMF1250. O.M.B and W.P. acknowledge the support of KAUST. J.B.N. and S.G.L. were supported by the U.S. Department of Energy, Director, Office of Science, Office of Basic Energy Sciences, Materials Sciences and Engineering Division, under Contract No. DEAC0205CH11231, through the Theory FWP (KC2301) at Lawrence Berkeley National Laboratory (LBNL). This research used resources of the National Energy Research Scientific Computing Center, a DOE Office of Science User Facility supported by the Office of Science of the U.S. Department of Energy under Contract No. DEAC02 05CH11231.
Author information
Affiliations
Contributions
M.B. conceived the project. K.F. and M.B. carried out the calculations and prepared the manuscript with input from coauthors. W.P. and O.B. synthesized and provided the samples. T.D., J.H., and D.H. designed and conducted the experiments. B.B. contributed to the calculations. J.B.N. and S.G.L. contributed to conceive the theory and calculations. All authors edited the manuscript.
Corresponding authors
Ethics declarations
Competing interests
The authors declare no competing interests.
Additional information
Publisher's note: Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons license, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons license and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/.
About this article
Cite this article
Frohna, K., Deshpande, T., Harter, J. et al. Inversion symmetry and bulk Rashba effect in methylammonium lead iodide perovskite single crystals. Nat Commun 9, 1829 (2018). https://doi.org/10.1038/s4146701804212w
Received:
Accepted:
Published:
Further reading

Electronic and optical properties of vacancy ordered double perovskites A2BX6 (A = Rb, Cs; B = Sn, Pd, Pt; and X = Cl, Br, I): a first principles study
Scientific Reports (2021)

Minimizing hydrogen vacancies to enable highly efficient hybrid perovskites
Nature Materials (2021)

Lowfield onset of WannierStark localization in a polycrystalline hybrid organic inorganic perovskite
Nature Communications (2021)

Implications of the size variation on the local structure and polarized emission of CsPbBr3 quantum dots
Journal of Materials Science (2021)

Enhanced stability in CH3NH3PbI3 hybrid perovskite from mechanochemical synthesis: structural, microstructural and optoelectronic characterization
Scientific Reports (2020)
Comments
By submitting a comment you agree to abide by our Terms and Community Guidelines. If you find something abusive or that does not comply with our terms or guidelines please flag it as inappropriate.