Molecular Origin of Properties of Organic ... - ACS Publications

Apr 11, 2016 - ... of Properties of Organic–Inorganic Hybrid Perovskites: The Big Picture ... We show that the electronic properties, including the ...
0 downloads 0 Views 921KB Size
Subscriber access provided by LAURENTIAN UNIV

Letter

Molecular Origin of Properties of Organic-Inorganic Hybrid Perovskites: the Big Picture from Small Clusters Hong Fang, and Puru Jena J. Phys. Chem. Lett., Just Accepted Manuscript • DOI: 10.1021/acs.jpclett.6b00435 • Publication Date (Web): 11 Apr 2016 Downloaded from http://pubs.acs.org on April 12, 2016

Just Accepted “Just Accepted” manuscripts have been peer-reviewed and accepted for publication. They are posted online prior to technical editing, formatting for publication and author proofing. The American Chemical Society provides “Just Accepted” as a free service to the research community to expedite the dissemination of scientific material as soon as possible after acceptance. “Just Accepted” manuscripts appear in full in PDF format accompanied by an HTML abstract. “Just Accepted” manuscripts have been fully peer reviewed, but should not be considered the official version of record. They are accessible to all readers and citable by the Digital Object Identifier (DOI®). “Just Accepted” is an optional service offered to authors. Therefore, the “Just Accepted” Web site may not include all articles that will be published in the journal. After a manuscript is technically edited and formatted, it will be removed from the “Just Accepted” Web site and published as an ASAP article. Note that technical editing may introduce minor changes to the manuscript text and/or graphics which could affect content, and all legal disclaimers and ethical guidelines that apply to the journal pertain. ACS cannot be held responsible for errors or consequences arising from the use of information contained in these “Just Accepted” manuscripts.

The Journal of Physical Chemistry Letters is published by the American Chemical Society. 1155 Sixteenth Street N.W., Washington, DC 20036 Published by American Chemical Society. Copyright © American Chemical Society. However, no copyright claim is made to original U.S. Government works, or works produced by employees of any Commonwealth realm Crown government in the course of their duties.

Page 1 of 25

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

The Journal of Physical Chemistry Letters

Molecular Origin of Properties of Organic-Inorganic Hybrid Perovskites: the Big Picture from Small Clusters Hong Fang, Puru Jena* Department of Physics, Virginia Commonwealth University, 701 West Grace Street, 23284, VA, United States

Abstract We show that the electronic properties, including the band gap, the gap deformation potential and the exciton binding energy, as well as the chemical stability of organic-inorganic hybrid perovskites can be traced back to their corresponding molecular motifs. This understanding allows one to quickly estimate the properties of the bulk semiconductors from their corresponding molecular building blocks. New hybrid perovskite admixtures are proposed by replacing halogens with superhalogens having compatible ionic radii. Mechanism of the boronhydride based hybrid perovskite reacting with water is investigated by using a cluster model.

TOC

1

ACS Paragon Plus Environment

The Journal of Physical Chemistry Letters

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

Endowed with high power-conversion efficiencies and low-cost manufacturing, organicinorganic hybrid perovskites with the general formula AMX3-xYx (A = organic cation; M = Pb, Sn, Ge; X,Y = I, Br, Cl) have emerged as the next-generation photovoltaic materials [1-14]. Their inorganic relatives CsSnX3 (X = I, Br, Cl) with Cs+ replacing the large organic cation have also attracted considerable interest due to their possible application as an excellent solar-cell absorber [15-18]. These materials show similar electronic and crystal structures, which usually makes them a ‘prototype’ in modeling the hybrid perovskites. Numerous studies on the bulk phases of these hybrid perovskites have been carried out to engineer their band gaps, understand the nature of their photoexcitation and improve their physical and chemical stability [1-14]. Most recently, based on the concept of Goldschmidt's tolerance factor, hundreds of possible hybrid perovskites composed by various organic cations A+, divalent metals M2+ and anionic X− moieties are proposed to exist [19]. Wide range of intriguing and diverse properties and functionalities are expected to emerge from these materials. However, given the large number of compositional possibilities, a fundamental understanding of the electronic structure and stability of hybrid perovskites is needed to facilitate the property prediction and design. One of the primary goals of cluster science has been to develop bulk properties from the molecular building blocks. Early works on semiconductors, such as Si, ZnS, GaAs, CdS, CdSe, W2O3, and PbS have shown continuous evolution of the structure and electronic properties of the bulk material from the molecular clusters [20-24]. Recognition of the electronic properties being modulated by the polymerization of the molecular motifs to different sizes and dimensions has led to the discovery of the quantum confinement effects and nonlinear optical properties [20]. Inspired by these early studies, we show in this work that the electronic, optical, and excitonic properties as well as the chemical stability of the hybrid perovskites stem from their respective 2

ACS Paragon Plus Environment

Page 2 of 25

Page 3 of 25

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

The Journal of Physical Chemistry Letters

molecular motifs. This observation allows it possible to understand and predict the properties of the bulk hybrid perovskites by using the simple molecular models. We have recently introduced the concept that the hybrid perovskites A[MX3] can be viewed as super alkali halides composed of superalkalis and superhalogens [25]. Their molecular motifs consist of the A+ cation and the [MX3]− anion bound by an ionic bond. As shown in the Supporting Table S1 [29], the organic cations A+ (as listed in Ref. 19) have lower ionization potentials (defined as the energy necessary to remove an electron from a neutral species) than lithium in the alkali group and hence can be regarded as superalkalis [26]. On the contrary, the anions [MX3]− have electron affinities (defined as the energy gained in adding an electron to a neutral species) higher than the highest value of the halogen group (chlorine) and are, therefore, named as superhalogens [27]. The nomenclature of these special ions suggests that they mimic the chemistry of the elementary halogens and alkalis, yet have distinctive properties (such as extreme size) by being atomic clusters [25].

Fig. 1 (a) Optimized molecular CH3NH3GeI3 and (b) a supercell of the experimental crystal structure of CH3NH3GeI3 from Ref. [28] using polyhedral representations of the organic cation and the [GeI3]− anion. The yellow-highlighted cation-anion couple in the supercell corresponds to the molecular motif. 3

ACS Paragon Plus Environment

The Journal of Physical Chemistry Letters

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

Table I Optimized (Calc.) average bond lengths of the isolated molecules compared to the available experimental data (Exp.) [28-29, 30-33] in the bulk crystals of the ambient phase. MA+ and FA+ denote the superalkalis [CH3NH3]+ and [HC(NH2)2]+, respectively. The bond lengths of the superalkalis in isolated molecular forms are found to be almost the same as those in the crystalline environment. The bond lengths of superhalogens coupled to different cations (given in parenthesis) show that the cations have very little influence in the geometry of the anions. For the hyperhalogens [M(BH4)3]− (M = Pb, Sn, Ge), the bond length given is that between M-B. The B-H bond length is always between 1.2 and 1.3 Å. For the hyper halogen [Ge(HCOO)3]−, the bond length is between Ge-O. Superalkali Bond (Calc.) Bond (Exp.) Superalkali Bond (Exp.) Bond (Calc.) (Å) MA+ (Å) (Å) FA+ (Å) C-N 1.49 1.49 C-N 1.31 1.31 C-H 1.09 1.09 C-H 1.09 1.09 N-H 1.04 1.04 N-H 1.01 1.02 Superhalogen Superhalogen Bond (Calc.) Bond (Exp.) Bond (Calc.) Bond (Exp.) (Å) (Å) (Å) (Å) [MX3]− [MX3]− + 2.72 (MA ) Pb-I 3.04 3.20a Sn-Br 2.89e 2.72 (Cs+) Pb-Br 2.83 2.98b Sn-Cl 2.57 2.64c Pb-Cl 2.66 2.84c Pb-(BH4) 2.77 + 2.72 (MA ) Ge-I 2.79 2.77f Sn-(BH4) 2.70(Cs+) g Ge-Br 2.55 2.53 Ge-(BH4) 2.60 Ge-Cl 2.39 2.33c Pb-(BH4) 2.77 2.99 (MA+) 3.05a Sn-I Ge-(HCOO) 1.98 2.96 (Cs+) 3.11d a Ref. [32]. b Data at 11 K from Ref. [30]. c Data of MAPbCl3, MASnCl3 and MAGeCl3 from Ref. [31] d Data of CsSnI3 from Ref. [32]. e Data of MASnBr3 from Ref. [11]. f Data of MAGeI3 from Ref. [30]. Data for CsGeI3 and FAGeI3 are 2.75 and 2.73 Å, respectively. g Data of CsGeBr3 from Ref. [33].

The optimized structures of the studied molecules are given in Figure 1 and in the Supporting Information (SI). One would imagine that the bond lengths in the isolated molecules may deviate significantly from those of the corresponding motifs in the bulk crystals, given the apparent difference between the vacuum and the crystal field. However, the calculated bond lengths of these molecules agree with the experimental data of the bulk crystals within 5% -- as shown in Table I. Such level of agreement suggests that both the superhalogens and the superalkalis (as 4

ACS Paragon Plus Environment

Page 4 of 25

Page 5 of 25

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

The Journal of Physical Chemistry Letters

identified above), to large extent, maintain their individual identities from isolated clusters to the composite ions in the bulk materials [25].

Fig. 2 (a) Onsets of the calculated UV-vis spectra for the admixture molecules CH3NH3PbI3−xBrx. The inset shows the first peaks of the spectra. x = 0-3 correspond to red, light red, cyan and blue, respectively. (b) The CCSD(T) (purple) and DFT-B3LYP (blue) calculated HOMO-LUMO gaps of the molecules compared to the experimental values (red) [6]. (c) The relative gap by shifting the x = 0 value at origin. The solid line shows the least-square fitting to the CCSD(T) relative gap by using Eq. (7). The calculated exciton binding energies of the molecules are shown as open stars.

For the hybrid perovskites, including CsSnI3−xBrx, CH3NH3PbI3-xBrx, CH3NH3SnI3−xBrx, CH3NH3PbI3−xClx, CH3NH3PbBr3−xClx and HC(NH2)2PbI3−xBrx, whose experimental band gaps are available, we calculated the UV-vis spectra of the their corresponding molecules, as shown in Figure 2a and Figure S2 in Supporting Information (SI). In each case, the onset of the spectrum is blue shifted with increasing x from 0 to 3. The same trend has also been observed in the available experimental spectra for the bulk materials [6, 11, 34-35]. The shift indicates that there is an increase of the optical gap upon gradual replacement of iodine by bromine or bromine by chlorine in both molecules and the bulk crystals. Such trend is connected to the increase of the 5

ACS Paragon Plus Environment

The Journal of Physical Chemistry Letters

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

fundamental gap -- either the HOMO-LUMO gap in the molecule or the band gap in the bulk material, as shown in Figure 2b and Figure S3 in SI. By comparing the dependence of the HOMO-LUMO gap and the band gap on x, it is interesting to see that the curves are almost parallel to each other. If we shift the lowest energy gap (when x = 0) in each case to the origin, as shown in Figure 2c and Figure S3 in SI , we can see that the curve of the HOMO-LUMO gap in the cluster is in good agreement with that of the bulk material (the largest deviation being about 0.2 eV). The CCSD(T) calculated HOMO-LUMO gaps and the TDDFT+B3LYP calculated optical gaps of the studied molecules are given in Table S2 in SI. Another exotic, yet fundamental property of these hybrid perovskites is that their band gap deformation potentials, ag =

∂Eg

∂ ln ( a )

(1)

are positive [15-16], i.e. the band gap Eg decreases with the lattice constant a, which is opposite to the behavior in most semiconductors [36]. To see if such property has a molecular origin, we calculated the change of HOMO-LUMO gaps of CsSnI3 and CsSnBr3 by changing the interatomic distance between Cs and Sn. The distance change is realized in three ways (one can refer to the schematic Figure S4 of SI for visualization): the first is to reduce the distance by evenly drawing the two atoms closer relative to the mid-point of the connecting line (Em); the second is to reduce the distance by only moving the cation along the line connecting Cs and Sn (Ec); and the third is to reduce the distance by only moving Sn along the line connecting Cs and Sn (Es). We choose CsSnI3 and CsSnBr3 here as prototype models of the hybrid perovskites. Note that changing the distance in the above fashion retains all the symmetry operations of the

6

ACS Paragon Plus Environment

Page 6 of 25

Page 7 of 25

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

The Journal of Physical Chemistry Letters

system as it would be the case when the material is under hydrostatic compression without any phase transition. The respective gap deformation potentials of the molecules are computed using Eq. (1) and are compared to those of the bulk crystals in Table II. The three calculated gaps against distance are given in Figure S4 of SI. Both Em and Es show positive band gap deformation potentials as does the band gap of the bulk crystal, while Ec shows negative ag. NBO analysis reveals that the HOMO orbital consists of Sn s and Br/I p states, and the LUMO orbital consists of Sn p and Br/I s states. The compositions of these states are found to be the same for the top of the valance and the bottom of the conduction bands, respectively, of the bulk crystals. Both Em and Es involve more overlapping atomic states of Sn and Br/I while Ec does not, which suggests that this is the reason for the positive ag of these materials. Indeed, according to the wavefunction analysis, when Sn is moved 0.08 Å towards the plane formed by the three Br (Figure S4), the total wave function becomes more blended: the contribution to the Sn-Br bond from Sn (s, p) increases by 0.08 %, while the contribution from Br (s, p) decreases the same amount. The larger ag of CsSnBr3 compared to CsSnI3 in the molecules as well as in the crystals should be due to the less diffusive electrons of Br compared to I, making the overlap between the atomic states more sensitive to the distance change. The calculated ag of Es is comparable to the gap deformation due to the change of lattice constant in the cubic hybrid perovskites [37-38]. The negative ag of Ec suggests that the bonding between Cs+ and the super halogen is similar to the ionic bonding in a CsCl molecule whose HOMO-LUMO gap will increase upon shortening the interatomic distance between Cs and Cl.

7

ACS Paragon Plus Environment

The Journal of Physical Chemistry Letters

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

Table II Calculated bandgap deformation potentials (in eV) for the Em, Ec and Es (see text for definition) of the studied molecules. These values are compared to the values of the bulk crystals obtained from the quasi-particle self-consistent GW calculations [16]. ag(Em) ag(Ec) ag(Es) ag(bulk) [15-16] CsSnI3 1.87 6.42 7.35 −0.44 10.44 8.99 CsSnBr3 2.47 −5.11

The above studies suggest that the electronic properties of the hybrid perovskites are largely decided by the bonding within the molecular motifs while the neighboring molecules should have less influence. This means that a tight-binding model [39] would be appropriate for these materials with the molecules as the basic building blocks. In such a model, the valence band of the bulk material originates from the HOMO orbital of the corresponding molecule as   () =  −  −   () ,

(2)

where EH(k) is the band energy at a wave factor k. EH is the HOMO energy of the molecule.  measures the energy change under the influence of the potential ∆V from all the other molecular motifs in the crystal,  = −  ()|∆V|  (),

(3)

where  () is the orbital of the molecular motif placed at the origin. Given that ∆V should be negative, the minus sign in Eq. (3) makes  positive.  depends on the overlap between orbitals of the central molecular motif and the neighboring ones,  = −  ()|∆V( − )|  ( − ), (4) where R is the lattice vector. As in Eq. (3), the minus sign here makes  a positive value.  () is a periodic function only related to the lattice constant and structure of the crystal. Its value is bounded by a positive upper limit FU and a negative lower limit of −FL (FL > 0). The conduction band originates from the LUMO orbital of the central molecular motif. Thus, we have an equation for the band energy EL(k), 8

ACS Paragon Plus Environment

Page 8 of 25

Page 9 of 25

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

The Journal of Physical Chemistry Letters

  () =  −  −   (),

(5)

where all the terms have the same meaning as in Eq. (2), except that they are all built from the LUMO orbital. Using Eq. (2) and (5), the band gap of the bulk crystal is written as,  = ( −  −  FL ) − ( −  +  FU ) =  HL −  −  ,

(6)

where  HL =  −  is the HOMO-LUMO gap. Both  =  −  and  =  FL +  FU should be positive, since the unoccupied orbitals are more diffusive and will be more influenced by the potential from the other molecular motifs resulting in bigger . In the crystal structures of the hybrid perovskites, if we pair the (organic) cation and the superhalogen anion to make a lattice, it will almost look like a simple cubic structure, where FL is 0 and FU is 12. By comparing the values of HOMO-LUMO gaps of the molecules and the band gaps of the crystal in Figure 2b, the FL (FU ) values suggest that the values of  ( ) and  (  ) are of the order of eV, which is normally expected from a tight-binding model [39]. By decreasing the interatomic distance, the integrals in Eq. (3) and (4) will grow bigger (larger  and  ) and therefore will further reduce the band gap by adding to the already decreased HOMO-LUMO gap  HL as shown in Table II. Referring to Eq. (6), the level of agreement between the change of the molecular gap ( HL ) with x and the crystal band gap ( ) with x, as shown in Figure 2c, suggests that the potential exerted on the central molecule from all the other molecular motifs (measured by ) and the overlap integral (  ) change little upon the change of halogens. Our calculations on the molecules do not consider the spin-orbit coupling (SOC). It has been found that, both from our calculations and others [31, 40-41], the SOC contributions to the band gap from the metal Pb, Sn and Ge are in the order of 1.0, 0.3 and 0.1 eV, respectively. For different x but the same metal, this contribution will be subtracted. The SOC contributions from the halogens will only be 9

ACS Paragon Plus Environment

The Journal of Physical Chemistry Letters

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

Page 10 of 25

partially subtracted for different x. However, these values are found to be quite small (with the largest value of about 0.1 eV for iodine) [31]. Therefore, the error introduced by the lack of SOC in these molecular calculations should be small. Another limitation of the current study is that materials with metal cation mixing cannot be studied by using the molecular model, since each motif only contains a single metal cation. Equipped with the above findings, we can now start to study other molecular admixtures and ‘extrapolate’ the properties of corresponding bulk materials from these motifs. Lead-free analogues of the hybrid perovskites are desired because of environmental considerations. Therefore, we first study the systems with lead replaced by germanium. The calculated UV-vis spectra of the molecules CH3NH3GeI3−xBrx, HC(NH2)2GeI3−xBrx and CH3NH3GeBr3−xClx show the blue shift of the onsets as shown in Figure S2 of SI, suggesting that the band gap of these crystals increases as iodine is replaced by bromine or bromine is replaced by chlorine. As already shown in the above, for the hybrid perovskites, the relative HOMO-LUMO gap of the molecular motif against x (by putting the gap for x = 0 at the origin) serves as a good estimate of the relative band gap of the bulk crystal against x. We can expand the band gap of the crystal as a second-order polynomial [36],  = ! + "# + $# %

(7)

with ! as the band gap for x = 0. Such an expansion has been successfully used to describe the experimental band gaps of the admixtures CH3NH3PbI3−xBrx [6]. Coefficients a and b in Eq. (7) can be obtained by polynomial fitting to the relative HOMO-LUMO gaps of the molecules against x, as shown in Figure 2c and Figure S3 of SI. With the fitted values of these two coefficients given in Table III, we can then use Eq. (7) to estimate the band gaps of the bulk

10

ACS Paragon Plus Environment

Page 11 of 25

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

The Journal of Physical Chemistry Letters

series CH3NH3GeI3−xBrx, HC(NH2)2GeI3−xBrx and CH3NH3GeBr3−xClx for any value of x (0 < x ≤ 3). Table III Parameters of Eq. (7) for various admixtures obtained by fitting the relative HOMO-LUMO gap of the molecules against x obtained from the CCSD(T) calculations. The x = 0 value of ! is from the available experiments. The calculated energy gaps at x = 3 using the empirical relation of Eq. (7) are compared to the experimental band gaps of the crystals. MA and FA denote for CH3NH3 and HC(NH2)2, respectively.

! (eV)

a (eV)

b (eV)

0.078 1.3a CsSnI3−xBrx b 1.57 0.220 MAPbI3−xBrx c 1.3 0.326 MASnI3−xBrx 2.15c 0.057 MASnBr3−xClx 1.57b MAPbI3−xClx 0.197 b e 2.29 /2.35 0.167 MAPbBr3−xClx 1.48f FAPbI3−xBrx 0.168 1.9g 0.300 MAGeI3−xBrx 2.63 0.095 MAGeBr3−xClx 2.2g FAGeI3−xBrx 0.291 1.57b MAPbI3−x(BH4)x −0.150 f 1.48 FAPbI3−x(BH4)x 0.129 1.3c MASnI3−x(BH4)x 0.178 g 1.9 0.040 MAGeI3−x(BH4)x g 2.2 0.442 FAGeI3−x(BH4)x 1.9g 0.125 MAGeI3−x(HCOO)x g 2.2 FAGeI3−x(HCOO)x −0.103

0.027 0.025b −0.035 0.051 0.102 0.020 0.034 −0.019 0.015 −0.029 0.248 0.171 0.089 0.150 0.022 0.197 0.295

Admixture

 (Calc.) (eV) 1.78 2.45 1.96 2.78 3.08 2.97/3.03 2.29 2.63 3.05 2.81 3.35 3.41 2.63 3.37 3.72 4.05 4.55

a

 (Exp.) (eV) 1.8a 2.29b 2.15c 3.69d 3.1e 3.1e 2.23f 3.72h -

Eb(x=0) Eb(x=3) (meV) (meV) 3.89 8.51 6.10 22.62 3.89 11.20 14.83 35.94 6.10 50.94 18.11 44.51 5.26 18.11 10.23 28.81 28.81 49.11 15.94 36.40 6.10 70.08 5.26 75.06 3.89 28.81 10.23 71.70 15.94 105.71 10.23 148.95 15.94 240.70

Experimental value from Ref. [16]. Experimental value from Ref. [6]. For MASnI3−xBrx, the fitted bowing factor b = 0.025 eV here from the molecular HOMO-LUMO gaps is equivalent in value to 0.025×9 = 0.225 eV, agreed well to the fitted bowing factor of 0.30 eV from the experimental data of band gaps [6]. c Experimental value from Ref. [31]. d Experimental value from Ref. [42]. e Experimental value from Ref. [35]. f Experimental value from Ref. [34]. g Experimental value from Ref [28]. h Derived from the PBE value and the GW value of the material at its high-temperature cubic phase (Pm−3m) [31], by assuming that the GW correction to the band gaps obtained from PBE is the same for different phases of the material [16]. b

11

ACS Paragon Plus Environment

The Journal of Physical Chemistry Letters

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

Page 12 of 25

There has been a recent report of successful synthesis of more than 30 new stable perovskite materials where halogen ions (Cl, Br, I) were replaced by [BH4]− superhalogen [43]. The effective ionic radius of the super ion [M(BH4)3]− is found to be very similar to that of [MBr3]− in the hybrid perovskites [25]. Based on these, we set out to study the properties of the admixtures of AMI3−x(BH4)x (A = CH3NH3, HC(NH2)2; M = Pb, Sn, Ge). There is a universal blue-shift of the onsets of the calculated UV-vis spectra upon the replacement of iodine by [BH4]−. What is interesting here is that the change of the HOMO-LUMO gap and the optical gap is small with partial replacement, however, the change jumps with full replacement in AM(BH4)3. These are shown in Figure 3 for CH3NH3PbI3−x(BH4)x and in Figure S2 of SI for the other admixtures. We expect the same behavior for the band gaps of the corresponding crystals according to our earlier discussions.

Fig. 3 (a) Onsets of the calculated UV-vis spectra for the admixture molecule CH3NH3PbI3−x(BH4)x. The inset shows the first peaks of the spectra. x = 0-3 correspond to red, light red, cyan and blue, respectively. (b) The CCSD(T) (purple) and DFT-B3LYP (blue) calculated HOMO-LUMO gaps of the molecules. (c) The relative gap by putting the x = 0 value at origin. The solid line shows the least square fitting to the CCSD(T) relative gap by using Eq. (7). The calculated exciton binding energies of the molecules are shown as open stars.

12

ACS Paragon Plus Environment

Page 13 of 25

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

The Journal of Physical Chemistry Letters

The reason for the small change of energy gaps on partial replacement can be found using the NBO analysis of the admixtures. In the extreme cases where we have AMI3 on one end and AM(BH4)3 on the other, the HOMO and LUMO of the former are mainly from the I(5p) and M(4p~6p) electronic states, respectively, while these orbitals of the latter are from the much lessenergetic s and p states of hydrogen and boron. This explains why the difference in the gap between these two cases is as large as shown in Figure 3 (about 2 eV). However, with only partial replacement of iodine by [BH4]−, the HOMO-LUMO gap will still be decided by the I(5p) and M(4p~6p) states of the electron due to the hybridization of these more energetic and more diffusive states with the s and p states of hydrogen and boron, respectively. The same argument goes for the crystal band structures where the valence and conduction bands originating from the I(5p) and M(4p~6p) orbitals will be much wider than the bands originating from the s and p orbitals of hydrogen and boron. This is because the electronic wave functions of the former are much more diffusive, resulting in much larger overlap-integrals in Eq. (4) and wider band according to Eq. (2). The band gap will also be much smaller in such cases according to Eq. (6). Upon partial replacement of iodine by [BH4]−, however, the top of the valence band will still be determined by the I(5p) and M(4p~6p) orbitals, since these orbitals will broaden out the top of the valence band and the bottom of the conduction band (see Figure S5 of SI). Such observation suggests that, although crystals AM(BH4)3 may show large band gaps, the admixed crystals AMI3−x(BH4)x with small x should have band gaps within the visible range based on the known small band gaps (~1.5 eV) of AMI3 [6, 34]. The fitted parameters in Eq. (7) of these new admixtures together with those of the ones studied before are summarized in Table III. It is now possible to derive the band gap of the 13

ACS Paragon Plus Environment

The Journal of Physical Chemistry Letters

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

Page 14 of 25

admixture material with any 0 < x ≤ 3 based on the band gap value for x = 0. The band gap calculated as such at x = 3 is given for each set of the admixtures. The predicted values agree well with the values from the available experiments, as shown in Table III. According to our prediction, band gaps of AM(BH4)3 should be similar to those of AMCl3. The TDDFT method together with the hybrid functional B3LYP is quite successful in the description of the optical excitation of finite systems of molecules and quantum dots [44-48]. The optical gap for each molecule is identified as the first electronic transition with non-zero oscillator strength, as shown by the first peak in Figure 2a and Figure S2 of SI . The exciton binding energy of the molecule is then the difference between the HOMO-LUMO gap and the optical gap [44]. The calculated exciton binding energies of the studied molecules are given in Table S2 of SI. The Wannier exciton binding energy is inversely proportional to the square of the relative dielectric constant which becomes smaller for larger band gaps [47]. The exciton binding energy of a single molecule is also related to its HOMO-LUMO gap [48]. By fitting to the data of a number of semiconductors (see Table S3 of SI) and the studied molecules in this work, we find a relation between the fundamental gap (Eg) and the exciton binding energy (Eb),

log10  Eb ( meV )  = 4.4538 − 4.9568 × exp  −0.1917 × Eg ( eV )  ,

(8)

as shown in Figure S4 of SI. Since the studied hybrid perovskites are all direct bandgap semiconductors, the semiconductors used here are of the same type. With such empirical relation, we can use the predicted band gaps of the hybrid perovskites from the molecular model (as given in Table III) to estimate the exciton binding energy of these materials, as shown in the last two columns of Table III. For example, the estimated exciton binding energy of the hybrid 14

ACS Paragon Plus Environment

Page 15 of 25

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

The Journal of Physical Chemistry Letters

perovskites CH3NH3PbI3 and CH3NH3PbCl3 are 6.10 and 50.94 meV, respectively, which agree well with the measured values in these materials [49-50]. According to the estimation, most studied mixed-halide perovskites are non-excitonic, while the ones with [BH4]− should have larger binding energy. We now address the stability of the new hybrid perovskites AM(BH4)3 against moisture. It is observed in the experiments [6, 50-51] that the hybrid perovskites, especially CH3NH3PbI3, are hygroscopic and readily decompose into CH3NH2, HI and PbI2 when exposed to a trace amount of water and/or intense light. In a proposed reaction of such decomposition, the water molecule serves as a catalyst [50]. It has been shown that the degradation mechanism can be successfully modeled by using merely a single molecular motif [25]. The newly invented admixtures here with [BH4]− are also expected to be hygroscopic due to ionic nature of the bond between the superalkali and superhalogen. Calculations using the molecular model reveal that the water molecule is likely to be trapped between the organic cation and the inorganic anion of [MX3]−. The binding energy between a water molecule and the molecular CH3NH3Pb(BH4)3 is 0.59 eV which is significantly smaller (12 %) than the binding energy of 0.67 eV between a water molecule and CH3NH3PbI3 [25]. It is well known that boron hydrides are easily hydrolyzed. The reaction mechanism between water and [BH4]− is as follows [52]. The first step of the reaction will be H2O + BH4− = H3B(OH)− + H2.

(9)

The reaction is based on the fact that H is negative in BH4− and, in H2O, O is strongly negative while H is positive. Therefore, H in BH4− and H in H2O will first attract each other to form H2. Meanwhile, the positive B will attract negative O to form a complex H3BOH−. There are a few 15

ACS Paragon Plus Environment

The Journal of Physical Chemistry Letters

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

Page 16 of 25

possible reaction routes afterwards between H2O and the complex H3BOH− to complete the hydrolysis process, but the first key step is to have H3BOH−. However, in the case of the hybrid perovskites, such CH3NH3[Pb(BH4)3], the reaction of the super halogen [Pb(BH4)3]− with H2O will not be the same, given the different charge states of atoms (from NBO analysis). Pb is strongly positive (+1.31e), B is strongly negative (−0.64e), and, since there are twelve H, each H will be almost neutral (−0.02e), i.e. ten times smaller than the charge on H in BH4−. Thus, the attraction between the positive H in the H2O and the H in the molecule CH3NH3[Pb(BH4)3] will be tiny and the strongly negative B in this case is unlikely to be attracted to the strongly negative O in H2O. Also, with the presence of CH3NH3+, O in H2O would rather stay with the organic cation. All these make the hydrolysis reaction in (9) difficult to happen.

16

ACS Paragon Plus Environment

Page 17 of 25

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

The Journal of Physical Chemistry Letters

Fig. 4 Calculated reaction pathway of CH3NH3Pb(BH4)3 with a water molecule (energy difference with the transition state vs. the H-O distance). The ball-and-stick plots show the reactant, transition and product states, respectively. The corresponding data points of these states are numbered accordingly. The values in these plots show the bond length in Å.

Based on the molecular model, we have been able to identify the reaction between CH3NH3Pb(BH4)3 and the water molecule as a Grotthuss mechanism [53], where one proton of NH3 in the organic cation hops onto the water molecule and one original proton of the water molecule on the other side is transferred to [Pb(BH4)3]−. Such process, including the reactant, transition and product states, is shown schematically in Figure 4. The energy barrier between the reactant and the transition state is 0.83 eV while the barrier between the product and the transition state is negligibly small, namely, 0.04 eV. Such behavior together with the relatively

17

ACS Paragon Plus Environment

The Journal of Physical Chemistry Letters

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

Page 18 of 25

small binding energy with water suggests that the studied admixtures with [BH4]− would be more resistant to moisture than the regular boronhydrides. We have shown that the properties of the organic-inorganic hybrid perovskites, including the fundamental gap, the gap deformation potential, the exciton binding energy, and the hygroscopicity, all originate from their corresponding molecular moieties. Therefore, it is possible to make quick estimations for the properties of new hybrid perovskites by using simple molecular models, which can serve as a preliminary screening process. Such process is useful given that there are hundreds of hybrid perovskites predicted to exist based on the Goldschmidt's tolerance factor [19] and the accurate bulk calculations based on DFT-HSE and relativistic effects, especially for those materials with heavy metals (such as lead), will be quite expensive. For example, one can easily extend the current study to new hybrid perovskites composed by superhalogens other than [BH4]−. [HCOO]− has one of the largest radii among the candidates of superhalogens [19]. According to the calculated tolerance factor, CH3NH3Ge(HCOO)3 and HC(NH2)2Ge(HCOO)3 can form the perovskite structures. By using the molecular model, we can readily predict the band gaps and the exciton binding energy for the admixtures MAGeI3−x(HCOO)x and FAGeI3−x(HCOO)x, as shown in Table III. A combination of the use of Goldschmidt's tolerance factor [19], the molecular-model method (as demonstrated here), and finally the accurate bulk calculation can be applied to identify colorful hybrid perovskites with improved stability and functionality.

18

ACS Paragon Plus Environment

Page 19 of 25

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

The Journal of Physical Chemistry Letters

METHODS Our calculations on clusters are carried out using the GAUSSIAN03 package [54]. The hybrid density functional theory (DFT) method with Becke three parameter Lee-Yang-Parr (B3LYP) [55-56] prescription for the exchange-correlation energy is used. The fundamental gaps of the molecules are calculated using the coupled cluster method CCSD(T). The basis set used for halogen ions is 6-31+G*. For Ge, Sn and Pb, aug-cc-pVDZ basis sets with the effective core potential are used. For each admixture of AMX3−xYx (x = 0-3), the corresponding molecular motif is optimized to its ground state without any imaginary frequency. The ultra-violet-visible (UV-vis) spectrum of the molecule is then computed from 100 singlet states by using the timedependent DFT (TD-DFT) [57]. Natural bond orbital (NBO) analysis is used to obtain the atomic charge state and the bond composition of the molecule. The binding energy between two molecules is calculated using the counterpoise correction to avoid the basis set superposition error. DFT calculations for the crystalline phase are performed with the projector augmented wave (PAW) [58] and the Perdew-Burke-Ernzerhof generalized gradient approximation (PBEGGA) [59] using the VASP package [60-61]. The electronic structure calculations are performed with and without the spin-orbit coupling (SOC). The plane-wave cutoff energy is set to 600 eV. Monkhorst-Pack 6 × 6× 6 k-meshes [62] are adopted.

19

ACS Paragon Plus Environment

The Journal of Physical Chemistry Letters

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

Page 20 of 25

ASSOCIATED CONTENT Supporting Information Calculated ionization potential and electron affinity for organic cations and inorganic anions in the hybrid perovskites, calculated UV-vis spectra for the studied molecular admixtures, calculated HOMO-LUMO gaps and the exciton binding energyies for the studied molecules, data of calculated HOMO-LUMO gaps and optical gaps, schematic plots showing the three ways to change the atomic distance, calculated gap deformation potentials, DFT-PBE calculated electronic density of states (DoS) for the crystalline admixtures MASnI2(BH4) and MAGeI2(BH4), experimental data of band gaps and exciton binding energy of direct bandgap semiconductors, and collected data of direct bandgap semiconductors and the studied molecules.

AUTHOR INFORMATION Corresponding Author *E-mail: [email protected] Notes The authors declare no competing financial interest.

ACKNOWLEDGEMENTS This work is supported in part by the U.S. Department of Energy, Office of Basic Energy Sciences, Division of Materials Sciences and Engineering under Award # DE-FG02-96ER45579. Resources of the National Energy Research Scientific Computing Center supported by the Office of Science of the U.S. Department of Energy under Contract No. DE-AC02-05CH11231 is acknowledged. We also thank the National Institute for Computational Sciences for letting us use their computational resources (Darter).

REFERENCES (1) Mitzi, D. B.; Field, C. A.; Harrison, W. T. A.; Guloy, A. M. Conducting Tin Halides with a Layered Organic-Based Perovskite Structure. Nature 1994, 369, 467-469. (2) Lee, M. M.; Teuscher, J.; Miyasaka, T.; Murakami, T. N.; Snaith, H. J. Efficient Hybrid Solar Cells Based on Meso-Superstructured Organometal Halide Perovskites. Science 2012, 338, 643647. (3) Etgar, L.; Gao, P.; Xue, Z.; Peng, Q.; Chandiran, A. K.; Liu, B.; Nazeeruddin, M. K.; Graetzel, M. Mesoscopic CH3NH3PbI3/TiO2 Heterojunction Solar Cells. J. Am. Chem. Soc. 2012, 134, 17396-17399. (4) Hodes, G. Perovskite-Based Solar Cells. Science 2013, 342, 317-318. 20

ACS Paragon Plus Environment

Page 21 of 25

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

The Journal of Physical Chemistry Letters

(5) Loi, M. A.; Hummelen, J. C. Hybrid Solar Cells: Perovskites under the Sun. Nature Materials 2013, 12, 1087-1089. (6) Noh, J. H.; Im, S. H.; Heo, J. H.; Mandal, T. N.; Seok, S. I. Chemical Management for Colourful, Efficient, and Stable Inorganic-Organic Hybrid Nanostructured Solar Cells. Nano Letters 2013, 13, 1764-1769. (7) Snaith, J. H. Perovskites: the Emergence of a New Era for Low-Cost, High-Efficiency Solar Cells. J. Phys. Chem. Lett. 2013, 4, 3623-3630. (8) Park, N. G. Organometal Perovskite Light Absorbers Toward a 20% Efficiency Low-Cost Solid-State Mesoscopic Solar Cell. J. Phys. Chem. Lett. 2013, 4, 2423-2429. (9) Bisquert, J. The Swift Surge of Perovskite Photovoltaics. J. Phys. Chem. Lett. 2013, 4, 25972598. (10) Lotsch, B. V. New Light on an Old Story: Perovskites Go Solar. Angew. Chem. Int. Ed. 2014, 53, 635-637. (11) Hao, F.; Stoumpos, C. C.; Cao, D. H.; Chang, R. P. H.; Kanatzidis, M. G. Lead-Free SolidState Organic-Inorganic Halide Perovskite Solar Cells. Nature Photonics 2014, 8, 489-494. (12) Tan, Z.; Moghaddam, R. S.; Lai, M. L.; Docampo, P.; Higler, R.; Deschler, F.; Price, M.; Sadhanala, A.; Pazos, L. M.; Credgington, D.; et al. Bright Light-Emitting Diodes Based on Organometal Halide Perovskite. Nature Nanotechnology 2014, 9, 687-692. (13) Xing, G.; Mathews, N.; Lim, S. S.; Yantara, N.; Liu, X.; Sabba, D.; Graetzel, M.; Mhaisalkar, S.; Sum, T. C. Low-Temperature Solution-Processed Wavelength-Tunable Perovskites for Lasing. Nature Materials 2014, 13, 476-480. (14) Deschler, F.; Price, M.; Pathak, S.; Klintberg, L. E.; Jarausch, D.; Higler, R.; Huettner, S.; Leijtens, T.; Stranks, S. D.; Snaith, H. J.; et al. High Photoluminescence Efficiency and Optically Pumped Lasing in Solution-Processed Mixed Halide Perovskite Semiconductors. J. Phys. Chem. Lett. 2014, 5, 1421-1426. (15) Borriello, I.; Cantele, G.; Ninno, D. Ab Initio Investigation of Hybrid Organic-Inorganic Perovskites Based on Tin Halides. Phys. Rev. B 2008, 77, 235214. (16) Huang, L.; Lambrecht, W. R. L. Electronic Band Structure, Phonons, and Exciton Binding Energies of Halide Perovskites CsSnCl3, CsSnBr3, and CsSnI3. Phys. Rev. B 2013, 88, 165203. (17) Chen, Z.; Wang, J. J.; Ren, Y.; Yu, C.; Shum, K. Schottky Solar Cells Based on CsSnI3 Thin-Films. Appl. Phys. Lett. 2012, 101, 093901. 21

ACS Paragon Plus Environment

The Journal of Physical Chemistry Letters

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

Page 22 of 25

(18)Yu, C.; Chen, Z.; Wang, J. J.; Pfenninger, W.; Vockic, N.; Kenney, J. K.; Shum, K. Temperature Dependence of the Band Gap of Perovskite Semiconductor Compound CsSnI3. J. Appl. Phys. 2011, 110, 063526. (19) Kieslich, G.; Sun, S.; Cheetham, A. K. Solid-State Principles Applied to Organic–Inorganic Perovskites: New Tricks for an Old Dog. Chem. Sci. 2014, 5, 4712-4715; An Extended Tolerance Factor Approach for Organic–Inorganic Perovskites. Chem. Sci. 2015, 6, 3430-3433. (20) Steigerwald, M. L.; Brus, L. E. Semiconductor Crystallites: a Class of Large Molecules. Acc. Chem. Res. 1990, 23, 183-188. (21) Soloviev, V. N.; Eichhoefer, A.; Fenske, D.; Banin, U. Molecular Limit of a Bulk Semiconductor:  Size Dependence of the “Band Gap” in CdSe Cluster Molecules. J. Am. chem. Soc. 2000, 122, 2673-2674. (22) Kasuya, A.; Sivamoohan, R.; Barnakov, Y. A.; Dmitruk, I. M.; Nirasawa, T.; Romanyuk, V. R.; Kumar, V.; Mamykin, S. V.; Tohji, K.; Jeyadevan, B.; et al. Ultra-Stable Nanoparticles of CdSe Revealed from Mass Spectrometry. Nature Materials 2004, 3, 99-102. (23) Zhai, H. J.; Huang, X.; Cui, L. F.; Li, X.; Li, J.; Wang, L. S. Electronic and Structural Evolution and Chemical Bonding in Ditungsten Oxide Clusters: W2On−and W2On (n = 1−6). J. Phys. Chem. A 2005, 109, 6019-6030. (24) Chandler, R. R.; Bigham, S. R.; Coffer, J. L. Spectroscopic Analysis of Semiconductor Colloids: an Experiment in Materials Science for the Advanced Inorganic or Physical Chemistry Laboratory. J. Chem. Educ. 1993, 70, A7-A10. (25) Fang, H.; Jena, P. J. Super-Ion Inspired Colorful Hybrid Perovskite Solar Cells. J. Mater. Chem. A 2016, 4, 4621-4621. (26) Gutsev, G. L.; Boldyrev, A. I. DVM Xα Calculations on the Electronic Structure of “Superalkali” Cations. Chem. Phys. Lett. 1982, 92, 262-266. (27) Gutsev, G. L.; Boldyrev, A. I. DVM Xα Calculations on the Ionization Potentials of MX−(k+1) Complex Anions and the Electron-Affinities of MKk+1 Superhalogens. Chem. Phys. Lett. 1981, 56, 277-283. (28) Stoumpos, C. C.; Frazer, L.; Clark, D. J.; Kim, Y. S.; Rhim, S. H.; Freeman, A. J.; Ketterson, J. B.; Jang, J. I.; Kanatzidis, M. G.; Hybrid Germanium Iodide Perovskite Semiconductors: Active Lone Pairs, Structural Distortions, Direct and Indirect Energy Gaps, and Strong Nonlinear Optical Properties. J. Am. Chem. Soc. 2015, 137, 6804-6819. 22

ACS Paragon Plus Environment

Page 23 of 25

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

The Journal of Physical Chemistry Letters

(29) Stoumpos, C. C.; Malliakas, C. D.; Kanatzidis, M. G. Semiconducting Tin and Lead Iodide Perovskites with Organic Cations: Phase Transitions, High Mobilities, and Near-Infrared Photoluminescent Properties. Inorg. Chem. 2013, 52, 9019-9038. (30) Swainson, I. P.; Hammond, R. P.; Soulliere, C.; Knop, O.; Massa, W.; Phase Transitions in the Perovskite Methylammonium Lead Bromide, CH3ND3PbBr3. J. Sol. Stat. Chem. 2003, 176, 97-104. (31) Katan, C.; Pedesseau, L.; Kepenekian, M.; Rolland, A.; Even, J. Interplay of Spin-Orbit Coupling and Lattice Distortion in Metal Substituted 3D Tri-Chloride Hybrid Perovskites. J. Mater. Chem. A 2015, 3, 9232-9240. (32) Yamada, K.; Mikawa, K.; Okuda, T.; Knight, K. S. Static and Dynamic Structures of CD3ND3GeCl3 Studied by TOF High Resolution Neutron Powder Diffraction and Solid State NMR. J. Chem. Soc., Dalton Trans. 2002, 2112-2118. (33) Thiele, G.; Rotter, H. W.; Schmidt, K. D. Kristallstruckturen und Phasentransformationen von Caesiumtrihalogenogermanaten (II) CsGeX3 (X = Cl, Br, I). Z. Anorg. Allg. Chem. 1987, 545, 148. (34) Eperon, G. E.; Stranks, S. D.; Menelaou, C.; Johnston, M. B.; Herz, L. M.; Snaith, H. J. Formamidinium Lead Trihalide: a Broadly Tunable Perovskite for Efficient Planar Heterojunction Solar Cells. Energy Environ, Sci. 2014, 7, 982-988. (35) Kitazawa, N.; Watanabe, Y.; Nakamura, Y. Optical Properties of CH3NH3PbX3 (X = Halogen) and Their Mixed-Halide Crystals. J. Mater. Sci. 2002, 37, 3585-3587. (36) Piprek, J. Semiconductor Optoelectronic Devices Introduction to Physics and Simulation, Academic Press, 2002. (37) Filip, M. R.; Eperon, G. E.; Snaith, H. J.; Giustino, F. Steric Engineering of Metal-Halide Perovskites with Tunable Optical Band Gaps. Nat. Commun. 2014, 5, 5757. (38) Kim, J.; Lee, S-C.; Lee, S-H.; Hong, K-H. Importance of Orbital Interactions in Determining Electric Band Structures of Organo-Lead Iodide. J. Phys. Chem. C. 2015, 119, 4627-4634. (39) Omar, M. A. Elementary Solid State Physics: Principles and Applications, Addison-Wesley, 1994.

23

ACS Paragon Plus Environment

The Journal of Physical Chemistry Letters

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

Page 24 of 25

(40) Even, J.; Pedesseau, L.; Jancu, J.; Katan, C. Importance of Spin-Orbit Coupling in Hybrid Organic/Inorganic Perovskites for Photovoltaic Aapplications. J. Phys. Chem. Lett. 2013, 4, 2999-3005. (41) Even, J.; Pedesseau, L.; Dupertuis, M. A.; Jancu,J. M.; Katan, C. Electronic Model for SelfAssembled Hybrid Organic/Perovskite Semiconductors: Reverse Band Edge Electronic States Ordering and Spin-Orbit Coupling. Phys. Rev. B 2012, 86, 205301. (42) Chiarella, F.; Zappettini, A.; Licci, F. Combined Experimental and Theoretical Investigation of Optical, Structural, and Electronic Properties of CH3NH3SnX3 Thin Films (X = Cl, Br). Phys. Rev. B 2008, 77, 045129. (43) Schouwink, P.; Ley, M. B.; Tissot, A.; Hagemann, H.; Jensen, T. R.; Smrčok, L.; Černý, R. Structure and Properties of Complex Hydride Perovskite Materials. Nature Communication 2014, 5, 5706. (44) Garoufalis, C. S.; Zdetsis, A. D. High Level Ab Initio Calculations of the Optical Gap of Small Silicon Quantum Dots. Phys. Rev. Lett. 2001, 87, 276402. (45) Williamson, A. J.; Grossman, J. C.; Hood, R. Q.; Puzder A.; Galli, G. Quantum Monte Carlo Calculations of Nanostructure Optical Gaps: Application to Silicon Quantum Dots. Phys. Rev. Lett. 2002, 89, 196803. (46) Ullrich, C. A.; Yang, Z. Excitons in Time-Dependent Density-Functional Theory. Top. Curr. Chem. Springer International Publishing Switzerland, 2014, DOI: 10.1007/128_2014_610. (47) Haug, H.; Koch, S. W. Quantum Theory of the Optical and Electronic Properties of Semiconductors, World Scientific, 2009. (48) Knupfer, M. Exciton Binding Energies in Organic Semiconductors. Appl. Phys. A 2003, 77, 623-626. (49) D'Innocenzo, V.; Grancini, G.; Alcocer, M. J. P.; Kandada, A. R. S.; Stranks, S. D.; Lee, M. M.; Lanzani, G.; Snaith, H. J.; Petrozza, A. Excitons Versus Free Charges in Organo-Lead TriHalide Perovskites, Nature Commun. 2014, 5, 3586. (50) Frost, J. M.; Butler, K. T.; Brivio, F.; Hendon, C. H.; van Schilfgaarde, M.; Walsh, A. Atomistic Origins of High-Performance in Hybrid Halide Perovskite Solar Cells. Nano lett. 2014, 14, 2584-2590.

24

ACS Paragon Plus Environment

Page 25 of 25

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

The Journal of Physical Chemistry Letters

(51) Deretzis, I.; Alberti, A.; Pellegrino, G.; Smecca, E.; Giannazzo, F.; Sakai, N.; Miyasaka, T.; La Magna, A. Atomistic Origins of CH3NH3PbI3 Degradation to PbI2 in Vacuum. Appl. Phys. Lett. 2015, 106, 131904. (52) Zhang, Y.; Li, M.; Yang, B. Reaction Mechanism of Hydrogen Generation from Hydrolysis of BH4−: Quantum Chemistry Computing. Acta Chim. Sinica 2013, 71, 769-776. (53) Agmon, N. The Grotthuss Mechanism. Chem. Phys. Lett. 1995, 244, 456-462. (54) Frisch, M. J.; Trucks, G. W.; Schlegel, H. B.; Scuseria, G. E.; Robb, M. A.; Cheeseman, J. R.; Jr. Montgomery, J. A.; Treven, T.; Kudin, K. N.; Burant, J. C.; et al. Gaussian 03, revision B.03; Gaussian, Inc.: Wallingford, CT, 2003. (55) Becke, A. Density-Functional Thermochemistry. III. The Role of Exact Exchange. J. Chem. Phys. 1993, 98, 5648-5652. (56) Lee, C.; Yang, W.; Parr, R. G. Development of the Colle-Salvetti Correlation-Energy Formula into a Functional of the Electron Density. Phys. Rev. B 1988, 37, 785. (57) Runge, E.; Gross, E. K. U. Density-Functional Theory for Time-Dependent Systems. Phys. Rev. Lett. 1984, 52, 997. (58) Blöchl, P. E. Projector Augmented-Wave Method. Phys. Rev. B 1994, 50, 17953. (59) Perdew, J. P.; Burke, K.; Ernzerhof, M. Generalized Gradient Approximation Mad Simple. Phys. Rev. Lett. 1996, 77, 3865. (60) Kresse, G.; Furthmueller, J. Efficiency of Ab-Initio Total Energy Calculations for Metals and Semiconductors Using a Plane-Wave Basis Set. J. Comput. Mater.Sci. 1996, 6, 15-50. (61) Kresse, G.; Furthmueller, J. Efficient Iterative Schemes for Ab Initio Total-Energy Calculations Using a Plane-Wave Basis Set. Phys. Rev. B 1996, 54, 11169. (62) Monkhorst, H. J.; Pack, J. D. Special Points for Brillouin-Zone Integrations. Phys. Rev. B 1976, 13, 5188.

25

ACS Paragon Plus Environment