Conduction-band effective mass and bandgap of ZnSnN2 earth-abundant solar absorber

Pseudo III-V nitride ZnSnN2 is an earth-abundant semiconductor with a high optical absorption coefficient in the solar spectrum. Its bandgap can be tuned by controlling the cation sublattice disorder. Thus, it is a potential candidate for photovoltaic absorber materials. However, its important basic properties such as the intrinsic bandgap and effective mass have not yet been quantitatively determined. This paper presents a detailed optical absorption analysis of disordered ZnSnN2 degenerately doped with oxygen (ZnSnN2−xOx) in the ultraviolet to infrared region to determine the conduction-band effective mass (mc*) and intrinsic bandgap (Eg). ZnSnN2−xOx epilayers are n-type degenerate semiconductors, which exhibit clear free-electron absorption in the infrared region. By analysing the free-electron absorption using the Drude model, mc* was determined to be (0.37 ± 0.05)m0 (m0 denotes the free electron mass). The fundamental absorption edge in the visible to ultraviolet region shows a blue shift with increasing electron density. The analysis of the blue shift in the framework of the Burstein-Moss effect gives the Eg value of 0.94 ± 0.02 eV. We believe that the findings of this study will provide important information to establish this material as a photovoltaic absorber.

important for establishing ZnSnN 2 as a photovoltaic absorber. Insufficient knowledge of the m c * and E g values make it difficult to use this material for photovoltaic absorbers. Although optically determined E g values have been reported by some research groups, the values are scattered in the range of 1.0-2.0 eV (Table 1) [9][10][11][12] . We hypothesize that the reported values are not incorrect and their scattering originates from the blue shift of E g due to the conduction-band filling effect (the so-called Burstein-Moss shift), because the conduction-electron densities were largely different in each study. As will be described later, this hypothesis was justified, and thus the intrinsic E g value was successfully determined from the analysis in the framework of the Burstein-Moss effect. As for the reported m c * values, a large discrepancy between the theoretical and experimental values exists, which are summarized in Table 1. The theoretical study predicted small m c * values of 0.12m 0 (m 0 denotes the free-electron mass) 9 , whereas the experimental study yielded a value four times larger at 0.5m 0 . This large discrepancy motivated us to discover which value is correct. The experimental m c * value was derived from the analysis of the Burstein-Moss shift 11 . Strictly speaking, the effective mass determined by this method is not m c * but the 'reduced' effective mass (m vc * ), as the m vc * is a function of m c * and the valence-band effective mass (m v * ), as discussed later. The m vc * value becomes equal to m c * value only when the m v * value is much higher than that of m c * (m c * /m v * ≈ 0). Thus, it is likely that the reported experimental value deviates somewhat from the true value. To obtain the correct m c * value, we performed an analysis of the free-electron absorption/reflection in the optical transmittance and reflectance spectra, because the effective mass obtained by this analysis purely corresponds to m c * . In fact, m c * in many semiconductors has determined by this method.
In addition, it has been reported that ZnSnN 2 thin films grown by physical vapor deposition were heavily doped with oxygen (ZnSnN 2−x O x ), even when oxygen was not intentionally introduced into the growth chamber 4,10,11 . Since the incorporated oxygen acts as an electron donor 5,13 , unintentionally oxygen-doped ZnSnN 2 films behave as n-type degenerate semiconductors with high electrical conductivity and free-electron absorption/ reflection in the infrared region. Fioretti et al. have addressed the unintentional oxygen-doping problem to obtain nondegenerate ZnSnN 2 14 . We can also make use of oxygen-doped ZnSnN 2 to determine the conduction-band effective mass by analysing the free-electron absorption/reflection.
In this study, epilayers of disordered ZnSnN 2−x O x with various x values were grown. Detailed analysis of their optical absorption/reflection in the ultraviolet to infrared region were performed to reveal the m c * and intrinsic E g values. As a result, the m c * and intrinsic E g values were found to be (0.37 ± 0.05)m 0 and 0.94 ± 0.02 eV, respectively. Herein, we describe the analysis in detail.

Results and Discussion
Film Growth and Structure. According to recent literature 15 , disordered ZnSnN 2 epitaxially grows on (111) yttria-stabilized zirconia (YSZ(111)). Thus, we used YSZ(111) single-crystal substrates for film growth. The structure of the ZnSnN 2 film strongly depended on substrate temperatures (T s ) and N 2 partial pressure (P N 2 ). Figure 2 depicts the phase diagram for ZnSnN 2 thin films, as a function of T s and P N 2 . As can be seen in this figure, epitaxial films were obtained in a relatively narrow region (closed circles in the shadowed region). Figure 3a shows the out-of-plane X-ray diffraction (XRD) pattern of ZnSnN 2 epilayers. Although the 111, 222, and 333 diffraction peaks are present in the substrate, only the 0002, 0004, and 0006 diffraction peaks from ZnSnN 2 were observed, indicating that the (0001) planes are parallel to the (111) planes of the substrate. The ZnSnN 2 epilayers were confirmed to be in the disordered-phase by in-plane ϕ-scans (see Methods section). Figure 3b shows a typical in-plane XRD pattern (scanned in ϕ-2θ χ mode) of ZnSnN 2 epitaxial layers. The diffraction pattern shows that the (1120) planes of the film are aligned with the (110) planes of the substrate. As shown in the inset in Fig. 3b, six evenly spaced 1120 peaks were observed in the ϕ-scan pattern, illustrating the 6-fold rotational symmetry of the basal plane in the disordered-phase (see Fig. 1b). Therefore, it can be confirmed that the disordered ZnSnN 2 phase was grown on the YSZ(111) substrate, with the epitaxial relationship of ZnSnN 2 (0001)//YSZ(111) and ZnSnN 2 (1120)//YSZ(110). Figure 4a shows the topographic atomic force microscopy (AFM) image of the ZnSnN 2 epitaxial film. The image shows that the epitaxial film had a homogeneous compact structure comprising of small grains with lateral diameters of ~300 nm. The surface possessed a rough topography as seen from the three-dimensional image shown in Fig. 4b. The root-mean-square roughness in the 2 × 2 μm 2 area was 1.6 nm. The rough surface is clearly  illustrated in the cross-sectional image (Fig. 4c). Such rough surface may stem from the resputtering effect 16 , which is due to the bombardment of high-energy particles from the plasma. However, the surface roughness did not essentially affect the electrical and optical analyses described below.
Film Composition and Electrical Properties. The cationic composition, hereafter referred to as Zn/ (Zn + Sn), was slightly Zn-rich in composition with Zn/(Zn + Sn) of ~0.55, which was almost independent of the growth condition (Fig. 5a). We detected oxygen in all the films, indicating that unintentional oxygen-doping occurred during the film growth (ZnSnN 2−x O x ). In contrast to Zn/(Zn + Sn), the oxygen concentration, x in ZnSnN 2−x O x , depended on P N 2 during film growth, as shown in Fig. 5a. The x value decreased with increasing P N 2 . The oxygen contamination likely originated from the residual vapor in the growth chamber. The molar ratio of water to nitrogen is expected to decrease as P N 2 increases or vice versa. This may be one of the reasons for the x-P N 2 dependency.
All the epilayers in the present study showed n-type conductivity. As seen from Fig. 5b, the electron density (n e ) also showed P N 2 dependency, as with x. The trend of the n e -P N 2 dependency was very similar to that of the x-P N 2 , suggesting that the incorporated oxygen impurities act as electron donors similar to nitrides such as InN 17,18 , GaN 19,20 , AlN 21,22 , and Zn 3 N 2 23,24 . Chen et al. theoretically predicted that oxygen impurities will occupy the nitrogen sublattice in ZnSnN 2 and behave as single electron donors 13 . Since the n e values were on the order of 10 20 cm −3 , the conduction electrons were highly degenerate (Fig. S1, Supplementary Information).
Effective Mass. Figure 6 shows the optical transmittance (T) and reflectance (R) spectra for the ZnSnN 2−x O x epilayers in the UV to IR region. In this figure, free-electron absorption/reflection can be seen. We performed a fitting analysis of the free-electron absorption/reflection using a dielectric function model to determine the conduction-band effective mass (m c * ) of disordered ZnSnN 2 (see Methods section). The square and circle marks in Fig. 6a-e represent the experimental data, and the solid lines denote the best fit spectra obtained by the double Tauc-Lorenz (TL) dispersion model coupled with the Drude model. In our combined double TL-Drude model, four fitting parameters for each TL term and two for the Drude term were considered, and described in detail in the Methods section. As can be seen in Fig. 6a-e, the model sufficiently reproduced the spectra across the entire spectral region for all the films. The values of the TL and Drude parameters obtained from the best-fit curves are summarized in Table S1 (Supplementary Information). In addition, we also show the dielectric functions obtained from the fitting analysis in Fig. S2. Using plasma energy (E p ) and n e values in Table S1, the conduction-band effective mass, m c * , can thus be derived (see Methods).    First principle studies predicted m c * values to be as small as ~0.1m 0 in both ordered-and disordered ZnSnN 2 4,9,25,26 . The predicted m c * values are comparable to those in other III-V and Zn-based nitride semiconductors including InN (0.055m 0 ) 27 , GaN (0.2m 0 ) 28 , GaAs (0.083m 0 ) 29 , and Zn 3 N 2 (0.08m 0 ) 30 . Although small m c * values generally leads to high electron mobility (μ), μ values in the ZnSnN 2 epilayers have been reported to be less than or equal to 10 cm 2 V −1 s −1 9,10 . Indeed, the μ values of the ZnSnN 2 epitaxial films in this study were as low as 10-20 cm 2 V −1 s −1 . A recent experimental study by Fioretti et al. suggested that the m c * value in disordered ZnSnN 2 is roughly 4-5 times larger than the predicted values (~0.5m 0 ) 11 . Our results were consistent with their value. Namely, we revealed that the m c * values predicted from first-principle calculations are probably underestimated and the m c * value is (0.37 ± 0.05)m 0 . Therefore, small electron mobility in ZnSnN 2 can be partly explained by the somewhat large m c * . However, the conduction-band effective mass alone cannot fully account for the low mobility in disordered ZnSnN 2 . Further investigation of the electron transport properties is still required to understand the reasons for the low electron mobility.
Meanwhile, enhancements in m c * with increasing n e have been reported for various degenerately doped semiconductors such as GaN 31 , InN 32 , Zn 3 N 2 33 , GaAs 34 , InSb 35 , ZnO 36,37 , TiO 2 38 , and CdSnO 4 39 . Such behaviours have been interpreted as a result of the nonparabolic conduction band. In contrast, the m c * of ZnSnN 2 seems to be almost independent of n e (Fig. 7), suggesting that the conduction band of disordered ZnSnN 2 has a parabolic shape.
Bandgap. The optical absorption coefficients (α) of the ZnSnN 2 epilayers were calculated using the relation- where d is the thickness of film. Figure 8 shows the absorption coefficient plotted on a log scale as a function of photon energy (hν) for four ZnSnN 2 epilayers with different n e values. The α spectra showed V-shaped dependency on photon energy. The absorption in the high hν region (hν > 1.5 eV) is due to the fundamental absorption, while the absorption in the low hν region (typically hν < 1.2 eV) corresponds to the free-electron absorption. We defined the optical gap energy (E g opt ) as the photon energy where the extrapolated straight lines of the free-electron absorption and fundamental absorption intersect (see Fig. 8).
The E g opt varied from 1.25 to 1.51 eV as n e increased from 2.8 × 10 20 to 5.1 × 10 20 cm −3 . Therefore, the increase in n e caused a blue shift of the E g opt . The blue shift has frequently been observed in degenerate ZnSnN 2 and is interpreted as the Burstein-Moss (BM) shift, which is a consequence of the conduction band filling effect due to heavy oxygen-doping 4,[9][10][11] . For degenerate n-type semiconductors, n e -dependent E g opt due to the BM effect is expressed as where E g is the intrinsic bandgap energy and m vc * is the reduced effective mass given by m vc * denotes the valence-band effective mass). In heavily doped semiconductors, bandgap narrowing due to electron-electron and electron-impurity interactions also occurs: the bandgap narrowing partially compensates the bandgap widening due to the BM shift [40][41][42] . We calculated the bandgap narrowing due to electron-electron and electron-impurity interactions (ΔE g e,e and ΔE g e,i , respectively) in a manner similar to the procedure described in literature 43 . Both the ΔE g e,e and ΔE g e,i were quite small in comparison with ΔE g BM , and those exhibited roughly linear dependence on n e 2/3 . Consequently, the total bandgap shift, ΔE g = ΔE g BM − ΔE g e,e −ΔE g e,i , was almost proportional to n e 2/3 (see Fig. S3, Supplementary Information). This result justifies plotting E g opt as a function of n e 2/3 to derive the intrinsic bandgap, E g . The E g opt values for all films are plotted (closed circles) together with the values reported in literature for disordered ZnSnN 2 films (squares 11 and triangles 9 ) as a function of n e 2/3 in Fig. 9.  Table S1. The E g opt values of these additional films are also presented in Fig. 9 (open circles). As seen from Fig. 9, E g opt shows linear dependency on n e 2/3 , which is consistent with the trend shown in Fig. S3. It is striking that the E g values in this study correlate with the values reported in literature. Hence, a linear regression fit (least-square approach) was applied to all the data including the values reported in literature (Fig. 9), indicating universality within the used data, the intercept of the straight-line gives the intrinsic bandgap (E g ). Accordingly, the E g value of disordered ZnSnN 2 was determined to be 0.94 ± 0.02 eV. The E g value reasonably agrees with the value predicted by the first-principle study for fully disordered ZnSnN 2 (E g = 0.98 eV) 9 .
In the present study, we did not take the cation composition into account for the determination of E g . The cation composition may be relevant for the E g value, though it is not clear whether the effect is significant or not. The cation composition dependency of E g is currently a research focus for ZnSnN 2 .  . The circles denote the data obtained in the present study (closed circles for films in Table S1 and open circles for additional films). Data obtained from Fioretti et al. 11 and Veal et al. 9 are denoted as squares and triangles, respectively.

Conclusion
In summary, the conduction-band effective mass and intrinsic bandgap of disordered ZnSnN 2 were optically determined to get an insight into the fundamental properties. Heavily oxygen-doped ZnSnN 2−x O x epitaxial layers were grown with various electron densities. These epilayers exhibited free-electron absorption/reflection in the infrared region. The detailed analysis of the free-electron absorption/reflection revealed that the conduction-band effective mass is (0.37 ± 0.05)m 0 . To the best of our knowledge, this study is the first to quantitatively determine m c * . Furthermore, the blue shift of the optical bandgap was analysed in the framework of the Burstein-Moss effect to determine the intrinsic bandgap, which has been in debate to date. Our optical bandgap values and those obtained from literature showed a universal 2/3-power dependency on the electron density. From this dependency, the intrinsic bandgap of disordered ZnSnN 2 was determined to be 0.94 ± 0.02 eV. We believe that the findings in the present study provide important information to establish this material as a photovoltaic absorber.

Methods
Thin Film Growth. ZnSnN 2 thin films were epitaxially grown on YSZ(111) single-crystalline substrates by reactive radio-frequency (RF) magnetron sputtering using a Zn 0.5 Sn 0.5 alloy target (diameter of 10 cm and purity of 3 N), at T s ranging from 250 to 350 °C. The base pressure of ~2 × 10 −4 Pa was established prior to the filmgrowth. An RF power of 70 W was applied to the target. A mixture of Ar and N 2 gas with various N 2 /(N 2 + Ar) ≡ f(N 2 ) ratios was introduced into the chamber through two independent mass flow controllers with a total flow rate of 5 sccm. The working pressure in the chamber (P w ) was held at 2.0 Pa during film growth. Films were grown under P N 2 ranging from 1.0 to 2.0 Pa. The nitrogen partial pressure was defined as P N 2 = f(N 2 ) × P w . The growth time was adjusted to obtain films with thicknesses of 100-300 nm.

Characterization Methods.
A Rigaku ATX-G X-ray diffractometer with Cu Kα radiation was employed to perform out-of-plane (θ-2θ) and in-plane (ϕ-2θ χ ) scans to evaluate the structure of the films. To ensure that grown ZnSnN 2 epilayers were in the disordered-phase, we performed in-plane ϕ-scans. Provided that the ZnSnN 2 films are in ordered-phase with the orthorhombic symmetry, a characteristic 110 peak will be observed at an in-plane diffraction angle (2θ χ ) of 20.07° in the ϕ-scans (Fig. S4a in Supplementary Information). We carried out the in-plane ϕ-scans at the fixed 2θ χ angle of 20.07° several times, but no peak was observed (Fig. S4b in Supplementary Information). Even when the 2θ χ angle was slightly increased or decreased, the 110 peak did not appear. These results indicate that the ZnSnN 2 films epitaxially grown on YSZ(111) were in the disordered-phase, which is consistent with the results reported in literature 15 . The surface morphology was characterized with an atomic force microscope (Hitachi AFM5100N). The compositions (Zn/(Zn + Sn) and x) of the ZnSnN 2 films were examined by X-ray photoelectron spectroscopy (PHI Versa Probe), using monochromated Al Kα (hν = 1486.6 eV) radiation. The relative sensitivity factor (RSF) approach was exploited to determine the compositions. It was confirmed that the compositions determined by the RSF method were consistent with those determined by Rutherford backscattering spectrometry 33,44 . Hence, it was concluded that the compositions in this study are sufficiently reliable. The description of the RSF approach are detailed in the Supplementary Information. Electrical properties were determined by Hall-effect measurements in the van der Pauw configuration (Toyo Corp. Resitest 8200). Optical transmittance and reflectance were collected between 0.3 and 5.0 µm using a UV-Vis-NIR spectrophotometer (Shimadzu UV-3150) and FTIR spectrometer (Shimadzu IRAffinity-1).

Fitting Analysis of T and R Spectra.
Owing to the large electron density on the order of 10 20 cm −3 , free-electron absorption and reflection are clearly seen in the infrared region (see Fig. 6a−e). The well-known Drude dielectric model was employed to describe the optical response by the free electrons in the IR region. The Drude function is given by, where E is the photon energy, Г D represents the broadening parameter, and E P denotes the plasma energy, which can be expressed as Here, ħ, e, and ε 0 , denote the reduced Planck constant, elemental charge, and the static dielectric constant of free space, respectively. As seen from Equation 3, E p 2 is theoretically proportional to n e when the dielectric response in the IR region can be described by the Drude model. In fact, E p 2 showed a linear dependency on n e (Fig. S5, Supplementary Information), suggesting appropriate use of the Drude function.
In addition to the Drude model, the Tauc-Lorentz (TL) dispersion model was considered to describe the optical response across the whole spectral region (see Fig. 6). The explicit expression of the TL model is given in literature 45 . The TL model has proved to be useful for the description of the fundamental absorption of many semiconductors 36,[45][46][47] . Recently, Deng et al. demonstrated that the dielectric response of ZnSnN 2 in UV to visible region can be reproduced by double TL functions 48 . Hence, we modelled the dielectric function of degenerate ZnSnN 2 as a sum of double TL (ε TL1 (ω) and ε TL2 (ω)) and Drude functions, whereby ε(ω) = ε TL1 (ω) + ε TL2 (ω) + ε D (ω). In this study, ε TL1 (ω) and ε TL2 (ω) were assigned to the fundamental absorption of the bandgap and the transition that corresponds to the M 4 −M 3 transition in group-III nitrides, respectively 49,50 . Theoretical T and R spectra calculated via the Fresnel formulas combined with the ε(ω) function were fitted to the experimental spectra. In the Scientific RepoRts | 7: 14987 | DOI:10.1038/s41598-017-14850-7 fitting procedure, a set of {A TL , Г TL , E 0 , E T } for each TL term and {E p , Γ p } for the Drude term were selected as fitting parameters. Here, A TL , Г TL , E 0 , and E T in the TL term represent the oscillator strength, broadening parameter, resonance energy, and Tauc gap energy, respectively 45 .