Arxiv:1612.04781V1

Home , Lattice energy

arXiv:1612.04781v1 [cond-mat.mtrl-sci] 14 Dec 2016 ulydcmoeatrapro ftm htrne rmsev- from ranges month a that to time days of eral period Th a ev after perovskites decompose either. organic tually hybrid degradation encapsulated their of prevent layer active not does cells ovskite thermodynamic oeua yai simulations dynamic molecular n aieeeet,rsetvl.Tengtv ecine reaction negative The thalpy respectively. involved. elements, compounds halide solid and the of energy Here total the on based yeautn h nhlyo h reaction the of calculations enthalpy first-principle the of evaluating level by the at captured be can of iin o xml,terato nhlyfrtetragonal for enthalpy reaction the example, CH For sition. vktswspeiul icse rmai-aechemistr per acid-base halide from for discussed pathway previously degradation was ovskites a creating in ra- takes moisture etc.) UV of that oxygen, (moisture, separation atmospheric conditions phase ambient diation, a under to instantly prone place materia active are per- the perovskites of organic instability brid the hybrid is of cells obsta solar commercialization ovskite main the The hindering devices. temper- cle alternativ low photovoltaic low-cost a thin-film a on offering crystalline thus based 20% method, is over solution perovksite ature to organic hybrid increased of already tion have cells lar uaunit mula r ewe ekysal n ntbearen ihexper- with agreeing observations unstable imental and stable weakly between ary factor an instability unclear. remain intrinsic improvement pe its of halide for origin enues hybrid the of structures, formability ovskite predicting in calculations nrni ntblt fhbi aieprvkt structu perovskite halide hybrid of instability Intrinsic h fcece fhbi rai-nrai eosieso perovskite organic-inorganic hybrid of efficiencies The emtia atr uha h odcmd’ tolerance Goldschmidt’s the as such factors Geometrical ∆ 3 NH H A 18 ∆ r h oesal h tutr saantdecompo- against is structure the stable more the , 3 ersnsa rai cation, organic an represents n caerlfco ucsflyepanformabil- explain successfully factor octahedral and PbI H 7,12,15 r 3 niae tbepout.Telwrtevalue the lower The products. stable indicates fhbi aieprvkts rset ffidn e perov discussed. new are finding applications cr of photovoltaic additional for Prospects an aimed as perovskites. cage halide a hybrid and of molecule a of vario energy of ionization analysis vari The of components. enthalpy energy reaction molecularization the formabil partition capture to cycle fully Born-Haber not perovsk do halide factors of geometrical stability tahedral predicting obse for experimentally criteria evalua to simple linked for is used which are perovskites, halide calculations First-principle cations. swti h ag of range the within is hc edr h tutr ob ttebound- the at be to structure the renders which , 12 ntblt fhbi rai-nrai aieperovskit halide organic-inorganic hybrid of Instability esetvs naslto fteper- the of Encapsulation perspectives. AX .INTRODUCTION I. 13,14 7 ept h ucs ffirst-principle of success the Despite . + oiaineeg sasaiiycieinfrhld perov halide for criterion stability a as energy Ionization . BX 1 eateto aeil cec n niern,McMaster Engineering, and Science Materials of Department 2 9,10 → yrlssreaction hydrolysis , 20Mi tetWs,Hmlo,OtroLS48 Canada 4L8, L8S Ontario Hamilton, West, Street Main 1280 5–7 ABX − h ermna role detrimental The . B 0 . 1 and 3 . . . 0 X . hoZheng Chao 06 1–4 r h metal the are Fabrica- . Vprfor- per eV Dtd eebr1,2016) 15, December (Dated: 11 .Hy- l. av- d 15–17 to e and en- (1) res y 1, n- r- 8 e - - - , ∗ n lgRubel Oleg and tsaelcigsneGlshitstlrneadoc- and tolerance Goldschmidt’s since lacking are ites ino eopsto ecinetap fhybrid of enthalpy reaction decomposition a of tion t ftoeprvkts nti ae,w xedthe extend we paper, this In perovskites. those of ity ue h au ftetlrnefco o CH for factor tolerance the of value The ture. nefactor ance t fvrosiogncprvkt structures perovskite inorganic various of ity h est ucinlter DT.I ilb hw htin that shown be will It (DFT). I theory functional structures density perovskite the various of for components enthalpies energy reaction of the analysis the to cycle Born-Haber the vkt opud.Frisac,taiinliogncpe inorganic per- II traditional non-halide the of instance, to ovskites For compared as compounds. structure ovskite latti ionic electrostatic their low of relatively ergy a to perovskites halide h ntblt fhbi aieperovskites. expl halide to hybrid sufficient of not instability are the factors geometrical that suggests nrylaigteinzto nhlyt eemn h dir the determine to enthalpy ionization moleculariza the the leaving by energy cancelled largely is contribution ergy inetap,adltieenthalpy lattice and ioni enthalpy, enthalpy, tion atomization ener formation components: the following break the to into allows It enthalpies. formation stabil for criterion a as structures. ionic used of be cannot alone energy lattice of hi energy its of lattice despite decomposition, without melts material the t aeo CH of case t eog oteI the to belongs nry ti lokonta CsPbI that known also is It energy. X atc nryof energy lattice of energy tice − tutrsadpooe oadPuigsothda factor octahedral perovkit Pauling’s of add r formability to capture proposed fully and structures not does alone factor ad h xeietlrato nhlyfrPbTiO for enthalpy oth reaction the experimental On the stability. chemical hand, weaker thus and energy static htI that nedsal pt h eprtr f460 of temperature the to up stable indeed ihntealwberneof range allowable the within shnesterdvlpetfrpoooti appli- photovoltaic for development their hinders es vddgaaino eiecaatrsis However, characteristics. device of degradation rved u eosiesrcue noltie oiain and ionization, lattice, into structures perovskite ous − scnrbtost h ecinetap onst an to points enthalpy reaction the to contributions us B 0 0 = 0 = epciey sa diinlgoerclciein In criterion. geometrical additional an as respectively) , rs tal. et Frost h onHbrccei rdtoal sdfraayi of analysis for used traditionally is cycle Born-Haber The II /r . kt tutrswt mrvdceia stability chemical improved with structures skite tro o rdcigceia rnsi stability in trends chemical predicting for iterion 38 − X − . VII . 8 eV 91 II ( r − − 20 B 3 22 VII 3 0 hc swti h ag facpal values acceptable of range the within is which , rai n nrai eosie h atc en- lattice the perovskites inorganic and organic hc sodr fmgiuels hnislattice its than less magnitude of orders is which , t and NH . 95 esrscmateso h eosiestruc- perovskite the of compactness measures 1 3 8 − 3 19 eosie aeitisclylwrelectro- lower intrinsically have perovskites r − − trbtdteisaiiyo yrdorganic hybrid of instability the attributed − − PbI 119 X ie al. et Li . 27 II IV 28 − r h oi ai fcation of radii ionic the are 3 − V hs bevtosidct htthe that indicate observations These eV. V hsvlei uhlwrta the that lower much is value This eV. VfrCH for eV VII University, h caerlfactor octahedral the VI 3 3 aiy ..PbTiO e.g. family, aiy hsagmn suggests argument This family. 19 skites 0 one u httetolerance the that out pointed . 414 3 NH 24 nti ae eextend we paper this In . − 3 3 PbI eosiesrcueis structure perovskite 0 . 732 3 ◦ eokie which perovksite, C r 21 Pb 23 hsanalysis This . 19 3 /r bv which above , aetelat- the have , h toler- The . 3 B I NH 0 = n anion and 3 3 PbI sonly is een- ce . using 54 gher tion 3 the ain za- ec- ity gy er 21 is is r- e 2 tion of the reaction. The instability of hybrid organic lead- III. COMPUTATIONAL DETAILS iodine perovskites can be attributed to the high energy associ- – ated with ionization of organic molecules and [PbI3] . Electronic structure calculations have been performed in the framework of DFT27 and Perdew-Burke-Ernzerhof generalized gradient approximation28 (GGA-PBE) for the II. BASIC CONCEPTS exchange-correlation functional. Total energies of all compounds were obtained using the Vienna ab initio simula- The Born-Haber cycle was originally proposed by Max tion program (VASP) and projector augmented-wave (PAW) Born and Fritz Haber as a way to measure formation ener- potentials29–31. gies of ionic structures24. The cycle also provides a method All crystal structures of compounds studied here are to determine the lattice energy of the structures, which oth- taken at their most stable polymorph at ambient conditions. erwise cannot be directly measured experimentally. Here we Among perovskite structures, CH3NH3PbI3 adapts a tetrag- will explain the essence of the Born-Haber cycle and its uti- onal β-phase at the ambient temperature, CH3NH3PbBr3 32,33 lization for analysis of reaction enthalpy components using and CH3NH3PbCl3 have a cubic phase , CN3H6PbI3 34–36 the CH3NH3PbI3 perovksite structure as an example. and (CH3)4NPbI3 favor hexagonal structures . CsPbI3, The formation process of CH3NH3PbI3 from solid CsPbBr3, and CsPbCl3 prefer an orthorhombic (Pnma) δ- 37,38 CH NH I and PbI compounds can be subdivided into sev- phase . The crystal structure of CH3NH3I, CH3NH3Br, 3 3 2 ′ eral consecutive steps illustrated in Fig. 1. and CH3NH3Cl organic salts correspond to α -tetragonal The initial step—molecularization (similar to the atomiza- (P4/nmm) phase at room temperature39–41. Szafrański tion in the original Born-Haber cycle)—involves breaking and Jarek 42 reported the structures of guanidinium iodide CN H I, and the structure of tetramethylammonium iodine the CH3NH3I and PbI2 lattice structures and formation of 3 6 (CH ) NI was obtained using (CH ) NAu43 as a parent struc- CH3NH3 and PbI3 molecules 3 4 3 4 ture followed by full relaxation of their structural parameters. ∆Hmo Cubic crystal structures of CsI, CsBr, CsCl and NaCl as well CH3NH3I(s)+ PbI2(s) −−−→ CH3NH3(g)+ PbI3(g) . (2) as hexagonal PbI2 and orthorhombic PbBr2 were taken from Graystone and Wyckoff 44, Gerlach 45 . The crystal structure of The rational for using CH3NH3 and PbI3 molecules as the smallest units in the Born-Haber cycle is justified by the exis- orthorhombic PbCl2 was derived from the structure of PbBr2. × × tence of the corresponding free standing ions25,26, and will be For reciprocal space integration, 4 4 4 Monkhorst-Pack 46 × × discussed in section IV. grid was used for cubic phases, 3 3 2 were used for tetrag- × × × × The next step is the ionization of CH NH molecule onal phases, 4 4 3 for hexagonal phases and 3 6 2 for 3 3 × × orthorhombic CsPbX3 phases and 4 8 4 for ohthorhombic H PbBr and PbCl . The cutoff energy for a plane wave expan- −−−−→∆ ion,1 + 2 2 CH3NH3(g)+PbI3(g) CH3NH3 (g)+PbI3(g) , (3) sion was set at 400 eV. The lattice constant and atomic posi- tions were optimized such that residual forces acting on atoms followed by the ionization of PbI3 did not exceed 2 meV/A,˚ and the residual hydrostatic pressure was less than 50 MPa. ∆H + −−−−→ion,2 + − + – CH3NH3 (g)+ PbI3(g) CH3NH3 (g)+ PbI3 (g) . Gaseous phases, such as Cs, [CH3NH3] , [PbI3] , were (4) modelled as an individual atom/molecule surrounded by 20 A˚ It can be seen from the diagram in Fig. 1 that the formation of of vacuum. All calculations related to gaseous phases were + [CH3NH3] ion is an endothermic process, whereas the ion- performed in conjunction with optimization of internal de- ization of PbI3 is an exothermic process. The resultant ioniza- grees of freedom. Only Γ-point was used in the Brillouin tion energy is an additive of two enthalphies zone. The ionization energy of positively charged ions was calculated by subtracting the total energy of cations (e.g. Cs+, + + ∆Hion = ∆Hion,1 + ∆Hion,2 . (5) [CH3NH3] , [CN3H6] ) from the energy of neutral atoms or molecules (e.g. Cs, CH3NH3, CN3H6). Similarly, the electron + – Finally, electrically charged [CH3NH3] and [PbI3] com- affinity of negatively charged ions was modelled by adding plex ions are combined to form CH3NH3PbI3 crystalline one electron to PbCl3, PbBr3, or PbI3 molecules to form – – – structure [PbCl3] , [PbBr3] , and [PbI3] anions. The electron affinity of these ions was represented as an energy difference be- + − ∆Hlatt CH3NH3 (g)+ PbI3 (g) −−−→ CH3NH3PbI3(s) . (6) tween negatively charged complex ions and neutral species. Monopole, dipole and quadrupole corrections implemented in 47,48 The amount of energy ∆Hlatt released in this reaction is called VASP were used for eliminating leading errors and ac- the lattice energy of the hybrid organic perovskite structure. quiring accurate total energies of all charged ions. 49 This concludes the Born-Haber cycle of CH3NH3PbI3. The VESTA 3 package was used to visualize crystal structures total reaction enthalpy is compiled from enthalpies of individ- and for computing the Madelung electrostatic energy using ual steps of the cycle oxidation state as formal charges. In these calculations, the radius of ionic sphere and the reciprocal-space range were set −1 ∆Hr = ∆Hmo + ∆Hion + ∆Hlatt . (7) at 1 Aand˚ 4 A˚ , respectively. 3

+ – FIG. 1. Born-Haber cycle of hybrid halide perovskites: Methylammonium (MA) lead iodide obtained with [CH3NH3] and [PbI3] ions as elementary species.

IV. RESULTS AND DISCUSSION trend line. This suggests that formation of A+ cations and – [BX3] complex anions is a result from melting of the per- A. Lattice energies of halide perovskites ovskite structures.

Calculation of individual energies associated with various steps in the Born-Haber cycle requires subdivision of the ionic B. Stability analysis of hybrid organic halide perovksites solid in question into elementary species. In the case of alkali halides (such as NaCl, CsCl, etc.), the atomization is an appar- Now we will utilize the Born-Haber cycle in order to − ent choice. Following the same strategy, Cs+, Pb2+, and I evaluate components of the reaction enthalpy of hybrid ions can be used to calculate the lattice energy, which yields halide perovskits. The lattice energies of CH3NH3PbCl3, ≃− ∆Hlatt 29 eV (Fig. 2). CH3NH3PbBr3 and β-CH3NH3PbI3 perovskites are listed in This value agrees well with the Madelung energy of Table I. All three compounds have similar values of the lattice −27 eV obtained from the point charge model. Gopal 50 no- energies (∼10% max-min difference). However, their stabil- 7 ticed existence of a trend between the lattice energy ∆Hlatt ity characteristics are quite different. Buin et al. demon- and the melting point Tm of alkali halides with the propor- strated that under ambient conditions CH3NH3PbCl3 and − ≈ · −3 tionality factor of ∆Hlatt/Tm 7.4 10 eV/K. Assuming CH3NH3PbBr3 do not undergo a phase separation, unlike β- that the same proportionality holds for perovskite structures, CH3NH3PbI3. Both CH3NH3PbCl3 and CH3NH3PbBr3 re- − − the melting point of I II VII3 perovskites would be near main stable up to the temperature of approximately 520 K, 3900 K, which is an order of magnitude greater than the ac- above which they decompose20. Lattice energies of the − tual values of 733 888 K for group-I lead halide perovskites corresponding inorganic perovskites (CsPbI3, CsPbBr3 and 51,52 (CsPbI3, CsPbBr3, and CsPbCl3) . CsPbCl3) are very similar to their organic counterparts. In Alternatively, we can separate CsPbI3 perovskite structure fact, these inorganic perovskites are chemically stable under + – into two ions Cs and [PbI3] . The existence of the corre- the ambient environment. Remarkably, the lattice energy of β- 25,26 sponding free-standing ions was veriﬁed experimentally . CH3NH3Pbl3 and δ-CsPbCl3 structures are identical, in spite Using this approach we re-evaluated the lattice energy of of the distinct stability characteristics. Therefore, we can con- − CsPbI3 as 5.55 eV using the Born-Haber cycle similar to clude that the lattice energy cannot be used as a criteria to that shown in Fig. 1. This result translates into a substantially predict the chemical stability of compounds. lower melting point of approximately 750 K, which is remark- The analysis of various contributions to the reaction en- ably close to the experimental value of 749 K. thalpies of hybrid halide perovskites (Table I) shows that Similar calculations of the lattice energy were performed the molecularization and lattice energies largely cancel each − − − − for other inorganic I II VII3 and II IV VI3 perovksites. other. The ionization energy is the remaining contribution to Results are summarized in Table I. the reaction enthalpy in Eq. (7) that ultimately controls the The plot of the melting point vs the lattice energy of those balance of the reaction. The lower ∆Hion is, the more stable compounds is shown in Fig. 3. the compound. From this ﬁgure, we can see that the melting point of dif- Let us examine the chemical trends in ionization energy of ferent ionic structures including alkali halides follows a linear various perovskites. The total ionization energy (Eq. 5) com- 4

FIG. 2. Born-Haber cycle of inorganic halide perovskites: Caesium lead iodide obtained with Cs+, Pb2+, and I− ions as elementary species.

TABLE I. Components (eV) of the reaction enthalpies extracted from Born-Haber cycle as well as the melting temperature and stability against spontaneous decomposition for halide perovskites and other ionic structures.

Compounds ∆Hmo ∆Hlatt ∆Hion ∆Hr Tm (K) Stability 52 δ-CsPbCl3 6.23 −5.95 −0.67 −0.39 888 Y − − − 51 δ-CsPbBr3 5.80 5.72 0.46 0.39 840 Y − − − 23 δ-CsPbI3 5.72 5.55 0.34 0.17 749 Y − − − · · · CH3NH3PbCl3 6.20 6.03 0.32 0.15 Y CH3NH3PbBr3 5.80 −5.81 −0.11 −0.11 · · · Y − · · · β-CH3NH3PbI3 5.66 5.61 0.01 0.06 N − − − · · · CN3H6PbI3 5.43 5.43 0.38 0.39 Y − − − · · · (CH3)4NPbI3 5.47 4.79 1.05 0.37 Y CsCl 2.59 −6.64 0.14 −3.91 91853 Y NaCl 3.01 −8.22 1.46 −3.76 107754 Y prises of two components: the ionization energy for the cation side of the favourable range for single-junction solar cells. + + – (Cs or [CH3NH3] ) and that for the complex ion ([PbI3] , Since caesium has the lowest ionization energy in the pe- – – [PbBr3] , or [PbCl3] ). Caesium has a lower ionization en- riodic table, it is a challenging task to find molecules with ergy than CH3NH3 (Table II), which explains trends in the smaller or similar ionization energy. Among the variety of or- higher chemical stability of Cs-based perovskites as compared + + ganic cations listed in the Table II, [CN3H6] and [(CH3)4N] to their CH NH -based counterparts. + 3 3 have the ionization energies lower than that for [CH3NH3] Switching halides in the complex ions from PbI3 to PbCl3 cation making them favourable candidates for perovskites lowers their electron affinity (Table II) and, thus, leads to with improved stability. However, the size of CN3H6 and the lower total ionization energy. This explains increase of (CH3)4N molecules is significantly greater than CH3NH3, the chemical stability when changing the inorganic cage from which raises the tolerance factor above the upper formability PbI3 to PbBr3 and PbCl3. limit of 0.95 (Table III). In order to achieve a chemically stable hybrid halide per- From two structures, CN3H6PbI3 and (CH3)4NPbI3 shows ovskite structures, the necessary requirements are favourable favourable reaction enthalpies of −0.39 eV and −0.37 eV, re- geometrical factors (t-factor and octahedral factor) in conjunc- spectively(Table I). A large size of the organic molecule hin- tion with the low ionization energy (∆Hion . 0 eV). Two ders formability of CN3H6PbI3 and (CH3)4NPbI3 perovksite strategies can be used to achieve this goal: (i) find a cation structures. They both prefer hexagonal structures at ambi- with the low ionization energy or (ii) select an inorganic cage ent temperature35,36. Marco et al. 60 successfully synthesized with the low electron affinity. The second avenue is not very and characterized CN3H6PbI3 perovskite solar cells. It was promissing, since the band gap of PbBr3- and PbCl3-based hy- found that CN3H6PbI3 solar cell is also unstable under the brid perovskites (2.3 eV55 and 2.9 eV56, respectively) is out- ambient environment, which is evident from degradation of 5

go through phase separation over time. The ionization energies of onium ions in Table II corre- late with the proton affinity of the corresponding molecules62. Molecules with the low ionization energy exhibit strong proton affinity and vice versa. For instance, the proton affinity of PH3 is 785 kJ/mol, which is much lower than 901 kJ/mol for CH3NH2. It turns out that methylamine has one of the strongest proton affinity among organic compounds. There very few organic molecules (including (CH3)2NH studied here) with stronger proton affinity than CH3NH2, but none of them have a size compatible with the PbI3 cage.

V. CONCLUSIONS

The Goldschmidt’s tolerance and octahedral geometrical FIG. 3. Correlation between the lattice energy and melting temperature of ionic compounds. The line is a guide to the eye. factors do not fully capture prerequisites for formability of hybrid halide perovskites. Here we used DFT calculations in conjunction with a Born-Haber cycle to evaluate contributions TABLE II. Ionization energies (eV) of atoms and molecules. of the lattice, ionization and molecularization energies to the decomposition reaction enthalpy of hybrid halide perovskites. ∆ Ions Hion,1/2 It was previously assumed that the instability of halide per- [(CH ) N]+ 3.15 3 4 ovskite is due to a lower lattice energy of their ionic struc- Cs+ 3.85 + ture. We observe a correlation between the lattice energies [CN3H6] 3.81 + and melting temperatures, but not with reaction enthalpies [CH3NH3] 4.20 + that are ultimately linked to the chemical instability of the [HCNH2NH2] 4.30 + [NH4] 4.78 perovskites. Analysis of Born-Haber cycle components sug- Na+ 5.17 gests that the reaction enthalpy of hybrid halide perovskites is + [CH3PH3] 5.20 governed by the sum of ionization energies of a cation, e.g., + + – [CH3SH2] 5.30 [CH3NH3] , and an anion, e.g., [PbI3] . The lower total ion- + [PH4] 5.36 ization energy, the more stable is the structure, provided the + [HCNH2PH2] 8.36 geometrical conditions are fulfilled (the tolerance and octa- + [CH3] 10.0 hedral factors). This explains chemical trends in stability of – − [PbCl3] 4.52 – − hybrid and inorganic halide perovskites. For instance, the rel- [PbBr3] 4.31 – − atively high stability of CH3NH3PbCl3 is attributed to a lower [PbI3] 4.19 – ionization energy of [PbCl3] complex ion, whereas the sta- + bility of CsPbI3 is due to the lower ionization energy of Cs . The ionization energy of organic cations correlates with their the power conversion efficiency over time. Interestingly, the proton affinity. In the search for hybrid perovskite with im- rate of the efficiency decay is slower for CN3H6PbI3 as com- proved chemical stability and the band gap suitable for pho- 61 pared to CH3NH3PbI3. Szafrański found that CN3H6PbI3 tovoltaic applications, several cations were investigated. The crystals transform from orange-reddish phase to yellow phase + + promising candidates are [CN3H6] and [(CH3)4N] with the after several hours at ambient pressure and temperature. This ionization energies even lower than Cs+. The corresponding color changing demonstrates that the bandgap increases dur- CN3H6PbI3 and (CH3)4NPbI3 structures have the decomposi- ing phase transformation. From the reaction enthalpy of tion reaction enthalpy approximately 0.3 eV more favourable CN H PbI (Table I), we concludethat the drop of power con- 3 6 3 than CH3NH3PbI3. However, these ions has a prohibitively version efficiency of CN3H6PbI3 photovoltaic device is due to large size that translates into a large band gap. It is the fact the phase transformation, and the CN H PbI structure won’t 3 6 3 that CH3NH2 has the highest proton affinity among molecules of comparable size. It makes challenging to find a cation suitable for PbI3 cage as a stable activate layer for photovoltaics. TABLE III. Size of organic cations, the tolerance factor, volume of the unit cell and the band gap of selected perovskites. ACKNOWLEDGMENTS Perovskite Cation radius Tolerance Volume Bandgap (pm) factor36,57 (A˚ 3/f.u.) (eV) 36,57 Funding was provided by the Natural Sciences and En- β-CH3NH3PbI3 217 0.91 262 1.62 36,57 gineering Research Council of Canada under the Discov- CN3H6PbI3 278 1.04 326 3.38 58,59 ery Grant Program RGPIN-2015-04518. The work was per- (CH3)4NPbI3 320 1.15 361 3.30 formed using computationalresourcesof the ThunderBay Re- 6 gional Research Institute, Lakehead University, and Compute Canada(CalculQuebec).

∗ E-mail: [email protected] 34 G. Giorgi, J.-I. Fujisawa, H. Segawa, and K. Yamashita, J. Phys. 1 M. M. Lee, J. Teuscher, T. Miyasaka, T. N. Murakami, and H. J. Chem. C 119, 4694 (2015). Snaith, Science 338, 643 (2012). 35 L. Dimesso, A. Quintilla, Y.-M. Kim, U. Lemmer, and 2 N.-G. Park, J. Phys. Chem. Lett. 4, 2423 (2013). W. Jaegermann, Materials Science and Engineering: B 204, 27 3 H. S. Jung and N.-G. Park, Small 11, 10 (2015). (2016). 4 W. S. Yang, J. H. Noh, N. J. Jeon, Y. C. Kim, S. Ryu, J. Seo, and 36 G. Liu, J. Liu, Z. Sun, Z. Zhang, L. Chang, J. Wang, X. Tao, and S. I. Seok, Science 348, 1234 (2015). Q. Zhang, Inorg. Chem. 55, 8025 (2016). 5 J. A. Christians, P. A. M. Herrera, and P. V. Kamat, J. Am. Chem. 37 S. Sharma, N. Weiden, and A. Weiss, Z. Phys. Chem. 175, 63 Soc. 137, 1530 (2015). (1992). 6 S. Wozny, M. Yang, A. M. Nardes, C. C. Mercado, S. Ferrere, 38 M. Hidaka, Y. Okamoto, and Y. Zikumaru, Phys. Status. Solidi.A M. O. Reese, W. Zhou, and K. Zhu, Chem. Mater. 27, 4814 79, 263 (1983). (2015). 39 H. Ishida, H. Maeda, A. Hirano, T. Fujimoto, Y. Kubozono, 7 A. Buin, R. Comin, J. Xu, A. H. Ip, and E. H. Sargent, Chem. S. Kashino, and S. Emura, Z. Naturforsch. A. 50, 876 (1995). Mater. 27, 4405 (2015). 40 E. J. Gabe, Acta Cryst. 14, 1296 (1961). 8 J. M. Frost, K. T. Butler, F. Brivio, C. H. Hendon, M. Van Schilf- 41 E. W. Hughes and W. N. Lipscomb, J. Am. Chem. Soc. 68, 1970 gaarde, and A. Walsh, Nano Lett. 14, 2584 (2014). (1946). 9 E. Mosconi, J. M. Azpiroz, and F. De Angelis, Chem. Mater. 27, 42 M. Szafrański and M. Jarek, CrystEngComm 15, 4617 (2013). 4885 (2015). 43 P. D. Dietzel and M. Jansen, Chem. Commun. pp. 2208–2209 10 L. Zhang and P. H.-L. Sit, J. Phys. Chem. C 119, 22370 (2015). (2001). 11 J. Zhao, B. Cai, Z. Luo, Y. Dong, Y. Zhang, H. Xu, B. Hong, 44 R. W. Graystone and Wyckoff, Crystal structures, vol. 2 (Inter- Y. Yang, L. Li, W. Zhang, et al., Sci. Rep. 6, 21976 (2016). science New York, 1960). 12 E. Tenuta, C. Zheng, and O. Rubel, Sci. Rep. 6, 37654 (2016). 45 W. Gerlach, Z. Phys. A-Hadron. Nucl. 9, 184 (1922). 13 J. Burschka, N. Pellet, S.-J. Moon, R. Humphry-Baker, P. Gao, 46 H. J. Monkhorst and J. D. Pack, Phys. Rev. B 13, 5188 (1976). M. K. Nazeeruddin, and M. Grätzel, Nature 499, 316 (2013). 47 G. Makov and M. C. Payne, Phys. Rev. B 51, 4014 (1995). 14 Y. Han, S. Meyer, Y. Dkhissi, K. Weber, J. M. Pringle, U. Bach, 48 J. Neugebauer and M. Scheffler, Phys. Rev. B 46, 16067 (1992). L. Spiccia, and Y.-B. Cheng, J. Mater. Chem. A 3, 8139 (2015). 49 K. Momma and F. Izumi, J. Appl. Crystallogr. 44, 1272 (2011). 15 Y.-Y. Zhang, S. Chen, P. Xu, H. Xiang, X.-G. Gong, A. Walsh, 50 R. Gopal, J. Indian Chem. Soc. 30 (1953). and S.-H. Wei, arXiv preprint arXiv:1506.01301 (2015). 51 C. C. Stoumpos, C. D. Malliakas, J. A. Peters, Z. Liu, M. Sebas- 16 A. M. Ganose, C. N. Savory, and D. O. Scanlon, J. Phys. Chem. tian, J. Im, T. C. Chasapis, A. C. Wibowo, D. Y. Chung, A. J. Lett. 6, 4594 (2015). Freeman, et al., Cryst. Growth Des. 13, 2722 (2013). 17 F. El-Mellouhi, E. T. Bentria, S. N. Rashkeev, S. Kais, and F. H. 52 F. Fayon, Y. Vaills, P. Simon, P. Echegut, C. Bessada, J. Emery, Alharbi, Sci. Rep. 6, 30305 (2016). and J.-Y. Buzare, Ferroelectrics 185, 201 (1996). 18 V. M. Goldschmidt, Naturwissenschaften 14, 477 (1926). 53 J. Johnson, P. Agron, and M. Bredig, J. Am. Chem. Soc. 77, 2734 19 C. Li, K. C. K. Soh, and P. Wu, J. Alloy. Compd. 372, 40 (2004). (1955). 20 G. Nagabhushana, R. Shivaramaiah, and A. Navrotsky, P. Natl. 54 L. Hunter, Phys. Rev. 61, 84 (1942). Acad. Sci. USA. 113, 7717 (2016). 55 S. Ryu, J. H. Noh, N. J. Jeon, Y. C. Kim, W. S. Yang, J. Seo, and 21 L. Pauling, J. Am. Chem. Soc. 51, 1010 (1929). S. I. Seok, Energy Environ. Sci. 7, 2614 (2014). 22 M. V. Rane, A. Navrotsky, and G. A. Rossetti, J. Solid. State. 56 L. Dimesso, M. Dimamay, M. Hamburger, and W. Jaegermann, Chem. 161, 402 (2001). Chem. Mater. 26, 6762 (2014). 23 S. Sharma, N. Weiden, and A. Weiss, Z. Phys. Chem. 175, 63 57 G. Kieslich, S. Sun, and A. K. Cheetham, Chem. Sci. 5, 4712 (1992). (2014). 24 R. S. Treptow, J. Chem. Educ. 74, 919 (1997). 58 J. Palomo and P. N. Pintauro, J Membrane. Sci. 215, 103 (2003). 25 P. Cremaschi and M. Simonetta, J. Mol. Struct. 29, 39 (1975). 59 S. Garde, G. Hummer, and M. E. Paulaitis, J. Chem. Phys. 108, 26 O. E. Lanford and S. J. Kiehl, J. Am. Chem. Soc. 63, 667 (1941). 1552 (1998). 27 W. Kohn and L. J. Sham, Phys. Rev. 140, A1133 (1965). 60 N. D. Marco, H. Zhou, Q. Chen, P. Sun, Z. Liu, L. Meng, E.-P. 28 J. P. Perdew, K. Burke, and M. Ernzerhof, Phys. Rev. Lett. 77, Yao, Y. Liu, A. Schiffer, and Y. Yang, Nano Lett. 16, 1009 (2016). 3865 (1996). 61 M. Szafrański, Thermochim. Acta 307, 177 (1997). 29 G. Kresse and J. Furthmüller, Phys. Rev. B 54, 11169 (1996). 62 A. L. East, B. J. Smith, and L. Radom, J. Am. Chem. Soc. 119, 30 G. Kresse and D. Joubert, Phys. Rev. B 59, 1758 (1999). 9014 (1997). 31 P. Blöchl, Phys. Rev. B 50, 17953 (1994). 32 A. Poglitsch and D. Weber, J. Chem. Phys. 87, 6373 (1987). 33 N. Onoda-Yamamuro, T. Matsuo, and H. Suga, J. Phys. Chem. Solids. 51, 1383 (1990).