Mechanisms of charge transfer and electronic properties of Cu2ZnGeS4 from investigations of the high-field magnetotransport

Recent development of the thin film solar cells, based on quaternary compounds, has been focused on the Ge contain compounds and their solid solutions. However, for effective utilization of Cu2ZnGeS4, deeper investigations of its transport properties are required. In the present manuscript, we investigate resistivity, ρ (T), magnetoresistance and Hall effect in p-type Cu2ZnGeS4 single crystals in pulsed magnetic fields up to 20 T. The dependence of ρ (T) in zero magnetic field is described by the Mott type of the variable-range hopping (VRH) charge transfer mechanism within a broad temperature interval of ~100–200 K. Magnetoresistance contains the positive and negative components, which are interpreted by the common reasons of doped semiconductors. On the other hand, a joint analysis of the resistivity and magnetoresistance data has yielded series of important electronic parameters and permitted specification of the Cu2ZnGeS4 conductivity mechanisms outside the temperature intervals of the Mott VRH conduction. The Hall coefficient is negative, exhibiting an exponential dependence on temperature, which is quite close to that of ρ(T). This is typical of the Hall effect in the domain of the VRH charge transfer.

extensively for both types. On the other hand, the electronic transport of CZGeS, which is important for any material utilization in optoelectronic devices, is still investigated quite insufficiently. There exist only a few articles, where the room temperature resistivity 16,17,20 , and its behavior in the high temperature range of 450-550 K 22,29 have been observed, without any analysis of the data. Probably, the only exclusion is an explicit analysis of the resistivity and magnetoresistance in the kesterite type Cu 2 ZnSn x Ge 1-x Se 4 solid solutions 11 .
In the present work, the temperature dependence of the zero-field resistivity, ρ(T), and magnetotransport in CZGeS single crystals with wurtzstannite structure are investigated within a wide interval below 300 K. This permits identification of the conductivity mechanisms in different temperature ranges and determination of important electronic parameters, which are required to estimate efficiency of the material in optoelectronic devices. In particular, the joint analysis of the resistivity and magnetoresistance data has been widely demonstrated to be an effective method for such a purpose (see e. g. refs 30-38 and references therein), which has been already utilized recently in various kesterite type quaternary compounds 11,39 .

Results and Discussion
Experimental results. As can be seen in the top panel of Fig. 1, ρ(T) of all samples is similar, exhibiting a strongly activated behavior within the whole temperature range, available for the resistivity measurements.
As follows from Fig. 2, magnetoresistance (MR) contains both the positive (pMR) and the negative (nMR) contributions within the whole temperature range of ~50-300 K. In particular, the relative MR, Δρ/ρ ≡ [ρ(B) − ρ(0)]/ρ(0), is negative almost at any B and T, excluding only the cases of T = 70 K in #1, and T = 50 K and 300 K in #3 at B > 18 T. The behavior of MR with increasing temperature is noticeable, having different temperature dependences within two different temperature intervals, addressed to the lower and higher T, as well as various shapes of the curves in Fig. 2 in the corresponding intervals. Namely, between T = 70-100 K, 60-80 K and 50-70 K in #1, #2 and #3, respectively, the MR effect increases with increasing T. The decay of MR with T, conventional for semiconductors, starts only with further increase of the temperature beyond the intervals above. In addition, such a transition is accompanied by a gradual flattening of the MR shape, showing different relations between temperature dependences of the nMR and pMR contributions in each of the temperature intervals, defined above.
In Fig. 3 is evident, that the Hall resistivity, ρ H , of sample #1 is negative at any T between 70-300 K, despite of the positive sign of the thermopower and the p-type conduction of our CZGeS samples. Note, that the hole conduction is typical of all Cu based quaternary compounds of this family of materials, although possessing the kesterite structure (see e. g. refs 11, 39 and references therein). One cannot observe distinguishable deviations of the ρ H (B) function from linearity, which, however, may be connected partially to a rather high scattering of the data points in Fig. 3. On the other hand, a strong dependence of ρ H (T) between T = 70 and 230 K is evident at any field, weakening substantially with further increase of T up to 300 K. Temperature dependence of resistivity at B = 0. In the quaternary Cu based kesterite materials the hopping charge transfer, mainly the Mott variable-range hopping (VRH), has been reported for different chemical content, stoichiometry and macroscopic state, including single crystals 11,[39][40][41][42] , thin polycrystalline films [43][44][45][46][47][48] and powder samples 49,50 . The temperature border of the hopping conduction is usually very high, stretching even up to the room temperature already in single crystals, whereas width of the defect acceptor band, W, vary from ~10 to ~100 meV and even higher [39][40][41][42][43][44][45][46][47][48][49][50] . Both these features are quite uncommon for conventional semiconductors with shallow hydrogenic impurities, exhibiting usually W between ~1-10 meV and the hopping charge transfer at temperatures well below 300 K [32][33][34][35][36][51][52][53][54] . The point is that the energy of the main (most stable) acceptor level, E A , connected to Cu Zn antisites [55][56][57][58][59][60] , in such materials is rather high (e.g. in Cu 2 ZnSnS 4 and Cu 2 ZnSnSe 4 the values of E A ≈ 120-140 meV have been reported 39,[55][56][57][58], which favors the hopping charge transfer in general. On the other hand, domination of the Mott VRH conduction is connected with a high intrinsic lattice disorder, caused by formation of the disordered kesterite phase [61][62][63] . In turn, the disorder can be effectively varied experimentally, as has been demonstrated by observations of considerable changes in W of the Cu 2 ZnSnS 4 films after thermal annealing 44 . At this point, the VRH conduction is quite expectable in CZGeS, too, supported by E A = 140 meV 64 and by a complicated crystal structure [23][24][25][26][27][28] permitting similar reasons for the microscopic disorder, as in kesterites [61][62][63] . Therefore, it is reasonable to start the analysis of ρ(T) in our samples by searching the possible interval of the Mott VRH charge transfer, given by the resistivity law Here, D is the VRH prefactor constant and is the density of the localized states (DOS) at the Fermi level, μ, and a is the localization radius of charge carriers, scaling the space decay of the impurity wave functions [51][52][53][54] . Generally, the VRH conduction sets in, when it is energetically more profitable for a carrier to jump beyond nearest centers, with levels with large energy difference, towards more distant centers with a smaller energy difference with respect to the initial center 52 . Therefore, the factors favoring the VRH charge transfer are connected to decreasing temperature and increasing degree of the disorder. Formally, the VRH conduction of Eq. (1) is expected at ε m (T) < W, where ε m (T) ≈ 0.5 k B (T 3 T 0 ) 1/2 is the mean hopping energy 40,53 . In addition, Eqs (1) and (2) have been obtained in refs 53 within a percolation approach for the shallow hydrogenic levels, satisfying the condition of E A /E g <<1. With the values of E g ~2.1-2.3 eV and E A = 140 meV, given above, one finds E A /E g ~0.06-0.07, in agreement with the condition above. Indeed, the plots in the middle panel of Fig. 1 exhibit broad linear intervals, pertinent to the Mott VRH conduction according to Eq. (1). The values of T 0 have been found with the slopes of these plots. The temperature intervals, ΔT v , of the Mott VRH charge transfer, have been obtained with the linear intervals of the corresponding plots. The data of T 0 and ΔT v are collected in Table 1. In addition, W in Table 1 has been evaluated with the expression of W ≈ 0.5 k B (T v 3 T 0 ) 1/4 40, 53 , where T v is the Mott VRH onset temperature on cooling (i. e. the right border of the intervals ΔT v in Table 1).

Density of the acceptor states.
In order to proceed the analysis of the resistivity and MR, below we assume a conventional DOS of the Anderson type 51,52 in the acceptor band, exhibited in Fig. 4. It is considered to be symmetrical, centered at E A = 0 and is characterized by a finite total width 2W. The parameters −E c and E c in Fig. 4 are the mobility edges, separating the delocalized acceptor states from the localized states of the acceptor band, which are hatched. Finally, a small overlap of the acceptor band with the states of the valence band, taking place for #1 (in Table 1, W = 146 meV exceeds slightly E A = 140 meV 64 ) can be neglected in the analysis below. Indeed, the position of μ is expected to lie among the localized states, because otherwise the conductivity would be metallic 51,52 contradicting to the activated behavior of ρ(T) in Fig. 1. Such position should be chosen near the left edge of the DOS, as shown in Fig. 4. This corresponds to a sufficiently high degree of compensation, where N A is the acceptor concentration and N D is the concentration of the compensating donors. Therefore, the details of the energy spectrum on the right edge of the DOS, which is close to the valence band edge, are unimportant. In turn, this position of μ is connected with the high W values, which are comparable with the E A value (see Table 1). This excludes the possibility of μ lying among the right-hand interval of the localized states in Fig. 4, because in this case the conductivity, connected to the activation of the holes from μ to the valence band states, would dominate at any T due to a small difference between W and E A . Instead, another type of the resistivity 51, 52 , a c has been predicted for the position of μ shown in Fig. 4. Namely, as outlined in Fig. 4, the charge transfer connected to Eq. (3) is determined by activation of the localized holes from μ to −E c , or, more strictly, from the localized states below the Fermi level into the interval of delocalized states of the impurity band 51,52 . In turn, the temperature interval of Eq. (3), ΔT a , may lie either above or below that of the VRH interval ΔT v , depending on the distance between μ and -E c according to Eqs (3) and (4). In particular, both cases of the position of ΔT a with respect to ΔT v have been observed in kesterites 39,41,50 .

Analysis of magnetoresistance. In doped semiconductors pMR in the Mott VRH conduction regime is
connected mainly with the shrinkage of the impurity wave functions by the magnetic field 53 . In weak magnetic fields of λ(B)>> a, where λ(B) = [ħ/(eB)] 1/2 is the magnetic length (here, ħ is the Planck constant and e is the elementary charge), pMR of this mechanism can be expressed as 53  at sufficiently low temperatures, where the inelastic collision frequency is expected to be also low, providing conservation of a carrier phase coherence during multiple scattering 65 . However, in strongly disordered systems the constraint above may be weakened, which expands the temperature interval, where the interference effects are important. The point is that a higher disorder, being eventually a main reason of the interference effects [65][66][67][68][69] , would lead to their stronger contribution to MR at the same T. Therefore, one can expect a corresponding increase of the upper temperature limit, where the interference effects still survive, when the disorder is increased. Indeed, nMR in the VRH regime, attributable to the interference mechanism above, has been observed at T up to 300 K in disordered polyaniline silicon nanocomposites 70 and highly disordered pregraphitic carbon nanofiber 71 , as well as up to 280 K in Cu 2 Zn(Sn x Ge 1-x )Se 4 single crystals 11 . As mentioned above, CZGeS and related quaternary compounds are characterized, generally, by high intrinsic lattice disorder. In particular, such disorder leads to a considerable broadening of the impurity band, which exceeds substantially typical values of W in doped semiconductors. As can be seen in Table 1, large values of W characterize our CZGeS samples, too. This permits us to use tentatively the orbital interference mechanism [65][66][67][68][69] to interpret the nMR effect in our material.
For the weak scattering, given by the relation of 68 , and below a critical field B c , nMR can be expressed as Δρ/ρ = −a 2 (T)B 2 , whereas above B c it is given by the expression and J 0 is the prefactor of the overlap integral 68 .
357 is a constant 72 ), the nMR law of Δρ/ρ ∝ B 1/2 has been predicted 69 . First, we estimate the value of B c ′ for our CZGeS samples. As evident in the top panel of Fig. 1, the strongly activated behavior of ρ(T) is consistent with the state of all the samples lying relatively far from the metal-insulator transition, meaning that the value of a should exceed only slightly the value of the Bohr radius, a B = ħ 2 κ 0 / (me 2 ). For evaluation of a B , the theoretical values of the dielectric constant far from the metal-insulator transition, κ 0 = 6.8 73 , and of the mean hole effective mass of a stannite CZGeS, m = 0.48m 0 (where m 0 is the free electron mass) 74 can be used, yielding a B = 7.5 Å. This result is supported completely by the value of E A = ħ 2 / (2ma B 2 ) = 140 meV, coinciding with its value observed experimentally in ref. 64. Hence, taking e. g. a = 10 Å, we obtain with the direct expression of B c ′ = ħ/(a 2 eγ 4 ) × (T/T 0 ) 1/4 , the values of B c ′ ≈ 100-140 T at T = 100 K and 110-160 T at T = 200 K. These data permit to exclude a possibility of the square-root dependence of nMR on B above from further consideration, meaning that B c ′ exceeds considerably the maximum applied field of B = 20 T irrespective the possible error of the a value used above. Similarly the magnetic field B 2 , addressed to violation of Eq. (5) and defined by the condition of λ(B 2 ) ~ a, i. e. B 2 ~ ħ/(ea 2 ), can be estimated to lie close to B 2 ~600 T. This guarantees applicability of Eq. (5) for our samples, too.
Eventually, the data of A s ≈ 0.040-0.049, 0.12-0.14 and 0.055-0.064 can be evaluated for #1, #2 and #3, respectively, with the expression of A s , given below Eq. (7) directly, using E A = 140 meV and the values of T 0 and ΔT v from Table 1. This guarantees the regime of the weak scattering in agreement with Eq. (7). However, the value of the critical field B c dividing the field intervals of the quadratic and linear nMR dependences above cannot be estimated beforehand requiring the value of the parameter J 0 . On the other hand, as evident in Fig. 2, the field interval of Δρ/ρ ∝ B 2 in our samples looks to be quite narrow and is shifted to the low fields (lying probably below B ~1-2 T), constituting only a negligible part of the investigated field interval, where the MR data are less confident than at higher B. Therefore, it is more consistent to concentrate the analysis of MR on the higher field interval of B > 10 T, where the linear field dependence of nMR and the quadratic dependence of the pMR on B are more expected to persist. Hence, taking into account the consideration of nMR above, Eq. (5) for pMR and the evident fact that the expression of ln[ρ (T, B)/ρ (T, 0)] ≈ Δρ/ρ is fulfilled with a high accuracy, provided that the relative MR does not exceed 1% anywhere (Fig. 2), we can use for the analysis of MR the expression The values of a 1 (T) are found by plotting the left-hand side of Eq. (9) vs. B 2 and choosing a 1 (T) to obtain the best linearization of the plots. This is done by minimizing the standard deviation (SD), provided that all the plots should pass through the origin. Finally, the slope of these plots in Fig. 5 gives A 4 (T).
As can be seen in Fig. 5, the procedure above can be done for the majority of the MR plots in Fig. 2, including those lying within the VRH interval ΔT v (some of the data are not shown only for convenience), whereas the onset of the linear dependence, lying between B ons ~5-12 T, can be found only approximately. However, the interval of linearity looks sufficient for a reasonably accurate determination of both parameters, a 1 (T) and A 4 (T), which are plotted in the top and middle panels of Fig. 6, respectively, as functions of T −3/4 . One can see a good linearity of both functions within a whole temperature interval ~100-200 K of the Mott VRH conduction regime (cf. Table 1), deviating from such behavior only below 80 K. The values of A 04 were obtained with linear fits of the plots in the middle panel of Fig. 6 according to the expression of A 4 (T), given below Eq. (5). They are collected in Table 1.
Finally, existence of other nMR mechanisms, probably less universal than that discussed above, should be mentioned for the sake of completeness. In particular, we have tested the nMR models based on the Zeeman effect 75 and on a possible sensitivity of DOS to the magnetic field 76 , as well as the Khosla-Fischer model, taking into account both nMR and pMR contributions and addressing nMR to the scattering of charge carriers by the localized magnetic moments 77 . Concerning the first two models 75,76 , their application could not reproduce correctly either the field or the temperature dependences of nMR in our samples. Formally, our MR data can be reproduced with the last out of the models above 77 . However, this model requires a degeneracy of the carrier gas, which obviously does not take place in our samples contradicting to the strongly activated behavior of ρ(T) in Fig. 1. In addition, no evidence for existence of the localized magnetic moments (or any types of paramagnetic centers) have been obtained in CZGeS in any previous investigations. In particular, those of the magnetization 78,79 and the EPR 79 measurements have established only a weak diamagnetic response, which does not permit to consider any model of nMR, based on localized magnetic moments, to be substantiated sufficiently.
Determination of microscopic parameters. The first pair of the hole parameters, namely those of a and g(μ), can be found directly from the pair of Eqs (2) and (6). This is done using the data of T 0 and A 04 (which have been obtained from the slope of the plots in the middle panel of Fig. 6) collected in Table 1, irrespective to the details of the DOS model shown in Fig. 4. The values of a and g(μ) are collected in Table 2.
However, further analysis designed for determination of such details of the hole spectrum, as positions of μ and E c , as well as those of N A , requires utilization of the DOS model, although without a detailed knowledge of its shape. At this point, the analysis above is valid for both possible positions of μ, namely, near −W and W, since a symmetrical DOS has been proposed. On the other hand, the arguments towards the strong compensation, formulated in Section "Density of the acceptor states", still dictates the position of μ as shown in Fig. 4. Further argument for such a choice will be obtained below.
The general expression of the localization radius can be written in a form where ν ≈ 1 is the critical index of the correlation length and N c is the critical acceptor concentration [51][52][53][54] . The latter is connected to the Bohr radius with the universal Mott criterion, ≈ .
obtained with Eq. (2) and (10). The data of N A are collected in Table 2.
On the other hand, according to ref. 51, another expression of a, conformed with the energy scale of Fig. 4, can be written as The values of μ and E c are obtained with Eqs (12)(13)(14), using the value of κ 0 = 6.8 cited above, the value of a B = 7.5 Å and those of N A in Table 2, by variation of μ to fit the data of a in Table 2 with Eq. (6) explicitly for each sample. The data of W in Table 1 and N A in Table 2 permitted evaluation of the average DOS values, g av ≡ N A /(2 W), which are collected also in Table 2. For convenience, the values of μ and E c are displayed in a separate Table 3 below (see section "Discussion"). Finally, the data of E a have been calculated with Eq. (4), using those of μ and E c (Table 3), and are collected in Table 3, too.
Hall coefficient. The Hall effect in the domain of the VRH conduction has different nature, than in the case of the band conductivity over the delocalized states (or drift conduction), connected to the Lorenz force. Therefore, the sign of the Hall resistivity, ρ H , and of the Hall coefficient, R H , should not be obligatory the same as that of the charge carriers. Namely, it can be opposite to the sign of the thermopower, which has been observed in the amorphous, as well as in the doped crystalline semiconductors 54, 80-85 , including our case of CZGeS (see Fig. 3).
The early theories of the Hall effect under VRH conduction predicted the behavior of R H (T) similar to ρ(T), given by Eq. (1), but with another (and much smaller) value of the characteristic temperature T 0H ≈ 0.15T 0 due to the influence of the exponential dependence of the Hall mobility, μ H (T) = R H (T)/ρ(T). However, this prediction has not obtained a comprehensive experimental support, meaning observations of R H (T) with T 0H both smaller than T 0 and close to T 0 (see ref. 54 and references therein). Moreover, it has been even suggested in investigations of CuInSe 2 , that the dependence of R H (T) is not exponential at all, whereas the resistivity-like behavior with the same value of T 0 should be addressed only to the Hall mobility 81,82 .
On the other hand, further theoretical work, based on an explicit percolation approach to the treatment of the Hall effect 80,[83][84][85] . has given the evidence for only a weak (power-law) temperature dependence of μ H (T) and for the exponential dependence of R H (T), same as that of ρ(T). Namely, according to the results of ref. 80, the temperature dependence of the Hall coefficient can be presented in the form where T 0H = T 0 and the dependence of C(T) is quite weak (logarithmic). As can be seen in Fig. 3, some scattering of the data points does not permit to establish deviations of ρ H (B) in # 1 from linearity unambiguously, excluding probably only the case of low fields at 70 K. Therefore, the plots in Fig. 3 have been fitted with the linear function within the whole interval of B to obtain R H (T), and the plots of ln(−R H ) vs. T −1/4 , neglecting the very weak dependence of C(T), are presented in the bottom panel of Fig. 6. The linear fit of the data yields the value of T 0H = (1.5 ± 0.4)×10 7 K.

Discussion
The good linearity of the plots in Fig. 5, accompanied with the dependences of a 1 (T) and A 4 (T) in the top and middle panels of Fig. 6, respectively, provide a sufficient evidence of the interpretation of both nMR and pMR contributions to the total MR of our CZGeS samples, proposed in Section "Analysis of magnetoresistance". In addition, the analysis given in previous Section permits a more straightforward evaluation of the critical field B c according to Eq. (8), using the explicit values of a in Table 2, those of T 0 in Table 1 and J 0 = 3.0, 2.4 and 2.5 eV for #1, #2 and #3, respectively, following from Eq. (14). So, one can obtain the values of B c ≈ 3-9 T, which are comparable but smaller than those of B ons ≈ 5-12 T marking the onset of the linear behavior in Fig. 5. Such a relation between B c and B ons is in a reasonable agreement with the linear contribution of nMR, as discussed in Section "Analysis of magnetoresistance". On the other hand, the deviations from linearity of the plots of A 4 vs. T −3/4 and especially those of a 1 vs. T −3/4 (see the top and middle panels of Fig. 6) already close to the lower border of the VRH conduction (i. e. below ~80 K) look too strong. This requires a special discussion. At this point, only the orbital interference nMR mechanism may be insufficient to account for all details of the temperature dependence of nMR in our samples. Namely, the violation of the linearity of the plots in Fig. 6 with lowering T may be connected also to importance of the spin disorder and spin correlation effects of the hopping electrons with decreasing temperature [86][87][88] . Indeed, the temperature dependence of nMR due to the spin effects above is determined  ), and the interval of activated conduction (ΔT a ).
by the ratio of τ/τ s , where τ and τ s are the characteristic hopping time and the spin relaxation time, respectively, having non-trivial and quite different dependences on T 87,88 . Hence, the issue above is attributable also to the contribution of spin-disorder and spin-correlation effects, which are more important when T is decreased [86][87][88] . A similar situation with the anomalous temperature behavior of the nMR contribution has been observed recently in the magnetotransport of Cu 2 Zn(Sn x Ge 1-x )Se 4 single crystals 11 . For the next, as follows from the data in Table 2, the small values of the ratio of N A /N C ~0.2-0.4 agree completely with those of another ratio, a/a B ~ 1.3-1.7 exceeding unity only slightly. As can be seen with Eq. (10), the values of both ratios above mean, that all the investigated CZGeS material lies relatively far from the metal-insulator transition. In addition, one can see in Table 2, that the values of g(μ) are comparable with the average DOS values, g av , but are systematically smaller. The comparability of g(μ) and g av supports the consistence of our analysis in general, whereas the smallness of g(μ) with respect to g av is evident from Fig. 4, provided that the position of μ is shifted towards one of the DOS edges.
Eventually, the shift of μ towards the left DOS edge, or −W, suggests that the degree of the compensation in our samples is sufficiently strong. The arguments for such a situations have been already given in Section "Density of the acceptor states". In addition, the values of E a have been evaluated in Section "Determination of microscopic parameters" above (see Table 3). Such values of E a suggest existence of the intervals ΔT a of the activated ρ(T) behavior, given by Eq. (3), which is not connected to the Mott VRH behavior within the corresponding interval ΔT v , but is due to the activation of the holes into the interval (−E c , E c ) of the delocalized states of the acceptor band (see Section "Density of the acceptor states", inset in the top panel and the bottom panel of Fig. 1). Namely, the relatively high value of E a for #1 should shift ΔT a above ΔT v , whereas small values of E a for #2 and #3 suggest that ΔT a should lie below ΔT v . The corresponding intervals ΔT a have been found by the linearization of the plots in the inset to the top panel of Fig. 1 and in the bottom panel in Fig. 1, and are given in Table 3. The values of E a (ex) , obtained from the slopes of the linear fits of these plots, are displayed in Table 3, too. One can see that, indeed, ΔT a for #1 lies above Δ T v , while those for #2 and #3 lie below ΔT v . Such behavior is accompanied with the experimental values, E a (ex) , lying quite close to those of E a , evaluated in previous Section (see Table 3). This supports completely the position of μ in Fig. 3, as well as clarifies the nature of the charge transfer in CZGeS outside the Mott VRH interval, indicating a sufficiently strong degree of the compensation K in our material, as has been supposed above. On the other hand, the too strong compensation does not favor the nMR effect in general 39,65 , which in turn may limit the value of K from above.
Finally, one can see that the Hall coefficient (bottom panel of Fig. 6) vary with T sufficiently close to the law of Eq. (15). In addition, the value of T 0H = (1.5 ± 0.4) × 10 7 K is close to that of T 0 = 1.43 × 10 7 K for #1, following from the ρ(T) dependence of Eq. (1). This suggests a reasonable agreement with the percolation model of the Hall effect in the domain of the Mott VRH conduction 80 . However, some deviations from the model of ref. 80, including a possible exponential contribution to the Hall mobility, cannot be excluded due to a rather high error of T 0H .
From the presented results, it is evident that the charge transport mechanisms in CZGeS have a complicated nature. In particular, this leads to observations of the VRH conductivity at unusually high temperatures, including the nMR and the pMR contributions in non-zero magnetic field, as well as to the non-conventional Hall effect. All these findings, including the variation of the macroscopic and microscopic parameters from sample to sample, should be considered anyway during the production of CZGeS based optoelectronic devices. Particularly for the solar cell use, the key problem of this material is its high tolerance to the intrinsic defect formation. Confirmed deep position of defect acceptor levels, as well as their broadening into the quite wide acceptor bands, may have an important detrimental influence on the photo carriers recombination and the final device efficiency. Therefore, the growing conditions of the CZGeS absorber layer for the thin film solar cells should be optimized to prevent or at least minimize the formation of the deep defects. In addition, the non-conventional Hall effect, which does not have any direct connections to the hole concentration, should be taken into account with a certain cautions for the characterization of CZGeS materials.

Conclusions
We have investigated the resistivity, the magnetoresistance and the Hall effect in p-Cu 2 ZnGeS 4 single crystals in pulsed magnetic fields up to 20 T. The dependence of ρ(T) in zero magnetic field is described by the Mott VRH charge transfer mechanism within a broad temperature interval of ~100-200 K. Magnetoresistance contains the positive and negative components, which are interpreted by the common reasons of doped semiconductors, including the shrinkage of the impurity wave functions by the magnetic field and the damping of the electron interference effects in the VRH hopping conduction regime, respectively. On the other hand, the joint analysis of the resistivity and MR data has yielded a series of important electronic parameters and permitted specification of the CZGeS conductivity mechanisms outside the temperature intervals of the Mott VRH conduction. The Hall coefficient is negative, exhibiting an exponential dependence on temperature, which is quite close to that of ρ(T), as typical of the Hall effect in the domain of the VRH charge transfer.

Materials and Methods
Single crystals of CZGeS were grown by a chemical vapor transport using iodine as a transport agent. The growing process was performed in the vertical two-zone furnace with 850 °C in the evaporation zone and 800 °C in the growth zone. More details of the crystal growth process could be found elsewhere 20 .
Three CZGeS samples were selected for detailed investigations by a most convenient shape for electrical measurements. Chemical composition of the samples was determined by the X-ray fluorescence method, performed in three different points of each sample. The obtained data do not vary significantly from point to point, and the mean values are presented in Table 4. It can be seen, that all samples have the composition close to stoichiometry, exhibiting a slight Cu and Zn excess for #2 and some Cu and Zn deficient for #3 and #1. A part of the obtained material was grounded for the X-ray diffraction analysis, which showed the wurtzstannite structure of all the samples 20 .
The hot point probe method, addressed to the thermopower measurements, showed p-type conductivity in all samples. Six indium contacts were made in each sample. The resistivity, ρ(T), was measured with a standard dc method, and MR was investigated in the pulsed magnetic field with B up to 20 T (see ref. 39 for the device details). The sample was fixed on the holder and inserted inside a cryostat. A filling helium cryostat was used to vary the sample temperature in the range of 20-300 K. A finger of the cryostat, containing the sample, was installed inside a pulsed solenoid. The main sample plate was oriented perpendicular to the direction of the magnetic field. A measurement at a certain T was performed by increasing and decreasing B from 0 up to the values between ~0.5-20 T, with subsequent changing of the field polarity. The results with the opposite polarities were averaged with a special software, to obtain the parallel (ρ) and the perpendicular (ρ H ) components and to avoid their mutual influences. Eventually, series of above measurements was integrated to obtain the dependence of ρ(B) and ρ H (B) over the whole field diapason. Finally, it should be noted that the measurements at lower temperatures were hindered by the high sample resistance, which exceeded somewhat the installation limits.