Impact Factor 3.994

Frontiers reaches 6.4 on Journal Impact Factors

This article is part of the Research Topic

Solar Energy Harvesting

Original Research ARTICLE

Front. Chem., 26 January 2018 |

Effect of Magnesium Incorporation on Solution-Processed Kesterite Solar Cells

Raquel Caballero1*, Stefan G. Haass2, Christian Andres2, Laia Arques3, Florian Oliva3, Victor Izquierdo-Roca3 and Yaroslav E. Romanyuk2
  • 1Departamento de Física Aplicada, Universidad Autónoma de Madrid, Madrid, Spain
  • 2Laboratory for Thin Films and Photovoltaics, Empa- Swiss Federal Laboratories for Materials Science and Technology, Dübendorf, Switzerland
  • 3Catalonia Institute for Energy Research (IREC), Barcelona, Spain

The introduction of the alkaline-earth element Magnesium (Mg) into Cu2ZnSn(S,Se)4 (CTZSSe) is explored in view of potential photovoltaic applications. Cu2Zn1−xMgxSn(S,Se)4 absorber layers with variable Mg content x = 0…1 are deposited using the solution approach with dimethyl sulfoxide solvent followed by annealing in selenium atmosphere. For heavy Mg alloying with x = 0.55…1 the phase separation into Cu2SnSe3, MgSe2, MgSe and SnSe2 occurs in agreement with literature predictions. A lower Mg content of x = 0.04 results in the kesterite phase as confirmed by XRD and Raman spectroscopy. A photoluminescence maximum is red-shifted by 0.02 eV as compared to the band-gap and a carrier concentration NCV of 1 × 1016 cm−3 is measured for a Mg-containing kesterite solar cell device. Raman spectroscopy indicates that structural defects can be reduced in Mg-containing absorbers as compared to the Mg-free reference samples, however the best device efficiency of 7.2% for a Mg-containing cell measured in this study is lower than those frequently reported for the conventional Na doping.


Kesterite-type material Cu2ZnSn(S,Se)4 (CZTSSe) has been recognized as a promising candidate for low-cost thin-film solar cells due to its large absorption coefficient, tunable band-gap Eg between 1.0 and 1.5 eV adjusted via S/Se-ratio, low toxicity and earth-abundant nature. This technology has achieved a 12.6% maximum performance (Wang C. et al., 2014), still far away from that of 22.6% for Cu(In,Ga)Se2 solar cells (Jackson et al., 2016). The main performance limitation of kesterite-based solar cells is the open circuit voltage deficit (Eg/q − Voc). One of the reasons is the non-optimal quality of kesterite absorber and the presence of secondary phases (Siebentritt and Schorr, 2012). Another reason can be the unfavorable alignment of the conduction band minimum (CBM) at the CZTSSe/CdS interface (Platzer-Björkman et al., 2015). Gokmen et al. (2013) pointed out that the [CuZn + ZnCu] defect cluster could be the origin of electrostatic potential fluctuation in the CZTSSe absorber layer. That fluctuation is probably a significant factor that decreases the photovoltaic (PV) device performance. It was suggested that the substitution of Cu or Zn by other elements as Magnesium (Mg) could suppress the antisite defects CuZn and/or ZnCu formation that limit kesterite solar cells efficiency (Zhong et al., 2016).

Several studies about the Mg incorporation into kesterite have recently been reported, but the observed effects of Mg are contradictory. CIGSe bulk material doped with Mg was deposited by liquid phase sintering method measuring a decreased hole concentration with the increase in Mg content (10 at %, Mg/(In+Mg) = 0.1), which was attributed to the MgCu donor defect (Monsefi and Kuo, 2014). At the same time, an increased hole concentration was observed for 5 at % Mg explained by the substitution of In with Mg ion to form MgIn acceptor defect. Cu2MgSnS4 thin films grown by ultrasonic co-spray pyrolysis showed p-type conductivity and band-gap energy of 1.76 eV (Guo et al., 2016). In contrary, n-type conductivity was estimated for (Cu2xMgx)ZnSnSe4 bulk materials with x = 0.1−0.4, which was attributed to the formation of the donor-type MgCu antisite defects (Kuo and Wubet, 2014). The formation of stable Cu2MgSn(S,Se)4 was calculated based on density functional theory (Zhong et al., 2016), whereas a complete phase separation was predicted by Wang W. et al. (2014).

The purpose of this work is to study the effect of Mg addition in various concentrations to CZTSSe solar cell absorbers. Mg is incorporated to the absorber thin films by adding a magnesium salt to the precursor solution. Two experimental series are presented. The first experiment involves the replacement of Zn for Mg concentrations in order to evaluate possible alloying effects on the absorber band-gap for Cu2Zn1−xMgxSn(S,Se)4 thin films, however a complete phase separation is observed for x = 0.55 and 1. For a lower Mg content of x = 0.04 the kesterite phase is obtained, and respective absorbers and solar cells are studied to reveal any structural and electronic effects of Mg-containing sample as compared to a nominally Mg-free one.

Materials and Methods

Absorber Preparation

The Cu2Zn1−xMgxSn(S,Se)4 thin films absorbers were deposited from the precursor solution with dimethyl sulfoxide (DMSO) as the solvent onto Mo/SiOx/soda-lime glass (SLG) substrates with a subsequent selenization using the methodology described in the previous work (Haass et al., 2015). 1 mm-thick SLG was cleaned in three different supersonic baths at a temperature of 80°C. The first bath consisted of salt-free water with soap (Borer Deconex), the second bath consisted of a weak acetic solution (~5%) and the final bath contained only de-ionized water (18 MΩ·cm). The residual water was blown off with nitrogen and the substrates were dried in vacuum prior to the subsequent layer deposition. The SiOx barrier layer was deposited at 200°C substrate temperature on SLG substrate by sputtering of a Si target in an Ar/O2 atmosphere. The layer thickness of SiOx is approximately 200-300 nm. The SiOx barrier layer was used to reduce the alkali elements out-diffusion from the SLG substrate although the presence of Na could not be eliminated since it can also be transported via the gas phase during the annealing step (Abzieher et al., 2016). The CZTSSe precursor solution contained thiourea (99%+, Sigma-Aldrich), SnCl2·2H2O (98%, Sigma-Aldrich), ZnCl2 (99.99%, Alfa Aesar) and CuCl2 (98%+, Alfa Aesar) dissolved in DMSO (99.9%, Alfa Aesar). Mg was introduced by adding Mg(CH3COO)2·2H2O to the precursor solution. Table 1 shows the nominal Mg composition and metal ratios in the precursor solutions.


Table 1. Mg content and nominal metal ratios in the precursor solution (sol) and selenized absorbers (abs) for all the samples.

The solution was spin-coated onto a Mo/SiOx-coated and dried on a hotplate at 320°C in air. The spin-coating and drying steps were repeated 12 times in order to obtain the desired precursor film thickness of approximately 1.5–2 μm (Haass et al., 2015). Samples were annealed in a rapid thermal annealing (RTP) reactor (AS-ONE 150 from Annealsys) inside a silica-coated graphite box with selenium pellets (0.8 g). A three- stage temperature gradient was employed for annealing samples with holding temperatures at 300, 500, and 550°C for 30, 45, and 5 min, respectively. After selenization the absorbers were immersed into a 10 wt% KCN solution for 30 s in order to remove any copper-rich phases and clean the surface from contaminations and oxides (Haass et al., 2015).

Device Preparation

A 50–70 nm thick CdS buffer layer was deposited by chemical bath deposition, and 70 nm/250 nm i-ZnO/Al:ZnO bilayer was sputtered. A Ni/Al top grid and an antireflection coating of MgF2 were deposited by e-beam evaporation. Individual solar cells were mechanically scribed to an area of 0.30 ± 0.02 cm2 (Haass et al., 2015).

Materials Characterization

Inductively coupled plasma mass spectrometry (ICP-MS) measurements were carried out to determine the Mg content in the selenized absorber layers, which is listed in Table 1. The ICP-MS measurements were done on the full cells (without antireflection coating) for samples A and B with prior etching of the window layer by 5% acetic acid for 60 s and rinsing in distilled H2O. For the sample C with much lower Mg concentration we used the ICP-MS results of the absorber layer on SLG (no Mo) without any etching. The Mg content measured by ICP-MS for sample C is 0.04, which is lower than the nominal content of 0.17, indicating that Mg is lost during the absorber fabrication, and the reason needs to be investigated. Secondary Ion Mass Spectrometry (SIMS) measurements were recorded on a TOF-SIMS system from ION-TOF using O2+ primary ions with 2 keV of ion energy, a current of 400 nA and a raster size of 400 × 400 μm2. An area of 100 × 100 μm2 in the case of depth profiles was analyzed using Bi+ ions with 25 keV of ion energy. Energy dispersive X-ray (EDX) analysis and X-ray fluorescence (XRF) were used to quantify the composition of matrix elements. The XRF measurements were done on selenized absorber layers without any etching. Scanning Electron Microscopy (SEM) and EDX measurements were done on a Hitachi S-4800 electron miscroscope. X-ray diffraction (XRD) patterns were recorded in 2θ\θ scan mode using a Bruker D8 diffractometer with CuKα radiation (λ = 1.5418 Å̀, beam voltage: 40 kV, beam current: 40 mA, calibrated using Si100 and Si111 single crystals), a step size of 0.04° and a scan rate of 0.5 s/step (Haass et al., 2015). Grazing incidence XRD (GIXRD) measurements were collected with a PaNAlytical X'Pert Pro MPD diffractometer using Cu Kα radiation (λ = 1.5418 Å̀, beam voltage: 40 kV, beam current: 40 mA), a step size of 0.02°, a scan rate of 2 s/step and a multilayer mirror (Caballero et al., 2015). Detector scans with incident angles of 1°, 3°, and 5° were carried out. Raman scattering measurements were performed in back scattering configuration using a highly sensitive Raman apparatus developed at IREC consisting in a Horiba Jobin Yvon iRH320 spectrometer coupled with a low noise CCD detector cooled at −70°C. In this system, excitation and light collection were made through a macro optic system with a laser spot diameter of about 70 μm. Back-scattering measurements were performed under 633 and 785 nm excitation wavelengths by focusing laser spot directly onto the layer surface which allows the assessment of the absorber layer without any contribution from the upper layers. In the case of back surface measurements, laser spot was directly focused on the back surface after mechanical lift-off process. Excitation power was kept below 26 W/cm2 in order to avoid presence of thermal effects in spectra. The first-order Raman spectrum of monocrystalline silicon (Si) was measured as a reference before and after each Raman spectrum acquisition, and spectra were corrected by imposing Si first order at 520 cm−1 (Oliva et al., 2017).

Device Characterization

The J-V characterization of solar cells was performed under standard test conditions (100 mWcm−2, 25°C, AM1.5G) using a solar simulator calibrated with a certified Si diode. External Quantum Efficiency (EQE) spectra were recorded using a chopped white light source (900 W halogen lamp) with a LOT MSH-300 monochromator, which was calibrated with Si and Ge photodiodes. The illuminated area on the sample was 0.1 cm2 including grid lines. Photoluminescence (PL) spectra were measured on a FT300 fluorescence lifetime spectrometer from PicoQuant with a 639 nm pulsed diode laser as excitation source (pulse width 90 ps, repetition rate 10 MHz) and a thermoelectric cooled Hamamatsu NIR-PMT module H10330A-45 (rise time 0.9 ns, transit time spread 0.4 ns). Capacitance-Voltage (C-V) room temperature measurements were carried out with a LCR-meter from Agilent (E4990A) with an AC-voltage of 30 mV (Haass et al., 2015). JV-T was carried out in the temperature range from 138 to 298 K. For temperature dependent current-voltage measurements the solar cell was placed on temperature controlled Cu stage inside an evacuated cryostat cooled with liquid nitrogen and illuminated by a 900 W halogen lamp. The sample temperature was measured by thermocouples and regulated by a PID controller. The intensity of the incident light was varied by 2 orders of magnitude from approximately 1 - 140 mWcm−2 using neutral density filters.

Results and Discussion

Alloying with High Concentration of Mg: Cu2Zn1−xMgxSn(S,Se)4

The first experiment series involved heavy alloying of Zn with Mg in Cu2Zn1−xMgxSn(S,Se)4 (CZMTSSe) thin films. Figure 1A shows XRD spectra of samples A (x = 1) and B (x = 0.5). One can observe a complete phase separation, identifying the reflexes corresponding to MgSe2, MgSe, SnSe2, and Cu2SnSe3 in both samples. The presence of CZTSSe and ZnSe phase is difficult to exclude due to the closeness of their diffraction peaks with those of Cu2SnSe3. Already during the deposition of the Cu2Zn1−xMgxSnS4 (CZMTS) precursor layer, it was observed that the solution was inhomogeneous with some precipitates.


Figure 1. (A) GIXRD patterns of samples A and B exhibiting reflexes of Mo, Cu2SnSe3, SnSe, SnSe2, MgSe, and MgSe2. The PDF files number 04-012-4693, 04-009-2277, 01-089-2939, 04-018-3911, and 04-004-2933 have been used to identify Cu2SnSe3, SnSe, SnSe2, MgSe, and MgSe2 phases respectively. The dashed lines correspond to Cu2SnS3 phase. (B,C) Surface morphology of the samples A and B respectively. Different grains are formed, which are characteristic of the different secondary phases identified by EDX.

Figures 1B,C illustrate SEM surface morphology of samples A and B. Different types of grains are observed. For sample A (x = 1), SnSe2 is observed on top of smaller grains with the presence of Cu (15.1 at %)-Mg (16.9 at %)-Sn (10.1%)-Se (44.7 at%)-S (13.2 at %) (see Figure 1B). For sample B with x = 0.5, larger grains corresponding to SnSe2 and the smallest grains consistent with MgSe2 and Cu2SnSe3 are detected as shown in Figure 1C.

The observed phase separation is in agreement with the predictions of Wang et al. who describes the instability of the Cu2MgSnS4 by the following reaction:

Cu2MgSnS4MgS+Cu2SnS3    (1)

The calculated energy change of this reaction is exothermic (ΔE = −0.01 eV < 0), meaning that the phase separation of Cu2MgSnS4 proceeds spontaneously, in accordance with the disappearance of the stable region in the chemical potential space. The instability for this compound is directly related to the fact that for group IIA element, Cu2MgSnS4 is more stable in the ionic rocksalt structure with S than in the more covalent tetrahedral environment as in kesterite. When an I2-II–IV–VI4 compound is unstable, the corresponding elements, in this case Mg, can only be incorporated into a stable I2-II–IV–VI4 with low concentration (Wang W. et al., 2014). No working solar cells could be obtained using samples A and B as a consequence of the phase separation.

Effect of Low Mg Concentration on Kesterite Solar Cells

The effect of the addition of low Mg concentration was studied by comparing sample C with a Mg content of x = 0.04 (measured by ICP-MS in Table 1) to the nominally Mg-free, sample D. The distribution of Mg inside sample C measured by SIMS is not homogeneous (Figure 2A), with a higher Mg signal in the lower part of the absorber. This can be correlated to the morphology of samples C and D visible in the cross-section. Both absorbers exhibit a bi-layer structure with large-grain material on top and smaller grains close to the back contact as reported for the solution-processed kesterite layers (Haass et al., 2015). We assume that a higher Mg content in the lower part of sample C is due to a larger number of grain boundaries which can accommodate more Mg in comparison to large grains of the top crust. It is worth mentioning that a significant Na signal can be detected by SIMS in both samples C and D despite the fact that no Na was intentionally added to the precursor solution and a SiOx barrier was employed. ICP-MS measurements were carried out to quantify the Na concentration incorporated into the absorber layer. A Mg-free reference sample had a sodium concentration of only 70 ppm, while samples without SiOx barrier layer were reported with around 2,000 ppm (Sutter-Fella, 2014). This fact shows the effectiveness of the alkali barrier layer. This confirms that Na can also be transported via the gas phase during selenization (Abzieher et al., 2016).


Figure 2. (A) SIMS depth profile of CdS/CZTSSe/Mo corresponding to device C. (B) Cross-sectional SEM micrographs of the completed devices C and D. In the device C the positions 1 and 2 correspond to the penetration depths where Raman measurements were carried out (see Figure 3). A dash line indicates the border where a clear change of grain size and crystallinity of the bilayer structure occur.

Figure 3 displays GIXRD diffractograms of solar cells corresponding to samples C and D. Reflexes at 15.7°, 17.5° and 22.2° confirm unambiguously the CZTSSe phase. Moreover, reflexes at 31.8° and 56.6° appear, indicating the presence of the Mo(S,Se)2 layer formed during the selenization. No other secondary phases can be detected, although the presence of Zn(S,Se) and Cu2Sn(S,Se)3 phases cannot be ruled out since their diffraction reflexes coincide with those of CZTSSe. Reflexes of ZnO are also visible due to the window layer of i-ZnO/AZO in complete devices. A zoom of the 112 Bragg reflex of the kesterite phase measured at different grazing incidence angles does not show significant differences between both samples at different depths.


Figure 3. GIXRD patterns of samples C and D exhibiting reflexes of CZTSSe, Mo and Mo(S,Se)2. A zoom of the 112 Bragg peak measured at different grazing incidence angles do not show any significant difference between both samples. The PDF files number 04-010-6295 corresponding to CZTSe (solid line) and number 01-075-4122 assigned to CZTS (dash line) are also plotted.

Raman spectra measurements of the completed devices corresponding to samples C and D are shown in Figure 4 for two wavelengths of 633 and 785 nm. At these excitation wavelengths the contribution from the ZnO and CdS layers do not interfere with the absorber thin film signal allowing the direct characterization of the surface (penetration of < 70 nm for the 633 nm) and subsurface regions (< 150 nm for the 785 nm) (Dimitrievska et al., 2016b; Oliva et al., 2016). In order to access CZTSSe bulk, additional Raman measurements under 633 nm excitation were performed at the back of the Mg-containing CZTSSe after its mechanical delamination (lift-off process), corresponding to interface regions 1 and 2 (shown as red dotted line in Figure 2B). All Raman spectra exhibit characteristic features of highly Se-rich CZTSSe solid solution ([S]/([S]+[Se] ≈ 3% evaluated by Raman spectroscopy) yet. A detailed analysis of the Raman spectra of samples C and D can bring more insights about the impact of low Mg concentration on the CZTSSe absorber. Raman related parameters are summarized in Tables S1–S3 (See Supporting Information). The analysis of the Raman shift (RS) and the full-width-at-half-maximum (FWHM) of the Se-Se A' mode of kesterite phase (peak located at 196 cm−1) from Mg-containing/Mg-free samples at surface/sub-surface regions show that the differences are within the experimental error. The apparent FWHM increases between the surface and the sub-surface measurements (around 2 cm−1) (see Tables S1, S2 in Supporting Information) is attributed to the activation by Raman pre-resonant process of weak contributions overlapped with the 196 cm−1 peak. The absence of significant variations of the characteristic properties of the Se-Se peak suggests that both samples (Mg-containing/Mg-free) present a similar crystalline quality at the surface and sub-surface regions.


Figure 4. Raman spectra of the samples C and D acquired under 633 and 785 nm wavelengths (corresponding to different penetration depths, thus surface and sub-surface regions) and from the bulk of the absorber layer after lift-off process. (1) And (2) sides are the Raman spectra obtained at the regions (1) and (2) indicated in Figure 2B. The region centered at 330 cm−1 has been magnified x3 to clarify the contribution S-S vibration that indicates the S presence in the kesterite material.

For the Mg-containing sample C an extended analysis has been performed to compare the surface/bulk and bulk regions. The results show a red-shift of the main peak as well as an increase of the FWHM. These changes are an indication of crystal quality degradation, in agreement with the smaller grains at the back of the absorber layer observed by cross-sectional SEM images (see Figure 2B). Additionally, this red-shift of the main peak could be compatible with the replacement of Cu or Zn by a lighter atom as Mg. As shown in Figure 2A, a minimum Cu SIMS-signal coincides with a maximum Mg SIMS-signal and is located near the border between the big and small grains shown in Figure 2B. These results seem to suggest that a replacement of Cu by Mg could take place in that region, which could reduce the CuZn antisite defects formation.

The analysis of the changes at 175 and 242 cm−1 regions is interesting in order to evaluate the concentration variation of the defect clusters [ZnCu+VCu] and [2ZnCu+ZnSn] (Dimitrievska et al., 2015). A significant variation of the 175 cm−1 peak intensity between the samples C and D at the surface and the sub-surface is observed. The 175 and 196 cm−1 peak intensity ratio has been correlated with the concentration of the [ZnCu+VCu] (A-type) defect cluster and with the Voc of the device (Dimitrievska et al., 2015). For the Mg-containing samples an increase of this ratio is observed. This suggests that with the inclusion of Mg, a reduction of Vcu and/or ZnCu antisite formation is promoted. It is in agreement with the hypothesis that for low Mg concentration in CZTSSe samples, Mg replaces Cu and/or the Mg allows the inhibition of ZnCu antisites formation during CZTSSe synthesis (Zhong et al., 2016). Additionally the analyses of the bulk region show similar values of the 175 cm−1 peak intensity than those observed for the surface indicating a similar defect property at the surface and bulk of the absorber layer. Furthermore, the variation of the PV parameters of both samples with the A[175 cm−1]/([A175 cm−1)+A[196 cm−1]) ratio (see Figure 5) shows a clear correlation between the improvement of Voc, short circuit current density Jsc, fill factor FF and efficiency η with the increase of the [ZnCu+VCu] defect cluster concentration (reduction of the A[175 cm−1]/([A175 cm−1)+A[196 cm−1]) ratio), similarly as what has already been reported in previous works for pure CZTSe (Dimitrievska et al., 2016a,b). Therefore, optoelectronic parameters increase with the increase of the defect cluster concentrations as it is shown in Figure 5 for both samples, Mg-containing and Mg-free. This fact suggests that Mg improves material properties allowing better device performances, but reduces the [ZnCu+VCu] defect clusters density.


Figure 5. Correlation of the photovoltaic parameters of devices C and D with the A[175 cm−1]/(A[175 cm−1] + A[196 cm−1]) Raman intensity ratio.

The evolution of the peak intensity at 242 cm−1 shows that for the Mg-containing and Mg-free samples, the surface and sub-surface regions are similar while a clear increase is observed for the spectra acquired from the points (1) and (2) of the lift-off for sample C, regions with smaller grain size. This suggests that for this region, the [2ZnCu+ZnSn] defect cluster is promoted inducing the degradation of the crystal quality observed by changes of the FWHM and Raman shift (Dimitrievska et al., 2015).

An average efficiency of 5.8%, Voc = 405 mV, Jsc = 36.3 mA cm−2 and FF = 39.2% was measured for nine Mg-containing solar cells on sample C, whereas the nominally Mg-free sample D yielded an average efficiency of 4.2%, Voc = 376 mV, Jsc = 29.0 mA cm−2 and FF = 38.2% for nine cells. J-V curves and EQE spectra for best cells from sample C and D are shown in Figure 6. It appears that the performance of Mg-containing cells is higher than that for nominally Mg-free devices, in particular by improving the charge collection in the long-wavelength part of the spectrum (Figure 6B). However, no reduction of the Voc deficit is measured after adding Mg, being of 0.59 V for both best cells. The band gap energy Eg determined from the inflection point of the EQE spectrum is 1.01 eV and 0.97 eV for samples C and D, respectively. The comparison of reverse-biased EQE measurement with the EQE without bias is also plotted in Figure 6B for the device C. Although the Mg-containing device exhibits an enhanced carrier collection, it can be further improved under negative bias (−1 V). This bias dependence is signature of a poor carrier collection toward the back of the absorber layer (Scheer and Schock, 2011), which can be related to the smaller grain size next to the back contact.


Figure 6. (A) J-V characteristic of the best devices C and D (dark curve—dash line, light curve—solid line). Photovoltaic parameters are also shown and refer to total area measurements. (B) External quantum efficiency of devices C and D. EQE spectrum for device C under a reversed bias of −1 V is also displayed as well as the ratio between voltage biased (−1 V) and unbiased EQE.

Figure 7A shows PL spectra acquired at room temperature of devices C and D. The PL spectra exhibit a broad peak for both samples. The broad shape indicates the presence of tail states and/or potential fluctuations (Gokmen et al., 2013). Red-shifted PL peak maxima compared to Eg, ΔEEg−PL, by 20 meV and 25 meV are determined for solar cells C and D, respectively. These values are in the range obtained for CZTSe grown by vacuum-based methods (20 meV) (Lee et al., 2015).


Figure 7. (A) PL spectra for samples C and D. The PL maxima occur at a smaller energy value than the band gap (as determined from the EQE curve) for both samples. (B) Temperature dependent J-V measurement (dark curve—dotted line, light curve—solid line) of device C. The inset shows a linear fit of Voc that can be extrapolated to an intersection value EA = 0.95 eV, which is close to the estimated Eg = 1.01 eV, as determined by EQE. (C) NCV vs. x for the devices C and D extracted from C-V measurements. The values of depletion width xd and NCV at 0 V are shown.

Figure 7B displays the JV-T measurement of the device C. No roll-over effect is observed even at the lowest temperature of 138 K. The temperature dependence of the Voc extrapolated to T = 0 K provides an activation energy EA = 0.95 eV, which can correspond to the main recombination channel. Since this EA value is very near the Eg of 1.01 eV determined by EQE, it appears that the dominant recombination path is located within the bulk of the kesterite absorber layer rather than at the CdS/CZTSSe interface.

Furthermore, capacitance-voltage (C-V) measurement was carried out for both solar cells (see Figure 7C). The depletion width (xd) of the Mg-containing CZTSSe device is approximately of 0.2 μm at a bias of 0 V. A carrier concentration NCV of 1 × 1016 cm−3 can be extracted for the Mg-containing absorber, whereas a lower carrier concentration NCV of 4 × 1015 cm−3 was measured for the Mg-free sample, indicating that the presence of Mg increases acceptor density. The doping level is in the same range as for Na-doped CZTSe solar cells fabricated using a NaF precursor layer and by diffusion from SLG substrates (xd = 0.23 μm and NCV = 1 × 1016 cm−3) (Lee et al., 2015).


The incorporation of Mg into Cu2Zn1−xMgxSn(S,Se)4 layers leads to the complete phase separation for x = 0.5…1, whereas a lower content of x = 0.04 can preserve the kesterite phase. As compared to the nominally Mg-free sample, low quantities of Mg can improve the grain growth, reduce the number of structural defects and increase the acceptor density in the solar cell absorber. In this respect Mg appears to behave similar to the conventional alkali dopants such as Na or K. The first Mg-containing kesterite cell has been fabricated with the highest conversion efficiency of 7.2%. This value is lower than the efficiency of 11–12% obtained for Li, K or Na-doped devices (Haass et al., 2017), indicating that the addition of Mg is less effective.

Author Contributions

RC, SH, and YR designed the research and experiments. RC and SH fabricated solar cells, characterized layers and solar cells. CA assisted with the analysis. LA and FO carried out Raman spectroscopy measurements. VI-R assisted with Raman spectra analysis. RC, SH, FO, and YR wrote the paper. All authors contributed with discussions.

Conflict of Interest Statement

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.


This work was supported by the Framework 7th program under the project KESTCELLS (FP7-PEOPLE-2012-ITN-316488), Spanish Ministry of Education, Culture and Sport within the José Castillejo program (CAS 15/00070) and MINECO project WINCOST (ENE2016-80788-C5-2-R). RC acknowledges financial support from Spanish MINECO within the Ramón y Cajal program (RYC-2011-08521). The authors would also like to thank the whole team of the Laboratory for Thin Films and Photovoltaics at Empa.

Supplementary Material

The Supplementary Material for this article can be found online at:


Abzieher, T., Schnabel, T., Hetterich, M., Powalla, M., and Ahlswede, E. (2016). Source and effects of sodium in solution-processed kesterite solar cells. Phys. Status Solidi A 213, 1039–1049. doi: 10.1002/pssa.201532619

CrossRef Full Text | Google Scholar

Caballero, R., Cano-Torres, J. M., Garcia-Llamas, E., Fontané, X., Pérez-Rodriguez, A., Greiner, D., et al. (2015). Towards the growth of Cu2ZnSn1-xGexS4 thin films by a single-stage process: effect of substrate temperature and composition. Sol. Energy Mater. Sol. Cells 139, 1–9. doi: 10.1016/j.solmat.2015.03.004

CrossRef Full Text | Google Scholar

Dimitrievska, M., Fairbrother, A., Saucedo, E., Pérez-Rodríguez, A., and Izquierdo-Roca, V. (2015). Influence of compositionally induced defects on the vibrational properties of device grade Cu2ZnSnSe4 absorbers for kesterite based solar cells. Appl. Phys. Lett. 106:073903. doi: 10.1063/1.4913262

CrossRef Full Text | Google Scholar

Dimitrievska, M., Fairbrother, A., Saucedo, E., Pérez-Rodríguez, A., and Izquierdo-Roca, V. (2016a). Secondary phase and Cu substitutional defect dynamics in kesterite solar cells: impact on optoelectronic properties. Sol. Energy Mater. Sol. Cells 149, 304–309. doi: 10.1016/j.solmat.2016.01.029

CrossRef Full Text | Google Scholar

Dimitrievska, M., Giraldo, S., Pistor, P., Saucedo, E., Pérez-Rodríguez, A., and Izquierdo-Roca, V. (2016b). Raman scattering analysis of the surface chemistry of kesterites: impact of post-deposition annealing and Cu/Zn reordering on solar cell performance. Sol. Energy Mater. Sol. Cells 157, 462–467. doi: 10.1016/j.solmat.2016.07.009

CrossRef Full Text | Google Scholar

Gokmen, T., Gunawan, O., Todorov, T. K., and Mitzi, D. B. (2013). Band tailing and efficiency limitation in kesterite solar cells. Appl. Phys. Lett. 103:103506. doi: 10.1063/1.4820250

CrossRef Full Text | Google Scholar

Guo, Y., Cheng, W., Jiang, J., Zuo, S., Shi, F., and Chu, J. (2016). The structural, morphological and optical-electrical characteristic of Cu2XSnS4 (X: Cu,Mg) thin films fabricated by novel ultrasonic co-spray pyrolysis. Mater. Lett. 172, 68–71. doi: 10.1016/j.matlet.2016.02.088

CrossRef Full Text | Google Scholar

Haass, S. G., Diethelm, M., Werner, M., Bissig, B., Romanyuk, Y. E., and Tiwari, A. N. (2015). 11.2% efficient solution processed kesterite solar cell with a low voltage deficit. Adv. Energy Mater. 5:1500712. doi: 10.1002/aenm.201500712

CrossRef Full Text | Google Scholar

Haass, S. G., Andres, C., Figi, R., Schreiner, C., Bürki, M., Romanyuk, Y. E., et al. (2017). Complex interplay between absorber composition and alkali doping in high-efficiency kesterite solar cells. Adv. Energy Mater. doi: 10.1002/aenm.201701760. [Epub ahead of print].

CrossRef Full Text | Google Scholar

Jackson, P., Wuerz, R., Hariskos, D., Lotter, E., Witte, W., and Powalla, M. (2016). Effects of heavy alkali elements in Cu(In,Ga)Se2 solar cells with efficiencies up to 22.6%. Phys. Status Solid RRL 10, 583–586. doi: 10.1002/pssr.201600199

CrossRef Full Text | Google Scholar

Kuo, D. H., and Wubet, W. (2014). Mg dopant in Cu2ZnSnSe4: an n-type former and a promoter of electrical mobility up to 120 cm2V−1s−1. J. Solid Stat. Chem. 215, 122–127. doi: 10.1016/j.jssc.2014.03.034

CrossRef Full Text | Google Scholar

Lee, Y. S., Gershon, T., Gunawan, O., Todorov, T. K., Gokmen, T., Virgus, T., et al. (2015). Cu2ZnSnSe4 thin-film solar ells by thermal co-evaporation with 11.6% efficiency and improved minority carrier diffusion length. Adv. Energy Mater. 5:1401372. doi: 10.1002/aenm.201401372

CrossRef Full Text | Google Scholar

Monsefi, M., and Kuo, D. H. (2014). Influence of Mg doping on electrical properties of Cu(In,Ga)Se2 bulk materials. J. Alloys Comp. 582, 547–551. doi: 10.1016/j.jallcom.2013.08.101

CrossRef Full Text | Google Scholar

Oliva, F., Arques, L., Acebo, L., Guc, M., Sánchez, Y., Alcobé, X., et al. (2017). Characterization of Cu2SnS3 polymorphism and its impact on optoelectronic properties. J. Mater. Chem. A 5, 23863–23871. doi: 10.1039/C7TA08705E

CrossRef Full Text | Google Scholar

Oliva, F., Kretzschmar, S., Colombara, D., Tombolato, S., Ruiz, C. M., Redinger, A., et al. (2016). Optical methodology for process monitoring of chalcopyrite photovoltaic technologies: application to low cost Cu(In,Ga)(S,Se)2 electrodeposition based processes. Sol. Energy Mater. Sol. Cells 158, 168–183. doi: 10.1016/j.solmat.2015.12.036

CrossRef Full Text | Google Scholar

Platzer-Björkman, C., Frisk, C., Larsen, J. K., Ericson, T., Li, S. Y., Scragg, J. J. S., et al. (2015). Reduced interface recombination in Cu2ZnSnS4 solar cells with atomic layer deposition Zn1–xSnxOy buffer layers. Appl. Phys. Lett. 107:243904. doi: 10.1063/1.4937998

CrossRef Full Text | Google Scholar

Scheer, R., and Schock, H. W. (2011). Chalcogenide Photovoltaics. Weinheim: Wiley-VCH Verlag.

Google Scholar

Siebentritt, S., and Schorr, S. (2012). Kesterites- a challenging material for solar cells. Prog. Photovoltaics Res. Appl. 20, 512–519. doi: 10.1002/pip.2156

CrossRef Full Text | Google Scholar

Sutter-Fella, C. M. (2014). Solution-Processed Kesterite Absorbers for Thin Film Solar Cells. Ph.D. thesis, ETH Zurich.

Wang, C., Chen, S., Yang, J. H., Lang, L., Xiang, H. J., Gong, X. G., et al. (2014). Design of I2-II-IV-VI4 semiconductors through element substitution: the thermodynamic stability limit and chemical trend. Chem. Mater. 26, 3411–3417. doi: 10.1021/cm500598x

CrossRef Full Text | Google Scholar

Wang, W., Winkler, M. T., Gunawan, O., Gokmen, T., Todorov, T., Zhu, Y., et al. (2014). Device characteristics of CZTSSe thin-film solar cells with 12.6% efficiency. Adv. Energy Mater. 4:13011465. doi: 10.1002/aenm.201301465

CrossRef Full Text | Google Scholar

Zhong, G., Tse, K., Zhang, Y., Li, X., Huang, L., Yang, C., et al. (2016). Induced effects by the substitution of Zn in Cu2ZnSnX4 (X = S and Se). Thin Solid Films 603, 224–229. doi: 10.1016/j.tsf.2016.02.005

CrossRef Full Text | Google Scholar

Keywords: thin film solar cells, kesterite, Mg, solution processing, structural defects

Citation: Caballero R, Haass SG, Andres C, Arques L, Oliva F, Izquierdo-Roca V and Romanyuk YE (2018) Effect of Magnesium Incorporation on Solution-Processed Kesterite Solar Cells. Front. Chem. 6:5. doi: 10.3389/fchem.2018.00005

Received: 11 December 2017; Accepted: 10 January 2018;
Published: 26 January 2018.

Edited by:

Alberto Mittiga, National Agency for New Technologies, Energy and Sustainable Economic Development, Italy

Reviewed by:

Maarja Grossberg, Tallinn University of Technology, Estonia
Simona Binetti, Università degli Studi di Milano Bicocca, Italy

Copyright © 2018 Caballero, Haass, Andres, Arques, Oliva, Izquierdo-Roca and Romanyuk. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.

*Correspondence: Raquel Caballero,