Using simultaneous magnetic force microscopy and transport measurements of a graphene spin valve, we correlate the non-local spin signal with the magnetization of the device electrodes. The imaged magnetization states corroborate the influence of each electrode within a one-dimensional spin transport model and provide evidence linking domain wall pinning to additional features in the transport signal.

Electrical injection and detection of spin have become common techniques to study spin transport and relaxation in metals and semiconductors1–6 and are attractive for technological implementation of spin-based logic. Graphene has demonstrated great promise as the spin transport channel in such devices.4,7,8 Electrical transport measurements of graphene non-local spin valves have been used to demonstrate spin lifetimes of 500 ps and spin diffusion lengths of 2.4 μm at room temperature9 (up to 2 ns and 12 μm for devices encapsulated in hexagonal boron nitride10). However, theoretical expectations for spin lifetime in graphene are orders of magnitude longer.11,12 Optimization of spintronic devices for application-readiness requires understanding all aspects of device behavior and noise, suggesting the need for measurement tools and techniques beyond electrical transport alone. Considering the great insight into charge transport and electronic behavior gained by combining scanning probe microscopy (SPM) with charge transport measurements,13–18 similar advances should be possible by integration of scanned probes with spin transport experiments. The measurement technique demonstrated in this letter—simultaneous magnetic force microscopy (MFM) and spin transport—provides independent and complementary probes of magnetization, and facilitates a more complete understanding of the interaction between spin and magnetization in nanoscale devices. This approach provides a useful tool for studying complex magnetization configurations and their relationship to transport in other technologically relevant devices.19 Furthermore, extending the combination of transport and SPM to spin transport systems could illuminate spin behavior by enabling spatially resolved studies not possible with transport measurements alone.20,21 While optical techniques provide a route to imaging spin transport in optically active systems,22 the addition of SPM opens the door for optically inaccessible samples such as silicon and graphene.

Here we report direct imaging of electrode magnetization of an operating non-local graphene spin valve. These images allow unambiguous determination of the overall magnetic state of the device and can be correlated with simultaneously acquired non-local magnetoresistance. Using this correlation, we verify that a one-dimensional spin diffusion model,23 combined with electrode magnetization switching, quantitatively captures the primary behavior of the device—high, low, and intermediate resistance states.

Direct correlation of imaging and transport also provides insight into the origins of various observed features in the transport signal beyond simple binary switching of electrode magnetization. We observe domain wall pinning and de-pinning in the ferromagnetic electrodes concurrent with otherwise obscure transitions in the non-local signal. Without imaging, the origins of these occasionally detected transitions4 would remain unclear. Our observation adds to a physical understanding of these sources of noise-like variations and aligns with the previous observation that ferromagnetic contacts can have a strong influence on spins in the transport channel via stray field effects.24 

The graphene spin valves used in this study consist of exfoliated single layer graphene with Co/MgO ferromagnetic contacts (see Ref. 25 for details). Using a combination of two-point and four-point measurements, we find the particular device described in this manuscript exhibits a graphene resistance of 1.46 kΩ and an average contact resistance of 6.13 kΩ.

Spin transport measurements are performed using the traditional four-terminal non-local scheme3,4 (see device schematic in Fig. 1(a)). All measurements are performed at room temperature and in vacuum. We source 10 μA of current at I+ (drain at I) at 11 Hz, and measure the non-local voltage VNL=V+V with a lock-in amplifier. Results are plotted as RNL=VNL/I. An external magnetic field is swept in-plane, parallel to the long axis of the Co electrodes, in order to obtain the non-local magnetoresistance (NLMR) signal (Fig. 1(b)). From a Hanle measurement of this device (with magnetic field oriented perpendicular to both the substrate plane and the injected spin orientation), we obtained a spin lifetime τs = 476 ps and spin diffusion length λs = 3.10 μm by using the measured device resistances and the fit algorithm of Ref. 26. These parameters are used as inputs to the spin diffusion model to produce the spin electrochemical potential (ECP) profiles shown in Figs. 2, 3, and 4 and to fit our NLMR data (Fig. 1, solid lines).

FIG. 1.

(a) Experimental schematic showing non-local circuit configuration of a graphene spin valve and the scanned magnetic probe positioned above the device. The device contains 7 magnetic Co electrodes, 4 of which are used for transport measurements. One possible magnetic configuration is shown with colored arrows representing magnetization direction. (b) Electrically detected non-local spin valve signal. At magnetic field values indicated by numbered triangles, we acquired MFM images. A fit to the data is performed using the MFM-determined magnetic configuration of each state (solid lines).

FIG. 1.

(a) Experimental schematic showing non-local circuit configuration of a graphene spin valve and the scanned magnetic probe positioned above the device. The device contains 7 magnetic Co electrodes, 4 of which are used for transport measurements. One possible magnetic configuration is shown with colored arrows representing magnetization direction. (b) Electrically detected non-local spin valve signal. At magnetic field values indicated by numbered triangles, we acquired MFM images. A fit to the data is performed using the MFM-determined magnetic configuration of each state (solid lines).

Close modal
FIG. 2.

Left: MFM images acquired during magnetic field up sweep. Color scale represents measured cantilever frequency shift, caused by tip-sample forces. Images correlate with the RNL measurements in Fig. 1(b) indicated by black numbered triangles. Right: 1-D spin diffusion calculation, showing spin-resolved electrochemical potential as a function of position in the graphene channel for up (blue) and down (red) spins. Positions of the four circuit electrodes (V,V+,I+, and I) are indicated with dashed lines, and the overall device state (magnetization direction of electrodes) is indicated with up (blue) and down (red) arrows in the inset. The measured voltage difference VNL is proportional to the spin density difference V+V highlighted in these 1-D calculations.

FIG. 2.

Left: MFM images acquired during magnetic field up sweep. Color scale represents measured cantilever frequency shift, caused by tip-sample forces. Images correlate with the RNL measurements in Fig. 1(b) indicated by black numbered triangles. Right: 1-D spin diffusion calculation, showing spin-resolved electrochemical potential as a function of position in the graphene channel for up (blue) and down (red) spins. Positions of the four circuit electrodes (V,V+,I+, and I) are indicated with dashed lines, and the overall device state (magnetization direction of electrodes) is indicated with up (blue) and down (red) arrows in the inset. The measured voltage difference VNL is proportional to the spin density difference V+V highlighted in these 1-D calculations.

Close modal
FIG. 3.

Sequence of MFM images and spin ECP profiles correlated with RNL measurements at specific field values (red numbered triangles in Fig. 1(b)) acquired on magnetic field down sweep. The MFM tip magnetization has reversed prior to imaging state 6, inverting the frequency shift color contrast relative to Fig. 2. The electrode magnetizations are indicated by colored arrows in the spin density profile insets using the convention established in Fig. 2.

FIG. 3.

Sequence of MFM images and spin ECP profiles correlated with RNL measurements at specific field values (red numbered triangles in Fig. 1(b)) acquired on magnetic field down sweep. The MFM tip magnetization has reversed prior to imaging state 6, inverting the frequency shift color contrast relative to Fig. 2. The electrode magnetizations are indicated by colored arrows in the spin density profile insets using the convention established in Fig. 2.

Close modal
FIG. 4.

(a) Down-sweep NLMR signal. Green triangle indicates field value at which the MFM image shown in (b) was acquired. In this MFM image, a pinned domain wall can be clearly seen on the electrode between V and V+. (c) Modeled spin ECP with (VNL,DW, green) and without (VNL, black) the pinned domain wall. The stray field from the domain wall causes local dephasing and a reduction in VNL (green curve). (d) Up-sweep NLMR signal. Green triangle indicates field value at which the MFM image shown in (e) was acquired. In this MFM image, the outermost electrode has reversed (relative to the state 1). Immediately before this electrode's magnetization is reversed, the NLMR signal was seen to increase, and then return to baseline upon reversal.

FIG. 4.

(a) Down-sweep NLMR signal. Green triangle indicates field value at which the MFM image shown in (b) was acquired. In this MFM image, a pinned domain wall can be clearly seen on the electrode between V and V+. (c) Modeled spin ECP with (VNL,DW, green) and without (VNL, black) the pinned domain wall. The stray field from the domain wall causes local dephasing and a reduction in VNL (green curve). (d) Up-sweep NLMR signal. Green triangle indicates field value at which the MFM image shown in (e) was acquired. In this MFM image, the outermost electrode has reversed (relative to the state 1). Immediately before this electrode's magnetization is reversed, the NLMR signal was seen to increase, and then return to baseline upon reversal.

Close modal

We have previously described the custom-built scanning probe microscope which enables simultaneous acquisition of force microscopy and transport measurements.27 This system utilizes frequency-shift detection of an oscillating cantilever to image atomic, electrostatic, and magnetic forces. The MFM signature of a bar magnet (which adequately approximates the Co electrodes) creates opposite frequency shifts at opposing ends of the bar, where the field gradient is strongest. This MFM fingerprint enables clear determination of the magnetization orientation of each electrode.

At specific fields during the NLMR measurement field sweep (indicated by numbered triangles in Fig. 1(b)), we acquired MFM images of the graphene flake and the nearby ends of all seven Co electrodes. Figs. 2 and 3 show the sequence of MFM images, constituting direct observation of magnetization orientation and switching sequence in an operational non-local spin valve. These magnetization states can be explicitly correlated with the measured NLMR.

To begin the measurement sequence, the field is first ramped to −250 mT to initialize all electrodes to the parallel state. The up sweep non-local voltage data are then recorded starting at −50 mT. The first MFM image is acquired on the up sweep at 0 mT applied field, verifying the all-parallel configuration (Fig. 2, subpanel 1). The graphene flake is outlined in dashed lines for reference. It is faintly observable in the force detection images owing to the periodic 11 Hz voltage applied to drive the injection current. This sinusoidal voltage is experienced by the scanned cantilever as a periodic frequency shift due to the capacitive interaction between the grounded, conducting cantilever and the biased sample.28 

On the up sweep, 5 distinct magnetic configurations are observed via MFM, corresponding to four unique RNL levels (the all-parallel configurations in subpanels 1 and 5 exhibit the same resistance). From this sequence of images, and correlation with the measured RNL levels at the corresponding field, it is clear that all four circuit electrodes play critical roles in determining the measured non-local voltage. The injector (I+) and extractor (I) contacts determine the steady-state spin polarization in the channel, while the detector (V+) and reference (V) contacts probe the spin chemical potential set by their respective magnetic orientations.

A similar set of states are observed on the down sweep, as shown in Fig. 3. The electrode magnetizations are again initialized to an all-parallel state with a large +250 mT field. This field reverses the MFM tip magnetization, which causes an inversion of the frequency shift contrast in Fig. 3 relative to Fig. 2. Therefore, the electrode magnetizations in states 5 and 6 are identical.

With knowledge of the total magnetization state determined by the MFM images, we can model the spin ECP in the channel to understand the magnitude and sign of the measured non-local voltage. The spin ECP spatial profile for each of the 9 observed electrode configurations is calculated using a 1-D spin diffusion model.23 These are shown immediately adjacent to the MFM images in Figs. 2 and 3 in order to highlight the correlation between the imaged magnetization orientation and the corresponding spin profile in the graphene. Inputs to the model include the geometry of our device, the switching fields of each electrode as determined by MFM and transport, and the spin diffusion length (determined from a Hanle measurement). The 1-D model does not include contact-induced spin relaxation29 (i.e., once injected, graphene spins cannot escape through the metallic electrodes). We find this to be justifiable given measured device resistances (Rcontact>Rgraphene) and the good agreement between experiment and fit. Using this model, we can replicate the 6 unique voltage levels obtained in the electrically detected non-local spin signal as shown from the fit in Fig. 1 (solid lines). There are no free parameters in the fit aside from a scaling of the overall magnitude, and a DC offset to match the measured data. The relative voltage levels of the various resistance states are fixed by the experimentally determined model inputs discussed above (i.e., the voltage level of each state cannot float freely with respect to the others).

The MFM images also uncover device behavior beyond the model of single domain switching. We can correlate several transitions in the NLMR signal (Figs. 4(a) and 4(d), green triangles) with interesting features in the MFM images obtained at the same field. This direct imaging provides evidence linking magnetization domain wall pinning and de-pinning to the observed changes in the NLMR signal. Uncovering this connection between domain wall motion and spin transport is a key new capability made possible by simultaneous imaging and correlated transport.

Fig. 4(b) shows a domain wall pinned directly above the graphene channel on the Co electrode between V and V+ at a field of −13.6 mT. This domain wall produces a large, but localized, stray field experienced by spins in the underlying graphene. These spins are dephased by the field, causing a change in RNL indicated by the green triangle in Fig. 4(a). Following the methods of Ref. 20, we can model spin transport in the presence of such a localized, inhomogeneous field to obtain the new steady-state spin density. Qualitatively, we reproduce the observed reduction in spin signal. We find that the dephasing has a greater impact on spins below V+ than those under V, thereby reducing the voltage difference VNL (see spin density profiles, Fig. 4(c)). Growth of Co on graphene can result in enhanced grain and domain formation of Co due to graphene's high surface diffusion,25,30,31 and domains could be pinned above the graphene by these grain boundaries.

A similar feature is observed during the up sweep in Fig. 4(d), between 8 and 9.4 mT. At 9.4 mT, we acquired the MFM image in (e) showing that the leftmost electrode (external to the non-local circuit) had fully reversed relative to the 0 mT image of Fig. 2, subpanel 1. Prior to this reversal, a domain wall necessarily propagated across the graphene channel which would have produced a stray field experienced by the graphene spins. Modeling spin diffusion in the presence of this localized field reproduces the qualitative behavior of the electrically detected NLMR: the non-local voltage is enhanced by the dephasing effect of the stray field (since spin accumulation near V is primarily affected). This enhancement persists until the electrode has completely reversed at 9.4 mT and the stray field has been eliminated.

In conclusion, we have utilized simultaneous transport and magnetic force microscopy to correlate the complete magnetization state of a non-local graphene spin valve with its spin signal. We find good agreement between measured voltages and a 1-D spin diffusion model informed by the MFM images. Such direct correlation studies also uncover device behavior that would be impossible to understand with transport measurements alone—namely, spin transport sensitivity to pinning and de-pinning of ferromagnetic domains in the contact electrodes.

These results point to possible studies and applications in which spins interact with a much richer set of states than the two allowed orientations of a monolithically magnetized ferromagnet. Understanding the interactions between spin transport and mobile domain walls or other magnetic textures in nanoscale geometries is an area of intense current research.19 In-operando MFM of functioning devices is well suited for such studies. For example, the fast and sensitive force detection provided by our microscope27 is able to detect domain wall transit below the scanned tip, which could be correlated with high-bandwidth features in transport data. Furthermore, the local spin dephasing due to stray fields from the contact electrodes is identical to the spin imaging mechanism used in scanning spin precession microscopy,20,21 and suggests that spins in graphene should be amenable to such an imaging technique. Overall, this approach shows clear promise for the precise characterization of spintronic device performance, and the development of reliable, application-ready devices.

Funding for this research was primarily provided by the Center for Emergent Materials, an NSF MRSEC at The Ohio State University (Award No. DMR-1420451). V.B. acknowledges the support from the Army Research Office under Grant No. W911NF-12-1-0587.

1.
M.
Johnson
and
R. H.
Silsbee
,
Phys. Rev. Lett.
55
,
1790
(
1985
).
2.
A. T.
Hanbicki
,
B. T.
Jonker
,
G.
Itskos
,
G.
Kioseoglou
, and
A.
Petrou
,
Appl. Phys. Lett.
80
,
1240
(
2002
).
3.
X.
Lou
,
C.
Adelmann
,
S. A.
Crooker
,
E. S.
Garlid
,
J.
Zhang
,
K. S. M.
Reddy
,
S. D.
Flexner
,
C. J.
Palmstrom
, and
P. A.
Crowell
,
Nat. Phys.
3
,
197
(
2007
).
4.
N.
Tombros
,
C.
Jozsa
,
M.
Popinciuc
,
H. T.
Jonkman
, and
B. J.
van Wees
,
Nature
448
,
571
(
2007
).
5.
I.
Appelbaum
,
B.
Huang
, and
D. J.
Monsma
,
Nature
447
,
295
(
2007
).
6.
S. P.
Dash
,
S.
Sharma
,
R. S.
Patel
,
M. P.
de Jong
, and
R.
Jansen
,
Nature
462
,
491
(
2009
).
7.
W.
Han
,
K.
Pi
,
W.
Bao
,
K. M.
McCreary
,
Y.
Li
,
W. H.
Wang
,
C. N.
Lau
, and
R. K.
Kawakami
,
Appl. Phys. Lett.
94
,
222109
(
2009
).
8.
W.
Han
,
K.
Pi
,
K. M.
McCreary
,
Y.
Li
,
J. J. I.
Wong
,
A. G.
Swartz
, and
R. K.
Kawakami
,
Phys. Rev. Lett.
105
,
167202
(
2010
).
9.
W.
Han
and
R. K.
Kawakami
,
Phys. Rev. Lett.
107
,
047207
(
2011
).
10.
M.
Guimares
,
P.
Zomer
,
J.
Ingla-Ayns
,
J.
Brant
,
N.
Tombros
, and
B.
van Wees
,
Phys. Rev. Lett.
113
,
086602
(
2014
).
11.
D.
Huertas-Hernando
,
F.
Guinea
, and
A.
Brataas
,
Eur. Phys. J. Spec. Top.
148
,
177
(
2007
).
12.
C.
Ertler
,
S.
Konschuh
,
M.
Gmitra
, and
J.
Fabian
,
Phys. Rev. B
80
,
041405
(
2009
).
13.
M. A.
Topinka
,
B. J.
LeRoy
,
S. E. J.
Shaw
,
E. J.
Heller
,
R. M.
Westervelt
,
K. D.
Maranowski
, and
A. C.
Gossard
,
Science
289
,
2323
(
2000
).
14.
M. A.
Topinka
,
B. J.
LeRoy
,
R. M.
Westervelt
,
S. E. J.
Shaw
,
R.
Fleischmann
,
E. J.
Heller
,
K. D.
Maranowski
, and
A. C.
Gossard
,
Nature
410
,
183
(
2001
).
15.
J.
Berezovsky
,
M. F.
Borunda
,
E. J.
Heller
, and
R. M.
Westervelt
,
Nanotechnology
21
,
274013
(
2010
).
16.
J.
Martin
,
N.
Akerman
,
G.
Ulbricht
,
T.
Lohmann
,
J. H.
Smet
,
K.
von Klitzing
, and
A.
Yacoby
,
Nat. Phys.
4
,
144
(
2008
).
17.
R.
Jalilian
,
L. A.
Jauregui
,
G.
Lopez
,
J.
Tian
,
C.
Roecker
,
M. M.
Yazdanpanah
,
R. W.
Cohn
,
I.
Jovanovic
, and
Y. P.
Chen
,
Nanotechnology
22
,
295705
(
2011
).
18.
V.
Panchal
,
A.
Lartsev
,
A.
Manzin
,
R.
Yakimova
,
A.
Tzalenchuk
, and
O.
Kazakova
,
Sci. Rep.
4
,
5881
(
2014
).
19.
M.
Hayashi
,
L.
Thomas
,
R.
Moriya
,
C.
Rettner
, and
S. S. P.
Parkin
,
Science
320
,
209
(
2008
).
20.
V. P.
Bhallamudi
,
A. J.
Berger
,
D. E.
Labanowski
,
D.
Stroud
, and
P. C.
Hammel
,
J. Appl. Phys.
111
,
013902
(
2012
).
21.
V. P.
Bhallamudi
,
C. S.
Wolfe
,
V. P.
Amin
,
D. E.
Labanowski
,
A. J.
Berger
,
D.
Stroud
,
J.
Sinova
, and
P. C.
Hammel
,
Phys. Rev. Lett.
111
,
117201
(
2013
).
22.
M.
Furis
,
D. L.
Smith
,
S.
Kos
,
E. S.
Garlid
,
K. S. M.
Reddy
,
C. J.
Palmstrom
,
P. A.
Crowell
, and
S. A.
Crooker
,
New J. Phys.
9
,
347
(
2007
).
23.
M. V.
Costache
,
M.
Zaffalon
, and
B. J.
van Wees
,
Phys. Rev. B
74
,
012412
(
2006
).
24.
S. P.
Dash
,
S.
Sharma
,
J. C.
Le Breton
,
J.
Peiro
,
H.
Jaffrs
,
J.-M.
George
,
A.
Lematre
, and
R.
Jansen
,
Phys. Rev. B
84
,
054410
(
2011
).
25.
W.
Han
,
K. M.
McCreary
,
K.
Pi
,
W. H.
Wang
,
Y.
Li
,
H.
Wen
,
J. R.
Chen
, and
R. K.
Kawakami
,
J. Magn. Magn. Mater.
324
,
369
(
2012
).
26.
E.
Sosenko
,
H.
Wei
, and
V.
Aji
,
Phys. Rev. B
89
,
245436
(
2014
).
27.
A. J.
Berger
,
M. R.
Page
,
J.
Jacob
,
J. R.
Young
,
J.
Lewis
,
L.
Wenzel
,
V. P.
Bhallamudi
,
E.
Johnston-Halperin
,
D. V.
Pelekhov
, and
P. C.
Hammel
,
Rev. Sci. Instrum.
85
,
123702
(
2014
).
28.
29.
T.
Maassen
,
I. J.
Vera-Marun
,
M. H. D.
Guimares
, and
B. J.
van Wees
,
Phys. Rev. B
86
,
235408
(
2012
).
30.
A. J.
Berger
,
W.
Amamou
,
S. P.
White
,
R.
Adur
,
Y.
Pu
,
R. K.
Kawakami
, and
P. C.
Hammel
,
J. Appl. Phys.
115
,
17C510
(
2014
).
31.
Z.
Zhou
,
F.
Gao
, and
D. W.
Goodman
,
Surf. Sci.
604
,
L31
(
2010
).