1. Introduction
One of the factors contributing to global warming is the burning of fossil fuels, which increases the concentration of carbon dioxide in the atmosphere, warming the planet and changing the climate [
1,
2,
3]. The Earth absorbs a significant portion of the light emitted by the Sun when it reaches its surface. As a result of heating the planet with this energy, its surface emits infrared radiation. The greenhouse effect, which causes global warming, occurs when carbon dioxide in the atmosphere absorbs much of the incoming thermal radiation and reflects it back to the Earth’s surface [
4,
5]. It is known that electricity generation emits about 40% of greenhouse gases into the atmosphere [
6,
7]. However, the world now obtains 80% of its energy from fossil fuels (oil, coal, and gas), which also contribute to environmental pollution and the daily depletion of these natural supplies [
8]. Thus, in light of the anticipated need for electricity worldwide and the 21st - century green energy agenda, the issue of the development of alternative and renewable energy sources is brought up, and the potential for transforming non-traditional energy into electricity is examined [
9].
Alternative sources such as windmills, moisture-to-electricity converters, thermoelectric (TE) generators, photovoltaic (PV) converters (solar panels), solar thermophotovoltaic converters and the use of geothermal waters have a positive effect on the air environment, but they are not available everywhere. One of the effective natural sources of alternative energy is the synthesis and optimization of the properties of new nanomaterials to create solar energy converters into electricity, which is very effective and environmentally friendly compared to many other methods [
1,
2,
3,
4,
5,
6,
7,
8,
9]. The great growth potential of this alternative energy industry is due to such global factors as the need to ensure national energy security and the rising cost of fossil energy sources. Alternative energy has other unique advantages: for example, solar energy is available to everyone, free, practically inexhaustible, and the process of converting it into electrical energy does not have a negative impact on the environment [
10,
11,
12,
13].
In recent years, photovoltaics has been rapidly developing for the developed of which, in addition to silicon and perovskites, many types of crystalline materials have been developed and proposed [
14,
15,
16,
17,
18,
19,
20,
21]. One of these promising and promising materials is kesterite materials with the general formula A
2BCD
4 (A = Cu, Ag; B = Zn, Cd, Ni, Mg, ...; C = Sn, Ge, Si; D = S, Se) [
22]. Since their emergence as so-called “Third generation generators,” photovoltaic (PV) systems have been hailed as an environmentally and economically viable alternative to conventional technologies for solving the world’s energy, safety and environmental problems [
23,
24,
25]. However, despite the obvious achievements in this area, the development and research of the fundamental properties of potentially new kesterite photovoltaic materials is of great importance for improving the performance of devices based on them.
Because of their direct bandgap in the 1.0–2.5 eV region, the most extensively used kesterite crystals based on the Cu
2ZnSnS
4 and Cu
2ZnSnSe
4 (CZTSSe) system are utilized in the industrial manufacture of solar panels [
26]. Kesterite’s range of use is increased by the adjustable band gap, which also enables task-specific adjustment for ideal spectral matching [
27]. Furthermore, they are excellent candidates for the production of solar energy due to their p-type conductivity and high optical absorption coefficient (>10
4 cm
-1) [
28]. Additionally, this allows for the kesterite film to be thinned down to a point where the solar panel’s cost can be decreased without sacrificing efficiency [
28]. The efficiency of kesterite-based solar cells has recently risen from 12.6% to 13.6% [
29,
30,
31]. However, the performance of solar cells based on this material is still far from the theoretical limit, indicating that the efficiency potential of kesterite is still little exploited. The reason for the poor performance of CZTSSe is mainly due to their high open circuit voltage deficit, which has been repeatedly reported in research papers [
32,
33,
34,
35]. This is due to fluctuations in the band gap and potential induced by crystalline disorder between elements A and B sites of kesterite which occurs at the structural and electronic levels [
32,
33,
34,
35]. Existing problems force researchers to develop new analogues of CZTSSe, however, obtaining new perovskites including some members of the Cu
2NiXS
4 family (X = Si, Ge, Sn) is labor-intensive work due to the complex single-phase growth of kesterite while obtaining a homogeneous and high-quality layer free of secondary phases [
34,
35]. Sometimes the final synthesis of kesterite materials results in many undesirable solid solutions, which complicate the work [
36], and some resulting materials with a kesterite structure may have undesirable properties. Therefore, in recent years, preliminary prediction of properties applied to the synthesis of materials has become an integral tradition among the solid state community.
In this regard, recently the properties of kesterites have also been studied by various theoretical methods, as a result of which the efficiency of solar cells based on them is constantly increasing [
25,
26,
27,
28,
29,
30,
31,
32,
33,
34,
35,
36]. Density functional theory (DFT) is a potent theoretical approach that has gained significant traction in the last ten years as a major tool for the theoretical study of solid materials. Its potent approach accounts for the behavior of electrons in all atomic-molecular environments and offers a highly accurate reformulation of quantum mechanical calculations of solids. This is because contemporary computing clusters can solve the Kohn-Sham equations efficiently [
37,
38,
39,
40,
41,
42,
43,
44,
45]. However, these formulas are predicated on a single estimate, that of the exchange-correlation energy, which accounts for the precision of quantum computations. Here, many of the basic characteristics of compounds based on the Cu
2NiXS
4 (X = Si, Ge, Sn) family still remain poorly studied and correspond to the current topic requiring in-depth research, despite the annual increase in publications devoted to the study of the properties of kesterites.
In this work, using quantum chemical calculations within the framework of density functional theory, the electronic and optical properties of kesterites of the Cu2NiXS4 (X = Si, Ge, Sn) family are studied, a detailed study and disclosure of which is important for the appropriate selection of the synthesized material for specific applications.
3. Results and Discussion
The energy band distribution diagram, the band energy’s dependence on the density of electronic states (DOS), and the band gap’s numerical values are used to evaluate the electronic characteristics of the Cu
2NiXS
4 (X = Si, Ge, and Sn) system.
Table 1 compares the bandgap values we calculated within DFT - HSE06 with the results of experimental measurements by independent authors.
As shown in
Table 1, the results of bandgap calculations of Cu
2NiSnS
4, Cu
2NiGeS
4 crystals obtained from the HSE06 functional are in good agreement with experiment.
The total density of electronic states of Cu
2NiSnS
4, Cu
2NiGeS
4, and Cu
2NiSiS
4 semiconductor crystals were then analyzed using mBJ calculations. The results of calculations of the total density of electronic states for Cu
2NiSnS
4, Cu
2NiGeS
4, and Cu
2NiSiS
4 crystals are shown in
Figure 1a–c.
The results shown in
Figure 2 demonstrate that the band gap lowers and the Fermi levels move towards the valence band when Si is substituted with Ge and Sn. Conversely, it is evident that the density of states drops when Sn takes the place of Si. In this case, new electronic states are formed in the energy gap of Cu
2NiGeS
4 and Cu
2NiSnS
4, which are important from the point of view of using the material in electronic devices.
To justify the change in the bandgap width, it is necessary to analyze the spectra of the optical properties of the materials under study. Calculated optical properties of materials, including their absorption coefficient and refractive index, provide information about what type of response these materials will exhibit when photons are incident on them [
57]. The optical properties of the Cu
2NiXS
4 (X = Si, Ge, Sn) system were investigated based on the calculation of their real (ɛ
1) and imaginary (ɛ
2) parts of the dielectric functions. The real part shows the energy-saving capacity of a material, which is something that is assumed to be inherent in all materials at zero energy or zero frequency limit.
Figure 2a,b shows the curves of ɛ
1 and ɛ
2 versus the energy of incident photons for kesterites of the Cu
2NiXS
4 (X = Si, Ge, Sn) family.
From
Figure 2a it can be seen that at the highest photon energies, all these materials, namely kesterite containing silicon, exhibit metallic behavior. That is, a negative value of the real part indicates the possession of a metallic nature at high energies. This makes it possible to estimate the metallicity fractions of materials using a real function, which shows feedback from the optical band gap. For solar devices, the behavior of these materials indicates the energy gain, whereas the imaginary component of the dielectric function indicates the compound’s absorptive capacity. This provides information about how the material reacts when exposed to electromagnetic radiation [
58,
59,
60,
61]. According to
Figure 2b, the replacement of Si with Ge and Sn leads to an increase in the absorption coefficient of the materials under study in the IR and visible radiation range, which is important for the use of materials in solar panels.
The results displayed in
Figure 3 make it evident that adding germanium and tin in place of silicon raises the Cu
2NiXS
4 (X = Si, Ge, Sn) system’s refractive index (n). In some energy ranges, the refractive index falls drastically below unity after reaching its maximum value. Furthermore, we can deduce from the expression n=c/ϑ that a refractive index value less than one means that the incident radiation’s phase velocity is greater than c, which allows the incident rays to pass through the material and turn it transparent to incoming radiation [
62,
63].
Table 2 shows the static values of ε
1x(0), ε
2z(0) and
n according to the DFT-mBJ-WIEN2k calculations.
Beyond the Fermi level, it is known that photon absorption excites occupied states toward unoccupied states. It is known as “optical conduction” when the photons cross the interband transition and as “interband absorption” when they are absorbed. When light is subjected to an electric field, conductivity is known as optical conductivity.
Figure 4 and
Figure 5 display the computed spectra of the real and imaginary components of the optical conductivity of the systems that are being studied. It is also evident from these spectra that adding Ge and Sn to Si results in improved photoconductivity.
Judging by
Figure 4 and
Figure 5, optical conductivity for the Cu
2NiGeS
4 and Cu
2NiSnS
4 system begins at energies lower than in the case of Cu
2NiSiS
4. From the data shown in
Figure 5, it is clear that all these materials actively absorb light even at low photon energies, namely Cu
2NiSnS
4 is sensitive even to rays with an energy of 0.7 eV. Cu
2NiSiS
4 also begins to be activated in the energy range above 1.3 eV, having the highest photoconductivity when absorbing short-wavelength radiation [
64,
65,
66,
67,
68,
69].