# Influence of Secondary Phases in Kesterite-Cu2ZnSnS4 Absorber Material Based on the First Principles Calculation

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

---

## Abstract

The influence of secondary phases of ZnS and Cu2SnS3 (CTS) in Cu2ZnSnS4 (CZTS) absorber material has been studied by calculating the band offsets at the CTS/CZTS/ZnS multilayer heterojunction interfaces on the basis of the first principles band structure calculation. The ZnS/CZTS heterointerface is of type I and since ZnS has a larger band gap than that of CZTS, the ZnS phase in CZTS is predicted to be resistive barriers for carriers. The CTS/CZTS heterointerface is of type I; that is, the band gap of CTS is located within the band gap of CZTS. Therefore, the CTS phase will act as a recombination site in CZTS.


---

## Body

## 1. Introduction

As candidates for low-cost and nontoxic thin-film solar-cell absorbers, quaternary semiconductor Cu2ZnSnS4 (CZTS) has attracted much attention recently as an absorber layer material in a heterojunction solar cell [1–6]. It is a direct-band gap p-type semiconductor having a 1.5 eV band gap, which is close to the optimum band gap for solar cells [7–9].Because of the increased number of elements in CZTS, their synthesis is more difficult than for binary and ternary semiconductors. The secondary phases in CZTS, such as ZnS, Cu2S, and Cu2SnS3 (CTS), could form depending on the processing conditions [10–12]. Experimentally it is found that working devices are obtained for “Cu-poor” and “Zn-rich” compositions [13]. It has been theoretically shown that ZnS and CTS have relatively small formation enthalpies among the above secondary phases [14].Presence of the secondary phases could affect the photovoltaic performance. It has been speculated that the phase separation of CZTS into ZnS and Cu-Sn-S compounds is one of the factors limiting the efficiency of CZTS solar cells [12, 15–19]. However, the effects of the secondary phases on solar cell performance can in fact be significantly different depending on the band offset with CZTS. There are several possible band alignments at heterointerfaces between CZTS and secondary phases as shown in Figure 1. If the band alignment is of type I, the conduction band minimum (CBM) of the secondary phase is higher (or lower) than that of CZTS, and the valence band maximum (VBM) of the secondary phase is lower (or higher) than that of CZTS, as shown in Figures 1(a) and 1(b), depending on whether the band gap of the secondary phase is larger or smaller than that of CZTS. If the band gap is smaller, the secondary phase will act as a sink for both the minority and majority carriers and facilitate the carrier recombination. Thus the carrier separation is hampered, and the carrier collection efficiency of the cell would be reduced. When the secondary phase has a larger band gap (Figure 1(b)), the band offset should become a barrier for both of the minority and majority carriers in CZTS. Then, the secondary phase just reduces the conductivity but does not enhance carrier recombination rate in the absorber layer. On the other hand, if the band alignment is of type II, both of CBM and VBM for the secondary phase are lower or higher than that of CZTS, as shown as Figures 1(c) and 1(d). In those cases, the secondary phase acts as a sink for only one type of carrier; that is, electrons and holes tend to be separated, and thus the carrier recombination will not be facilitated.Band alignment of CZTS/secondary phase. (a) Type I heterointerface withEg of the secondary phase smaller than that of CZTS. (b) Type I heterointerface with Eg of the secondary phase larger than that of CZTS. (c) Type II heterointerface with CBM and VBM in the secondary phase lower than that of CZTS. (d) Type II heterointerface with CBM and VBM in the secondary phase higher than that of CZTS.
(a)
(b)
(c)
(d)Thus, to discuss the effects of the secondary phase on the solar cell performance, it is essential to determine whether the interface with the secondary phase is type I or II. However, the band alignment at the interface with the secondary phase has never been investigated either experimentally or theoretically so far. In this work, the influence ofsecondary phases of ZnS and CTS in CZTS absorber material is studied by calculating the band offsets at the CTS/CZTS/ZnS multilayer heterojunction interfaces on the basis of the first principles band structure calculation.


## 2. Calculation

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 Tokyo [20]. The valence electron configurations used in our calculations are S (3s2, 3p4), Cu (3d10, 4s1), Zn (3d10, 4s2), and Sn (4d10, 5s2, 5p2). We used the generalized gradient approximation (GGA) for the exchange-correlation interaction. For a given atomic arrangement, the lattice constants and atom positions were optimized to minimize the total energy.


### 2.1. Crystal Structure of Cu2SnS3

The most stable structure of CTS is reported to be the monoclinic structure withCc symmetry [21], as shown in Figure 2. However, because of the lower symmetry, one cannot construct (monoclinic-CTS)/(kesterire-CZTS)/(zincblende-ZnS) supercell. CTS can also have a cubic or tetragonal structure, similar to the zincblende structure [21]. In this case, there is no difficulty due to symmetry mismatch in constructing a hetero supercell with CZTS and ZnS. However, in the zincblende-like structure of CTS, arrangement of Cu and Sn atoms is not perfectly ordered [21]. Since perfect periodicity of atom arrangement is needed for the band structure calculation, we consider the unit cells shown in Figure 3. The (001) layer of the metal sublattice alternately consists of (Cu, Sn), (Cu, Sn), and (Cu,Cu). The positions of Cu and Sn atoms are exchanged or not exchanged between the two consecutive equivalent layers, and thus varieties of unit cells shown in Figure 3 can be considered.Figure 2
Monoclinic structure of Cu2SnS3.Figure 3
Unit cells for Cu2SnS3 with different position of Cu and Sn atoms in the (Cu, Sn) layers.First we calculated the lattice constant for these CTS structures, and the calculated lattice constants are the same for all; that is,a=0.55 nm and c=1.644 nm. On the other hand, the total energy is different, lowest for structure-2. If we take the total energy zero for structure-2, the total energies are 2.5 meV/atom and 1.1 meV/atom for structure-1 and structure-3, respectively. The energy difference is small, and considering limited accuracy of calculation, we could not definitively conclude which is the most stable structure.


### 2.2. Density of State

Figure4 shows the calculated density of states (DOS) for (a) structure-2 (tetragonal structure) of CTS, (b) for zincblende ZnS, and (c) for kesterite CZTS. In DOS for both CTS and CZTS, the states near VBM are derived from interaction of S 3p and Cu 3d orbital, and the lowest conduction band is composed of S 3p and Sn 5s orbitals [21]. The reported value of band gap is 0.9 eV for CTS [21] and 3.7 eV for ZnS [22].Density of states (a) for tetragonal structure of CTS, (b) for ZnS, and (c) for kesterite CZTS.
(a)
(b)
(c)


## 2.1. Crystal Structure of Cu2SnS3

The most stable structure of CTS is reported to be the monoclinic structure withCc symmetry [21], as shown in Figure 2. However, because of the lower symmetry, one cannot construct (monoclinic-CTS)/(kesterire-CZTS)/(zincblende-ZnS) supercell. CTS can also have a cubic or tetragonal structure, similar to the zincblende structure [21]. In this case, there is no difficulty due to symmetry mismatch in constructing a hetero supercell with CZTS and ZnS. However, in the zincblende-like structure of CTS, arrangement of Cu and Sn atoms is not perfectly ordered [21]. Since perfect periodicity of atom arrangement is needed for the band structure calculation, we consider the unit cells shown in Figure 3. The (001) layer of the metal sublattice alternately consists of (Cu, Sn), (Cu, Sn), and (Cu,Cu). The positions of Cu and Sn atoms are exchanged or not exchanged between the two consecutive equivalent layers, and thus varieties of unit cells shown in Figure 3 can be considered.Figure 2
Monoclinic structure of Cu2SnS3.Figure 3
Unit cells for Cu2SnS3 with different position of Cu and Sn atoms in the (Cu, Sn) layers.First we calculated the lattice constant for these CTS structures, and the calculated lattice constants are the same for all; that is,a=0.55 nm and c=1.644 nm. On the other hand, the total energy is different, lowest for structure-2. If we take the total energy zero for structure-2, the total energies are 2.5 meV/atom and 1.1 meV/atom for structure-1 and structure-3, respectively. The energy difference is small, and considering limited accuracy of calculation, we could not definitively conclude which is the most stable structure.


## 2.2. Density of State

Figure4 shows the calculated density of states (DOS) for (a) structure-2 (tetragonal structure) of CTS, (b) for zincblende ZnS, and (c) for kesterite CZTS. In DOS for both CTS and CZTS, the states near VBM are derived from interaction of S 3p and Cu 3d orbital, and the lowest conduction band is composed of S 3p and Sn 5s orbitals [21]. The reported value of band gap is 0.9 eV for CTS [21] and 3.7 eV for ZnS [22].Density of states (a) for tetragonal structure of CTS, (b) for ZnS, and (c) for kesterite CZTS.
(a)
(b)
(c)


## 3. Results and Discussion

To evaluate the band offset, we consider a supercell consisting one unit cell of CTS, one unit cell of CZTS, and two unit cells of ZnS. For CTS, the (001) layer of the metal sublattice consists of (Cu, Sn) or (Cu, Cu). For kesterite CZTS, the (001) layer of the metal sublattice consists of (Cu, Sn) or (Cu, Zn). Thus, there are six different CTS/CZTS/ZnS supercells with different arrangements of the interface layer atoms as shown in Table1. Interface 1 is between CTS and CZTS, interface 2 between CZTS and ZnS, and interface 3 between ZnS and CTS. The lattice constant parallel to the interface is assumed to be equal to that of the dominant material CZTS; that is, a=0.551 nm. The lattice spacing in the perpendicular direction is assumed to be the same as in the respective bulk materials.Table 1
Valence band offsets (eV) at each interface and the total energies (meV/atom) for six supercells with different interface layer atoms. The interface between CTS and CZTS is interface 1,  interface between CZTS and ZnS is interface 2, and interface between ZnS and CTS is interface 3.
CTS-CZTS
CZTS-ZnS
ZnS-CTS
Δ
E
v
1
Δ
E
v
2
Δ
E
v
3
Δ
E
total
Interface 1
Interface 2
Interface 3
(1)
(Sn, Cu)-(Zn, Cu)
(Sn, Cu)-(Zn, Zn)
(Zn, Zn)-(Cu, Cu)
0.5
0.3
0.8
14.4
(2)
(Sn, Cu)-(Sn, Cu)
(Zn, Cu)-(Zn, Zn)
(Zn, Zn)-(Cu, Cu)
−0.1
0.1
0
27.9
(3)
(Sn, Cu)-(Zn, Cu)
(Sn, Cu)-(Zn, Zn)
(Zn, Zn)-(Sn, Cu)
0.4
1.4
1.8
0
(4)
(Sn, Cu)-(Sn, Cu)
(Zn, Cu)-(Zn, Zn)
(Zn, Zn)-(Sn, Cu)
0.1
0.8
0.9
5.2
(5)
(Cu, Cu)-(Zn, Cu)
(Sn, Cu)-(Zn, Zn)
(Zn, Zn)-(Sn, Cu)
0
1.8
1.8
27.7
(6)
(Cu, Cu)-(Sn, Cu)
(Zn, Cu)-(Zn, Zn)
(Zn, Zn)-(Sn, Cu)
−0.2
0.8
0.6
13.2Table1 lists the calculated valence band offsets and total energies for the six supercells with different arrangements of the interface layer atoms. The VBM energy averaged in each layer was used to calculate the band offset. ΔEv1 is the band offset at CTS/CZTS interface, and a positive value indicates that VBM of CTS is higher than that of CZTS. ΔEv2 is the band offset at CZTS/ZnS interface, and a positive value indicates that VBM of CZTS is higher than that of ZnS. ΔEv3 is the band offset at ZnS/CTS interface, and a positive value indicates that VBM of CTS is higher than that of ZnS.The total energy per atom is given with the energy of supercell (3) taken as zero. From Table1, one can see that the total energies are relatively small for supercells (3) and (4). In these cells, the average value of valence electrons is near two for the four adjacent metal atoms at all of the three interfaces, and thus the octet rule is almost fulfilled. In contrast, in supercell (5), the (Cu, Cu)-(Zn, Cu) interface has a high energy, and in supercell (2), both the (Sn, Cu)-(Sn, Cu) and (Zn, Zn)-(Cu, Cu) interfaces have a relatively high energy, because of violation of the octet rule.Because the total energy is relatively lower for supercells (3) and (4) as listed in Table1, the band alignment for supercell (3) and supercell (4) is shown in Figure 5. VBM of CZTS is below that of CTS (ΔEv=0.4 and 0.1 eV), and VBM of ZnS is below that of CZTS (ΔEv=1.4 and 0.8 eV). Considering that the band gaps are 3.7 eV for ZnS, 0.9 eV for CTS, and 1.5 eV for CZTS, CBM of CTS is below that of CZTS (ΔEc=0.2 and 0.5 eV), and CBM of CZTS is below that of ZnS (ΔEc=0.8 and 1.4 eV).Figure 5
Predicted band alignments for supercell (3) and supercell (4) of CTS/CZTS/ZnS.Barkhouse et al. measured the band offset at ZnS/CuZnSn(S0.4Se0.6)4 by laser pump/probe ultraviolet photoemission spectroscopy [23]. Cu2ZnSn(S0.4Se0.6)4 has a band gap of 1.2 eV, smaller than that of CZTS by 0.3 eV, and the measured value of ΔEv is 1.3 eV. The present results for ZnS/CZTS seem not inconsistent with their results.Thus, the band alignment of ZnS/CZTS interface is of type I, and there is large spike (>0.8 eV) for both the conduction and valence bands. Such a large spike will create a high resistance barrier to carrier flows. Thus, the ZnS phase may hinder the flow of photo-excited carriers (electrons), reducing the photo current, and impede the flow of majority carriers (holes), increasing the series resistance [17, 23, 24]. On the other hand, since ZnS can expel both the carriers, the recombination rate in ZnS can be low. In addition, since the lattice constant of ZnS (0.541 nm) is close to that of CZTS (a=0.551 nm), the interface will not have high density of defects. Therefore, if the surface or boundary of CZTS grains is covered with ZnS, recombination at the surface/boundary will be suppressed. In fact, Mendis et al. found that the ZnS/CZTS has a low recombination velocity [17].As shown in Figure5, the band alignment of CTS/CZTS interface is also of type I. CTS can be regarded as a cluster of (SnZn+2CuZn) defect pairs in CZTS. The formation of (SnZn+2CuZn) defect pairs can shift VBM upward because the formation of CuZn enhances the p-d repulsion and shift CBM downward because the formation of SnZn makes CBM wavefunction more localized on the more electronegative Sn sites [18, 21]. Therefore, CBM and VBM of CTS are located within the forbidden energy region of CZTS. In this case, since CTS has a smaller band gap, both of minority and majority carriers will flow to CTS. Thus the CTS phase in CZTS will act as a recombination center, lowering both output current and voltage. If the CZTS surface is covered with CTS, CTS is expected to act as a sink for carriers and enhance the surface recombination velocity. However, according to Mendis et al., the CTS/CZTS interface also has a low recombination velocity [17]. Thus, their results cannot be interpreted on the basis of the band alignment but could be related with other factors such as space charge distribution near the interface.The present results show that the CTS phase is definitely detrimental for the solar cell performance. Thus, Zn-rich composition is preferable to prevent formation of CTS, although excessively Zn-rich composition will result in dominance of the highly resistive ZnS phase and reduce the solar cell efficiency.In this work, we consider only the (001) interface. In fact, the band offset may depend on the interface orientation [15]. The actual CZTS layers in solar cells are polycrystalline, and thus the interface with the secondary phases can have various orientations. Therefore, for more accurate prediction of the band offset, the calculation had better be done for the other orientations. However, for the orientations other than (001), we need to consider a much larger supercell. For example, the supercell needed for the calculation of the (010) interface is six times larger than that for the (001) orientation. Since we cannot perform the calculation for such a large cell because of the limited computational power, we consider only the (001) orientation. In general, the band offset is expected to be determined by the two factors. One is the position of VBM with respect to an absolute reference level, and the other is the interface charge-dipole, which causes discontinuity in the electrostatic potential. The second one will depend on the interface orientation, but the first factor is common for any interface orientation. Thus, we may expect that the calculation for only one orientation is still useful to qualitatively predict the band offset.


## 4. Conclusions

In this study, we have calculated the band offsets at the CTS/CZTS/ZnS heterojunction interfaces on the basis of the first principles band structure calculation. The ZnS/CZTS heterointerface is of type I, and since ZnS has a larger band gap than that of CZTS, the ZnS phase in CZTS is predicted to be resistive barriers for carriers. The CTS/CZTS heterointerface is also of type I, and CBM and VBM of CTS are located within the band gap of CZTS. Therefore, the CTS phase will act as a recombination site in CZTS. It should be noted that comparison to experiment is limited by the limits of computing power and that there may be some discrepancies due to inability to carry out more generalized computational modeling.



---

*Source: 592079-2015-01-06.xml*