Diwen Liu1,2, Wenying Zha3, Rusheng Yuan3, Benyong Lou4, Rongjian Sa1. 1. Institute of Oceanography, Fujian Key Laboratory of Functional Marine Sensing Materials, Minjiang University Fuzhou Fujian 350108 China rjsa@mju.edu.cn. 2. College of Materials and Chemical Engineering, Pingxiang University Pingxiang 337055 China. 3. State Key Laboratory of Photocatalysis on Energy and Environment, College of Chemistry, Fuzhou University Fuzhou 350108 P. R. China. 4. Research and Testing Center of Pharmaceutical Formulations, Ocean College, Minjiang University Fuzhou Fujian 350108 P. R. China lby@mju.edu.cn.
Lead-based halide perovskites have attracted great attention for optoelectronic applications in the past ten years.[1-3] The power conversion efficiency (PCE) of lead-based halide perovskite solar cells increased to a high value of 25.2% in 2019.[4] However, the toxicity of lead and poor stability still remain significant challenges for commercial applications. Extensive efforts have been paid to exploring lead-free or lead-less and stable halide perovskite materials for solar cells. Sn2+ and Ge2+ metal ions have been expected to replace Pb2+ ions in perovskites.[5-9] Sn2+ ion is highly unstable and can be easily oxidized to Sn4+,[10] which results in decreased photovoltaic performance. To date, the highest PCE of Ge-based halide perovskites is only 4.94%.[11] Two dimensional (2D) layered perovskites have been also explored such as A3B2X9 (A = Cs, Rb; B = Bi, Sb; X = I, Br, Cl).[12-17] However, these materials have large band gaps over 2.0 eV, which are not suitable for perovskite solar cells.In recent years, halide double perovskites A2B+B3+X6 have been developed as stable and non-toxic alternatives to lead-based perovskite materials, where two toxic lead ions are substituted by using monovalent B+ and trivalent B3+, such as Ag+, Bi3+, and Sb3+. Although most of the reported double perovskites show excellent stability, their band gaps are too large for application in single-junction solar cells.[18-21] A low band gap (1.66 eV) of double perovskites Cs2NaBiI6 has been synthesized, which possesses excellent stability and strong light absorption performance.[22] However, its PCE is only 0.42%, which is much lower than those of lead-based hybrid perovskites.The vacancy-ordered double perovskites with formula of A2BX6 can be regarded as a derivative structure of the traditional ABX3, where two toxic lead ions are substituted by one tetravalent ion to keep charge conservation as the conventional perovskites. Cs2SnI6 is one of the most representative perovskite materials, which belongs to a cubic structure with space group Fm3̄m.[23] The previous study classified Cs2SnI6 as zero dimensional (0D) structure.[24] However, this compound still maintains three dimensional (3D) structure and presents 3D materials' properties, so it is still considered as 3D material.[25] Cs2SnI6 exhibits air and moisture stability, and strong visible light absorption, which are advantageous properties for photovoltaic applications.[26,27] The first tested solar cell using Cs2SnI6 as absorber material in air showed a small PCE of 0.96%.[26] Meanwhile, MA2SnI6 was found to have a direct band gap of 1.81 eV with a strong absorption coefficient of ∼7 × 104 cm−1.[28] The material properties of Cs2Sn1−TeI6 were investigated to explore structure–property relationships with potential applications in photovoltaics.[29] In 2017, a promising family of all-inorganic based-Ti double perovskites were predicted to possess suitable band gaps, excellent optical absorption, and high stability.[30] In particular, the band gap of Cs2TiI6−Br can be tuned from 1.38 to 1.78 eV for single-junction and tandem solar cells.[30] A comprehensive and systematic study of halide perovskites has been performed by using density functional theory (DFT) calculations.[31] The theoretical study on the material properties of Cs2TiX6 (X = I, Br) has been investigated.[32-34] The PCE of 3.28% can be obtained for Cs2TiBr6 when it is incorporated into planar-heterojunction PSCs.[35] In 2018, Zhou et al. synthesized lead-free all-inorganic Cs2PdBr6 perovskite nanocrystals with a narrow band gap (1.69 eV) for the first time.[36] A series of lead-free Te-based double perovskites A2TeX6 (A = MA, FA, or BA; X = Br− or I−; MA = CH3NH3; FA = CH(NH2)2; BA = benzylamine) were reported to be potential materials for optoelectronic devices.[37] These Te-based double perovskites exhibited a tunable band gap, a low trap density, and a high mobility. Recently, the structural stability and electronic and optical properties of Te-based double perovskite materials were investigated by using first-principles calculations.[38] Therefore, the design of new non-/low-toxic stable halide perovskites for solar cells is inevitable.In this work, the structural stability, mechanical properties, electronic and optical properties of double perovskites Cs2Te1−TiI6 were investigated for the first time. The results show that Cs2Te0.50Ti0.50I6 has excellent stability, suitable direct band gap, and strong optical absorption, which is a promising candidate material for optoelectronic device. The mechanical properties of these materials are discussed in detail, and the results indicate that they are ductile materials except for Cs2Te0.50Ti0.50I6.
Computational details
First-principles calculations were performed by using Vienna ab initio simulation package (VASP).[39] The interaction between the core and valence electrons was described by the projector-augmented wave (PAW) method.[40] The generalized gradient approximation (GGA) of Perdew–Burke–Ernzerhof (PBE) was employed to describe the exchange–correlation functional.[41] A plane-wave energy cutoff of 500 eV was used. A 4 × 4 × 4 k-point grid was adopted for geometric optimization until the convergence criteria of energy and forces were set to be 10−5 eV and 0.01 eV Å−1, respectively. The structures were fully relaxed without any constraints. A denser 8 × 8 × 8 k-point grid was employed for electronic and optical properties calculations. Four k-points X (0.5, 0, 0), R (0.5, 0.5, 0.5), M (0.5, 0.5, 0), and G (0, 0, 0) were selected for the bandgap calculations. In order to gain accurate lattice parameters of the studied systems, the van der Waals (vdW) interactions were adopted within vdW-DFT in all calculations.[42] For Cs2TiI6, considering the underestimation of the band gap in the standard DFT calculation, we used the DFT+U method[43] to obtain accurate electronic structure. Various values of U have been tested for Ti 3d orbitals of Cs2TiI6. U = 2 eV was finally chosen for the electronic and optical calculations for the Cs2Te1−TiI6 solid solutions.
Results and discussion
Structure properties
The crystal structures of both Cs2TeI6 and Cs2TiI6 belong to the cubic phase with space group Fm3̄m.[29,30] The calculated lattice constant of Cs2TeI6 is 11.66 Å, which is in good agreement with experimental value (11.71 Å).[29] The corresponding lattice constant of Cs2TiI6 is predicted to be 11.38 Å, which is much lower than its experimental value (11.67 Å).[30] The ion radii of Te4+ and Ti4+ are 0.97 Å and 0.61 Å, respectively. Interestingly, the experimental lattice constants of Cs2TeI6 and Cs2TiI6 are almost the same. It should be noted that the lattice parameter of Cs2TiI6 is estimated to be 11.67 Å according to the experimental XRD results. The optimized structures of the Cs2Te1−TiI6 solid solutions are depicted in Fig. 1. x is the concentration of Ti4+. The x values of 0.00, 0.25, 0.50, 0.75, and 1.00 are adopted in this study. For the species with x = 0.25, 0.50, and 0.75, all the possible arrangements are considered in this work. The results show that different structures with the same doping percentage have the same stability, with calculated energy differences within less than 2 meV. The lattice constant of Cs2Te1−TiI6 decreases gradually when the concentration of Ti4+, as shown in Fig. 2. Moreover, the bond length Te–I (2.96 Å) of Cs2TeI6 is larger than the bond length Ti–I (2.76 Å) of Cs2TiI6.
Fig. 1
The optimized structures of the Cs2Te1−TiI6 solid solutions.
Fig. 2
The lattice constants of the Cs2Te1−TiI6 solid solutions.
The tolerance factor (t) has been widely used to predict the structural stability of perovskite material. The formula is defined as follows:[44]where RA, RB, and RX are the effective ionic radii for A, B, and X, respectively. In general, the perovskite structure is stable when the t value is in the range of 0.81–1.11.[45] Our calculated values of t are 0.86 and 0.97 for Cs2TeI6 and Cs2TiI6, respectively. These results confirm that the substitution of Ti4+ for Te4+ in Cs2TeI6 is feasible. To further assess the thermodynamic stability of the Cs2Te1−TiI6 solid solutions, we have calculated their formation energies. The formation energy can be obtained according to the following equation:where E(Cs2Te1−TiI6), E(CsI), E(TeI4), and E(TiI4) are the total energy of Cs2Te1−TiI6, CsI, TeI4, and TiI4, respectively. The calculated formation energies of the Cs2Te1−TiI6 solid solutions are demonstrated in Fig. 3. Our calculated result shows that the ΔH value of MAPbI3 is 0.01 eV per f.u., which is in good agreement with the previous reported data.[46] These results indicate that MAPbI3 is marginally stable. It can be seen that all the compounds exhibit good stability with larger negative ΔH values compared to MAPbI3. Moreover, it is observed that the structural stability of Cs2Te1−TiI6 decreases gradually with the increasing concentration of Ti4+. In addition, the phonon spectrum of the mixed-metal double perovskites are also studied. As shown in Fig. 4, it can be seen that none of the imaginary phonon mode exists for Cs2Te0.75Ti0.25I6 and Cs2Te0.25Ti0.75I6, which indicate that they are kinetically stable. Interestingly, Cs2Te0.50Ti0.50I6 has very small imaginary frequencies. According to the previous study,[47] all small imaginary frequencies will disappear when it is at room temperature. Therefore, Cs2Te0.50Ti0.50I6 is considered to be a stable structure at room temperature.
Fig. 3
Calculated formation energies of the Cs2Te1−TiI6 solid solutions.
Fig. 4
Phonon spectrum of (a) Cs2Te0.75Ti0.25I6, (b) Cs2Te0.50Ti0.50I6, and (c) Cs2Te0.25Ti0.75I6.
Mechanical properties
To predict the mechanical stability of double perovskite Cs2Te1−TiI6, the corresponding elastic constants were calculated using the finite strain theory.[48] For cubic structure, the Born stability criteria is given as follows:[49]The calculated results of three independent elastic constants are listed in Table 1. All the compounds satisfy the above Born stability criteria, indicating that they are mechanical stable. The bulk modulus (B) and shear modulus (G) can be obtained by using Hill approximations.[50] The Young's modulus (E) and Poisson ratio (ν) can be calculated from the values of bulk modulus and shear modulus by applying the following formula:[51]
The calculated elastic constants C (GPa), bulk modulus B (GPa), shear modulus G (GPa), Young's modulus (E), Pugh's ratio (B/G), and Poisson's ratio ν of Cs2Te1−TiX6
C11
C12
C44
B
G
E
B/G
ν
Cs2TeI6
19.14
10.24
8.78
13.21
6.68
17.16
1.98
0.28
Cs2Te0.75Ti0.25I6
20.20
9.80
9.86
13.27
7.63
19.20
1.74
0.26
Cs2Te0.50Ti0.50I6
17.51
8.39
9.98
11.43
7.29
18.03
1.57
0.24
Cs2Te0.25Ti0.75I6
19.14
11.08
10.27
13.77
7.06
18.08
1.95
0.28
Cs2TiI6
17.85
10.29
10.70
12.81
7.05
17.88
1.82
0.27
In order to investigate the ductility and brittleness of the studied double perovskites, we calculated the Pugh's and Poisson ratios. The Pughs's ratio is the ratio between bulk modulus and shear modulus (B/G), which is proposed by Pugh.[52] The material is considered as ductile when the value of Pugh's ratio is larger than 1.75, otherwise it is brittle. The Poisson's ratio can also be used to separate the ductile materials from brittle materials. The ductile materials have the Poisson's ratio greater than 0.26, and it is smaller than 0.26 for brittle materials.[53] From Table 1, it can be seen that Cs2Te0.50Ti0.50I6 is a brittle material, while other four compounds are ductile materials. It should be noted that the Poisson's ratio values of the studied double perovskites are close to the required lowest limit. The smaller the value of Young's modulus, the better flexibility of a material. The results show that Young's modulus of Cs2TeI6 is the smallest among the five structures, which indicates that it is the most flexible. Moreover, it is beneficial to form high quality film for ductile materials.
Electronic properties
The band structures calculated along high symmetry directions in the Brillouin zone are shown in Fig. 5. Cs2TeI6 has an indirect band gap between G and R point, and a direct band gap at R point. The indirect and direct band gaps of Cs2TeI6 are 1.12 and 1.48 eV, respectively. The previous reported band gap of Cs2TeI6 is about 1.5 eV,[54,55] and the recent experimentally observed optical gap is 1.59 eV.[29] The band gap of Cs2TiI6 with the PBE method is 0.75 eV, which is lower than its experimental value (1.02 eV).[30] Therefore, the DFT+U method is employed to obtain accurate band gap of Cs2TiI6. The result shows that Cs2TiI6 possesses a direct band gap with 0.98 eV at the G point, which is in good agreement with its experimental value.[30] The Spin–orbit coupling (SOC) effect on the electronic structures of Te- and Ti-based perovskites is small.[30,37] It is observed that Cs2Te0.50Ti0.50I6 is a direct band gap semiconductor with the value of 1.09 eV, while Cs2Te0.75Ti0.25I6 is an indirect band gap semiconductor. The direct band gap of Cs2Te0.25Ti0.75I6 is 1.05 eV. The trend in band gap is displayed in Fig. 6. The band gap of Cs2TeI6 reduces gradually when the ratio of Ti4+ increases. The band gaps of the two calculations with PBE and PBE+U show the same tendency. The variation in band gap is only 0.11 eV when the concentration of Ti4+ changes from 0.50 to 1.00. It should be noted that the direct band gap of Cs2TeI6 is close to the experimental value. The direct band gap of Cs2Te0.75Ti0.25I6 is 1.31 eV, which is a potential candidate for single-junction solar cells according to the Shockley–Queisser theory.[56] An indirect-to-direct band gap transition can be observed when the doping content of Ti4+ increases from 0.25 to 0.50. Moreover, the direct band gap is significantly reduced by 0.22 eV. It can be seen that the band gaps of the Cs2Te1−TiI6 solid solutions are apparently underestimated at the PBE method, as shown in Fig. 6.
Fig. 5
Calculated band structures of double perovskites along the X–R–M–G–R path for (a) Cs2TeI6, (b) Cs2Te0.75Ti0.25I6, (c) Cs2Te0.50Ti0.50I6, (d) Cs2Te0.25Ti0.75I6, and (e) Cs2TiI6.
Fig. 6
Experimental and theoretical band gap trends with different methods of Cs2Te1−TiI6.
Fig. 7 shows the calculated density of states of the Cs2Te1−TiI6 solid solutions. For Cs2TeI6, the valence band maximum (VBM) is mainly contributed by the I-p orbitals, while the conduction band minimum (CBM) is dominated by the Te-p orbitals along with the I-p orbitals. For Cs2TiI6, the VBM is mainly composed of the I-p orbitals, whereas the CBM is composed of the Ti-d orbitals. The similar situation is observed for the mixed Te–Ti double perovskites with respect to Cs2TiI6. It can be seen that Cs+ does not contribute to the band edge states.
Fig. 7
Density of states of (a) Cs2TeI6, (b) Cs2Te0.75Ti0.25I6, (c) Cs2Te0.50Ti0.50I6, (d) Cs2Te0.25Ti0.75I6, and (e) Cs2TiI6.
Optical properties
The calculated optical absorption spectrum is given in Fig. 8. Cs2TeI6 shows a large absorption coefficient in the visible light region. The absorption coefficient of Cs2TiI6 is over 10−5 cm−1 in the whole region, but it is lower than that of Cs2TeI6 in the range of 300–700 nm. For Cs2TiI6, two strong absorption peaks are located at about 450 and 800 nm, respectively. The optical absorption ability gradually increases when the proportion of Ti4+ decreases. In general, the Cs2Te1−TiI6 solid solutions have strong and wide optical absorption in the visible spectrum. Cs2Te0.50Ti0.50I6 is the most promising candidate due to its superior stability, suitable direct band gap, and strong optical absorption.
Fig. 8
The optical spectra of the Cs2Te1−TiI6 solid solutions.
Conclusions
In summary, the structural, mechanical, electronic and optical properties for the Cs2Te1−TiI6 solid solutions have been investigated by using first-principles calculations for the first time. The calculated results show that Cs2TeI6 is an indirect band gap semiconductor, while Cs2TiI6 is a direct band gap semiconductor with the value of 0.98 eV. These compounds show good structural stability compared to CH3NH3PbI3. An indirect–direct band gap transition can be observed when the doping concentration of Ti4+ is 0.50. Cs2Te0.50Ti0.50I6 has a suitable direct band gap (1.09 eV), which is considered as a promising photovoltaic material for single-junction solar cells. Further analysis reveals that all the structures are ductile materials except for Cs2Te0.50Ti0.50I6. Moreover, these materials have large optical absorption coefficients in the visible light region. We expect that this study can provide insights into developing the stable, non-toxic, and high-efficiency perovskite materials for optoelectronic devices.
Authors: George Volonakis; Amir Abbas Haghighirad; Rebecca L Milot; Weng H Sio; Marina R Filip; Bernard Wenger; Michael B Johnston; Laura M Herz; Henry J Snaith; Feliciano Giustino Journal: J Phys Chem Lett Date: 2017-01-31 Impact factor: 6.475
Authors: Jaya Pal; Suman Manna; Anirban Mondal; Shyamashis Das; K V Adarsh; Angshuman Nag Journal: Angew Chem Int Ed Engl Date: 2017-10-09 Impact factor: 15.336
Authors: Annalise E Maughan; Alex M Ganose; Mitchell M Bordelon; Elisa M Miller; David O Scanlon; James R Neilson Journal: J Am Chem Soc Date: 2016-07-05 Impact factor: 15.419