Abstract
The magnetic proximity effect can induce a spin dependent exchange shift in the band structure of graphene. This produces a magnetization and a spin polarization of the electron/hole carriers in this material, paving the way for its use as an active component in spintronics devices. The electrostatic control of this spin polarization in graphene has however never been demonstrated so far. We show that interfacing graphene with the van der Waals antiferromagnet CrSBr results in an unconventional manifestation of the quantum Hall effect, which can be attributed to the presence of counterflowing spin-polarized edge channels originating from the spin-dependent exchange shift in graphene. We extract an exchange shift ranging from 27 – 32 meV, and show that it also produces an electrostatically tunable spin polarization of the electron/hole carriers in graphene ranging from − 50% to + 69% in the absence of a magnetic field. This proof of principle provides a starting point for the use of graphene as an electrostatically tunable source of spin current and could allow this system to generate a large magnetoresistance in gate tunable spin valve devices.
Similar content being viewed by others
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.
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.
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.
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:
so that the voltage at contacts in which no net charge current flows (Ik = 0) will be a weighted average of its nearest neighbours:
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:
and the transverse resistance:
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.
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.
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:
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.
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.
Data availability
All data supporting this study are available at https://doi.org/10.5281/zenodo.10654565.
References
Himpsel, F. J., Ortega, J. E., Mankey, G. J., & Willis, R. F. Magnetic nanostructures. Adv. Phys. 47, 511–597 (1998)
Žutić, I., Fabian, J. & Sarma, S. D. Spintronics: Fundamentals and applications. Rev. Mod. Phys. 76, 323 (2004).
Diény, B. Giant magnetoresistance in spin-valve multilayers. J. Magn. Magn. Mater. 136, 335–359 (1994).
Butler, W. H. Tunneling magnetoresistance from a symmetry filtering effect. Sci. Technol. Adv. Mater. 9, 014106 (2008).
Dieny, B. et al. Opportunities and challenges for spintronics in the microelectronics industry. Nat. Electron. 3, 446–459 (2020).
Manipatruni, S., Nikonov, D. E. & Young, I. A. Beyond CMOS computing with spin and polarization. Nat. Phys. 14, 338–343 (2018).
Novoselov, K. S. et al. Electric field in atomically thin carbon films. Science 306, 666–669 (2004).
Sierra, J. F., Fabian, J., Kawakami, R. K., Roche, S. & Valenzuela, S. O. Van der Waals heterostructures for spintronics and opto-spintronics. Nat. Nanotechnol. 16, 856–868 (2021).
Piquemal-Banci, M. et al. Spin filtering by proximity effects at hybridized interfaces in spin-valves with 2d graphene barriers. Nat. Commun. 11, 5670 (2020).
Geim, A. K. & Grigorieva, I. V. Van der Waals heterostructures. Nature 499, 419–425 (2013).
Yang, H. et al. Two-dimensional materials prospects for non-volatile spintronic memories. Nature 606, 663–673 (2022).
Bora, M. & Deb, P. Magnetic proximity effect in two-dimensional Van der Waals heterostructure. J. Phys. Mater. 4, 034014 (2021).
Huang, B. et al. Emergent phenomena and proximity effects in two-dimensional magnets and heterostructures. Nat. Mater. 19, 1276–1289 (2020).
Yang, H.-X. et al. Proximity effects induced in graphene by magnetic insulators: first-principles calculations on spin filtering and exchange-splitting gaps. Phys. Rev. Lett. 110, 046603 (2013).
Lee, J. & Fabian, J. Magnetotransport signatures of the proximity exchange and spin-orbit couplings in graphene. Phys. Rev. B 94, 195401 (2016).
Leutenantsmeyer, J. C., Kaverzin, A. A., Wojtaszek, M. & Van Wees, B. J. Proximity induced room temperature ferromagnetism in graphene probed with spin currents. 2D Mater. 4, 014001 (2016).
Wei, P. et al. Strong interfacial exchange field in the graphene/EuS heterostructure. Nat. Mater. 15, 711–716 (2016).
Wu, Y. F. et al. Magnetic proximity effect in graphene coupled to a BiFeO3 nanoplate. Phys. Rev. B 95, 1–7 (2017).
Asshoff, P. et al. Magnetoresistance of vertical Co-graphene-NiFe junctions controlled by charge transfer and proximity-induced spin splitting in graphene. 2D Mater. 4, 031004 (2017).
Wu, Y. et al. Large exchange splitting in monolayer graphene magnetized by an antiferromagnet. Nat. Electron. 3, 604–611 (2020).
Tang, C., Zhang, Z., Lai, S., Tan, Q. & Gao, W.-B. Magnetic proximity effect in graphene/CrBr3 Van der Waals heterostructures. Adv. Mater. 32, 1908498 (2020).
Wu, Y. et al. Magnetic exchange field modulation of quantum Hall ferromagnetism in 2D Van der Waals CrCl3/graphene heterostructures. ACS Appl. Mater. Interfaces 13, 10656–10663 (2021).
Kaverzin, A. A., Ghiasi, T. S., Dismukes, A. H., Roy, X. & Van Wees, B. J. Spin injection by spin–charge coupling in proximity induced magnetic graphene. 2D Mater. 9, 045003 (2022).
Zatko, V. et al. Artificial graphene spin polarized electrode for magnetic tunnel junctions. Nano Lett. 23, 34–41 (2023).
Song, Y. & Dai, G. Spin filter and spin valve in ferromagnetic graphene. Appl. Phys. Lett. 106, 223104 (2015).
Song, Y., Liu, Y., Feng, X., Yan, F. & Zhang, W. Spin-selectable, region-tunable negative differential resistance in graphene double ferromagnetic barriers. Phys. Chem. Chem. Phys. 20, 1560–1567 (2018).
Solis, D., Hallal, A., Waintal, X. & Chshiev, M. Proximity magnetoresistance in graphene induced by magnetic insulators. Phys. Rev. B 100, 104402 (2019).
Zhai, X. & Blanter, Y. M. Proximity-induced diversified magnetic states and electrically controllable spin polarization in bilayer graphene: towards layered spintronics. Phys. Rev. B 106, 075425 (2022).
Ghiasi, T. S. et al. Electrical and thermal generation of spin currents by magnetic bilayer graphene. Nat. Nanotechnol. 16, 788–794 (2021).
Lee, K. et al. Magnetic order and symmetry in the 2D semiconductor CrSBr. Nano Lett. 21, 3511–3517 (2021).
Telford, E. J. et al. Layered antiferromagnetism induces large negative magnetoresistance in the van der waals semiconductor CrSBr. Adv. Mater. 32, 2003240 (2020).
Telford, E. J. et al. Coupling between magnetic order and charge transport in a two-dimensional magnetic semiconductor. Nat. Mater. 21, 754–760 (2022).
Hu, J. et al. Tunable spin-polarized states in graphene on a ferrimagnetic oxide insulator. Adv. Mater. 36, 2305763 (2023).
Novoselov, K. S. et al. Unconventional quantum Hall effect and Berry’s phase of 2π in bilayer graphene. Nat. Phys. 2, 177–180 (2006).
Jiang, Z., Zhang, Y., Tan, Y.-W., Stormer, H. & Kim, P. Quantum Hall effect in graphene. Solid State Commun. 143, 14–19 (2007).
Roth, A. et al. Nonlocal transport in the quantum spin Hall state. Science 325, 294–297 (2009).
Akiho, T., Irie, H., Onomitsu, K. & Muraki, K. Counterflowing edge current and its equilibration in quantum Hall devices with sharp edge potential: roles of incompressible strips and contact configuration. Phys. Rev. B 99, 121303 (2019).
Büttiker, M. Absence of backscattering in the quantum Hall effect in multiprobe conductors. Phys. Rev. B 38, 9375–9389 (1988).
Couëdo, F., Irie, H., Suzuki, K., Onomitsu, K. & Muraki, K. Single-edge transport in an InAs/GaSb quantum spin Hall insulator. Phys. Rev. B 94, 1–5 (2016).
Barlas, Y., Yang, K. & MacDonald, A. Quantum Hall effects in graphene-based two-dimensional electron systems. Nanotechnology 23, 052001 (2012).
Cooper, D. R. et al. Experimental review of graphene. Int. Sch. Res. Not. 2012, 1–56 (2012).
Young, A. F. et al. Spin and valley quantum Hall ferromagnetism in graphene. Nat. Phys. 8, 550–556 (2012).
Young, A. F. et al. Tunable symmetry breaking and helical edge transport in a graphene quantum spin Hall state. Nature 505, 528–532 (2014).
Julliere, M. Tunneling between ferromagnetic films. Phys. Lett. A 54, 225–226 (1975).
Hallal, A., Ibrahim, F., Yang, H., Roche, S. & Chshiev, M. Tailoring magnetic insulator proximity effects in graphene: first-principles calculations. 2D Mater. 4, 025074 (2017).
Zollner, K., Gmitra, M., Frank, T. & Fabian, J. Theory of proximity-induced exchange coupling in graphene on hBN/(Co, Ni). Phys. Rev. B 94, 155441 (2016).
Cardoso, C., Costa, A.T., MacDonald, A. H. & Fernández-Rossier J. Strong magnetic proximity effect in Van der Waals heterostructures driven by direct hybridization. Phys. Rev. B 108, 184423 (2023).
Galceran, R. et al. Control of spin–charge conversion in Van der Waals heterostructures. APL Mater. 9, 100901 (2021).
Safeer, C. et al. Room-temperature spin Hall effect in graphene/MoS2 Van der Waals heterostructures. Nano Lett. 19, 1074–1082 (2019).
Acknowledgements
We acknowledge the financial support of the Zernike Institute for Adv. Mater. and the European Union’s Horizon 2020 research and innovation program under Grant Agreement No. 785219 and No. 881603 (Graphene Flagship Core 2 and Core 3). This project is also financed by the NWO Spinoza prize awarded to B.J. van Wees. by the NWO and has received funding from the European Research Council (ERC) under the European Union’s 2DMAGSPIN (Grant agreement No. 101053054). Synthesis of the crystals was supported by the National Science Foundation (NSF) through the Columbia University Materials Research Science and Engineering Center (MRSEC) on Precision-Assembled Quantum Materials DMR-2011738. Electrical transport measurements at Columbia were supported as part of Programmable Quantum Materials, an Energy Frontier Research Center funded by the US Department of Energy, Office of Science, Basic Energy Sciences, under award DE-SC0019443.
Author information
Authors and Affiliations
Contributions
B.J.v.W. and B.Y. conceived the experiment. B.Y and B.B. fabricated the device, performed the measurements and developed the model. B.J.v.W., B.Y., B.B. and M.C.C. analyzed the data and discussed the experimental results. B.J.v.W., B.Y. and M.C.C. wrote the manuscript. D.G.C., E.J.T. and X.R. performed the growth of the CrSBr crystals. All authors discussed the results and commented on the manuscript. F.I. and M.C. performed the DFT calculations.
Corresponding authors
Ethics declarations
Competing interests
The authors declare no competing interests.
Peer review
Peer review information
Nature Communications thanks the anonymous reviewer(s) for their contribution to the peer review of this work. A peer review file is available.
Additional information
Publisher’s note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Supplementary information
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.
About this article
Cite this article
Yang, B., Bhujel, B., Chica, D.G. et al. Electrostatically controlled spin polarization in Graphene-CrSBr magnetic proximity heterostructures. Nat Commun 15, 4459 (2024). https://doi.org/10.1038/s41467-024-48809-w
Received:
Accepted:
Published:
DOI: https://doi.org/10.1038/s41467-024-48809-w
This article is cited by
-
A seamless graphene spin valve based on proximity to van der Waals magnet Cr2Ge2Te6
Nature Electronics (2024)