Insights into Intrinsic Defects and the Incorporation of Na and K in the

Jan 6, 2016 - IBM, Hechtsheimer Str. 2, 55131 Mainz, Germany. § Max Planck Institute for Chemical Physics of Solids, Nöthnitzer Str. 40, 01187 Dresd...
2 downloads 4 Views 690KB Size
Subscriber access provided by UNIV OF NEBRASKA - LINCOLN

Article

Insights into Intrinsic Defects and the Incorporation of Na and K in the CuZnSnSe Thin-Film Solar Cell Material from Hybrid-Functional Calculations 2

4

Elaheh Ghorbani, Janos Kiss, Hossein Mirhosseini, Markus Schmidt, Johannes Windeln, Thomas D. Kuehne, and Claudia Felser J. Phys. Chem. C, Just Accepted Manuscript • DOI: 10.1021/acs.jpcc.5b11022 • Publication Date (Web): 06 Jan 2016 Downloaded from http://pubs.acs.org on January 11, 2016

Just Accepted “Just Accepted” manuscripts have been peer-reviewed and accepted for publication. They are posted online prior to technical editing, formatting for publication and author proofing. The American Chemical Society provides “Just Accepted” as a free service to the research community to expedite the dissemination of scientific material as soon as possible after acceptance. “Just Accepted” manuscripts appear in full in PDF format accompanied by an HTML abstract. “Just Accepted” manuscripts have been fully peer reviewed, but should not be considered the official version of record. They are accessible to all readers and citable by the Digital Object Identifier (DOI®). “Just Accepted” is an optional service offered to authors. Therefore, the “Just Accepted” Web site may not include all articles that will be published in the journal. After a manuscript is technically edited and formatted, it will be removed from the “Just Accepted” Web site and published as an ASAP article. Note that technical editing may introduce minor changes to the manuscript text and/or graphics which could affect content, and all legal disclaimers and ethical guidelines that apply to the journal pertain. ACS cannot be held responsible for errors or consequences arising from the use of information contained in these “Just Accepted” manuscripts.

The Journal of Physical Chemistry C is published by the American Chemical Society. 1155 Sixteenth Street N.W., Washington, DC 20036 Published by American Chemical Society. Copyright © American Chemical Society. However, no copyright claim is made to original U.S. Government works, or works produced by employees of any Commonwealth realm Crown government in the course of their duties.

Page 1 of 20

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

Insights Into Intrinsic Defects and the Incorporation of Na and K in the Cu2ZnSnSe4 Thin-Film Solar Cell Material from Hybrid-Functional Calculations Elaheh Ghorbani,†,‡ Janos Kiss,¶ Hossein Mirhosseini,∗,¶ Markus Schmidt,‡ Johannes Windeln,§ Thomas D. Kühne,k and Claudia Felser¶ †Institute of Inorganic and Analytical Chemistry, Johannes Gutenberg University of Mainz, Staudinger Weg 9, 55122 Mainz, Germany ‡IBM, Hechtsheimer Str. 2, 55131 Mainz, Germany ¶Max Planck Institute for Chemical Physics of Solids, Nöthnitzer Str. 40, 01187 Dresden, Germany §Wilhelm Büchner Hochschule, Ostendstraße 3, 64319 Pfungstadt, Germany kDepartment of Chemistry and Institute for Lightweight Design with Hybrid systems, University of Paderborn, Warburger Str. 100, 33098 Paderborn, Germany E-mail: [email protected] Abstract We have performed density functional theory calculations using the HSE06 hybrid functional to investigate the energetics, atomic, and electronic structure of intrinsic defects as well as Na and K impurities in the kesterite structure of the Cu2 ZnSnSe4 (CZTSe) solar cell material. We found that both Na and K atoms prefer to be incorporated into this material as substitutional defects in the Cu sublattice. At this site

1

ACS Paragon Plus Environment

The Journal of Physical Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

highly stable (Na-Na), (K-K), and (Na-K) dumbbells can form. While Na interstitial defects are stable in CZTSe, the formation of K interstitial defects is unlikely. In general, the calculated formation energies for Na-related defects are always lower compared to their K-related counterparts. Based on thermodynamic charge transition level calculations we can conclude that the external defects are harmless except NaSn and KSn . These defects induce gap states that might be detrimental for the device performance. Keywords: CZTSe, Thin-film solar cell, Group-I dopant

Introduction The global energy demand is constantly on the rise. 17 If renewable-energy resources like solar cells are to make a significant contribution to this growing demand, they should be able to produce energy in the terawatt scale. However, current solar cells with high conversion efficiency such as single-crystal Si cells and thin-film solar cells based on CuIn1−x Gax Se2 (CIGSe) have natural limitations that prevent them reaching the terawatt production scale. Namely, for CIGSe, with the world record efficiency of 21.7%, 13,19 the low attainability and high expense of In and Ga make the large-scale production of CIGSe-based cells unpractical. For the manufacturing of single-crystal Si cells, tremendous amount of energy and natural resources are needed as well. 2 Hence, in recent years there is a constantly growing interest and ongoing scientific discussion to look for alternative photovoltaic materials containing earth-abundant elements. A promising alternative to the CIGSe ternary compound is quaternary CZTSe. Instead of In and Ga, CZTSe contains zinc (Zn) and tin (Sn) that are earth-abundant, inexpensive, and non-toxic elements. Although so far CZTSe-based solar cells are not as efficient as CIGSe-based solar cells, the conversion efficiency of 12.6% for CZTSe-based devices 37 shows a great potential of CZTSe. To further boost the efficiency of CZTSe cells, it is crucial to understand the effects of various defects on this material. By widening our knowledge about the defect physics, we can assess which defects are detrimental for the performance of the 2

ACS Paragon Plus Environment

Page 2 of 20

Page 3 of 20

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

cell, so that the formation of such harmful defects could be prevented during the synthesis. For CIGSe cells it is already well established that the incorporation of a small amount of Na improves the p-type conductivity and conversion efficiency. 32,38 It was commonly accepted that doping with Na is the best choice for improving the performance of CIGSe cells. However, very recently it was found that K impurities also have a benign effect on the CIGSe-device performance. 7,33,34,39 In addition to the diffusion of Na and K from the soda-lime glass substrate into CIGSe, alkaline metals can be incorporated into the absorber layers via the so-called post deposition treatment (PDT). During the PDT a small amount of NaF and/or KF salt is deposited onto the CIGSe layer after its growth has been finished. The PDT unveiled the remarkable role of K in improving the quality of the pn-junction by reducing the electron-hole recombination at the CdS/CIGSe interface. 7,31 Also, the PDT with additional KF shortens the deposition period of CdS, which results in a thinner CdS buffer layer without weakening the cell performance. 7,9,35 Due to the structural similarity between CIGSe and CZTSe (see Fig. 1) one expects that the conclusions obtained from vastly investigated CIGSe devices can also be transferred to CZTSe cells. Considering these regards, one expects that sodium and potassium play a favorable role in solar cells based on CZTSe as well. Therefore, in this work we study the formation energies of Na and K impurities in kesterite CZTSe 5 via hybrid functional calculations. Based on our results, we discuss the structure and energetics of the most likely and unlikely positions in CZTSe for the incorporation of Na and K extrinsic defects. Beside the Na and K point defects, we have also investigated the formation of defect complexes, such as (Na-Na), (Na-K) and (K-K) dumbbells.

3

ACS Paragon Plus Environment

The Journal of Physical Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Computational methodology Computational details In the present work, all calculations have been performed within the framework of density functional theory (DFT), using the VASP 22,23 code. We employed the projector augmented wave method 4,24 with a plane-wave cut-off energy of 400 eV. The exchange-correlation potential has been treated with the HSE06 hybrid functional, 16 which provides a rather good description of the formation energies, atomic and electronic structure of semiconductors. There is a recent theoretical work 27 that tackles the issue of Li, Na and K extrinsic defects in CZTSe via GGA (PBE) calculations. 30 It is known that GGA functionals severely underestimate the band gap of CZTSe. The results achieved by PBE, therefore, are affected by the shortcomings of this functional. Regarding the experimental band gap of CZTSe, there are different values published in the literature. The discrepancy between the values are either due to the different methods that have been used to determine the gap or due to the sample preparation, i. e. samples might contain secondary phases like ZnSe and Cu2 SnSe3 . 1 One group of experimental studies reports a direct band gap in the range of 1.4-1.6 eV for CZTSe. 3,20,26,28 Other works report an experimental gap in the order of 0.8-1.0 eV. 1,10,14 In our work, the contribution of the

Cu

In, Ga

Se

Sn

Zn

Figure 1: Comparison between the conventional tetragonal unit cell of CIGSe in the chalcopyrite structure (left panel) and CZTSe in the kesterite structure (right panel). 4

ACS Paragon Plus Environment

Page 4 of 20

Page 5 of 20

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

Hartree-Fock exchange in HSE06 calculations has been set to 27%, so that the calculated band gap of 1.0 eV is in good agreement with the latter experimental findings and previous HSE06 hybrid functional calculations 5 The fully ordered crystal structure of the conventional tetragonal cell of CZTSe is shown in Fig. 1. From experimental works it is established that there is intermixing between Cu and Zn atoms in CZTSe. 8 In our calculations, however, we have used the fully ordered structure, in which this disorder was not taken into account. The CZTSe bulk was represented by employing a large 216-atom supercell built based on the primitive cell. In these supercells a point defect is separated by about 21 Å from its own periodic replica within periodic boundary conditions. The Brillouin zone integration was performed on a Γ-centered 2×2×2 mesh of k-points. To calculate the atomic and electronic structure of various defects, first we have relaxed the atomic positions and equilibrium volume of the defect-free CZTSe bulk. The formation energies of the point defects and impurities have been calculated using the theoretically obtained equilibrium volume for the bulk. In all calculations convergence was assumed when the largest residual force component dropped below 0.1 eV/Å. The calculated a = 5.73 Å, c = 11.42 Å cell parameters, tetragonal elongation η= c/2a= 0.996 and band gap of Eg = 1.0 eV are in good agreement with the experimental results. 11

Formation energy calculation Defect formation energies are calculated via the supercell approach. In this approach defects are surrounded by a large number of CZTSe atoms, and the whole cell is repeated periodically. 36 The intrinsic point defects were created by removing a Cu, Zn, Sn or Se atom. The Na and K substitutional defects were created when an atom on one atomic site was substituted by a Na or a K impurity. The interstitial defects were created by placing one impurity atom on an interstitial site.

5

ACS Paragon Plus Environment

The Journal of Physical Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 6 of 20

The defect formation energies are calculated as q Ef [dq ] = E[dq ] − E[p] + Ecorr

(1)

− Σi ni µi + q[EVBM + µe + ∆v0/b ], where dq stands for a defect in charge state q. E[dq ] denotes the total energy of the supercell containing the impurity or defect d and E[p] is the total energy of the pure defect-free charge neutral bulk CZTSe. ni represents the number of atoms of type i (host atoms or impurity atoms) that have been added to (ni > 0) (like Na and K) or removed from (ni < 0) (like Cu, Zn, Sn and Se) the supercell when a defect or an impurity was created, and µi is the equivalent chemical potential of the elements in their native elemental state. EVBM is the valence band maximum (VBM) of the defect-free bulk and µe represents the position of the Fermi energy relative to EVBM . The thermodynamic charge transition levels ǫd (q/q′) were computed as the Fermi energy in Eq. 1 at which the charge state of defect d transforms from q to q′. Due to the finite-size of the supercells the formation energies are corrected by two terms: q q∆v0/b and Ecorr , where ∆v0/b is the difference between the average electrostatic potential of q the bulk and the supercell containing the associated neutral defect. 21 Ecorr accounts for the

interaction of a charged defect with its periodic replicas. In the current work we used the q approach proposed by Lany and Zunger 25 to calculate Ecorr .

RESULTS and DISCUSSION Intrinsic defects in CZTSe The removal of a Cu, Zn, Sn or Se atom from the CZTSe lattice creates intrinsic point defects. The formation of such defects in CZTSe and the isostructural sulfide compound (Cu2 ZnSnS4 ) have been studied. 6,15 In agreement with the aforementioned work, 6 our results show that

6

ACS Paragon Plus Environment

Page 7 of 20

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

among intrinsic point defects in CZTSe the copper vacancy VCu has the lowest formation energy for the whole range of the chemical potential of the electrons, i. e. under both p- and n-type conditions (see Fig. 2). By computing the thermodynamic charge transition levels, we found that the only stable charge state for VCu is -1, i. e. this vacancy does not induce any defect level within the band gap and is not detrimental for the carrier transport in CZTSe. In contrast to VCu , the tin vacancy VSn creates three levels: ǫ(-2/-) at EVBM +0.08 eV, ǫ(-3/-2) at EVBM +0.27 eV, and ǫ(-4/-3) at EVBM +0.7 eV. As it is shown in Fig. 2, the formation energy of VSn falls rapidly by increasing µe , i. e. when CZTSe becomes n-type. Considering that VSn has the highest formation energy among native vacancies (4.21 eV for charge state -1 at the VBM), such defects should have a negligible concentration in a p-type CZTSe material. Although the creation of the ǫ(-4/-3) and ǫ(-3/-2) level of VSn could have negative impact on the efficiency of CZTSe thin-film solar cells, this detrimental effect is controlled by the conditions of the CZTSe deposition process, where films are grown p-type and the probability of the formation of Sn vacancies is low. The Zn vacancy VZn , similar to VCu , does not induce transition levels within the band gap. The formation energy of VZn is rather high compared to VCu but decreases by increasing µe (see Fig. 2). For the anion vacancy VSe we have calculated the formation energy of 1.91 eV for charge state +1 at the VBM. The presence of VSe in the material induces a charge transition level ǫ(0/+) which is located at EVBM +0.05 eV and is not detrimental for the carrier transport.

Formation energy of Na and K defects After investigating the intrinsic defects, now we turn our attention toward the discussion on the incorporation of Na and K impurities into CZTSe. To assess the relative stability of Na- and K-related defects, we have computed the formation energies of various defects in different charge states using Eq. 1. Interstitial and substitutional defects are formed by inserting Na and K into the CZTSe bulk and substituting the constituent atoms with Na 7

ACS Paragon Plus Environment

The Journal of Physical Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

and K, respectively. Our results are summarized in Fig. 3 that shows the formation energies as a function of the chemical potential of the electrons µe . Na and K atoms can occupy two different interstitial positions with different chemical environments: a Na or K impurity in an interstitial position is either octahedrally coordinated (oc) by three Cu, two Sn, and one Zn cations (see the left panel in Fig. 4), or tetrahedrally coordinated (tc) by two Cu, one Sn, and one Zn atom (as shown in the right panel of Fig. 4). Our results reveal that for both Na and K impurities, it is energetically more preferable to be located in the tc position rather than in the oc site. For Na impurities in CuInSe2 similar results have been found. 12,29 Both Na and K interstitials have the charge state +1 over the whole range of µe . In p-type CZTSe the respective formation energy for Na+ oc and Na+ tc are -1.14 eV and -1.41 eV, so they are both likely to form. Conversely, due to the larger + size of the K impurity, K+ oc and Ktc interstitials have considerably high formation energies

of 1.08 eV and 0.57 eV, respectively. This suggests that in contrast to Na, K interstitials are highly unlikely to be formed in CZTSe. Next, we have studied the formation of Na and K substitutional defects. On the basis of our results we can state that both Na and K prefer to be incorporated into CZTSe in the

Figure 2: The calculated formation energies of native vacancies in Cu2 ZnSnSe4 as a function of the chemical potential of the electrons µe , where the left (µe =0) and right side (µe =1.0 eV) of the figure corresponds to p-type and n-type conditions, respectively. The positive and negative numbers next to the lines indicate the charge state.

8

ACS Paragon Plus Environment

Page 8 of 20

Page 9 of 20

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

copper sublattice as Na0Cu and K0Cu substitutional defect. The defects are charge neutral, since both Na and K are monovalent, and they are substituting a monovalent Cu atom. For Na0Cu and K0Cu , our computed formation energies are -1.56 eV and -0.23 eV, respectively. The negative formation energies of these defects indicate that at low concentrations it is energetically favorable to incorporate Na and K impurities into VCu defects. Moreover, since the formation energy of Na0Cu is more than 1.3 eV lower than K0Cu , if a Cu site is already occupied by Na, then K most likely cannot kick out Na to take its place. For the CuInSe2 light absorber we have obtained qualitatively similar results. 12 Given the high formation energy of 1.58 eV for K− Zn , the formation of a K substitutional defect in a Zn position is not likely. In contrast to K, the formation energy of 0.34 eV computed for Na− Zn suggests that such substitutional defects might form. The Sn-related

Figure 3: The calculated formation energies of Na- and K-related defects in Cu2 ZnSnSe4 as a function of the chemical potential of the electrons µe . 9

ACS Paragon Plus Environment

The Journal of Physical Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Figure 4: Atomic structure showing the local coordination of octahedrally coordinated interstitial (left panel) and tetrahedrally coordinated interstitial (right panel). The Na or K interstitial atom is represented by a red sphere, and for the rest of the atoms the same color scheme has been used as in Fig. 1. For clarity, only half of the conventional tetragonal cell is shown. defects can have a wide range of stable charge states and their formation energy drops rapidly as µe increases. Since the light absorber layers are grown under p-type conditions, the role of NaSn and KSn should be negligible in the films used for actual solar cells. Regarding the Se site, our data shows that Na and K-related substitutional defects at the Se site in CZTSe have similar formation energies computed for substitutional defects in the CuInSe2 thin film solar cell material. 12,29 The qualitatively and quantitatively similar results clearly show that Se-related substitutional defects are energetically rather unfavorable both in CZTSe and CuInSe2 . In addition to the point defects, we have also looked into the formation of three types of dumbbells in CZTSe, namely (Na-Na), (Na-K) and (K-K) dumbbells in a copper site. + + Our data presented in Fig. 3 show that the formation of (Na-Na)+ Cu , (Na-K)Cu , and (K-K)Cu

dumbbells with the respective formation energies of -3.18 eV, -1.95 eV and -0.66 eV is highly probable in CZTSe, similar to CuInSe2 . 12,29 The (Na-Na) dumbbells are the most stable defects over the whole range of the chemical potential of the electrons µe . For (Na-K) and (K-K) as the system becomes more n-type, the Na0Cu and K0Cu substitutional defects will become more stable. Comparing the formation energies of a (Na-Na)+ Cu dumbbell (-3.18 eV) with two Na in two distinct and spatially well-separated copper positions (-2.84 eV) unveils that for the p-type system it is more favorable to cluster two Na atoms together, probably to lower the strain in the material. A comparison between the formation energy of a (K-K)+ Cu dumbbell and two 10

ACS Paragon Plus Environment

Page 10 of 20

Page 11 of 20

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

KCu point defects reveals that K atoms behave in the same way as Na: K impurities also have the preference to form dumbbells instead of occupying two distinct substitutional Cu positions. Concerning the charge state of the Na- and K-related defects, as one would expect, 3− acceptor defects such as Na− Zn and NaSn form more easily in a n-type material, where µe is

at the CBM. Conversely, donor defects such as Na+ tc form more easily in a p-type material where µe is at the VBM. We note in passing, that regarding the most stable charge state of dumbbells, we found a qualitative difference between CuInSe2 and CZTSe. Namely, while in the former compound dumbbells with the charge state +2 are stable 12 at the VBM, in CZTSe the +1 charged dumbbells are stable for the whole range of µe .

Band structure and thermodynamic charge transition levels In this section we discuss the effect of defects upon the electronic structure of CZTSe. Since the band structure calculations with hybrid functionals are computationally extremely demanding, these calculations have been performed for a smaller 64-atom supercell employing the same methodology as described in the computational details. Fig. 5 shows the band structure of the defect-free CZTSe bulk and native defects. While the Sn vacancy creates a gap state, other intrinsic defects induce no gap state. The charge transition level calculations for the large 216-atom supercell show that Se and Sn vacancies create gap states but Cu and Zn vacancies do not induce gap states. We note in passing that the discrepancy between band structure calculations and charge transition levels can have several origins. One reason is the difference between the supercell sizes. Band structure calculations are performed for 64-atom supercells, where there is a considerable interaction between a defect and its periodic replica. For comparison, the charge transition levels are calculated for supercells with 216 atoms. The other reason has to do with the position of the orbitals that might be placed wrong by DFT. 18 3− 3− − Fig. 6 depicts the band structures calculated for the Na0Cu , K0Cu , Na− Zn , KZn , NaSn , KSn ,

11

ACS Paragon Plus Environment

The Journal of Physical Chemistry

Cu2ZnSnSe4

-

-

0

VCu

+

VSn

VZn

VSe

4 3 Energy (eV)

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 12 of 20

2 1 0 -1 -2 Z

Γ

XZ

Γ

XZ

Γ

X Z

Γ

XZ

Γ

X

Figure 5: Band structure for the bulk CZTSe and native defects in CZTSe. The occupied bands are drawn in black and the unoccupied band are shown in blue. The red dashed line represents the defect level. The dotted lines indicate the position of the Γ point in the Brillouin zone. 0

-

NaCu

-

3-

NaZn

+

NaSe

NaSn

Natc

4 3 2 1 0 -1 -2 Z

Γ

XZ

0 KCu

Γ

X Z

Γ

XZ

3KSn

KZn

Γ

XZ

KSe

Γ

X

+ Ktc

4 3 2 1 0 -1 -2 Z

Γ

XZ

Γ

XZ

Γ

XZ

Γ

XZ

Γ

X

Figure 6: Band structure for Na- and K-related defects in CZTSe. + + − Na− Se , KSe , Natc and Ktc defects. The band structures show that NaSe and KSe together with

NaSn and KSn are responsible for introducing a defect level within the band gap. However, we have to mention that due to the high formation energy of such defects, we expect that Na and K impurities in Se and Sn substitutional positions play a negligible role in the efficiency of solar cells based on CZTSe. The charge transition levels calculated for the large 216-atom supercell reveal that only Sn-related defects have transition level located close to the middle of the band gap. Therefore, they could act as a trap for carriers in CZTSe and become detrimental for the device performance. NaSn exhibit three transition levels: the ǫ(0/-) level around EVBM +0.17 eV, the (-/-2) 12

ACS Paragon Plus Environment

Page 13 of 20

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

level close to EVBM +0.55 eV that is a deep level, and the (-2/-3) level near EVBM +0.87 eV. KSn also shows three transition levels within the band gap: (0/-), (-/-2) and (-2/-3) which are located at EVBM +0.39 eV, EVBM +0.45 eV and EVBM +0.9 eV, respectively. Other Na- and K-related defects do not show any charge transition levels and therefore they do not act as carrier traps.

SUMMARY AND CONCLUSIONS We have carried out hybrid functional computations on native point defects as well as Na and K extrinsic defects in several substitutional and interstitial positions in CZTSe. According to our results the most favourable way to incorporate Na and K impurities is to form substitutional defects in the copper sublattice, similar to the CuInSe2 . The incorporation of Na in a tetrahedrally coordinated interstitial site is also highly probable, whereas the formation of K interstitials is not very likely, which can be understood due to the large size of the K impurity. Comparing the formation energies of two Na in two separate Cu sites with a (Na-Na) dumbbell in a Cu site reveals that Na atoms prefer to assemble as dumbbells in this site. The formation energy of (Na-Na) dumbbell in a Cu site is -3.18 eV, so such dumbbells are highly stable. Similarly, K atoms prefer to assemble as (K-K) dumbbells in Cu sites with a formation energy of -0.66 eV. If CZTSe is exposed to both Na and K simultaneously, then under p-type conditions (Na-K) mixed dumbbells can be formed as well. As a general remark, the formation energy of Na-related defects is always lower than their respective K-related counterparts in all studied substitutional and interstitial sites. Therefore, it is unlikely that K atoms replace Na from those sites that are already occupied by Na. Also, one has to keep in mind that most of the time Na diffuses out from the sodalime glass into CZTSe and K is incorporated into CZTSe via the PDT after the growth of the absorber. From the thermodynamic charge transition level calculations we can conclude that all Sn

13

ACS Paragon Plus Environment

The Journal of Physical Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 14 of 20

and Se vacancies create charge transition levels. While The gap state induced by Se vacancies is harmless for the carrier transport, the gap states created by Sn vacancies might be detrimental. The NaSn and KSn defects introduce transition levels within the band gap that might be detrimental for the device performance. Other extrinsic defects do not introduce gap levels.

Acknowledgment Financial support is gratefully acknowledged from the German Bundesministerium für Wirtschaft und Energie (BMWi) for the comCIGS II project (0325448C). E.Gh. would like to acknowledge IBM for the work proposal and for the provided computing resources.

References (1) Ahn, S.; Jung, S.; Gwak, J.; Cho, A.; Shin, K.; Yoon, K.; Park, D.; Cheong, H.; Yun, J. H. Determination of Band Gap Energy (Eg ) of Cu2 ZnSnSe4 Thin Films: On the Discrepancies of Reported Band Gap Values. Appl. Phys. Lett., 2010, 97, 021905 (2) Andersson, B. A. Materials Availability for Large-Scale Thin-Film Photovoltaics. Prog. in Photovoltaics: Research and Applications, 2000, 8, 61–76 (3) Babu, G. S.; Kumar, Y. K.; Bhaskar, P. U.; Vanjari, S. R. Effect of Cu/(Zn+Sn) Ratio on the Properties of Co-Evaporated Cu2 ZnSnSe4 Thin Films. Sol. Energy Mat. and Sol. Cells, 2010, 94, 221–226 (4) Blöchl, P. E. Projector Augmented-Wave Method. Phys. Rev. B, 1994, 50, 17953–17979 (5) Chen, S.; Gong, X. G.; Walsh, A.; Wei, S.-H. Crystal and Electronic Band Structure of Cu2 ZnSnX4 (X=S and Se) Photovoltaic Absorbers: First-Principles Insights. Appl. Phys. Lett., 2009, 94, 041903 14

ACS Paragon Plus Environment

Page 15 of 20

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

(6) Chen, S.; Gong, X.-G.; Walsh, A.; Wei, S.-H. Structural, Electronic and Defect Properties of Cu2 ZnSn(S,Se)4 Alloys. Computational Semiconductor Materials Science Symposium, 2011, 1370 (7) Chirilă, A.; Reinhard, P.; Pianezzi, F.; Bloesch, P.; Uhl, A.; Fella, C.; Kranz, L.; Keller, D.; Gretener, C.; Hagendorfer, H.; Jaeger, D.; Erni, R.; Nishiwaki, S.; Buecheler, S.; Tiwari, A. N. Potassium-Induced Surface Modification of Cu(In, Ga)Se2 Thin Films for High-Efficiency Solar Cells. Nat. Mat., 2013, 12, 1107–1111 (8) Choubrac, L.; Lafond, A.; Paris, M.; Guillot-Deudon, C.; Jobic, S. The Stability Domain of the Selenide Kesterite Photovoltaic Materials and NMR Investigation of the Cu/Zn Disorder in Cu2 ZnSnSe4 (CZTSe). Phys. Chem. Chem. Phys., 2015, 17, 15088–15092 (9) Contreras, M. A.; Romero, M. J.; To, B.; Hasoon, F.; Noufi, R.; Ward, S.; Ramanathan, K. Optimization of CBD CdS Process in High-Efficiency Cu(In, Ga)Se2 -Based Solar Cells. Thin Solid Films, 2002, 403-403, 204 – 211 (10) Friedlmeier, T. M.; Wieser, N.; Walter, T.; Dittrich, H.; Schock, H. W. Heterojuncitons Based on Cu2 ZnSnS4 and Cu2 ZnSnSe4 Thin Films. 14th European Conference of Photovoltaic Science and Engineering and Exhibition, 1997, 1242 (11) Gao, F.; Yamazoe, S.; Maeda, T.; Wada, T.

Structural Study of Cu-Deficient

Cu2(1−x) ZnSnSe4 Solar Cell Materials by X-Ray Diffraction and X-Ray Absorption Fine Structure. Jpn. J. Appl. Phys., 2012, 51, 10NC28 (12) Ghorbani, E.; Kiss, J.; Mirhosseini, H.; Roma, G.; Schmidt, M.; Windeln, J.; Kühne, T.; Felser, C. Hybrid-Functional Calculations on the Incorporation of Na and K Impurities Into the CuInSe2 and CuIn5 Se8 Solar-Cell Materials. J. of Phys. Chem. C, 2015, 119, 25197-25203

15

ACS Paragon Plus Environment

The Journal of Physical Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 16 of 20

(13) Green, M. A.; Emery, K.; Hishikawa, Y.; Warta, W.; Dunlop, E. D. Solar Cell Efficiency Tables (Version 42). Progress in Photovoltaics: Research and Applications, 2013, 21, 827–837 (14) Grossberg, M.; Krustok, J.; Timmo, K.; Altosaar, M. Radiative Recombination in Cu2 ZnSnSe4 Monograins Studied by Photoluminescence Spectroscopy.

Thin Solid

Films, 2009, 517, 2489 – 2492 (15) Han, D.; Sun, Y.; Bang, J.; Zhang, Y.; Sun, H.-B.; Li, X.-B.; Zhang, S.; Deep Electron Traps and Origin of p-Type Conductivity in the Earth-Abundant Solar-Cell Material Cu2 ZnSnS4 . Phys. Rev. B, 2013, 87, 155206 (16) Heyd, J.; Scuseria, G. E.; Ernzerhof, M. Hybrid Functionals Based on a Screened Coulomb Potential. J. Chem. Phys., 2003, 118, 8207–8215 (17) Hoffert, M. I.; Caldeira, K.; Jain, A. K.; Haites, E. F.; Harvey, L. D. D.; Potter, S. D.; Schlesinger, M. E.; Schneider, S. H.; Watts, R. G.; Wigley, T. M. L.; Wuebbles, D. J. Energy Implications of Future Stabilization of Atmospheric Co2 content. Nature 1998, 395, 881–884 (18) Ivády, V.; Abrikosov, I. A.; Janzén, E.; Gali, A. Role of Screening in the Density Functional Applied to Transition-Metal Defects in Semiconductors. Phys. Rev. B, 2013, 87, 205201 (19) Jackson, P.; Hariskos, D.; Wuerz, R.; Kiowski, O.; Bauer, A.; Friedlmeier, T. M.; Powalla, M. Properties of Cu(In, Ga)Se2 Solar Cells With New Record Efficiencies up to 21.7%. phys. status solidi (RRL), 2015, 9, 28–31. (20) Katagiri, H.; Saitoh, K.; Washio, T.; Shinohara, H.; Kurumadani, T.; Miyajima, S. Development of Thin Film Solar Cell Based on Cu2 ZnSnS4 Thin Films. Sol. Energy Mat. and Sol. Cells, 2001, 65, 141 –148

16

ACS Paragon Plus Environment

Page 17 of 20

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

(21) Komsa, H.-P.; Rantala, T. T.; Pasquarello, A. Finite-Size Supercell Correction Schemes for Charged Defect Calculations. Phys. Rev. B, 2012, 86 ,045112 (22) Kresse, G.; Furthmueller, J. Efficiency of Ab-Initio Total Energy Calculations for Metals and Semiconductors Using a Plane-Wave Basis Set. Comp. Mat. Sci., 1996, 6, 15–50 (23) Kresse, G.; Furthmueller, J. Efficient Iterative Schemes For Ab Initio Total-Energy Calculations Using a Plane-Wave Basis Set. Phys. Rev. B, 1996, 54, 11169–11186 (24) Kresse, G.; Joubert, D. From Ultrasoft Pseudopotentials to the Projector AugmentedWave Method. Phys. Rev. B, 1999, 59, 1758–1775 (25) Lany, S.; Zunger, A. Assessment of Correction Methods For the Band-Gap Problem and for Finite-Size Effects in Supercell Defect Calculations: Case Studies for ZnO and GaAs. Phys. Rev. B, 2008, 78, 235104 (26) Liu, M.-L.; Huang, F.-Q.; Chen, L.-D.; Chen, I.-W. A Wide-Band-Gap p-Type Thermoelectric Material Based on Quaternary Chalcogenides of Cu2 ZnSnQ4 (Q=S, Se). Appl. Phys. Lett., 2009, 94, 202103 (27) Maeda, T.; Kawabata, A.l; Wada, T. First-Principles Study on Alkali-Metal Effect of Li, Na, and K in Cu2 ZnSnS4 and Cu2 ZnSnSe4 . phys. status solidi (c) ,2015, 12, 631–637 (28) Matsushita, H.; Maeda, T.; Katsui, A.; Takizawa, T. Thermal Analysis and Synthesis From the Melts of Cu-Based Quaternary Compounds Cu–iii–iv–vi4 and Cu2 –ii–iv–vi4 (ii=Zn, Cd; iii=Ga, In; iv=Ge, Sn; vi=Se). J. Crys. Growth, 2000, 208, 416–422 (29) Oikkonen, L.; Ganchenkova, M.; Seitsonen, A.; Nieminen, R. Effect of Sodium Incorporation Into CuInSe2 From First Principles. J. Appl. Phys., 2013, 114, 083503 (30) Perdew, J. P.; Burke, K.; Ernzerhof, M. Generalized Gradient Approximation Made Simple. Phys. Rev. Lett., 1996, 77, 3865–3868 17

ACS Paragon Plus Environment

The Journal of Physical Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

(31) Pianezzi, F.; Reinhard, P.; Chirilă, A.; Bissig, B.; Nishiwaki, S.; Buecheler, S.; Tiwari, A. Unveiling the Effects of Post-Deposition Treatment with Different Alkaline Elements on the Electronic Properties of CIGS Thin Film Solar Cells. Phys. Chem. Chem. Phys., 2014, 16, 8843–8851 (32) Rau, U.; Schock, H. Electronic Properties of Cu(In, Ga)Se2 Heterojunction Solar Cells– Recent Achievements, Current Understanding, and Future Challenges. Appl. Phys. A, 1999, 69, 131–147 (33) Reinhard, P.; Chirilă, A.; Blösch, P.; Pianezzi, F.; Nishiwaki, S.; Buechelers, S.; Tiwari, A. Review of Progress Toward 20% Efficiency Flexible CIGS Solar Cells and Manufacturing Issues of Solar Modules. IEEE J. of Photovoltaics, 2013, 3, 572–580 (34) Reinhard, P.; Pianezzi, F.; Bissig, B.; Chirilă, A.; Blösch, P.; Nishiwaki, S.; Buecheler, S.; Tiwari, A. Cu(In, Ga)Se2 Thin-Film Solar Cells and Modules- A Boost in Efficiency Due to Potassium. IEEE J. of Photovoltaics, 2015, 5, 656–663 (35) Steirer, K. X.; Garris, R. L.; Li, J. V.; Dzara, M. J.; Ndione, P. F.; Ramanathan, K.; Repins, I.; Teetera, G.; Perkins, C. L. Co-Solvent Enhanced Zinc Oxysulfide Buffer Layers in Kesterite Copper Zinc Tin Selenide Solar Cells. Phys. Chem. Chem. Phys., 2015, 17, 15355–15364 (36) Van de Walle, C. G.; Denteneer, P. J. H.; Bar-Yam, Y.; Pantelides, S. T. Theory of Hydrogen Diffusion and Reactions in Crystalline Silicon. Phys. Rev. B, 1989, 39, 10791–10808 (37) Wang, W.; Winkler, M. T.; Gunawan, O.; Gokmen, T.; Todorov, T. K.; Zhu, Y.; Mitzi, D. B. Device Characteristics of CZTSSe Thin-Film Solar Cells with 12.6% Efficiency. Adv. Energy Mat., 2014, 4, 1301465 (38) Wei, S.-H.; Zhang, S. B.; Zunger, A. Effects of Na on the Electrical and Structural Properties of CuInSe2 . J. Appl. Phys., 1999, 85, 7214–7218 18

ACS Paragon Plus Environment

Page 18 of 20

Page 19 of 20

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

(39) Wuerz, R.; Eicke, A.; Kessler, F.; Paetel, S.; Efimenko, S.; Schlegel, C. CIGS ThinFilm Solar Cells and Modules on Enamelled Steel Substrates. Sol. Energy Mat. and Sol. Cells, 2012, 100, 132–137

19

ACS Paragon Plus Environment

The Journal of Physical Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

20

ACS Paragon Plus Environment

Page 20 of 20