Introduction

The 3d transition metal ferromagnets possess a large exchange interaction of the order of 1 eV, which has the effect of shifting the energy of their 3d electronic band structure depending on its spin1. This exchange shift leads to the presence of a magnetization, and of a spin polarization of the density of states at the Fermi level, which induces a spin polarization of the conductivity2. These have led to large magnetoresistance effects in spin valve structures3 which exploit the spin polarization of the conductivity through the giant magnetoresistance (GMR), as well as the tunnel magnetoresistance (TMR)4 which relies on the spin polarization of the density of states. GMR and TMR allowed the development of spintronics devices where information is stored in the parallel/antiparallel orientation of the magnetic layers, while these magnetoresistance effects provide effective means to readout the magnetic configuration. Writing the magnetic information however remains challenging, as it relies on energy consuming processes, either using an external magnetic fields, spin transfer or spin orbit torques involving energy dissipative currents, to reverse the magnetization direction5.

The possibility to electrostatically reverse the spin polarization would lift this limitation of spintronic devices6. This electrostatic tunability of the magnetization direction is however not possible in conventional metallic ferromagnets owing to their large carrier density. In contrast, graphene is an air stable metallic material, with a low carrier concentration which allows the electrostatic control of its transport properties7. In addition, the possibility to modify its band structure through proximity effects8, and its spin filtering abilities9, has led graphene to find its way into van der Waals heterojunctions10,11. However, owing to its lack of magnetic properties, graphene is not able to detect or generate a spin current by itself and remains a passive component of spintronic devices.

The magnetic proximity effect (MPE)12,13 has changed this paradigm. It has indeed been predicted14,15 and demonstrated12,16,17,18,19,20,21,22,23 that having graphene in contact with a magnetic material can induce an exchange energy shift in its band structure. It has been experimentally shown24 that this shift leads to the appearance of a large spin polarization in graphene. This type of system could therefore be used in the predicted spin valves devices25,26,27,28.

A spin polarization with a value of 14% has been reported in graphene interfaced with CrSBr29, an air stable and insulating layered type A antiferromagnetic van der Waals material30,31,32. It is expected that a proper electrostatic adjustment of the Fermi energy could lead to an electrostatically sign reversible ± 100% spin polarization in this system owing to the electron/hole symmetry of its band structure. Such an effect has however not been reported so far. Indeed, while large exchange shifts have been reported in proximity-magnetized graphene33, the region close to the charge neutrality Dirac point has not been explored. This is important since, as we will show, an efficient electrostatic control and sign reversal of the spin polarization can only be achieved when the Fermi energy lies in between the exchange shifted Dirac points.

In this paper we provide a direct demonstration of the tunability of the spin polarization of the carrier density, which reaches values of + 69% to − 50%. We use high field magnetotransport measurements to extract the exchange energy shift of the graphene band structure. The presence of the exchange shift modifies the well established energy spectrum of the graphene Landau Levels (LLs)34, resulting in an unusual manifestation of the electronic transport in the Quantum Hall Effect (QHE) regime due to spin polarized counterflowing electron and hole edge channels. These high field quantum transport measurements allow us to obtain the specific values of the exchange shift using a self-consistent model of the spin polarized QHE edge channels, and of the spin polarization at zero magnetic field in proximitized graphene without relying on ferromagnetic contacts. Our results confirm the possibility to electrostatically control the sign and amplitude of the spin polarization in graphene.

The device discussed in this work is shown in Fig. 1(a). It consists of a Hall bar-shaped bilayer graphene exfoliated on SiO2, and on top of which a 100-nm thick CrSBr flake has been deposited [Fig. 1(b)]. The CrSBr consists of antiferromagnetically coupled layers, such that the layer in contact with graphene induces an uniform MPE. Magnetotransport measurement were performed at a temperature of 20 K and under the application of an out of plane magnetic field B. Additional magnetotransport measurements were performed at 180 K, above the Néel temperature of CrSBr. Details on the sample preparation and magnetotransport measurement procedure can be found in the “Method" section and in the Supplementary Information.

Fig. 1: Device structure.
figure 1

a Optical micrograph of the bilayer graphene-CrSBr device. The white outline is the area covered by the exfoliated graphene flake, the red outline is the graphene Hall bar geometry after excess graphene was etched away. The green-yellow bar in the centre is the CrSBr flake, its non-uniform color is due to the irregular height of its surface. The large number of contacts allows to perform magnetotransport measurements on various electrically independent areas of the device to confirm the reproducibility of the magnetic proximity effect. b Cross section of the device presented in (a) showing the doped Si backgate separated from the device by an insulating SiO2 layer.

Results and Discussion

Theory

We performed four probe measurements of the magnetic field and gate dependence of longitudinal [Fig. 2(a)] and transverse (Hall) resistances [Fig. 2(b) and (c)]. The MPE in graphene is expected to lift the spin degeneracy17,18,20,22,29, causing a modification of the conventional QHE measured in pristine graphene35. In this section we provide a theoretical framework which describes the effect of the exchange energy shift on the high field magnetotransport in bilayer graphene.

Fig. 2: Magnetotransport measurements.
figure 2

a Rxx measured between contacts 9 and 6 as a function of the gate voltage. This quantity is a square sheet resistance. b Rxy measured between contacts 10 and 12 (dataset A). c Rxy measured between contacts 8 and 14 (dataset B). Rxx and Rxy have been symmetrized and anti-symmetrized, respectively, with respect to the magnetic field to remove a small cross contamination between Rxx and Rxy due to a possible misalignment of the voltage contacts (see Supplementary Section VII for the raw data).

Figure 3a qualitatively shows the effect of the MPE on the graphene density of states. It induces an exchange energy shift \(\pm \frac{\Delta }{2}\) of the electronic states depending on their spin. The electron and hole carrier densities can be indexed in separate populations according to their spin state. In the following they will be denoted nσ and pσ for electrons and holes with spin σ =  or . The density of these spin-polarized carriers depends on the position of the Fermi energy (EF). When EF lies in between the two exchange shifted Dirac points, the system is in a two carrier regime with n = 0 and p = 0, so that there is a one-to-one correspondence between the electron (hole) states and the spin down (up) states. We call the position for which n = p the charge compensation point (CCP), which corresponds to EF = 0 meV.

Fig. 3: Exchange shifted densities of state, Landau level energies and counterflowing edge channels in proximity-magnetized graphene.
figure 3

a Left: Spin-up polarized holes (red) and spin-down polarized electrons (blue) densities of state of proximity-magnetized bilayer graphene at B = 0 T and Δ = 27 meV. The exchange shifted Dirac points are marked as DP. EF (grey dashed line) is located in the 2 carrier regime, corresponding to Vg − Vg0 = − 1.3 V. Right: LL energies in bilayer graphene at B = 2.2 T with Δ = 27 meV. Blue and red lines corresponds to electron and hole LLs, respectively. Solid (dashed) lines are occupied (unoccupied) energy levels. The arrows indicate the spin state. The numbers indicate the indices of the spin polarized LLs, with the energy of the degenerate n = (0, 1) LLs at the exchange shifted Dirac points and independent on the magnetic field. b LL spatial dispersion and edge state chemical potentials of the sample section located between contacts 3 and 2, and 5 and 6. The Fermi energy is the same as in (a). c Coexisting spin-down electron (blue) and spin-up hole (red) quantum Hall edge channels. Channels originating from the source contact (4) are in solid lines, while those originating from the drain (1) are in dashed lines. The μi with i = 1 to 6 are the contact electrochemical potentials.

The application of a strong magnetic field causes the formation of coexisting spin-polarized electron and hole LLs [Fig. 3(a)]. As depicted in Fig. 3(b), this results in the formation of edge channels with opposite flow direction owing to their different spatial dispersion relations close to the edges. These counterflowing edge channels will be populated up to the electrochemical potentials of the contacts from which they originate, and equilibrate in the voltage contacts36,37. In contrast to the case of the conventional QHE, the electrochemical potential of the side contacts will therefore not be a direct copy of the source or drain contacts. This will result in a non-zero longitudinal resistance and in a non-conventional transverse resistance as observed in Fig. 2.

We made a self-consistent model to describe the electronic transport (Supplementary Section V). The carrier density of the thermally-broadened exchange shifted LLs is first self-consistently evaluated as a function of the gate voltage. We thus obtain the Fermi energy as a function of the gate voltage and magnetic field. It is then possible to evaluate the finite temperature conductances GN(P) = (2e2/h) N(P) of the spin-polarized electron (N) and holes (P) counter flowing edge channels between adjacent contacts. N and P correspond to the (thermal) occupation of the spin-polarized edge channels. We applied a gate voltage ranging from 0 V–7 V, for which the system was found to be in the spin-polarized two carrier regime. Finally, the voltage Vk at the contact numbered k is evaluated in the framework of the Landauer-Buttiker formalism38. The net charge current Ik flowing in or out of the contact k is:

$${I}_{k}=\frac{2{e}^{2}}{h}\left[-(N+P){V}_{k}+N{V}_{k+1}+P{V}_{k-1}\right]$$
(1)

so that the voltage at contacts in which no net charge current flows (Ik = 0) will be a weighted average of its nearest neighbours:

$${V}_{k}=\frac{N{V}_{{{{{{{{\rm{k}}}}}}}}+1}+P{V}_{{{{{{{{\rm{k}}}}}}}}-1}}{N+P}$$
(2)

This model relies on the following assumptions: (1) The counterflowing edge channels with opposite spins do not equilibrate between the contacts. (2) There is no bulk transport pathway which allows current flow other than through the edge channels. (3) The counter-flowing edge channels are fully absorbed and the chemical potentials carried by them are equilibrated at the voltage contacts36,39.

We now apply our model to the example diagram with six contacts shown in Fig. 3(c), in which the LLs and Fermi energies are given in Fig. 3(a), with their spatial dispersion relation in Fig. 3(b) for an illustrative set of parameters. The longitudinal resistance is (at T = 0 K) is:

$${R}_{32}=\frac{{V}_{3}-{V}_{2}}{I}=\frac{h}{2{e}^{2}}\frac{14}{63}$$
(3)

and the transverse resistance:

$${R}_{26}=\frac{{V}_{2}-{V}_{6}}{I}=\frac{h}{2{e}^{2}}\frac{1}{3}$$
(4)

Eq. (3) shows that in contrast to the conventional QHE, the longitudinal resistance is non zero. Eq. (4) predicts the appearance of transverse resistance plateau-like features at values which are not expected for the conventional QHE in non-magnetic graphene40,41. These features originate from the specific ways in which the top and bottom edge channel electrochemical potentials change when moving between contacts from the current source to the drain, and this is directly linked to the presence of counter propagating edge channels. This spatial dependence of the contact electrochemical potential for the device presented in Fig. 1(a) is reported in Supplementary Section V. An important consequence of this model is that for N = P the electrochemical potential difference between the top and bottom contacts vanishes, leading to a zero transverse voltage. Furthermore the number of occupied electron/hole edge channels depends crucially on Δ, which is the relevant parameter we will extract from the model.

Comparison with measurement results

In this section we discuss the results of the magnetotransport measurements, and their interpretation in the framework of the model presented in the Theory section. The bias current flows between contact 11 and 1 [see Fig. 1], while the longitudinal sheet resistance Rxx was measured between contacts 9 and 6 [Fig. 2(a)]. The transverse resistance Rxy was measured between two pairs of contacts: 10 and 12 (dataset A) [Fig. 2b], 8 and 14 (dataset B) [Fig. 2c].

We first focus on the low field results of dataset A. Shifting Vg from 0 – 7 V changes the low field slope of Rxy from negative to positive, with a sign change at 3.1 V, which is defined as Vg0, corresponding to the CCP. For Vg < Vg0, the majority of the carriers are spin-up holes. For Vg > Vg0 the majority of the carriers are spin-down electrons. A similar behavior was observed in dataset B where Vg0 ≈ 3.3 V. The low field value of Rxx is weakly gate voltage dependent. This feature is interpreted in the framework of our model as a consequence of the exchange shift, which creates a gate voltage-independent total carrier density close to the charge compensation point (see below).

In contrast to the usual longitudinal resistance in bilayer graphene subjected to a large magnetic field, Rxx does not go to zero [Fig. 2(a)]. It instead displays a minimum and then a sharp upturn (at 4.5 T for Vg = 0). The minimum progressively becomes less pronounced when increasing Vg and the magnetic field at which it occurs decreases down to 3.3 T at Vg = Vg0, then increases again when Vg > Vg0. At a given gate voltage, this minimum of Rxx corresponds to an extremum of Rxy which displays a similar voltage dependence [Fig. 2b, c]. The gate dependence of the magnetic field at which this feature occurs is reported in the Supplementary Section V. This behavior has not been observed in the high field magnetotransport of pristine graphene35. In the next section, the observed features are interpreted in the framework of the counterpropagating spin up and spin down channels model described in the Theory Section.

The carrier mobility μ was extracted from the low field longitudinal resistance [Fig. 2a] using 1/ρ = (n + p)eμ with ρ the resistivity, −e the charge of the electron. Finally, n = n + n and p = p + p are the total electron and hole densities, the determination of which as a function of the applied gate is described in the Supplementary Section II. The mobility depends weakly on the gate voltage, with values ranging from 1.42 – 1.72 m2 V−1 s−1 (see Supplementary Section III for the gate dependence).

The main feature to be analyzed using our model is the non-monotonic behavior of the transverse resistance and its dependence on the gate voltage. In the following we illustrate the effect of the exchange energy shift on the magnetotransport in two carrier regime ranging from electron (Vg > Vg0) to hole (Vg < Vg0) majority. We then show that the observed features are well reproduced by our model calculations. An analysis similar to the one developed below is given for other sets of data in Supplementary Section VI.

In Fig. 4a, b, we show the longitudinal resistance from Fig. 2a and transverse resistance extracted from the dataset A [Fig. 2b] at Vg = 5 V. The system is in a spin down-electron majority regime (n > p) with a carrier density estimated in the Supplementary Section II. We used the model described in the Theory Section to calculate the transverse resistance, with the energy shift Δ as a fitting parameter. The calculated Δ value which gives a Rxy best reproducing the measured one displayed in Fig. 4b is 27 meV. In Fig. 4c, d, we show the LL positions calculated using this parameter at B = 2.8 and 6.2 T, respectively. We also show the position of the Fermi energy at Vg = 5 V. Fig. 4c shows the LL energies calculated with a magnetic field of 2.8 T. Here the only populated LLs are the spin down n = (0, 1) and n = 2, and the spin up n = (0, 1) (here n refers to the bilayer graphene LL number). When the magnetic field goes above 6.2 T, as shown in Fig. 4d, the spin down n = 2 LL has moved above the Fermi energy and became depopulated. The occupation of spin up (P) and down (N) channels therefore becomes approximately the same, causing Rxy to decrease. Since the measurements have been carried out at 20 K, the spin up n = 2 channel is still partially populated for a magnetic field of 6.5 T, so that a complete cancellation of Rxy is not reached. The longitudinal resistance in Fig. 4(a) could not be reproduced using our model owing to the possible presence of bulk transport paths. We show in the Supplementary Section VI that our model nevertheless qualitatively predicts the upturn of Rxx at 6.2 T.

Fig. 4: Comparison of the model prediction and the transport measurements.
figure 4

a Sheet resistance Rxx and b Rxy from dataset A. The gate voltage is Vg = 5 V. The dashed line corresponds to the Rxx calculated using Δ = 27 meV. c LLs energies calculated using Δ = 27 meV and B = 2.8 T, corresponding to the first vertical dashed line on the left. d same as c with B = 6.2 T, corresponding to the second vertical dashed line on the left. e Rxx and f Rxy from dataset B, Vg = 1.8 V. The dashed line corresponds to the transverse resistance using Δ = 32 meV. g LLs energies calculated using Δ = 32 mev and B = 3 T.

The longitudinal and transverse resistance from dataset B [Fig. 2c] at Vg = 1.8 V, close to the charge compensation point (Vg0 = 2.2 V) in this region of the sample, are shown in Fig. 4e, f, respectively. At this gate voltage there is a spin up-hole majority (p > n). We used the procedure described above and find Δ = 32 meV. Here, although our model qualitatively reproduces the features of Rxy, the calculated amplitude is much larger than the measured one. We believe that this discrepancy originates from our simplifying assumption of edge channel-only electronic transport. The calculated results are therefore normalized for the experimental maximum in order to qualitatively compare the main features of Rxy. Both Rxx and Rxy display an extremum at B ≈ 3 T. The corresponding Fermi and LLs energies are shown in Fig. 4g. At this field, the n = − 1 spin up and n = 2 spin down LLs have moved respectively below/above the Fermi energy, such that only the n = (0, 1) LLs for these spin states remain occupied, and N and P both come close to 1. This same number for the electron and hole edge channels causes the transverse resistance to further decrease when the magnetic field goes above 6.5 T. Another noteworthy feature in Fig. 4(f) is the small shoulder at B = 1.7 T, this is caused by the n = − 2 spin up channel moving below EF at this field value. This feature is also reproduced by our model. When moving the Fermi level from the electron to the hole-majority regime, we observed that the features described in Fig. 4 are anti-symmetric in terms of relative gate voltage Vg − Vg0 for the same magnetic fields (Supplementary Section VI). This indicates that the sign and magnitude of the exchange shift is independent on the Fermi energy in the range of gate voltages used in this study.

We observed a variation of Δ extracted from different areas of the device, ranging from 27 – 32 meV. This indicates that the MPE is not fully homogeneous across our device. These values are however in agreement with the 20 meV extracted previously on a similar system using a different method29. The energy shift corresponds to an effective magnetic field Beff with Δ = gμBBeff22, in which g ≈ 2 is the Landé factor and μB the Bohr magneton. The resulting effective fields vary from 231 to 276 T. During the measurements our applied magnetic field never exceeds 7 T, such that the Zeeman and valley splittings can be neglected42,43.

The exchange energy shift originates from the MPE and is therefore also present at low and zero magnetic field. Figure 5a shows the calculated electron/hole density as a function of the gate voltage. The presence of an exchange shift implies the existence of a two carrier regime for Vg − Vg0 < 4.5 V.

Fig. 5: Extraction of the spin polarization of the carrier density as a function of the gate voltage.
figure 5

a Calculated electron and hole densities. b Calculated (dashed line) and the measured (red dots) Hall coefficient RH from dataset A as a function of the gate voltage, and c calculated (colored line) and measured (black line) spin polarization. A temperature of 20 K and an exchange shift Δ = 27 meV were used for the calculations.

The Hall coefficient RH for two carrier n and p with the same mobilities is expressed as eRH = (n − p)/(n + p)2 (see Supplementary Section II) and corresponds to the slope of the transverse resistance measured at low field. The Hall coefficient is expected to display a smooth crossing and a sign change at the CCP when an exchange energy shift is present in the graphene band structure. It then reaches a maximum and decreases when the system enters the single carrier regimes (n = 0 or p = 0). In Fig. 5b we plot the Hall coefficient calculated using the n and p extracted from the model described in the Supplementary Section II, and using Δ = 27 meV, extracted from the high field magnetotransport measurements. Figure 5b also shows the measured low field Hall coefficient. We find a good agreement between our model and the experimental values close to the charge compensation point, while a deviation occurs away from the CCP.

The exchange shift induces an equilibrium magnetization in graphene, expressed as MGr = gμB(n − p + p − n). In bilayer graphene, MGr = gμB(meff/π2)Δ, which is independent of the Fermi energy since the density of states is approximately energy independent in this system. We however expect this magnetization to be energy dependent in monolayer graphene. This offers the possibility for the electrostatic manipulation of the magnetization. Another relevant quantity is the spin polarization, which is defined as:

$$P=\frac{{n}_{\uparrow }+{p}_{\uparrow }-{n}_{\downarrow }-{p}_{\downarrow }}{{n}_{\uparrow }+{p}_{\uparrow }+{n}_{\downarrow }+{p}_{\downarrow }}(\times 100\%)$$
(5)

In the following we use the extracted Δ = 27 meV to evaluate P.

In the two carrier regime, there is a one to one correspondence between the electron (hole) and spin down (up) densities, such that P = (p − n)/(p + n). Here, Vg < Vg0 corresponds to a hole-spin up majority regime, which translates into a positive spin polarisation [Fig. 5c]. Similarly, Vg > Vg0 corresponds to a negative spin polarisation at the Fermi energy. At T = 0 K, the spin polarization can reach ± 100% at the spin down (up) Dirac point since n = 0 (p = 0). The polarization is zero at the charge compensation point where n = p, i.e. for Vg = Vg0. When moving outside the two carrier regime for Vg − Vg0 > 4.5 V, the spin up and down density of both electrons and holes increase, causing a decrease of the spin polarization at large positive and negative gate voltage [Fig. 5c]. In the two carrier regime, the Hall coefficient writes eRH = − P/(n + p), hence providing a direct measurement of the spin polarization. In Fig. 5c we show the spin polarization evaluated using the measured Hall coefficient and the calculated n and p. We observe a polarization ranging from + 69% at Vg = 0 V to − 50% Vg = 7 V. At zero magnetic field the CrSBr magnetization is transverse to the flake long axis in Fig. 1a, and sets the direction of the graphene spin polarization29.

The conductivity of spin () carriers is expressed as σ() = e(n() + p())μ(). The symmetry of the quantum Hall magnetoresistance when moving from the electron to hole-dominated regime indicates that the exchange shift modifies the electron and hole LL energies in the same manner. We therefore expect the mobility to be independent of the spin of the carrier, such that μ = μ. In this case, the spin polarization of the carrier density and the spin polarization of the conductivity are directly related.

It is now possible to evaluate the GMR of a lateral all-graphene electrostatically-modulated heterostructure. We use a spin polarization of  ± 60%, leading to a gate-tunable GMR of up 50% for systems with lateral lengths comparable to the spin diffusion length. For a system with an energy-dependent density of states, such as monolayer graphene, it is also possible to view the polarization in our system as a spin polarization of the density of states at the Fermi energy. In this case, one may expect a tunneling magnetoresistance (TMR)44 of 112%. Here we propose as a proof of concept the magnetic tunnel junction made from proximity-magnetized graphene shown in Fig. 6, in which the switching between “parallel" and “anti-parallel" states can be achieved by applying different gate voltages to the graphene layers.

Fig. 6: Proposed electrostatically-controlled tunnel junction made of proximity-magnetized graphene.
figure 6

The high resistance ("anti-parallel") case is on the left while the low resistance ("parallel") case is on the right. The thickness of the arrows represent the magnitude of the spin up (red) and down (blue) tunneling current. An electrostatic bias eV is applied between the graphene layers to achieve the high and low resistance states.

We performed low and high field magnetotransport measurement in proximity-magnetized graphene, and extracted an exchange shift of 27 – 32 meV, and showed that it results in a spin polarization of the carrier density. The spin polarization is tunable electrostatically, and reached values between +69 and −50%.

Although our modelling matched closely, both qualitatively and quantitatively, with our magnetotransport measurements, discrepancies between the model and the experiment remain. The most striking one is the deviation of the measured Hall coefficient from its predicted values when moving away from the charge neutrality point. Our model assumes that the only effect of the MPE on the band structure of graphene is to cause a constant exchange energy shift, which does not depend on the carrier type or on the electronic state energy. This assumption might be too restrictive to fully account for the details of our data in view of theoretical works predicting a carrier dependent exchange shift, as well as a possible gap opening at the graphene Dirac point45. Furthermore, an energy dependence of the exchange shift can be expected in this type of systems46,47, and the MPE strongly depends on the stacking between graphene and CrSBr as shown in Supplementary Section II.D. Further refinement of our model is possible to test these theoretical predictions.

Overall, the proximity-magnetized graphene system appears in this work as promising for the development of spintronic devices based on the electrostatic tunability of magnetic properties. It indeed displays possibilities for the realization of gate-controllable spin valves and spin filters, and will be of great significance in the exploitation of spin-charge interconversion phenomena in van der Waals heterostructures8,48,49 owing to the optimized spin injection efficiency offered by the combination of a large electrostatically-tunable spin polarization with the long spin relaxation length of graphene29. We also predict that proximitized graphene hosts an equilibrium magnetization, which could be modified upon the application of an a.c. or d.c. gate voltage. We show that our results can be used to electrostatically generate the GMR and TMR effects.

Methods

The device presented in Fig. 1 was fabricated dry-transferring a CrSBr flake (thickness ~ 100 nm) on top of a bilayer graphene flake which was exfoliated on a Si/SiO2(285 nm) substrate. Ti/Au contacts were deposited on the graphene in a Hall bar pattern, and the excess graphene outside of the Hall bar was etched away. We used conventional lock-in measurement techniques, with a bias current of 10 – 100 nA and frequency 13.7 Hz. A gate voltage Vg ranging from 0 to 7 V was applied between the p++ doped Si substrate and graphene during the measurements.