Low Density of Conduction and Valence Band ... - ACS Publications

Jan 4, 2017 - School of Chemistry, The University of Melbourne, Parkville, Victoria 3010, Australia. ‡. School of Materials Science and Engineering,...
0 downloads 0 Views 690KB Size
Subscriber access provided by Fudan University

Article

The Low Density of Conduction and Valence Band States Contribute to the High Open Voltage in Perovskite Solar Cells Yecheng Zhou, and Guankui Long J. Phys. Chem. C, Just Accepted Manuscript • DOI: 10.1021/acs.jpcc.6b10914 • Publication Date (Web): 04 Jan 2017 Downloaded from http://pubs.acs.org on January 4, 2017

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 C 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 27

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

The Low Density of Conduction and Valence Band States Contribute to the High Open Voltage in Perovskite Solar Cells Yecheng Zhou† and Guankui Long∗,‡ †School of Chemistry, The University of Melbourne, Parkville, VIC 3010, Australia ‡School of Materials Science and Engineering, Nankai University, Tianjin 300071, China E-mail: [email protected]

1

ACS Paragon Plus Environment

The Journal of Physical Chemistry

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

Abstract Hybrid perovskites are widely used for high performance solar cells. Large diffusion lengths and long charge carrier lifetimes are considered as two main factors for their high performance. Here, we argue that not only large diffusion lengths and long carrier lifetimes, but also the low densities of the conduction and valence band states (Nc , Nv ) contribute to highperformance perovskite solar cells. We estimated Nc and Nv of silicon, CdTe, and typical perovskites with two different methods. It was found that the Nc and Nv of perovskites and CdTe are much lower than that of silicon. Using numerical models, we found that the solar cell of a material with same characteristics but lower Nc and Nv can realize a higher open-circuit voltage (Voc ) and higher power conversion efficiency (PCE). We put forward and proved that the low Nc and Nv in hybrid perovskite is one of the factors for its high performance. This provides a new guideline for finding and developing new photovoltaic candidate materials.

2

ACS Paragon Plus Environment

Page 2 of 27

Page 3 of 27

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

Introduction Open-circuit voltage (Voc ) of a solar cell is always lower than the absorber band gap. The difference between them is called open voltage deficit. A large Voc deficit means large potential energy lost and leads to a lower performance. The effort to reduce the open voltage deficit never stop. Most of this kind of work is focusing on fabrication technics and methods, 1–3 few of them focus on the influence from the intrinsic properties of perovskites. Here, we discuss the open voltage deficit from the aspect of material itself and propose that we can reduce the open voltage deficit by using materials with low DOS and small effective mass. Perovskite solar cells have realized Power Conversion Efficiencies (PCEs) up to 22 %. 4–9 Excellent electron and hole diffusion ability 10,11 and a very wide absorption wavelength range with high IPCE 12–14 are reckoned as two main factors for perovskite solar cell’s high performance. Both theoreticians 15–19 and experimentalists 10,11,20,21 are trying to find more fundamental reasons. They attribute the high performance of perovskite solar cell to its proper band gap and large static dielectric constant. 16,19,22,23 Others believe that polarization of perovskite benefits to its high PCEs through grain boundary 24–27 or screening effect. 23 If we examine it carefully, it will be found that the reason of excellent electronic properties found by both experimental and theoretical works only benefits to the high output current. Whereas the maximum Jsc (21 mA/cm2 ) of perovskite solar cells 8 is much lower than the Jsc (30 mA/cm2 ) of CdTe solar cells, 28 which has similar band gaps. Therefore, there must be other reasons. The other most possible reason is the high Voc , which hasn’t been fully discussed for perovskite solar cells. Before our investigation, we collected experimental band gaps and Voc related data of different type solar cells. For amorphous silicon, its band gap ranges from 1.55 to 2.10 eV, 29 whereas its highest Voc is only about 0.98V. 30 Its voltage deficit is 570 mV. For multi-crystalline silicon solar cells, the highest Voc is 660 mV, which is 460 mV lower than its band gap at 1.12 eV. 31 The maximum Voc of CdTe solar cells is 1.10 V, 32–34 which is 400 mV lower than its band gap of 1.5 eV at 300K. 35 However, the high Voc in CdTe 3

ACS Paragon Plus Environment

The Journal of Physical Chemistry

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

solar cell is obtained by scarifying its PCE. For CdTe solar cells with Voc larger than 1.0 V, the maximum PCE is 13.6%. 34 For the best performance CdTe solar cell the voltage deficit is 630 meV. 31 For a perovskite solar cell, its Voc reaches 1.14 V corresponding to the band gap of 1.57 V. 8 The voltage deficit is 430 mV. It is obvious that the PCE is higher when its voltage deficit is lower. Therefore, we are going to discuss some factors related to Voc and check whether there is a rule governing solar cell’s Voc and performance. Open-circuit voltage deficits were discussed in organic solar cells, 36–39 quantum dot solar cells, 2 a-Si:H/c-Si heterojunction solar cells. 40,41 In quantum dot solar cells, it was found that sub-bandgap states are the most likely origin of the high Voc deficit. 2 For bulk heterojunction solar cells, theoretical works showed that the Voc depends on (1) the donor-acceptor energy gap; (2) charge-carrier recombination rates; (3) illumination intensity; (4) the contact work functions (if not in the pinning regime); and (5) the amount of energetic disorder. 40,41 In 2013, Edri et al. showed that perovskite solar cells can realize high Voc up to 1.3 V. While, this high Voc is based on its large band gap of 2.3 eV. Its open-circuit voltage deficit is as large as 1.0 V. Shao’s experiment showed that reducing the energy disorder in the fullerene electron transport layer through a simple solvent annealing process increases the Voc from 1.04 V to 1.13 V. 42 It was also found that the tuning of the work function can improve Voc in the devices. 43 Chen et al. demonstrated that the limiting of Voc of CH3 NH3 PbBr3 perovskite solar cells is due to the interface loss induced by the charge extraction layer rather than by bulk dominated recombination losses. 44 In theory, Yang et al. illustrated that the reduced bimolecular recombination is the origin of high Voc . 45 With numerical simulations, we also showed the inert interface defects, long charge carrier lifetime and strong incident light intensity increase Voc . 46 Kirchartz et al. showed that photon recycling enhances the Voc in the range of 10-100 mV. 47 The upper limit of Voc was estimated to be 1.24-1.32 V based on different recombination mechanisms. 48 For organic solar cells, Garcia-Belmonte et al. found that small tails of DOS distribution will lead to high Voc . 36 Sulas et al. also showed that the increased density of charge-transfer

4

ACS Paragon Plus Environment

Page 4 of 27

Page 5 of 27

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

states causes Voc loss. 37 The same conclusion was observed by Collins et al.. 39 So the question is whether the DOSs in perovskite contribute to its solar cell high Voc ? Especially, there are some works reported that perovskites exhibit low DOSs at the edges of the valence and conduction band. 49,50 In this paper, we are going to show that besides the defects, charge carrier lifetime and incident light intensity, there are other parameters that affects Voc . They are Nc and Nv . With two different methods based on ab initio calculations, we calculated the effective conduction band and valence band DOSs (Nc and Nv ) of silicon, CdTe and typical hybrid perovskites, such as CH3 NH3 PbI3 51,52 and HC(NH2 )2 PbI3 , 53–57 which have been widely used in solar cells. It was found that Nc and Nv of hybrid perovskites and CdTe are several times lower than that of silicon. The low Nc and Nv are in good agreement with the low DOS observed by Endres et al. 49,50 By implementing these Nc and Nv into the numerical model we developed for perovskite solar cell, 46 we found that the low Nc and Nv in perovskite materials could contribute 100 mV higher Voc compared to the ‘perovskite solar cell’ with silicon’s Nc and Nv .

Methods DFT calculation details In this paper, we calculated Nc and Nv of typical perovskites, silicon and CdTe in order to find out DOS effect on their power conversion performance. Perovskite structures came from Ref. 58. All of the calculations are performed by VASP. 59,60 Silicon band structure based on GGA-PS functional, 61 GW method 62,63 and HSE functional 64,65 are carried out and compared in order to find the optimum method. It is found that the GW method and HSE functional do give correct band gaps of 1.10 eV and 1.14 eV comparing to 1.12 eV of experiments. Whereas, their effective masses are very similar to the masses from GGA-PS functional. Effective mass calculations in MAPbI3 have been carried out with HSE functional 5

ACS Paragon Plus Environment

The Journal of Physical Chemistry

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 6 of 27

and GW methods. 18,66 The harmonic mean effective electron and hole masses of HSE+SOC calculations are 0.18 m0 and 0.22 m0 respectively. 66 Their average effective masses are 0.19 m0 and 0.25 m0 for the method of SOC+GW. 18 Both of these effective masses close to 0.1920.199 m0 23 and 0.21 m0 (Shown in Table S6) in our GGA-PS calculations. The effective mass error between our GGA-PS method and these two highly accurate methods is smaller than the error between these two accurate methods themselves. This confirms that GW method and HSE functional do give correct band gaps while their effective masses are the same with GGA-PS functional. For the sake of computational time limitation and little accuracy sacrifice, we use GGA-PS functional in this work. The effective masses are the quadratic coefficients of the fitting quadratic function with five points near the high-symmetry k-point. Fittings are performed by R code. 67 Detailed processes and methods are described in SI.

Nc and Nv calculation from DOS(E) When the charge concentration is low, the interaction between charges is weak and can be neglected. We assume there is no interaction between charge carriers, the charge carrier density follows the Maxwell-Boltzmann distribution. It is expressed by: f (E) = f0 exp(−

E−Ef ), kT

where, E is the energy level of charges, Ef is the Fermi energy level, and f0 is the charge density at Fermi level. Thus, the electron density in the conduction band is given by

n= =

Z



c ZE∞

DOS(E) × f (E)dE = DOS(E − Ec )e

E−Ec −kT

Z



DOS(E) × e

E−Ef n −kT

dE/Vcell

Ec

dE/Vcell × e

Ec −Ef n −kT

= Nc × e

Ec −Ef n −kT

(1)

Ec

where, Ec and Ef n are the lowest energy of the conduction band and quasi Fermi level, Vcell is the volume of the unit cell, T is the temperature and k is the Boltzmann constant, DOS(E)

6

ACS Paragon Plus Environment

Page 7 of 27

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

is shown in Figure 1. Nc is calculated by:

Nc =

Z

∞ Ec

DOS(E − Ec ) E−E c e −kT dE Vcell

(2)

The integration of DOS(E) is very sensitive to calculation methods, in this work we used tetrahedron method with Blochl corrections (ISMEAR=-5) as recommended by VASP for accurate electronic DOS calculation.

Nc and Nv calculation from effective mass Theoretically, the density-of-states function is given by 1 3 DOS(E) = 2 2π

Z

dSk , |∇k E|

(3)

where, 2 is for the spin, dSk = 4πk 2 for the integral of the energy surface. For free electrons, its energy is calculated by E =

¯ 2 k2 h . 2m0

Thus ∇k E =

¯2k h . m0

Substituting them into Equation 3,

we obtain

DOS(E) =

4π  1  4πk 2 (2m0 )3/2 E 1/2 . = 2 3 4π 3 h h ¯ k/m0

(4)

Equation 4 is based on free electrons with parabolic energy band. This equation can also be used for semiconductors with parabolic energy bands just by replacing the free electron mass with the effective mass of electron or hole. 68 Choosing the conduction band edge as the reference level, Equation 4 can be expressed by

DOS(E − Ec ) =

4π  (2m∗ )3/2 (E − Ec )1/2 . h3

7

ACS Paragon Plus Environment

(5)

The Journal of Physical Chemistry

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 8 of 27

Finally, the electron density can be calculated by 68,69 ∞ 1 f (E)DOS(E − Ec )dE = DOS(E − Ec )dE (Ef −E)/kB T ] Ec [1 + e Ec Z ∞ (E − Ec )1/2 dE 4π  ∗ 3/2 = 3 (2m kB T ) (Ef −E)/kB T ] h Ec [1 + e

n=

Z



Z

(6)

= Nc F1/2 (η), where T m∗c 3/2 ) 300m0 R∞ is the effective density of the conduction band states; F1/2 (η) = √2π Ec Nc = 2(2πm∗ kB T /h2 )3/2 = 2.5409 × 1019 (

(7) (E−Ec )1/2 dE [1+e

(Ef −E)/kB T

]

is

the Fermi integral.

For the non-spherical constant-energy surface, the effective mass is replaced by the reduced effective mass for DOS, which is given by

m∗ =

p 3 g 2 ml mt mt

(8)

where g is the number of equivalent band minima, ml and mt are the longitudinal and transverse effective masses respectively. Equation 7 is calculated for the DOS calculation of materials with parabolic bands. The Nc and Nv are dominated by the bottom of the conduction band or the top of the valence band, the bottom and the top can be fitted by parabolic curves without bringing about significant errors. For hybrid perovskites, the conduction band bottom is strongly affected by MA+ orientations. 23 Under certain cases, 70 hybrid perovskites can show two valleys. We think the Nc and Nv with and without disordered MA+ should be similar. On the other hand, in this paper, all of the conduction bands show one bottom. All of the fittings are reasonable without significant variation, as shown in the SI. For a material has several conduction bands, the total Nc is the summation of these conduction band Nc if their energy levels are the same. If their energy levels have Ei −E0 P small difference, the total DOS is calculated by N = ni=0 Ni × e −kT , where Ni and Ei are 8

ACS Paragon Plus Environment

Page 9 of 27

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

the effective DOS and the energy level of the i th band.

Results Known factors influence on Voc Open circuit voltages are always lower than their band gaps, which is due to charge carrier densities are lower than the effective DOSs of conduction band and valence band. In silicon solar cells, it was already found that Voc is related to the charge density. 71–73 We pointed out that the output voltage of a solar cell is determined by its charge densities at two electrodes and Nc and Nv : 46,48 p|x=d n|x=0 ) − (Ev − kT ln( )) Nc Nv n|x=0 p|x=d = Ebgap + kT ln( ) + kT ln( ) Nc Nv

eV = Ef n − Ef p = Ec + kT ln(

(9)

= Ebgap + kT ln(n|x=0 p|x=d ) − kT ln(Nc Nv ) According to Equation (9), there are three approaches to reduce the voltage deficit. The first is operating solar cells at low temperature. If the working temperature close to zero Kelvin, then kT ≈ 0, which leads to V ≈ Ebgap . But keeping solar cells work under low temperature may cost more. This conclusion agrees with experiment and theory done by Sachenko et al.. 74 The second is increasing the charge carrier densities at two boundaries (n|x=0 p|x=d ). To do so, both increasing the generation rate by using concentrated incident light and efficient charge transport are feasible. Using high-density incident light source to generate high density charge carriers is already applied in concentrator solar cells that higher concentrated solar cells have higher Voc . 75–77 We also verified that with numerical simulations. 46 The other method is to improve charge transfer efficiency by either reducing recombination or increasing carriers lifetime and mobility. Numerical simulations showed that longer lifetime is, higher Voc outputs. 46 We also demonstrated that the presence of 9

ACS Paragon Plus Environment

The Journal of Physical Chemistry

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 27

interface recombination lowers Voc , 46 which was verified by recent publications that Voc is primarily limited by the interface losses. 44,45 Here, we found that besides the light intensity, charge carries lifetime and mobility and recombination, there is another parameter controls Voc .

Nc and Nv from DOS(E) The last approach to improve Voc is using materials with low Nc and Nv . Here, we calculated DOS(E) of typical perovskites, CdTe and silicon by VASP. DOS(E) are displayed in Figure 1. Nc and Nv were then calculated by integrating these DOS(E), as illustrated by Equation (2). Table 1 shows Nc and Nv of these materials. It was found that all the perovskites have similar Nc and Nv , with product around 6 × 1036 cm−6 , which is about one third of silicons. CdTe has a very low Nc of 3.0×1016 cm−3 as a result of the very sharp DOS(E) at bottom of conduction (Figure 1D). While, its Nv is the largest. According to Equation (9), only the product of Nc and Nv influences Voc as the output potential is the relative potential difference between quasi hole and electron Fermi levels. Among them, Si has the highest Nc and Nv product, which is up to 1037 cm−6 . The lower Nc and Nv in perovskite and CdTe are supposed to have higher Voc than silicon under same parameters and conditions. The Nc and Nv calculated from DOS(E) integration is very sensitive to calculation methods, although we used recommended settings for accurate electronic DOS(E), the integration of discrete DOS(E) brings about errors as well. The accuracy of Nc and Nv calculated from DOS(E) integral here is poor. Only the relative magnitude of the calculated Nc and Nv are reliable. Table 1: Nc and Nv calculated from integration of conduction and valence band DOS(E) of unit cells through Equation (2). Compound α-FAPbI3 MAPbCl3 MAPbI3 CdTe Si

V (˚ A3 ) 765.85 716.61 984.66 271.72 40.05

Nc -integral 4.55×10−3 8.56×10−3 1.03×10−2 8.19×10−6 4.56×10−4

Nv -integral 8.99×10−4 3.11×10−4 6.26×10−4 5.53×10−3 6.17×10−5

10

Nc (cm−3 ) 5.94×1018 1.19×1019 1.04×1019 3.00×1016 1.14×1019

ACS Paragon Plus Environment

Nc (cm−3 ) 1.17×1018 4.34×1017 6.36×1017 2.03×1019 1.54×1018

Nc Nv (cm−6 ) 6.97×1036 5.18×1036 6.62×1036 6.09×1035 1.75×1037

Page 11 of 27

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

A)

B)

D)

E)

C)

Figure 1: Density of states of A), α-FAPbI3 ; B), MAPbCl3 ; C), MAPbI3 ; D), CdTe and E), Si respectively.

Nc and Nv from effective mass Table 2: Nc , Nv and simulated solar cell performance parameters for α-FAPbI3 , MAPbCl3 , MAPbI3 , CdTe and silicon. Compound α-FAPbI3 MAPbCl3 MAPbI3 CdTe Si Si 78

Nv (cm−3 ) 2.90×1018 3.64×1018 2.49×1018 4.54×1019 2.32×1019 2.85×1019

Nc (cm−3 ) 1.20×1019 1.36×1019 6.98×1018 4.25×1017 3.35×1019 1.85×1019

Nc N v 3.47×1037 4.94×1037 1.74×1037 1.93×1037 7.76×1038 5.27×1038

Voc (mV) 1132.02 1123.32 1149.14 1146.46 1055.70 1065.21

PCEs (%) 19.70 19.53 20.04 19.99 18.18 18.37

Nc and Nv calculated from effective mass are expected to be more accurate. The error comes from the evaluation of effective mass. Here, we calculated both Nc and Nv from effective electron mass and effective hole mass respectively. The effective mass is calculated according to its definition and estimated by fitting band curves near conduction band bottoms and valence band tops. Details about the calculation method can be found in Ref. 11

ACS Paragon Plus Environment

The Journal of Physical Chemistry

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

79 and 23. Related band structures and fittings are available in SI. The calculated reduced effective mass of MAPbI3 for conductivity with the same method is 0.09 m0 , 23 similar to experimental values ranging from 0.09-0.15. (0.09-0.12 m0 , 80 0.12 m0 , 81 0.11 m0 , 82 0.10 m0 83 and 0.15 m0 . 84 ) This agreement confirms that our calculations are reliable. The estimated silicon hole effective masses are 0.12-0.76 m0 , and effective electron mass is 1.20 m0 . Both the calculated Nc and Nv are comparable to reference data. 78 Effective electron masses of the lowest conduction bands are 0.31-0.70 m0 and 0.17-0.53 m0 , for MAPbCl3 and MAPbI3 . Both of them are similar to the SOC-GW calculated 0.28-0.40 m0 and 0.18-0.30 m0 . 85 These comparisons suggest that our method is reliable. According to Equation (9), the product of Nc and Nv determines the voltage deficit from DOS. Therefore, Nc , Nv and their products are calculated in Table 2. The product sequence is Si > MAPbCl3 ≈ α-FAPbI3 ≈ MAPbI3 ≈ CdTe, which is in concert with the result calculated from integration. The absolute value difference between the Nc and Nv from effective mass and the Nc and Nv from the integral is due to the poor accuracy of Nc and Nv from DOS(E) integral. Although CdTe has a very small Nc , its Nv is large as it has three flat valence bands. The Nc and Nv products of CdTe and perovskites are in the same scale, which is much smaller than that of silicon. Considering Equation (9), the voltage deficits of perovskites should be about 100 mV smaller than silicon’s due to the products of hybrid perovskites are about 45 times lower than that of silicon.

Numerical simulations In order to figure out how the Nc and Nv influences solar cell performance and eliminate other parameters influences, such as band gap and carrier mobility influences, we assumed all of these materials have same parameters except the Nc and Nv . Several other parameters are correlated with Nc and Nv , such as effective mass and absorption. A low Nc or Nv is related to the small effective electron or hole mass, which also contributes to the high performance. As the small effective mass produces high charge carrier mobility and small 12

ACS Paragon Plus Environment

Page 12 of 27

Page 13 of 27

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

Figure 2: J-V curves of various materials with same electronic parameters except for Nc and Nv . The thickness of simulated solar cell is 350 nm. Diffusion coefficients of the perovskite were assumed to be 0.017 cm2 s−1 and 0.011 cm2 s−1 . 11 The band gap is set as 1.55 eV. 8 The interface recombination width is 5 nm and its lifetime is 5 ns. The charge carrier lifetime is assumed to be 736 ns. 8 Cells operate at 300 K and 1 sun (1.5AM). More simulated details can be found in Ref. 46. exciton separation energy, 23 mobilities and IPCEs are assumed the same in order to avoid the influence from the small effective mass. The low Nc and Nv may lead to low absorption. The absorption influence was excluded as we use same absorption coefficients. Using the numerical model we have developed for perovskite solar cells, 46 we simulated the working J-V curves of these materials. All of the simulated solar cells are ‘perovskite’ solar cells. For the simulated silicon and CdTe solar cells, they are not real silicon solar cells or CdTe solar cells, they are perovskite solar cells just with the Nc and Nv of silicon and CdTe. Figure 2 shows performances of their solar cells with different Nc and Nv . For the solar cell with silicon Nc and Nv , its Nc × Nv is 7.76×1038 cm−3 , the corresponding simulated Voc is 1055.70 mV. For the MAPbI3 , its Voc is calculated to be 1149.14 mV, which is about 93 mV higher than that of the silicon solar cell. If we lower the silicon Nc and Nv to the Nc and Nv of MAPbI3 , its PCEs will be improved from 18.18 % to 20.04 %. The relative improvement is 13

ACS Paragon Plus Environment

The Journal of Physical Chemistry

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

about 10 %.

Discussion Since the effective mass is calculated from the second derivative of bands, and number of K point in the band structure calculation are limited, this second derivative is not accurate. The error in effective mass evaluation is approximately 10 % in our calculations. , and the error for the reduced mass (m∗ ) is estimated to be approximately 15 %. The error in Nc or Nv is therefore approximately 23 %. The error in these two methods is obvious, but the sequence of Nc and Nv magnitude calculated from these two methods are the same. This confirms the low Nc and Nv in perovskite. An important question is does low Nc and Nv lead to low charge density near electrodes (n|x=0 p|x=d )? We find there will be no enhancement if the changing of n|x=0 p|x=d is the same as the changing of Nc ×Nv . Charge carriers are generated and then migrate to the electrodes. Hence, charge carrier densities at the two ends should increase with generation rate and mobilities, and decrease with thickness and recombination rate. The DOS has no effect on the recombination rate as the recombination rate is mostly dependent on the charge carrier density itself, defect density and defect energy levels according to the Shockley-Read-Hall Model. 68,86–88 Low Nc and Nv are related with small effective masses according to Equation 7, and produce high mobilities. The mobility effect was excluded in our simulations by using same mobilities. If take this into consideration, Voc will be further improved. The generation rate and thickness greatly influence each other. Theoretically, the genRdRλ RdRλ eration rate is calculated from G = 0 0 0 G(λ, x)dλ = 0 0 0 ηg [1 − R(λ)] × I(λ) × α(λ) × Rλ e−α(λ)x dλdx = 0 0 I(λ) × IP CE(λ, d)dλ, where α is the absorption coefficient and I(λ) is

incident photon density; R(λ) is the reflect coefficient and ηg is the generation quantum effi-

ciency; IP CE(λ, d) is the quantum efficiency of a whole device, which is also called External Quantum Efficiency (EQE). 89 In some experiments, solar cells show IP CE(λ) of approxi-

14

ACS Paragon Plus Environment

Page 14 of 27

Page 15 of 27

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

mately 80% in the visible light range. 6,12,90–93 For a small α material, a thicker thin film is required to reach the same IP CE. The low Nc and Nv in perovskite might lead to low absorption, which will lower solar cell performance. As shown in Ref. 94, saturable absorption happens at extreme incident light intensities. The intensity to exhibit saturable absorption for Bi2 Se3 is about 41-53 MW/cm2 , 95,96 which is about 108 times strong than 100mW/cm2 (the light intensity of AM1.5). There is no experimental saturable absorption data for perovskite materials, but the intensity should be of similar magnitude. Hence, the Nc and Nv ’s influence on absorption of sunlight is negligible for normal solar cells. Furthermore, the absorption of MAPbI3 is higher than silicon as shown in Ref. 97. As such, to reach the same IP CE, perovskite solar cell could be thinner. This thin film also leads to high n|x=0 p|x=d . In summary, if we consider the Nc and Nv related parameter influences, higher n|x=0 p|x=d are expected. The high Voc calculated from numerical models will be further improved.

Conclusion In conclusion, using two different methods based on ab initio results, we estimated Nc and Nv of silicon, CdTe and typical perovskites. It was found that the DOSs of perovskites and CdTe are much lower than that of silicon. The low Nc and Nv consist with the small effective mass, which directly contributes to high mobility. To avoid the small effective mass induced improvement, the mobility in these simulations is fixed. Analytic and numerical results show that the material with lower Nc and Nv can realize higher open circuit voltage and higher PCEs. With systematic discussions, we determined how the Nc and Nv of material works on its solar cell performance. The lower Nc and Nv of perovskites lead to about 100 mV higher Voc and 10 % (relative) higher PCEs. Considering low Nc and Nv consist with small effective masses, the performance of solar cells with low Nc and Nv will be further improved. In this paper, we show that besides inert interface defects, 44,46 long charge carrier lifetime, 46 and reduced bimolecular recombination, 45 that low Nc and Nv in perovskites contribute to the

15

ACS Paragon Plus Environment

The Journal of Physical Chemistry

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 27

high Voc and high performance of perovskite solar cells. Therefore, the materials with low Nc and Nv should be good candidates for solar cells.

Acknowledgments Y. Zhou is grateful to the National Computational Infrastructure (NCI, Australia) for High Performance Computing resources. Y. Zhou thanks Dr. Angus Gray-Weale for discussion.

ASSOCIATED CONTENT Supporting Information The Supporting Information is available free of charge on the ACS Publications website at DOI: 10.1021/acs.jpcc.00000. Method choosing, calculation details of effective masses and DOSs of silicon, FAPbI3 , MAPbCl3 , MAPbI3 , and CdTe.

References (1) Sites, J.; Pan, J. Strategies to Increase CdTe Solar-Cell Voltage. Thin Solid Films 2007, 515, 6099–6102. (2) Chuang, C.-H. M.; Maurano, A.; Brandt, R. E.; Hwang, G. W.; Jean, J.; Buonassisi, T.; Bulovi´c, V.; Bawendi, M. G. Open-Circuit Voltage Deficit, Radiative Sub-Bandgap States, and Prospects in Quantum Dot Solar Cells. Nano Lett. 2015, 15, 3286–3294. (3) Bourdais, S.; Chon´e, C.; Delatouche, B.; Jacob, A.; Larramona, G.; Moisan, C.; Lafond, A.; Donatini, F.; Rey, G.; Siebentritt, S. et al. Is the Cu/Zn Disorder the Main Culprit for the Voltage Deficit in Kesterite Solar Cells? Adv. Energy Mater. 2016, 6, 1502276. 16

ACS Paragon Plus Environment

Page 17 of 27

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

(4) Docampo, P.; Ball, J. M.; Darwich, M.; Eperon, G. E.; Snaith, H. J. Efficient Organometal Trihalide Perovskite Planar-Heterojunction Solar Cells on Flexible Polymer Substrates. Nat. Comm. 2013, 4, 2761. (5) Hodes, G. Applied Physics. Perovskite-Based Solar Cells. Science 2013, 342, 317–8. (6) 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, 643–7. (7) Liu, M.; Johnston, M. B.; Snaith, H. J. Efficient Planar Heterojunction Perovskite Solar Cells by Vapour Deposition. Nature 2013, 501, 395–8. (8) Zhou, H.; Chen, Q.; Li, G.; Luo, S.; Song, T.-B.; Duan, H.-S.; Hong, Z.; You, J.; Liu, Y.; Yang, Y. Interface Engineering of Highly Efficient Perovskite Solar Cells. Science 2014, 345, 542–546. (9) National Renewable Energy Laboratory (NREL), Research Cell Efficiency Records. http://www.nrel.gov/ncpv/images/efficiency_chart.jpg, Accessed: 2016-08-21. (10) Xing, G.; Mathews, N.; Sun, S.; Lim, S. S.; Lam, Y. M.; Gr¨atzel, M.; Mhaisalkar, S.; Sum, T. C.; Gratzel, M. Long-Range Balanced Electron- and Hole-Transport Lengths in Organic-Inorganic CH3 NH3 PbI3 . Science 2013, 342, 344–347. (11) Stranks, S. D.; Eperon, G. E.; Grancini, G.; Menelaou, C.; Alcocer, M. J. P.; Leijtens, T.; Herz, L. M.; Petrozza, A.; Snaith, H. J. Electron-Hole Diffusion Lengths Exceeding 1 Micrometer in an Organometal Trihalide Perovskite Absorber. Science 2013, 342, 341–4. (12) Burschka, J.; Pellet, N.; Moon, S.-J.; Humphry-Baker, R.; Gao, P.; Nazeeruddin, M. K.; Gr¨atzel, M. Sequential Deposition as a Route to High-Performance Perovskite-Sensitized Solar Cells. Nature 2013, 499, 316–9. 17

ACS Paragon Plus Environment

The Journal of Physical Chemistry

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

(13) Kim, H.-S.; Lee, C.-R.; Im, J.-H.; Lee, K.-B.; Moehl, T.; Marchioro, A.; Moon, S.-J.; Humphry-Baker, R.; Yum, J.-H.; Moser, J. E. et al. Lead Iodide Perovskite Sensitized All-Solid-State Submicron Thin Film Mesoscopic Solar Cell with Efficiency Exceeding 9%. Sci. Rep. 2012, 2, 591. (14) Ku, Z.; Rong, Y.; Xu, M.; Liu, T.; Han, H. Full Printable Processed Mesoscopic CH3 NH3 PbI3 /TiO2 Heterojunction Solar Cells with Carbon Counter Electrode. Sci. Rep. 2013, 3, 3132. (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) Brivio, F.; Walker, A. B.; Walsh, A. Structural and Electronic Properties of Hybrid Perovskites for High-Efficiency Thin-Film Photovoltaics from First-Principles. APL Materials 2013, 1, 042111. (17) Even, J.; Pedesseau, L.; Jancu, J.-M.; Katan, C. Importance of Spin-Orbit Coupling in Hybrid Organic/Inorganic Perovskites for Photovoltaic Applications. J. Phys. Chem. Lett. 2013, 4, 2999–3005. (18) Umari, P.; Mosconi, E.; De Angelis, F. Relativistic GW calculations on CH3 NH3 PbI3 and CH3 NH3 SnI3 Perovskites for Solar Cell Applications. Sci. Rep. 2014, 4, 4467. (19) Wang, Y.; Gould, T.; Dobson, J. F.; Zhang, H.; Yang, H.; Yao, X.; Zhao, H. Density Functional Theory Analysis of Structural and Electronic Properties of Orthorhombic Perovskite CH3 NH3 PbI3 . Phys. Chem. Chem. Phys. 2013, 16, 1424–9. (20) Snaith, H. J. Perovskites: The Emergence of a New Era for Low-Cost, High-Efficiency Solar Cells. J. Phys. Chem. Lett. 2013, 4, 3623–3630. (21) Baikie, T.; Fang, Y.; Kadro, J. M.; Schreyer, M.; Wei, F.; Mhaisalkar, S. G.; Graetzel, M.; White, T. J. Synthesis and Crystal Chemistry of the Hybrid Perovskite 18

ACS Paragon Plus Environment

Page 18 of 27

Page 19 of 27

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

(CH3 NH3 )PbI3 for Solid-State Sensitised Solar Cell Applications. J. Mater. Chem. A 2013, 1, 5628. (22) Du, M. H. Efficient Carrier Transport in Halide Perovskites: Theoretical Perspectives. J. Mater. Chem. A 2014, 2, 9091. (23) Zhou, Y.; Huang, F.; Cheng, Y.-B.; Gray-Weale, A. Photovoltaic Performance and the Energy Landscape of CH3 NH3 PbI3 . Phys. Chem. Chem. Phys. 2015, 17, 22604–22615. (24) Frost, J. M.; Butler, K. T.; Walsh, A. Molecular Ferroelectric Contributions to Anomalous Hysteresis in Hybrid Perovskite Solar Cells. APL Mater. 2014, 2, 081506. (25) 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–90. (26) Quarti, C.; Mosconi, E.; De Angelis, F. Interplay of Orientational Order and Electronic Structure in Methylammonium Lead Iodide: Implications for Solar Cells Operation. Chem. Mater. 2014, 26, 6557–6569. (27) Stroppa, A.; Quarti, C.; De Angelis, F.; Picozzi, S. Ferroelectric Polarization of CH3 NH3 PbI3 : A Detailed Study Based on Density Functional Theory and Symmetry Mode Analysis. J. Phys. Chem. Lett. 2015, 6, 2223–2231. (28) Green, M. A.; Emery, K.; Hishikawa, Y.; Warta, W.; Dunlop, E. D. Solar Cell Efficiency Tables (version 42). Prog. Photovoltaics Res. Appl. 2013, 21, 827–837. (29) Fukutani, K.; Kanbe, M.; Futako, W.; Kaplan, B.; Kamiya, T.; Fortmann, C.; Shimizu, I. Band Gap Tuning of a-Si:H from 1.55 eV to 2.10 eV by Intentionally Promoting Structural Relaxation. J. Non. Cryst. Solids 1998, 227-230, 63–67. (30) Liang, J.; Schiff, E. A.; Guha, S.; Yan, B.; Yang, J. Light-Soaking Effects on the OpenCircuit Voltage of a-Si:H Solar Cells. MRS Proceedings 2005, 862, A13.6. 19

ACS Paragon Plus Environment

The Journal of Physical Chemistry

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

(31) Green, M. A.; Emery, K.; Hishikawa, Y.; Warta, W.; Dunlop, E. D. Solar Cell Efficiency Tables (version 47). Prog. Photovoltaics Res. Appl. 2016, 24, 3–11. (32) Kuciauskas, D.; Dippo, P.; Kanevce, A.; Zhao, Z.; Cheng, L.; Los, A.; Gloeckler, M.; Metzger, W. K. The Impact of Cu on Recombination in High Voltage CdTe Solar Cells. Appl. Phys. Lett. 2015, 107, 0–5. (33) Zhao, Y.; Boccard, M.; Liu, S.; Becker, J.; Zhao, X.-H.; Campbell, C. M.; Suarez, E.; Lassise, M. B.; Holman, Z.; Zhang, Y.-H. Monocrystalline CdTe Solar Cells with OpenCircuit Voltage over 1 V and Efficiency of 17%. Nat. Energy 2016, 16067. (34) Burst, J. M.; Duenow, J. N.; Albin, D. S.; Colegrove, E.; Reese, M. O.; Aguiar, J. A.; Jiang, C.-S.; Patel, M.; Al-Jassim, M. M.; Kuciauskas, D. et al. CdTe Solar Cells with Open-Circuit Voltage Breaking the 1 V Barrier. Nat. Energy 2016, 1, 16015. (35) Fonthal, G.; Tirado-Meja, L.; Marn-Hurtado, J.; Ariza-Calder´on, H.; MendozaAlvarez, J. Temperature Dependence of the Band Gap Energy of Crystalline CdTe. J. Phys. Chem. Solids 2000, 61, 579–583. (36) Garcia-Belmonte, G.; Boix, P. P.; Bisquert, J.; Lenes, M.; Bolink, H. J.; La Rosa, A.; Filippone, S.; Martin, N. Influence of the Intermediate Density-of-States Occupancy on Open-Circuit Voltage of Bulk Heterojunction Solar Cells with Different Fullerene Acceptors. J. Phys. Chem. Lett. 2010, 1, 2566–2571. (37) Sulas, D. B.; Yao, K.; Intemann, J. J.; Williams, S. T.; Li, C.-Z.; Chueh, C.-C.; Richards, J. J.; Xi, Y.; Pozzo, L. D.; Schlenker, C. W. et al. Open-Circuit Voltage Losses in Selenium-Substituted Organic Photovoltaic Devices from Increased Density of Charge-Transfer States. Chem. Mater. 2015, 27, 6583–6591. (38) Elumalai, N. K.; Uddin, A. Open Circuit Voltage of Organic Solar Cells: an in-Depth Review. Energy Environ. Sci. 2016, 9, 391–410.

20

ACS Paragon Plus Environment

Page 20 of 27

Page 21 of 27

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

(39) Collins, S. D.; Proctor, C. M.; Ran, N. A.; Nguyen, T.-Q. Understanding Open-Circuit Voltage Loss through the Density of States in Organic Bulk Heterojunction Solar Cells. Adv. Energy Mater. 2016, 6, 1501721. (40) Zhong, C.-L.; Geng, K.-W.; Luo, L.-E.; Yang, D.-W. An Analytical Model to Explore Open-Circuit Voltage of a-Si:H/c-Si Heterojunction Solar Cells. J. Cent. South Univ. 2016, 23, 598–603. (41) Blakesley, J. C.; Neher, D. Relationship between Energetic Disorder and Open-Circuit Voltage in Bulk Heterojunction Organic Solar Cells. Phys. Rev. B 2011, 84, 075210. (42) Shao, Y.; Yuan, Y.; Huang, J. Correlation of Energy Disorder and Open-Circuit Voltage in Hybrid Perovskite Solar Cells. Nature Energy 2016, 1, 15001. (43) Yan, W.; Li, Y.; Ye, S.; Li, Y.; Rao, H.; Liu, Z.; Wang, S.; Bian, Z.; Huang, C. Increasing Open Circuit Voltage by Adjusting Work Function of Hole-transporting Materials in Perovskite Solar Cells. Nano Research 2016, 9, 1600–1608. (44) Chen, S.; Hou, Y.; Chen, H.; Richter, M.; Guo, F.; Kahmann, S.; Tang, X.; Stubhan, T.; Zhang, H.; Li, N. et al. Exploring the Limiting Open-Circuit Voltage and the Voltage Loss Mechanism in Planar CH3 NH3 PbBr3 Perovskite Solar Cells. Adv. Energy Mater. 2016, 1600132. (45) Yang, W.; Yao, Y.; Wu, C.-Q. Origin of the High Open Circuit Voltage in Planar Heterojunction Perovskite Solar Cells: Role of the Reduced Bimolecular Recombination. J. Appl. Phys. 2015, 117, 095502. (46) Zhou, Y.; Gray-Weale, A. A Numerical Model for Charge Transport and Energy Conversion of Perovskite Solar Cells. Phys. Chem. Chem. Phys. 2016, 18, 4476–4486. (47) Kirchartz, T.; Staub, F.; Rau, U. Impact of Photon Recycling on the Open-Circuit Voltage of Metal Halide Perovskite Solar Cells. ACS Energy Letters 2016, 1, 731–739. 21

ACS Paragon Plus Environment

The Journal of Physical Chemistry

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

(48) Staub, F.; Hempel, H.; Hebig, J.-C.; Mock, J.; Paetzold, U. W.; Rau, U.; Unold, T.; Kirchartz, T. Beyond Bulk Lifetimes: Insights into Lead Halide Perovskite Films from Time-Resolved Photoluminescence. Physical Review Applied 2016, 6, 044017. (49) Endres, J.; Egger, D. A.; Kulbak, M.; Kerner, R. A.; Zhao, L.; Silver, S. H.; Hodes, G.; Rand, B. P.; Cahen, D.; Kronik, L. et al. Valence and Conduction Band Densities of States of Metal Halide Perovskites: A Combined ExperimentalTheoretical Study. J. Phys. Chem. Lett. 2016, 7, 2722–2729. (50) Kawai, H.; Giorgi, G.; Marini, A.; Yamashita, K. The Mechanism of Slow Hot-Hole Cooling in Lead-Iodide Perovskite: First-Principles Calculation on Carrier Lifetime from ElectronPhonon Interaction. Nano Letters 2015, 15, 3103–3108. (51) Lin, Q.; Armin, A.; Chandra, R.; Nagiri, R.; Burn, P. L.; Meredith, P. Electro-optics of Perovskite Solar Cells. Nat. Photonics 2014, 9, 106–112. (52) Xu, W.; Liu, L.; Yang, L.; Shen, P.; Sun, B.; McLeod, J. A. Dissociation of Methylammonium Cations in Hybrid OrganicInorganic Perovskite Solar Cells. Nano Lett. 2016, 16, 4720–4725. (53) Lee, J.-W.; Seol, D.-J.; Cho, A.-N.; Park, N.-G. High-Efficiency Perovskite Solar Cells Based on the Black Polymorph of HC(NH2 )2 PbI3 . Adv. Mater. 2014, 6, 4991–4998. (54) Zhou, Y.; Yang, M.; Kwun, J.; Game, O. S.; Zhao, Y.; Pang, S.; Padture, N. P.; Zhu, K. Intercalation Crystallization of Phase-Pure α-HC(NH2 )PbI3 Upon Microstructurally Engineered PbI2 Thin Films for Planar Perovskite Solar Cells. Nanoscale 2016, 8, 6265–6270. (55) Wang, Z.; Zhou, Y.; Pang, S.; Xiao, Z.; Zhang, J.; Chai, W.; Xu, H.; Liu, Z.; Padture, N. P.; Cui, G. Additive-Modulated Evolution of HC(NH2 )2 PbI3 Black Polymorph for Mesoscopic Perovskite Solar Cells. Chem. Mater. 2015, 27, 7149–7155.

22

ACS Paragon Plus Environment

Page 22 of 27

Page 23 of 27

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

(56) Song, J.; Hu, W.; Wang, X.-F.; Chen, G.; Tian, W.; Miyasaka, T. HC(NH2 )2 PbI3 as a Thermally Stable Absorber for Efficient ZnO-based Perovskite Solar Cells. J. Mater. Chem. A 8435–8443. (57) Kim, J.; Yun, J. S.; Wen, X.; Soufiani, A. M.; Lau, C. F. J.; Wilkinson, B.; Seidel, J.; Green, M. A.; Huang, S.; Ho-Baillie, A. W. Y. Nucleation and Growth Control of HC(NH2 )2 PbI3 for Planar Perovskite Solar Cell. J. Phys. Chem. C 2016, 120, 11262– 11267. (58) 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–38. (59) Kresse, G.; Furthm¨ uller, J. Efficiency of ab-initio Total Energy Calculations for Metals and Semiconductors Using a Plane-Wave Basis Set. Comput. Mater. Sci. 1996, 6, 15– 50. (60) Kresse, G.; Furthmuller, J.; Furthm¨ uller, J. Efficient Iterative Schemes for ab initio Total-Energy Calculations Using a Plane-Wave Basis Set. Phys. Rev. B 1996, 54, 11169–11186. (61) Perdew, J. P.; Ruzsinszky, A.; Csonka, G. I.; Vydrov, O. A.; Scuseria, G. E.; Constantin, L. A.; Zhou, X.; Burke, K. Restoring the Density-Gradient Expansion for Exchange in Solids and Surfaces. Phys. Rev. Lett. 2008, 100, 136406. (62) Aryasetiawan, F.; O. Gunnarsson, The GW Method. Reports Prog. Phys. 1998, 61, 237–312. (63) Hedin, L. New Method for Calculating the One-Particle Green’s Function with Application to the Electron-Gas Problem. Phys. Rev. 1965, 139, A796–A823.

23

ACS Paragon Plus Environment

The Journal of Physical Chemistry

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

(64) Peralta, J. E.; Heyd, J.; Scuseria, G. E.; Martin, R. L. Spin-orbit Splittings and Energy Band Gaps Calculated with the Heyd-Scuseria-Ernzerhof Screened Hybrid Functional. Phys. Rev. B 2006, 74, 073101. (65) Heyd, J.; Scuseria, G. E.; Ernzerhof, M. Hybrid Functionals Based on a Screened Coulomb Potential. J. Chem. Phys. 2003, 118, 8207. (66) Men´endez-Proupin, E.; Palacios, P.; Wahn´on, P.; Conesa, J. C. Self-Consistent Relativistic Band Structure of the CH3 NH3 PbI3 Perovskite. Phys. Rev. B 2014, 90, 045207. (67) R Core Team, R: A Language and Environment for Statistical Computing. R Foundation for Statistical Computing: Vienna, Austria, 2016. (68) Sheng, L. Semicond. Phys. Electron., 2nd ed.; Springer: Berlin, 2006. (69) Simon M. Sze, Physics of Semiconductor Devices, 2nd ed.; Wiley: New York, 1981. (70) Brivio, F.; Butler, K. T.; Walsh, A.; van Schilfgaarde, M. Relativistic Quasiparticle SelfConsistent Electronic Structure of Hybrid Halide Perovskite Photovoltaic Absorbers. Phys. Rev. B 2014, 89, 155204. (71) G´alvez, F. E.; Kemppainen, E.; M´ıguez, H.; Halme, J. Effect of Diffuse Light Scattering Designs on the Efficiency of Dye Solar Cells: An Integral Optical and Electrical Description. J. Phys. Chem. C 2012, 116, 11426–11433. (72) Soedergren, S.; Hagfeldt, A.; Olsson, J.; Lindquist, S.-e. Theoretical Models for the Action Spectrum and the Current-Voltage Characteristics of Microporous Semiconductor Films in Photoelectrochemical Cells. J. Phys. Chem. 1994, 98, 5552–5556. (73) Halme, J.; Vahermaa, P.; Miettunen, K.; Lund, P. Device Physics of Dye Solar Cells. Adv. Mater. 2010, 22, E210–34. (74) Sachenko, A. V.; Kryuchenko, Y. V.; Kostylyov, V. P.; Bobyl, A. V.; Terukov, E. I.; Abolmasov, S. N.; Abramov, A. S.; Andronikov, D. A.; Shvarts, M. Z.; Sokolovskyi, I. O. 24

ACS Paragon Plus Environment

Page 24 of 27

Page 25 of 27

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

et al. Temperature Dependence of Photoconversion Efficiency in Silicon Heterojunction Solar Cells: Theory vs Experiment. J. Appl. Phys. 2016, 119, 225702. (75) King, R. R.; Law, D. C.; Edmondson, K. M.; Fetzer, C. M.; Kinsey, G. S.; Yoon, H.; Sherif, R. A.; Karam, N. H. 40% Efficient Metamorphic GaInPGaInAsGe Multijunction Solar Cells. Appl. Phys. Lett. 2007, 90, 183516. (76) Sinton, R.; Young Kwark,; Gan, J.; Swanson, R. 27.5-Percent Silicon Concentrator Solar Cells. IEEE Electron Device Lett. 1986, 7, 567–569. (77) Bi, D.; Yang, L.; Boschloo, G.; Hagfeldt, A.; Johansson, E. M. J. Effect of Different Hole Transport Materials on Recombination in CH3 NH3 PbI3 Perovskite-Sensitized Mesoscopic Solar Cells. J. Phys. Chem. Lett. 2013, 4, 1532–1536. (78) Principles of Electronic Devices: Effective Mass in Semiconductors. http://ecee. colorado.edu/~bart/book/effmass.htm, Accessed: 2016-08-21. (79) Zhou, Y.-C.; Zhang, H.-L.; Deng, W.-Q. A 3N Rule for the Electronic Properties of Doped Graphene. Nanotechnology 2013, 24, 225705. (80) Galkowski, K.; Mitioglu, A.; Miyata, A.; Plochocka, P.; Portugall, O.; Eperon, G. E.; Wang, J. T.-W.; Stergiopoulos, T.; Stranks, S. D.; Snaith, H. J. et al. Determination of the Exciton Binding Energy and Effective Masses for Methylammonium and Formamidinium Lead Tri-halide Perovskite Semiconductors. Energy Environ. Sci. 2016, 9, 962–970. (81) Hirasawa, M.; Ishihara, T.; Goto, T.; Uchida, K.; Miura, N. Magnetoabsorption of the Lowest Exciton in Perovskite-Type Compound CH3 NH3 PbI3 . Phys. B Condens. Matter 1994, 201, 427–430. (82) Ishihara, T. Optical Properties of PbI-based Perovskite Structures. J. Lumin. 1994, 60-61, 269–274. 25

ACS Paragon Plus Environment

The Journal of Physical Chemistry

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

(83) Miyata, A.; Mitioglu, A.; Plochocka, P.; Portugall, O.; Wang, J. T.-W.; Stranks, S. D.; Snaith, H. J.; Nicholas, R. J. Direct Measurement of the Exciton Binding Energy and Effective Masses for Charge Carriers in OrganicInorganic tri-halide Perovskites. Nature Physics 2015, 11, 582–587. (84) Tanaka, K.; Takahashi, T.; Ban, T.; Kondo, T.; Uchida, K.; Miura, N. Comparative Study on the Excitons in Lead-Halide-Based Perovskite-Type Crystals CH3 NH3 PbBr3 CH3 NH3 PbI3 . Solid State Commun. 2003, 127, 619–623. (85) Mosconi, E.; Umari, P.; De Angelis, F. Electronic and Optical Properties of MAPbX3 Perovskites (X = I, Br, Cl): a Unified DFT and GW Theoretical Analysis. Phys. Chem. Chem. Phys. 2016, 18, 27158–27164. (86) Shockley, W.; Read, W. Statistics of the Recombinations of Holes and Electrons. Phys. Rev. 1952, 87, 835–842. (87) Hall, R. Electron-Hole Recombination in Germanium. Phys. Rev. 1952, 87, 387–387. (88) Altermatt, P. P. Models for Numerical Device Simulations of Crystalline Silicon Solar Cellsa Review. J. Comput. Electron. 2011, 10, 314–330. (89) Solanki, C. S. Solar Photovoltaics: Fundamentals, Technologies And Applications, 3rd ed.; PHI Learning, 2015; p 108. (90) Zhang, Q.; Kan, B.; Liu, F.; Long, G.; Wan, X.; Chen, X.; Zuo, Y.; Ni, W.; Zhang, H.; Li, M. et al. Small-Molecule Solar Cells with Efficiency over 9%. Nat. Photonics 2014, 9, 35–41. (91) Zhang, X.; Chen, L.; Li, X.; Mao, J.; Wu, W.; ˚ Agren, H.; Hua, J. Photovoltaic Properties of bis(octyloxy)benzo-[c][1,2,5]thiadiazole Sensitizers Based on an N,Ndiphenylthiophen-2-amine Donor. J. Mater. Chem. C 2014, 2, 4063.

26

ACS Paragon Plus Environment

Page 26 of 27

Page 27 of 27

The Journal of Physical Chemistry

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 ACS Paragon Plus Environment