Electronic structure and optical properties of Cs2AX2′X4 (A=Ge,Sn,Pb; X′,X=Cl,Br,I) Electronic structure and optical properties of (A=Ge,Sn,Pb; X′,X=Cl,Br,I) Guangtao Wang, Dongyang Wang, and Xianbiao[.]
Electronic structure and optical properties of X′,X=Cl,Br,I) Guangtao Wang, Dongyang Wang, and Xianbiao Shi Citation: AIP Advances 5, 127224 (2015); doi: 10.1063/1.4939016 View online: http://dx.doi.org/10.1063/1.4939016 View Table of Contents: http://aip.scitation.org/toc/adv/5/12 Published by the American Institute of Physics (A=Ge,Sn,Pb; AIP ADVANCES 5, 127224 (2015) Electronic structure and optical properties of Cs2AX′2X4 (A=Ge,Sn,Pb; X′,X=Cl,Br,I) Guangtao Wang,a Dongyang Wang, and Xianbiao Shi College of Physics and Information Engineering, Henan Normal University, Xinxiang, Henan 453007, People’s Republic of China (Received 25 August 2015; accepted December 2015; published online 21 December 2015) We studied the crystal structures, electronic structures and optical properties of Cs2AX2′ X4 (A=Ge,Sn,Pb; X′, X=Cl, Br, I) compounds using the first-principles calculation Our optimized structures agree well with experimental and theoretical results Band structure calculations, using the modified Becke-Johnson (mBJ) potential method, indicate that these compounds (with the exception of Cs2PbX2′ I4) are semiconductors with the direct band gap ranging from 0.36 to 4.09 eV We found the compounds Cs2GeBr2I4, Cs2GeCl2I4, Cs2GeI2Br4, Cs2SnI6, and Cs2SnBr2I4 may be good candidates for lead-free solar energy absorber materials C 2015 Author(s) All article content, except where otherwise noted, is licensed under a Creative Commons Attribution 3.0 Unported License [http://dx.doi.org/10.1063/1.4939016] I INTRODUCTION As the global demand for energy grows inexorably, solar energy is attracts an increasing amount of attention owing to its abundant, clean, and renewable characteristics The pursuit of high efficiency solar cells has not slowed since the first generation reported The first generation crystalline silicon based solar cells reached a high efficiency, but they were limited by a low optical absorption coefficient and high cost of fabrication.1–3 Cu-based selenide and sulfide compounds are considered to be a promising class of materials for solar cell applications, with the higher power conversion efficiency (PCE) to date, about 23.3%.4–8 Recently, the hybrid organic-inorganic perovskites metal halides have emerged as promising light-harvesting materials for the next generation of solar cells, with advantageous properties including direct band gaps, large absorption coefficients, and high carrier mobilities.9–13 Miyasaka et al.10 pioneered the incorporation of the hybrid organic-inorganic perovskite halides into solar cells, and found that CH3NH3PbI3 based solar cells had a 3.8% PCE, which was improved to 6.5% by Lee et al.14 A breakthrough occurred in 2012 with the realization of the long term stability of all-solid-state solar cells using this materials; these had a PCE of 9.7%.15 Chung et al.16 improved the PCE up to 10.2% in a new type of all-solid-state inorganic solar cell, which consisted of the p-type direct band gap semiconductor CsSnI3 and n-type TiO2 with dye N719 The work of Liu et al.11 indicates that the PCE could reach over 15% using an absorbing layer of the vapor-deposited perovskite CH3NH3PbI(3−x)Cl x A solution-processed solar cell using CH3NH3PbI3 as the absorption material could achieve a PCE of 20.1%.8 More recently, Lee et al.17 designed a new type of solar cell by replacing CsSnI3 with the more stable Cs2SnI6 as the absorbing layer and achieved a PCE of about 8% Unlike the tin-based perovskite compound CsSnI3, which has Sn2+ and must be handled in an inert atmosphere when fabricating the solar cells, Cs2SnI6 contains Sn4+, making it more stable under exposure to air and moisture This characteristic is vital for the realization of the low-cost, stable, lead-free, and environmentally friendly next generation of solid-states solar cells In addition to the experimental results as mentioned above, there are many first-principle calculations about various perovskite solar materials.18–22 While searching for lead-free, high efficiency photovoltaic materials, Lang et al.23 studied the optical properties of a Electronic mail: wangtao@henannu.edu.cn 2158-3226/2015/5(12)/127224/7 5, 127224-1 © Author(s) 2015 127224-2 Wang, Wang, and Shi AIP Advances 5, 127224 (2015) MAB3, where M=Cs, CH3NH3, NH2CHNH2; A=Sn, Pb; B=Cl, Br, I However, some of these MAB3 compounds are unstable in the air and moist environments In this paper, we report the crystal structures, electronic structures and absorption coefficients of Cs2AX2′ X4 (A=Ge,Sn,Pb; X′,X=Cl,Br,I) materials determined through the first-principles calculations We found Cs2GeBr2I4, Cs2GeCl2I4, Cs2GeI2Br4, Cs2SnI6, and Cs2SnBr2I4 may be good candidates for lead-free solar energy absorber materials These materials will provide more options in the pursuit of high PCE solar cell absorber materials II METHODOLOGY The structure optimization calculations are based on the density functional theory (DFT) method, as implemented in the Vienna Ab-initio Simulation Package (VASP).30 The generalized gradient approximation (GGA)31 formulated by Perdew-Burke-Ernzerhof (PBE)32 is employed for exchange-correlation functional The project-augmented wave (PAW) potential is employed to describe the electron-ion potential The wave functions are expanded in plane waves with a cut off energy of 500 eV The convergence criterion for the electronic self-consistent iteration and the Hellmann-Feynman forces are set to 10−5 eV and < 0.01 eV Å−1, respectively The Brillouin-zone integration is performed with 11 × 11 × 11 and 13 × 13 × 11 Γ-centered Monkhorst-Pack (MP)33 grids for cubic and tetragonal phase, respectively The lattice constants of the analyzed materials are listed in Table I Taking the optimized structures listed in Table I, we use the WIEN2k34 code with the PBE and modified Becke-Johnson (mBJ) potential35 to study the electronic structure and optical properties of the materials The mBJ method successfully in overcomes the problem of underestimating the band gap, since it considers the core electrons as fully relativistic and the valence electrons are treated semi-relativistically To achieve the energy eigenvalue convergence, the basis function are expanded up to R MT × Kmax = 7, where RMT is the minimum radius of muffin-tin spheres and Kmax is the magnitude of the largest K vector in the plane wave expansion The valence wave functions inside the spheres are expanded up to l max = 10, and the charge density is expanded as a Fourier series up to Gmax = 12 The self-consistent potentials are calculated on 1000 k-points in the first Brillouin zone, and the convergence criterion is set to 0.0001 Ry Denser (8000) k-points in the irreducible wedge of the first Brillouin zone are used for calculations of the optical properties TABLE I Optimized lattice constant (in Å) of Cs2AX′2X4 The available experimental results are shown in square brackets The predicted lattice constant from the formula in Ref 24 is shown in parentheses A Cs2AI6 Cs2ABr2I4 Cs2ACl2I4 Cs2AI2Br4 Cs2ABr6 Ge a=11.804 (a=11.317) a=8.153 c=11.792 a=8.163 c=11.538 a=8.162 c=10.958 a=11.081 (a=10.633) Sn a=12.032 [a=11.65]26 a=8.369 c=11.771 a=8.335 c=11.563 a=8.251 c=11.259 a=11.243 [a=10.770]27 Pb a=12.066 (a=11.553) a=8.210 c=12.246 a=8.146 c=12.069 a=8.378 c=11.089 a=11.327 (a=10.869) A Cs2ACl2Br4 Cs2AI2Cl4 Cs2ABr2Cl4 Cs2ACl6 Ge a=7.687 c=11.037 a=8.047 c=10.331 a=7.678 c=10.536 a=10.560 [a=10.23]25 Sn a=7.805 c=11.196 a=8.125 c=11.563 a=7.787 c=10.714 a=10.724 [a=10.355]28 Pb a=7.825 c=11.376 a=8.269 c=10.422 a=7.871 c=10.709 a=10.821 [a=10.416]29 127224-3 Wang, Wang, and Shi AIP Advances 5, 127224 (2015) FIG Crystal structure of (a) the cubic phase Cs2AX6 and (b) the tetragonal phase Cs2AX′2X4 III RESULTS AND DISCUSSION ¯ (No.225) space group, as Cs2AX6 compounds crystallize in a cubic structure with the Fm3m shown in Fig 1(a) In the case of Cs2AX2′ X4, shown in Fig 1(b), the X ion at the apex of the octahedron along the z-axis is replaced by the other halogen ion X′ The symmetry of this new crystal is lowered to I4/mmm (No.139) The optimized lattice constants of Cs2AX2′ X4 for cubic and tetragonal phase are presented in Table I The results are in good agreement with the available experiment data25–29 and the theoretical prediction reported by Brik et al.24 The lattice constants and bond length in the octahedron of these crystals are shown in Fig 2(a), 2(b), 2(c) and Fig 2(d), 2(e), 2(f) The changes in lattice parameters and bond length can be ascribed to radius changes among A and the halogen ions When X′ and X are fixed, the lattice parameter a of Cs2AX2′ X4 in the tetragonal phase decreases with the radius of the A ion from Pb, Sn, to Ge, with a corresponding increase of the lattice parameter c This leads to a distortion of the octahedron and lowers the symmetry of the crystal, as shown in Fig When the A and X ions are fixed, the bond length of A-X′ along the z axis decreases from I, Br, to Cl FIG The lattice constant (solid-line for a, dashed-line for c, Å), bond length (BL in Å, solid-line for BL in x-y plane, dashed-line for BL along c-axis), band gap (E g in eV), and octahedron ionicity (Fo ) of Cs2AX′2X4 are presented as a function of X′ atoms type in the (a,b,c), (d,e,f), (g,h,i) and (j,k,l), respectively The color line in √ black, red, and green for A = Ge, Sn, and Pb, respectively For the tetragonal phase, the lattice constant was defined to be a = 2a 0, where a is the lattice constant of a conventional cell and a is the super-cell lattice constant 127224-4 Wang, Wang, and Shi AIP Advances 5, 127224 (2015) FIG (a) Band structure with mBJ (blue lines) and PBE (red lines) methods and (b) projected density of state with mBJ method in cubic phase Cs2PbBr6 The zero energy is at the VBM We calculate the projected density of states (PDOS) and the band structure of Cs2AX2′ X4 with mBJ and PBE methods, and summarize their band gaps with mBJ method in Fig 2(g), 2(h), 2(i) Our results indicate that Cs2PbX2′ I4 (X′ = Cl,Br,I) show metallic characteristic, and all the other compounds behave as semiconductors with a band gap (Eg ) ranging from 0.36 to 4.09 eV, as shown in Fig 2(g), 2(h), 2(i) We consider the example of Cs2PbBr6 for further investigation; its band structure and PDOS are shown in Fig The mBJ method significantly improves the band gap and overcomes the underestimation of PBE method In Fig 3(a), we note the direct band gap of 1.26 eV located at Γ(0, 0, 0) point This Eg is smaller than the direct band gap of CH3NH3PbBr3 (2.2 eV), located at R(0.5, 0.5, 0.5).36 The conduction band maximum (CBM) only consists of one band derived from the hybridization of Pb-6s and Br-4p states, while the valence band minimum (VBM) is derived from the Br-5p states This is significantly different from the CBM of CsPbBr3, which is threefold degenerated at the R point and derived from the Pb-6p state.23 For CH3NH3PbBr3, the calculated band gap with spin-orbital coupling (SOC) is 1.0 eV smaller than that without SOC, but the SOC effect in the Cs2PbBr6 has little impact on the Eg (a change of ∼ 0.2 eV) This is because that the CBM in Cs2PbBr6 is derived from the hybridization of Pb-6s and Br-4p states, while the CBM of CsPbBr3 is derived from the Pb-6p state and the SOC have a significant impact for this compound All the tetragonal phase compounds have a similar electronic structure, except for a difference in band gap, as shown in Fig 2(g), 2(h), 2(i) Therefore, we show results of one representative compound Cs2PbI2Cl4 The band structure and PDOS are presented in Fig The band structure in Fig 4(a) also shows that the mBJ method can effectively overcome the underestimate of band gap with PBE method We note that both the CBM and VBM lie at the Γ(0, 0, 0) point, with a direct band gap Eg = 0.96 eV The CBM is derived from the mixture of Pb-6s, I-5p, and Cl-3p states The VBM is mainly derived from the I-5p and Cl-3p states; the band center of I-5p is higher than the Cl-3p states, as shown in Fig 4(b) Since the Pb-6p state has no contribution to the CBM, the SOC effect does not need to be considered This is also provides a evidence for the existence of Pb4+ in the tetragonal phase The band gaps of all the compounds are plotted in Fig 2(g), 2(h), 2(i) The trend of Eg can be understood by considering the Pauling electronegativity ( χ A) of the A ion and the octahedron ionicity (Fo ) We defined the octahedron ionicity (Fo ) as Fo = (2 × FAX ′ + × FAX )/6, where the FAX = − e−(χ A−χ X ) /4 is the Pauling ionicity, and χ A and χ X refer to the electronegativity of A and X atoms When X′ and X are fixed, the band gap Eg increases with the decreasing of χ A, from Pb, 127224-5 Wang, Wang, and Shi AIP Advances 5, 127224 (2015) FIG (a) Band structure with mBJ (blue lines) and PBE (red lines) methods and (b) projected density of state with mBJ method in tetragonal phase Cs2PbI2Cl4 The zero energy is at the VBM Ge, to Sn ( χ Ge = 2.01, χ Sn = 1.96, χ Pb = 2.33) The octahedron ionicity Fo indicates the electron transfers from X (and/or X′) to A in the octahedron In Fig 2(g), 2(h), 2(i), it shows that the larger Fo , the larger band gap The parameter Fo can help us to search for new compounds with suitable band gap, before doing first-principle calculations Semiconductors with Eg between 1.0 eV and 1.5 eV are generally well accepted to have the greatest potential in the application of high efficiency solar cells With this in mind, we predict that Cs2GeBr2I4, Cs2GeCl2I4, Cs2GeI2Br4, Cs2SnI6, Cs2SnBr2I4, Cs2SnCl2I4, Cs2PbBr6, Cs2PbI2Cl4, Cs2PbCl2Br4, and Cs2PbBr2Cl4 will be more promising candidates for the solar cell absorbers Now, we focus on these materials and compare their optical properties In Fig 5, we present the absorption coefficient of Cs2AX2′ X4 in the range of 300 to 800 nm It is seen from figure that the optical absorption coefficient of tetragonal phase exhibit strong anisotropy with respect to the relative orientation of the polarization of light and crystal lattice This polarization anisotropy results from anisotropy of the tetragonal phase dielectric tensor ϵ(ω) We denote the two components of the electric field (E) relative to the optical axis The dielectric component for E⊥c is ϵ ⊥ = ϵ xx = ϵ yy and that for E∥c is ϵ ∥ = ϵ zz, So we label the optical absorption coefficient for the light polarization in the ab plane as α⊥ and along c axis as α ∥ for the tetragonal phase The absorption coefficient α⊥ of Cs2GeCl2I4 has a high and sharp peak around 400 nm (Fig 5(a)), indicating a strong absorption of purple light The absorption coefficients α⊥ of Cs2GeBr2I4, Cs2GeCl2I4, Cs2SnBr2I4, Cs2SnCl2I4, Cs2PbI2Cl4, and Cs2PbCl2Br4 are larger than those α ∥, while the cases are reversed for Cs2GeI2Br4 and Cs2PbI2Cl4 The absorption coefficients of the Pb-based compounds are larger than those of the Ge- and Sn-based compounds In Fig 6, we define the Ideal Power Absorption Coefficient IPAC = I0(E) ∗ a(E), where the I0(E) is the standard AM 1.5G solar spectrum39 and a(E) is the photon absorptivity The absorptivity is given by a(E) = − e−2α(E)×L , which depends on the absorption coefficient α(E) and absorber thickness L.40,41 Here we used typical absorber thickness of L = 500 nm, because that such absorber thickness is often used in the experiments.13,17,42 For α(E), it is derived from weighted average of absorption coefficient as α(E) = (α⊥ + 2α ∥)/3 The area enclosed by the IPAC curves and the abscissa axis is equal to the absorption power of the materials We integrate the IPAC curves from 300 to 800 nm and present the resulting absorption power values in Fig In Fig 6(a), the absorption power of Cs2GeBr2I4 (718.21 W/m2) is larger than the other Ge-based compounds In Fig 6(b), the absorption power of Cs2SnCl2I4 (719.78 W/m2) is higher than the other Sn-based compounds Both of them are compatible to the largest value of 748.79 W/m2 derived from Cs2PbI2Cl4 The Pb-based compounds 127224-6 Wang, Wang, and Shi AIP Advances 5, 127224 (2015) FIG Absorption coefficient of Cs2AX′2X4, where A ion is (a) Ge, (b) Sn, and (c) Pb, respectively The optical absorption coefficient37,38 for the light polarization in the a-b plane (solid line) note as α ⊥ and along c axis (dashed line) note as α ∥ for the tetragonal phase, respectively Cs2PbBr6 and Cs2PbI2Cl4 have larger absorption powers than the Ge-based and Sn-based compounds, but they are not environmentally friendly The environmentally friendly compounds Cs2GeBr2I4, Cs2GeCl2I4, Cs2GeI2Br4, Cs2SnI6, and Cs2SnBr2I4 may be good candidates for solar energy absorber materials FIG Ideal Power Absorption Coefficient (IPAC) of Cs2AX′2X4 under the standard AM 1.5G solar spectrum,39 as a function of the wavelength Panels (a), (b), and (c) correspond to A=Ge, Sn and Pb, respectively IPAC = I0[1 − e (2L×(2α ⊥+α ∥)/3], where the I0 is the standard AM 1.5G solar spectrum, α ⊥ and α ∥ are the absorption coefficient and the L is the absorber thickness.40,41 An arrow points from each curve to a number that represents the result of the integration of that curve; this corresponds to the absorption power (in W/m2) 127224-7 Wang, Wang, and Shi AIP Advances 5, 127224 (2015) IV SUMMARY In this work, we have studied the electronic structure and optical properties of Cs2AX2′ X4 (A=Ge, Sn, Pb; X′, X=Cl, Br, I) compounds using first-principles calculations The optimized crystal structures we determined are in good agreement with the experimental data We found that the trends of band gaps are consistent with the octahedron ionicity The Pb-based compounds have larger absorption powers than the Ge-based and Sn-based compounds, but they are not environmentally friendly We found that Cs2GeBr2I4, Cs2GeCl2I4, Cs2GeI2Br4, Cs2SnI6, and Cs2SnBr2I4 might be environmentally friendly solar energy absorber materials ACKNOWLEDGMENTS The authors acknowledge support from the NSF of China (No.11274095, No.10947001) and the Program for Science and Technology Innovation Talents in the Universities of Henan Province (No.2012HASTIT009, No.104200510014, and No.114100510021) This work is supported by The High Performance Computing Center of Henan Normal University J Nelson, The physics of solar cells (Imperial College Press, London, 2003) V Fthenakis et al., Renew Sust Energy Rev 13, 2746 (2009) M.A Green, Prog Photovoltaics 9, 123 (2001) T Gokmen, O Gunawan, and D B Mitzi, Appl Phys Lett 105, 033903 (2014) H Cui, X Liu, F Liu, X Hao, N Song, and C Yan, Appl Phys Lett 104, 041115 (2014) Xiaolei Liu, Hongtao Cui, Charlie Kong, Yidan Huang, Fangyang Liu, Ning Song, Gavin Conibeer, and Martin Green, Appl Phys Lett 106, 131110 (2015) M A Green, K Emery, Y Hishikawa, W Warta, and E D Dunlop, Prog Photovolt: Res Appl 23, (2015) http://www.nrel.gov/ncpv/images/efficiency_chart.jpg NREL research solar cell power conversion efficiency records chart (accessed July 2015) L Etgar, Gao, P Xue, Z Peng, Q Chandiran et al., J Am Chem Soc 134, 17396 (2012) 10 A Kojima, K Teshima, Y Shirai, and T Miyasaka, J Am Chem Soc 131, 6050 (2009) 11 M Liu, M B Johnston, and H J Snaith, Nature 501, 395 (2013) 12 J Burschka, N Pellet, S.-J Moon et al., Nature 499, 316 (2013) 13 G Hodes, Science 342, 317 (2013) 14 J.-H Im, C.-R Lee, J.-W Lee et al., Nanoscale 3, 4088 (2011) 15 H.-S Kim, C.-R Lee, J.-H Im et al., Sci Rep 2, 591 (2012) 16 I Chung, B Lee, J He, and R P H Chang, Nature 485, 486 (2012) 17 B Lee, C.C Stoumpos, N.J Zhou, F Hao et al., J Am Chem Soc 136, 15379 (2014) 18 F Chiarella, A Zappettini, and F Licci, Phys Rev B 77, 045129 (2008) 19 I Borriello, G Cantele, and D Ninno, Phys Rev B 77, 235214 (2008) 20 Y.H Chang, C.H Park, and K Matsuishi, J Korean Phys Soc 44, 889 (2004) 21 G Murtaza and I Ahmad, Physica B 406, 3222 (2011) 22 E Mosconi, A Amat, M.K Nazeeruddin, M Gratzel, and F.D Angelis, J Phys Chem C 117, 13902 (2013) 23 L Lang, J.H Yang et al., Phys Lett A 378, 290 (2014) 24 M G Brik and I V Kityk, J Phys Chem Solids 72, 1256 (2011) 25 A.W Laubengayer, O.B Billings, and A.E Newkirk, J Am Chem Soc 62, 546 (1940) 26 W Werker, Re Trav Chim Pays-Bas Belg 58, 257 (1939) 27 J.A.A Ketelaar, A.A Rietdijk, A.P van, and Staveren, Re Trav Chim Pays-Bas Belg 56, 907 (1937) 28 T.B Brill, R.C Gerhart, and W.A Welsh, J Magn Reson 13, 27 (1974) 29 G Engel, Z Kristallogr 90, 341 (1935) 30 G Kresse and J Furthmller, Comput Mater Sci 6, 15 (1996) 31 G Kresse and D Joubert, Phys Rev B 59, 1758 (1999) 32 J P Perdew, K Burke, and M Ernzerhof, Phys Rev Lett 77, 3865 (1996) 33 H.J Monkhorst and J.D Pack, Phys Rev B 13, 5188 (1976) 34 P Blaha, K Schwarz, G Madsen, D Kvasnicka, and J Luitz, Wien2k package, http://www.wien2k.at/ 35 F Tran and P Blaha, Phys Rev Lett 102, 226401 (2009) 36 Jun Hong Noh, Sang Hyuk Im, and Jin Hyuck Heo, NanoLett 13, 1764 (2013) 37 J Sun and H Tian Wang, Phys Rev B 71, 125132 (2005) 38 P E Blochl, O Jepsen, and O K Andersen, Phys Rev B 49, 16223 (1994) 39 AM 1.5G http://rredc.nrel.gov/solar/spectra/am1.5/ 40 L.-J Yu and A Zunger, Phys Rev Lett 108, 068701 (2012) 41 X Huang, T.R Paudel, S Dong, and E.Y Tsymbal, Phys Rev B 92, 125201 (2015) 42 M.M Lee, J Teuscher, T Miyasaka et al., Science 338, 6107 (2012) ... ADVANCES 5, 127224 (2015) Electronic structure and optical properties of Cs2AX′2X4 (A=Ge,Sn,Pb; X′,X=Cl,Br,I) Guangtao Wang,a Dongyang Wang, and Xianbiao Shi College of Physics and Information Engineering,... Republic of China (Received 25 August 2015; accepted December 2015; published online 21 December 2015) We studied the crystal structures, electronic structures and optical properties of Cs2AX2′ X4 (A=Ge,Sn,Pb;. .. some of these MAB3 compounds are unstable in the air and moist environments In this paper, we report the crystal structures, electronic structures and absorption coefficients of Cs2AX2′ X4 (A=Ge,Sn,Pb;