Devendra
Tiwari
a,
Tristan
Koehler
b,
Reiner
Klenk
b and
David J.
Fermin
*a
aSchool of Chemistry, University of Bristol, Cantocks Close, Bristol BS8 1TS, UK. E-mail: David.Fermin@bristol.ac.uk
bHelmholtz-Zentrum Berlin für Materialien und Energie, Hahn-Meitner-Platz 1, D-14109 Berlin, Germany
First published on 27th March 2017
High quality microcrystalline tetragonal p-Cu2SnS3 (CTS) films are prepared by spin-coating a single precursor of metal salts and thiourea on to 5 × 5 cm2 Mo substrates. The best of the solar cells completed on these films with a substrate configuration: glass/Mo/CTS/CdS/i-ZnO/Al:ZnO/Ni–Al and a total area of 0.5 cm2 displays an open-circuit voltage of 206 mV, a fill factor of 34.5%, a short circuit current density of 27.6 mA cm−2 and a power conversion efficiency of 1.9% under simulated AM1.5 illumination. This the best performance reported for such solar architectures obtained by solution processing, with dispersion below 20% for 24 devices. For the first time, the key factors limiting CTS device efficiency are quantitatively established based on temperature dependent current–voltage curves and impedance spectroscopy, namely: (i) carrier recombination at the CTS/CdS interface, (ii) MoS2 non-ohmic back contact, formed due to sulfurization of the top Mo layer, with a barrier height of up to 216 ± 14 meV and (iii) the presence of two trap levels with activation energies 41 ± 0.4 meV and 206 ± 7 meV. The shallower trap is linked Cu vacancies, while the deeper trap is associated with Sn in Cu antisite defects based on DFT supercell calculations.
These studies have reinvigorated research in similar materials from the Cu–Sn–S family. The Cu–Sn–S phase diagram is rich with many stable phases: Cu2SnS3, Cu3SnS4, Cu4SnS4, Cu2Sn3S7, and Cu5Sn2S7; among which only Cu2SnS3 (CTS) has shown photovoltaic potential.9 Theoretically, CTS is shown to have a wide stability range and be devoid of Fermi level pinning, while the other forms exhibit poor hole-mobility, Fermi level pinning or high conductivity.10 Cu2SnS3 has a band gap of approximately 1 eV, an absorption coefficient above 105 cm−1 and conductivity between 0.5 and 10 S cm−1, which correspond to a hole concentration of ∼1018 cm−3 and mobility between 1 and 80 cm2 V−1 s−1.11,12 The present state of the art includes solar cells, featuring materials obtained by physical vapor deposition, with efficiency up to 4.63% for pure CTS and up to 6% for Ge-alloyed CTS.13,14 So far, reports have primarily focused on the material synthesis and assessment of PV properties of CTS and very little knowledge has been gathered with regard to the phenomena limiting the power device conversion.
A variety of physical deposition methods have been employed to prepare CTS including sputtering,11,14–16 thermal or e-beam evaporation,11,13 and pulsed laser deposition.11,17 Wet chemical methods, involving single or sequential steps, have also been implemented such as spin/dip coating,12,18,19,22–25 electrodeposition20,21 and nanoparticle inks.18,19,22 Solution based methods are rather appealing from the manufacturing point of view, particularly in terms of the PV module payback time.7,8,26–29 In the context of CTS based devices, the vast majority of the work has focused on physical vapor deposition, with only two studies employing solution processing generating efficiencies up to 2.1%.23,24
In this contribution, we assessed the PV performance of phase pure CTS films obtained by spin coating a single solution based precursor. The interactions between the metal and sulfur components of the precursor are investigated by IR spectroscopy, while thermally induced crystallization and stability of the CTS phase are probed by thermo-gravimetric analysis and differential scanning calorimetry. X-ray diffraction and Raman spectroscopy confirm the generation of polycrystalline CTS films with a tetragonal phase. Solar cells in the superstrate configuration: glass/Mo/CTS/CdS/i-ZnO/Al:ZnO/Ni–Al were fabricated featuring power conversion efficiencies close to 2%, which in itself is the highest reported for solution processed solar cells in the industrially adopted substrate configuration. Temperature dependent J–V and impedance measurements show that carrier recombination at the CTS/CdS junction is a dominant factor in the cell performance, as well as the presence of two different bulk defect states with activation energies 41 ± 0.4 meV and 206 ± 7 meV. Employing DFT supercell calculations of defect formation energies identifies Sn on Cu antisites as the key bulk recombination state.
The enthalpy of formation of different point defects is calculated employing DFT in a 96 atom supercell. Calculations are performed with the CASTEP code employing a DFT pseudopotential approach. A generalized gradient approximation functional: PBESOL with ultrasoft pseudopotentials, an energy cutoff of 500 eV and a Monkhorst–Pack grid with a spacing of ∼0.02 Å is used for geometry optimization and energy calculations. The atomic positions are optimized using a BFGS scheme with convergence tolerances of 1.0 × 10−7 eV per atom for energy, 0.01 eV Å for maximum force, 0.02 GPa for maximum stress, and 5.0 × 10−4 Å for maximum displacement.
Fig. 1 (a) FT-IR spectra of pure thiourea (black) and the precursor (red); (b) X-ray diffraction with Rietveld fitting and (c) Raman spectrum of the CTS film on a Mo substrate. |
A Mo coated glass substrate is spin coated with the precursor solution and heated at 300 °C. After the desired thickness is achieved by repeating the last steps, the films are annealed in a S atmosphere at 550 °C. Further details of deposition and processing are provided in the Experimental section. Thermogravimetric analysis of the precursor is shown in Fig. S1,† featuring a weight loss between 180 °C and 320 °C which corresponds to the decomposition of the metal-thiourea precursor to the corresponding sulphide. A small mass loss is also observed between 380 °C and 600 °C, most probably associated with SnS or S. Based on these results, a thermal procedure is established involving a heating step at 300 °C after spin-coating followed by annealing at 550 °C.
Fig. 1b shows a characteristic X-ray diffraction (XRD) pattern of the CTS film on a Mo substrate. CTS has been synthesized in different polymorphs, mostly monoclinic, cubic and tetragonal.15,31–37 Cubic CTS has a zinc blende type of structure. The monoclinic form is essentially a superstructure of the cubic arrangement, while the tetragonal structure derives from the cubic form in case there is a random distribution of Cu and Sn. This canonical relation leads to a similar diffraction pattern for these different CTS phases.37 In order to establish the type of crystal structure, analysis of diffraction patterns employing Rietveld refinement is performed using the Fullprof suite.38 Values of the correlation coefficients Rwp and Rp of 6.48 and 8.26, respectively, indicate a high quality fit. The results suggest that CTS crystallizes in a stannite-like tetragonal polymorph (I2m) with lattice parameters a = 5.426(3) Å and b = 10.8375(1) Å.‡
A typical Raman spectrum of the CTS thin film is displayed in Fig. 1c. Raman spectroscopy is particularly useful for distinguishing the various polymorphs of CTS, with characteristic bands reported at 290 cm−1 and 352 cm−1 for monoclinic; at 303 cm−1 and 355 cm−1 for cubic and at 336 cm−1 and 351 cm−1 for the tetragonal forms.15,39,40 The spectrum of the thin film shows a broad peak with a shoulder between 300 and 380 cm−1, which can be deconvoluted into two Voigt functions centered at 334 and 352 cm−1. These two modes are consistent with a tetragonal lattice. Based on the XRD and Raman analysis, a schematic of the refined unit cell is shown in Fig. S2.† Each metal site is tetrahedrally coordinated to S atoms and vice versa with absence of a direct S–S bond. The Cu–S distance is estimated to be 2.322 Å while the distances between mixed metal atom sites M(I) and M(II), i.e. Wyckoff positions 4d and 2b, and S are 2.354 Å and 2.364 Å, respectively. The bond angles S–Cu–S, S–M(I)–S and S–M(II)–S vary between −0.62° and +1.16° around the tetrahedral angle. These results are in agreement with previous studies.22,23,32,35 The metal atom distribution on site M(I) is dominated by Cu (56.39% occupancy) while the M(II) site is Sn dominated (52.53% occupancy). The overall unit cell composition reflects a Cu/Sn ratio of 2.12, which is in agreement with the elemental ratio obtained from EDAX.
XPS analysis of the Cu 2p, Sn 3d and S 2p core levels is shown in Fig. 2. The binding energies of Cu 2p3/2 and Cu 2p1/2 are observed to be 933.2 and 953.1 eV, with a full width at half maximum of 1.90 eV and 2.24 eV, respectively, consistent with a Cu+ oxidation state.22,23,32,35 The absence of any satellite or shake-up peak around 942 eV further confirms the absence of the Cu2+ state. In sulfides, the tetragonal coordination of Cu typically promotes the +1 oxidation state, as the Cu d10s1 configuration is the most stable for this geometry. The binding energies of Sn 3d5/2 and 3d3/2 are 486.47 and 494.99 eV, respectively, closely matching those reported for SnS2 (i.e. +4 oxidation state).22,23,32,35 The S 2p3/2 and 2p1/2 are closely located at 161.26 and 162.37 eV, respectively, which is similar to the spectra of CuFeS2. The XPS data allow estimating the relative atomic ratios of Cu/Sn/S of 30.6%/14.7%/50.6%, while the ratios of C and O are 1.5% and 2.6%, respectively. These compositions are consistent with EDX analysis as well as the results from the XRD refinement.
The morphology of films as probed through scanning electron microscopy of the CTS film on a Mo coated substrate is presented in Fig. 3a. The films are homogenous and compact with grain sizes between 500 and 1100 nm, which are appropriate for device fabrication. Fig. 3b shows the cross-sectional image of the films which reveals a uniform growth of the adherent film with a thickness of 1.2 μm. The contrast between CTS and the Mo layer indicates the partial sulfurization of the Mo layer at the interface. Formation of this MoS2 would implicate the device performance which will be discussed in the later section. Fig. 3c shows the diffuse reflectance spectrum of the CTS film, featuring a large change in reflectance between 1050 and 1350 nm corresponding to the band-to-band transition. The inset in Fig. 3c corresponds to a modified Tauc's plot using the Kubelka–Munk transformation of the diffuse reflectance, from which an optical band gap of 1.1 eV can be estimated. This value is consistent with previous optical studies of tetragonal CTS.15,23,40
Four probe conductivity and Hall measurements of the CTS films revealed p-type conductivity with a resistivity and hole mobility of 2.226 Ω cm and 4.581 cm2 V−1 s−1, respectively. These values allow estimating a room temperature hole concentration of 6.12 × 1017 cm−3, which is appropriate for heterojunction solar cell applications. Many of the Cu–Sn–S phases suffer from high carrier concentration limiting their utility for PV devices. Theoretical studies have predicted CTS doping densities between 1015 and 1018 cm−3 due to native vacancies, although other types of defects can contribute as well.10 A combinatorial study of a wide range [Cu]/([Cu] + [Sn]) ratio concluded that values between 0.6 and 0.72 lead to a hole density close to 1018 cm−3.41 This is fully consistent with our findings, considering that the [Cu]/([Cu] + [Sn]) ratio is 0.67.
The performance of 0.5 cm2 devices with the substrate architecture glass/Mo/CTS/CdS/i-ZnO/Al:ZnO/Ni–Al, with no antireflective coating, under a standard simulated AM1.5 spectrum is displayed in Fig. 4. The CdS layer (70 nm) was grown by chemical bath deposition, while the i-ZnO and Al–ZnO layers (400 nm total) were deposited by RF sputtering. The J–V characteristics of the best cell in the dark and under illumination are shown in Fig. 4a, featuring a power conversion efficiency of 1.92%, an open-circuit voltage (VOC) of 206 mV, a fill factor (FF) of 34.5% and a short circuit current density, JSC of 27.6 mA cm−2. These figures of merit are the highest reported for solution processed CTS on a substrate architecture. The relatively narrow dispersion of key parameters for 24 solar cells is summarized in Table S2.†
The spectral response of the best cell is shown in Fig. 4b, showing a maximum external quantum efficiency (EQE) of 70% at around 560 nm. At shorter wavelengths the device performance is restricted by the CdS (540 nm edge) and ZnO layers (400 nm edge). A band gap of 1.1 eV can be estimated from the onset of the EQE spectrum, which corroborates the value estimated from diffuse reflectance (Fig. 3c). The integrated value of the photocurrent over the entire spectrum is found to be 25.8 mA, falling close to the JSC calculated from J–V curves.
Fig. 4c shows the temperature dependence of the key device performance in the range of 300 to 80 K. As the temperature increases, the VOC linearly decreases over large portions of the temperature range, while the JSC slightly increases reaching a maximum at 240 K. The non-monotonic temperature dependence of JSC is a manifestation of a non-linear series resistance, most probably associated with a MoS2 rectifying back contact.42 The device efficiency increases with decreasing temperature reaching a maximum of 4.8% at 120 K. The temperature dependence of the VOC can be expressed in terms of:42–44
(1) |
The origin of the interfacial recombination is most likely connected to the cliff like band alignment between CTS and CdS which has been measured using photoelectron spectroscopy.45 Recombination can also be linked to clusters of cation disorder and stacking faults, which has been recently theoretically postulated and experimentally observed employing electron microscopy.46,47 Further studies are required in order to clearly identify the nature of the interfacial recombination site associated with the VOC deficiency shown in Fig. 4c.
As discussed in the ESI,† quantitative analysis of carrier dynamics can be extracted from the temperature dependence of the device impedance spectra (Fig. S3†). A systematic analysis using different equivalent circuits reveals contributions from two RC time constants associated with the CdS/CTS and CTS/MO interfaces, as well as the dynamic responses of two defect sites (Fig. S3a–c†). The adopted equivalent circuit has been used in numerous impedance studies featuring Cu2ZnSnS4, CdTe, CuInSe2 and Cu(In,Ga)Se2 absorbers.42,43,48–50
The temperature dependent relaxation frequency (ωT = 1/RC) of the two defect states measured by impedance spectroscopy is displayed in Fig. 5a, based on the following expression:
(2) |
Both frequencies show a clear Arrhenius behavior over the entire temperature range investigated, providing activation energy values of EA,D1 = 41 ± 0.4 meV and EA,D2 = 206 ± 7 meV. Interestingly, EA,D1 is very similar to the values associated with a Cu vacancy in chalcopyrites and kesterite cells.42,43,48,50 The nature of EA,D2 is discussed further below. The temperature dependence of the back contact Rb is displayed in Fig. 5b. In this case, the temperature dependence can be described in terms of,42–44
(3) |
In order to identify the nature of the deep defect state exhibiting an activation barrier of 216 meV, we have performed DFT calculations on a 96 atom monoclinic CTS supercell displayed in Fig. 6a. The use of a monoclinic unit cell as a first approximation can be justified considering that the tetragonal polymorph of CTS is a special case of the cubic form with a random cation distribution, while the cubic form is in itself a superstructure of the monoclinic polymorph.37 Furthermore, the latter is the ground state structure of CTS at 0 K. Calculations on a tetragonal unit cell would require extensive computational resources which are beyond the scope of this work. Details of calculations are included in the Experimental section.
The formation energies of various types of neutral point defects are calculated at five different points in the Cu–Sn–S ternary phase diagram schematically shown in Fig. 6b. The CTS phase lies at the center of the pentagon, while the five corners represent the points at which CTS is in equilibrium with stable metal, binary or ternary phases. The point defects investigated were Cu (VCu), Sn (VSn), and S (VS) vacancies as well as Cu on Sn (CuSn) and Sn on Cu (SnCu) antisites. The monoclinic structure consists of two different types of Cu sites and three S sites, thus defect formation energies of all of the various sites are calculated.
Fig. 6c shows the formation energies of the various defect sites, revealing that VCu has the lowest formation energy. Consequently, VCu is considered as the main acceptor state of CTS responsible for its p-type conductivity. On the other hand, the key donor state consists of a SnCu defect, exhibiting the second lowest formation energy. VSn and CuSn states appear thermodynamically unfavorable based on these calculations. In view of these findings, the deep state associated with EA,D2 = 206 meV can be attributed to SnCu antisite defects. Interestingly, recent TEM studies with atomic-scale resolution have shown clear evidence of the presence of SnCu antisite domain boundaries in Cu2ZnSnS4 kesterite nanoparticles.51 Consequently, optimization of the CTS preparation is required in order to simultaneously control the doping density via VCu and suppress the SnCu bulk recombination sites. Solution based precursor methods as described in this work are uniquely suited to explore these conditions by: (i) controlling the composition ratio of the elemental precursors and (ii) introduction of dopants which can decrease atomic disorder.30
Footnotes |
† Electronic supplementary information (ESI) available: This includes film characterization details, thermogravimetric analysis of the precursor, schematic of the refined crystal structure, typical impedance spectra of devices with trial circuits and temperature dependence of R1, C1, R2 and C2. See DOI: 10.1039/c7se00150a |
‡ Further details of the crystal structure investigation(s) may be obtained from the Fachinformationszentrum Karlsruhe, 76344 Eggenstein-Leopoldshafen (Germany), on quoting the depository number CSD-432354. |
This journal is © The Royal Society of Chemistry 2017 |