A Strategy to Stabilize Kesterite CZTS for High-Performance Solar Cells

Apr 3, 2015 - Cu2ZnSnS4–xSex (CZTS) is an important semiconductor with significant potential for applications in the next generation of solar cells...
0 downloads 0 Views 531KB Size
Subscriber access provided by SETON HALL UNIV

Article

A Strategy to Stabilize Kesterite CZTS for High Performance Solar Cells Kuang Yu, and Emily A. Carter Chem. Mater., Just Accepted Manuscript • DOI: 10.1021/acs.chemmater.5b00172 • Publication Date (Web): 03 Apr 2015 Downloaded from http://pubs.acs.org on April 6, 2015

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.

Chemistry of Materials 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 24

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

Chemistry of Materials

A Strategy to Stabilize Kesterite CZTS for High Performance Solar Cells Kuang Yu1 and Emily A. Carter1,2,* 1

Department of Mechanical and Aerospace Engineering and 2Program in Applied and Computational

Mathematics and Andlinger Center for Energy and the Environment, Princeton University, Princeton, New Jersey 08544-5263

ABSTRACT Cu2ZnSnS4-xSex (CZTS) is an important semiconductor with significant potential for applications in the next generation of solar cells. CZTS has an optimal band gap (~1.5 eV) and contains no expensive or toxic elements. However, CZTS-based solar cells suffer from low efficiency due to poor crystal quality, which is partly caused by secondary phase formation during synthesis. We use density functional theory + U calculations to systematically investigate the stabilities of three CZTS phases: kesterite, stannite, and wurtzite. In agreement with previous experiment and theory, we confirm that these three phases have very similar formation energies. This finding is consistent with the known difficulties in synthesizing pure kesterite CZTS, the phase desirable for photovoltaic applications. To overcome this problem, we characterize surfaces and interfaces of CZTS and are able to identify certain “beneficial surfaces” that could be exploited to potentially provide extra stability for the kesterite phase.

* Author to whom correspondence should be addressed. Electronic mail: [email protected].

1

ACS Paragon Plus Environment

Chemistry of Materials

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 2 of 24

We propose the zinc blende ZnS (001) surface as a substrate to induce formation of these beneficial surfaces and to stabilize the kesterite phase, thereby serving as an effective crystallization template for the fabrication of high performance CZTS solar cells.

Interest in developing clean and sustainable sources of electricity continues to rise, with unprecedented levels of investments in renewable energy technologies. Among these technologies, photovoltaics (PVs) represent a particularly promising and fast growing industry, as the sun provides a near-inexhaustible supply of free, clean, non-CO2-producing energy. To date, multicrystalline silicon (Si) PVs, with ~20% efficiency at best,1 dominate the commercial market, mainly due to the abundance of Si and the mature fabrication techniques provided by the existing electronics industry. Bulk Si features an indirect band gap, which on the one hand prevents electron-hole recombination, while on the other hand reduces light absorption. Therefore, thick Si layers are necessary to enhance sunlight absorption, increasing materials and manufacturing energy costs. In order to overcome this disadvantage, thin-film technologies, which reduce the thickness of the absorption layer to ~micrometers, have been devised based on various semiconductors with direct band gaps, including GaAs, CdTe, CuInxGa1-xSe2 (CIGS), and more recently organic-inorganic hybrid perovskites. Within this category, a quaternary compound Cu2ZnSnS4-xSex (CZTS or CZTSSe) is of particular interest, as it has an almost perfect band gap (1.0~1.5 eV with different Se concentrations) and does not contain any expensive metals (such as Ga in GaAs and In in CIGS), or toxic elements (such as As in GaAs, Cd in CdTe, and Pb in the perovskites).2 Despite their significant economic and environmental advantages, CZTS-based PVs suffer from relatively low efficiencies. Compared to other thin film semiconductors, such as GaAs (~28.8%),3 CdTe 2

ACS Paragon Plus Environment

Page 3 of 24

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

Chemistry of Materials

(~20.4 %),1 CIGS (~23.3%),1 or the fast developing hybrid perovskite (~19.3%),4 the efficiency record for CZTS remains at 11% thus far,5 which is similar to organic solar cells3 but about half of the efficiency of state-of-the-art Si modules. Previous photoluminescence (PL) experiments6, 7 indicate that the minority carrier lifetime in CZTS is mainly limited by the nonradiative Shockley-Read-Hall (SRH) recombination process, which is a major bottleneck for performance improvement. In the SRH process, electron-hole pair recombination is mediated by deep mid-gap levels, which are typically associated with defects, impurities, or unsaturated dangling bonds residing at interfaces or phase boundaries. Thus, in order to further improve the performance of CZTS solar cells, it is crucial to develop better strategies to synthesize high quality CZTS single crystals, diminishing the concentration of recombination centers.

Figure 1. CZTS (pure sulfide) crystal structures for: (a) kesterite (space group I 4 ), (b) stannite (space group I 42m ), and (c) wurtzite (space group Pc ) phases. Formation of secondary phases during crystal growth prevents the synthesis of high quality CZTS crystals. These secondary phases are either secondary compounds, such as ZnS, Cu2S, and Cu2SnS3, or multiple phases of CZTS itself, namely kesterite, stannite, or wurtzite (see Figure 1). Similar to CIGS and CdS, the structures of the first two phases (kesterite and stannite) are both based on the cubic zinc blende topology, and the only difference lies in the stacking pattern of the Cu/Zn sub-lattice. The wurtzite phase of CZTS has a lower symmetry unit cell that resembles the wurtzite phase of ZnS. 3

ACS Paragon Plus Environment

Chemistry of Materials

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 24

Neutron diffraction experiments8 and various computational studies that use pure or hybrid density functionals9-11 identify the kesterite phase as most stable. The kesterite phase is also the most desirable structure since both stannite and wurtzite phases have ~100 meV lower band gaps compared to kesterite 12, 13

. Hence, the admixture of multiple phases would lower the final output voltage and be detrimental to

solar cell efficiency. Previous computations9-11 indicate that the energy difference between the kesterite and the stannite phases is less than 3.2 kJ/mol/formula-unit, so stannite is likely to form even though it is slightly less stable than kesterite. Partial disorder of the Cu/Zn sub-lattice was indeed observed in the neutron diffraction and electron energy loss spectroscopy studies.8,

14

This implies that the ordered

Cu/Zn sub-lattice is not stable, which is consistent with the small energy difference between the kesterite and the stannite phases. As part of efforts to improve CZTS crystal quality, new strategies are needed to further stabilize the kesterite phase and minimize secondary wurtzite, stannite or other disordered phases. In this work, we use the two ordered secondary phases (stannite and wurtzite) as a comparison to search for strategies to stabilize both the tetragonal cell and the Cu/Zn stacking pattern of the kesterite phase. Disordered phases are not explicitly considered here as they would require formidably large simulation cells to model realistically. It is known that in very thin films, surface properties play a much more pronounced role than in bulk samples. For example, surface energies play a critical role in the phase ordering of nanostructured zirconia: calculations show that for low-index cleavage planes, the tetragonal phase is energetically favored over the monoclinic phase because of significantly lower surface energies, despite a lower energy of monoclinic ZrO2 in the bulk.15 This finding is in accord with previous experimental work demonstrating that tetragonal rather than monoclinic ZrO2 forms in nanolaminates.16 Similarly, even though bulk FeO is thermodynamically unstable under ambient conditions,17 it has been observed in both nanofilms18 and nanometer-sized islands.19 Small differences in surface energies are also able to reverse the relative stabilities of the rutile and anatase phases of nanocrystalline TiO2.20 And a strong 4

ACS Paragon Plus Environment

Page 5 of 24

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

Chemistry of Materials

size-dependence was observed for wurtzite/rock-salt phase transitions in CdSe nanocrystals,21 confirming that this effect is not just limited to oxides. Inspired by these findings, we propose to use a similar strategy for CZTS to attempt to enhance the stability of kesterite. Most previous computational studies of CZTS used density functional theory (DFT) with generalized gradient approximation (GGA) exchange-correlation functionals,9,

22-25

such as Perdew-Wang 91

(PW91)26 or Perdew-Burke-Ernzerhof (PBE),27 which are fast to evaluate but potentially problematic for first-row transition-metal-containing compounds. The approximate exchange in these local functionals tends to poorly describe the interactions between tightly held 3d electrons in these materials. Some other calculations used hybrid functionals (such as Heyd-Scuseria-Ernzerhof 06, HSE0628),10, 29, 30 which are quite accurate for calculating energetics but are also orders of magnitude more expensive and therefore cannot be readily applied to phonon, extended surface, or interface calculations requiring large cells. In this work, we instead use the DFT+U formalism31, 32 with ab initio-derived U parameters,33, 34 which is more accurate than pure DFT-GGA while retaining high computation throughput. With this methodology, we compute the relative stabilities of all three CZTS phases, as well as other possible secondary phases. We calculate cleavage energies of different surfaces of CZTS and thereby identify “beneficial surfaces” that might be exploitable to stabilize the desirable kesterite structure. Finally, we offer evidence that the corresponding ZnS (zinc blende) surface has potential to act as a growth template to synthesize high-quality, single-crystal, kesterite CZTS films.

METHODS AND COMPUTATIONAL DETAILS The VASP 5.3.3 program35-38 was used to perform all-electron projector-augmented-wave (PAW)39 DFT+U calculations within the frozen-core approximation, in order to self-consistently optimize the valence/outer core electron distributions. PAW potentials acted on the 4s/3d electrons of Cu and Zn, the 5s/5p/4d electrons of Sn, the 4s/4p electrons of Se and the 3s/3p electrons of S, while all other (core) 5

ACS Paragon Plus Environment

Chemistry of Materials

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 24

electrons were kept frozen .40 For sulfur, the softer PAW potential with a larger augmentation region was used to accelerate basis set convergence (vide infra). In combination with PBE exchangecorrelation,27 we utilized the U correction developed by Dudarev et al.32 for the 3d electrons in Zn(II), Cu(I), and the 4d electrons in Sn(IV), described further below. We also found that it is critical to include long-range dispersion to obtain reasonable cell parameters for non-bonded layered structures like SnS2. Including dispersion also significantly improves bulk moduli predictions compared to experimental values (vide infra). We therefore applied Grimme’s a posteriori D2 correction41 to account for this interaction and denote our overall computational model then as PAW-PBE+U+D2. The bulk moduli were derived via fitting pressure-volume curves to the Murnaghan equation of state.42 A planewave basis with a 700 eV kinetic energy cutoff was used in all our PAW-PBE+U+D2 calculations, which is sufficient to converge the total energies of all systems examined to within 1 meV/atom. The Γ-point-centered Monkhorst-Pack scheme43 was used to sample the Brillouin zone, with the k-point meshes dense enough to converge all our energies to within 1 meV/atom. The k-point meshes used for each periodic cell model are given in the Supporting Information. The electron occupancy in reciprocal space was determined using the Gaussian smearing method, with a 0.05 eV smearing width. The effective U parameters for Zn(II) and Cu(I) were taken from reference 44 (4.5 eV and 3.6 eV, respectively), and were derived from ab initio unrestricted Hartree-Fock computations on electrostatically embedded metal oxide clusters.33 This method for determining the effective U (actually U-J, where U and J are respectively the effective intra-atomic Coulomb and exchange energies) requires no experimental data as input. The U and J parameters were derived fully ab initio by computing the average onsite Coulomb and exchange interactions between the UHF orbitals. The same method was used here to obtain the U parameter for Sn(IV). All embedded cluster models were constructed from the experimental SnO2 crystal structure,45 and we used the 6-31G(d) basis set for O,46 in conjunction with 6

ACS Paragon Plus Environment

Page 7 of 24

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

Chemistry of Materials

the large core Stuttgart-Koel effective core potential (ECP28MDF) and valence double zeta basis set47, 48 for Sn. All clusters were embedded in a point charge array with ~4000 point charge sites to reproduce the Madelung potential of the ionic crystal lattice, with the cluster atoms replacing point charges in the center of the array. The Sn(IV) sites directly bonded to the terminating oxygen ions in the clusters were represented not by point charges but rather by capping Ge4+ ECPs (LANL2DZ ECP49) to approximately account for short range repulsive interactions between the oxygen ions and the Sn(IV) core electrons. A Ge4+ ECP was used since Ge4+ has a similar valence structure and size as Sn (IV), hence reproducing similar interactions with the terminal anions. The U value was converged to within 0.5 eV with respect to cluster size; the resulting U value for Sn(IV) is 4.8 eV. Earlier work has shown that changing U by up to 0.5 eV has little effect on properties predicted by DFT+U.50 Using the above computational setups, the bulk temperature-dependent free energies were computed by adding zero-point energy and thermal corrections to the electronic total energy:

F (T )  Eelec  Fthermal (T ) Fthermal 

1   (q)  kBT q ln(1   (q) / kBT ) 2 q

Here, the thermal correction Fthermal (T ) was evaluated within the harmonic approximation,51 using the q-dependent phonon frequencies (  (q) ). Force constants were calculated utilizing supercells of at least 10 Å length in each dimension, which are large enough to converge the final thermal corrections within 0.1 meV/atom (see Supporting Information for details of supercells used), and the q-dependent phonon frequencies  (q) were obtained using the PHONOPY code suite,52 using the force constants computed by density functional perturbation theory.53 All calculations were carried out at 670 K, which is within the relevant annealing temperature range during synthesis.54

7

ACS Paragon Plus Environment

Chemistry of Materials

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 8 of 24

Slab models of surfaces were constructed by inserting a 15 Å vacuum gap into the supercell. Cleavage energies were computed by simply taking the difference between the energy of slab with the vacuum gap present and the energy of the corresponding supercell without vacuum. The cleavage energy is the sum of the surface energies of the two newly formed, fully relaxed surfaces (which may or may not be the same) and can be considered a metric for the stability of the surfaces created. Dipole energy corrections (all with magnitudes less than 5 meV/atom) were applied for polar surfaces and all cleavage energies were converged to within 0.01 J/m2 with respect to the slab thickness. Se addition effects on cleavage energies were investigated for kesterite and stannite phases containing 50% Se. First, a 16-atom unit cell containing four Se and four S sites was constructed utilizing the special quasirandom structures (SQS) method, which approximately reproduces correlations of the Se positions in a random solid solution55,

56

(the pair correlation function of the SQS cell matches the

random alloy within 5 Å). Then supercells and slabs were built based on the SQS unit cell analogous to the pure sulfide case, and the cleavage energies were then computed. The CZTS/ZnS interface models were constructed by adding the CZTS slab on top of the ZnS(001) slab, followed by insertion of a 15 Å vacuum gap. Since we mean to simulate a CZTS thin film growing on a ZnS substrate, we scaled the CZTS slab structure to adopt the ZnS cell parameter, inducing 0.2% lateral strain. All interface structures were fully relaxed using the same level of theory and numerical parameters introduced above. A 0.01 eV/Å convergence threshold was used for all forces on atoms during geometry optimizations. The relative energies ( E ) of CZTS slabs deposited on ZnS (001) were calculated using the bare ZnS substrate and ground state bulk CZTS (i.e., the kesterite phase) as common reference states. This relative energy metric is given by:

E  ( Einter  Esub  ekesterite nczts ) / nczts 8

ACS Paragon Plus Environment

Page 9 of 24

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

Chemistry of Materials

Here, Einter is the energy of the entire structurally relaxed CZTS (kesterite or stannite)/ZnS interface slab, E sub is the energy of the structurally relaxed, isolated ZnS substrate slab maintaining the same lattice vectors as the interface slab, ekesterite is the ground-state equilibrium bulk CZTS (kesterite phase) energy per formula unit, and nczts is the number of formula units in the CZTS layer. Essentially, this relative energy metric contains two contributions: the cleavage energy to create a kesterite slab or a stannite slab; and the adhesion energy between CZTS slab and the ZnS substrate. Both terms will contribute to the relative stability of the interface structures and should be included in the metric. In this equation, it is important to note that we always adopt a bare ZnS substrate surface and bulk kesterite CZTS as our reference states, regardless of which phase of CZTS we are studying in the interfaces. By using consistent reference states, the relative energies of different interface structures containing different phases of CZTS are directly comparable. Using this definition, we are able to directly quantify the relative stability of a CZTS layer (for different phases) coated on a ZnS substrate. To gain deeper insight into the contributions of the cleavage energy and the adhesion energy to the interface energy, we also compute the adhesion energy between the two slabs:

Ead  ( Einter  Esub  Eczts  slab ) / S Here, Eczts  slab is the energy of the optimized free-standing CZTS slab with bulk CZTS lattice constants and S is the surface area. Note that for different interfaces, the reference free-standing CZTS slab is different, so the Ead metric does not directly reflect the relative stability of a particular interface system. In principle, this metric only represents the adhesion contributions to the relative energy, but does not include the cleavage energy contribution.

RESULTS AND DISCUSSION 9

ACS Paragon Plus Environment

Chemistry of Materials

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 10 of 24

Crystal Structures In this work, we aim to explore the phase stabilities of CZTS, so we will benchmark our methodology using mainly ground state calculations, such as crystal parameters and thermodynamic properties; we will not study excited states properties such as the optical band gap, which requires a higher level of theory (such as GW). To verify the validity of the Sn(IV) U parameter and our methodology, we first compare the lattice parameters and bulk moduli of SnO2 and kesterite CZTS computed using PBE+U(+D2) and the pure PBE functional with experimental data (see Table 1). Table 1. SnO2 and pure sulfide CZTS lattice parameters (a and c) and bulk modulus (B).

SnO2

CZTS (kesterite)

a

a (Å)

c (Å)

B (GPa)

Expa

4.737

3.185

212.3

PBE

4.827

3.243

168.7

PBE+U

4.698

3.137

183.7

PBE+U+D2

4.662

3.13

205.8

Expb

5.426-5.435

10.81-10.848

-

PBE

5.475

10.938

69.7

PBE+U

5.444

10.866

66.5

PBE+U+D2

5.348

10.681

78.6

References 45, 57; b Reference 58

Simple DFT-PBE overestimates the cell parameters for both materials and greatly underestimates the bulk modulus of SnO2, while the predictions of PBE+U are superior to PBE for both materials. The D2 correction slightly overcorrects the cell parameters compared to the PBE+U method, but greatly improves the bulk modulus of SnO2. We note that due to the large nuclear charge of Sn(IV), the d-orbitals of Sn(IV) are more contracted, leading to stronger intraatomic (on-site) Coulomb and exchange interactions and larger Ueff compared to 10

ACS Paragon Plus Environment

Page 11 of 24

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

Chemistry of Materials

Cu(I) and Zn(II).44 So the U correction leads to much larger effects (2.7-3.3%) on the lattice constants of SnO2 compared to Cu2O (0.4%) and ZnO (1.3-1.4%). The lattice constants of CZTS are mainly determined by the sizes of the largest ions in the lattice, which are Cu and S. Therefore, the CZTS lattice constants only change by 0.7%, similar to the Cu2O case. We found that methods missing dispersion lead to completely incorrect geometries for layered materials, such as SnS2, since the interactions between layers are dominated by dispersion. Hence, it is critical to include the D2 correction so that all compounds considered can be treated on the same footing. Therefore, all subsequent calculations that follow here use the PBE+U+D2 method. Phase Stability To investigate the thermodynamics of CZTS synthesis, we calculated the free energies of the three phases (kesterite, stannite, and wurtzite) of bulk CZTS, as well as various relevant secondary compounds, namely, ZnS, CuS, Cu2S, SnS, SnS2, and Cu2SnS3. All results are listed in Table 2. Table 2. Thermally corrected free energies (kJ/mol/formula-unit) of CZTS and relevant secondary phases. All values are computed at T = 670 K. Compound

Phase

Eelec

Fthermal

F

SnS2

Berndtite

-1410

-45

-1454

SnS

Herzenbergite

-945

-44

-989

CuS

Covellite

-737

-33

-770

Cu2S

Low-chalcocitea

-1025

-64

-1089

ZnS

Zinc blende

-684

-30

-714

ZnS

Wurtzite

-684

-31

-716

Cu2SnS3

Monoclinic

-2474

-106

-2581

Cu2ZnSnS4

Kesterite

-3171

-135

-3306

Cu2ZnSnS4

Stannite

-3169

-135

-3303

11

ACS Paragon Plus Environment

Chemistry of Materials

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

Cu2ZnSnS4

Wurtzite

-3166

-135

Page 12 of 24

-3301

a

We used the low-temperature phase because high-chalcocite has partially occupied sites, leading to a more complicated unit cell which is too large to model. First, we confirm that the most stable CZTS phase is the kesterite phase. However, the free energy differences are only ~3 and ~5 kJ/mol/formula-unit compared to the stannite and wurtzite phases, respectively. Hence, it is not surprising that it is a challenge to obtain pure kesterite in experiments as the other phases are almost as energetically favorable to form. The data listed in Table 2 allow us to also evaluate the thermodynamics of different CZTS formation reactions, all of which are predicted to be exoergic at 670 K: the F of the reaction of

is -47 kJ/mol, the F

is -61 kJ/mol, and the F of

is -9

kJ/mol. We used wurtzite ZnS data in these calculations since wurtzite ZnS is predicted to be the ground state in our calculations. The experimental ground state for ZnS is the zinc blende phase, contrary to the calculations. However, the predicted energy difference between the two phases of ZnS is only ~2 kJ/mol, smaller than the typical error of DFT approximations. Therefore all the thermodynamic predictions will lead to essentially the same result no matter which phase of ZnS we use. Thus there is a strong thermodynamic driving force for CZTS formation, which should aid its synthesis. Previous computations23 reported that thermal corrections Fthermal (T ) significantly affect defect formation energies. However, we find that the thermal corrections for all three phases of CZTS are identical (see Table 2), so thermal effects are not essential to our analysis here and thus are neglected in the rest of this study for simplicity. Similar data have been reported by several other groups,22, 25 computed using either the PW91 or PBE GGA exchange-correlation functionals without any thermal corrections. While roughly generating the similar trends, the pure GGA functionals yield qualitatively incorrect predictions of the thermodynamics of certain reactions. For example, according to the data reported in reference 25, the reaction free energy for

is 13.5 kJ/mol, which is endothermic. The 12

ACS Paragon Plus Environment

Page 13 of 24

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

Chemistry of Materials

experimental standard reaction enthalpy is exothermic, -26.3 kJ/mol. Notably, our calculation predicts a qualitatively correct value of -16.1 kJ/mol; the error of about 10 kJ/mol is to be expected considering the typical accuracy of density functional approximations. The same reaction energy calculated using the hybrid HSE06 functional is -25.1 kJ/mol,30 in excellent agreement with the experimental value, albeit with a two orders of magnitude higher computational expense. Both the HF terms in HSE06 and the U correction terms we adopt here alleviate the self-interaction error and achieve better accuracy. Thus this reaction enthalpy comparison highlights the importance of self-interaction error corrections for first-row transition metal sulfides and validates our choice of methodology. To summarize this section, we computed thermally corrected free energies of different bulk CZTS phases and various secondary compounds. Calculated relative stabilities agree with previous studies, showing that the kesterite phase is the most stable, even though it is only favored by a few kilojoules per mole. This result indicates that even though the formation of CZTS is exothermic, we need special strategies to further stabilize the kesterite phase such that crystals of better quality can be obtained. CZTS Cleavage Energies As stated in the introduction, we aim to exploit surface properties to help stabilize the kesterite phase in nanoscale structures. In order to find surfaces that may facilitate formation of kesterite CZTS, we evaluated the cleavage energies of all symmetry-unique, low-index surfaces for all three phases of CZTS. The structures of all calculated surfaces can be found in Figure 2 and all cleavage energies are listed in Table 3 and Table 4. Table 3. Cleavage energies for kesterite and stannite phases of CZTS, in J/m2. Pure sulfide

50% Se

Surface

Kesterite

Stannite

Kesterite

Stannite

(001)-aa

1.22

1.49

0.92

1.28

13

ACS Paragon Plus Environment

Chemistry of Materials

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 24

(001)-b

1.61

1.29

1.27

1.15

(100)

1.22

1.39

1.07

1.16

(110)

0.78

0.79

0.66

0.65

(112)

0.78

0.77

0.70

0.69

(102)

0.77

0.68

0.69

0.65

a

Two different cleavage planes exist in this particular surface normal direction, denoted by letters (a) and (b), respectively (see Figure 2).

Table 4. Cleavage energies for the wurtzite phase of pure sulfide CZTS, in J/m2. Surface

Pure Sulfide Wurtzite

(001)

0.78

(010)-aa

0.36

(010)-b

0.95

(100)

0.78

(210)

0.68

a

Two different cleavage planes exist in this particular surface normal direction, denoted by letters (a) and (b), respectively (see Figure 2).

14

ACS Paragon Plus Environment

Page 15 of 24

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

Chemistry of Materials

Figure 2. Structures for: (a) surfaces for the kesterite and stannite phases, (b) surfaces for the wurtzite phase and (c) optimized structure of the interface between the metal-terminated ZnS (001) substrate and the four bilayers of the kesterite (001)-a slab model. Red lines indicate cleavage planes identified with Miller indices in black. By comparing Table 3 and Table 4, we observe that the wurtzite phase of CZTS tends to have lower cleavage energies. Especially the (010)-a surface of wurtzite CZTS is much more stable than all lowindex surfaces of kesterite and stannite. Thus, when the crystal particle is small enough that the surface energy dominates the energetics, CZTS will favor formation of the wurtzite phase. This result agrees with the fact that most wurtzite CZTS samples that have been made are predominantly nanocrystalline.59-68 To suppress formation of wurtzite will require some means of bypassing formation of nanocrystals or to convert them once formed. We consider the former in what follows. For both kesterite and stannite, the (110), (112), and (102) surfaces have significantly lower cleavage energies than the (100) and (001) surfaces. Indeed, in a typical powder X-ray diffraction spectrum of CZTS,69, 70 the peaks corresponding to the interplanar distances of (110), (112), and (102) are much 15

ACS Paragon Plus Environment

Chemistry of Materials

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 24

more distinctive than the (100)/(001) surfaces, suggesting a possible preference for growth along those directions. Subsequent experimental studies71-80 confirmed that the most preferable growth direction for the CZTS films is always [112] on either Mo or glass substrates, in qualitative agreement with our calculated cleavage energies. For these three low-energy surfaces, we find no essential difference in stability between the kesterite and stannite phases. On the other hand, the kesterite (001)-a and (100) cleavage energies are up to 0.27 J/m2 lower than the same surfaces of stannite. However, these two surfaces have higher energies compared to the other low-index surfaces and therefore are disfavored during growth. Therefore, to exploit surface thermodynamics to stabilize pure CZTS in the kesterite phase, we need to find a way to stabilize the “beneficial (001) and (100) surfaces” that favor kesterite over stannite, while disfavoring wurtzite nanocrystal formation. In other words, we need to find the means to encourage crystal growth - of thin films rather than nanocrystals - along the [001] and [100] directions via providing proper crystallization templates. We return to this point momentarily. We next investigated the effects of Se additions on the cleavage energies of kesterite and stannite CZTS to see if adding Se might be a way to stabilize kesterite; the results are shown in Table 3. Se alloying is predicted to stabilize all surfaces of both phases, while the effect is most prominent for the (001) and (100) surfaces. With 50% Se, the energies for the (001) and (100) surfaces are lowered by ~0.2-0.3 J/m2 whereas for the other three directions only a ~0.1 J/m2 decrease is predicted. Se additions therefore could be effective at stabilizing the “beneficial surfaces”, although this effect is not strong enough to alter the trend of the most stable surfaces being (110)/(102)/(112) rather than (100)/(001). Based on these results, we instead explore whether the stability of (100) and (001) kesterite CZTS surfaces can be enhanced at interfaces with appropriate substrates that might be used as growth templates in the synthesis process. CZTS/ZnS(001) Interfaces 16

ACS Paragon Plus Environment

Page 17 of 24

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

Chemistry of Materials

We consider zinc blende ZnS, which topologically is a parent compound for both the kesterite and stannite CZTS structures. The computed cell parameter for ZnS is 5.34 Å, which is very close to the kesterite cell parameter of 5.35 Å. ZnS has the same stacking pattern as kesterite/stannite CZTS and very minor lattice mismatch, making its (001) surface a perfect crystallization template for CZTS (001) and (100) surfaces. We consider here only the two tetragonal phases and not the wurtzite phase for two reasons: 1) wurtzite CZTS has highest bulk free energy and 2) the hexagonal lattice of wurtzite CZTS does not match well the cubic lattice of the zinc blende ZnS (001) substrate. Thus it is highly unlikely that wurtzite will be a competitive phase on this particular substrate. Which phase is more stable on the substrate, kesterite or stannite? If the kesterite phase is more stable, is there any extra stability provided by the substrate and the thin film geometry, compared to the bulk situation? To answer these questions, we calculated the relative energy of CZTS overlayers deposited on the ZnS (001) surface. The ZnS substrate was built with four ion bilayers, with each of the bilayers composed of one layer of cations and one layer of anions. We tested two models of the CZTS overlayer with different thicknesses, containing two (corresponding to half a unit cell in the (001) direction) or four bilayers (corresponding to one unit cell in the (001) direction) respectively, representing different stages of crystal growth. Both models maintained the correct stoichiometry of the bulk CZTS. The interfaces can be further divided into two different types: sulfur-terminated (S-term) CZTS in contact with metal-terminated (M-term) substrate, or M-term CZTS with S-term substrate. The relative energies computed for all these systems are summarized in Table 5. Table 5. Interface relative energies (in kJ/mol/formula-unit). Overlayer phase and orientation

Kesterite (001)-a

S-terminated substrate

Zn-terminated substrate

2-bilayer

4-bilayer

2-bilayer

4-bilayer

48.96

43.12

-129.76

-49.20

17

ACS Paragon Plus Environment

Chemistry of Materials

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 24

Kesterite (001)-b

79.28

39.44

-0.40

17.04

Stannite (001)-a

118.48

83.68

-98.48

13.44

Stannite (001)-b

69.76

35.92

-79.52

22.40

Kesterite (100)

52.00

24.00

-77.76

-40.88

Stannite (100)

9.12

6.16

-73.44

-38.08

For the S-terminated substrate, the most stable configuration is the stannite (100) surface in both 2bilayer and 4-bilayer cases, which is not desirable. However, for the M-terminated substrate, kesterite (001)-a becomes the most stable geometry and is energetically favored by ~30 kJ/mol/formula-unit (in the 2-bilayer case) and ~11 kJ/mol/formula-unit (in the 4-bilayer case) compared to the most stable stannite configuration. Considering the bulk energy difference between the two phases is only 3 kJ/mol/formula unit, we can conclude that the zinc-terminated ZnS (001) substrate provides significant extra stability for the kesterite phase, especially in the early stage of crystallization (the 2-bilayer case). Furthermore, the adhesion energies between the kesterite(001)-a slab and the Zn-terminated ZnS slab are very attractive (-3.09 J/m2 for the 2-bilayer case and -3.02 J/m2 for the 4-bilayer case), showing that the two slabs adhere well to each other and formation of the interface is naturally favored. We thus predict that use of a zinc-terminated ZnS (001) substrate may aid in synthesis of pure kesterite CZTS crystals. It is interesting to further compare the adhesion energies of the ZnS/kesterite(001)-a interface with the ZnS/stannite(001)-a system. For the 2-bilayer case, we obtained ΔEad = -3.09 J/m2 for kesterite(001)-a and -3.44 J/m2 for stannite(001)-a. Apparently the stannite phase should be favored over the kesterite phase from the perspective of adhesion energies, in contrast with the trend observed in the total relative energy metric. Therefore, the extra stability of the kesterite(001)-a interface is controlled by the lower cleavage energies of the kesterite(001) slab (see Table 3). The only advantage of the ZnS(001) substrate is to provide an appropriate template for the beneficial orientations, and we do not have to rely on any 18

ACS Paragon Plus Environment

Page 19 of 24

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

Chemistry of Materials

special chemical affinities between the ZnS substrate and the kesterite film. Hence, in later stages of growth, we expect the extra stability of the kesterite phase to remain as long as the orientation of the film surface is unchanged. We can consider the entire growth process to be a combination of many elementary thin-layer deposition steps. Within each elementary step, the formation of a new kesterite layer is always strongly favored over a new stannite layer. Therefore, the lower cleavage energy of the kesterite(001) surface provides an advantage for the kesterite phase during the entire deposition process. CONCLUSIONS In this study, we investigated the bulk, surface, and interface stabilities of three phases of CZTS (kesterite, stannite, and wurtzite) using PBE+U+D2 level of theory. The results show that the formation of CZTS is thermodynamically favored and that the different phases of CZTS are very close in energy, in accordance with previous studies. This near-degeneracy is likely responsible for secondary phase formation in CZTS films, introducing SRH recombination centers, lowering the band gap and the open circuit voltage, ultimately decreasing PV cell efficiency. The stannite phase is particularly problematic for the desired kesterite phase, as these two phases share the same lattice topology. Partial Cu/Zn disorder has been observed experimentally, indicating the instability of the ordered Cu/Zn sub-lattice, which is consistent with the quasi-degeneracy of the kesterite and stannite phases. Via computer simulations, we explored the feasibility of a strategy that exploits the surface and interface properties of CZTS to enhance kesterite’s relative stability. Our calculations indicate that the wurtzite phase of CZTS generally has the lowest cleavage energies and thus will be the most favored phase in nanocrystals where surface energetics dominates, consistent with recent measurements. No essential difference in stability was found between kesterite and stannite for all low energy surfaces that naturally emerge in current synthesis procedures. However, two surfaces ((001) and (100)) with slightly higher energies may provide extra stability for kesterite and thus might be exploitable as “beneficial surfaces”. Se additions could lower the formation energies of these “beneficial surfaces” but not very 19

ACS Paragon Plus Environment

Chemistry of Materials

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 20 of 24

significantly. Instead, we assessed whether a zinc blende ZnS (001) surface could serve as a successful template to induce selective growth of kesterite CZTS (001)/(100) surfaces. This substrate has two important features: 1) its lattice parameters are almost identical to the corresponding surfaces of kesterite/stannite CZTS and formation of wurtzite CZTS will be inhibited due to significant lattice mismatch and 2) we predict that the kesterite/ZnS(001) interface indeed has significantly higher stability than the stannite/ZnS(001) interface, provided the ZnS(001) substrate is zinc-terminated. Compared to the bulk, the energy difference between the two phases is significantly enhanced in the interface structure, especially in the initial stage of crystallization when surface properties play a more important role. We therefore propose that zinc-terminated zinc blende ZnS (001) or other similar substrates that stabilize “beneficial surfaces” during synthesis should be explored to see if indeed better quality CZTS crystals are formed.

Supporting Information Available Table S1 contains the sizes and k-point setups of the simulation cells we used for all secondary phases, surfaces and interfaces.

ACKNOWLEDGMENTS E.A.C. thanks the U.S. Department of Energy, Office of Science, Basic Energy Sciences under Award No. DE-SC0002120 for funding this project.

REFERENCES 1.

NREL research cell efficiency records chart. http://www.nrel.gov/ncpv/images/efficiency_chart.jpg.

2.

Delbos, S. EPJ Photovoltaics 2012, 3, 35004.

3.

Green, M. A.; Emery, K.; Hishikawa, Y.; Warta, W.; Dunlop, E. D. Prog. photovolt. 2014, 22, 1-9.

4.

Service, R. F. Science 2014, 344, 458.

20

ACS Paragon Plus Environment

Page 21 of 24

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

5.

Chemistry of Materials

Todorov, T. K.; Tang, J.; Bag, S.; Gunawan, O.; Gokmen, T.; Zhu, Y.; Mitzi, D. B. Advanced Energy Materials 2013, 3,

34-38. 6.

Halliday, D. P.; Claridge, R.; Goodman, M. C. J.; Mendis, B. G.; Durose, K.; Major, J. D. J. Appl. Phys. 2013, 113,

223503. 7.

Phuong, L. Q.; Okano, M.; Yamada, Y.; Nagaoka, A.; Yoshino, K.; Kanemitsu, Y. Appl. Phys. Lett. 2014, 104, 081907.

8.

Schorr, S. Sol Energ Mat Sol C 2011, 95, 1482-1488.

9.

Chen, S. Y.; Gong, X. G.; Walsh, A.; Wei, S. H. Appl. Phys. Lett. 2009, 94, 041903.

10. Paier, J.; Asahi, R.; Nagoya, A.; Kresse, G. Phys. Rev. B 2009, 79, 115126. 11. Persson, C. J. Appl. Phys. 2010, 107, 053710 12. Siebentritt, S.; Schorr, S. Prog. photovolt. 2012, 20, 512-519. 13. Lu, X. T.; Zhuang, Z. B.; Peng, Q.; Li, Y. D. Chem. Commun. 2011, 47, 3141-3143. 14. Mendis, B. G.; Shannon, M. D.; Goodman, M. C.; Major, J. D.; Claridge, R.; Halliday, D. P.; Durose, K. Prog Photovoltaics 2014, 22, 24-34. 15. Christensen, A.; Carter, E. A. Phys. Rev. B 1998, 58, 8050-8064. 16. Scanlan, C. M.; Gajdardziska-josifovska, M.; Aita, C. R. Appl. Phys. Lett. 1994, 64, 3548-3550. 17. Darken, L. S.; Gurry, R. W. J. Am. Chem. Soc. 1945, 67, 1398-1412. 18. Galloway, H. C.; Benitez, J. J.; Salmeron, M. Surf. Sci. 1993, 298, 127-133. 19. Wang, W.; Zhang, H.; Wang, W. H.; Zhao, A. D.; Wang, B.; Hou, J. G. Chem. Phys. Lett. 2010, 500, 76-81. 20. Gribb, A. A.; Banfield, J. F. Am. Mineral. 1997, 82, 717-728. 21. Tolbert, S. H.; Alivisatos, A. P. Science 1994, 265, 373-376. 22. Nagoya, A.; Asahi, R.; Wahl, R.; Kresse, G. Phys. Rev. B 2010, 81, 113102. 23. Kosyak, V.; Amiri, N. B. M.; Postnikov, A. V.; Scarpulla, M. A. J. Appl. Phys. 2013, 114, 124501. 24. Chen, S. Y.; Walsh, A.; Gong, X. G.; Wei, S. H. Adv. Mater. 2013, 25, 1522-1539. 25. Chen, S. Y.; Yang, J. H.; Gong, X. G.; Walsh, A.; Wei, S. H. Phys. Rev. B 2010, 81, 245204 26. Perdew, J. P.; Wang, Y. Phys. Rev. B 1992, 45, 13244-13249. 27. Perdew, J. P.; Burke, K.; Ernzerhof, M. Phys Rev Lett 1996, 77, 3865-3868. 28. Krukau, A. V.; Vydrov, O. A.; Izmaylov, A. F.; Scuseria, G. E. J. Chem. Phys. 2006, 125, 224106. 29. Nagoya, A.; Asahi, R.; Kresse, G. J. Phys.- Condens. Matter 2011, 23, 404203. 30. Han, D.; Sun, Y. Y.; Bang, J.; Zhang, Y. Y.; Sun, H. B.; Li, X. B.; Zhang, S. B. Phys. Rev. B 2013, 87, 155206. 31. Anisimov, V. I.; Aryasetiawan, F.; Lichtenstein, A. I. J. Phys.- Condens. Matter 1997, 9, 767-808. 32. Dudarev, S. L.; Botton, G. A.; Savrasov, S. Y.; Humphreys, C. J.; Sutton, A. P. Phys. Rev. B 1998, 57, 1505-1509. 33. Mosey, N. J.; Liao, P.; Carter, E. A. J. Chem. Phys. 2008, 129, 014103. 34. Mosey, N. J.; Carter, E. A. Phys Rev B 2007, 76, 155123. 35. Kresse, G.; Hafner, J. Phys. Rev. B 1993, 47, 558-561. 36. Kresse, G.; Hafner, J. Phys. Rev. B 1994, 49, 14251-14269. 37. Kresse, G.; Furthmuller, J. Comput. Mater. Sci. 1996, 6, 15-50. 38. Kresse, G.; Furthmüller, J. Phys Rev B 1996, 54, 11169-11186.

21

ACS Paragon Plus Environment

Chemistry of Materials

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 22 of 24

39. Blochl, P. E. Phys. Rev. B 1994, 50, 17953-17979. 40. Kresse, G.; Joubert, D. Phys. Rev. B 1999, 59, 1758-1775. 41. Grimme, S. J. Comput. Chem. 2006, 27, 1787-1799. 42. Murnaghan, F. D. Proc. Natl. Acad. Sci. U.S.A. 1944, 30, 244-247. 43. Monkhorst, H. J.; Pack, J. D. Phys. Rev. B 1976, 13, 5188-5192. 44. Yu, K.; Carter, E. A. J. Chem. Phys. 2014, 140, 121105. 45. Baur, W. H. Acta Crystallogr 1956, 9, 515-520. 46. Hehre, W. J.; Ditchfie.R; Pople, J. A. J. Chem. Phys. 1972, 56, 2257. 47. Metz, B.; Stoll, H.; Dolg, M. J. Chem. Phys. 2000, 113, 2563-2569. 48. Peterson, K. A. J. Chem. Phys. 2003, 119, 11099-11112. 49. Hay, P. J.; Wadt, W. R. J Chem Phys 1985, 82, 270-283. 50. Liao, P. L.; Keith, J. A.; Carter, E. A. J Am Chem Soc 2012, 134, 13296-13309. 51. McQuarrie, D. A. Statistical Mechanics. University Science Books: 2000. 52. Togo, A.; Oba, F.; Tanaka, I. Phys. Rev. B 2008, 78, 134106. 53. Gonze, X.; Lee, C. Phys. Rev. B 1997, 55, 10355-10368. 54. Jiang, M.; Yan, X. Cu2ZnSnS4 Thin Film Solar Cells: Present Status and Future Prospects. In Solar Cells - Research and Application Perspectives, Morales-Acevedo, A., Ed. InTech: 2013. 55. van de Walle, A.; Tiwary, P.; de Jong, M.; Olmsted, D. L.; Asta, M.; Dick, A.; Shin, D.; Wang, Y.; Chen, L. Q.; Liu, Z. K. Calphad 2013, 42, 13-18. 56. Zunger, A.; Wei, S. H.; Ferreira, L. G.; Bernard, J. E. Phys. Rev. Lett. 1990, 65, 353-356. 57. Chang, E.; Graham, E. K. J. Geophys. Res. 1975, 80, 2595-2599. 58. Mitzi, D. B.; Gunawan, O.; Todorov, T. K.; Wang, K.; Guha, S. Sol. Energy Mater. Sol. Cells 2011, 95, 1421-1436. 59. Zou, Y.; Su, X.; Jiang, J. J. Am. Chem. Soc. 2013, 135, 18377-18384. 60. Zhao, Y.; Qiao, Q.; Zhou, W. H.; Cheng, X. Y.; Kou, D. X.; Zhou, Z. J.; Wu, S. X. Chem. Phys. Lett. 2014, 592, 144148. 61. Thompson, M. J.; Ruberu, T. P. A.; Blakeney, K. J.; Torres, K. V.; Dilsaver, P. S.; Vela, J. J. Phys. Chem. Lett. 2013, 4, 3918-3923. 62. Liao, H. C.; Jao, M. H.; Shyue, J. J.; Chen, Y. F.; Su, W. F. J. Mater. Chem. A 2013, 1, 337-341. 63. Li, M.; Zhou, W. H.; Guo, J.; Zhou, Y. L.; Hou, Z. L.; Jiao, J.; Zhou, Z. J.; Du, Z. L.; Wu, S. X. J. Phys. Chem. C 2012, 116, 26507-26516. 64. Li, C. Y.; Ha, E. N.; Wong, W. L.; Li, C. L.; Ho, K. P.; Wong, K. Y. Mater. Res. Bull. 2012, 47, 3201-3205. 65. Kong, J.; Zhou, Z. J.; Li, M.; Zhou, W. H.; Yuan, S. J.; Yao, R. Y.; Zhao, Y.; Wu, S. X. Nanoscale Res. Lett. 2013, 8, 464. 66. Kang, C. C.; Chen, H. F.; Yu, T. C.; Lin, T. C. Mater. Lett. 2013, 96, 24-26. 67. Jao, M. H.; Liao, H. C.; Wu, M. C.; Su, W. F. Jpn. J. Appl. Phys. 2012, 51, 10NC30. 68. Cattley, C. A.; Cheng, C.; Fairclough, S. M.; Droessler, L. M.; Young, N. P.; Warner, J. H.; Smith, J. M.; Assender, H. E.; Watt, A. A. R. Chem. Commun. 2013, 49, 3745-3747.

22

ACS Paragon Plus Environment

Page 23 of 24

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

Chemistry of Materials

69. Guo, Q. J.; Hillhouse, H. W.; Agrawal, R. J. Am. Chem. Soc. 2009, 131, 11672-+. 70. Tanaka, K.; Moritake, N.; Uchiki, H. Sol. Energy Mater. Sol. Cells 2007, 91, 1199-1201. 71. Sun, L.; He, J.; Kong, H.; Yue, F. Y.; Yang, P. X.; Chu, J. H. Sol. Energy Mater. Sol. Cells 2011, 95, 2907-2913. 72. Arba, Y.; Rafi, M.; Tchognia, H.; Hartiti, B.; Ridah, A.; Thevenin, P. J. Optoelectron. Adv. Mater. 2013, 15, 1200-1203. 73. Byeon, M. R.; Chung, E. H.; Kim, J. P.; Hong, T. E.; Jin, J. S.; Jeong, E. D.; Bae, J. S.; Kim, Y. D.; Park, S.; Oh, W. T.; Huh, Y. S.; Chang, S. J.; Lee, S. B.; Jung, I. H.; Hwang, J. Thin Solid Films 2013, 546, 387-392. 74. Chalapathi, U.; Uthanna, S.; Raja, V. S. J. Renew. Sustain. Ener. 2013, 5, 031610. 75. Kheraj, V.; Patel, K. K.; Patel, S. J.; Shah, D. V. J. Cryst. Growth 2013, 362, 174-177. 76. Seboui, Z.; Cuminal, Y.; Kamoun-Turki, N. J. Renew. Sustain. Ener. 2013, 5, 023113. 77. Wang, J. S.; Li, S.; Cai, J. J.; Shen, B.; Ren, Y. P.; Qin, G. W. J. Alloys Compd. 2013, 552, 418-422. 78. Patel, M.; Mukhopadhyay, I.; Ray, A. Semicond. Sci. Technol. 2013, 28, 055001. 79. Zhang, J.; Long, B.; Cheng, S. Y.; Zhang, W. B. Int. J. Photoenergy 2013, 2013, 986076. 80. Jiang, X.; Shao, L. X.; Zhang, J.; Li, D.; Xie, W.; Zou, C. W.; Chen, J. M. Surf. Coat. Technol. 2013, 228, S408-S411.

23

ACS Paragon Plus Environment

Chemistry of Materials

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

TOC graphic:

24

ACS Paragon Plus Environment

Page 24 of 24