First-Principle Calculation of High Absorption-TlGaTe2 for Photovoltaic Application.

We use first-principle calculations based on hybrid functional and the Bethe-Salpeter equation method to investigate the electronic and optical properties of dichalcogenide TlGaTe2. Based on theoretical studies, TlGaTe2 has until recently been considered as an indirect band gap material, however; by employing more accurate hybrid functional model, we showed that although TlGaTe2 has an indirect band gap of 1.109 eV, it also exhibits a fundamental direct band gap of 1.129 eV. Our finding was further confirmed by the optical studies on TlGaTe2, which show that the absorption peak is registered at a photon energy of 1.129 eV. It was also shown that TlGaTe2 has high optical absorption peaks in the visible region. Based on phonon and elastic constant calculations, it was shown that TlGaTe2 is dynamically and mechanically stable. Our findings show that TlGaTe2 is a potential candidate for photovoltaic application.

1. Introduction

Photovoltaics (PV) show great potential in meeting the needs of the world’s energy consumption. Silicon-based solar cell technologies have over 90% of the PV market today, but the research community has been yearning for new materials that are cheap, flexible and efficient in solar cells. Band gap value and absorption coefficient of the material are the two major factors that decide whether the material can be of interest for photovoltaic (PV) applications. We have earlier reported several promising non-silicon materials with intermediate band gaps for photovoltaic applications [1,2]. There is an increasing interest in the family of ternary compounds with the general formula III-III-VI2 (where III = Al, Ga, In, Tl; VI = Se, Te, S)—known as dichalcogenides [3]. In a recent published article, we carried out in-depth analysis of TlBiS2 (thallium-based) and calculated electronic and optical properties for the photovoltaic applications [4]. Thallium-bearing ternary materials have been studied in optoelectronics [5], but not extensively investigated in photovoltaics due to the fact that thallium is toxic. Researchers have shown interest to the Hall effect analysis, the photoconductivity and recombination kinetics in dichalcogenide-TlGaTe2 material [6,7]. From our earlier study on non-silicon solar cell materials [8], we identified that TlGaTe2 has a strong absorption coefficient in the visible region, but a complete study of TlGaTe2 is needed to confirm that it is potential material for PV applications.

In this article, we study structural, electronic, spectroscopic and optical properties of TlGaTe2 dichalcogenide using first-principle calculations. As the density functional theory (DFT) is considered as one of the most effective and accurate computational methods, we investigate the electronic band structure of TlGaTe2 employing DFT analysis. We report accurate electronic band structure of TlGaTe2 based on a hybrid density functional (HSE06) calculations. The calculated HSE06 band structure for this compound helps to identify the type of band gap TlGaTe2 possess and its value, authenticating the suitability for the solar cell application. We also study the structural stability and mechanical stability of TlGaTe2. To verify the Raman active (R-active) and infrared active (IR-active) modes with previously published experimental results [9], we simulated the Raman spectra and infrared spectra of TlGaTe2.

2. Computational Details

Total energies of TlGaTe2 were determined by the projected augmented plane-wave (PAW) utilization of the Vienna ab initio simulation package (VASP) [10]. We used the Perde-Burke-Ernzerhof (PBE) of the generalized gradient approximation (GGA) to achieve the structure relaxations [11]. The ionic coordinates are fully relaxed until the total energy is changed 10−6 eV per atom. Ground-state geometries were calculated by minimizing stresses and Hellman-Feynman forces using the conjugate-gradient algorithm with a force convergence threshold of 10−3 eV Å−1. We used the screened-exchange Heyd-Scuseria-Ernzerhof hybrid function (HSE06) to study the electronic properties of TlGaTe2. The screened parameter was set to 0.2 Å−1 and 30% of the screened Hartree-Fock (HF) exchange was mixed with the PBE exchange functional in the HSE06 method [12]. We applied a 10 × 10 × 10 Γ-centered Monkhorst-pack k-point mesh, which is used in HSE06 self-consistent calculations for TlGaTe2. For the calculation of vibrational properties, the software PHONOPY along with VASP, was utilized [13,14]. For vibrational properties, we used 2 × 2 × 2 supercell (64 atoms). Solving the Casida’s equation is the best approach for calculating the dielectric function [15]. We used a 16 × 16 × 16 -centered grid for both GW (green function with Coulomb interaction) and the Bethe-Salpeter equation (BSE). For all these calculations, a plane-wave cut off of 500 eV was used. The Raman activity was obtained by calculating the polarizability for the vibrational modes [16]. The IR spectra is simulated with the analytic formula of Baroni et al. (using the Born effective charge tensor) [17].

3. Results and Discussions

TlGaTe2 crystallizes in a tetragonal structure (I4/mcm, space group No. 140) as shown in Figure 1a. The calculated Ga-Te distance is 2.69 Å, and Tl-Tl bond length is smaller than the separation between Te and Tl atoms, which make up the GaTe4 tetrahedra (Figure 1a) [18]. The open-circuit voltage and short-circuit current of the solar cells are highly influenced by the band gap of the photoactive semiconductors. The detailed-balance limiting efficiency of an ideal solar cell is 32% for a material with the optimal band gap of 1.4 eV [19]. Band gap calculation using electronic band structure provides a promising opportunity to identify suitable PV materials. The band structure and site-projected density of states of dichalcogenide TlGaTe2 calculated with HSE06 are shown in Figure 1b,c, respectively. The valence band maximum (VBM) is located at the Z-point and Z1-point. The first conduction band minimum (CBM) is located along M-P directions and the second CBM is located at Z-point. In two different experimental studies, it was reported that the band gap energy of TlGaTe2 is 0.84 eV [20] and 1.2 eV [21]. The reason for this deviation is not known. Our calculations show that TlGaTe2 has an indirect band gap of 1.109 eV, which is higher than the mBJ calculation (0.94 eV) [22] but has a good agreement with an experimentally verified band gap value of 1.2 eV [21]. Also, we observe that TlGaTe2 exhibits a direct band gap of 1.129 eV at the Z-point and the difference between direct and indirect band gap is 20 meV. This phenomenon is verified further by a study on the optical properties of TlGaTe2 later in this article. The site-projected density of states shows that valence band maximum (VBM) is mainly due to Tl-s state with little contribution of Te-p state, while conduction band maximum (CBM) is mainly contributed by Ga-s state and Tl-d state in Figure 1c.

Figure 1

(a) Crystal structures of TlGaTe2. The illustration shows the legends for the different type of atoms; (b) calculated HSE06 electronic band structure of TlGaTe2; (c) total and site projected density of states of TlGaTe2. The Fermi level is set to zero. (colour code: red-Ga, green-Te, blue-Tl).

Since TlGaTe2 exhibits narrow indirect band gap at the Z1-point (Figure 1b), the VBM is flatter than the first CBM along M-P directions. Note that around Z-point, the VBM is somewhat flatter than the second CBM. Hence, it would appear that electrons are lighter than holes. The effective masses of holes and electrons have been calculated by the finite difference method [23]. The effective masses of TlGaTe2 calculated by HSE06 are listed in Table 1. For the indirect electronic transition, we found that the effective mass of electrons is 0.196 me (me–a mass of an electron) along the M|X-P direction and the effective mass of holes is 0.736 me along the Z1-M direction. It is interesting to note that Qasrawi et al. used the Hall effect measurement and reported the hole effective mass of TlGaTe2 to be 0.730 me [7], which is in good agreement with our calculated hole effective mass (0.736 me). For the direct transition, we noted that the EM of electrons is 0.378 me and EM of holes is 0.736 me at Z-point.

Table 1

The calculated effective masses of TlGaTe2, where m*hh and m*e are the effective masses of heavy holes and electrons respectively. me is mass of the electron.

Compounds	Effective Masses
	m*hh me	m*e me
TlGaTe2 (indirect electronic transition 1.109 eV)	0.736	0.196
TlGaTe2 (direct electronic transition 1.129 eV)	0.736	0.378

For comparison, the effective mass of electron such as silicon (Si) and gallium arsenide (GaAs) are 0.26 and 0.12 me [24], respectively. The effective mass (EM) of the electron is low for TlGaTe2 compared to Si. Thus, TlGaTe2 has high electron mobility than silicon, which would be useful for an efficient solar cells [25].

In order to understand the dynamical stability of the TlGaTe2, we carried out the phonon calculations. The calculated vibrational spectrum curves together with the corresponding site projected density of states are plotted in Figure 2. The positive phonon frequencies suggest that TlGaTe2 is dynamically stable. The phonon frequencies in the TlGaTe2 have a separation between the low and high-frequency modes (the so-called phonon band gap). The vibrational modes are spread over from 0 to 1735 cm−1. As expected, based on the relationship between masses and frequencies given by equation [26], the heavier element thallium is normally connected with the lower frequencies. The lattice vibrational modes for Tl, Ga, and Te are present in the 0 to 450 cm−1 range, 1167 to 1735 cm−1 range and 0 to 767 cm−1 range, respectively.

Figure 2

Calculated phonon spectra and site projected phonon density of states for TlGaTe2.

To validate the mechanical stability of tetragonal-TlGaTe2, we calculated the single-crystal elastic constant tensor using the finite strain technique. This structure has six independent elastic constants namely C11 = 32.205, C12 = 15.894, C13 = 11.959, C33 = 58.188, C44 = 10.237, and C66 = 12.009. The mechanical stability criteria for the tetragonal (type I) phase are given by [27]:C 2(C C C

All four conditions for the mechanical stability of tetragonal-TlGaTe2 given in Equations (1)–(4) are satisfied, and this finding indicates that tetragonal-TlGaTe2 phase is mechanical stable. From the calculated elastic constants, the bulk (Bv, BR) and shear moduli (Gv, GR) are calculated from Voigt–Reuss–Hill approximations [28,29]. Our calculated values of bulk moduli (B) and shear moduli are 22 GPa and 20.143 GPa from the elastic constants, respectively. The critical G/B value that differs the ductile and brittle materials is 0.5 [30]. As known, when the G/B value of materials is lesser than 0.5, then those materials are ductile—otherwise, they are brittle. The value of G/B for TlGaTe2 is 0.916, which shows that TlGaTe2 is brittle.

According to the crystal symmetry, we have the following modes and the selection rules for IR-active and R-active modes, Γ = A1g (R) + 2A2u (IR) + B1g (R) + B1u + 2B2g (R) + 3Eu (IR) + 3Eg (R). As seen in Figure 3a and Table 2, we observed seven R-active modes in the Raman spectra for TlGaTe2, whereas Gasanly et al. [9] found only three out of seven predicated R-active modes. Our calculated R-active modes confirm the existence of three Eg mode with frequencies of 1Eg= 19 cm−1 (16 cm−1 [31]), 2Eg = 60 cm−1 (67 cm−1 [9]) and 3Eg = 157 cm−1 (165 cm−1 [9]). Our calculated R-active modes confirm the existence of one A1g mode with a frequency of A1g = 128 cm−1 (135 cm−1 [9]) and three B1g mode with frequencies of B1g = 76 cm−1 (76 cm−1 [31]), B1g = 99 cm−1 (99 cm−1 [31]) and 2B2g = 216 cm−1 (210 cm−1 [31]).

Figure 3

(a) Calculated Raman spectra of TlGaTe2; (b) calculated infrared spectra of TlGaTe2.

Table 2

Calculated phonon frequencies ωcal (in this work and previous work [31]) and their experimental values ωexp [9].

TlGaTe2-Phonon Activity	Phonon Frequencies (in cm−1)
	ωcal [31]	ωcal (This Work)	ωexp [9]
R-active 1Eg	16	19	-
R-active 2Eg	60	60	67
R-active 3Eg	152	157	165
R-active B1g	76	76	-
R-active A1g	125	128	135
R-active 1B2g	99	99	-
R-active 2B2g	210	216	-
IR-active 1A2u	8	18	27
IR-active 2A2u	161	170	175
IR-active 1E2u	31	34	44
IR-active 2E2u	77	78	88
IR-active 3E2u	182	188	192

From Figure 3b and Table 2, we observed five IR-active modes in infrared spectra calculation [17]. Our calculated IR-active modes confirm the existence of two A2u mode with frequencies of 1A2u = 18 (27 cm−1 [9]) and 2A2u = 170 (175 cm−1 [9]), along with three Eu mode with frequencies of 1Eu = 34 (44 cm−1 [9]), 2Eu = 78(88 cm−1 [9]) and 3Eu = 188(192 cm−1 [9]). In Table 2, we presented our calculated values for R-active and IR-active modes along with results published earlier in [31] and the experimental values [9]. Note that our calculated R-active and IR-active modes are in better agreement with the experimental results obtained in [9] than [31]. However, the calculated phonon frequencies of TlGaTe2 are lower than the corresponding experimental phonon frequencies. The exchange correlation functional used in this study is not enough to estimate the transition pressure. On the other hand, it is very difficult to obtain a pure form of TlGaTe2 at ambient conditions due to experimental limitations and metastability.

The optical properties have major influence on the solar cell materials; thus, we study the absorption properties of TlGaTe2 material by employing optical dielectric function ε(ω) = ε1(ω) + iε2(ω). The dielectric function is contingent on the wavelength of the electromagnetic radiation, and it is related to the interaction between photons and electrons. The absorption coefficient of the material is highly influenced by the imaginary part ε2(ω), and it can be determined from the inter-band optical transitions by summing over the unoccupied states using the equation [32]. The real part ε1(ω) of dielectric function can be calculated from the ε2(ω) by Kramer-Kronig relationship [33]. The calculated imaginary part of the dielectric function and absorption coefficients of the TlGaTe2 compound are presented in Figure 4. According to the directional dependency of ε1(ω) and ε2(ω), we can conclude that TlGaTe2 is an anisotropic material to the optical properties. Therefore, we present both the real and imaginary parts of the dielectric function of TlGaTe2 in Figure 4 along x and z-directions.

Figure 4

(a) Calculated dielectric function of TlGaTe2; (b) absorption coefficient of TlGaTe2.

In Figure 4a, both the real and the imaginary part of the dielectric function are plotted against the photon energy. The absorption coefficients (α) for TlGaTe2 are calculated numerically with the inclusion of excitonic effects treated within the Bethe-Salpeter equation (BSE) method, which gives numerical results that are in better agreement with the experimental absorption spectra [34,35,36,37]. To achieve BSE calculation, we start from the DFT level, which permits one to obtain a ground state energy of the TlGaTe2. In the second stage, GW (G–Green function and W–Coulomb interactions) approximations are used to calculate the quasi-particle energies. The GW function fails to calculate the excitonic effect between electron-hole interactions. The electron-hole interactions are then included beyond GW approximations, and the BSE method is employed in order to get improved dielectric function for TlGaTe2. To shed light on the accuracy of the BSE method, we calculated the absorption coefficient for silicon and compared with the corresponding experimental results. In Figure 4b, both the numerical and experimental absorption coefficients for Si and TlGaTe2 are depicted in Figure 4b. Note that the calculated absorption coefficient based on the BSE method for silicon is almost similar to the experimental values obtained from [36].

From Figure 4a,b, it is clearly seen that TlGaTe2 exhibits absorption peaks at 1.13 eV along the x-direction, which has to be due to the direct band gap of 1.129 eV at Z-point (from electronic band structure). Figure 4b shows a high absorption coefficient in the visible region. Absorption peaks of TlGaTe2 are observed at 1.129, 2.02, 2.5 and 2.7 eV along both x and z-directions. When the photon energy is about 2.7 eV along the x-direction, and 2.5 eV along the z-direction, the absorption coefficient of TlGaTe2 reached a maximum in the visible region. The valence band maximum is mainly attributed to the states Tl-s, and thallium is a very good absorbing material in the visible region. Consequently, we observed that TlGaTe2 exhibits high absorption in the visible region. Due to the indirect band gap of silicon, the absorption occurs for the photon energies higher than 2.5 eV as shown in Figure 4b. The difference between the indirect and direct band gap in silicon is about 1.40 eV, which results in low absorption coefficients for silicon in the visible region. Our results show that TlGaTe2 has better absorption coefficients compared to silicon in the visible region.

4. Concluding Remarks

In summary, we can mention that the electronic band structure and optical properties of TlGaTe2 with accurate first-principle calculations were investigated. The calculated HSE06 electronic structures confirmed that TlGaTe2 exhibits an indirect band gap, with a value of 1.109 eV, at the Z1-point. However, TlGaTe2 also exhibits a direct band gap, with a value of 1.129 eV at the Z-point. The first absorption peak at a photon energy of 1.129 eV confirms this finding. Until now, numerical limitations in calculating the band gap energies in the range of a few meV difference have been leading to claims stating that the fundamental gap of TlGaTe2 is indirect. Our accurate calculations from employing BSE revealed that TlGaTe2 possesses an additional direct band gap of 1.129 eV with a difference of 20 meV compared to the indirect band gap. Our comparison of the calculated absorption coefficient of silicon and TlGaTe2 shows that TlGaTe2 has a better absorption coefficient than silicon in the visible region.

Our phonon calculations show that TlGaTe2 is dynamically stable as no imaginary frequency was observed. Our elastic constant calculations illustrate that TlGaTe2 is mechanically stable. Raman and IR calculations are carried out for TlGaTe2 and compared with previous experimental and calculated results. Our results seem to be in better agreement with the experimental results than any other calculated results published previously. Our detailed studies of electronic and optical properties of the ternary dichalcogenide material suggest that TlGaTe2 is a potential candidate for photovoltaic applications.