2 Calculations Combined with a ... - ACS Publications

Jul 13, 2019 - maximum absorbance with a limit for the red shift at x = 1.0, the cubic alloy sets a .... ionized (ae-1/2) atoms, trimmed by a step fun...
1 downloads 0 Views 2MB Size
Letter Cite This: J. Phys. Chem. Lett. 2019, 10, 4245−4251

pubs.acs.org/JPCL

Relativistic DFT-1/2 Calculations Combined with a Statistical Approach for Electronic and Optical Properties of Mixed Metal Hybrid Perovskites Diego Guedes-Sobrinho,*,†,‡ Ivan Guilhon,*,‡ Marcelo Marques,*,‡ and Lara K. Teles*,‡ †

São Carlos Institute of Physics, University of São Paulo, P.O. Box 369, 13560-970 São Carlos, SP, Brazil Grupo de Materiais Semicondutores e Nanotecnologia, Instituto Tecnológico de Aeronáutica, DCTA, 12228-900 São José dos Campos, São Paulo, Brazil

Downloaded via UNIV OF SOUTHERN INDIANA on July 18, 2019 at 09:52:09 (UTC). See https://pubs.acs.org/sharingguidelines for options on how to legitimately share published articles.



S Supporting Information *

ABSTRACT: We overcome the great theoretical computational challenge of mixed perovskites, providing a rigorous and efficient model by including quasiparticle, spin−orbit coupling, and disorder effects. As a benchmark, we consider the mixed MAPb1−xSnxI3 perovskites. The calculations are based on the generalized quasichemical approach and the DFT-1/2 approximated quasiparticle correction. Both cubic and tetragonal structures are investigated. By mapping the entire range of compositions, we correctly describe the bowing-like behavior for the energy gaps with 1.24 eV as the minimum value at x = 0.70, in very good agreement with the experimental data. Furthermore, while the tetragonal alloy reaches the maximum absorbance with a limit for the red shift at x = 1.0, the cubic alloy sets a maximum absorbance/red shift for the optimal composition at x = 0.70.

M

optical properties.14−16 As such, the theoretical prediction of electronic properties of MHP alloys requires that we overcome an additional challenge, that is the bowing-like curve with respect to the composition for the energy gap of MHP alloys. The challenge of rigorous first-principles calculation of a MHP mixture is threefold: (i) the inclusion quasiparticle (QP) effects, (ii) strong spin−orbit coupling (SOC), and (iii) disorder effects due to alloying. Therefore, standard density functional theory (DFT) calculations are incapable of providing reliable electronic and optical properties across the range of compositions.17−20 An alternative would be the relativistic GW approximation, which is accurate for electronic structure description,21,22 but its application in a statistical approach is infeasible due to the extremely high computational cost. In this sense, an efficient and accurate theoretical approach for obtaining the electronic and optical properties is highly desirable and could open an avenue for a deeper understanding of mixed MHP, such as MAPb1−xSnxI3. In this Letter, we present a protocol for calculating accurate band gap energies and absorption coefficients for perovskite alloys throughout the entire composition range, which is also applied by regarding different symmetries for the structure, e.g., cubic and tetragonal. Our framework combines a computationally efficient first-principles method based on the DFT-1/2 method for approximated quasiparticle corrections,

etal halide perovskites (MHP) as photovoltaic materials, e.g., the methylammonium (MA) lead iodide, CH3NH3PbI3 (MAPbI3), with efficiencies in excess of 22%, have attracted a significant amount of interest in the photovoltaic solar cell community.1−3 Replacing or complementing Si-based devices with MHP offers intriguing possibilities with respect to higher energy conversion efficiencies and lower processing costs.4,5 However, significant additional effort has been spent on finding alternatives for toxic Pb by preserving the material stability against moisture and air, among others.6−8 An alternative is to replace Pb with other nontoxic or less toxic cations. Among the possible candidates, Sn is one of the most viable, because it provides fewer lattice distortions. Moreover, the mixed MAPb1−xSnxI3 perovskites have been extensively investigated to extend the lightharvesting region. Experiments have revealed the stabilization of MHP mainly as tetragonal (I4/mcm) and pseudocubic (P4mm) structures, so that the increase in symmetry from I4/mcm to P4mm (second-order phase transition)9,10 is correlated with the increase in temperature. For instance, MAPbI3 is stable as tetragonal (lattice parameters given by a = 8.84 Å, b = 12.58 Å, and c = 8.55 Å) between 163 and 328 K,11,12 moving toward pseudocubic (a = 6.31 Å) for temperatures above 328 K.13 Nevertheless, for MHP alloys, e.g., the Pb/Sn mixture, the structural motif is also dependent on the composition. For MAPb1−xSnxI3 over a few compositions, experiments have revealed the I4/mcm to P4mm phase transition for x > 0.5, promoting a deviation from linearity of its electronic and © XXXX American Chemical Society

Received: May 24, 2019 Accepted: July 13, 2019 Published: July 13, 2019 4245

DOI: 10.1021/acs.jpclett.9b01499 J. Phys. Chem. Lett. 2019, 10, 4245−4251

Letter

The Journal of Physical Chemistry Letters including SOC, with a statistical treatment named generalized quasichemical approximation (GQCA), which describes the statistical ensemble of a solid solution of MHP alloys with a cluster expansion approach.23 This approach yielded electronic and optical parameters in excellent agreement with those from experiments with MAPb 1−x Snx I3 alloys for the whole composition range and an arbitrary growth temperature. All of the total energy calculations and structure relaxation were DFT-based24,25 within the Perdew−Burke−Ernzerhof (PBE)26 approximation for the electronic exchange-correlation functional (Exc). We also considered SOC for the valence electrons, which is an important relativistic treatment for nonspherical atomic orbitals in Pb-based perovskites. Thus, SOC affects the directionality of the bonds,27 and yields the Rashba−Dresselhaus splitting in band edges.17,28,29 We used the VASP30,31 code to solve the Kohn−Sham (KS) equations, from which a plane-wave cutoff energy of 500 eV was used, and to integrate over the Brillouin zone was considered a Monkhorst−Pack k-mesh of 4 × 4 × 4 (3 × 6 × 4) for a cubic structure. Shape, lattice, and atomic positions were fully optimized via minimization of Hellmann−Feynman forces on every atom, for values of 0.5. In view of the electronic properties of the MAPb1−xSnxI3 alloys determined via the DFT-1/2+SOC approach, we turn to the average optical absorption within our statistical treatment at some compositions aiming to observe its correspondence with the lowest Eg(x). Figure 5a shows the absorbance spectra

Figure 4. Energy gaps (solid lines) and their standard deviation (shaded regions) calculated within the GQCA statistical approach for the MAPb1−xSnxI3 alloy as cubic (tetragonal), employing the standard (a) DFT and (b) SOC-DFT-1/2 approaches. The empty diamonds are the experimental values (diamond from ref 14 and both square and circle from refs 15 and 16) varying as bowing-like for x = 0.0, 0.25, 0.50, 0.75, and 1.0 compositions.

DFT also fails in describing the empirically observed bowing, from which the minimum of the energy gap lies in the window Eg(min) = 1.17−1.23 eV for x = 0.75.14−16 Conversely, the relativistic interaction through the DFT+SOC enlarges the difference between the average Eg values for the cubic and tetragonal structures (Figure S4a), reproducing the bowing shape for the Eg curves. However, DFT+SOC strongly reduces the Eg values over a range of compositions, such that the magnitude of the energy gap is very far from the experimental results. Via application of the quasiparticle correction method, here as DFT-1/2+SOC in Figure 4b, average Eg values for the MAPb1−xSnxI3 alloys, connecting the MAPbI3 and MASnI3 pure perovskites (from which the accuracy has already been discussed in Figure 1), are in good agreement with the experimental values. DFT-1/2+SOC describes the bowing observed by the experiments, wherein we find a minimum Eg(min) of 1.24 ± 0.21 eV, i.e., average and standard deviation (shaded regions in Figure 4), through the GQCA approach for the cubic structure. This result lies at x = 0.70, differing between 0.81 and 5.98% in excellent agreement with the experimental result (1.17−1.23 eV) at x = 0.75.14−16 On the basis of that, we estimated a bowing parameter (b) for the energy gap dependency with composition as Eg(x) = xEg,MASnI3 + (1 − x)Eg,MAPbI3 − bx(1 − x). We found for the cubic (tetragonal) form a bowing value of b = 0.69 eV (0.54 eV), from which the error does not exceed 0.007 eV (0.006 eV). Therefore, the relativistic and quasiparticle approach to our statistical method shows an excellent agreement with empirical

Figure 5. Average optical absorbance for MAPb1−xSnxI3 alloys in (a) cubic and (b) tetragonal forms for the x = 0.0−1.0 composition range within the GQCA statistical approach at 300 K. The small charts expand the photon energy region between 1.2 and 1.7 eV, which permit the observation of the maximum red shift for the cubic form at x = 0.7 (due to the minimum energy gap from the bowing-like curve) and the tetragonal form at x = 1.0 (due to the gradual increase in the optical absorption from x = 0.0 to x = 1.0).

for the cubic structure, where notably one observes a maximum red shift for 0.5 < x < 0.8 compositions in close correlation with their correspondent lowest Eg(x) values. Our optical absorption results computed via DFT-1/2+SOC are in agreement with empirical investigations14−16 through transformed Kubelka−Munk treatment from near-infrared absorption spectra. Because the lowest Eg(x) values from the bowinglike curve are closer to the x = 1.0 value for the tetragonal structure (panel b), the maximum red shift lies on the solid solution compositions with a large excess of Sn, nearly as the MASnI3 pure perovskite. Finally, to provide an overview of all of the DFT-1/2+SOC results for the statistical ensemble of the MAPb1−xSnxI3 alloys, we correlate the energy gap, lattice parameter, and total absorbance average in Figure 6 as a complete and accurate representation of the systems. The cubic structure, panel a, depicts the coincident maximum red shift with the minimum energy gap from the bowing-like curve (Eg = 1.24 eV), from which the lattice constant around 6.31 Å correlates with the x 4248

DOI: 10.1021/acs.jpclett.9b01499 J. Phys. Chem. Lett. 2019, 10, 4245−4251

Letter

The Journal of Physical Chemistry Letters

skite alloy at x = 0.70, i.e., composition with the smallest amount of Pb presenting the best photovoltaic and optical performance, suggests the importance of seeking eco-friendly materials from the reduction of toxic chemical species. Therefore, our protocol for mixed perovskites can assist theorists and experimentalists dedicated to the design of mixed perovskites immersed in solar cells aiming to achieve high photovoltaic performance. The approach presented here can be directly applied to similar complex systems, such as 3D and 2D mixed metal halide perovskites.



ASSOCIATED CONTENT

S Supporting Information *

The Supporting Information is available free of charge on the ACS Publications website at DOI: 10.1021/acs.jpclett.9b01499.



Additional results and methodological details about GQCA (PDF)

AUTHOR INFORMATION

Corresponding Authors

*E-mail: *E-mail: *E-mail: *E-mail:

Figure 6. Correlation between the energy gap and the total average for the lattice constant calculated within GQCA at 300 K for the MAPb1−xSnxI3 alloy in the (a) cubic and (b) tetragonal structures. The color map corresponds to the total absorbance average obtained from GQCA for the photon energy interval of 0.0−3.5 eV.

[email protected]. [email protected]. [email protected]. [email protected].

ORCID

Diego Guedes-Sobrinho: 0000-0002-3313-2822 Ivan Guilhon: 0000-0001-8882-4490 Lara K. Teles: 0000-0002-8713-2094

= 0.70 composition (as shown in Figure 3). For the tetragonal form, even though its energy gap has a minimum at 8.60 Å for the lattice constant (correlated with the x = 0.90 composition), one observes that its spectral absorption increases gradually from pure MAPbI3 to MASnI3 (point of the maximum red shift). In summary, we performed a systematic and rigorous investigation via first-principles calculations of the MAPb1−xSnxI3 perovskite alloys within the DFT-1/2 approximated quasiparticle correction, considering the spin−orbit coupling and addressing the disorder effects with a cluster expansion approach. The spin−orbit coupling is mandatory for the correct description of the bowing parameter of the energy gap with respect to the composition and for distinguishing the electronic properties of cubic and tetragonal structures. In MHP, the spin−orbit coupling is very strong, resulting in a significantly smaller energy band gap. Therefore, the inclusion of QP correction is indispensable for a proper description of the electronic band structure. The DFT-1/2+SOC approach provides precise results for energy band gaps of the pure compound and mixed compositions, as well as its bowing-like behavior, wherein the minimum of the cubic configurations at x = 0.70 under 300 K is in excellent agreement with experimental findings. Our computed minimum energy gap of 1.24 eV differs by only 0.81% from the experimental value. This allows the calculation of reliable composition-dependent absorption spectra. We identify the correlation among the lattice parameter, energy gaps, and the spectral absorption over the entire composition range. Our results indicate for the cubic MAPb1−xSnxI mixed perovskite at x = 0.70 (optimal composition for a minimum energy gap) a maximum optical absorption/red shift, which is a crucial parameter for photovoltaic performance. As such, the MAPb1−xSnxI3 perov-

Notes

The authors declare no competing financial interest.



ACKNOWLEDGMENTS The authors thank the Brazilian funding agency Coordination for Improvement of Higher Level Education (CAPES) (PVE Grants 88881.068355/2014-01 and 88887.145962/2017-00), the National Council for Scientific and Technological Development (CNPq) (Grants 308742/2016-8 and 306322/ 2017-0), and the São Paulo Research Foundation (FAPESP) (Grants 2006/05858-0 and 2012/50738-3) for financial support. The authors also thank the Scientific Computation National Laboratory for providing Santos Dumont Supercomputer resources for the performance of some of the calculations.



REFERENCES

(1) Grätzel, M. The Light and Shade of Perovskite Solar Cells. Nat. Mater. 2014, 13, 838−842. (2) Saparov, B.; Mitzi, D. B. Organic−Inorganic Perovskites: Structural Versatility for Functional Materials Design. Chem. Rev. 2016, 116, 4558−4596. (3) Ali, R.; Hou, G.-J.; Zhu, Z.-G.; Yan, Q.-B.; Zheng, Q.-R.; Su, G. Predicted Lead-Free Perovskites for Solar Cells. Chem. Mater. 2018, 30, 718−728. (4) Lee, J.-W.; Kim, S.-G.; Bae, S.-H.; Lee, D.-K.; Lin, O.; Yang, Y.; Park, N.-G. The Interplay Between Trap Density and Hysteresis in Planar Heterojunction Perovskite Solar Cells. Nano Lett. 2017, 17, 4270−4276. (5) Lee, W.; Li, H.; Wong, A. B.; Zhang, D.; Lai, M.; Yu, Y.; Kong, Q.; Lin, E.; Urban, J. J.; Grossman, J. C.; et al. Ultralow Thermal Conductivity in All-Inorganic Halide Perovskites. Proc. Natl. Acad. Sci. U. S. A. 2017, 114, 8693−8697. 4249

DOI: 10.1021/acs.jpclett.9b01499 J. Phys. Chem. Lett. 2019, 10, 4245−4251

Letter

The Journal of Physical Chemistry Letters (6) Berry, J.; Buonassisi, T.; Egger, D.; Hodes, G.; Kronik, L.; Loo, Y.; Lubomirsky, I.; Marder, S. R.; Mastai, Y.; Miller, J. S.; et al. Hybrid Organic−Inorganic Perovskites (HOIPs): Opportunities and Challenges. Adv. Mater. 2015, 27, 5102−5112. (7) Lee, S.-W.; Kim, S.; Bae, S.; Cho, K.; Chung, T.; Mundt, L. E.; Lee, S.; Park, S.; Park, H.; Schubert, M. C.; et al. UV Degradation and Recovery of Perovskite Solar Cells. Sci. Rep. 2016, 6, 38150. (8) Zhang, Y.-Y.; Chen, S.; Xu, P.; Xiang, H.; Gong, X.-G.; Walsh, A.; Wei, S.-H. Intrinsic Instability of the Hybrid Halide Perovskite Semiconductor CH3NH3 PbI3. Chin. Phys. Lett. 2018, 35, No. 036104. (9) Brivio, F.; Frost, J. M.; Skelton, J. M.; Jackson, A. J.; Weber, O. J.; Weller, M. T.; Goñi, A. R.; Leguy, A. M. A.; Barnes, P. R. F.; Walsh, A. Lattice Dynamics and Vibrational Spectra of the Orthorhombic, Tetragonal, and Cubic Phases of Methylammonium Lead Iodide. Phys. Rev. B: Condens. Matter Mater. Phys. 2015, 92, 144308. (10) Jaffe, A.; Lin, Y.; Karunadasa, H. I. Halide Perovskites Under Pressure: Accessing New Properties Through Lattice Compression. ACS Energy Lett. 2017, 2, 1549−1555. (11) 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 (CH3NH3)PbI3 for Solid-State Sensitised Solar Cell Applications. J. Mater. Chem. A 2013, 1, 5628− 5641. (12) G, S.; Mahale, P.; Kore, B. P.; Mukherjee, S.; Pavan, M. S.; De, C.; Ghara, S.; Sundaresan, A.; Pandey, A.; Guru Row, T. N.; Sarma, D. D. Is CH3NH3 PbI3 Polar? J. Phys. Chem. Lett. 2016, 7, 2412− 2419. (13) 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. (14) Hao, F.; Stoumpos, C. C.; Chang, R. P. H.; Kanatzidis, M. G. Anomalous Band Gap Behavior in Mixed Sn and Pb Perovskites Enables Broadening of Absorption Spectrum in Solar Cells. J. Am. Chem. Soc. 2014, 136, 8094−8099. (15) Ogomi, Y.; Morita, A.; Tsukamoto, S.; Saitho, T.; Fujikawa, N.; Shen, Q.; Toyoda, T.; Yoshino, K.; Pandey, S. S.; Ma, T.; et al. CH3NH3 Snx Pb(1‑x)I3 Perovskite Solar Cells Covering up to 1060 nm. J. Phys. Chem. Lett. 2014, 5, 1004−1011. (16) Hu, H.; Dong, B.; Zhang, W. Low-Toxic Metal Halide Perovskites: Opportunities and Future Challenges. J. Mater. Chem. A 2017, 5, 11436−11449. (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) Brivio, F.; Butler, K. T.; Walsh, A.; van Schilfgaarde, M. Relativistic Quasiparticle Self-Consistent Electronic Structure of Hybrid Halide Perovskite Photovoltaic Absorbers. Phys. Rev. B: Condens. Matter Mater. Phys. 2014, 89, 155204. (19) Geng, W.; Zhang, L.; Zhang, Y.; Lau, W.; Liu, L. FirstPrinciples Study of Lead Iodide Perovskite Tetragonal and Orthorhombic Phases for Photovoltaics. J. Phys. Chem. C 2014, 118, 19565−19571. (20) Jiang, L.; Wu, T.; Sun, L.; Li, Y.-J.; Li, A.; Lu, R.; Zou, K.; Deng, W.-Q. First-Principles Screening of Lead-Free Methylammonium Metal Iodine Perovskites for Photovoltaic Application. J. Phys. Chem. C 2017, 121, 24359−24364. (21) Umari, P.; Mosconi, E.; De Angelis, F. Relativistic GW calculations on CH3NH3 PbI3 and CH3NH3 SnI3 Perovskites for Solar Cell Applications. Sci. Rep. 2015, 4, 4467. (22) Mosconi, E.; Umari, P.; De Angelis, F. Electronic and Optical Properties of Mixed Sn−Pb Organohalide Perovskites: a First Principles Investigation. J. Mater. Chem. A 2015, 3, 9208−9215. (23) Brivio, F.; Caetano, C.; Walsh, A. Thermodynamic Origin of Photoinstability in the CH3NH3 Pb(I1‑xBrx) 3 Hybrid Halide Perovskite Alloy. J. Phys. Chem. Lett. 2016, 7, 1083−1087. (24) Hohenberg, P.; Kohn, W. Inhomogeneous Electron Gas. Phys. Rev. 1964, 136, B864−B871.

(25) Kohn, W.; Sham, L. J. Self-Consistent Equations Including Exchange and Correlation Effects. Phys. Rev. 1965, 140, A1133− A1138. (26) Perdew, J. P.; Burke, K.; Ernzerhof, M. Generalized Gradient Approximation Made Simple. Phys. Rev. Lett. 1996, 77, 3865−3868. (27) Manser, J. S.; Christians, J. A.; Kamat, P. V. Intriguing Optoelectronic Properties of Metal Halide Perovskites. Chem. Rev. 2016, 116, 12956−13008. (28) Even, J.; Pedesseau, L.; Dupertuis, M.; Jancu, J.; Katan, C. Electronic Model for Self-Assembled Hybrid Organic/Perovskite Semiconductors: Reverse Band Edge Electronic States Ordering and Spin-Orbit Coupling. Phys. Rev. B: Condens. Matter Mater. Phys. 2012, 86, 205301. (29) Even, J.; Pedesseau, L.; Katan, C.; Kepenekian, M.; Lauret, J.; Sapori, D.; Deleporte, E. Solid-State Physics Perspective on Hybrid Perovskite Semiconductors. J. Phys. Chem. C 2015, 119, 10161− 10177. (30) Kresse, G.; Hafner, J. Ab-Initio Molecular Dynamics for OpenShell Transition Metals. Phys. Rev. B: Condens. Matter Mater. Phys. 1993, 48, 13115−13118. (31) Kresse, G.; Furthmüller, J. Efficient Iterative Schemes for AbInitio Total-Energy Calculations Using a Plane-Wave Basis Set. Phys. Rev. B: Condens. Matter Mater. Phys. 1996, 54, 11169−11186. (32) Adolph, B.; Gavrilenko, V. I.; Tenelsen, K.; Bechstedt, F.; Del Sole, R. Nonlocality and Many-Body Effects in the Optical Properties of Semiconductors. Phys. Rev. B: Condens. Matter Mater. Phys. 1996, 53, 9797−9808. (33) Gajdoš, M.; Hummer, K.; Kresse, G.; Furthmü ller, J.; Bechstedt, F. Linear Optical Properties in the ProjectorAugmented Wave Methodology. Phys. Rev. B: Condens. Matter Mater. Phys. 2006, 73, No. 045112. (34) Ferreira, L. G.; Marques, M.; Teles, L. K. Approximation to Density Functional Theory for the Calculation of Band Gaps of Semiconductors. Phys. Rev. B: Condens. Matter Mater. Phys. 2008, 78, 125116. (35) Matusalem, F.; Marques, M.; Teles, L. K.; Bechstedt, F. Stability and Electronic Structure of Two-Dimensional Allotropes of Group-IV Materials. Phys. Rev. B: Condens. Matter Mater. Phys. 2015, 92, No. 045436. (36) Ataide, C. A.; Pelá, R. R.; Marques, M.; Teles, L. K.; Furthmü ller, J.; Bechstedt, F. Fast and Accurate Approximate Quasiparticle Band Structure Calculations of ZnO, CdO, and MgO Polymorphs. Phys. Rev. B: Condens. Matter Mater. Phys. 2017, 95, No. 045126. (37) Lucatto, B.; Assali, L. V. C.; Pela, R. R.; Marques, M.; Teles, L. K. General Procedure for the Calculation of Accurate Defect Excitation Energies from DFT1/2 Band Structures: The Case of the NV− Center in Diamond. Phys. Rev. B: Condens. Matter Mater. Phys. 2017, 96, No. 075145. (38) Guilhon, I.; Koda, D. S.; Ferreira, L. G.; Marques, M.; Teles, L. K. Approximate Quasiparticle Correction for Calculations of the Energy Gap in Two-Dimensional Materials. Phys. Rev. B: Condens. Matter Mater. Phys. 2018, 97, No. 045426. (39) Tao, S. X.; Cao, X.; Bobbert, P. A. Accurate and Efficient Band Gap Predictions of Metal Halide Perovskites Using the DFT-1/2 Method: GW Accuracy with DFT Expense. Sci. Rep. 2017, 7, 14386. (40) Rodrigues Pela, R.; Gulans, A.; Draxl, C. The LDA1/2 Method Applied to Atoms and Molecules. J. Chem. Theory Comput. 2018, 14, 4678−4686. (41) Stoumpos, C. C.; Cao, D. H.; Clark, D. J.; Young, J.; Rondinelli, J. M.; Jang, J. I.; Hupp, J. T.; Kanatzidis, M. G. Ruddlesden−Popper Hybrid Lead Iodide Perovskite 2D Homologous Semiconductors. Chem. Mater. 2016, 28, 2852−2867. (42) Ferreira, L. G.; Marques, M.; Teles, L. K. Slater HalfOccupation Technique Revisited: the LDA-1/2 and GGA-1/2 Approaches for Atomic Ionization Energies and Band Gaps in Semiconductors. AIP Adv. 2011, 1, No. 032119. 4250

DOI: 10.1021/acs.jpclett.9b01499 J. Phys. Chem. Lett. 2019, 10, 4245−4251

Letter

The Journal of Physical Chemistry Letters (43) Slater, J. C.; Johnson, K. H. Self-Consistent-Field Xα Cluster Method for Polyatomic Molecules and Solids. Phys. Rev. B 1972, 5, 844−853. (44) Hao, F.; Stoumpos, C. C.; Cao, D. H.; Chang, R. P. H.; Kanatzidis, M. G. Lead-Free Solid-State Organic-Inorganic Halide Perovskite Solar Cells. Nat. Photonics 2014, 8, 489−494. (45) Yamada, Y.; Nakamura, T.; Endo, M.; Wakamiya, A.; Kanemitsu, Y. Near-Band-Edge Optical Responses of SolutionProcessed Organic−Inorganic Hybrid Perovskite CH3NH3 PbI3 on Mesoporous TiO 2 electrodes. Appl. Phys. Express 2014, 7, No. 032302. (46) Leguy, A. M. A.; Azarhoosh, P.; Alonso, M. I.; Campoy-Quiles, M.; Weber, O. J.; Yao, J.; Bryant, D.; Weller, M. T.; Nelson, J.; Walsh, A.; van Schilfgaarde, M.; et al. Experimental and Theoretical Optical Properties of Methylammonium Lead Halide Perovskites. Nanoscale 2016, 8, 6317−6327. (47) Phuong, L. Q.; Yamada, Y.; Nagai, M.; Maruyama, N.; Wakamiya, A.; Kanemitsu, Y. Free Carriers versus Excitons in CH3NH3 PbI3 Perovskite Thin Films at Low Temperatures: Charge Transfer from the Orthorhombic Phase to the Tetragonal Phase. J. Phys. Chem. Lett. 2016, 7, 2316−2321. (48) Fujiwara, H.; Collins, R. W. Spectroscopic Ellipsometry for Photovoltaics: Vol. 1: Fundamental Principles and Solar Cell Characterization; Springer Series in Optical Sciences; Springer, 2018. (49) Teles, L. K.; Furthmüller, J.; Scolfaro, L. M. R.; Leite, J. R.; Bechstedt, F. First-Principles Calculations of the Thermodynamic and Structural Properties of Strained InxGa1−xN and AlxGa1−xN Alloys. Phys. Rev. B: Condens. Matter Mater. Phys. 2000, 62, 2475−2485. (50) Guilhon, I.; Teles, L. K.; Marques, M.; Pela, R. R.; Bechstedt, F. Influence of Structure and Thermodynamic Stability on Electronic Properties of Two-Dimensional SiC, SiGe, and GeC Alloys. Phys. Rev. B: Condens. Matter Mater. Phys. 2015, 92, No. 075435.

4251

DOI: 10.1021/acs.jpclett.9b01499 J. Phys. Chem. Lett. 2019, 10, 4245−4251