High-throughput screening of stable and efficient double inorganic halide perovskite materials by DFT

Abstract

Perovskite solar cells have become the most promising third-generation solar cells because of their superior physical–chemical properties and high photoelectric conversion efficiency. However, the current obstacles to commercialization of perovskite solar cells are their poor stability and harmful elements. How to find high-efficiency, high-stability and non-toxic perovskite materials from thousands of possible perovskite crystals is the key to solve this problem. In this paper, the inorganic halide double perovskite A₂BX₆ and its crystal structure are considered, and the data mining algorithm in informatics is introduced into the high-throughput computing data to analyze various elements in nature to study the perovskite materials that can meet the requirements of high performance. The photoelectric conversion properties and stability of 42 inorganic double perovskite materials are studied based on density functional theory (DFT). The results show that the tolerance factors of 39 crystals are between 0.8 and 1.10, indicating that these crystals have stable perovskite structure. In addition, the dielectric function, PDOS, elastic modulus, shear modulus and poison’s ratio of these crystals are analyzed. According to the above theoretical simulation results, three candidate materials for ideal light absorption are presented. This can provide a theoretical basis for the industrial application of perovskite solar cells.

Introduction

Halide perovskite materials have shown great potential as light absorbers in perovskite solar cells and which have been aroused great enthusiasm by researchers^1,2,3,4. Because they have many excellent properties such as tunable band gap^5,6,7, strong light absorption^8,9, long carrier diffusion length^10,11, and small exciton binding energy^12,13. Since the first application in 2009 for perovskite materials, it has developed rapidly in just a few years^{14,15,16,17,18}.The photoelectric conversion efficiency (PCE) of metal halide perovskite solar cells rose from 3.8% for dye-sensitized solar cell configurations to 25.7% for planar heterojunction cells which is certified by National Renewable Energy Laboratory (NREL)^19,20,21,22. However, the current obstacle to commercialization is the stability of perovskite solar cells and its harmful elements^23,24,25,26. The specific solution is to improve humidity stability, allowing modules to operate in outdoor environments without the need for expensive packaging, as well as replacing the harmful element lead (Pb) with other non-toxic metal halides.

At present, these problems have not been well resolved. How to construct perovskite crystal materials that meet the above requirements is the key to solving the problem. As we know that there are 118 elements in the periodic table. It is difficult to screen out high-efficiency, high-stability and non-toxic perovskite materials from thousands of possible combinations for perovskite crystals. Fortunately, the situation would improve if we expand our search from ternary A⁺B²⁺X₃ perovskites to quaternary A₂⁺B⁴⁺X₆ double perovskites. The combined advantages of the halide double perovskite Cs₂TiX₆ with low cost, high efficiency, nontoxicity, stability, and tunable bandgap have also been recently reported^{27,28,29,30,31,32}. On the other hand, many studies have shown that replacing MA⁺ or FA⁺ with cesium ion (Cs⁺) or rubidium ion (Rb⁺) can significantly improve the thermal stability of materials¹³. These perovskite materials are considered suitable as carrier transport materials (CTM) for CsPbBr₃ and Cs₂TiBr₆ in light-collecting or light-emitting devices. Most of double inorganic halide perovskite materials are direct band gap semiconductor materials, and they have strong visible light absorption³³. There are many tetravalent ions in nature, such as Ge⁴⁺, Zr⁴⁺, Sn⁴⁺, Hf⁴⁺, Se⁴⁺, Te⁴⁺ and Pd⁴⁺, which can combine with monovalent metal ions Cs⁺ and Rb⁺ to form double perovskite structures. These materials have suitable electronic configurations and are very stable in air exposure.

Herein, the properties of inorganic metal halide perovskite materials are comprehensively discussed based on DFT. And the organization of this paper mainly includes four parts. First of all, we describe the calculation method and specific parameters used in this work. Then, we consider inorganic halide double perovskite A₂BX₆ and its crystal structure with the method of high-throughput screening. Thirdly, we further simulation calculated the parameters of these materials, such as band structure (the position of conduction band and valence band), electronic density of states, absorption spectrum, etc. Finally, the thermodynamic stability of these crystals is analyzed. The results show that Cs₂SeI₆, Rb₂SeI₆, Cs₂SeBr₆ and Rb₂SeBr₆ can be used as ideal candidates for light absorption materials. This work can not only guide experiments and design experiments rationally, but also reduce research and development costs, shorten research time, and improve the success rate of material design.

Computational methods

In this work we discuss the crystal structure, electronic band structure, elastic modulus and optical properties of perovskite crystals in detail. The calculations in this paper are mainly performed by the plane wave pseudopotentials method as implemented in Cambridge serial total energy package (CASTEP) in Material Studio^33,34,35,36, and the conduction band and valence band edge positions are calculated with the DMOL³ module^37,38. The theoretical basis of both computational modules is density functional theory (DFT). Further, the generalized gradient approximation (GGA) within the Perdew–Burke–Ernzerhof (PBE) functional is used for the evaluation of exchange–correlation interactions for two computational modules. The lattice parameters of the crystal are shown in Fig. 4, and k pionts of supercell are set as 2 × 2 × 2. The calculation of the Dmol³ software package adopts the numerical basis set, which is set as 4.4 Double Numerical plus d-functions basis (DND). The accuracy and qualities are selected as "fine". SCF tolerance is 10^–5 eV/atom for perovskite crystals. All the other sets are the default. While the calculation of CASTEP uses plane waves basis set. The electronic configurations 5s¹, 5s²5p², 3s²3p⁵, 4s²4p⁵ and 5s² 5p⁵ are considered in valence for Rb, Sn, Cl, Br and I atoms, respectively. And the cut-off energy is set to 500 eV. The structure optimization process stops until the atomic Hellmann–Feynman force is less than 0.01 eV/Å. The energy convergence criterion is set to 10^–5 eV.

Model for perovskite solar cells

In this section, we consider inorganic halide double perovskite A₂BX₆ and its crystal structure (shown in Fig. 1), where A is a monovalent cation, B is a tetravalent cation, and X is a halide anion) type perovskite. For perovskite compounds, the octahedral factor is used to predict the formation of BX₆ octahedra, and the tolerance factor is used to predict calcium formation and distortion of the titanite structure empirically. Likewise, in the perovskite crystal-derived structure A₂BX₆, we can combine the octahedral factor and the radius ratio to predict the formation and deformation of the structure, where A = Cs, Rb, B = Ge, Zr, Sn, Hf, Se, Te and Pd, and the halide anions X = I, Cl, Br. High-throughput screening methods and properties of important elements is shown in Fig. 2. The small octahedral factor indicates that it is not favorable for the formation of BX₆ octahedra. Small radius ratios lead to evacuated cavities and lower structural symmetry, even different connectivity of fully octahedral networks. Each crystal structure parameters of these compounds are obtained from the International Crystal Structures of Database (ICSD) synthesized and characterized by experiments^39,40.

Figure 1

The crystal structure model of A₂BX_6, which is double perovskites have high-symmetric cubic.

Figure 2

(a) High-throughput screening methods for network chain structure diagram and (b) properties of important elements germanium (Ge), selenium (Se), zirconium (Zr), palladium (Pd), tin (Sn), hafnium (Hf), tellurium (Te).

To evaluate the structural stability of materials in perovskite structures, we further explore the empirical tolerance factor of them, which is an important parameter for evaluating whether cations in perovskite materials can be substituted. This factor has been widely used by previous researchers⁴¹ in the current double perovskite A₂BX₆, and the effective^42,43 tolerance factor t could be defined as

$$ t = \frac{{R_{A} + R_{X} }}{{\sqrt 2 \left[ {R_{B} + R_{X} } \right]}}. $$

(1)

In the formula, if there are multiple ions in A and B position, the average radius can be taken. It is a semi-empirical formula that can roughly describe the stability of the perovskite structure. Studies have shown that the tolerance factor of structurally stable perovskite compounds are generally between 0.78 and 1.10. In general: (1) When tolerance factor t is close to 1.0, the compound has an equiaxed Pm3m structure. (2) If tolerance factor t deviates greatly from 1.0, other structures would be formed, in which R_Cs = 167 pm, R_Rb = 152 pm, R_I = 220 pm, R_Cl = 181 pm, R_Br = 196 pm. The calculation result of tolerance factor t is shown in Fig. 3.

Figure 3

Tolerance factor of (a) Cs₂BX₆ and (b) Rb₂BX₆.

The results show that only three tolerance factors of Rb₂GeI₆, Rb₂SeI₆, Cs₂TeI₆ are (1.14, 1.15, 1.14) slightly larger than 1.10, and the tolerance factors of other crystals are all between 0.8–1.10. The results show that these combinations are in line with the properties of perovskite compound materials. Next, the lattice parameters of these crystals are analyzed, and all possible combinations of crystal compounds are constructed based on the structure of the crystals by using MS software, and the crystal parameters are optimized which is shown in Fig. 4, with the unit (Å) for crystal bond length a, b, or c.

Figure 4

Lattice parameters of A₂BX₆.

The value of the band gap of the crystal material is important factor affecting the light absorption efficiency of the perovskite layer. As we know that the smooth transition of electrons from the valence band to the conduction band can be achieved only when the photon energy $h\nu$ is greater than the band gap of the crystal. In order to analyze the optical properties of these materials, we calculated the band gaps of 42 kinds of materials with CASTEP and DMOL³ modules, as shown in Fig. 5 refer to the calculation results of some general perovskite materials in other papers ref. Table 1, our calculated band gap results are well agree with them.

Figure 5

The value of bandgaps for A₂BX₆.

Table 1 Some known band gap reference values of perovskite crystals.

It can be seen from Fig. 5 that the value of band gap for 14 kinds of crystals are between 1.33 eV and 2.40 eV, such as Cs₂ZrI₆, Cs₂HfI₆, Cs₂SeI₆, Cs₂TeI₆, Rb₂ZrI₆, Rb₂HfI₆, Rb₂SeI6, Rb₂TeI₆, Cs₂GeCl₆, Rb₂GeCl₆, Cs₂SnBr₆,Cs₂SeBr₆, Rb₂HfBr₆, Rb₂SeBr₆, they are more suitable as light-absorbing layer materials. The band gap size of the iodine-containing halide perovskite material shows relatively excellent visible light absorption performance, and its performance is better than that of the halide perovskite material containing Br and Cl.

Electronic properties of A₂BX₆

As a photoelectric conversion material for solar cells, the electronic structure of A₂BX₆ is a key factor affecting spectral absorption. Therefore, we further analyze the energy band structures of these materials, and we compared the electric potential of the back electrode of the hole transport layer in the conducting glass in perovskite solar cells. Figure 6 shows that the conduction band position of the perovskite is higher than the conduction band potential value of the electron transport layer, and the valence band position is lower than the potential value of the hole transport layer, which is more conducive to the separation of electrons and holes and forms a stable potential difference external power supply.

Figure 6

The device configuration and working principle of perovskite solar cells.

The perovskite light-absorbing layer absorbs solar photons, and electrons in the valence band obtained energy transitions to the conduction band. And holes are left in the valence band, so the electrons and holes are separated. The maximum value of the valence band of the perovskite material is lower than that of the holes transport material, while the minimum value of its conduction band is higher than that of the electron transport material. In this way, smooth separation of carriers can be easily achieved.

We calculate the valence band maxima (VBM) and conduction band minima (CBM) of 42 compounds and compared them with the band gaps of two efficient organic/inorganic perovskite materials (CH₃NH₃PbI₃) MAPbI₃, (NH₂CH = NH₂PbI₃) FAPbI₃ form. It can be seen from Fig. 7 that the maximum value of the valence band of these materials and the size of the valence band of the electron transport layer material, the hole transport layer and the back electrode material are helpful to analyze the efficiency of carrier migration of these materials. We preferentially analyze the above 14 materials with band gaps between 1.33 and 2.40 eV. Considering the energy band value of ITO, Ag and PEDOT:PSS, it requires the minimum edge of the conduction band is greater than − 4.6 eV, and the maximum value of the valence band is lower than − 5.4 eV. There are 11 kinds of crystals that satisfy the above conditions as follows: Cs₂ZrI₆, Cs₂HfI₆, Cs₂SeI₆, Cs₂TeI₆, Rb₂HfI₆, Rb₂SeI₆, Rb₂TeI₆, Rb₂GeCl₆ Cs₂SeBr₆ , Rb₂SeBr₆ and Rb₂HfBr₆ . And we found that the conduction band and valence band edge value of the four crystals such as Cs₂SeI₆, Rb₂SeI₆, Cs₂SeBr₆, Rb₂SeBr₆ are relatively close to the values of (CH₃NH₃PbI₃) MAPbI₃ and (NH₂CH=NH₂PbI₃) FAPbI_3.

Figure 7

The valence band maxima (VBM) and conduction band minima (CBM) of 42 crystals, the energy band value of ITO, Ag and PEDOT:PSS is shown in Figure with unit (eV) .

Optical properties of A₂BX₆

The dielectric function ε(ω) is used to describe the linear response properties of photovoltaic materials to electromagnetic radiation, and is one of the effective indicators to reflect the spectral properties of photovoltaic materials. Its optical response properties can be described by the following complex dielectric function ε(ω)⁴⁸

$$ \varepsilon (\omega ) = \varepsilon_{1} (\omega ) + i\varepsilon_{2} (\omega ) $$

(2)

where $\varepsilon_{1} (\omega ) = n^{2} - k^{2} ,$ $\varepsilon_{2} (\omega ) = 2nk$ (n and k are the reflection coefficient and extinction coefficient, respectively)ε₁ (ω) and ε₂ (ω) are the real and imaginary parts of the dielectric function, respectively, depending on the optical frequency. The characteristic curve of the real part of the dielectric function obtained by calculation is shown in Fig. 8, and the characteristic curve of the imaginary part is shown in Fig. 9. ω represents the frequency of light. The absorption coefficient is directly related to the band gap of the material. Optical properties can be obtained from complex dielectric functions. The absorption coefficient I (ω) can be derived from ε₁(ω) and ε₂(ω). The specific formula of the equation is shown as following⁴⁸:

$$ I(\omega ) = \sqrt 2 \omega \left[ {\sqrt {\varepsilon_{1} (\omega ) + \varepsilon_{2} (\omega )} - \varepsilon_{1} (\omega )} \right]^{{{1 \mathord{\left/ {\vphantom {1 2}} \right. \kern-\nulldelimiterspace} 2}}} $$

(3)

Figure 8

Real part characteristic curve of dielectric.

Figure 9

Imaginary part characteristic curve of dielectric.

Next, the optical absorption coefficients of these 42 compounds are calculated by the PBE method and are shown in Fig. 10. Absorption is seen at energies less than the bandgap, which is due to an error in the calculation method and did not affect our analysis of the results. The absorption coefficient is as high as 10⁵ cm^−1, which is mainly due to the absorption between the p orbitals of halide ions Cl, Br, and I in the valence band and the s-orbital of Sn in the conduction band. The optical absorption coefficients of the three structures are shown in Fig. 10. When X changed from I^- to Cl^-, the absorption spectrum of Cs₂TeX₆ is blue-shifted due to the band gap increase. Comparatively, A₂XI₆ has better absorption than A₂XCl₆ and A₂XBr₆. The light absorption curves of seven tetravalent ions (Ge, Zr, Sn, Hf, Se, Te, Pd) compounds can be seen, their light absorption peaks are concentrated in the ultraviolet region, and those with better absorption effect in the visible light region have Se and Pd-containing materials.

Figure 10

Absorption spectra of 42 crystalline materials.

In contrast, the light absorption behavior of Rb₂XCl₆ in the visible region is poor, but the light absorption in the ultraviolet region is also considerable. Others with better light absorption are halide perovskite materials containing Hf and Zr. From the absorption spectrum, it can be seen that Cs₂SeI₆, Rb₂SeI₆, Cs₂SeBr₆, Rb₂SeBr₆, Rb₂PdBr₆, Cs₂PdBr₆ are potential candidates for the absorption layer of solar cells. Comparing with the Figs. 5 and 7, we can find that the conclusion obtained according to the energy band theory is in good agreement with the absorption spectrum Fig. 10.

Electronic properties of A₂BX₆

Among halogen elements, as the atomic number increases, the ability of atoms to attract the outermost electrons gradually decreases, so the ability to obtain electrons decreases with the increase of atomic number. It can be seen from the PDOS curve in Fig. 11a–c that the curve peaks of Cl, Br, and I proceed in the direction of decreasing energy, and the energy of some electrons concentrate in the energy band decreases continuously.

Figure 11

(a–c) show the PDOS curve of Halogen element chlorine (Cl), bromine (Br), iodine (I) respectively.

Figure 12a,b illustrate the partial density of states (PDOS) for Cesium (Cs) and Rubidium (Rb). And the contribution to Rb mainly came from 5 s-orbital and 4p-orbital electrons, the element Cs mainly comes from 6 s-orbital electrons and 5d-orbital electrons. It can be seen from the Fig. 12a,b that they have no band gap, which shows that they have obvious metallicity. From the position in the periodic table of elements, Rb and Cs are at the same group (IA). The difference of them is that Rb located in the 5th periods and Cs in 6th periods. There are 55 electrons in outer cores for Cs, and 37 electrons in outer Rb cores. But their outermost layer has only one electron, which is easy to lose and become a positive monovalent ion. In contrast, the compositions of PDOS curves for the two elements are very similar. Nevertheless the electronic density of states of rubidium is more concentrated. Further analysis shows that the 6 s-orbital electrons of Cs and the 5 s-orbital electrons of Rb did not contribute to their low energy bands. Above appearance may help us to comprehend the fact of difference in the effective mass between hole and electron.

Figure. 12

PDOS curves of (a) Cs and (b) Pd.

Since selenium(Cs)-containing perovskite compounds have good light absorption properties and band gap structures, this paper uses selenium-containing compounds as an example to calculate the PDOS curves of six crystals, involving electrons 3d¹⁰ 4s² 4p⁶. It can help us to analyze the ability of Cs to lose electrons when it bonds with different halogens. The PDOS curves of Fig. 13a–f suggest that with the increase of atomic radii of Cl, Br and I, the number of exonuclear electrons increases correspondingly, and the binding ability of the nucleus and electrons decreases. Further, we found that the ability to obtain electrons decreases and their oxidizing property decreases. Longitudinal comparison we found that when Se combined with Rb to form a crystal, the area covered by the density of states would be widen.

Figure 13

PDOS curves of Se-containing perovskite materials (a) Cs₂SeI₆ (b) Cs₂SeCl₆ (c) Cs₂SeBr₆ (d) Rb₂SeI₆ (e) Rb₂SeCl₆ (f) Cs₂SeBr₆.

Stability of perovskite crystals A₂BX₆

Thermodynamical stability of A₂BX₆

To explore the thermodynamic stability of these crystals, we calculate their formation energies. For perovskite materials, the formation energy can be calculated using the following expression⁴⁹:

$$ \Delta E_{1} = E(A_{2} BX_{6} ) - 2E(AX) - E(BX_{4} ) $$

(4)

In general, the initial state of compound formation is the elemental state of A and B, and the bond between A and B atoms need to be broken first, which requires absorbing a certain amount of energy. If we want to form compound AB again, energy need to be released. The energy releases and absorbs by the difference between the two is the heat of generation, which is equivalent to the reaction heat of a chemical reaction. The more negative the formation energy, the more stable the compound. The calculation results are shown in Fig. 14a,b. The formation of compounds containing Cs and Rb elements is roughly the same, and the perovskite crystals containing Cs are more stable. Among these elements with positive tetravalent ion, there are also Zr and Hf perovskite. The crystal structure is more stable. If the band gap size and light absorption properties of Cs₂SeI₆, Rb₂SeI₆, Cs₂SeBr₆ and Rb₂SeBr₆ are considered comprehensively, they can be considered as ideal candidates for light absorption materials. In addition, Cs₂ZrI₆, Cs₂HfI₆, Rb₂ZrI₆, Rb₂HfI₆, Rb₂HfBr₆ are also superior in stability and band gap structure. Their light absorption behavior in the visible region is poor, but their light absorption in the ultraviolet region is also considerable.

Figure 14

The curve trend Figure of formation energy (a) Cs₂X (I, Cl, Br)₆, (b) Rb₂X (I, Cl, Br)₆.

Mechanical properties of A₂BX₆

Then, the mechanical stability of these materials is predicted from the calculated elastic moduli, which is calculated using finite strain theory⁵⁰. The cubic crystal has three independent elastic constants (C₁₁, C₁₂ and C₄₄).The natural mechanical stability criterion is given by: C₁₁–C₁₂ > 0, C₁₁ + 2C₁₂ > 0, C₄₄ > 0. Then we calculate elastic constants of the 42 inorganic double perovskites are shown in Tables 1, 2 and 3. According to the above criteria, these compounds are mechanically stable. To differentiate between ductile and brittle materials, Pugh's ratio (B/G) and Poisson's ratio (v) are used. Poisson's ratio can be evaluated using the following formula:

$$ \nu = (3B - 2G)/(6B + 2G) $$

(5)

where B and G are the bulk modulus and shear modulus, respectively. The critical B/G and values are 1.75 and 0.26 respectively^38,44. In other words, the composite is considered to be a ductile material when the B/G value is greater than 1.75 or the υ value is greater than 0.26, otherwise it is considered a brittle material. The calculation results are shown in Tables 2, 3 and 4. The calculated value of elastic constants C_ij, bulk modulus B (GPa), shear modulus G (GPa), Young’s modulus E (GPa), Pugh’s ratio K/G, and Poisson ratio (ν) of cubic perovskite, Averaged sound velocity V(m/s).

Table 2 Elastic Modulus of Iodide Inorganic Perovskites.

Table 3 Elastic Modulus of Chloride Inorganic Perovskites.

Table 4 Elastic moduli of bromide inorganic perovskites.

In order to analyze the stability of the crystal in more detail, we plot three elastic constants C₁₁, C₁₂ and C₄₄ sufficient to explain the complete mechanical behavior of cubic symmetry as shown in Fig. 15. It is easy to conclude that they meet the criteria of C₁₁-C₁₂ > 0, C₄₄ > 0, and C₁₁ + 2C₁₂ > 0. It indicates that their mechanical stability is relatively good.

Figure 15

The intuitive graphics of the three components of elastic modulus are (a) C_44, (b) C₁₁, C₁₂ respectively, and they can show the elastic modulus along different axes of the crystal.

The critical limit of 1.75 for the Pugh ratio (B/G) distinguishes the brittle (B/G < 1.75) and ductile (B/G > 1.75) behavior of the studied material. The values report in Tables 1, 2 and 3 show the ductility behavior, the Poisson's ratio (υ) further ensures ductility, and the critical limit for ductile materials is υ > 0.26^50,51. Figure 16 is a graph showing its Poisson's ratio and Pugh ratio, and the comparison shows that Rb₂TeI₆, Cs₂TeI₆, Rb₂HfCl₆, Cs₂PdBr₆, Cs₂SeBr₆, Rb₂TeBr₆, Rb₂SeBr₆ are more ductile, while others are more brittle.

Figure 16

(a) Intuitive plot of Poisson's ratio for 42 kinds of inorganic double perovskite materials, the figure shows that the Poisson ratio of Cs₂TeCl₆ is the largest, (b) shows the B/G of Rb₂TeI₆ is the 42 kinds of them.

Conclusion

Here we have reported the photoelectric conversion properties and stability of 42 inorganic double perovskite materials based on density functional theory. Starting from the construction of double perovskite crystal structures, the data mining algorithm in informatics is introduced into high-throughput computing data analysis. First of all we discuss the valence band maximum (VBM) and conduction band minimum (CBM) of 42 compounds, and compare the band gaps of two efficient organic/inorganic perovskite materials. Among them, there are 39 kinds of crystals with their tolerance factors which are between 0.8 and 1.10. Then, the important parameters such as dielectric function, PDOS curve, elastic modulus, shear modulus and Poisson's ratio of the crystals are analyzed. Moreover, we investigate the maximum value of the valence band of these materials and the size of the valence band of the electron transport layer material, hole transport layer and back electrode material. These results are helpful to analyze the efficiency of carrier migration of these materials. We further screen out 14 kinds of crystalline perovskite crystals with energy band gaps between 1.33 eV and 2.40 eV. Then we analyze the minimum edge of the conduction band of these materials (Cs₂ZrI₆, Cs₂HfI₆, Cs₂SeI₆, Cs₂TeI₆, Rb₂HfI₆, Rb₂SeI₆, Rb₂TeI₆, Rb₂GeCl₆ Cs₂SeBr₆ , Rb₂SeBr₆ and Rb₂HfBr₆) is larger than − 4.6 eV, and the maximum value of the valence band lower than − 5.4 eV, refer to Fig. 7. Considering the above factors, we have reduced the range of crystals screened.

From the absorption spectrum, the light absorption curves of Cs₂SeI₆, Cs₂SeBr₆ Rb₂SeI₆, Rb₂SeBr₆ and Rb₂SeCl₆ fall in a relatively large range in the visible light region. The conclusion obtained according to the energy band theory is in well agreement with the absorption spectrum. The stability factor should be further considered, and the formation energy is an important factor for the thermodynamic stability of the reaction. The results show that the formation of compounds containing Cs and Rb elements is roughly the same, and the perovskite crystals containing Cs are more stable. Among several positive tetravalent valence elements, the perovskite crystal structure containing Zr and Hf is more stable. However, its light absorption efficiency is not very high. Taking into consideration with its band gap size and light absorption properties Cs₂SeI₆, Cs₂SeBr₆ and Rb₂SeBr₆ can be used as ideal candidates for light absorption materials. This work can not only guide experiments and design experiments rationally, but also reduce research and development costs, shorten research time, and improve the success rate of material design. Above all, it can provide an important theoretical basis for the industrialization of perovskite solar cells.