# First-Principles Study on Influences of Crystal Structure and Orientation on Band Offsets at the CdS/Cu2ZnSnS4 Interface

Authors: Wujisiguleng Bao; Masaya Ichimura
Journal: International Journal of Photoenergy
 (2012)

---

## Abstract

Cu2ZnSnS4 (CZTS) has attracted much attention recently as an absorber layer material in a heterojunction solar cell. Using the first-principles method, we calculate the band offsets for the CdS/CZTS heterojunction. The valence band offset ΔEv is 1.2 eV for the (001) CdS/CZTS heterointerface and 1.0 eV for the (010) heterointerface, when CZTS is considered to crystallize in the kesterite structure. When CZTS is considered to crystallize in the stannite structure, ΔEv=1.1 eV for the (001) heterointerface and ΔEv=1.3 eV for the (010) heterointerface. In any case, the conduction band minimum of CZTS is higher than that of CdS, and the conduction band offset ΔEc is in a range between 0.1 and 0.4 eV.


---

## Body

## 1. Introduction

Cu2ZnSnS4 (CZTS) is a p-type semiconductor with a band gap of 1.5 eV, which is close to the optimum band gap for solar cells and has high absorption coefficient (>104 cm−1) in the visible region. In addition, the earth abundance and the nontoxic nature of the elements eliminate environmental concerns. Thus, CZTS has attracted much attention recently as an absorber layer material in a heterojunction solar cell [1, 2]. So far, CdS has been selected as a partner of the heterojunction [2, 3]. CdS is an n-type semiconductor with a band gap of 2.4 eV and has also been used for CuInSe2 and CdTe-based solar cells. The solar cells based on CZTS absorber layers and CdS window layers have been developed rapidly, and the highest efficiency obtained so far is about 8.4% [4].Electrical transport properties of heterojunction solar cells are greatly influenced by the band offsets at the heterointerface. The difference in the band gap across the interface is divided into the valence band and the conduction band offsets. These offsets may dramatically influence both the photocurrent and dark current of heterojunction solar cells. The valence band offset for CdS/CZTS has been studied both theoretically and experimentally. Chen et al. obtained the valence band offsetΔEv of 1.01 eV on the basis of the first-principles calculation [5]. On the other hand, Nagoya et al. obtained a significantly smaller value of valence band offset, 0.3~0.57 eV on the basis of the first-principles calculation, considering that the lattice of CdS is strained to fit the CZTS lattice [6]. Haight et al. experimentally obtained a similar value, 0.54 eV, using ultraviolet photoelectron spectroscopy (UPS) [7]. If ΔEv=1.01 eV, the conduction band minimum (CBM) of CZTS is higher than that of CdS; that is, the band alignment is of type II. On the other hand, if ΔEv=0.54 eV, CBM of CZTS is lower than that of CdS and the band alignment is of type I. In a previous work [8], we also carried out the first-principles calculation for the CdS/CZTS heterostructure and obtained ΔEv of 1.2 eV, similar to the results of Chen et al. In the previous calculation, we assumed that the heterointerface is (001) oriented and that CZTS crystallizes in the kesterite structure. It is generally believed that the kesterite is the most stable structure for CZTS, but several theoretical calculations predict that the energy difference from the stannite structure is small compared with thermal energy at the material growth temperature [9, 10]. Thus, we may expect that the actual CZTS film is a mixture of the kesterite and stannite phases. In addition, both the CZTS and CdS films in the solar cell are in fact polycrystalline, and thus the interfaces of various different crystallographic orientations should be formed. Therefore, in this study, we attempt to predict the influences of crystal structure and interface orientation on the band offset by calculating the band offsets for both the kesterite and stannite structures with (001) and (010) interface orientations.


## 2. Calculation

Figure1 shows the supercells of CdS/CZTS heterojunctions considered in this study. CdS is assumed to have the cubic crystal structure. Figure 1(a) shows the (001) CdS/kesterite-CZTS supercell. In the supercell along the (001) orientation, two unit cells of CdS and one unit cell of CZTS are stacked alternately. There are two inequivalent interfaces between CdS and CZTS, and the (Cu, Zn) layer in CZTS is adjacent to one interface, while the (Cu, Sn) layer is adjacent to the other interface. Figure 1(b) shows the (001) CdS/stannite-CZTS supercell. The (Sn, Zn) layer of CZTS is adjacent to one interface and the (Cu, Cu) layer is adjacent to the other interface. Figure 1(c) shows the (010) supercell of CdS/kesterite-CZTS. The two atomic layers in CZTS adjacent to the interfaces contain the same number of atoms (2Cu, Zn, and Sn), and thus the two interfaces are considered to be effectively equivalent. The two interfaces in the (010) CdS/stannite-CZTS (not shown here) are also considered to be equivalent.(a) Supercell of (001)-oriented CdS/kesterite-CZTS. (b) Supercell of (001)-oriented CdS/stannite-CZTS. (c) Supercell of (010)-oriented CdS/kesterite-CZTS.
(a)
(b)
(c)In the calculation, the average value of lattice constants of CdS and CZTS bulk materials was used for the lattice constant along the heterointerface. The lattice spacing in the perpendicular direction was assumed to be the same as in the respective bulk material. The CZTS and CdS layers in the superstructure are considered to simulate the respective bulk material. Thus, the position of each atom was not optimized; that is, the relative positions of atoms in the CdS and CZTS unit cells are the same as in the bulk material. On the other hand, the S atoms at the interface were allowed to move to minimize the total energy.The calculation of density of states (DOS) was performed on the basis of the first-principles, density-functional, pseudopotential method using the PHASE code developed by Institute of Industrial Science, University of Tokyo [11]. The generalized gradient approximation (GGA) is developed by Perdew-Burke-Ernzerhof exchange-correlation energy [12].We first calculated the energy level difference between the valence band maximum (VBM) and reference core levels for CdS and CZTS then calculated the core levels for the CdS/CZTS supercell. The band offsets are given by [13, 14]
(1)ΔEv=ΔEVBM-1CZTS-ΔEVBM-2CdS-ΔEcoreCZTS/CdS,(2)ΔEc=|ΔEg-ΔEv|,
where ΔEVBM-1CZTS(ΔEVBM-2CdS) is the energy separation between the reference core levels and VBM for the CZTS (CdS) bulk, and ΔEcoreCZTS/CdS is the difference in the core level in the CZTS/CdS supercell. This equation is based on the idea that the energy difference between the core level and VBM in the respective bulk material is conserved in the heterostructure. ΔEg is the band gap difference between CdS and CZTS. Since the band gap is generally underestimated in the calculation based on GGA, we used the experimentally determined band gap values (2.4 eV for CdS and 1.5 eV for CZTS).The energy separation between VBM and the reference core level is obtained by calculating the DOS for CZTS and CdS. For CZTS, the S 3s, Sn 4d, and Zn 3d energy levels form a narrow band well separated from the upper continuous valence bands as shown in Figure2(a) and thus may be regarded as the reference core level. For CdS, the S 3s and Cd 4d levels may be regarded as the reference core level, as shown in Figure 2(b). Thus, there are several different choices of the core levels in the above equation, and all of them are examined in the calculation. Both CZTS and CdS have the same anion S, and the S 3p orbital composes the VBM states. The energy of VBM relative to S 3s is 12.8 eV in CZTS and 11.9 eV in CdS. The larger separation between the S 3s level and VBM in CZTS is caused by the repulsion between the S 3p and Cu 3d states, which moves VBM upwards [15].Density of states for (a) Cu2ZnSnS4 and (b) CdS.
(a)
(b)


## 3. Results and Discussion

The calculated energies of the reference core levels and obtained VBM for the respective layer in the CdS/CZTS supercell are shown in Figure3. Figure 3(a) shows the results of (001) CdS/kesterite-CZTS. Layer 0 (or 16) corresponds to the S layer at the (Zn, Cu) interface, and layer 8 corresponds to the S layer at the (Sn, Cu) interface. As shown in this figure, the energy level of a certain kind of atom orbital is not constant but dependent on the position (layer number). Polarization field due to the different atom configurations in the two interface layers, as discussed below, may cause this. As noted above, (1) implies that the energy difference between VBM and each reference level remains the same as in the respective compound. The “Ev” energies denoted by the stars in Figure 3(a) show the predicted VBM positions for both the CdS and CZTS layers. The VBM position of layer 1 was obtained based on the Zn 3d level in that layer, and the VBM position of layer 2 based on the S 3s level, and so on. As noted above, there are two inequivalent interfaces, and for the (Zn, Cu) interface, the valence band offset was estimated as ΔEv=1.1 eV from the VBM positions at layer 1 and layer 15. For the (Sn, Cu) interface, ΔEv was estimated as ΔEv=1.6 eV from the VBM positions of layer 7 and layer 9. The corresponding conduction band offsets are ΔEc=0.2 eV for the (Zn, Cu) interface and ΔEc=0.7 eV for the (Sn, Cu) interface. If the VBM energy is averaged within CdS and CZTS parts, then ΔEv=1.2 eV and ΔEc=0.3 eV.Calculated reference core levels and VBM for the respective layer in the (001)-oriented CdS/CZTS supercell. (a) CdS/kesterite-CZTS. (b) CdS/stannite-CZTS.
(a)
(b)The results for the (001) CdS/stannite-CZTS structure are shown in Figure3(b). Layer 0 (or 16) corresponds to the S layers at the (Zn, Sn) interface and layer 8 to the S layer at the (Cu, Cu) interface. As in Figure 3(a), the energy level of a certain kind of atom orbital is dependent on the position (layer number). The energy level tends to be lower near the (Zn, Sn) interface (layer 0 or 16). The valence band offset was estimated as ΔEv=1.5 eV for the (Zn, Sn) interface and ΔEv=0.6 eV for the (Cu, Cu) interface. The corresponding conduction band offsets are ΔEc=0.6 eV for the (Zn, Sn) interface and ΔEc=-0.3 eV for the (Cu, Cu) interface. If the VBM energy is averaged within each parts, then ΔEv=1.1 eV and ΔEc=0.2 eV.The results for (010) CdS/kesterite-CZTS are shown in Figure4(a). In contrast to the (001) heterostructures (Figure 3), the energy level of each atom orbital is only weakly dependent on the layer number. As noted above, in the (010) structure, the two interfaces are effectively equivalent, and thus no polarization field is formed in each layer. The VBM energy is averaged within each parts, the band offsets are estimated as ΔEv=1.0 eV and ΔEc=0.1 eV. Figure 4(b) shows the results for (010) CdS/stannite-CZTS. The obtained band offsets are ΔEv=1.3 eV and ΔEc=0.4 eV.Calculated reference core levels and VBM for the respective layer in the (010)-oriented CdS/CZTS supercell. (a) CdS/kesterite-CZTS. (b) CdS/stannite-CZTS.
(a)
(b)We obtained the charge of each atom for the respective layers in the (001) supercells of both kesterite and stannite CZTS, to discuss the position dependence of the core levels. The Sn atom has the largest positive charge among the constituent atoms. In the (001) CdS/kesterite-CZTS structure, the net charge is more positive near the (Sn, Cu) interface than near the (Zn, Cu) interface as shown in Figure5(a), which leads to lowering of the VBM position near the (Sn, Cu) interface (layer 8). For the (001) CdS/stannite-CZTS structure, the net charge is more positive near the (Zn, Sn) interface than near the (Cu, Cu) interface, and thus the VBM position is lower near the (Zn, Sn) interface (layer 0 or 16). The difference in ΔEv between the two interfaces may be caused by electric dipole at the interface [16]. For the (001) supercells, since the two interfaces are not equivalent, the dipole charge may be different between the two interfaces. However, we are unable to identify the dipole charge in the analysis so far.Obtained charge of each atom in the (001)-oriented supercell. (a) CdS/kesterite-CZTS. (b) CdS/stannite-CZTS.
(a)
(b)The band alignments obtained in this work are summarized in Figure6. The obtained ΔEv values are in a range of 1.0 eV ~ 1.3 eV and ΔEc 0.1 eV ~ 0.4 eV. Thus, there will be some difference in the band offsets between the different crystal structures and the different orientations. In addition, for a certain crystal structure and orientation, the band offsets may be different depending on the local atom arrangement at the interface, as can be seen from the difference in ΔEv between the two inequivalent interfaces in Figures 3(a) and 3(b). Therefore, the local band offset may fluctuate in a CdS/CZTS heterostructure. However, in any case, CBM of CZTS is higher than that of CdS; that is, the CdS/CZTS heterojunction is of type II. Therefore, there will be no barrier for the flow of photogenerated electrons from CZTS to CdS, and thus the quantum efficiency can potentially be high. However, since VBM of CZTS is close to CBM of CdS, the recombination probability of majority carriers will be high, which will lead to reduced output voltage.Figure 6
Band alignments of CdS and different crystal structures of CZTS (kesterite and stannite) with (001) and (010) orientations.


## 4. Conclusion

We have calculated the band offsets for the CdS/CZTS heterojunction by using the first-principles, density-functional, pseudopotential method for two different crystal structures of CZTS (kesterite and stannite) and two different interface orientations ((001) and (010)). For CdS/kesterite-CZTS, the obtained valence band offsetΔEv=1.2 eV for the (001) orientation and ΔEv=1.0 eV for the (010) orientation. For CdS/stannite-CZTS, ΔEv=1.1 eV for the (001) orientation and ΔEv=1.3 eV for the (010) orientation. In any case, the CdS/CZTS heterojunction is of type II; that is, CBM of CZTS is higher than that of CdS.



---

*Source: 619812-2012-12-03.xml*