1D ferromagnetic edge contacts to 2D graphene/h-BN heterostructures

We report the fabrication of one-dimensional (1D) ferromagnetic edge contacts to two-dimensional (2D) graphene/h-BN heterostructures. While aiming to study spin injection/detection with 1D edge contacts, a spurious magnetoresistance signal was observed, which is found to originate from the local Hall effect in graphene due to fringe fields from ferromagnetic edge contacts and in the presence of charge current spreading in the nonlocal measurement configuration. Such behavior has been confirmed by the absence of a Hanle signal and gate-dependent magnetoresistance measurements that reveal a change in sign of the signal for the electron- and hole-doped regimes, which is in contrast to the expected behavior of the spin signal. Calculations show that the contact-induced fringe fields are typically on the order of hundreds of mT, but can be reduced below 100 mT with careful optimization of the contact geometry. There may be an additional contribution from magnetoresistance effects due to tunneling anisotropy in the contacts, which needs further investigation. These studies are useful for optimization of spin injection and detection in 2D material heterostructures through 1D edge contacts.


Introduction
Graphene has been shown to be a promising material for spin-polarized electron transport due to its low spin-orbit coupling and negligible hyperfine interactions [1][2][3]. Recently, long-distance spin transport and electrical control over spin signal and lifetimes have been achieved in graphene at room temperature [1,2,[4][5][6][7][8]. These experiments have shown that the performance of graphene spintronic devices is significantly affected by the quality of the contacts to the graphene channel [1]. Conventional spin transport experiments use top ferromagnetic metal/tunnel barrier contacts of micrometer width on graphene channels for the purpose of spin injection and detection [1][2][3]. However, the fabrication of atomically smooth oxide tunnel barriers is challenging, as they usually contain pinholes, roughness, and defects [9][10][11]. The use of such wide contacts is known to limit the device performance due to inhomogeneous injection and detection of spins, and due to interface-induced spin dephasing under the contacts [9][10][11]. Additionally, at the nanoscale graphene edges can become important and are predicted to host spin-polarized edge states, giving rise to spin filtering [12,13]. Such spintronic properties are also predicted to be tunable by external electric fields [12,13]. If realized, this would add the spin degree of freedom to graphene-based devices and circuits, where spin currents can be generated and injected from zigzag nanoribbons to graphene without the need of ferromagnetic spin injectors.
Here we present the fabrication and characterization of a device that utilizes 1D ferromagnetic edge contacts to graphene/hexagonal boron nitride (h-BN) heterostructures. The non-local spin-valve measurements are found to show spurious magnetoresistance effects due to local Hall effects. Our detailed measurements in different geometries and gate voltages, as well as calculations show that the ferromagnetic contactinduced fringe fields give rise to such a signal.

Results and discussion
A schematic representation and optical microscopy picture of the device are shown in figures 1(a) and (b), respectively. The devices with 1D edge contacts to graphene/h-BN heterostructures were fabricated using the process steps presented in figure 1(c). The shape anisotropy of the ferromagnetic contacts (different widths) was used to achieve different switching fields. The details about the fabrication steps of heterostructure and 1D edge contacts are described in Methods section.
The basic characterization of the fabricated graphene/h-BN heterostructure is depicted in figure 2. The Raman spectrum of the heterostructure is shown in figure 2(a), where the G and 2D peaks related to graphene [14] are visible at positions 1568 cm −1 and 2632 cm −1 respectively, and the h-BN peak [15] is at 1343 cm −1 . The 1D edge contact resistances to graphene were found to be in the range of 4-16 kΩ. The electrical properties of the graphene channel were characterized by gate voltage (V g ) dependence in both the local (figure 2(b)) and nonlocal (R NL = V NL /I bias , figure 2(c)) configurations. The charge neutrality (Dirac) point of graphene was found to be in the range of V D ~ −5 to +26 V for different channels measured in local or nonlocal configurations between different contacts due to variations of doping levels within the graphene sheet. The presence of a measurable nonlocal resistance R NL in the detection circuit indicates a charge current spreading outside the bias current circuit, giving rise to an Ohmic resistance contribution R s ·e −πL/W , where R s is graphene sheet resistance, and L and W are the length and width of the graphene channel, respectively [16].
Next, nonlocal magnetoresistance measurements were carried out in the 1D edge contact devices as shown in figure 3(a). We observed a nonlocal volt age V NL with single switching and hysteresis behavior while sweeping the in-plane magnetic field aligned with the contacts at fixed bias currents ( figure 3(b)). By changing the current direction, a similar switching signal with opposite sign was observed. The complete bias dependence of the signal is presented in figure 3(c), which shows a linear dependence in the measured bias range. The temperature dependence of the magnetoresistance signal V NL was measured at a constant current bias of I = +15 µA (figure 4). We observed a decay of the switching amplitude of V NL (ΔV ↑↓ ) with increasing temperature, which could be measured up to 200 K.  The presence of only one step of magnetoresistance for each sweep direction is not typical for spin signals that arise from spin injection and detection [6,11,[17][18][19][20]. For comparison, a typical spin-valve signal for top ferromagnetic contacts to graphene is shown in supplementary figure S1(b) (stacks.iop.org/ TDM/5/014001/mmedia). In the latter case, at least two steps are usually visible for each sweep direction, when both injector and detector contacts switch their magnetization direction [6,11,[17][18][19][20]. The coercivity values of the ferromagnetic contacts used here are within the typical sweeping field range, which we have also verified with higher field sweep ranges (supplementary figure S2(b)). A single-step switching signal can arise if the graphene edge itself generates a spin current; however, such effects are only expected in graphene nanoribbons and not in the micrometer-scale devices used here [12,13]. The quantum spin Hall effect is also unlikely to be responsible for the observed signal considering the negligible spin-orbit coupling in graphene. Additional confirmation of the absence of spin transport in our 1D contact devices comes from the out-of-plane field sweeps, where no Hanle spin precession signal is observed (supplementary figure S2(c)). The continuous linear change of magnetoresistance here (inset in supplementary figure S2(c)) could be due to Hall effect in the presence of continuous sweeps of external perpendicular magnetic field, in contrast to abrupt changes of stray fields as e.g. in figure 3. Therefore, the observed single-switching magnetoresistance with 1D ferromagnetic contacts could be due to the local Hall effect in the graphene in the presence of stray magnetic fields emanating from the edges of ferromagnetic contacts [10,[21][22][23][24].
In order to further clarify the origin of the magnetoresistance, gate-dependent measurements were carried out, where we tune the concentration and type of the charge carriers in graphene. Figure 5(a) shows the magnetoresistance switching with the application of V g = ±30 V, where a change in the sign of the ΔV ↑↓ is observed due to electron or hole conduction of the graphene channel. The complete gate dependence of ΔV ↑↓ , along with the channel resistance is shown in figure 5(b), revealing a correlation between the sign of ΔV ↑↓ and type of charge carrier in graphene. The absolute value of ΔV ↑↓ is found to have two maxima  near the charge neutrality point in graphene, where the charge density n 2D is minimal. These results support the argument in favor of a local Hall effect-dominated magnetoresistance switching. At the same time, this gives additional evidence to rule out any spin-related nature of the observed switching, since the spin signal should not change sign with a change of charge carrier type [17].
Next, regular Hall measurements were performed (figure 6(a)) on the same device at 75 K to correlate with the magnetoresistance measurements. Figure 6(b) shows a change in the sign of the slope of the Hall response V H in the graphene channel by the application of a gate voltage V g = ±40 V, due to a change in conduction from electrons to holes. The full gatedependent Hall measurements were performed by sweeping V g from electron to hole conduction across the Dirac point in the presence of different perpendicular magnetic fields (figure 6(c)). We observe a change in the amplitude and sign of V H with the applied V g . A clear similarity is also observed between the V g dependence of the magnetoresistance signal ΔV ↑↓ (figure 5(b)) and the regular Hall voltage V H (figure 6(c)). This similarity is justified due to the similar origin of the local Hall magnetoresistance and the regular Hall effects, which are due to stray magnetic fields from the ferromagnetic contacts, or the Lorentz force acting on moving charges in the presence of a perpendicular external magnetic field, respectively.
In order to quantify the impact of the local Hall effect on the magnetoresistance signal, we use the Biot-Savart law to numerically calculate the stray magnetic fields that can arise from the 1D edge contacts. The contact geometry we consider is shown in figure 7(a). To describe the contact we assume a surface current density given by the saturation magnetization of cobalt, M s = 1.42 MA m −1 , which is parallel to the +x (−x) axis on the top (bottom) surfaces of the side view, and parallel to the +z (−z) axis on the left (right) surfaces of the top view. Away from the sloped region, this magnetizes the cobalt along the y-axis. In the sloped region, this surface current also yields an out-of-plane component to the polarization. For metal surfaces parallel to the principal axes, the stray magnetic field can  be calculated analytically [26], while the sloped portion of the contact requires a numerical solution.
In figure 7(b) we show the spatial distribution of the perpendicular component of the stray magnetic field B z within the graphene layer, using the experimental contact geometry (with t h-BN,lower = 0 and the rest of the parameters listed in the caption of figure 7). The black dashed lines indicate the contact metal edges. Here we see that the stray field can be quite strong, reaching more than 600 mT at the graphene edge. In addition, the end of the top portion of the ferromagnetic contact also induces a strong perpendicular stray field on the order of 350 mT. Near the 1D edge contact to graphene, the stray field decays from 600 mT to zero over a distance of 50 nm. This average field of 300 mT corresponds to an average magnetic length of ~50 nm, indicating that the stray field can indeed induce a significant Hall effect at the charge injection/detection point.
In figure 7(c) we plot B z along the middle of the contact, indicated by the white dashed line in figure 7(b). This plot shows that by increasing the thickness of the bottom h-BN layer, the magnitude of the stray field at the graphene edge can be reduced by nearly a factor of three for a typical set of experimental parameters. This reduction occurs because the stray fields due to the top and bottom surfaces of the contact tend to cancel one another as the vertical position of the graphene layer increases. While the perpendicular stray field at the injection edge is still relatively large, ~250 mT, further optimization of the contact geometry can reduce this. For example, tuning the ratio of t h-BN,upper /t h-BN,lower can reduce B z at the injection edge ( figure 7(d)), but this can result in deeper penetration of stray fields into the graphene sheet, and in sign changes of B z at different distances from the graphene edge (figure 7(e)). Generally, with thicker top and bottom h-BN layers, the stray field at the edge can be reduced to below 100 mT (figure 7(f) and supplementary figure S3(a)). A shallower etching angle θ could also reduce the stray fields at the graphene edge (supplementary figure S3(b)), but this is difficult to tune experimentally. Additionally, thicker top h-BN layer can significantly reduce the stray fields at the contact edge corresponding to y = 300 nm (supplementary figure S3(c)). It should also be noted that this calculation overestimates the stray field by a factor M s /M r , where M r is the remanent magnetization of the contact.
A similar magnetoresistance switching effect can also arise due to tunneling anisotropic magnetoresistance (TAMR) [27][28][29][30][31][32]. In magnetic tunnel junctions, the TAMR signal strongly depends on the orientation of the magnetization with respect to the current direction, crystallographic axes, spin-orbit interaction (SOI) and density of states (DOS) anisotropies in the materials. In such devices, TAMR signals were measured at low temperatures and the magnitude is found to enhance in the presence of heavy elements. However, the TAMR signal is generally observed only at very low temperatures below 100 K due to sampling a wider region of the tunneling DOS at higher temperatures. In contrast, our magnetoresistance signal in graphene devices is persistent up to 200 K and vanishes only due to increased noise level of the junctions at higher temperatures. The TAMR effects in literature are also known to have a characteristic bias, temperature and angle dependence, and are strongly dependent on the interface SOI. Further detailed investigations are required to identify the contribution of TAMR in the observed magnetoresistance signal (in addition to the stray Hall effect) in our graphene devices with 1D ferromagnetic contacts.

Conclusions
We have investigated the operation of a graphene device with 1D ferromagnetic edge contacts in a nonlocal configuration. The analysis of the measurements indicates a local Hall effect-dominated phenomenon in the presence of charge current spreading outside the injection circuit, resulting in a single-step magnetoresistance switching. This local Hall magnetoresistance can arise due to the Lorentz force acting on moving charges in the presence of ferromagnetic stray fields from the 1D edge contacts. This origin of the signal is also supported by the absence of a Hanle signal and the change in sign of the magnetoresistance with the change in carrier type from electron to hole in the graphene transport channel. Our calculations reveal that with proper optimization of the contact geometry, the magnitude of the perpendicular stray field can be reduced to below 100 mT at the contact edge. Further investigations are required to identify other contributions to the measured signal such as TAMR in our graphene devices. These findings provide insight into the phenomena that hinder the operation of nonlocal spin valve devices with 1D edge contacts, and they lay the foundation for further developments.

Methods
To fabricate h-BN encapsulated graphene devices, first CVD graphene was wet-transferred onto a Si/SiO 2 substrate (from Graphenea). After cleaning and annealing in Ar/H 2 atmosphere, h-BN flakes were dry-transferred using the scotch-tape method. The graphene was patterned by means of oxygen plasma etching where h-BN flakes were used as a mask. Electron beam lithography and electron beam evaporation of metals were utilized to fabricate the contacts (1 nm TiO 2 /65 nm Co), which formed a 1D connection to encapsulated graphene. Electronic and magnetotransport measurements were performed in a cryostat by means of a Keithley 6221 current source and a Keithley 2182A nanovoltmeter. Gate voltage was applied using a Keithley 2612 source metre.