Scale-invariant large nonlocality in polycrystalline graphene
Mário Ribeiro, Stephen R. Power, Stephan Roche, Luis. E. Hueso, Fèlix Casanova
SScale-invariant large nonlocality in polycrystalline graphene
Mário Ribeiro
1, † , Stephen R. Power
2, 3 , Stephan Roche
2, 4 , Luis E. Hueso
1, 5, *, Fèlix Casanova
1, 5, *. CIC nanoGUNE, 20018 Donostia-San Sebastian, Basque Country, Spain Catalan Institute of Nanoscience and Nanotechnology (ICN2), CSIC and The Barcelona Institute of Science and Technology, Campus UAB, 08193 Bellaterra, Catalonia, Spain Universitat Autònoma de Barcelona, 08193 Bellaterra, Catalonia, Spain ICREA – Institució Catalana de Recerca i Estudis Avançats, 08010 Barcelona, Catalonia, Spain IKERBASQUE, Basque Foundation for Science, 48013 Bilbao, Basque Country, Spain †Present address: Center for Quantum Nanoscience, Institute for Basic Science (IBS), Seoul 03760, Republic of Korea and Department of Physics, Ewha Womans University, Seoul 03760, Republic of Korea. *Email: [email protected] , [email protected]
Abstract
The observation of large nonlocal resistances near the Dirac point in graphene has been related to a variety of intrinsic Hall effects, where the spin or valley degrees of freedom are controlled by symmetry breaking mechanisms. Engineering strong spin or valley Hall signals on scalable graphene devices could stimulate further practical developments of spin- and valleytronics. Here we report on scale-invariant nonlocal transport in large-scale chemical vapour deposition graphene under an applied external magnetic field. Contrary to previously reported Zeeman spin Hall effect, our results are explained by field-induced spin-filtered edge states whose sensitivity to grain boundaries manifests in the nonlocal resistance. This phenomenon, related to the emergence of the quantum Hall regime, persists up to the millimeter scale, showing that polycrystalline morphology can be imprinted in nonlocal transport. This suggests that topological Hall effects in large-scale graphene materials are highly sensitive to the underlying structural morphology, limiting practical realizations.
Introduction
In recent years, there has been a continuous effort to harvest the spin and valley degrees of freedom of charge carriers for developing innovative information processing as an alternative to complementary metal-oxide-semiconductor (CMOS) technologies . The high-mobility, low spin-orbit coupling (SOC), and linear energy dispersion of graphene has made it a core platform in such quest. In this context, graphene has been explored in magnetic-element-free nonlocal transport using Hall-bar geometries. The large nonlocal signals are observed close to the Dirac point, and have been tentatively related to topological Hall effects such as the spin Hall effect (SHE) , or long-range chargeless topological valley Hall currents for gapped graphene/h-BN heterostructures . Spin signals in graphene using nonlocal approaches free of magnetic elements have been reported by applying external magnetic fields , by using extrinsic sources of spin-orbit coupling , and by proximity effect to ferromagnetic insulators . Theoretical predictions have mainly explored spin diffusion in the spin Hall regime, where spin currents are respectively generated and detected by means of the SHE and the inverse spin Hall effects ISHE) . Nonlocal signatures related to valley Hall currents have also been predicted to emerge from the different diffraction features of the two valleys in graphene . The theoretical frameworks in the spin Hall regime consider extrinsic sources of SOC, Zeeman interaction (Zeeman spin Hall effect), and proximity-induced strong-exchange bias as possible mechanisms driving such enhancements of nonlocal signals close to the Dirac point . More recently, spin generation and detection via extrinsic sources of SOC have been questioned by Wang et al. , Kaverzin et al. , and Cresti et al. . All current demonstrations however rely on microscale Hall bars of pristine graphene obtained from micromechanical exfoliation techniques, with natural flakes with dimensions ranging several micrometers, or on microscale Hall bars patterned from chemical vapor deposition (CVD) graphene . The micrometer spin and valley relaxation lengths therefore restrict the experimental realizations to such scale. Core in the analysis of nonlocal signals, these approaches rely on comparing the nonlocal signal to the background arising from the classical current spreading (Ohmic) contribution and studying its dependence with the channel geometry to sustain the claim of other-than-charge sources for the nonlocality. Here, we demonstrate the persistence of large nonlocal signals at the Dirac point of CVD graphene devices from the micrometer up to the millimeter scale in presence of an external magnetic field applied perpendicular to the surface. These signals exhibit a similar dependence with the magnetic field to those induced via Zeeman spin Hall effect in microscale high-quality pristine graphene devices, and closely follow the same dependence with the device aspect ratio for device lengths that differ by two orders in magnitude. We however exclude an origin based on pure spin and valley currents due to the length scales involved in the transporting channel of the macroscale devices, and conclude that the origin of the intriguing nonlocality for both the microscale and macroscale stems from the same fundamental mechanism. Considering the microscopic details of the fabricated samples and the polycrystalline morphology of the CVD graphene samples, the large nonlocal signals are consistent with a dissipative quantum Hall regime driven by Zeeman-split counter-propagating edge states . Additionally, the grain boundaries shunt the insulating bulk induced by the strong magnetic field, and generate a strong asymmetry of the nonlocal magnetoresistance with external magnetic field, a striking feature already observed, but not yet understood . Results
Device fabrication and measurement configuration.
The availability of mm-sized continuous CVD graphene films enables the fabrication of Hall bars in a wide range of scales, and the study of nonlocal signals for device geometries not achievable using standard micromechanical exfoliation techniques . The devices studied in this work were fabricated following standard electron-beam lithography procedures on monolayer CVD graphene films wet-transferred onto 1×1 cm Si/SiO (300 nm) substrates (see Methods). Standard local electrical characterization was first employed to ensure the quality of the metal contacts, and that the graphene samples exhibited characteristic transport properties unique to graphene, as the half-integer quantum Hall effect (see Supplementary Notes 1 and 2). Figure 1a shows a sketch of a Hall bar depicting the measurement configuration of the nonlocal transport, and the two sources of nonlocal signals mainly discussed in this paper. Figure 1b shows an optical microscope picture of the macroscale sample. Figure 1 | Device structure and nonlocal measurements. a) Sketch of the Hall-bar-shaped channel for nonlocal measurements. The nonlocal measurement consists of driving a current, I, (black arrow) across contacts 1 and 2, and measuring the resulting voltage signal between probes 3 and 4. The nonlocal resistance is defined as 𝑅 𝑁𝐿 = 𝑉 𝐼 ⁄ . In the spin Hall regime, a spin current transverse to the injected charge current is generated via SHE (spin “down”, blue arrow; spin “up, red arrow), and diffuses through the channel of length L until it reaches electrodes 3 and 4, where via ISHE the spin current generates a charge current. The black dashed lines represent the van der Pauw contribution to the nonlocal resistance detected between probes 3 and 4. This background signal originates from the current spreading from terminals 1 and 2, and depends only on the resistivity of the conducting channel multiplied by an exponential decaying geometrical factor. In our devices, the current is always injected in the left arm, and detected at the right arm, with the external magnetic field applied perpendicular to the surface. b) Optical microscope image of the macroscale CVD graphene sample. The channel width, W, is 500 μm, with the electrodes distance L being 500, 700, 900, 1200 and 1500 μm. The width of the contacts is defined to be 1/10 th of the width of the channel. White scale bar is 500 μm. In a concise manner, nonlocal measurements in the context of spin Hall currents consist of driving a current in one arm of an H-bar shaped channel and detecting the voltage drop at the other arm (see Fig. 1a). Using the terminal notation shown in Fig. 1a, the nonlocal resistance is determined as 𝑅 NL = 𝑉 𝐼 ⁄ . If the channel length matches the spin diffusion length in the non-magnetic material, this simple device scheme and measurement setup makes possible the detection of spin-related signals in the spin Hall regime. The Hall-bar geometry allows for the generation of a transverse spin current from the charge current via SHE, which will then diffuse across the channel and be converted back into a charge current via ISHE. The magnitude of the effect will depend on the efficiency of the spin-to-charge conversion, and on the spin relaxation length. These two quantities can be determined from a transmission line method measurement of the nonlocal resistance , where the nonlocal resistance is related to the length, L , and width, W , of the transporting channel as 𝑅 NL = 𝜃 SH2 𝜌 𝑥𝑥 𝑊𝜆 s exp (− 𝐿 𝜆 s ⁄ ) , where 𝜃 SH is the spin Hall angle, 𝜌 𝑥𝑥 the sheet resistance, and 𝜆 s the spin diffusion length. In graphene, 𝑅 NL in the spin Hall regime is greatly enhanced at the Dirac point, requiring the use of a gate voltage to sweep the Fermi level of graphene to the charge neutrality point (CNP) . Although nonlocal measurements are used to probe non-charge related transport, there are classical, charge-related sources of nonlocality that can contribute to the signal detected between terminals 3 and 4 (see Fig. 1a). A robust source of nonlocal signals is the classical van der Pauw current spreading from the injecting terminals . By injecting a current on the left arm, a net current will reach the detection terminals, with magnitude decreasing exponentially with the distance to the injecting terminals. This Ohmic contribution in the device scheme resented is determined using the formula 𝑅 NL,Ohmic = 𝜌 𝑥𝑥 𝜋 ln [ cosh(𝜋𝐿 𝑊⁄ )+1cosh(𝜋𝐿 𝑊⁄ )−1 ] , which, for cases where L > W, is usually approximated as 𝑅 NL,Ohmic ≈ 𝜌 𝑥𝑥 exp (−𝜋 𝐿𝑊 ) . The sheet resistance is determined by performing a four-probe measurement of the respective channel, injecting a current between electrodes S and D, and measuring the voltage drop between electrodes 2 and 4, or 1 and 3. By comparing the detected nonlocal resistance versus the expected nonlocal Ohmic contribution, we can evaluate the emergence of signals not related to this classical source. In our study, the nonlocal resistance at the CNP of graphene is measured for different applied magnetic fields ( B ), and then compared to the Ohmic contribution evaluated under the same conditions. Throughout the manuscript, we keep the same relative orientation of the perpendicular B and of the arms injecting the current and detecting the voltage signal. We explored CVD graphene Hall bars with channel 5, 50, and 500 µm wide, with center-to-center distance between terminals maintaining similar aspect ratios ( 𝐿 𝑊⁄ ) of typically 1, 1.4, 1.8, 2.4 and 3.2. The width of the terminals is 1/10 th or less of the channel width (see Methods). We report mainly on the results obtained for the extreme cases of 5 and 500-µm-wide samples, using the data of the 50-µm-wide sample to extend the discussion. Further details on the fabrication and electrical measurements are provided in Methods. Nonlocality in macroscale devices.
Figure 2 shows the results obtained for the 500-µm-wide macroscale CVD graphene sample, for a channel length of 1500 μm, at cryogenic temperatures.
Figure 2 | Nonlocality in macroscale CVD graphene devices at low temperatures. a) Nonlocal resistance as a function of back-gate voltage for a 500-μm-wide and 1500-μm-long Hall bar. b) Nonlocal resistance as a function of the applied out-of-plane magnetic field, at the CNP of graphene, for different emperatures. c) Transmission line plot of the nonlocal resistance with channel length, at the CNP. In all figures, dashed lines represent the determined Ohmic contributions from four-probe measurements under equivalent conditions.
In Fig. 2a, a sweep of the gate voltage, V G , in the absence of magnetic field reveals the carrier density dependence of graphene, with CNP at V G = 34 V. The Ohmic background closely follows the nonlocal resistance and, at the CNP, matches the detected nonlocal resistance. When B is applied, a significant increase of R NL at the Dirac point occurs, together with a narrowing of the curve (see Supplementary Note 3 for the temperature dependence). This effect is not reproduced by the Ohmic contribution. By fixing the gate voltage at the CNP and sweeping B (Fig. 2b), one verifies two fundamental features. Firstly, the nonlocal magnetoresistance is asymmetric, with significantly higher nonlocal values for negative values of B . Secondly, for the side with higher signals, the difference between the nonlocal resistance and the Ohmic contribution increases with increasing B . Furthermore, the Ohmic signal shows a more symmetric magnetoresistance. Expanding the study to the different channel lengths and focusing on the side with highest nonlocal signal, Fig. 2c shows the dependence of R NL and of the Ohmic contribution at 2 K as a function of the channel length, at the CNP, with ( B = -6 T) and without ( B = 0) applied magnetic field. Without magnetic field, the expected Ohmic contribution and the measured R NL coincide, following the same dependence with the channel length. Fitting the dependence to an exponential decay, the resulting decay 𝜆 = 159.4 ± 1.6 µm closely matches the expected 𝑊 𝜋⁄ from the Ohmic term exp (−𝜋 𝐿𝑊 ) , calculated to be ~159.2 μ m. Upon switching on the magnetic field, R NL follows a similar exponential decay with channel length, but the magnitude of the resistance is greatly enhanced, in some cases being one order of magnitude higher than the Ohmic signal for similar conditions. If the signal were to be interpreted in terms of charge neutral spin currents, (or, equivalently, valley currents), a fitting to the expression 𝑅 NL = 𝜃 SH2 𝜌 𝑥𝑥 𝑊𝜆 s exp (− 𝐿 𝜆 s ⁄ ) would yield a spin relaxation length of 𝜆 s = 163 ± 19 μm and 𝜃 SH = 1.6 ± 0.4. Such values are clearly unreasonably high for such a disordered polycrystalline graphene sample supported onto an oxide substrate. Actually, the spin Hall angle determined is one order of magnitude higher than what is typical of metals with strong SOC, such as platinum ( 𝜃 SH ~ 0.2 ) . Nonlocality in microscale devices.
Figure 3a shows the nonlocal resistance measured for the 5- μ m-wide microscale CVD graphene sample, for a channel length of 16 μ m, at cryogenic temperatures. Figure 3 | Nonlocality in microscale CVD graphene devices at low temperatures. a) Nonlocal resistance as a function of back-gate voltage for a 5-μm-wide and 16-μm-long Hall bar. b) Nonlocal resistance as a function of the applied out-of-plane magnetic field, at the CNP of graphene, for different temperatures. c) Transmission line plot of the nonlocal resistance with channel length, at the CNP. In all figures, dashed lines represent the determined Ohmic contributions from four-probe measurements under equivalent conditions.
Similarly to what was found before, in the absence of an applied magnetic field ( B = 0), and by sweeping the gate voltage, R NL is well described by the Ohmic background, with a decay 𝜆 = 1.58 ± 0.03 µm. However, in this case, under an applied magnetic field ( B = 6 T), the enhancement of the nonlocal signal is one order of magnitude higher than in the macroscale sample for an equivalent aspect ratio, 𝐿 𝑊⁄ = 3, with R NL two orders of magnitude higher than the Ohmic background (see Supplementary Note 3 for the temperature dependence). Strikingly, fixing the gate voltage at the CNP and sweeping B (Fig. 3b), the magnetoresistance of the graphene becomes greatly asymmetric, but this time the higher nonlocal values happen for positive values of B . A strong asymmetry of the nonlocal resistance with B was also reported in other studies of nonlocality, but such effects were not analyzed . Evaluating R NL at the CNP for the different channel lengths at 2 K with ( B = 6 T) and without ( B = 0) applied magnetic field, Fig. 3c exhibits a similar trend to that reported for the macroscale sample, with a clear agreement between the nonlocal resistance and the Ohmic contribution at B = 0, following the same dependence with the channel length. In the presence of B , the nonlocal resistance dependence is also clearly described by a similar exponential decay with channel length, but the magnitude of the signals is even further enhanced, from one to two orders of magnitude larger than the expected Ohmic background. An analysis in terms of non-charge currents yields 𝜆 s = .7 ± 1.2 µm, and 𝜃 SH = 1.3 ± 0.6. We note that this relaxation length is similar to reports that place 𝜆 s between 1 and 5 μ m in CVD graphene on SiO . Discussion
So far, we have demonstrated that CVD graphene samples that differ in dimensions by two orders of magnitude can show similar large nonlocal signals close to the Dirac point, which dominate over the conventional Ohmic contribution. In the absence of magnetic field, these samples exhibit nonlocal resistance in perfect agreement with what is expected from the van der Pauw background. The source of the strong nonlocal signal when the magnetic field is applied is less clear and more complex. In the macroscale devices, an origin related to spin diffusion within the spin Hall regime would require the spins generated in the injecting terminal to survive a disordered 1.6- mm -long channel, and convert back into a charge current, which disagrees with all the estimations of spin diffusion lengths reported to date in CVD graphene samples . This strongly suggests that the signal visible in the macroscale sample is not related to a pure spin transport mechanism within a SHE and ISHE process. At the microscale sample, the signal increases even further, by one order of magnitude. Again, the dependence with channel length indicates that the signal decays exponentially similarly to that predicted from the Ohmic contribution of 1.59 μ m. Besides the origins related to spin transport and van der Pauw backgrounds, thermoelectric effects have also been proposed as a possible mechanism driving nonlocality at the CNP of graphene, in particular the Ettingshausen-Nernst effect . In this picture, the nonlocal signal would be generated in two steps: first, under an applied perpendicular magnetic field, the current being driven would lead to a transverse heat flow (Ettingshausen effect); second, the thermal gradient across the detection terminal under B would generate the nonlocal voltage (Nernst effect). Since thermoelectric coefficients in graphene show a strong increase around the Dirac point, the gate dependence should manifest an enhancement at the CNP. But considering both Ettingshausen and Nernst effects, the dependence of 𝑅 NL with B should be quadratic, 𝑅 NL ∝ 𝐵 . This seems not to be the case in our samples, neither in the microscale sample (Fig. 3b) nor in the macroscale one (Fig. 2b), where the dependence of the nonlocal resistance with magnetic field is strongly asymmetric. From this discussion, the nonlocal signals in our experiment do not seem to originate from spin or valley Hall effects, thermoelectric effects, or purely Ohmic contributions. To further clarify the nature of the observed features, we repeated the experiments for a sample with dimensions in between 5 and 500 μ m. With a 50- μ m-wide channel, and equivalent aspect ratios between 1 and 3.2, the sample is on the macroscale side, and should exclude spin transport. Figure 4a summarizes the nonlocal resistance at the CNP of graphene taken at | B | = 6 T as a function of the channel aspect ratio, for all three samples (5, 50 and 500- μ m-wide channel). Figure 4 | Origin of nonlocal resistance for the CVD graphene samples. a) Nonlocal resistance as a function of aspect ratio of the Hall bars for the 5-, 50-, and 500-μm-wide samples, at CNP and 2 K. b) Optical image of the microscale CVD graphene sample during the fabrication marking the visible grain boundaries and bilayer islands. White scale bar is 5 μm. c) At the CNP, and under an external field, the counter-propagating edge currents are coupled to the charge of the carrier via Zeeman interaction. The presence of hole/electron asymmetric grain boundaries shunting the edge states through the insulating bulk leads to nonlocal measurements that are sensitive to the sample topology and sign of the magnetic field.
The nonlocal resistance for the 50- μ m-wide sample remains between those obtained for the two other samples, with signals 1 to 2 orders of magnitude higher than expected from the van der Pauw background. Interestingly, the signal magnitude is similar to the 5- μ m-wide sample, where spin signals could conceivably survive a 16- μ m-long channel. In this case, with a 160- μ m-long channel, we exclude a spin transport origin. The mechanism driving the nonlocal response must be able to generate signals dependent with the magnitude and sign of the external magnetic field, seemingly invariant with scale, and that follow a dependence with channel length similar to the Ohmic background. To explain all these anomalous features, we propose an interpretation based on counter-propagating edge states shunted by grain boundaries (GBs). When an out-of-plane magnetic field is applied to graphene, the bulk gradually becomes insulating, while charge current increasingly flows at the edges, resulting in a well-defined quantum Hall effect for large fields (see Supplementary Note 2). The uniqueness of the electronic properties of graphene give rise to a zeroth Landau level (LL0) at the CNP which is populated by both electrons and holes . The Zeeman splitting of the LL0 couples the charge carrier type to the spin state, so that electrons and holes are now associated with opposite spin olarizations, which triggers a dissipative quantum Hall effect near the CNP, driven by spin-polarized counter-propagating edge states (as shown in Fig. 4c) . A key feature of this effect is that the longitudinal resistivity 𝜌 𝑥𝑥 at the CNP is now predominantly driven by edge, and not bulk, channels . Under these conditions, and considering a 4-point contact configuration similar to that explored in our work, a nonlocal signal emerges at the CNP solely driven by potential drops at the contacts (see Supplementary Note 4 for a multi-terminal Landauer-Büttiker formulation of the device detailing the origin of the nonlocality). A similar fingerprint (peak in R NL at the CNP) has already been found in a nonlocal measurement on two-dimensional system with a simultaneous presence of electrons and holes in a 20 nm HgTe quantum well . In our experimental data, the presence of a strong asymmetric behavior of R NL with the direction of B is a salient observation that demands interpretation. The strong nonlocal signal can be almost entirely quenched by changing the sign of B . The quenching is consistent with bulk conduction, meaning that a shunting of the edge currents that drive the nonlocal signal enter into play. Indeed, previous studies associate conducting bulk states with a decay of the 𝜌 𝑥𝑥 peak away from the CNP . In CVD graphene, line defects (i.e., GBs) are the most likely source of this conduction through the bulk, and as shown in Fig. 4b, our nonlocal transport geometry intercepts many GBs along the transport channel. Importantly, such defects can display strong electron-hole asymmetry in their transmission properties due to local sublattice symmetry breaking or by intercepting electron-hole puddles across the device. This manifests as an asymmetry with respect to magnetic field direction due to the Zeeman splitting, which associates each charge carrier type with a different spin orientation. In Supplementary Note 4, we explicitly demonstrate how nonlocal resistances develop an asymmetry with respect to the sign of an applied magnetic field in the presence of GBs. It is important to note that in the context of quantum Hall experiments a strong dependence of 𝜌 𝑥𝑥 with the sign of B or a signature of edge transport in nonlocal magnetotransport measurements does not necessarily require a graphene-specific model. In particular, an asymmetry of 𝜌 𝑥𝑥 with B has been observed in well-behaved semiconductor two-dimensional electron gases , with mechanisms such as conduction paths in the bulk , carrier density gradients , hybrid constrictions , and in-plane electric fields driving the effect. Our proposed mechanism, however, differs in that it can explain not only the asymmetry of the nonlocal resistance, but also the large nonlocal signals beyond the Ohmic contribution at the CNP. The physics of the LL0 in graphene, with a coexistence of electrons and holes (and thus a sensitivity to electron-hole asymmetries), has no equivalent in other low-dimensional systems. To demonstrate that GBs are at the origin of the B asymmetry of the nonlocal signal, we perform a tight-binding simulation in the nonlocal geometry with a single 558-type line defect connecting the edges (see Methods) of a 40-nm-wide and 85-nm-long graphene channel. As seen in Fig. 5a, a highly asymmetric nonlocal peak at the CNP emerges when a Zeeman splitting is introduced in the presence of a line defect. This follows directly from the electron-hole asymmetry of transmission through the defect, discussed in further detail in Supplementary Notes 4 and 5. The similarity between the modeled system and the experimental results shown in Fig. 5b for the device with 𝐿 𝑊⁄ = 3.2 of the 5- μ m-wide sample strongly suggests counter-propagating edge states shunted by GBs to be at the root of the large, asymmetric nonlocality observed in the samples. An extended study over all devices with different aspect ratios of the 5- μ m-wide sample shows a strong asymmetry of the nonlocality with the sign of the magnetic field for all devices, further corroborating our analysis (see Supplementary Note 6). While all the previous discussed origins could not address this strong asymmetry, our proposed mechanism explains it. Furthermore, our simulations capture additional peaks associated with the first Landau level (LL1). A finite R NL is here associated ith conduction through bulk states and not to counter-propagating states, which only occur near the CNP. Higher LL peaks in R NL follow 𝜌 𝑥𝑥 , a feature previously demonstrated in other works . Figure 5 | Nonlocal resistance dependence with magnetic field direction in a graphene Hall bar device . a) Tight-binding simulation in the nonlocal geometry with a single 558-type line defect connecting the edges of a Hall-bar graphene channel. Nonlocal signals versus position of the Fermi level relative to the CNP for opposite directions of the applied perpendicular magnetic field. b) Nonlocal resistance as a function of back-gate voltage for a 5-μm-wide and 16-μm-long graphene Hall bar at 2 K for opposite directions of the applied perpendicular magnetic field.
Due to the random formation of inhomogeneities and GBs during the growth of CVD graphene, samples with individual, random line defects can show nonlocality preferentially for either positive or negative magnetic fields. Electron-hole asymmetry in GB transmissions can arise due to both intrinsic sublattice effects and external doping effects. Other considerations, including the positioning of GBs relative to the probes, the rate of backscattering between counter-propagating states, and the density of electron-hole puddles, can play an additional role in determining the exact nonlocal signature. However, the qualitative effect (a nonlocal peak at the CNP whose bulk-mediated suppression is dependent on the sign of B ) is very general and, in principle, independent of sample size. However, larger samples contain more grains across their width and the asymmetric effects of individual GBs will tend to be averaged out. This will lead to smaller asymmetries with the direction of B , but also smaller nonlocal signals due to a greater number of bulk conduction channels. This picture fits convincingly with the range of effects exhibited by our measured graphene samples from the microscale to the macroscale. Importantly, the appearance of a strong nonlocal signal is independent of spin or valley transport mechanisms arising from, for example, long-ranged coherent spin transport and the SHE/ISHE effects. In summary, we identify strong nonlocal signals in both microscale and macroscale CVD graphene Hall-bar devices, emerging when a perpendicular magnetic field is applied. The observed signals share many similarities to experiments relying on spin Hall or valley Hall effect mechanisms, but our control experiments at the macroscale rule out both long range spin and valley polarized transport. The similarity of the nonlocal phenomenon across different scales suggests that a different mechanism is at the origin of such nonlocality. We propose a mechanism driven by field-induced spin-split edge states. The sensitivity of these states to electron-hole asymmetric transport, induced by features such as GBs, strongly influences the nonlocal resistance profiles. The persistence of this behavior to millimeter length scales shows that defect-induced contributions to nonlocal transport could emerge in a wide range of easurements, and may indeed dominate over more exotic sources of nonlocality in practical devices made from scalable graphene materials. Methods
Device fabrication.
The same fabrication procedures were followed for the fabrication of the 5-, 50- and 500- μ m-wide CVD graphene samples. 1×1 cm Si(n++ doped)/SiO (300nm) chips with monolayer CVD graphene grown on copper foils were acquired from a commercially available supplier, Graphenea S.A. All samples reported here come from the same batch. The samples where first spin coated with double-layer of Poly(methyl methacrylate) (PMMA) (495/950kDa) and then baked at 180 ºC for 90 s. To define the Hall-bar-shaped graphene, the areas surrounding the Hall bar were exposed with electron-beam lithography, and etched with a chemistry of Ar/O (80/5 sccm) in a capacitive coupled plasma reactive-ion etching (RIE) Oxford PlasmaLab 80 equipment. The remaining resist leftovers were striped in acetone bath at room temperature for 4 hours, immersed in isopropanol, and dried with a nitrogen gun. For the definition of the electrical contacts, the samples were again spin coated with double layer PMMA, and the metal electrodes defined over the Hall bar. The metallization was done on ultra-high vacuum at a base pressure of 10 -9 mbar, using electron-beam deposition of Ti (5 nm)/Au (40 nm) at a rate of 0.5 and 1.5 Ås -1 . The lift-off was performed with acetone at room temperature. One important feature of the device design was the terminals width of 1/10 th or less of the channel width. This comes from a straightforward analysis of the van der Pauw expression for charge diffusive backgrounds, where if the width of the contact is on the same order of the channel width it can lead to edge-to-edge signals ~20 times in magnitude larger than the center-to-center signal. Electrical characterization.
The devices were characterized in a Quantum Design physical property measurement system (PPMS) using standard four-probe direct current measurement methods. The measurements were performed using a Keithley 2182A as current source and a Keithley 6221 as voltmeter. The gate voltage was applied using one channel of the Keithley 2636. In all measurements, the excitation current was 10 μ A. Before measuring the devices, we performed an in-situ annealing at 400 K for three hours, with helium flushing cycles to release the sample chamber of evaporated water. The samples would then be cooled down at the maximum rate to 2 K, and any temperature dependent study performed for increasing temperatures.
Tight-binding simulation.
The system is described by a standard nearest-neighbor tight-binding model with the magnetic field incorporated using the Peierl’s phase approach
𝐻 = ∑ 𝑡 𝑒 ∫ 𝐴 𝑟𝑗𝑟𝑖 (𝑟)⋅𝑑𝑟 𝑐^ 𝑖† 𝑐^ 𝑗 , where 𝐵 = 𝛻 × 𝐴 . The transmissions between each set of probes is given by the Caroli formula 𝑇 𝑝𝑞 =𝑇𝑟[𝐺 𝑅 𝛤 𝑞 𝐺 𝐴 𝛤 𝑝 ] with the required Green’s functions ( 𝐺 𝑅 , 𝐺 𝐴 ) calculated using efficient recursive techniques. The left and right leads are included via broadening terms ( 𝛤 ) calculated from the surface Green’s function of semi-infinite nanoribbons, whereas constant broadening terms, representing metallic contacts, are chosen for the top and bottom probes. The Zeeman term is included separately via equal and opposite energy shifts of the spinless transmissions by half the required splitting (10 -3 eV). The total transmission is then the sum of the two spin contributions (our model excludes spin-mixing terms). The nonlocal resistance is calculated from the potentials 𝑉 𝑝 and currents 𝐼 𝑝 at each probe, which emerge from solving the multi-terminal Landauer-Büttiker relation 𝐼 𝑝 = ∑ (𝑇 𝑞𝑝 𝑉 𝑝 − 𝑇 𝑝𝑞 𝑉 𝑞 ) 𝑞 with suitable boundary conditions. Additional potential terms are included to represent the presence of charge inhomogeneities (electron-hole puddles) in the system induced by the graphene/substrate interaction, and the random nature of experimental grain boundaries is accounted for by random local potentials in the vicinity of the grain boundary. Additional details and simulations are presented in Supplementary Note 5. Data availability . All relevant data are available from the authors.
References
1. Tombros, N., Jozsa, C., Popinciuc, M., Jonkman, H. T. & van Wees, B. J. Electronic spin transport and spin precession in single graphene layers at room temperature.
Nature
Nat. Phys. et al. Giant Nonlocality Near the Dirac Point in Graphene.
Science et al.
A two-dimensional spin field-effect switch.
Nat. Commun. et al. Large room temperature spin-to-charge conversion signals in a few-layer graphene/Pt lateral heterostructure.
Nat. Commun.
661 (2017). 6. Schaibley, J. R. et al.
Valleytronics in 2D materials.
Nat. Rev. Mater. et al. Toward Wafer Scale Fabrication of Graphene Based Spin Valve Devices.
Nano Lett. et al.
Giant spin Hall effect in graphene grown by chemical vapour deposition.
Nat. Commun. et al. Spin Hall Effect and Origins of Nonlocal Resistance in Adatom-Decorated Graphene.
Phys. Rev. Lett. et al.
Spin–orbit proximity effect in graphene.
Nat. Commun. et al. Detecting topological currents in graphene superlattices.
Science et al.
Edge currents shunt the insulating bulk in gapped graphene.
Nat. Commun. Phys. Rev. Lett. et al.
Strong interface-induced spin–orbit interaction in graphene on WS . Nat. Commun. et al. Strong interfacial exchange field in the graphene/EuS heterostructure.
Nat. Mater.
Phys. Rev. B
Phys. Rev. Lett.
Phys. Rev. B
Phys. E
Phys. Rev. Lett. et al.
Proximity Effects Induced in Graphene by Magnetic Insulators: First-Principles Calculations on Spin Filtering and Exchange-Splitting Gaps.
Phys. Rev. Lett.
2D Mater. Phys. Rev. B
Phys. Rev. B
Phys. Rev. B
Riv. del Nuovo Cim.
Phys. Rev. Lett. et al.
Dissipative Quantum Hall Effect in Graphene near the Dirac Point.
Phys. Rev. Lett.
Nat. Nanotechnol. Nat. Mater. J. Phys. Condens. Matter et al.
Charge Transport in Polycrystalline Graphene: Challenges and Opportunities.
Adv. Mater. et al.
Two-Dimensional Graphene with Structural Defects: Elastic Mean Free Path, Minimum Conductivity, and Anderson Transition.
Phys. Rev. Lett. et al.
Quantum Hall effect on centimeter scale chemical vapor deposited graphene films.
Appl. Phys. Lett.
Nature
Philips Res. Reports
Philips Tech. Rev. et al.
Tuning the spin Hall effect of Pt from the moderately dirty to the superclean regime.
Phys. Rev. B
Nat. Commun. et al. Extremely large magnetoresistance in few-layer graphene/boron–nitride heterostructures.
Nat. Commun. et al. Room-Temperature Quantum Hall Effect in Graphene.
Science
Phys. Rev. B et al.
Nonlocal Transport Near Charge Neutrality Point in a Two-Dimensional Electron-Hole System.
Phys. Rev. Lett.
Phys. Rev. B
Phys. Rev. B
Phys. Rev. Lett. et al.
Quantum Hall effect in a high-mobility two-dimensional electron gas on the surface of a cylinder.
Phys. Rev. B
Phys. Rev. B et al.
The effect of carrier density gradients on magnetotransport data measured in Hall bar geometry.
Solid State Commun.
Appl. Phys. Lett.
Phys. Rev. B
Phys. E Low-dimensional Syst. Nanostructures et al.
New resistivity for high-mobility quantum Hall conductors.
Phys. Rev. Lett. et al.
Nonlocal Shubnikov-de Haas oscillations through edge and bulk currents in GaAs/AlGaAs mesoscopic quantum wires.
Solid State Commun. et al.
Role of edge and bulk currents through a gate barrier in nonlocal resistance of GaAs/Al x Ga As.
Phys. Rev. B
Acknowledgements
This work was supported by the European Union 7th Framework Programme under the Marie Curie Actions (607904-13-SPINOGRAPH), by the Spanish MINECO (Project No. MAT2015-65159-R) and by the Regional Council of Gipuzkoa (Project No. 100/16). CIC nanoGUNE is supported by the Spanish MINECO under the Maria de Maeztu Units of Excellence Programme (MDM-2016-0618). S.R. acknowledges funding from Spanish MINECO (Project No. MAT2012-33911), the European Regional Development Fund (Project No. FIS2015-67767-P (MINECO/FEDER)), and the Catalan Government (Secretaria d’Universitats i Investigació del Departament d’Economia i Coneixement). S.P. and S.R. acknowledge funding from the European Union Horizon 2020 Programme under the Marie Sklodowska-Curie grant agreement No. 665919 and No. 696656 Graphene Flagship. ICN2 is supported by the Spanish MINECO under the Severo Ochoa Centres of Excellence Programme (SEV-2013-0295) and by the Catalan Government under the CERCA Programme.
Author Contributions
M.R. and F.C. designed the experiments. M.R. fabricated the samples, performed the measurements, and did the data analysis. S.R.P. did the tight-binding simulations. M.R. wrote the paper. All authors participated in the discussion of the results and contributed to the writing of the manuscript. L.E.H. and F.C. supervised the entire study.
Additional Information
Supplementary information
Competing financial interests:
The authors declare no competing financial interests upplementary Information
Supplementary Note 1: Preliminary electrical characterization of chemical vapour deposition graphene devices.
Standard electrical characterization was performed in the different chemical vapour deposition (CVD) graphene Hall bars. Supplementary Figure 1 shows I - V curves, transfer curves, temperature dependence of the sheet resistance ( 𝜌 𝑥𝑥 ), and the temperature dependence of the carrier density ( 𝑛 ) obtained from Hall measurements for the 500- μ m-wide sample. Supplementary Table 1 shows the mobility determined from the expression 𝜇 = 𝑥𝑥 𝑛 𝑒 , where 𝜌 𝑥𝑥 is taken for a doping of cm -2 , and 𝜌 𝑥𝑥 at the charge neutrality point (CNP) for the samples shown in the main manuscript. Supplementary Figure 6 | Typical electrical performance of the devices fabricated. a) Two-point output characteristics of the 500-μm-wide CVD graphene sample contacted with Ti (5 nm)/Au (40 nm) showing linear
I-V curves. b) Sheet resistance as a function of gate voltage (transfer characteristics) of the 500-μm-wide CVD graphene device with
L/W = 3 after the in-situ annealing at 400 K with helium flushes, at several temperatures. Thickness of the SiO dielectric is 300 nm. c) Temperature dependence of the sheet resistance of the 500-μm-wide CVD graphene device with L/W = 3 at V G = 0 V, with weak-localization effects emerging for the lowest temperatures. d) Temperature dependence of the carrier density determined from Hall measurements of the 500-μm-wide CVD graphene device with L/W = 3 at V G = 0 V. Supplementary Table 1 |
Sheet resistance at the CNP and mobility for the devices shown in the manuscript at a temperature of 2 K.
Channel width ( μ m) 500 50 5 μ (cm V -1 s -1 ) @ 2 K ρ xx @ 2 K @ CNP (Ω) upplementary Note 2: Magnetotransport measurements of chemical vapour deposition graphene devices. We performed magnetotransport measurements of the CVD graphene Hall bars in order to characterize their edge transport. Supplementary Figure 2 shows the half-integer quantum Hall effect (QHE) well developed for an external magnetic field of 9 T for the whole range of device dimensions explored in our study.
Supplementary Figure 2 | Half-integer quantum Hall Effect in the CVD graphene samples. a) Transverse conductivity 𝜎 𝑥𝑦 and sheet resistance 𝜌 𝑥𝑥 as a function of the gate voltage for the 500-μm-wide macroscale samples. b) 𝜎 𝑥𝑦 and 𝜌 𝑥𝑥 as a function of the gate voltage for the 5-μm-wide microscale samples. Dashed lines represent the quantized transverse conductivity according to 𝜎 𝑥𝑦 = ±4(𝑖 + 1/2) 𝑒 ℎ , with i = 0, 1, 2, … Supplementary Note 3: Temperature dependence of the nonlocality for CVD graphene devices.
While the study of the nonlocal resistance was mainly performed at the lowest temperatures available to our system ( T = 2 K), the effects are visible for the entire range of temperatures up to room temperature. Supplementary Figure 3 shows the nonlocal resistance as a function of the gate voltage for different temperatures between 2 K and 300 K, for both macroscale and microscale devices. Supplementary Figure 3 | Nonlocal resistance as a function of gate voltage under external magnetic field magnitude of 6 T for several temperatures. a) 500-μm-wide macroscale device with
L/W = 3. b) 5-μm-wide microscale device with
L/W = 3.2. In both cases, the shift in the CNP for the 300 K measurements is related to the adsorption of water (see Methods in the main manuscript). upplementary Note 4: Magnetic field asymmetry in nonlocal resistances
In this note, we provide a simple Landauer-Büttiker analysis of how nonlocal resistances can develop an asymmetry with respect to the sign of an applied magnetic field. We consider a four-probe device with three different device, grain boundary (GB) and backscattering configurations.
Quantum spin Hall limit.
Firstly, we consider the spin-split QHE induced by Zeeman splitting at the CNP in the absence of GBs or scattering. In this case, perfectly ballistic counter-propagating edge states allow unitary transmission between neighbouring probes and give rise to an effective quantum spin Hall effect (QSHE). This case, shown schematically in Supplementary Fig. 4, gives rise to a quantized nonlocal resistance. A nonlocal resistance arises solely due to potential drops at the probes, and the quantized R NL value is dependent on the device geometry and probe layout. Supplementary Figure 4| Sketch of ballistic, quantized, counter-propagating channels in Zeeman-split QHE.
Solid arrows show electron channels, colored by source probe, and dashed black arrows shows the hole channels. Arrow directions are inverted upon changing the sign of the external magnetic field B . At the CNP, up and down spins are carried by separate electron and hole channels. The transmission matrices for electrons and holes here are given by 𝑇 e = (0 0 0 11 0 0 00 1 0 00 0 1 0) and 𝑇 h = (0 1 0 00 0 1 00 0 0 11 0 0 0) respectively, where 𝑇 𝑝𝑞e describes transmission from probe q to probe p due to electrons. As the spin channels are independent, the total transmission is a simple sum of these terms, giving 𝑇 = (0 1 0 11 0 1 00 1 0 11 0 1 0) . The nonlocal resistance is determined by solving the multi-terminal Landauer formula 𝐼 𝑝 = 𝑒 ℎ ∑ (𝑇 𝑞𝑝 𝑉 𝑝 − 𝑇 𝑝𝑞 𝑉 𝑞 ) 𝑞 for the currents and potentials at each probe, and in this case is found to be 𝑅 NL =
14 ℎ𝑒 . Changing the sign of B changes the direction of current flow for both the electron and hole channels, and is equivalent to taking the transpose of 𝑇 . Since 𝑇 is symmetric, the nonlocal resistance is unchanged. Grain boundary with electron-hole asymmetric transmission.
We now consider the simple case where the GB has a transmission channel for electrons only. This is easily generalizable to he case where the GB transmits electrons and holes with different probabilities. We also assume that one end of the GB is located near one of the contacts (probe 3) so that direct transmission from the GB to the contact is possible – we will later explore the relaxation of this assumption. The setup, as shown in Supplementary Fig. 5, differs from the previous case only for the spin channel carried by electrons: e.g. a certain amount (‘ x ’) of electron current from probe 1 towards probe 2 is now deflected along the grain boundary, where it again splits with part (‘ y ’) travelling along the electron edge channel to probe 3, and the remainder (‘ x - y ’) coupling directly to probe 4. Electron current from probe 3 experiences a similar split. The transmission matrix here is 𝑇 = ( 0 1 0 11 − 𝑥 0 𝑥 + 1 0𝑥 − 𝑦 1 0 1𝑦 + 1 0 1 − 𝑦 0) , giving nonlocal resistances 𝑅 NL (𝐵 >0) = and 𝑅 NL (𝐵 < 0) = . Supplementary Figure 5 | Sketch for Zeeman-split QHE system with a grain boundary allowing transmission only for electrons.
As in Supplementary Fig. 4, but with a grain boundary connecting the two sides of the device and allowing electron transmission. The red and blue labels show the division of current in the electron channels emerging from leads 1 and 3 due to the grain boundary.
The difference between these values depends strongly on a direct coupling between the GB and one of the probes – as this is switched off (𝑦 → 𝑥) , the asymmetry vanishes despite changes in the individual potentials, and despite T not being a symmetric matrix. An asymmetric T is a necessary but not sufficient condition for a B -field asymmetry – in the QHE, T is asymmetric but the nonlocal resistance is zero for both charge carrier types. However, we also note that electron-hole asymmetry is essential to the appearance of a B -field asymmetry – equal transmission for elections and holes through the GB results in a symmetric T matrix where changing the sign of B merely swaps the roles of electrons and holes while conserving the total transmissions. So, even in the limit of a quasi-QSHE with perfectly ballistic edge channels, an asymmetric non-local response can emerge due to the presence of a GB supporting electron-hole asymmetric transmission and a direct coupling to one of the probes. In order to directly couple to a probe the separation between the GB and the probe should be on the order of the magnetic length 𝑙 𝐵 , which for the experimental case at 𝐵 = 6 T is ~10 nm. This length may be extended due to, e.g., electron-hole puddles which can open additional conducting channels near pn interfaces, but is most likely too small to account for all the experimental signatures, in particular for the largest systems. However, it is a good system to test electron-hole nonlocal asymmetries within a tight-binding framework, and we shall examine it further in Supplementary Note 5. Electron-hole asymmetric grain boundaries with backscattering.
We now move from the effective QSHE regime to consider the dissipative regime where the edge channels are no longer ballistic. As discussed in Abanin et al . (see Supplementary Ref. 1), this can occur by ackscattering along one edge or through bulk-mediated mechanisms. Backscattering between channels along the same edge requires mixing of the different spin channels, which can be mediated by spin-orbit coupling or magnetic disorder. In this situation, we consider a very simple example where backscattering reduces transmission along the top edge of the device by an amount a for both electrons and holes (see Supplementary Fig. 6). This can be easily generalised to cases with different backscattering rates for different edges and carrier types to account for, e.g., the distance dependence of transmissions through systems with scattering. The grain boundary once more transmits only for electrons, but this time does not have a direct coupling to any lead. We now have 𝑇 = ( 0 1 − 𝑎 0 11 − 𝑎 − 𝑥 0 𝑥 + 1 00 1 0 11 + 𝑥 0 1 − 𝑥 0) , 𝑅 NL (𝐵 > 0) = and 𝑅 NL (𝐵 < 0) = (1−𝑎)(1−𝑥)4−𝑎(𝑥+3) ℎ𝑒 . As before, the asymmetry between B and -B only arises when an electron-hole asymmetric channel opens between the two sides of the device. Supplementary Figure 6 | Sketch for Zeeman-split QHE system with backscattered channels and an electron-hole asymmetric grain boundary.
As in Supplementary Fig. 5, but now the grain boundary has no direct connection to any probe. Electrons and hole channels between probes 1 and 2 are now dissipative, with backscattering (e.g. the yellow region) reducing the total current in these channels to . Supplementary Note 5: Tight-binding calculations
The setup for the tight-binding simulation is shown in Supplementary Fig. 7a. We consider a ~ 40-nm-wide 85-nm-long armchair nanoribbon with four metallic probes attached as shown. A 558-grain boundary (GB) is placed between the current and voltage probes, and we displace probe 3 so that it can couple directly with the GB as discussed in the previous section. To model a disordered GB as expected in experiment, Anderson disorder is added in a narrow channel surrounding the GB (shown by red and blue circles in Supplementary Fig. 7a for positive and negative potentials, respectively). In addition, Gaussian-type electron-hole puddles are included throughout the device to simulate the effects of charged defects in the substrate. An example of the net potential profile in the system is shown by the color map in Supplementary Fig. 7b. The interprobe transmissions are calculated using the Landauer formula with recursive Green’s functions techniques. From these the nonlocal resistance is calculated. The two effects of an external, perpendicular magnetic field are considered separately within our simulations. Firstly, the onset of a stable quantum Hall regime is controlled by including Peierl’s phase factors in the tight-binding hopping parameters. Since our simulated systems have significantly narrower width than the experimental cases, we chose B = ±60 T which ensures that the edge states on either side of the device are sufficiently decoupled. A Zeeman splitting is then manually introduced to break the degeneracy between up- and down-spin states. We choose ΔZ ~ 10 -3 eV, agreeing with the expected values, including interaction effects, for an experimental 𝐵 ~ 10
T (see Supplementary Ref. 1). This is applied as a rigid shift of the two spin bands of ΔZ/2. Similar effects to those discussed in the main paper, and below, can be achieved with a variety of different B and ΔZ values. To test our simulation approach, we calculate the nonlocal resistance for a pristine ribbon without a GB or potential disorder. A nonlocal resistance peak at the CNP is observed with a height slightly below the value 𝑅 NL =
14 ℎ𝑒 expected in the QSHE regime. Supplementary Figure 7c shows this system with an exaggerated ΔZ to highlight that the peak width agrees exactly with the Zeeman splitting. The minor deviation from the quantized value is due to quantum tunneling processes through the bulk of the narrow ribbon simulated – this deviation decreases on increasing the magnetic field or the Zeeman splitting as expected. Furthermore, the curves for positive and negative B coincide for this system. Additional peak features at 𝐸 F =±0.09|𝑡| are due to the bulk states of the 1 st Landau level (LL1), which are broadened and split due to the Zeeman term. Features occur at this position for all our simulated systems, and at the expected LL1 positions in the experimental results for the microscale devices. We note that these peaks emerge due to a nonlocal signal driven by bulk states, unlike the edge state processes determining the signal at the CNP.
Supplementary Figure 7 | Tight-binding calculations . a) Sketch of the 4-probe device setup for the tight-binding simulation. A current is driven between probes 4 and 1, and the nonlocal resistance is measured between probes 3 and 2. A disordered 558-type grain boundary connects the two sides of the device. b) An example of the potential profile due to both electron-hole puddles and local Anderson disorder on the grain boundary sites. c-f) Examples of 𝑅 𝑁𝐿 calculated for different systems, with the width of the Zeeman splitting in each case shown by dashed lines. c) shows a pristine ribbon device, d) introduces a perfect 558-GB without Anderson disorder, but with electron-hole puddles. e) and f) show two examples with Anderson-disordered GBs and electron-hole puddles. In all cases, B is perpendicular to the sample plane and has magnitude of 60 T. Supplementary Figure 7d shows the dramatic effect that a highly conducting GB has on the nonlocal signal. The pristine GB introduces a direct channel between the current probes, hunting the edge state channels which underpin the 𝑅 NL peak at the CNP. The resulting signal is highly asymmetric for different signs of B and oscillates around zero. The range of non-zero values is broadened by the introduction of electron-hole puddles. Realistic CVD devices are unlikely to contain long regions of atomically pristine GBs, so we add a strong Anderson-type disorder of strength 3 t to atomic sites within 2 lattice constants of the grain boundary. This reduces the transmission efficiency of the bulk leakage channel, and allows us to tune the relative contributions of the counter-propagating edge states and grain boundary channels. Two examples with different background puddle configurations are shown in panels in Supplementary Fig. 7e and Supplementary Fig. 7f. These simulations replicate the main features of the experimental data – a large peak near the CNP for one sign of B and an almost complete suppression of the peak for the opposite sign. Furthermore, we note that the experimental features observed at LL1 in the microscale devices are also observed here. Comparing the two examples, we note that the sign of B displaying the peak depends on the underlying electron-hole puddle distribution. The distribution corresponding to Supplementary Fig. 7e is shown in Supplementary Fig. 7b, where we can expect that the transmission of the GB will be strongly affected by the two hole-doped regions through which is passes. Indeed, very large puddles capable of spanning the device can set up leakage channels through the formation of snake states along pn boundaries. However, this is unlikely in the experimental setup due to a device width much larger than the puddle length scale (~10 nm). upplementary Note 6: Dependence of the nonlocality with the sign of the magnetic field for different aspect ratios and sample widths. The strong asymmetry of the nonlocal signal with the sign of the magnetic field was present for all the devices considered in the study. For the 5- μ m-wide microscale sample, the suppression of the nonlocal signal at the CNP is more evident. Supplementary Figure 8 shows the nonlocal resistance as a function of the gate voltage for all the devices in the 5- μ m-wide microscale sample, and for one device in the 500- μ m-wide macroscale sample. Supplementary Figure 8 | Nonlocal resistance as a function of the gate voltage for opposite signs of an external magnetic field with magnitude of 6 T for several aspect ratios and channel widths.
Black solid lines: B = +6 T, Red solid lines: B = - 6 T. Panels a), b), c), d), and e) correspond to the 5- μ m-wide sample, with L/W = 1, 1.4, 1.8, 2.4, and 3.2, respectively. Panel f) corresponds to the 500- μ m-wide sample, with L/W = 2.4.
Supplementary References
1. Abanin, D. A., Lee, P. A. & Levitov, L. S. Spin-Filtered Edge States and Quantum Hall Effect in Graphene.
Phys. Rev. Lett.96,