Abstract
Van der Waals interactions with transition metal dichalcogenides were shown to induce strong spin-orbit coupling (SOC) in graphene, offering great promises to combine large experimental flexibility of graphene with unique tuning capabilities of the SOC. Here, we probe SOC-driven band splitting and electron dynamics in graphene on WSe2 by measuring ballistic transverse magnetic focusing. We found a clear splitting in the first focusing peak whose evolution in charge density and magnetic field is well reproduced by calculations using the SOC strength of ~ 13 meV, and no splitting in the second peak that indicates stronger Rashba SOC. Possible suppression of electron-electron scatterings was found in temperature dependence measurement. Further, we found that Shubnikov-de Haas oscillations exhibit a weaker band splitting, suggesting that it probes different electron dynamics, calling for a new theory. Our study demonstrates an interesting possibility to exploit ballistic electron motion pronounced in graphene for emerging spin-orbitronics.
Similar content being viewed by others
Introduction
The interfacial interactions with semiconducting transition metal dichalcogenides (TMDCs) have shown to be highly efficient in inducing strong spin-orbit coupling (SOC) in graphene1,2,3. For monolayer graphene, it was theoretically predicted to have two distinctive terms, one that couples out-of-plane spin and valley degrees of freedom (referred to as a spin-valley Zeeman term \({\tau }_{z}{s}_{z}\)) and another that couples in-plane spin and sublattice degrees of freedom, similar to the Rashba term \(({\tau }_{z}{\sigma }_{x}{s}_{y}-{\sigma }_{y}{s}_{x})\), as follows:
where \({H}_{0}\) is the graphene’s Dirac Hamiltonian, \({{{{{\boldsymbol{\sigma }}}}}}{{{{{\boldsymbol{=}}}}}}({\sigma }_{{x}},{\sigma }_{{y}},{\sigma }_{z})\) is a Pauli matrix vector that acts on the sublattice degree of freedom in graphene, \({{{{{\boldsymbol{s}}}}}}=({s}_{x},{s}_{y},{s}_{z})\) is a Pauli matrix vector that acts on spin, and \({\tau }_{z}=\pm 1\) identifies two different valleys in graphene (\(\Delta \approx 0\) due to a large lattice mismatch between the graphene and TMDCs)2,4,5,6,7. Both SOC terms induce band splitting in graphene, and their strengths (\(\lambda\) and \({\lambda }_{R}\)) can be further tuned by an electric field perpendicular to the layers8,9,10,11,12,13, twisting14,15, or by pressure16.
Combined with the high electron mobility and large experimental flexibility of graphene, such a strong interface-induced SOC makes the graphene on TMDCs ideal for ballistic spin-orbitronics where ballistic electron motion can be used to control or detect electron spin through the SOC17,18,19,20,21,22. It is particularly interesting as graphene has shown pronounced ballistic transport effects with large tunability, such as transverse magnetic focusing (TMF)23,24,25, Veselago lensing26,27, Fabry-Pérot interference28,29, and ballistic snake states30,31 among many others. They have also shown unique features originating from the relativistic nature of Dirac electrons23,24,25,26,27,28,29,30,31. However, a vast number of previous studies on graphene-TMDC heterostructures have focused on detecting spin relaxation due to electron scattering rather than the ballistic motion1,2,3,7,16,32,33,34,35,36,37. Moreover, only a few studies have found direct evidence for the SOC-induced band splitting by measuring beatings in Shubnikov-de Hass (SdH) oscillations (for both mono- and bilayer graphene)3,10,38 or tracing changes in quantum capacitance (for bilayer graphene only)11. Not only to understand the effect of the SOC on the electronic properties of the system, such as the band topology39,40 but also to exploit its full potential on ballistic spintronics17,18,19,20,21,22, it is therefore essential to demonstrate the ballistic transport in graphene on TMDCs while simultaneously probing their band structures and electron dynamics.
To fill this missing link, we employ TMF technique in monolayer graphene on WSe2 as it can not only probe the SOC-induced band splitting but also investigate electron dynamics simultaneously (see Fig. 1a–d). TMF occurs when ballistic carriers injected from a narrow aperture (“injector”) at the edge of the sample are subject to a small perpendicular magnetic field (\({{{{{\bf{B}}}}}}=B{{{{{\boldsymbol{z}}}}}}\))41,42,43. Owing to the Lorentz force, the carriers follow skipping cyclotron orbits and focus on another narrow aperture (“collector”) at a distance (\(L\)) that equals an integer multiple of \(2{r}_{c}\) with a cyclotron radius \({r}_{c}=\hslash {k}_{F}/{eB}\), where \(\hslash\) is the reduced Planck constant, \({k}_{F}\) is the Fermi momentum, and \(e\) is the elementary charge. Upon sweeping magnetic fields, the collector voltage will exhibit a set of resonance peaks at certain \(B\)-values determined by \({k}_{F}\),
where \(j\) is an integer and \(\pm\) represents electron and hole for the configuration shown in Fig. 1b. This enables the detection of Fermi surface configurations41,42,43. In the systems with multiple bands23, for instance, there will be multiple sets of resonance peaks at different \(B\)-values from which one can deduce their band structures. Moreover, TMF can also be used to study or control electron dynamics as charge carriers follow skipping cyclotron orbits during the process. In 2D electron gas (2DEG) systems with SOC, the TMF was indeed used to probe spin-orbit split bands and extract SOC strength by studying the separation of the peaks17,19,21, deduce spin polarization by comparing their heights17, or focus spin-polarized current by controlling ballistic electron motion17,19,20.
Here, we demonstrate that all these studies are possible in graphene-TMDC heterostructures. We first show a clear splitting of the first focusing peak in graphene on WSe2, whose evolution in charge density and magnetic field matches well with the theoretical calculations shown in Fig. 1c, d using the SOC strength of ~13 meV (see Supplementary Notes 1,2 for simulation details44,45) and no splitting in the second peak, which indicates a stronger Rashba SOC in the system. From the temperature dependence measurement, we also find a possible suppression of electron-electron scatterings that may originate from the dielectric scattering of WSe2 layers and/or the induced SOC. Further, we show that Shubnikov-de Haas oscillations exhibit a weaker band splitting, suggesting that it probes different electron dynamics. Interestingly, a similar behavior was found in studies on 2DEG with SOC19,21, indicating that this is universal and a new theory is needed to explain the phenomenon. Our study, therefore, places graphene-TMDC heterostructures as an interesting new material platform to explore the effect of SOC on ballistic electron transport in low-dimensional systems.
Results
Sample characterization
To study the ballistic TMF effect, we used a dry pick-up and transfer technique to assemble a stack of hexagonal boron-nitride (hBN), monolayer graphene, and multilayer WSe2 such that the graphene is protected from the harsh chemical environments in the following nanofabrication process46. Standard electron beam lithography and lift-off were carried out to make Hall bar devices on 285-nm-thick SiO2 substrate with doped silicon underneath used as a gate to control charge density \(n\) (see Fig. 1a, b and Methods for fabrication details). Electron transport through the fabricated devices was measured in a 1.5 K variable temperature insert with a 14-T superconducting magnet using a standard low-frequency AC lock-in technique (see Methods). Classical and quantum Hall effect measurements were used to estimate gate capacitance and confirm that the graphene flakes in our samples are monolayers (see Supplementary Fig. 1). Figure 1e shows the carrier density dependence of the four-probe resistances of two different samples 1 and 2 measured at 1.5 K, exhibiting high quality with carrier mobility of 200,000 ~ 400,000 cm2V-1s-1. Especially on the hole side, the symmetry broken quantum-Hall states were observed at magnetic fields as low as ~3 T (marked by black down-triangles in Fig. 1f). This indicates higher hole mobility in our sample than on the electron side. It is also consistent with the observation of the larger negative non-local Hall resistance \({R}_{{ae},{bf}}\) on the hole side (the inset of Fig. 1e) originating from the ballistic electron motion47 (\({R}_{\alpha \beta,\gamma \delta }\equiv {V}_{\gamma \delta }/{I}_{\alpha \beta }\) which refers to the resistance measured by sending a current from contact \(\alpha\) to \(\beta\) and measuring the voltage between contacts \(\gamma\) and \(\delta\)). Having such a high mobility—equivalently, a long mean free path—is important to resolve the small splitting of the focusing peak expected theoretically (Fig. 1d).
Transverse magnetic focusing spectra
The TMF signal (\({R}_{{nl}}={V}_{{cj}}/{I}_{{ai}}\)) is measured in a non-local configuration upon varying \(n\) and \(B\), as depicted in Fig. 1b. Figure 2a, b show the resulting maps of \({R}_{{nl}}\left(n,B\right)\) from sample 1 and 2, respectively. Both exhibit similar TMF spectra and their evolutions in \(n\) and \(B\). Overall, the positions of the \(j\)-th TMF peak in \(B\) follow Eq. (2) with \({k}_{F}=\sqrt{\pi \left|n\right|}\), as expected for monolayer graphene. However, on the hole side (where we found the higher sample quality), we can identify the splitting of the first focusing peak that evolves continuously in \(n\) and \(B\) and no splitting in the second (here, we only focus on the peaks that appear in all density range). Figure 2c, d further magnify the features by plotting 1D cuts \({R}_{{nl}}\left(B\right)\) of the map at different\(\,n\), which qualitatively matches the simulation result shown in Fig. 1d (see Supplementary Fig. 2a for the 1D cuts of the map on the electron side for comparison). Moreover, different from the TMF measurements on pristine monolayer graphene23, we found a large TMF signal constantly exceeding 100 \(\Omega\) around zero density (nearly two orders of magnitude larger than the values at finite densities; see the dark red bands near zero density in the color maps shown in Fig. 2a, b). All features found in the experiment (Fig. 2) match well with our expectations for the graphene with SOC and provide valuable insights about microscopic electron processes in the system, as discussed below.
Analysis of the first focusing peak
First, the first focusing peak splits due to the SOC-induced multiple bands \({S}_{+}\) and \({S}_{-}\) in our sample (Fig. 1c). Note that such a prolonged splitting in both \(n\) and \(B\) has not been observed in other graphene systems without SOC23,24,25. Moreover, we were able to fit the positions of the first focusing peaks with calculations using the SOC strength of \({\lambda }_{{SOC}}\equiv \sqrt{{\lambda }^{2}+{{\lambda }_{R}}^{2}}=13.9\) meV and \(12.1\) meV for samples 1 and 2 respectively as marked by black dotted lines in Fig. 2a, b. Figure 3a and Supplementary Fig. 2b further emphasize the accuracy of the fitting by plotting the average of the normalized difference between the data and the calculation \(\langle {\delta B}^{2}\rangle\) as a function of \(\lambda\) and \({\lambda }_{R}\) in a color scale for devices 1 and 2, respectively (a darker color indicates the smaller \(\langle {\delta B}^{2}\rangle\) so the better fitting; see the caption for more details). We note that the fitting works for any values of \(\lambda\) and \({\lambda }_{R}\) as long as they satisfy \({\lambda }_{{SOC}}=13.0\pm 4.7\) meV for both samples, indicating that the \(\lambda\) and \({\lambda }_{R}\) have a similar effect on the splitting of the first TMF peak—equivalently, the band splitting—in the density range we studied. This is expected as in the density range explored, the Fermi energy is larger than both \(\lambda\) and \({\lambda }_{R}\), leading to an identical splitting in momentum, \(\Delta k=2\sqrt{{\lambda }^{2}+{\lambda }_{R}^{2}}/\hslash {v}_{F}\) with Fermi velocity \({v}_{F}\approx {10}^{6}\,{{{{{{\rm{m}}}}}}\cdot {{{{{\rm{s}}}}}}}^{-1}\) for both \(\lambda\) and \({\lambda }_{R}\). This is consistent with the previous SdH oscillations measurements that probe electron bands3,38.
In addition, we found that the amplitude of the first split peak closer to zero \(B\) is always lower than that of the second (see, e.g., Fig. 2c, d). Interestingly, in 2DEG with Rashba SOC, such uneven heights of the split peaks have been used as a signature of the spin polarization in the system, even though more rigorous analysis and more controlled experiments are needed to identify the origin of the spin polarization17,48,49,50,51. For instance, it was suggested that adiabatic transition between the quantized sub-bands formed at the injector with width \(w\) could polarize electron spin48,51. Slightly modifying the condition derived for 2DEG51, we get \(\frac{\Delta k}{2{k}_{F}} > \frac{3}{16}{\left(\frac{{\lambda }_{F}}{w}\right)}^{2}\to {\lambda }_{R} > \frac{3\pi }{8}\frac{\hslash {v}_{F}{\lambda }_{F}}{{w}^{2}} \, \approx \, 0.25 \sim 0.35\) meV for our sample with Fermi wavelength \({\lambda }_{F} \, \approx \, 30 \sim 40\) nm, and \(w \, \approx \, 300\) nm. The SOC strengths estimated in our study are well in the range (Fig. 3a), but due to the absence of the quantum point contacts in our device, we cannot confirm the formation of the sub-bands at the injector. On the other hand, it was also shown that the relative heights of the split peaks could vary with the distance between the injector and detector due to the difference in scattering lengths for carriers at different spin-orbit-split bands50 and that the exact momentum dependence of the SOC should be considered to understand the spin-polarization49. Interestingly, in graphene on TMDCs, the presence of both spin-valley Zeeman and Rashba SOC terms was predicted to induce a characteristic spin winding of the spin-orbit-coupled bands that leads to a current-induced spin polarization52, which may result in the uneven heights of the peaks. For better understanding, we will need to conduct more sophisticated experiments, such as applying in-plane magnetic fields to control Zeeman energy17,49,53, using samples with various distances between the injector and detector50, or using ferromagnetic contacts for spin-sensitive detection54. Nevertheless, we note that all these previous studies have used the uneven heights of the split peaks as a signature of the spin polarization in the system17,48,49,50,51. Thus, our study shows a possibility of using TMF to detect spin polarization of the ballistic carriers in graphene-TMDC heterostructures.
Analysis of the second focusing peak
The second key finding of this study is the absence of the splitting in the second focusing peak (Fig. 2), which provides more information about the nature of the SOC in the system. It first can be interpreted as the scattering of charge carriers between the spin-orbit-coupled bands \({S}_{\pm }\) at the sample edge, which leads to a single peak as illustrated in the bottom inset of Fig. 2a. To confirm this origin, we have further calculated electron trajectories for the second peak with or without the inter-band transition in the sample in Supplementary Fig. 3. As expected, without the inter-band mixing, the second focusing peak also exhibits splitting. This confirms that the absence of the splitting in the second peak originates from the scattering between the bands \({S}_{\pm }\) at the sample edge.
Interestingly, from the behavior of the second peak, we can learn more about the relative strength of the spin-valley Zeeman and Rashba SOC terms, \(\lambda\) and \({\lambda }_{R}\), because the inter-band scattering at the edge depends sensitively on the spin textures of the split bands \({S}_{\pm }\). For more accessible discussion, let us first consider a single valley only and discuss the effect of the intervalley scattering later. As depicted in Fig. 3b, c, when only spin-valley Zeeman term exists (in other words, when \({\theta }_{{SOC}}=0\); see Fig. 3a), spins in \({S}_{+}\) (\({S}_{-}\)) band are aligned up (down) in z-direction. Thus, when backscattered at the edge, the electron at state A in the band \({S}_{+}\) will jump to state B in the same band unless there are enough magnetic impurities to flip the spin, which is unlikely in high-quality graphene samples like ours. This would lead to the splitting of the second focusing peak as illustrated in Fig. 3b. On the other hand, when the Rashba term dominates (i.e., when \({\theta }_{{SOC}}=\pi /2\)), \({S}_{\pm }\) bands have an opposite spin winding such that the electron at the state A in the band \({S}_{+}\) will jump to the state C in the opposite band, leading to the merging of the peaks as depicted in Fig. 3d, e. Therefore, we can estimate that in our system, the Rashba term dominates. Similarly, studies on 2DEG systems with Rashba SOC have indeed shown no splitting in the second focusing peak19,21,55.
To further confirm our analysis, we have simulated TMF spectra for different values of \({\theta }_{{SOC}}\) in Fig. 3f. As shown in the figure, the splitting in the second peak disappears rapidly as \({\theta }_{{SOC}}\) increases from zero and becomes nearly invisible as \({\theta }_{{SOC}}\gtrsim \pi /4\), consistent with our analysis above. It is, however, worth mentioning that in the simulation, we used an ideal edge, so we may have underestimated the intervalley scattering probabilities that occur in the real sample edge with atomic defects56. Although this does not influence the Rashba-dominating case as the spin winding direction remains the same for the \({S}_{\pm }\) bands in different valleys (see Fig. 3e), it can affect the result when the spin-valley Zeeman term dominates because the spin orientation for each \({S}_{\pm }\) band becomes opposite in different valleys (see Fig. 3c). Thus, the intervalley scattering can lead to the scattering from state A at K valley to state C’ at K’ valley, i.e. the backscattering between the \({S}_{\pm }\) bands at the edge, when the spin-valley Zeeman term dominates, suppressing the splitting in the second peak. Although more experimental and theoretical works are required for a complete understanding of this feature, we can roughly assume that our sample has disordered edges with atomic scale defects with resonance energy near the charge neutrality57, leading to the intervalley scattering rate close to or less than that of the intravalley scattering in the density range explored. In this case, we would still expect to see the splitting when \({\theta }_{{SOC}}=0\). Thus, we believe that the absence of the second peak splitting in Fig. 2 indicates the stronger Rashba SOC in our system.
Large non-local resistance and temperature dependence
We can also explain the observed large \({R}_{{nl}}\) near zero density (Fig. 2a, b) as the presence of the spin Hall effect (SHE)1 in the system. From the very weak temperature dependence of the conductance minimum at zero density (Fig. 4a), we first confirm that this is not from the gap opening at charge neutrality, which may have given a large TMF signal as found in the gapped trilayer graphene23. In contrast, we found a large non-local Hall signal \({R}_{{nl}}^{H}={R}_{{ae},{cg}}\) near the charge neutrality that exceeds the ohmic contribution by about 230 \(\Omega\) (Fig. 4b; here we used Hall probes that are further apart from those used in the inset of Fig. 1e to reduce the influence from the ballistic negative resistance). This is consistent with the SHE found in a similar system1 previously, except that the signal appears near zero density below \(2.0 \sim 3.0\times {10}^{11}\) cm−2 in our device. We attribute this to the crossover from the diffusive regime, where the SHE occurs, to the ballistic regime, where the TMF effect appears. In fact, we found that the density \(2.0 \sim 3.0\times {10}^{11}\) cm-2 coincides well with the value above which the ballistic negative non-local resistance (the inset of Fig. 1e) and TMF signals appear (Fig. 2 and Supplementary Fig. 4) within the resolution of our experiment. More in-depth studies on the crossover between the diffusive spin-Hall and ballistic TMF effects or their coexistence may lead to a better understanding of the charge transport in spin-orbit-coupled systems, and our study shows that it is possible in high-quality graphene-TMDC heterostructures.
Additionally, the observation of both diffusive SHE at low density and ballistic TMF peak splitting at higher density indicates that the SOC in our sample is induced by proximity with TMDC2,4,5,6,7 not by defects1 as the defects would have strongly suppressed ballistic transport in the sample. It confirms that in graphene-TMDC heterostructures—thanks to the atomically sharp interface—the atomic potentials generated by the TMDC can influence graphene band strongly to create an effective Hamiltonian with distinctive SOC terms shown in Eq. (1)2,4,5,6,7. This is similar to how the atomic potentials of boron and nitrogen atoms in hBN creates the moiré minibands in graphene-hBN moiré structure58,59,60,61. We note that such a proximity effect does not require charge carriers in graphene to fill the energy bands in TMDC (or hBN) and thus it can occur even when there are no defect sites in TMDC that can sink charge carriers from graphene and suppress the ballistic transport. Our observation also aligns well with other studies on similar graphene-TMDC heterostructures2,3,7,16,32,33,34,35,36,37,38,62.
We now examine the temperature dependence of the TMF spectra in Fig. 4c, d to study the electron dynamics in the system. Upon increasing temperature, we found that the amplitude of the TMF spectra decreases (Fig. 4c). This suggests enhanced electron scattering at high temperatures. To identify the main scattering mechanism in our system, we extracted the total area below the first focusing peaks \({A}_{1}\) at varying temperatures from 1.5 K to 300 K and calculated the relative scattering length from \({L}_{S}/{L}_{0}=\left(\right.{{{{{\mathrm{ln}}}}}[({A}_{1}(1.5K)/{A}_{1}(T)])}^{-1}\), which is proportional to the effective scattering time (\({L}_{0}\) is the length of the semi-circular electron trajectory corresponding to the first focusing peak)24,25. Figure 4d shows the result exhibiting a clear \({T}^{-1.8}\) dependence on both electron and hole side for different charge densities. This is between electron-phonon scattering (\({T}^{-1}\)) and electron-electron (e-e) scattering (\({T}^{-2}\)), indicating that although the e-e scattering is dominant, it is also slightly suppressed in our sample. In comparison, similar TMF studies on graphene-hBN heterostructures24,25 have shown \({T}^{-2}\) dependence in a wide range of temperatures and charge density. Thus, the \({T}^{-1.8}\) dependence found in our sample should be from the WSe2 itself and/or the induced SOC. Although the exact origin is unclear, we note that the WSe2 has a dielectric constant (\({\varepsilon }_{0}\approx 7.9\)) about twice larger than hBN used in the previous studies (\({\varepsilon }_{0}\approx 3.8\))63 that may induce more screening. In addition, recent studies64,65 have shown that the SOC can affect electronic interaction phenomena such as superconductivity in twisted or Bernal bilayer graphene. This suggests a possibility to measure ballistic transport effects in graphene-TMDC heterostructures to study the effect of SOC on e-e or electron-hole interaction phenomena, such as viscous charge transport66,67,68, and electron-hole collisions69, or superconductivity64,65.
Comparison with Shubnikov–de Haas (SdH) oscillations
To further elucidate the band splitting in our system, we have measured SdH oscillations at a higher magnetic field range. The results are summarized in Fig. 5. In all the density ranges including the electron side, we found beatings in the oscillations originating from the spin-orbit-coupled split bands (Fig. 5b). For quantitative analysis, we performed fast Fourier transforms (FFT) and extracted the frequency \(f\) at which the spectra exhibit a peak (Fig. 5a) which is directly connected to the area of the corresponding Fermi surface by \(f={nh}/2e\) and \(n={{k}_{F}}^{2}/2\pi\) assuming broken spin degeneracy due to SOC. It can therefore be used to estimate the SOC strengths independently. Figure 5c shows the result (we have selected the peaks that evolve continuously in density only), which offers a good agreement with the calculation using the SOC strengths along the circle \(\sqrt{{\lambda }^{2}+{{\lambda }_{R}}^{2}}=3.4\pm 0.7\) meV (Fig. 5d). Interestingly, the fitting is slightly better (i.e., \(\langle {\delta B}^{2}\rangle\) is smaller) near \(\lambda=0\) which indicates the larger Rashba SOC in the system, consistent with our estimation from the analysis on the second focusing peak (Fig. 3b–f). The absolute value is, however, about 4 times smaller than those extracted from the TMF data (Fig. 3a and Supplementary Fig. 2b), indicating that the TMF and SdH oscillations probe different electron dynamics.
To compare the results more directly, we have calculated \(\Delta {B}_{{SdH}}=\left(2\hslash /\left|e\right|L\right)\left|{k}_{F1}-{k}_{F2}\right|\), the expected size of the splitting of the first TMF peak in \(B\) using the \({k}_{F}\)-values extracted from the SdH oscillations (Fig. 5c), and plotted the splitting measured in TMF, \({\varDelta B}_{{TMF}}\) (extracted from Fig. 2), together in Fig. 5e after normalizing the values with the averaged peak positions of the two sub-peaks. As shown in the figure, \(\Delta {B}_{{TMF}}/{B}_{{TMF}}\) remains larger than \(\Delta {B}_{{SdH}}/{B}_{{SdH}}\) in all density range accessed in the experiment. Interestingly, we found that similar behavior was observed in 2DEG systems with SOC19,21, where it was suggested21 that there might be a SOC term, such as a linear-k term, that does not affect the total area of the Fermi surface by shifting a circle in one momentum direction. Since the SdH oscillation requires carriers to complete a full cyclotron motion while in TMF, they only make a half turn, it may be possible that one finds larger splitting in TMF than in SdH, as seen experimentally. However, it is also possible that the SdH oscillations require relatively large magnetic fields to form Landau levels which can induce non-negligible Zeeman energy and may affect the spin-valley Zeeman and Rashba terms in Eq. (1) differently (see Supplementary Note 3 for more discussions). Although more studies are required to understand this discrepancy, our measurement shows that the behavior occurs not only in semiconductor heterostructures-based 2DEG systems19,21 but also in graphene when SOC exists. Therefore, there is likely a fundamental origin behind this phenomenon.
Comparison with previous studies
In Fig. 5d, we also include all SOC strengths extracted from the previous measurements on monolayer graphene-TMDC heterostructures for comparison. Overall, the relaxation time analysis from weak anti-localization or spin-Hall effect measurements1,32,33,37 shows a considerable sample-to-sample variation (Fig. 5d). It can be from the fact that these measurements rely on the model to connect the spin relaxation process in the system with the SOC strength, which is sensitive to the sample-specific electron scattering process7,37. On the other hand, TMF and SdH oscillations directly probe the size of the Fermi surface, which can be compared with theoretically calculated band structures without considering details of the scattering processes. The recent study on SdH oscillation in monolayer graphene-WSe2 heterostructures38 has indeed shown a SOC strength \({\lambda }_{{SOC}}=2.51\) meV close to ours (a dotted circle in Fig. 5d). Interestingly, the study on Landau level splitting62, which is closely related to the SdH oscillations, also showed a similar SOC strength (a square in Fig. 5d). Moreover, in our TMF study, we found similar SOC strengths in two different samples (Fig. 3a and Supplementary Fig. 2b). This further elaborates the benefits of carrying out the (ballistic) transport spectroscopy on understanding electronic properties of the system with SOC.
Discussion
In summary, we have successfully demonstrated the ballistic electron motion in graphene on WSe2 by measuring TMF signals at different charge densities, magnetic fields, and temperatures. From the density and magnetic field dependence of the first focusing peaks (Fig. 2), we confirmed that there exist two split bands in the system as expected theoretically2,4,5,6,7 and estimated the SOC strength of \({\lambda }_{{SOC}}=13.0\pm 4.7\) meV (Fig. 3a and Supplementary Fig. 2b). More interestingly, by analysing the behavior of the second focusing peak that shows no splitting and by carrying out quantum transport simulations, we were able to learn that the Rashba SOC is likely dominant in our system (Fig. 3b-f). Both the presence of the band splitting and a stronger Rashba SOC are well reproduced in SdH oscillations measurements (Fig. 5) even though they showed a smaller SOC strength of \({\lambda }_{{SOC}}=3.4\pm 0.7\) meV. A similar discrepancy was found in other 2DEG with Rashba SOC19,21, indicating that there is a fundamental reason behind it. This calls for a new theory.
In addition to providing spectroscopic evidence of the spin-orbit-coupled bands, our work demonstrates that graphene on TMDCs can support ballistic transport that can be used not only to gain more insights into the microscopic electron process in the system but also to exploit various ballistic transport effects that are pronounced in graphene. It is particularly interesting as, in contrast to the existing studies on graphene spintronics7,22,70,71, TMF separates spin-up and spin-down carriers in real space. This enables the detection and measurement of both spins independently, instead of only the majority one injected from magnetic contacts. This, therefore, offers an alternative venue for graphene spintronic applications7,22,70,71. Similar strategies can also be used to study other 2D materials or heterostructures with strong SOC, such as bilayer graphene-TMDC heterostructures10,11,12,13, black phosphorus72, and more, which will offer new understandings about the effect of SOC in these material systems.
Methods
Sample fabrication
The WSe2, hBN, and graphene flakes were exfoliated from corresponding crystals onto silicon wafers and examined under an optical microscope. The flakes with suitable thicknesses and surfaces were selected and assembled onto highly doped silicon substrates with 285-nm-thick oxide, following the standard dry pick-up and transfer technique46. The WSe2 flakes used in this study are around 20 nm and 40 nm for samples 1 and 2, respectively. After the assembly, the stacks were annealed at 250 °C for 2 h in a tube furnace in Ar/H2 forming gas. 1D electrical contacts were fabricated on the annealed sample through a standard electron-beam lithography and reactive-ion etching (CF4/O2 mixture gas with flow rates of 5/25 sccm, RF power: 60 W), followed by electron beam evaporation of 5 nm Cr and 50 nm Au films. The devices were finally shaped into Hall bars by another electron-beam lithography and reactive-ion etching process.
Electrical measurement
Devices were measured in a 1.5 K cryogen-free variable temperature insert (VTI) with a superconducting magnet. The electrical signals were measured by applying a small low-frequency (17.777 Hz) AC current of 0.1–1 μA between the source and drain terminals and measuring the voltage drop between another two probes using a lock-in amplifier (Stanford Research SR830). The low-noise filters and amplifiers were used to detect small TMF signals. The back gate was controlled by Keithley 2400 source-meter.
Data availability
The data used in this study are freely available in the figshare database at https://doi.org/10.6084/m9.figshare.22644469.
Code availability
Codes to analyze the data and perform numerical calculations are available upon reasonable request.
References
Avsar, A. et al. Spin-orbit proximity effect in graphene. Nat. Commun. 5, 4875 (2014).
Wang, Z. et al. Strong interface-induced spin–orbit interaction in graphene on WS2. Nat. Commun. 6, 8339 (2015).
Wang, Z. et al. Origin and magnitude of ‘designer’ spin-orbit interaction in graphene on semiconducting transition metal dichalcogenides. Phys. Rev. X 6, 041020 (2016).
Gmitra, M. & Fabian, J. Graphene on transition-metal dichalcogenides: a platform for proximity spin-orbit physics and optospintronics. Phys. Rev. B 92, 155403 (2015).
Zubair, M., Vasilopoulos, P. & Tahir, M. Influence of interface induced valley-Zeeman and spin-orbit couplings on transport in heterostructures of graphene on WSe2. Phys. Rev. B 101, 165436 (2020).
Gmitra, M., Kochan, D., Högl, P. & Fabian, J. Trivial and inverted Dirac bands and the emergence of quantum spin Hall states in graphene on transition-metal dichalcogenides. Phys. Rev. B 93, 155104 (2016).
Garcia, J. H., Vila, M., Cummings, A. W. & Roche, S. Spin transport in graphene/transition metal dichalcogenide heterostructures. Chem. Soc. Rev. 47, 3359 (2018).
Khoo, J. Y., Morpurgo, A. F. & Levitov, L. On-demand spin–orbit interaction from which-layer tunability in bilayer graphene. Nano Lett. 17, 7003–7008 (2017).
Gmitra, M. & Fabian, J. Proximity effects in bilayer graphene on monolayer WSe2: Field-effect spin valley locking, spin-orbit valve, and spin transistor. Phys. Rev. Lett. 119, 146401 (2017).
Afzal, A. M. et al. Gate modulation of the spin-orbit interaction in bilayer graphene encapsulated by WS2 films. Sci. Rep. 8, 3412 (2018).
Island, J. O. et al. Spin–orbit-driven band inversion in bilayer graphene by the van der Waals proximity effect. Nature 571, 85–89 (2019).
Wang, D. et al. Quantum Hall effect measurement of spin–orbit coupling strengths in ultraclean bilayer graphene/WSe2 heterostructures. Nano Lett. 19, 7028–7034 (2019).
Ingla-Aynés, J., Herling, F., Fabian, J., Hueso, L. E. & Casanova, F. Electrical control of valley-Zeeman spin-orbit-coupling–induced spin precession at room temperature. Phys. Rev. Lett. 127, 047202 (2021).
David, A., Rakyta, P., Kormányos, A. & Burkard, G. Induced spin-orbit coupling in twisted graphene–transition metal dichalcogenide heterobilayers: Twistronics meets spintronics. Phys. Rev. B 100, 085412 (2019).
Veneri, A., Perkins, D. T. S., Péterfalvi, C. G. & Ferreira, A. Twist angle controlled collinear Edelstein effect in van der Waals heterostructures. Phys. Rev. B 106, L081406 (2022).
Fülöp, B. et al. Boosting proximity spin–orbit coupling in graphene/WSe2 heterostructures via hydrostatic pressure. npj 2D Mater. Appl. 5, 82 (2021).
Rokhinson, L. P., Larkina, V., Lyanda-Geller, Y. B., Pfeiffer, L. N. & West, K. W. Spin separation in cyclotron motion. Phys. Rev. Lett. 93, 146601 (2004).
Frolov, S. M. et al. Ballistic spin resonance. Nature 458, 868–871 (2009).
Lo, S.-T. et al. Controlled spatial separation of spins and coherent dynamics in spin-orbit-coupled nanostructures. Nat. Commun. 8, 15997 (2017).
Kohda, M., Okayasu, T. & Nitta, J. Spin-momentum locked spin manipulation in a two-dimensional Rashba system. Sci. Rep. 9, 1909 (2019).
Rendell, M. J. et al. Gate voltage dependent Rashba spin splitting in hole transverse magnetic focusing. Phys. Rev. B 105, 245305 (2022).
Avsar, A. et al. Colloquium: Spintronics in graphene and other two-dimensional materials. Rev. Mod. Phys. 92, 021003 (2020).
Taychatanapat, T., Watanabe, K., Taniguchi, T. & Jarillo-Herrero, P. Electrically tunable transverse magnetic focusing in graphene. Nat. Phys. 9, 225–229 (2013).
Lee, M. et al. Ballistic miniband conduction in a graphene superlattice. Science 353, 1526–1529 (2016).
Berdyugin, A. I. et al. Minibands in twisted bilayer graphene probed by magnetic focusing. Sci. Adv. 6, eaay7838 (2020).
Lee, G.-H., Park, G.-H. & Lee, H.-J. Observation of negative refraction of Dirac fermions in graphene. Nat. Phys. 11, 925 (2015).
Chen, S. et al. Electron optics with p–n junctions in ballistic graphene. Science 353, 1522–1525 (2016).
Miao, F. et al. Phase-coherent transport in graphene quantum billiards. Science 317, 1530–1533 (2007).
Young, A. F. & Kim, P. Quantum interference and Klein tunnelling in graphene heterojunctions. Nat. Phys. 5, 222–226 (2009).
Taychatanapat, T. et al. Conductance oscillations induced by ballistic snake states in a graphene heterojunction. Nat. Commun. 6, 6093 (2015).
Rickhaus, P. et al. Snake trajectories in ultraclean graphene p–n junctions. Nat. Commun. 6, 6470 (2015).
Völkl, T. et al. Magnetotransport in heterostructures of transition metal dichalcogenides and graphene. Phys. Rev. B 96, 125405 (2017).
Yang, B. et al. Strong electron-hole symmetric Rashba spin-orbit coupling in graphene/monolayer transition metal dichalcogenide heterostructures. Phys. Rev. B 96, 041409 (2017).
Ghiasi, T. S., Ingla-Aynés, J., Kaverzin, A. A. & van Wees, B. J. Large proximity-induced spin lifetime anisotropy in transition-metal dichalcogenide/graphene heterostructures. Nano Lett. 17, 7528–7532 (2017).
Benítez, L. A. et al. Strongly anisotropic spin relaxation in graphene–transition metal dichalcogenide heterostructures at room temperature. Nat. Phys. 14, 303–308 (2018).
Wakamura, T. et al. Strong anisotropic spin-orbit interaction induced in graphene by monolayer WS2. Phys. Rev. Lett. 120, 106802 (2018).
Zihlmann, S. et al. Large spin relaxation anisotropy and valley-Zeeman spin-orbit coupling in WSe2/graphene/h-BN heterostructures. Phys. Rev. B 97, 075434 (2018).
Tiwari, P. et al. Experimental observation of spin−split energy dispersion in high-mobility single-layer graphene/WSe2 heterostructures. npj 2D Mater. Appl. 6, 68 (2022).
Kane, C. L. & Mele, E. J. Z2 Topological order and the quantum spin Hall effect. Phys. Rev. Lett. 95, 146802 (2005).
Kane, C. L. & Mele, E. J. Quantum spin Hall effect in graphene. Phys. Rev. Lett. 95, 226801 (2005).
Tsoi, V. S. Focusing of electrons in a metal by a transverse magnetic field. JETP Lett. 19, 70–71 (1974).
van Houten, H. et al. Coherent electron focusing with quantum point contacts in a two-dimensional electron gas. Phys. Rev. B 39, 8556–8575 (1989).
Tsoi, V. S., Bass, J. & Wyder, P. Studying conduction-electron/interface interactions using transverse electron focusing. Rev. Mod. Phys. 71, 1641–1693 (1999).
Groth, C. W., Wimmer, M., Akhmerov, A. R. & Waintal, X. Kwant: a software package for quantum transport. New. J. Phys. 16, 063065 (2014).
Liu, M.-H. et al. Scalable tight-binding model for graphene. Phys. Rev. Lett. 114, 036601 (2015).
Wang, L. et al. One-dimensional electrical contact to a two-dimensional material. Science 342, 614 (2013).
Mayorov, A. S. et al. Micrometer-scale ballistic transport in encapsulated graphene at room temperature. Nano. Lett. 11, 2396–2399 (2011).
Reynoso, A., Usaj, G. & Balseiro, C. A. Detection of spin polarized currents in quantum point contacts via transverse electron focusing. Phys. Rev. B 75, 085321 (2007).
Chesi, S. et al. Anomalous spin-resolved point-contact transmission of holes due to cubic Rashba spin-orbit coupling. Phys. Rev. Lett. 106, 236601 (2011).
Rendell, M. J. et al. Spin polarization and spin-dependent scattering of holes observed in transverse magnetic focusing. Phys. Rev. B 107, 045304 (2023).
Eto, M., Hayashi, T. & Kurotani, Y. Spin polarization at semiconductor point contacts in absence of magnetic field. J. Phys. Soc. Jpn. 74, 1934–1937 (2005).
Offidani, M., Milletarì, M., Raimondi, R. & Ferreira, A. Optimal charge-to-spin conversion in graphene on transition-metal dichalcogenides. Phys. Rev. Lett. 119, 196801 (2017).
Bladwell, S. & Sushkov, O. P. Magnetic focusing of electrons and holes in the presence of spin-orbit interactions. Phys. Rev. B 92, 235416 (2015).
Raes, B. et al. Determination of the spin-lifetime anisotropy in graphene using oblique spin precession. Nat. Commun. 7, 11444 (2016).
Usaj, G. & Balseiro, C. A. Transverse electron focusing in systems with spin-orbit coupling. Phys. Rev. B 70, 041301 (2004).
Morpurgo, A. F. & Guinea, F. Intervalley scattering, long-range disorder, and effective time-reversal symmetry breaking in graphene. Phys. Rev. Lett. 97, 196804 (2006).
Halbertal, D. et al. Imaging resonant dissipation from individual atomic defects in graphene. Science 358, 1303–1306 (2017).
Yankowitz, M. et al. Emergence of superlattice Dirac points in graphene on hexagonal boron nitride. Nat. Phys. 8, 382–386 (2012).
Ponomarenko, L. A. et al. Cloning of Dirac fermions in graphene superlattices. Nature 497, 594–597 (2013).
Dean, C. R. et al. Hofstadter’s butterfly and the fractal quantum Hall effect in moiré superlattices. Nature 497, 598–602 (2013).
Hunt, B. et al. Massive Dirac fermions and Hofstadter butterfly in a van der Waals heterostructure. Science 340, 1427–1430 (2013).
Wang, D. et al. Spin-orbit coupling and interactions in quantum Hall states of graphene/WSe2 heterobilayers. Phys. Rev. B 104, L201301 (2021).
Laturia, A., Van de Put, M. L. & Vandenberghe, W. G. Dielectric properties of hexagonal boron nitride and transition metal dichalcogenides: from monolayer to bulk. npj 2D Mater. Appl. 2, 6 (2018).
Arora, H. S. et al. Superconductivity in metallic twisted bilayer graphene stabilized by WSe2. Nature 583, 379–384 (2020).
Zhang, Y. et al. Enhanced superconductivity in spin–orbit proximitized bilayer graphene. Nature 613, 268–273 (2023).
Kim, M. et al. Control of electron-electron interaction in graphene by proximity screening. Nat. Commun. 11, 2339 (2020).
Bandurin, D. A. et al. Negative local resistance caused by viscous electron backflow in graphene. Science 351, 1055–1058 (2016).
Crossno, J. et al. Observation of the Dirac fluid and the breakdown of the Wiedemann-Franz law in graphene. Science 351, 1058–1061 (2016).
Nam, Y., Ki, D.-K., Soler-Delgado, D. & Morpurgo, A. F. Electron–hole collision limited transport in charge-neutral bilayer graphene. Nat. Phys. 13, 1207–1214 (2017).
Han, W., Kawakami, R. K., Gmitra, M. & Fabian, J. Graphene spintronics. Nat. Nano. 9, 794–807 (2014).
Roche, S. et al. Graphene spintronics: the European flagship perspective. 2D Mater. 2, 030202 (2015).
Avsar, A. et al. Gate-tunable black phosphorus spin valve with nanosecond spin lifetimes. Nat. Phys. 13, 888–893 (2017).
Acknowledgements
The work is financially supported by the National Key R&D Program of China (2020YFA0309600) and by the University Grants Committee/Research Grant Council of Hong Kong SAR under schemes of Area of Excellence (AoE/P-701/20), ECS (27300819), and GRF (17300020, 17300521, 17309722). K.W. and T.T. acknowledge support from the JSPS KAKENHI (Grant Numbers 21H05233 and 23H02052) and World Premier International Research Center Initiative (WPI), MEXT, Japan. N.W. acknowledges support from William Mong Institute of Nano Science and Technology. W.-H.K. and M.-H.L. gratefully acknowledge National Science and Technology Council of Taiwan (grant numbers: MOST 109-2112-M-006-020-MY3 and NSTC 112-2112-M-006-019-MY3) for financial support and National Center for High-performance Computing (NCHC) for providing computational and storage resources.
Author information
Authors and Affiliations
Contributions
D.-K.K. conceived and supervised the project. M.-H.L. supervised the theoretical part carried out by W.-H.K. Q.R. fabricated the samples and performed the measurements, assisted by H.X. Some of the data was collected in NW’s group with help from Z.Y. and X.F. Q.R. and W.-H.K. analyzed the data, and D.-K.K. and M.-H.L. interpreted them with input from all authors. T.T. and K.W. synthesized the hBN crystals. Q.R., W.-H.K., M.-H.L., and D.-K.K. wrote the paper with input from all authors. All authors discussed the results.
Corresponding authors
Ethics declarations
Competing interests
The authors declare they have no competing interests.
Peer review
Peer review information
Nature Communications thanks Jonghwa Eom, Shun-Tsung Lo and the other, 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 license, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons license, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons license 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 license, visit http://creativecommons.org/licenses/by/4.0/.
About this article
Cite this article
Rao, Q., Kang, WH., Xue, H. et al. Ballistic transport spectroscopy of spin-orbit-coupled bands in monolayer graphene on WSe2. Nat Commun 14, 6124 (2023). https://doi.org/10.1038/s41467-023-41826-1
Received:
Accepted:
Published:
DOI: https://doi.org/10.1038/s41467-023-41826-1
This article is cited by
-
Designer spin-orbit superlattices: symmetry-protected Dirac cones and spin Berry curvature in two-dimensional van der Waals metamaterials
Communications Physics (2024)
-
Probing miniband structure and Hofstadter butterfly in gated graphene superlattices via magnetotransport
npj 2D Materials and Applications (2023)