The Scientific World Journal

Volume 2014 (2014), Article ID 737296, 8 pages

http://dx.doi.org/10.1155/2014/737296

Research Article

Weak Localization in Graphene: Theory, Simulations, and Experiments

Michael Hilke,¹ Mathieu Massicotte,² Eric Whiteway,¹ and Victor Yu¹

¹Department of Physics, McGill University, Montreal, QC, Canada, H3A 2T8

²Institute of Photonic Sciences, Castelldefels, 08860 Barcelona, Spain

Received 7 February 2014; Accepted 20 May 2014; Published 9 June 2014

Academic Editor: Jau-Wern Chiou

Copyright © 2014 Michael Hilke et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract.

We provide a comprehensive picture of magnetotransport in graphene monolayers in the limit of nonquantizing magnetic fields. We discuss the effects of two-carrier transport, weak localization, weak antilocalization, and strong localization for graphene devices of various mobilities, through theory, experiments, and numerical simulations. In particular, we observe a minimum in the weak localization and strong localization length reminiscent of the minimum in the conductivity, which allows us to make the connection between weak and strong localization. This provides a unified framework for both localizations, which explains the observed experimental features. We compare these results to numerical simulation and find a remarkable agreement between theory, experiment, and numerics. Various graphene devices were used in this study, including graphene on different substrates, such as glass and silicon, as well as low and high mobility devices.

1. Introduction

Graphene has attracted a considerable amount of attention due to the ease in isolating a single sheet of graphite via mechanical exfoliation [1, 2]. Despite the fact that it is only one- atom thick, exfoliated graphene has shown extraordinary transport properties and can be used as a novel material for many potential applications. Its unique band structure has led to many interesting phenomena such as tunable charge carriers densities [3], anomalous integer quantum Hall effect [1, 4], and ultrahigh mobilities at room temperature [5]. The most noteworthy property of the band structure is the existence of two degenerate Dirac cones [6], which leads to two degenerate valleys (K and K′). The existence of these valleys, each with a linear dispersion, is the main reason that graphene has transport properties which differ from most other semiconductors or semimetals. At high magnetic fields, a striking example is the quantum Hall effect, where the Hall conductance quantization occurs in steps of , because of the spin and valley degeneracy, and the Landau level quantization is proportional to , where is the perpendicular magnetic field component. This square root dependence leads to a very large lowest Landau level splitting, where the field quantization has been observed up to room temperature [7]. While the fourfold degeneracy can be lifted at very high fields in high mobility graphene, it is invisible at lower fields, which is the focus of this work.

The goal of this work is to present a comprehensive understanding of the magnetoresistance of monolayer graphene in a regime where the Landau quantization is not important; that is, , where is the quantum mobility. In this regime, the Landau quantization does not play any role, but the valley degeneracy, the linear dispersion, and the underlying disorder potential lead to interesting magnetotransport phenomena such as two-carrier transport (2CT), weak localization (WL), weak antilocalization (WAL), and strong localization (SL).

Most of the magnetotransport properties of graphene can be understood in terms of two types of scattering mechanisms: intervalley scattering, where electrons are scattered from one valley to the other (K to K′) with a rate hereby noted as , and intravalley scattering, where electrons scatter within a valley as described by a rate . In general, intervalley stems from short range scattering, such as atomic defects, including grain boundaries, whereas intravalley scattering is long range and is typically stronger and also includes large scale inhomogeneities and charged impurities in the substrate. Long range potential variations are responsible for the existence of both carriers' types (electrons and holes) simultaneously at the charge neutrality point (CNP). The carrier distribution can be understood by assuming that the total carrier density is given by [8], where , is the gate capacitance, the gate voltage, and the gate voltage at the charge neutrality point (CNP) with a total residual density (in Dorgan’s [8] work, was defined as half of the total residual density). is composed of the thermal carriers plus remaining electrons and holes due to sample inhomogeneities as a consequence of the underlying disorder potential. This disorder induced residual density dominates at low temperatures. and are the carrier densities of electrons and holes, respectively. This allows us to define electron and hole mobilities, and , where and are the corresponding electron and hole conductivities. Experimentally, and are often slightly different, which is likely due to the asymmetry of the scattering potential for electrons and holes.

2. Two-Carrier Transport

We start our magnetotransport analysis by discussing the simplest nontrivial contribution, which is due to the simultaneous existence of the two types of carriers near the CNP. As long as , this contribution can be evaluated using Drude's expression for two carriers, where and for the two types of carrier densities ( or ). The total resistivity is then given by , where are the single band conductivity matrices. This yields the following field dependence of the magnetoresistivity:

From this it follows that the relative field dependence can be written as

where

Expressing the conductivities in (3) in terms of mobilities, equating the zero field resistivity to , and identifying as the zero field resistivity at the CNP, we obtain

Hence, as a function of the 2CT gives rise to a parabolic positive relative magnetoresistance before saturating to when , which is valid as long as . Close to the CNP, that is, , is simply equal to the inverse median mobility and , which is the regime where the 2CT effect is the largest. Away from the CNP , we have a maximum , which decreases away from the CNP, since , which is the Hall density, increases. The 2CT will play an important role, as long as , to give rise to positive magnetoresistance due to large scale inhomogeneities. This effect is very important in graphene as compared to other two-dimensional systems because of the absence of a gap between the electron and hole carriers. The 2CT will serve as basis for the understanding of the classical contribution due to long range scattering, which will be present at all temperatures and even increases with temperature, since increases with temperature due to the thermal activation of the electron and hole carriers.

3. Theory of Weak Localization

Moving beyond the Drude and classical description of transport, we need to include quantum effects, of which coherent backscattering is the most important contribution. Coherent backscattering leads to WL [9, 10]. To obtain an expression for the WL correction, we have to evaluate the return probability of all possible trajectories [11]. At zero magnetic field, the coherent return probability can be expressed as [12]

where the first term in the integrant is the return probability () for diffusion constant , the second term represents the short time cut-off (), below which no elastic scattering occurs, and the third term is the phase coherence time () cut-off, beyond which phase coherence is lost. At low temperatures, where , this leads to the typical logarithmic WL correction of the conductivity .

To evaluate the magnetic field dependence of the return probability , we have to solve for the return probability according to the field dependent diffusion equation [11]:

where is the vector potential. The solution to this diffusion equation can be evaluated for and is given by

where . At zero field we recover . In the long time limit, where we have , this leads to the field induced destruction of the phase coherence when approximately one flux penetrates the area . The full magnetic field dependence is obtained by inserting into the expression of the coherent return probability and evaluating the integral; that is,

where and are the characteristic magnetic fields for and , whereas is the digamma function. The magnetic field correction to the coherent return probability is then given by

where , for , for , and is the Euler constant. The second term in (9) is responsible for the reduction of the return probability due to the presence of a magnetic field, which leads to the observed peak in resistance of width ~ around zero magnetic field, characteristic of weak localization in disordered systems [13]. The first term renormalizes the coherent return probability at large fields; that is, . This expression leads to the usual WL correction to the conductivity when expressing the conductivity correction in terms of the coherent return probability correction; that is, [11, 14].

In deriving (9) we assumed the existence of only two types of scattering mechanisms (phase breaking) and (elastic). However, in graphene, we have two very different types of elastic scattering mechanisms: those that scatter electrons within a valley (intravalley, ) and those that scatter between valleys (intervalley, K to K′, ). Here only intervalley scattering contributes to coherent backscattering; hence, the short time cut-off scattering time is provided by in graphene. Intervalley scattering is typically induced by short ranged impurities and can be due to grain boundaries or lattice defects. However, intravalley scattering is usually much stronger in graphene but does not contribute to coherent backscattering.

In an important work by McCann and coworkers [14], based on earlier work in honeycomb lattices [15], the authors have obtained a general expression for the WL and WAL correction specific to graphene, which determines the dependence of the magnetoresistivity as a function of involving the scattering parameters and explicitly. They obtained the following expression:

where and . This expression can also be obtained directly using (5) by replacing the short time cut-off in the magnetic field correction by . While the first two terms in (10) lead to the usual WL localization effect as in (9), the last term leads to WAL due to the presence of intravalley scattering as illustrated in Figure 1. This term becomes important when approaches and leads to a positive magnetoresistance contribution. In general, short ranged disorder will lead to WL corrections, due to intervalley scattering, while long range disorder leads to WAL [15]. However, onsite long ranged impurity potentials, like charge impurities, can induce both WL and WAL corrections [16, 17]. The scattering rates for electrons and holes could be different, and the expressions are valid for each type of carriers independently.

Figure 1: The relative return probability is shown as a function of the magnetic field for the different models. The green curve is from expression (9), the blue curve is from expression (10) with , and the red curve is obtained from expression (19) with and . For all curves we assumed that . The inset shows a zoom-in of the region close to zero magnetic field.

4. Strong Localization

In some graphene samples disorder can be very important, for instance, in intentionally disordered exfoliated graphene, which was shown to lead to strong localization [18]. Strong localization or Anderson localization [19] is obtained when the transmission is exponentially suppressed due to coherent backscattering. In the language of Anderson localization, the important parameter is the localization length , which measures the exponential increase of the resistance with size. In a system where the coherence length is infinite, the inverse localization length can be defined by

where is the length of the system and is the resistance. In order to incorporate strong localization in our discussion on weak localization, we can define the localization time as , which represents the time before the charge carrier is localized in a diffusive system. The coherent return probability then has to be modified according to

where and . The effect of strong localization is therefore simply to replace by . This is very intuitive, since if , will play the role of in determining the magnetic field correction, and no closed trajectory can occur beyond . The relative change in coherent return probability can then be written as

where the field dependence is mainly determined by . Since we define the localization length assuming a coherent system, the field induced relative change in is directly determined by the change in the relative coherent return probability for zero dephasing:

Equation (14) tells us how strong localization is connected to the coherent return probability and therefore to WL. For weak uncorrelated disorder, can be estimated using an extension of the Thouless expression in one dimension [20], which was recently used in the context of graphene [21], as , where is related to the quantum conductance. For small corrections in the conductivity, this expression yields . However, for strong disorder this expression is not valid and we will evaluate numerically below.

5. Graphene Nanoribbons

In many situations, both experimentally and numerically, one deals with a situation where the sample width is finite. This leads to a different field dependence, since the sample effectively becomes quasi-one-dimensional (Q1D) instead of 2D, when and is the width of the graphene ribbon. In this limit, the governing diffusion equation is no longer 2D but has to be replaced by the 1D analogue; that is, [12], where the first term represents 1D diffusion at zero field and the second term represents the field induced destruction of phase coherence in Q1D. is given as usual by , where scales with the flux quantum through the relevant characteristic area [14]. This leads to

When we obtain the following expression for the magnetic field correction of the coherent return probability:

Recalling that , we recover the expression for the WL correction to the conductivity obtained by McCann and coworkers [14].

For the strong disorder case, we have to again introduce the localization time , which has the effect of simply replacing by . We then obtain for the relative change in the coherent return probability

where . This is the result for the Q1D limit, when or . Here again we have for the relative change of the localization length.

6. Numerical Simulations

We now move to evaluate numerically for graphene assuming that , . There are two different limits: 2D and Q1D; that is, (the high disorder limit) and (the low disorder limit), respectively. These two limits are mainly determined by the amount of disorder in the system. For the numerical simulations, we will only consider short ranged disorder, since this is the relevant source for WL. The effect of long range disorder is to induce the coexistence of electrons and holes close to the charge neutrality point as well as WAL.

To model the system accurately, we consider a tight binding model with a honeycomb structure with two types of edges: armchair and zigzag. The hopping term is given by  eV. The two ends of the graphene device are assumed to be connected to wide ohmic contacts with a quadratic dispersion. The short range disorder is implemented by adding a random onsite potential on every atom site uniformly distributed with . The magnetic field is simply a phase factor in the hopping term, where corresponds to one flux quantum in a hexagon. The two-terminal transmission probability is then evaluated iteratively as a function of the length, , of the system, using an efficient iterative Green’s function method [22]. Similar models for disordered graphene have been used previously [23, 24]. The average two-terminal conductance , resistance , or logarithm conductance is obtained by averaging over many disorder configurations . For a given and at large enough , the system is always localized due to Anderson localization and can be extracted whenever .

We now move to compare the numerically obtained with the analytical expressions derived above for the relative change in , that is, testing (14). In the 2D limit, where , the field dependence of the coherent return probability is determined by (13), where in the coherent limit. The numerical field dependence of the relative is then rescaled by , where , and can be compared to the relative coherent return probabilities as shown in Figure 2. The agreement is quite remarkable with only two fitting parameters and . Away from the Dirac point we expect to have , where is approximately the number of quantum channels in a quasi-one-dimensional system. We indeed obtained for the comparison to the numerical data where in units of the lattice constant.

Figure 2: We compare the analytical expression (solid line) for from (13) with the numerically obtained relative change in the localization length using the parameter . The different dot colors (green, blue, and red) represent different widths (40, 80, 120) of the simulated graphene device, where is the lattice constant.

In the Q1D limit, where , the field dependence is rescaled by , where as determined from (17) and shown in Figure 3.

Figure 3: We compare the analytical expression (solid line) for from (17), relevant to the quasi-1D situation, with the numerically obtained relative change in the localization length using the parameter . The different dot colors (blue, red) represent different widths (80, 120) of the simulated graphene device, where is the lattice constant.

In both limits, the agreement between the numerically determined relative localization length and the relative coherent return probability is excellent, which confirms that indeed . Hence, the field dependence of the coherent return probability determines , which in turn determines , and therefore provides us with a way to extract the localization length simply from the magnetic field dependence.

It is now possible to make the connection between and the resistance, since . Therefore, small variations in lead to . This is assuming that . However, in experiments, where this is often not the case, plays the role of effective sample size () so that when we can write instead. This is often used in experiments to extract from the temperature dependence of the resistance, since also depends on temperature [25]. Here, we restrict ourselves to a fixed temperature, where is constant, but is magnetic field dependent. Hence, when we have, using (14),

For small variations , different contributions to the resistance are additive; hence the total effect on the resistance, including the two-carrier effect, yields

which allows the extraction of all relevant length scales from the magnetic field dependence alone, recalling that , , , and .

7. Experiments

While WL was observed in graphene in a number of previous experiments [26–29] and fitted using McCann’s WL expression [14], they were generally limited to low disorder. Here we will use expression (19) which is derived above and valid for all disorder strengths and fit it to the experimental data. We performed experiments on large scale graphene as well as lithographically defined Hall bars and graphene nanoribbons, in addition to large (over 100 m) single crystal graphene. Monolayers of graphene were grown by chemical vapor deposition (CVD) of hydrocarbons on 25 m thick commercial Cu foils. The CVD process used was similar to those described in previous works [30–33].

We start with low mobility samples, which cannot be fitted with the standard WL theory [14] alone because of the importance of strong localization in this regime. Instead we use the expression derived in (13), which incorporates the effect of strong localization. In Figure 4 we show the relative resistance change of a representative low mobility graphene Hall bar. Typical features of low mobility samples are a wide peak around zero magnetic field and an important relative resistance change. Using the best fit over the entire magnetic field range, we obtain  nm and  nm, which shows that the localization length is of the same order of magnitude as the phase coherence length. The agreement between the fit and the data is quite remarkable over the entire available magnetic field range and only relies on three fitting parameters, in this case , , and . The other parameters and are too large to significantly affect the fit. For even higher disorder, we expect to be smaller than and we would need to modify our assumption of the linear approximation used in (18). Instead, we expect the change in magnetoresistance to depend exponentially on and its associated field dependence. This is amplified further for long but narrow samples (nanoribbons).

Figure 4: The measured magnetoresistance of a low mobility ( cm²/Vs) graphene Hall bar device at 100 mK. The relative resistance (blues dots) and the fit using (19) with fitting parameters  mT,  T, and (red line) are shown. and are larger than 100 T and do not affect the fit.

The next step is to look at a large scale (cm) graphene sample grown by CVD deposited on a glass slide with evaporated gold as contacts. The magnetic field dependence is shown in Figure 5, which shows a relatively broad peak at zero magnetic field. There is also a parabolic increase in resistance at large field, which we attribute to 2CT. Using (19) to fit the data we extract the following length scales  nm, m, and  nm. Here again, the agreement between the fit and the data is quite remarkable, considering that the fit extends over the entire magnetic field range.

Figure 5: The measured magnetoresistance of a large scale graphene sample on glass (blue dots). The relative resistance and the fit (red line) using (19) with fitting parameters  mT, ,  T, and  T are shown.

The next sample we are considering is a large scale (cm) sized graphene sample grown by CVD and deposited on a standard SiO₂/Si substrate, where the doped silicon can be used as a back gate. The main difference is that this sample was grown using isotopically pure C¹³ methane gas instead of C¹², so that the Raman peaks are shifted by [34]. However, the different isotope does not affect the magnetotransport. The magnetic field dependence is shown in Figure 6, which shows a narrow weak localization peak at zero magnetic field. There is also a large parabolic increase in resistance at large field, which we attribute again to the 2CT. Using (19) to fit the experimental data we obtain  nm, m,  nm, and a mean field effect mobility of  m²/Vs, which is common for large scale CVD graphene, and a residual density of  cm⁻².

Figure 6: The measured magnetoresistance of a large scale graphene sample on SiO₂/Si (blue dots). The symmetrized relative resistance and the fit (red line) using (19) with fitting parameters  mT, ,  T,  mT, and  T and a zoom-in of the low field part in the inset are shown.

Interestingly, this is a prime example of a large parabolic-like positive magnetoresistance background, which cannot be fitted using WAL as per (10), since WAL gives negatively curved magnetoresistance at large fields. Here, we clearly demonstrate that the origin of the large positive magnetoresistance background is likely due to the 2CT because of the coexistence of both and type carriers. Strikingly, the quality of the fit is very remarkable over the entire available experimental magnetic field range.

We now turn to the higher mobility case, where we consider a graphene field effect device made out of a large (m diameter) single graphene crystal grown by CVD [35]. Because of its single crystal nature, this sample has a mobility comparable to exfoliated graphene on silicon oxide. From the gate voltage dependence we extracted a field effect mobility of about  m²/Vs for electrons and  m²/Vs for holes. The CNP was at  V. Using (19) to fit the experimental data shown in Figure 7 we can extract the various length scales and obtain  nm,  mm,  nm, and . For the median mobility, we used the experimental resistivity ratio of and found a remarkable agreement with the field effect mobility extracted from the gate voltage dependence. For the residual density we obtain  cm⁻² using .

Figure 7: The measured symmetrized magnetoresistance of a large single crystal graphene sample on SiO₂/Si measured at  V (blue dots) and at 0.3 K. The relative resistance and the fit (red line) using (19) with fitting parameters  mT, ,  T, and  mT are shown.

The overall experimental behavior can be summarized as follows: at low temperatures, very low mobility samples show a very wide peak in the resistance at zero field (as shown in Figure 4). In this regime it is important to include strong localization effects to understand the width of the peak, since is comparable to . The effective width is then determined by ~. With increasing mobility, the peak becomes increasingly sharper (Figures 5–7) and the relevant length scale becomes , whose inverse square determines the peak width. In addition to the WL effect (sharp negative magnetoresistance), which all samples show, there is often a parabolic positive magnetoresistance which is mainly due to the 2CT effect because of large scale inhomogeneities in addition to WAL in very clean samples.

8. Conclusion

We presented a comprehensive picture of the observed magnetoresistance in graphene spanning all disorder levels. At high disorder, localization becomes the dominant mechanism, where the field dependence can be understood in terms of the field dependent localization length, which interestingly follows the field dependence of the coherent return probability and hence the field dependence of the WL correction. In this picture, WL is simply a consequence of the field dependent localization length. Both WL and strong localization are due to short range scattering. Long range scattering causes a positive magnetoresistance effect on top of the negative WL magnetoresistance, which stems from WAL (due to intravalley scattering) and from the presence of two types of carriers close to the CNP (the 2CT). While WL is generic to all disordered systems, both WAL and 2CT are specific to graphene, where both play an important role in the understanding of magnetotransport.

Conflict of Interests

The authors declare that there is no conflict of interests regarding the publication of this paper.

Acknowledgments

The authors thank NSERC and FQRNT for financial assistance.

References

K. S. Novoselov, A. K. Geim, S. V. Morozov et al., “Electric field in atomically thin carbon films,” Science, vol. 306, no. 5696, pp. 666–669, 2004.

K. S. Novoselov, D. Jiang, F. Schedin et al., “Two-dimensional atomic crystals,” Proceedings of the National Academy of Sciences of the United States of America, vol. 102, no. 30, pp. 10451–10453, 2005.

D. R. Cooper, B. DAnjou, N. Ghattamaneni et al., “Experimental review of graphene,” ISRN Condensed Matter Physics, vol. 2012, Article ID 501686, 56 pages, 2012.

Y. Zhang, Y. W. Tan, H. L. Stormer, and P. Kim, “Experimental observation of the quantum Hall effect and Berry's phase in graphene,” Nature, vol. 438, no. 7065, pp. 201–204, 2005.

K. I. Bolotin, K. J. Sikes, Z. Jiang et al., “Ultrahigh electron mobility in suspended graphene,” Solid State Communications, vol. 146, no. 9-10, pp. 351–355, 2008.

P. R. Wallace, “The band theory of graphite,” Physical Review, vol. 71, p. 622, 1947.

K. S. Novoselov, Z. Jiang, Y. Zhang et al., “Room-temperature quantum hall effect in graphene,” Science, vol. 315, no. 5817, p. 1379, 2007.

V. E. Dorgan, M.-H. Baer, and E. Pop, “Mobility and saturation velocity in graphene on SiO2,” Applied Physics Letters, vol. 97, no. 8, Article ID 082112, 2010.

S. Hikami, A. I. Larkin, and Y. Nagaoka, “Spin-orbit interaction and magnetoresistance in the two dimensional random system,” Progress of Theoretical Physics, vol. 63, no. 2, pp. 707–710, 1980.

B. L. Altshuler, D. Khmel'nitzkii, A. I. Larkin, and P. A. Lee, “Magnetoresistance and Hall effect in a disordered two-dimensional electron gas,” Physical Review B, vol. 22, no. 11, pp. 5142–5153, 1980.

S. Chakravarty and A. Schmid, “Weak localization: the quasiclassical theory of electrons in a random potential,” Physics Reports, vol. 140, no. 4, pp. 193–236, 1986.

C. W. J. Beenakker and H. Van Houten, “Boundary scattering and weak localization of electrons in a magnetic field,” Physical Review B, vol. 38, no. 5, pp. 3232–3240, 1988.

G. Bergmann, “Weak localization in thin films. a time-of-flight experiment with conduction electrons,” Physics Reports, vol. 107, no. 1, pp. 1–58, 1984.

E. McCann, K. Kechedzhi, V. I. Falko, H. Suzuura, T. Ando, and B. L. Altshuler, “Weak-localization magnetoresistance and valley symmetry in graphene,” Physical Review Letters, vol. 97, Article ID 146805, 2006.

H. Suzuura and T. Ando, “Crossover from symplectic to orthogonal class in a two-dimensional honeycomb lattice,” Physical Review Letters, vol. 89, no. 26, Article ID 266603, pp. 1–4, 2002.

Y. Y. Zhang, J. Hu, B. A. Bernevig, X. R. Wang, X. C. Xie, and W. M. Liu, “Localization and the Kosterlitz-Thouless transition in disordered graphene,” Physical Review Letters, vol. 102, Article ID 106401, 2009.

F. Ortmann, A. Cresti, G. Montambaux, and S. Roche, “Magnetoresistance in disordered graphene: hhe role of pseudospin and dimensionality effects unraveled,” Europhysics Letters, vol. 94, no. 4, Article ID 47006, 2011.

J. Moser, H. Tao, S. Roche, F. Alzina, C. M. Sotomayor Torres, and A. Bachtold, “Magnetotransport in disordered graphene exposed to ozone: from weak to strong localization,” Physical Review B-Condensed Matter and Materials Physics, vol. 81, no. 20, Article ID 205445, 2010.

P. W. Anderson, “Absence of diffusion in certain random lattices,” Physical Review, vol. 109, p. 1492, 1958.

D. J. Thouless, “Localization distance and mean free path in one-dimensional disordered systems,” Journal of Physics C: Solid State Physics, vol. 6, no. 3, p. L49, 1973.

N. Leconte, A. Lherbier, F. Varchon, P. Ordejon, S. Roche, and J.-C. Charlier, “Quantum transport in chemically modified two-dimensional graphene: from minimal conductivity to Anderson localization,” Physical Review B-Condensed Matter and Materials Physics, vol. 84, no. 23, Article ID 235420, 2011.

S. Rotter, J. Z. Tang, L. Wirtz, J. Trost, and J. Burgdorfer, “Modular recursive Green’s function method for ballistic quantum transport,” Physical Review B, vol. 62, p. 1950, 2000.

A. Lherbier, B. Biel, Y.-M. Niquet, and S. Roche, “Transport length scales in disordered graphene-based materials: strong localization regimes and dimensionality effects,” Physical Review Letters, vol. 100, Article ID 036803, 2008.

A. Cresti, N. Nemec, B. Biel et al., “Charge transport in disordered graphene-based low dimensional materials,” Nano Research, vol. 1, no. 5, pp. 361–394, 2008.

X. Hong, S.-H. Cheng, C. Herding, and J. Zhu, “Colossal negative magnetoresistance in dilute fluorinated graphene,” Physical Review B, vol. 83, Article ID 085410, 2011.

S. V. Morozov, K. S. Novoselov, M. I. Katsnelson et al., “Giant intrinsic carrier mobilities in graphene and its bilayer,” Physical Review Letters, vol. 100, no. 1, Article ID 016602, 2008.

F. V. Tikhonenko, D. W. Horsell, R. V. Gorbachev, and A. K. Savchenko, “Weak localization in graphene flakes,” Physical Review Letters, vol. 100, no. 5, Article ID 056802, 2008.

D.-K. Ki, D. Jeong, J.-H. Choi, H.-J. Lee, and K.-S. Park, “Inelastic scattering in a monolayer graphene sheet: A weak-localization study,” Physical Review B, vol. 78, Article ID 125409, 2008.

Y.-F. Chen, M.-H. Bae, C. Chialvo, T. Dirks, A. Bezryadin, and N. Mason, “Magnetoresistance in single-layer graphene: weak localization and universal conductance fluctuation studies,” Journal of Physics Condensed Matter, vol. 22, no. 20, Article ID 205301, 2010.

A. Reina, X. Jia, J. Ho et al., “Large area, few-layer graphene films on arbitrary substrates by chemical vapor deposition,” Nano Letters, vol. 9, no. 1, pp. 30–35, 2009.

X. Li, W. Cai, J. An et al., “Large-area synthesis of high-quality and uniform graphene films on copper foils,” Science, vol. 324, no. 5932, pp. 1312–1314, 2009.

E. Whiteway, V. Yu, J. Lefebvre, R. Gagnon, and M. Hilke, “Magneto-transport of large CVD-grown graphene,” http://arxiv.org/abs/1011.5712

V. Yu, E. Whiteway, J. Maassen, and M. Hilke, “Raman spectroscopy of the internal strain of a graphene layer grown on copper tuned by chemical vapor deposition,” Physical Review B, vol. 84, Article ID 205407, 2011.

S. Bernard, E. Whiteway, V. Yu, D. G. Austing, and M. Hilke, “Probing the experimental phonon dispersion of graphene using ¹²C and ¹³C isotopes,” Physical Review B, vol. 86, Article ID 085409, 2012.

M. Massicotte, V. Yu, E. Whiteway, D. Vatnik, and M. Hilke, “Quantum Hall effect in fractal graphene: growth and properties of graphlocons,” Nanotechnology, vol. 24, no. 32, Article ID 325601, 2013.