Theoretical Investigation of Interaction of CuInSe2 Absorber Material

6 days ago - We performed ab-initio calculations to study oxygen and hydrogen point defects in CuInSe2 (CISe) solar–cell material. We found that H ...
7 downloads 0 Views 1MB Size
Subscriber access provided by University of South Dakota

C: Energy Conversion and Storage; Energy and Charge Transport 2

Theoretical Investigation of Interaction of CuInSe Absorber Material with Oxygen, Hydrogen, and Water Sudhir Kumar Sahoo, Ramya Kormath Madam Raghupathy, Thomas D. Kuehne, and Hossein Mirhosseini J. Phys. Chem. C, Just Accepted Manuscript • DOI: 10.1021/acs.jpcc.8b06709 • Publication Date (Web): 23 Aug 2018 Downloaded from http://pubs.acs.org on August 24, 2018

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 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 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.

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 29 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

Theoretical Investigation of Interaction of CuInSe2 Absorber Material with Oxygen, Hydrogen, and Water Sudhir K. Sahoo, Ramya Kormath Madam Raghupathy, Thomas D. K¨uhne, and Hossein Mirhosseini∗ Dynamics of Condensed Matter and Center for Sustainable Systems Design, Department of Chemistry, University of Paderborn, Warburger Str. 100, D–33098 Paderborn, Germany E-mail: [email protected] Abstract We performed ab–initio calculations to study oxygen and hydrogen point defects in CuInSe2 (CISe) solar–cell material. We found that H interstitial defects (when one H atom is surrounded by four Se atoms) and HCu (when a H atom is replacing a Cu atom) are the most stable defects. While these H substitutional defects remain neutral, H interstitial defects act as a donor defect and are detrimental for the cell performance. The incorporation of H2 into the CISe lattice, on the other hand, is harmless for the p–type conductivity. Oxygen atoms tend to either substitute Se atoms in the CISe lattice or form interstitial defects, though, the formation of substitutional defects are more favorable. All oxygen point defects have high formation energies, which results in the low concentration of these defects in CISe. However, the presence of oxygen in the system leads to the formation of secondary phases such as In2 O3 and InCuO2 . In addition to the point defects, we studied the adsorption of H2 O molecules on a defect–free surface and a surface with a (2VCu +InCu ) defect using the ab–initio

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

thermodynamics technique. Our results indicate that the dissociative water adsorption on the CISe surface is energetically unfavorable. Furthermore, in order to obtain a water–free surface, the surface with defects has to be calcined at a higher temperature compared to the defect–free surface.

Introduction Cu(In, Ga)Se2 (CIGSe) is considered to be the most promising absorber material for thin– film solar cells. 1–3 It is known that the properties of CIGSe films are very sensitive to the nature of defects, growing conditions, and heat treatment. 4–9 Point defects in CIGSe have attracted numerous interests in the scientific community in the hope of understanding the defects nature and improving the efficiency of CIGSe–based devices. Existing studies, so far, focused mainly on the understanding of the role of alkali–metal defects, which are found to alter the electronic properties of CIGSe films significantly. 10–21 A recent study shows that incorporating Rb into the CIGSe absorber improves the conversion efficiency of the cell up to 22.6 %, the world–record efficiency for thin–film solar cells. 22 Alkali atoms, however, are not the only impurities in CIGSe. The CIGSe thin films synthesized by coevaporation techniques in vacuum can be in contact with air after alkali– metal post deposition treatment when CIGSe films are rinsed with water to remove the excess alkali metals and secondary phases from the surface. 20,23–26 Oxygen and hydrogen atoms are also found in the CIGSe lattice after air annealing of the cells that improves the power–conversion efficiency. 7,8,27 Furthermore, oxygen atoms can diffuse from the back contact or the buffer layer (Zn(O, OH) layer for example) into CIGSe. 28 Water vapor can be deliberately introduced into the growth chamber during the growth of CIGS films to increase the hole carrier density and conductivity. 29 The coevaporation technique, usually employed for synthesizing effective CIGSe devices, requires high vacuum, which increases the cost for the thin–film synthesis. The solovothermal technique is one of the alternative non–vacuum methods to synthesize CIGSe. In this method, the annealing of the CIGSe thin 2

ACS Paragon Plus Environment

Page 2 of 29

Page 3 of 29 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

film is carried out in the presence of H2 gas to remove Cu2−x Se2 . 30 It is reported that excess of H2 gas causes etching of Se from the CIGSe surface that in turn can change the surface composition. 31 Oxygen- and hydrogen-based defects, despite their presence in the absorber layer 32,33 and their potential impacts on the cell performance, have not been characterized thoroughly. In this respect, it is necessary to understand the role of O- and H-point defects in CIGSe to prevent the formation of harmful defects which can hinder the further improvement of the cells efficiency. Absorption of water molecules on the CIGSe surface, on the other hand, could affect the post treatment of CIGSe films, CdS layer deposition for example, and should be considered for the optimization of the fabrication process. Hence, in this work we study H and O point defects in the CISe bulk and their effects on the electronic structure of CISe by performing density functional theory (DFT) based calculations. 34 Further, we probed the interactions of H2 O molecules with the CISe surface by employing ab–initio thermodynamics techniques. CIGSe is an alloy of ternary copper indium selenide (CISe) and copper gallium selenide (CGSe). It has been reported that the device efficiency of CIGSe–based solar cells is highest when the Ga/(In+Ga) ratio is nearly 0.30 and the band gap of the system is about 1.1–1.2 eV, which is close to that of CISe. 22,35–37 The atomic and electronic structures of CISe and CGSe are very similar. 38 It is expected that the study of a rather simpler system (CISe) can shed light on the properties of a more complicated system (CIGSe). Therefore, in this work we studied the CISe system. The paper is organized as follows: the computational details and the models employed in this work will be discussed in the next section followed by the results and discussion. A summary of salient results concludes this paper.

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

Page 4 of 29

Computational Details H and O Point Defects in CISe In order to model point defects in solids, total energy calculations of supercells with corresponding defects is well suited. 10,39 In these calculations, the size of the supercell is particularly important: it has to be large enough such that spurious interactions between a defect and its periodic images are minimal. 11,40 The bulk CISe was modelled by a supercell (Cu54 In54 Se108 ) of 216 atoms where the defect is created by either substituting or adding/removing an atom. The defect formation energy is computed as ∆Ef =Etot (D, q) − Etot (bulk) ±

X

ni µ0i + ∆µi

i

 (1)

  q + q EVBM + µe + ∆ν0/b + Ecorr , where Etot (D, q) and Etot (bulk) are the potential energies of the supercells with the defect (D) in the charge state q and the bulk supercell (without defect), respectively. In addition, ni is the number of atoms of species i that are either removed from or added to the system, whereas µ0i is the chemical potential of species i in the native elemental state, and ∆µi is the thermodynamic limit of the chemical potential in order to avoid the formation of secondary phases. The thermodynamic limits to the chemical potentials were computed by determining the stability region for all the competing structures with respect to the reference structure (see Figure 1 and Table 1 in the Supporting Information (SI)). The valence band maximum of the bulk is denoted as EVBM , µe is the Fermi energy position (electron chemical potential), which varies from 0 to the band gap of the system (1.04 eV for CISe), and ∆ν0/b represents the correction term for the electrostatic alignment that is calculated as the difference between the averaged electrostatic potential of the bulk and the host with a defect in the neutral q charge state. Moreover, Ecorr represents the correction term for the total energy when the

4

ACS Paragon Plus Environment

Page 5 of 29 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

system is charged which also depends on the size of the supercell. 41 All the calculations were performed using the projector augmented plane wave method 42 as implemented in the Vienna Ab initio Simulation Package (VASP). 43,44 A plane wave cutoff of 350 eV was used and all structures were optimized using the conjugate gradient method until the ionic forces were less than 0.01 eV/˚ A. The hybrid exchange correlation functional HSE06 45 was employed to calculate the electronic structures and defect formation energies. 46 In this work, the parameter α, which indicates the percentage of Hartree–Fock exchange to be included, is modified to 0.30 (mod–HSE06 hereafter). It is noted that the mod–HSE06 functional (α = 0.30) predicts the band gap of CISe (1.04 eV) very close to the experimental value. 11 To model the oxygen diffusion in the CISe lattice, we employed the climbing–image nudged elastic band (CI–NEB) method. 47 In this work, five equidistant images were created for each NEB calculation by a linear interpolation technique. The force constant for the string, which connects the images, was set to 5.0 eV/˚ A2 . It is computationally very expensive to carry out such calculations using hybrid functionals, which is why Hubbard–corrected DFT calculations with Hubbard on–site interaction parameter U = 5.0 eV for Cu 3d electrons has been employed here. 48 It is well known that PBE functional 49 underestimates the band gap of solids. By employing PBE+U functional the band gap and bulk structure of CISe improved (see Table 8 in the SI). It is to be noted that our results of the atomic structures, vacancy formation energies, as well as the barriers for diffusion of impurities in CISe are in a reasonable agreement with those calculated using the HSE06 hybrid functional. 50

H2 O on the CISe Surface The CISe bulk structure was initially optimized using a 1×1×1 supercell (Cu4 In4 Se8 ) where both the atoms and the lattice parameters were relaxed. The (110) surface was created from the optimized structure using the atomic simulation environment (ASE) program. 51 The surface is then translated along the X and Y directions, and a total of eight atomic layers (Cu64 In64 Se128 ) has been constructed. We employed the slab model for the surface where a 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 29

vertical vacuum of 14 ˚ A was created along the Z–direction. Hubbard–corrected DFT (PBE+U ) calculations with U = 5.0 eV for 3d electrons of Cu were carried out for all surfaces. The integration of the Brillouin zone was approximated by Γ–point only. All other computational details are the same as mentioned in the previous section. We employed the ab–initio thermodynamics technique 52–54 to investigate the structure of various H2 O/CISe interfaces at finite temperature and pressure. We first computed the surface free energy (γ) in order to find the water adsorption on the surface as

nH2 O(g) + CISe(110) → nH2 O/CISe(110),

(2)

where n, CISe(110), and nH2 O/CISe(110) are the number of H2 O molecules, clean surface (no water molecules on the surface), and surface with water, respectively. Specifically, γ is computed as

γ (T, p) =

1 [G (nH2 O/CISe(110)) − G (CISe(110)) − nH2 O µH2 O ] , A

(3)

where A is the area of the surface, G (nH2 O/CISe(110)) and G (CISe(110)) are the corresponding free energies, and µH2 O is the chemical potential of H2 O. The free energy of water absorption on the surface can be approximated to the potential energy difference of CISe(110) and nH2 O/CISe(110), 54,55 i.e.

γ (T, p) =

1 [E (nH2 O/CISe(110)) − E (CISe(110)) − nH2 O µH2 O ] , A

(4)

where E (nH2 O/CISe(110)) and E (CISe(110)) are the corresponding potential energies. In

6

ACS Paragon Plus Environment

Page 7 of 29 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

Eq. 4, only µH2 O is dependent on T and p, which is computed as µH2 O (T, p) =E (H2 O) + ∆µ (T, p) =E (H2 O) + ∆H

T, p0H2 O



− T ∆S

T, p0H2 O



 + kB T ln

pH 2 O p0H2 O

 ,

where E (H2 O) is the potential energy of a single H2 O molecule. Also, ∆H and ∆S are enthalpy and entropy differences which are taken from thermochemical table, 56 where p0H2 O , pH2 O , kB , and T are the standard pressure (1 bar), partial pressure of H2 O, Boltzmann constant, and temperature, respectively.

Results and Discussion H and O Point Defects in CISe The formation energies ∆Ef of various defects in their most stable state computed for point A in the CISe phase diagram (shown in Figure 1 in the SI) are given in Table 1. Point A in the CISe phase diagram 57 corresponds to the growing condition of the high–performance cells, which is Cu–poor and Se–rich. We note in passing that the chemical potentials of Cu, In, and Se for point A (Table 3 of the SI) are in an agreement with those calculated before. 46,58 The formation energies for all charge states are listed in Table 4 of the SI. For calculating the formation energies of O defects, the elemental chemical potential of O is modified to avoid the formation of the competing phases (see Table 2 of the SI). In the case of H defects, the chemical potential of H is taken from the H2 molecule. Our results show that among H substitutional point defects, HCu has the lowest ∆Ef . This is similar to other impurities such as Na and K, where NaCu and KCu substitutional defects were the most favourable defects in Cu–poor CISe. 11,59 The H atom in the HCu position forms only a covalent bond with the nearest Se atom owing to its small size (see Figure 2 in the SI). Similar to HCu , H in the HIn position forms

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

Table 1: The formation energies ∆Ef for various O- and H-point defects (in Cu–poor condition) in CISe in their most stable state. The two interstitial sites between four anions (Se atoms) and four cations (2Cu+2In) are denoted as int1 and int2, respectively. Defects HCu HIn HSe Hint1 Hint2 (H2 )Cu (H2 )int2 OCu OIn OSe Oint1 Oint2

q 0 0 1+ 1+ 1+ 10 0 0 0 0 0

∆Ef (eV) 0.34 2.08 2.02 0.29 0.80 2.08 0.85 3.34 4.29 1.74 1.89 2.85

a covalent bond with one of the Se atoms. But this defect has a larger ∆Ef compared to that of HCu . It is evident from Figure 1 that the HCu point defects remain charge neutral for almost all values of µe , i.e. they have no effect on carrier concentration. The H atom can also occupy the interstitial position and forms a bridging bond along the Cu–Se bond owing to its low value of ∆Ef . Moreover, the Hint1 defect is stable with 1+ charge (donor), which makes this defect detrimental for hole conductivity in CISe. HCu and Hint1 defects can form rather easily in CISe owing to their low formation energies. Therefore, one can expect that H atoms quench the VCu present in the top few layers of the CISe surface or form a H–Se bond at the surface. It is to be noted that our results are in a qualitatively agreement with those of Kilic and Zunger 59 that indicated H–related defects are either detrimental for the p–type conductivity or can neutralize the p typeness of CISe. We also considered the incorporation of H2 molecules into the CISe lattice. In particular, H2 in the Cu sublattice and H2 in the interstitial position (Figure 3 in the SI), where the H2 molecule remains intact. It is interesting to see that the nature of HCu and (H2 )Cu defects are very different from each other: (a) the defect formation energy is larger in the later

8

ACS Paragon Plus Environment

Page 8 of 29

Page 9 of 29

case and (b) (H2 )Cu enhances the hole concentration, while HCu is a neutral defect. It has been reported that VCu is stable in its 1− state 58,60 and this might be the reason for the enhanced p–type conductivity in Cu–poor CISe. Interestingly, when a H2 molecule occupies the Cu vacancy site, it remains negatively charged for all values of µe . This is to say that the occupation of Cu vacancies with H2 molecules does not decrease the hole concentration. In addition, when the H2 molecule occupies the interstitial position in the CISe lattice ((H2 )int2 ) as found in our calculations (see Table 1), it is stable in its charge neutral state. Based on our results, we can conclude that the incorporation of H2 molecules into the CISe lattice does not affect the p–type conductivity of CISe. In contrast, the presence of H atoms are detrimental either by the formation of the positively charged interstitial defects or by

2.5 2

∆Ef (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

The Journal of Physical Chemistry

1.5 1

0 HCu Hint1 OSe Oint1 (H2)Cu (H2)int2

0

1-

1-

0

0

1+

0.5

0 1-

0 0

0.2

0.4

0.6

0.8

1

µe (eV) Figure 1: Formation energies ∆Ef for various H and O point defects in CISe as a function of electron chemical potential (µe ). Point defects with low formation energies are shown here.

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

reducing the concentration of Cu vacancies. Regarding oxygen impurities, the most stable point defect forms when a Se atom is replaced by an O atom, as can be seen in Table 1. Similar to Se atoms in the CISe lattice, OSe forms four bonds with two Cu and two In atoms, though, the O–In bond length is smaller than that of Se–In (Figure 4 in the SI). Both O and Se belong to the same group in the periodic table, even though, O has a larger electronegativity. Therefore, OSe remains charge neutral and becomes negatively charged for µe larger than 0.85 eV. Considering interstitial defects, the O atoms occupy the bridging position along the In–Se bond. Oxygen interstitial defects remain charge neutral for the whole range of µe (see Figure 1). This means that Oint1 defects do not alter the carrier concentration in CISe. Considering the rather high values of ∆Ef of the oxygen defects, it can be concluded that these defects do not play a major role in the defect physics of CISe. The most striking effect of the presence of oxygen in the system is the formation of secondary phases: In2 O3 and InCuO2 are the most stable compounds that can coexist with CISe. 61 The band gap of In2 O3 is much larger than that of CISe 62 and the formation of this secondary phase could affect the efficiency of CIGSe–based devices. Various secondary phases that can form when O is present in the system are listed in Table 2 of the SI. Ishizuka et al. 29 have reported that the addition of water vapor during the thin–film deposition improves the crystal structure of CIGSe and the performance of the solar cell. In CuGaSe2 (CGSe), VSe defects act as donor defects 46 and are detrimental for the performance of the CGSe–based cell. The effects of VSe defects in CGSe can be compensated by the formation of OSe . 29 In CISe, however, the VSe defects are charge neutral. 46,60 Our results show that OSe is also charge neutral for p–type CISe. This means that in contrast to CGSe, the incorporation of O into the CISe lattice and formation of OSe cannot improve the p– conductivity. This is in agreement with experimental findings suggesting that the treatment of CIGSe with only O2 gas does not improve the cell performance. 29,63 We showed that hydrogen point defects are either detrimental or benign in CISe. Given these facts, we

10

ACS Paragon Plus Environment

Page 10 of 29

Page 11 of 29 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

might conclude that the cell performance improvement in the presence of water vapor in the fabrication chamber could be due to the enhanment of alkali metal effects. 29,63 In order to understand the diffusion of oxygen defects in the CISe lattice, we calculated the diffusion barriers for two specific paths: (a) diffusion of an O atom from one Se position to the nearest VSe site and (b) diffusion of an O atom from the most stable interstitial position to an equivalent interstitial position. We considered these two paths based on the corresponding low formation energies of these two point defects (see Table 1). The diffusion barriers obtained from CI–NEB calculations are given in Figures. 5 and 6 of the SI. The barrier for diffusion of an O atom from one Se lattice to the nearest VSe site is found to be 0.84 eV, whereas the barrier of migration of O atom from one interstitial site to another is found to be 1.15 eV. Considering the fact that ∆Ef and the diffusion barrier of OSe are smaller than those of Oint1 , we can conclude that OSe defects are more favorable than Oint1 and they can move in the CISe lattice by jumping from one VSe site to another.

H2 O on CISe Surface Surface Relaxation It has been reported that the (112) surface of CISe is the most exposed surface, however, this surface is a polar surface. 64,65 Among all non–polar surfaces, the (110) surface is the most stable one. 66 Interestingly, all atoms in the top layer of the (110) surface are tri–coordinated which is similar to the coordination number of atoms in the (112) surface. 66 Therefore, we believe that the chemical properties of these two surfaces are not much different. The atomic structure of the (110) surface before and after relaxation are shown in Figure 2. It can be seen that Cu and In atoms are moving downwards, while Se atoms are moving upwards during relaxation. The surface energy is calculated as

∆Esurf =

1 (Esurf − n × Ebulk ) , A 11

ACS Paragon Plus Environment

(5)

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

(a)

(b)

Figure 2: Structure of the CISe (110) surface before (a) and after (b) relaxation. Atoms color: Cu–cyan, In–pink and Se–yellow.

12

ACS Paragon Plus Environment

Page 12 of 29

Page 13 of 29 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

where A is the area of the surface and n is the number of CuInSe2 units present in the surface. Therein, Esurf and Ebulk are the potential energies of the relaxed surface and the bulk (of CuInSe2 unit), respectively. The values of ∆Esurf with respect to the number of layers that are relaxed indicate that mainly the top 3 layers of the (110) surface contribute to ∆Esurf (see Table 5 in SI). Therefore, in all following surface calculations only the top 4 layers are relaxed while all others are kept fixed. It is known that the topmost layers of the CISe surface are Cu poor. 67,68 It has been suggested that Cu–poor top layers of the CISe surface could be an ordered defect compound (ODC). 69,70 It has been also observed from theoretical studies that the introduction of (2VCu +InCu ) defect in CISe leads to the formation of ODC. 71 Hence, a (2VCu +InCu ) defect is created in the top layer of the (110) surface to represent the Cu–poor condition.

Single Water Molecule on the (110) Surface In order to find the structures of the H2 O/CISe interface, we first studied adsorption of a single water molecule on the CISe (110) surface. On this surface there are two types of unsaturated cations, Cu and In, where a water molecule could be adsorbed. It is to note that there is also a possibility of dissociative water adsorption on the surface, where a OH group and a H atom of the H2 O molecule can bind to one of the cations (either Cu or In) and Se atoms, respectively. The adsorption energy is calculated as

∆Eads = E (nH2 O/CISe(110)) − E (CISe(110)) − nH2 O E (H2 O) ,

(6)

where E (nH2 O/CISe(110)), E (CISe(110)), and E (H2 O) are the potential energies of nH2 O molecules adsorbed on the CISe (110) surface, the CISe (110) surface, and the H2 O molecule, respectively. The negative values of ∆Eads indicate that the adsorption of a H2 O molecule on the surface is thermodynamically favorable. Interestingly, from our calculations of ∆Eads , we found that the H2 O molecule is adsorbed

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

on the In atom instead of the Cu atom, as can be seen from Table 2 and Figure 7 of the SI. The In–O distance (2.39 ˚ A) is slightly larger than that of Cu–O (2.33 ˚ A). On the surface with the (2VCu +InCu ) defect, the water molecule is adsorbed on the InCu site (see Figure 3) and the corresponding ∆Eads value is found to be very large (−125.8 kJ mol−1 ). Moreover, ∆Eads of the H2 O molecule on the In atom close to the InCu site is larger than ∆Eads of the H2 O molecule on the In atom of the defect–free surface. This indicates that the presence of this defect changes the nature of water adsorption on the CISe surface. We also studied the dissociative water adsorption on both surfaces, defect–free and with defect. Interestingly, we found that ∆Eads is positive for both the surfaces (see Tables 2 and 3), which indicates that the dissociative water adsorption is energetically not . Recently, Senftle and Carter reported a study on band edge alignments at the CuInS2 (112) surface/water interface. 72 Interestingly, on the CuInS2 (112) surface the water molecule is adsorbed on In atoms rather than on Cu atoms. Moreover, dissociatively water adsorption is energetically not favorable. This indicates that our results on the CISe (110) surface are in the same line with the CuInS2 (112) surface. 72 Table 2: Adsorption energy ∆Eads for a water molecule on the defect–free (110) surface, where OIn and OCu represent the water molecule adsorbed on In and Cu atoms, respectively. Structure H2 OIn /CISe H2 OCu /CISe H2 Odiss In /CISe

∆Eads (kJ mol−1 ) −73.5 −32.5 17.3

Table 3: Adsorption energy ∆Eads for a water molecule on the (110) surface with (2VCu +InCu ) defect, where OInCu and OIn0 represent the water molecule adsorbed on InCu and the In atom close to InCu , respectively. Structure H2 OInCu /CISe H2 OIn0 /CISe H2 Odiss InCu /CISe

∆Eads (kJ mol−1 ) −125.8 −81.1 16.6

14

ACS Paragon Plus Environment

Page 14 of 29

Page 15 of 29 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

In

Cu

(a)

(b)

(c) Figure 3: Structure of H2 O molecule adsorbed on the InCu site (a), In atom close to InCu site (b) and dissociatively adsorbed over InCu site (c) on the CISe surface with a (2VCu +InCu ) defect. 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

Structure of H2 O on the (110) Surface In order to determine the structure of the H2 O/CISe interface, we optimized the structures where more than one water molecule was adsorbed on the surface. In the case of the defect– free surface, the second water molecule (when two water molecules are there) is also adsorbed on the In atom rather than the Cu atom. Similarly, other water molecules are adsorbed on In atoms only (Figure 8 in the SI). The computed adsorption energies per water molecule does not change much with respect to the number of H2 O molecules, which indicates that the water molecules on the surface are non interacting (Table 6 in the SI). In the case of a surface with defects, the second water molecule is adsorbed on the In atom close to InCu . The effect of temperature T and pressure pH2 O on the structure of CISe/H2 O interface is investigated using ab–initio thermodynamics, as explained in the previous section. We increased the number of water molecules one by one to the total number of In atoms on the surface (8 for the defect–free surface and 9 for the surface with defects). The computed free energies (γ) as a function of T (at pH2 O = 1 bar) are given in Figures 4 and 5 for the defect–free and the surface with defects, respectively. The negative and positive values of γ indicate adsorption and desorption of water on/from the surface, respectively. Our results show that at low temperatures, water molecules bind to the surface, while at higher temperatures water molecules tend to remain in the gas phase. We ascribe this trend to the fact that ∆Eads is more dominant at lower temperature, while at higher temperature the entropy of water is more dominant. In the case of the defect–free surface, the clean surface is obtained at 430 K and when the temperature is less than 418 K all In atoms are occupied with water molecules. For the surface with defects, the clean surface is obtained at a higher temperature (T > 682 K) and when the temperature is below 382 K all In atoms are occupied with water molecules. We find that this is due to the larger adsorption energy of the H2 O molecule on the surface with defects compared to the defect– free surface. It is evident from Figure 5 that there are many cross–overs in the diagram of the surface with a defect as compared to that of the defect–free surface. We ascribe this is 16

ACS Paragon Plus Environment

Page 16 of 29

Page 17 of 29

3 1H2O/CISe

2

-1

-2

γ (kJ mol Å )

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

1

2H2O/CISe 3H2O/CISe 4H2O/CISe 5H2O/CISe 6H2O/CISe 7H2O/CISe 8H2O/CISe

0 -1 -2 0

200

400

T (K)

600

800

1000

Figure 4: Surface free energy as a function of temperature for the defect–free CISe (110) surface.

17

ACS Paragon Plus Environment

The Journal of Physical Chemistry

3 1H2O/CISe

2

-1

-2

γ (kJ mol Å )

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 18 of 29

1 0

2H2O/CISe 3H2O/CISe 4H2O/CISe 5H2O/CISe 6H2O/CISe 7H2O/CISe 8H2O/CISe 9H2O/CISe

-1 -2 0

200

400

T (K)

600

800

1000

Figure 5: Surface free energy as a function of temperature for the CISe (110) surface with a (2VCu +InCu ) defect.

18

ACS Paragon Plus Environment

Page 19 of 29 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

due to different adsorption energies of water molecules on different In atoms in this surface. From our calculations, we could confirm that at ambient conditions (room temperature and 1 bar pressure) CISe surfaces are hydrated. Our results also suggest that in the presence of defects, it is required to calcine the CISe surface at a higher temperature in order to remove water molecules from the surface.

Conclusions We have studied the atomic and electronic structures of H- and O-point defects in CISe using DFT–based calculations. From the defect formation energies, it is evident that the H atom could form substitutional (HCu ) as well as interstitial (Hint1 ) defects. The incorporation of H atoms into CISe has detrimental effects on the performance of solar cell because (a) HCu could quench the excess hole generated by VCu , and (b) H interstitial defects act as a donor defect which could reduce the hole concentration. Interestingly, when H2 molecules are incorporated into the CISe lattice they have a benign effect on the carrier concentration, albeit, the formation energy is larger in comparison to H point defects. When oxygen is incorporated into CISe, the substitutional defect (OSe ) is more likely to form compared to the interstitial defect. All oxygen point defects remains charge neutral which means these defects do not alter the carrier concentration. The relatively low diffusion barriers of oxygen suggests that O atom could migrate from one Se site to another. We have also studied the structure of the H2 O/CISe interface using the ab–initio thermodynamics approach. We observe that water molecules tend to be adsorbed molecularly rather than dissociatively. On the defect–free surface, water molecules tend to absorb on top of the In atoms instead of the Cu atoms. In addition, on the surface with defects, water molecules tend to adsorb on the InCu site or on the In atoms next to these defects. The surface free energy calculations suggest that at higher temperature, the water molecule desorbs from the surface, while at lower temperature the surfaces are hydrated. Our results indicate

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

that to obtain a clean (water free) surface, the surface with defects has to be calcined at a higher temperature compared to the defect–free surface.

Acknowledgment The authors would like to acknowledge financial support from the German Bundesministerium f¨ ur Wirtschaft und Energie (BMWi) for the speedCIGS project (0324095C). The authors gratefully acknowledge the Gauss Centre for Supercomputing e.V. (www.gausscentre.eu) for funding this project by providing computing time on the GCS Supercomputer SuperMUC at the Leibniz Supercomputing Centre (www.lrz.de). The authors also acknowledge the Paderborn Center for Parallel Computing (PC2 ) (https://pc2.uni-paderborn.de/hpcservices/available-systems/oculus/) supercomputing time on OCuLUS where a part of the calculations were carried out.

Supporting Information The contents of the material supplied as Supporting Information: 1. Stability diagram for CISe 2. Calculated formation energies for different possible phases 3. Formation energies for various oxygen and hydrogen defects in CISe 4. Diffusion paths for O point defects in CISe 5. Structure of H2 O on the (110) surface 6. Comparison between the experimentally measured and theoretically calculated lattice parameters

20

ACS Paragon Plus Environment

Page 20 of 29

Page 21 of 29 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

References (1) Friedlmeier, T. M.; Jackson, P.; Bauer, A.; Hariskos, D.; Kiowski, O.; Menner, R.; Wuerz, R.; Powalla, M. High-Efficiency Cu(In,Ga)Se2 Solar Cells. Thin Solid Films 2017, 633, 13 – 17, E-MRS 2016 Spring Meeting, Symposium V, Thin-Film Chalcogenide Photovoltaic Materials. (2) Ramanujam, J.; Singh, U. P. Copper Indium Gallium Selenide Based Solar Cells - a Review. Energy Environ. Sci. 2017, 10, 1306–1319. (3) Chiril˘a, A.; Reinhard, P.; Pianezzi, F.; Bloesch, P.; Uhl, A. R.; Fella, C.; Kranz, L.; Keller, D.; Gretener, C.; Hagendorfer, H.; et al. Potassium-Induced Surface Modification of Cu(In,Ga)Se2 Thin Films for High-Efficiency Solar Cells. Nat. Mater. 2013, 12, 1107–1111. (4) Mccluskey, M. D.; Haller, E. E. Dopants and Defects in Semiconductors; CRS Press, Taylor & Francis Group: Boca Raton, 2012. (5) Kemell, M.; Ritala, M.; Leskel, M. Thin Film Deposition Methods for CuInSe2 Solar Cells. Crit. Rev. Solid State Mater. Sci. 2005, 30, 1–31. (6) Moons, E.; Gal, D.; Beier, J.; Hodes, G.; Cahen, D.; Kronik, L.; Burstein, L.; Mishori, B.; Shapira, Y.; Hariskos, D.; Schock, H.-W. Effect of Air Annealing on the Electronic Properties of CdS/Cu(In,Ga)Se2 Solar Cells. Sol. Energy Mater. Sol. Cells 1996, 43, 73–78. (7) Chung, Y.-D.; Cho, D.-H.; Park, N.-M.; Lee, K.-S.; Kim, J. Effect of Annealing on CdS/Cu(In,Ga)Se2 Thin-film Solar Cells. Curr. Appl. Phys. 2011, 11, S65–S67, International Conference on Electronic Materials. (8) Niu, X.;

Zhu, H.;

Liang, X.;

Guo, Y.;

21

Li, Z.;

ACS Paragon Plus Environment

Mai, Y. Air-Annealing of

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

Cu(In,Ga)Se2 /CdS and Performances of CIGS Solar Cells. Appl. Surf. Sci. 2017, 426, 1213–1220. (9) Garris, R. L.; Johnston, S.; Li, J. V.; Guthrey, H. L.; Ramanathan, K.; Mansfield, L. M. Electrical Characterization and Comparison of CIGS Solar Cells Made With Different Structures and Fabrication Techniques. Sol. Energy Mater. Sol. Cells 2018, 174, 77–83. (10) Malitckaya, M.; Komsa, H.-P.; Havu, V.; Puska, M. J. Effect of Alkali Metal Atom Doping on the CuInSe2 -Based Solar Cell Absorber. J. Phys. Chem. C 2017, 121, 15516–15528. (11) Ghorbani, E.; Kiss, J.; Mirhosseini, H.; Roma, G.; Schmidt, M.; Windeln, J.; K¨ uhne, T. D.; Felser, C. Hybrid-Functional Calculations on the Incorporation of Na and K Impurities into the CuInSe2 and CuIn5 Se8 Solar-Cell Materials. J. Phys. Chem. C 2015, 119, 25197–25203. (12) Hariskos, D.;

Powalla, M. Thermodynamic Limitations for Alkali Metals in

Cu(In,Ga)Se2 . J. Mater. Res. 2017, 32, 3789–3800. (13) Reinhard, P.; Bissig, B.; Pianezzi, F.; Avancini, E.; Hagendorfer, H.; Keller, D.; Fuchs, P.; D¨obeli, M.; Vigo, C.; Crivelli, P.; et al. Features of KF and NaF Postdeposition Treatments of Cu(In,Ga)Se2 Absorbers for High Efficiency Thin Film Solar Cells. Chem. Mater. 2015, 27, 5755–5764. (14) Hauschild, D.; Kreikemeyer-Lorenzo, D.; Jackson, P.; Friedlmeier, T. M.; Hariskos, D.; Reinert, F.; Powalla, M.; Heske, C.; Weinhardt, L. Impact of a RbF Postdeposition Treatment on the Electronic Structure of the CdS/Cu(In,Ga)Se2 Heterojunction in High-Efficiency Thin-Film Solar Cells. ACS Energy Lett. 2017, 2, 2383–2387. (15) Nicoara, N.; Kunze, T.; Jackson, P.; Hariskos, D.; Duarte, R. F.; Wilks, R. G.; Witte, W.; Br, M.; Sadewasser, S. Evidence for Chemical and Electronic Nonunifor-

22

ACS Paragon Plus Environment

Page 22 of 29

Page 23 of 29 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

mities in the Formation of the Interface of RbF-Treated Cu(In,Ga)Se2 with CdS. ACS Appl. Mater. Interfaces 2017, 9, 44173–44180. (16) Muzzillo, C. P. Review of Grain Interior, Grain Boundary, and Interface Effects of K in CIGS Solar Cells: Mechanisms for Performance Enhancement. Sol. Energy Mater. Sol. Cells 2017, 172, 18–24. (17) Avancini, E.; Carron, R.; Weiss, T. P.; Andres, C.; Brki, M.; Schreiner, C.; Figi, R.; Romanyuk, Y. E.; Buecheler, S.; Tiwari, A. N. Effects of Rubidium Fluoride and Potassium Fluoride Postdeposition Treatments on Cu(In,Ga)Se2 Thin Films and Solar Cell Performance. Chem. Mater. 2017, 29, 9695–9704. (18) Oikkonen, L. E.; Ganchenkova, M. G.; Seitsonen, A. P.; Nieminen, R. M. Effect of Sodium Incorporation into CuInSe2 from First Principles. J. Appl. Phys. 2013, 114, 083503. (19) Maeda, T.; Kawabata, A.; Wada, T. First-Principles Study on Alkali-Metal Effect of Li, Na, and K in CuInSe2 and CuGaSe2 . Jpn. J. Appl. Phys. 2015, 54, 08KC20. (20) Yuan, Z.-K.; Chen, S.; Xie, Y.; Park, J.-S.; Xiang, H.; Gong, X.-G.; Wei, S.-H. NaDiffusion Enhanced P-type Conductivity in Cu(In,Ga)Se2 : A New Mechanism for Efficient Doping in Semiconductors. Adv. Energy Mater. 2016, 6, 1601191–n/a, 1601191. (21) Reinhard, P.; Pianezzi, F.; Bissig, B.; Chiril˘a, A.; Blsch, P.; Nishiwaki, S.; Buecheler, S.; Tiwari, A. N. Cu(In,Ga)Se2 Thin-Film Solar Cells and Modules-A Boost in Efficiency Due to Potassium. IEEE J. Photovolt. 2015, 5, 656–663. (22) Jackson, P.; Wuerz, R.; Hariskos, D.; Lotter, E.; Witte, W.; Powalla, M. Effects of Heavy Alkali Elements in Cu(In,Ga)Se2 Solar Cells with Efficiencies up to 22.6%. Phys. Status Solidi RRL 2016, 10, 583–586.

23

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

(23) Forest, R. V.; Eser, E.; McCandless, B. E.; Chen, J. G.; Birkmire, R. W. Reversibility of (Ag,Cu)(In,Ga)Se2 Electrical Properties with the Addition and Removal of Na: Role of Grain Boundaries. J. Appl. Phys. 2015, 117, 115102. (24) Handick, E. et al. Formation of a KInSe Surface Species by NaF/KF Postdeposition Treatment of Cu(In,Ga)Se2 Thin-Film Solar Cell Absorbers. ACS Appl. Mater. Interfaces 2017, 9, 3581–3589, PMID: 28058843. (25) Forest, R. V.; McCandless, B. E.; He, X.; Rockett, A. A.; Eser, E.; Dobson, K. D.; Birkmire, R. W. Diffusion of Sodium in Single Crystal CuInSe2 . J. Appl. Phys. 2017, 121, 245102. (26) Khatri, I.; Shudo, K.; Matsuura, J.; Sugiyama, M.; Nakada, T. Comparative Study of Water and Ammonia Rinsing Processes of Potassium Fluoride-Treated Cu(In,Ga)Se2 Thin Film Solar Cells. Jpn. J. Appl. Phys. 2017, 56, 08MC12. (27) Joshi, R. A.; Taur, V. S.; Sharma, R. Effect of Annealing on Conversion Efficiency of Nanostructured CdS/CuInSe2 Heterojunction Thin Film Solar Cell Prepared by Chemical Ion Exchange Route at Room Temperature. Mater. Res. Bull. 2012, 47, 2206–2211. (28) Masanobu, I.; Shinya, S.; Tetsuhito, O.; Yuya, K.; Takeshi, M.; Hironori, K.; Hajime, S.; Shigeru, N. Structure of Chemically Deposited Zn(S,O,OH) Buffer Layer and the Effects on the Performance of Cu(In,Ga)Se2 Solar Cell. Prog. Photovolt: Res. Appl. 2015, 24, 397–404. (29) Ishizuka, S.; Shibata, H.; Yamada, A.; Fons, P.; Sakurai, K.; Matsubara, K.; Niki, S.; Yonemura, M.; Nakanishi, H. Control of the Thin Film Properties of Cu(In,Ga)Se2 Using Water Vapor Introduction During Growth. J. Appl. Phys. 2006, 100, 096106. (30) Wang, H.-C.; Wang, C.-C.; Feng, S.-W.; Chen, L.-H.; Lin, Y.-S. Synthesis of CIGS Thin Film by Solvothermal Route. Opt. Mater. Express 2013, 3, 54–66. 24

ACS Paragon Plus Environment

Page 24 of 29

Page 25 of 29 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

(31) Yakushev, M. V.; Ogorodnikov, I. I.; Volkov, V. A.; Mudryi, A. V. Incorporation of Hydrogen in CuInSe2 : Improvements of the Structure. J. Vac. Sci. Technol. A 2011, 29, 051201. (32) Noufi, R.; Powell, R.; Matson, R. On the Effect of Stoichiometry and Oxygen on the Properties of CuInSe2 Thin Films and Devices. Solar Cells 1987, 21, 55–63. (33) Niles, D. W.; Al-Jassim, M.; Ramanathan, K. Direct Observation of Na and O Impurities at Grain Surfaces of CuInSe2 Thin Films. J. Vac. Sci. Technol. A 1999, 17, 291–296. (34) Freysoldt, C.; Grabowski, B.; Hickel, T.; Neugebauer, J.; Kresse, G.; Janotti, A.; Van de Walle, C. G. First-Principles Calculations for Point Defects in Solids. Rev. Mod. Phys. 2014, 86, 253–305. (35) 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. (36) Dejene, F. The Structural and Material Properties of CuInSe2 and Cu(In,Ga)Se2 Prepared by Selenization of Stacks of Metal and Compound Precursors by Se Vapor for Solar Cell Applications. Sol. Energy Mater. Sol. Cells 2009, 93, 577 – 582. (37) Ishizuka, S.; Taguchi, N.; Nishinaga, J.; Kamikawa, Y.; Tanaka, S.; Shibata, H. Group III Elemental Composition Dependence of RbF Postdeposition Treatment Effects on Cu(In,Ga)Se2 Thin Films and Solar Cells. J. Phys. Chem. C 2018, 122, 3809–3817. (38) Belhadj, M.; Tadjer, A.; Abbar, B.; Bousahla, Z.; Bouhafs, B.; Aourag, H. Structural, Electronic and Optical Calculations of Cu(In,Ga)Se2 Ternary Chalcopyrites. Phys. Status Solidi (b) 2004, 241, 2516–2528.

25

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 26 of 29

(39) Nakamura, S.; Maeda, T.; Wada, T. First-Principles Study of Diffusion of Cu and In Atoms in CuInSe2 . Jpn. J. Appl. Phys. 2013, 52, 04CR01. (40) Komsa, H.-P.; Rantala, T. T.; Pasquarello, A. Finite-Size Supercell Correction Schemes for Charged Defect Calculations. Phys. Rev. B 2012, 86, 045112. (41) 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. (42) Bl¨ochl, P. E. Projector Augmented-Wave Method. Phys. Rev. B 1994, 50, 17953–17979. (43) Kresse, G.; Furthm¨ uller, 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. (44) Kresse, G.; Furthm¨ uller, J. Efficient Iterative Schemes for Ab Initio Total-Energy Calculations Using a Plane-Wave Basis Set. Phys. Rev. B 1996, 54, 11169–11186. (45) Heyd, J.; Scuseria, G. E.; Ernzerhof, M. Hybrid Functionals Based on a Screened Coulomb Potential. J. Chem. Phys. 2003, 118, 8207–8215. (46) Pohl, J.; Albe, K. Intrinsic Point Defects in CuInSe2 and CuGaSe2 as Seen via ScreenedExchange Hybrid Density Functional Theory. Phys. Rev. B 2013, 87, 245203. (47) Henkelman, G.; Uberuaga, B. P.; J´onsson, H. A Climbing Image Nudged Elastic Band Method for Finding Saddle Points and Minimum Energy Paths. J. Chem. Phys. 2000, 113, 9901–9904. (48) Dudarev, S. L.; Botton, G. A.; Savrasov, S. Y.; Humphreys, C. J.; Sutton, A. P. Electron-Energy-Loss Spectra and the Structural Stability of Nickel Oxide:

An

LSDA+U Study. Phys. Rev. B 1998, 57, 1505–1509. (49) Perdew, J. P.; Burke, K.; Ernzerhof, M. Generalized Gradient Approximation Made Simple. Phys. Rev. Lett. 1996, 77, 3865–3868. 26

ACS Paragon Plus Environment

Page 27 of 29 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

(50) Kiss, J.; Gruhn, T.; Roma, G.; Felser, C. Theoretical Study on the Diffusion Mechanism of Cd in the Cu-Poor Phase of CuInSe2 Solar Cell Material. J. Phys. Chem. C 2013, 117, 25933–25938. (51) Larsen, A. H. et al. The Atomic Simulation Environmenta Python Library for Working With Atoms. J. Phys.: Condens. Matter 2017, 29, 273002. (52) Reuter, K.; Scheffler, M. Composition, Structure, and Stability of RuO2 (110) as a Function of Oxygen Pressure. Phys. Rev. B 2001, 65, 035406. (53) Meyer, B. First-Principles Study of the Polar O-Terminated ZnO surface in Thermodynamic Equilibrium With Oxygen and Hydrogen. Phys. Rev. B 2004, 69, 045416. (54) Guhl, H.; Miller, W.; Reuter, K. Water Adsorption and Dissociation on SrTiO3 (001) Revisited: A Density Functional Theory Study. Phys. Rev. B 2010, 81, 155455. (55) Ghosh, T. K.; Nair, N. N. Nature of β-TaON Surfaces at Ambient Conditions. Surf. Sci. 2015, 635, 19–26. (56) Chase, M. W. J. Phys. Chem. Ref. Data. 1998, NIST-JANAF Thermochemical Tables. (57) Buckeridge, J. Chemical Potential Limits Analysis Program. https://github.com/ jbuckeridge/cplap. (58) Maria, M.; Hannu-Pekka, K.; Ville, H.; J., P. M. First-Principles Modeling of Point Defects and Complexes in Thin-Film Solar-Cell Absorber CuInSe2 . Adv. Electron. Mater. 2017, 3, 1600353. (59) Kilic, C.; Zunger, A. Doping of Chalcopyrites by Hydrogen. Appl. Phys. Lett. 2003, 83, 2007–2009. (60) Oikkonen, L. E.; Ganchenkova, M. G.; Seitsonen, A. P.; Nieminen, R. M. Vacancies in CuInSe 2 : New Insights From Hybrid-Functional Calculations. J. Phys.: Condens. Matter 2011, 23, 422202. 27

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

(61) Niki, S.; Kim, I.; Fons, P.; Shibata, H.; Yamada, A.; Oyanagi, H.; Kurafuji, T.; Chichibu, S.; Nakanishi, H. Effects of Annealing on CuInSe2 Films Grown by Molecular Beam Epitaxy. Sol. Energy Mater. Sol. Cells 1997, 49, 319–326. (62) King, P. D. C.; Veal, T. D.; Fuchs, F.; Wang, C. Y.; Payne, D. J.; Bourlange, A.; Zhang, H.; Bell, G. R.; Cimalla, V.; Ambacher, O.; Egdell, R. G.; Bechstedt, F.; McConville, C. F. Band Gap, Electronic Structure, and Surface Electron Accumulation of Cubic and Rhombohedral In2 O3 . Phys. Rev. B 2009, 79, 205211. (63) Ishizuka, S.; Sakurai, K.; Yamada, A.; Matsubara, K.; Shibata, H.; Yonemura, M.; Nakamura, S.; Nakanishi, H.; Kojima, T.; Niki, S. Water Vapor Introduction During Cu(In1−x Gax )Se2 Thin-Film Deposition and its Effect on Solar Cell Performance. 2006 IEEE 4th World Conference on Photovoltaic Energy Conference. 2006; pp 338–342. (64) Jaffe, J. E.; Zunger, A. Defect-Induced Nonpolar-to-Polar Transition at the Surface of Chalcopyrite Semiconductors. Phys. Rev. B 2001, 64, 241304. (65) Siebentritt, S.; Papathanasiou, N.; Albert, J.; Lux-Steiner, M. C. Stability of Surfaces in the Chalcopyrite System. Appl. Phys. Lett. 2006, 88, 151919. (66) Shigemi, A.; Wada, T. Surface Stabilities of Various Crystal Faces of CuInSe2 and Related Compounds by First-Principles Calculation. Jpn. J. Appl. Phys. 2012, 51, 10NC22. (67) Liao, D.; Rockett, A. Cu Depletion at the CuInSe2 Surface. Appl. Phys. Lett. 2003, 82, 2829–2831. ¨ ur, B.; Greiner, D.; Kauf(68) Br¨oker, S.; K¨ uck, D.; Timmer, A.; Lauermann, I.; Ums¨ mann, C. A.; M¨onig, H. Correlating the Local Defect-Level Density with the Macroscopic Composition and Energetics of Chalcopyrite Thin-Film Surfaces. ACS Appl. Mater. Interfaces 2015, 7, 13062–13072, PMID: 26010380.

28

ACS Paragon Plus Environment

Page 28 of 29

Page 29 of 29 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

(69) Schmid, D.; Ruckh, M.; Grunwald, F.; Schock, H. W. Chalcopyrite/Defect Chalcopyrite Heterojunctions on the Basis of CuInSe2 . J. Appl. Phys. 1993, 73, 2902–2909. (70) M¨onig, H.; Fischer, C.-H.; Caballero, R.; Kaufmann, C.; Allsop, N.; Gorgoi, M.; Klenk, R.; Schock, H.-W.; Lehmann, S.; Lux-Steiner, M.; Lauermann, I. Surface Cu Depletion of Cu(In,Ga)Se2 Films: An Investigation by Hard X-ray Photoelectron Spectroscopy. Acta Materialia 2009, 57, 3645 – 3651. (71) Zhang, S. B.; Wei, S.-H.; Zunger, A.; Katayama-Yoshida, H. Defect Physics of the CuInSe2 Chalcopyrite Semiconductor. Phys. Rev. B 1998, 57, 9642–9656. (72) Senftle, T. P.; Carter, E. A. Theoretical Determination of Band Edge Alignments at the Water-CuInS2 (112) Semiconductor Interface. Langmuir 2017, 33, 9479–9489.

Figure 6: TOC

29

ACS Paragon Plus Environment