Graphene, a sheet of sp²-bonded carbon atoms, for several years already continues to be in the focus of attention of physics community, due to its potentially transformative impact across a wide range of applications including advanced electronics and sensing. Carrier scattering in graphene can be due to all kinds of disorder, including ripples in the graphene layer, point defects and their associated short-range potentials, charged impurities residing in the supporting substrate, and adsorbed atoms on the surface. Investigation of the influence of disorder on the properties of graphene is attracting a tremendous interest due to possibility to modify this novel and promising material using weak or strong localization of charge carriers. By controllably introducing defects into graphene, one may be able to understand how these mechanisms limit transport.

Previously, there were observed separately either weak localization (WL) or different kinds of Variable- Range-Hopping (VRH) conductivity of strongly localized carriers in graphene samples disordered by different methods like doping, oxidation, ion irradiation (see, for example, [1] -[8] ). However, we are not aware of observations of all regimes of localization with gradual increase of disordering in graphene. In this paper we report the results of study of the localization process in monolayer graphene (MG) samples subjected by different doses of ion irradiation.

The initial large size (5 ´ 5 mm) specimens were supplied by “Graphenea” company. Monolayer graphene was produced by CVD on copper catalyst and transferred to a 300 nm SiO₂/Si substrate using wet transfer process. It is specified in the certificate, that the sheet resistivity of the specimen is 350 Ohm/sq. On one specimen, gold electrical contacts were deposited directly on the graphene surface. This sample was marked as 0. The resistivity of the sample 0 was 380 Ohm/sq which was close to the data in the certificate. On the surface of the other specimen, six groups of mini-samples (200 ´ 200 μm) were fabricated by means of electron-beam lithography (EBL) as well as electrical contacts (5 nm Ti and 45 nm Pd) for 2-probe measurements. The samples from the first group, marked as sample 1 were not irradiated, while 5 other groups were subjected to different doses (from up to cm⁻²) of irradiation by C⁺ ions with energy 35 keV. Ion irradiation was performed on HVEE-350 Implanter, C⁺ ions were obtained from the hollow cathode ion source by CO₂ decomposition.

The choice of the irradiation conditions is not occasional. They are the same as in [6] where the Hall Effect was measured in monolayer graphene subjected to ion irradiation. Resistances of the initial sample as in [6] and our initial sample 0 are the same. It is shown in [6] that both Hall coefficient and mobility are temperature independent up to high dose of irradiation. Moreover, it is shown that the sheet Hall concentration of charge carriers (about 10¹³ cm⁻²) does not depend on the dose of irradiation at all. We are based on these results in our decision to prefer the 2-probe geometry instead of the Hall bar geometry, because this allows us to measure more samples on the same sample holder. As a result, we are able to enhance the reproducibility of the obtained results. Comparison of the graphene resistance measured by 2-probe method and by 4-probe method presented in [6] shows that in this system, contact resistance is insignificant.

In our previous work [9] , a concentration of structural defects N_D was determined for each group of samples using measurements of the Raman scattering. These values are shown in the inset in Figure 1 (it turns out that sample 1 is also slightly disordered due to EBL process). Measurements of the current-voltage characteristics (I-V) for all samples, performed at room temperature, show that for samples 5 and 6, I-V is strongly non-linear even at very small current. That is why in this paper, the temperature dependences of resistance R(T) are shown only for samples 0 - 4.

The resistance is measured by two-probe method in helium cryostat down to 1.8 K in zero magnetic fields and in magnetic fields up to 4 Tesla. Figure 1 shows the general picture R(T) for all samples. The sample 0 shows typical metallic behavior, when R slightly decreases with decrease of T. For sample 1, R slightly increases with decreasing T, which is characteristic for “dirty” metals. For other samples R changes with T exponentially, which is characteristic for strongly localized carriers.

Figure 2 shows the experimental dependences of magnetoconductance (MC) of sample 1 in wide temperature interval, from 300 K down to 1.8 K. Plot of the temperature dependence of conductivity on the scale vs. ln T (Figure 3) shows the logarithmic temperature behavior of at low T, characteristic for regime of WL [10] , with tendency to saturation at very low temperatures.

WL regime of conductivity in monolayer graphene has important features due to the facts that charge carriers are chiral Dirac fermions, which are reside in two inequivalent valleys at the K and K' points of the Brillouin zone. Due to chirality, Dirac fermion acquires a phase of upon intravalley scattering, which leads to destructive interference with its time-reversed counterpart and weak antilocalization (WAL). Intervalley scattering leads to restoration of WL because fermions in K and K' valleys have opposite chiralities.

Quantum corrections to the conductivity of graphene have been intensely studied theoretically [11] -[19] . It is predicted that at relatively high temperatures WAL corrections will dominate, while with decreasing T the WL corrections will dominate. There were several experimental papers reporting logarithmic dependence of conductivity on temperature and magnetic field at low temperatures [1] -[5] [7] . However, in our sample 1, the logarithmic dependence is observed in wide temperature interval, starting from 300 K, which gives an opportunity to check in a very detailed way the theoretical predictions. For the MC the theory [15] predicts

where is the digamma function, is the coherence time, is the intervalley scattering rate, is the combined scattering rate of intravalley and intervalley scattering and of trigonal warping.

The result of fitting Equation (1) to experimental data for magnetoconductance of sample 1 at different temperatures is illustrated on Figure 2. In the process of fitting we are able to extract all three parameters entering the equation, which are shown in inset in Figure 3. It turns out, that the parameters are temperature-dependent which is not predicted by theory. Using these parameters we can calculate conductance at zero magnetic field according to Equation (10) of [15] , which can be rewritten in the form

where A is a constant, dependent upon the unit of magnetic field (chosen as 1T). We compare Equation (2) with the experimentally measured conductivity. The comparison is presented on Figure 3. The good agreement proves the correctness of the chosen parameters. The inset to Figure 3 shows that, apart from the lowest temperatures,. Mechanisms that can give such dependence are: electron-electron scattering in dirty limit [10] , electron-phonon scattering [20] [21] and electron-flexural phonon interaction [22] . The saturation of at low temperatures is well known in classical 2d systems and may be connected with existence of dephasing centers (for example, magnetic impurities) [23] [24] .Obtained value of allow us to determine the values of dephasing length. When the temperature decreases from 300 K to 3 K, the dephasing length increases from 7 nm to 70 nm and then saturates. The maximal value of allows us to estimate the density of dephasing centers as cm^-2.

Let’s discuss now samples 2 - 4 with pronounced insulating behavior. Plotting the data on the Arrhenius scale lnR vs. 1/T shows that energy of activation continuously decreases with decreasing T which is characteristic for the variable-range-hopping (VRH) conductivity [25] . There are two kinds of VRH depending on the structure of the density-of-states (DOS) in the vicinity of the Fermi level (FL) when, R(T) is described by the Mott law in the case of two-dimensional (2d) conductivity

Here is the numerical coefficient [25] , a is the radius of localization.

The Coulomb interaction between localized carriers leads to appearance of the soft Coulomb gap in the vicinity of FL which in the case of 2d has a linear form

where is the dielectric constant of the material. This leads to the Efros-Shklovskii (ES) VRH or law

where the numerical coefficient [25] .

Coulomb interaction can alter the DOS only near the FL. Far from FL, the DOS is restored to its initial value, which is approximately equal to, see inset in Figure 4. Denoting the half-width of the Coulomb gap as one can conclude, therefore, that at ES law has to be observed, while in the opposite case, the Mott law should dominate.

There are a number of reports about observation of either Mott or ES laws in different disordered graphene- based materials [26] -[29] . We show that in samples 3 and 4, both VRH laws are observable at different temperatures (For sample 2, the VRH regime will be observed at lower temperatures). In Figure 4 and Figure 5, log R is plotted versus and. At high temperatures, dependences R(T) are straightened on the scale, while at low temperatures they are straightened on the scale. The latter shows the approach to the ES law which should be observed at the lowest temperatures. These plots allow us to determine both parameters and (Table 1) and calculate the temperature of deviation from law to law in the case of 2d conductivity similarly to the calculation of for 3d conductivity [30] . In VRH, only localized states in an optimal band of width near the Fermi level are involved in the hopping process. The band becomes continuously narrower with decreasing. Hopping resistance is determined by the critical parameter

Here energy and temperature are measured in the same units, r is the mean distance of hopping. In the Mott regime, and, therefore, the total number of states in the optimal band is and the mean distance between states in two dimensions is. Substituting in Equation (6), one can find from the minimal value of determined from,

This gives the relationship between and the width of the optimal band:. At the crossover temperature T_c we get, so

Being inside the Coulomb gap, the crossover temperature can be determined from, which gives. Substituting into Equation (8) and using expressions for and, from Equation

(3) and Equation (5), we get

The functional proportionality of to the ratio has been obtained earlier [31] [32] , but with significantly different numerical coefficient. The latter is, however, crucial for comparison with experiment. The values of calculated from Equation (9) for samples 3 and 4 are given in Table 1 and shown as arrows in Figure 4. The good agreement shows the correctness of the obtained numerical coefficient. We can also estimate the width of the Coulomb gap. In calculation, the value for the monolayer graphene on the SiO₂ surface was used as, where (for SiO₂) and (for air). One can see that indeed, the ES law is observed when.

Comparison of samples 3 and 4 shows that increase of the density of defects leads to stronger localization which manifests in a significant decrease of a and increase of the energy band needed for hopping. This looks like increase of the amplitude of the random potential relief in classical semiconductors induced by randomly distributed positively charged donors and negatively charged acceptors in the case of compensation. Hence, we assume that structural defects in graphene are of amphoteric impurity action [33] , i.e. they can be either acceptors or donors and compensate each other. Thus, increase of leads to the random potential relief amplitude increase.

One can also assume that this phenomenon may be connected with the fact that some point defects (say, vacancies) may produce complex associative centers with other defects. These associations can show an amphoteric impurity action, whereas the individual components are not amphoteric [34] .

In conclusion, a gradual transformation of conductivity measured in a wide interval of temperatures (300 - 1.8 K) and magnetic fields (up to 4 T) was observed in a series of monolayer graphene samples subjected by different dose of ion irradiation. Increase of the density of structural defects induced by irradiation has led to change the mechanism of electron transport from metallic conductivity to the regime of WL and finally to the VRH of strongly localized carriers. It is shown that WL regime in slightly disordered sample starts from high (room) temperatures.

Comparison of experimental magnetoconductance curves with theory allows us to find the parameters which determine the logarithmic temperature dependence of conductivity in WL regime. In VRH, a crossover from the Mott law to the Efros-Shklovskii law is observed in the same samples with decreasing temperature. The calculated crossover temperatures are in good agreement with experimental values. It is suggested that strengthening of localization with increase of can be explained by amphoteric impurity action of graphene structural defects induced by ion irradiation.

ErezZion,AvnerHaran,AlexanderButenko,LeonidWolfson,YuriKaganovskii,TalHavdala,AmosSharoni,DoronNaveh,VladmirRichter,MosheKaveh,EugeneKogan,IssaiShlimak, (2015) Localization of Charge Carriers in Monolayer Graphene Gradually Disordered by Ion Irradiation. Graphene,04,45-53. doi: 10.4236/graphene.2015.43005