Abstract
Snake states are trajectories of charge carriers curving back and forth along an interface. There are two types of snake states, formed by either inverting the magnetic field direction or the charge carrier type at an interface. The former has been demonstrated in GaAsâAlGaAs heterostructures, whereas the latter has become conceivable only with the advance of ballistic graphene where a gap-less pân interface governed by Klein tunnelling can be formed. Such snake states were hidden in previous experiments due to limited sample quality. Here we report on magneto-conductance oscillations due to snake states in a ballistic suspended graphene pân junction, which occur already at a very small magnetic field of 20âmT. The visibility of 30% is enabled by Klein collimation. Our finding is firmly supported by quantum transport simulations. We demonstrate the high tunability of the device and operate it in different magnetic field regimes.
Similar content being viewed by others
Introduction
A magnetic field fundamentally modifies the transport properties of an electronic conductor by acting on its charge carriers via the Lorenz force. The most prominent magnetotransport effect is the quantum Hall effect in a two-dimensional electron gas. A strong perpendicular magnetic field forces the charge carriers into one-dimensional conduction channels in which they flow along the edges of a sample. At moderate magnetic fields, however, electron trajectories can be understood in a quasiclassical picture where the Lorenz force bends charge carriers into cyclotron orbits. In the bulk this leads to localization, whereas at the boundary charge carriers can propagate via so-called skipping orbits. Magnetic focusing experiments1,2 represent a direct proof of the skipping orbit picture. In such experiments, an increase of conductance is observed if the distance between two contacts is an integer multiple of the diameter of a cyclotron orbit. One condition for the observation of such trajectories is ballistic transport over the relevant device dimensions. This has limited the observation of skipping orbits to the cleanest available semiconductor samples.
Since 2004, graphene as a new two-dimensional conductor has moved into the focus of condensed-matter research and its behaviour in magnetic field has been intensely investigated. The quality of graphene devices has improved over the recent years and ballistic transport over distances of several microns have been demonstrated recently3,4. Since graphene is a gap-less semiconductor, it offers the possibility of creating internal interfaces with opposite charge carrier polarity. These so-called pân interfaces5,6,7 are formed by local electrostatic gating.
If electrons that propagate via skipping orbits encounter such a pân interface, they will turn into snake states. These states consist of alternating half circles with opposite chirality and they transport current along the interface. Similar snake states have first been realized in GaAs/AlGaAs two-dimensional electron gases by defining regions of alternating magnetic field direction8. These states share the condition of commensurability (similar to the above-described magnetic focusing experiments) with pân snake states but they do not propagate along a single and tunable interface. Snake states in graphene pân junctions were claimed to have been observed in disordered substrate-supported samples9 but the experiment lacked direct evidence for snaking trajectories.
In this article, we investigate ballistic transport across a graphene pân junction in different magnetic field regimes and identify magneto-conductance oscillations as a direct signature of snake states. These findings are supported by detailed tight-binding simulations that allow us to visualize the alternating cyclotron orbits of the snake states.
Results
Evolution of electron trajectories in graphene pân junction
In Fig. 1aâd, we illustrate schematically how trajectories under increasing perpendicular magnetic field evolve in such a device. Figure 1a describes the low-field situation where transport is still dominated by FabryâPérot oscillations with slightly bent trajectories10. As the field is increased, resonant scar states (Fig. 1b) may occur, as observed in semiconductor quantum dots11. At higher fields (Fig. 1c), snake states at the pân interface govern the electronic properties. Finally (Fig. 1d) the system enters the quantum Hall regime12,13 where transport is dominated by edge states and Landau level mixing at the pân interface can occur14,15. Even though this article is focusing on snake states, we will discuss the phenomenology of the mentioned magnetic field regimes. By doing so, we present an integral picture of graphene pân physics in magnetic field.
(a) A two-terminal graphene device consisting of a hole (blue) and an electron (red) cavity is sketched. By applying a weak field, the electron trajectories in the p- and n-cavities bend, leading to dispersing FabryâPérot resonances. (b) The field is increased until the cyclotron orbit becomes comparable to the cavity size, where resonant scar states can occur. (c) The field is further increased and transport is still described by quasiclassical cyclotron orbits. Snake states are formed along the pân interface. (d) Finally, quantum Hall edge states propagate in opposite directions in the p- and n-region at higher fields. (e) Three-dimensional design of the measured device. The SiO2 substrate is coloured in blue and the bottom gates in gold. The contacts, supported by the lift-off resist (green) are coloured in grey. (f) Scanning electron microscope image of a device similar to the measured sample. The graphene is coloured in blue and the bottom gates in gold. Scale bar, 1âμm.
Device architecture
Figure 1e shows the design of the measured device and Fig. 1f a scanning electron microscope picture of a similar sample. In a suspended 2 à 2âμm graphene sheet, a pân junction is formed by applying different voltages on the left (Vleft) and right (Vright) bottom gates, resulting in different charge carrier concentrations nleft and nright. The fabrication follows partly ref. 16, which we combined with a wet transfer process allowing us to align the graphene with the bottom gates. For details, see ref. 17 and Methods.
Measurements in the FabryâPérot and quantum Hall regime
In the following, we characterize the measured device in the zero, low- and high-field regimes. Figure 2a shows a two-dimensional colour map of the electrical conductance G(Vright,Vleft). As soon as a pân interface is formed, G is lowered drastically. Regular FabryâPérot resonances in the left/right cavity are visible as oscillations perpendicular to the zero density line in the left/right cavity (horizontal/vertical white dashed line in Fig. 2a), indicating ballistic transport3,4. The inset in Fig. 2a is a slice along the ppânn diagonal (blue dashed). We estimate the mobility to be at carrier density n=1.1 à 109âcmâ2, which we calculated by a simple parallel plate capacitor model. The mobility is mainly limited by scattering at the contacts3,18.
(a) Two-terminal conductance as a function of left and right gate voltage shows regular FabryâPérot oscillations at zero magnetic field. The inset reveals the narrow Dirac dip along the ppânn diagonal (blue dashed line) from which a mobility of μâ470 à 103âcm2âVâ1âsâ1 is deduced. (b) Cuts along the same diagonal at different magnetic field strengths exhibit the expected quantum Hall plateaus at 2, 6, 10âe2/h. The G=2e2/h plateau is already visible at 60âmT. (c) The colour plot as a function of Vleft and Vright at 200âmT shows quantum Hall plateaus in the unipolar region. (d) The numerical derivative (in arbitrary units) is recorded as a function of gate voltage
and B displaying the dispersion of the FabryâPérot oscillations.
is the magnitude of gate voltage in the situation of antisymmetric charge densities nleft=ânright (red dashed line in a, that is, the npâpn diagonal). The white dashed curve indicates the line along which the cyclotron radius Rc is equal to the cavity length L=0.8âμm, the region Rc<L is darkened and will be discussed in the main text. (e) The measured pattern is reproduced by a tight-binding quantum transport calculation based on the designed geometry of the measured device. The small inset shows a resonant electron trajectory at low magnetic field. Constructive interference occurs if the phase along this trajectory is an integer of 2Ï, leading to the numerical solution of the yellow lines (see Methods).
In Fig. 2b, we show cuts along the ppânn diagonal at different magnetic fields B and obtain quantum Hall plateaus at G=G0 à ν, where ν=2,6,10,... is the filling factor12,13 and G0=e2/h is the conductance quantum. Even at fields as low as 60âmT, the ν=2 plateau is visible. The colourscale map in Fig. 2c taken at 200âmT shows that plateaus develop in the unipolar region with conductance values given by the lowest number of edge modes in the left or right cavity, that is, G=G0 à min(νright,νleft)=2, 6, 10âe2/h. This observation compares well to experiments of ref. 14, even though we apply only 0.2âT instead of 4âT. In the bipolar region, the conductance stays well below 2e2/h due to the smoothness of the pân interface.
In a further step, we study the dispersion of the FabryâPérot interference pattern in low magnetic field10. Figure 2d shows the numerical derivative as a function of B and
, where
represents the magnitude of gate voltage in the situation of antisymmetric charge density (red dashed line in Fig. 2a, that is, the npâpn diagonal). In this configuration, the device consists of two FabryâPérot cavities of equal length Lâ0.8âμm. The darkened region where the cyclotron radius Rc<L will be discussed later. In Fig. 2e, a tight-binding transport calculation is shown, which reproduces the measured interference pattern very well (for details see ref. 19 and Methods). We highlight the quality of the measured graphene and the ability of the simulation to capture the complex oscillation pattern of this micron-sized system in the magnetic field. The dispersion of the FabryâPérot oscillations can be described by bent electron trajectories such as the one sketched in the inset of Fig. 2e (refs 3, 10). The condition for constructive interference is met if the accumulated phase along such a trajectory is a multiple of 2Ï (see âMethodsâ). The yellow lines in Fig. 2e are numerical solutions based on such a condition.
Measurement of snake states
We now discuss the regime where snake states emerge. In Fig. 3a, a snake state at a sharp pân junction is sketched. Consider a charge carrier trajectory starting at the grey cross with momentum k in âx direction. Due to the magnetic field, the trajectory is bent towards the pân interface within the cyclotron diameter . If the trajectory hits the pân interface, the hole will be transmitted to the n side with high probability due to Klein tunnelling20. At the upper edge of the sample, the snake trajectory scatters at the left side, resulting in a current towards the left contact. At lower n, Rc is reduced. As sketched in Fig. 3b, the snake trajectory scatters to the right at the upper edge, resulting in a net current towards the right contact. With this mechanism, one expects conductance oscillations that depend on B and n, and constant conductance along curves where Rc is constant. In Fig. 3c, we display calculated functions n(B) for constant Rc. Snake states occur once the cyclotron diameter 2Rc is smaller than the sample width W (green dashed curve) and can be described by quasiclassical trajectories as long as Rc is larger than the magnetic length
(red dashed line). In the regime W>2Rc>lB, additional parabolic lines show the condition for which the number of oscillations in the snake pattern is fixed and commensurate with W, that is, m à 2Rc=W with m=1,2,3,....
(a) Charge carrier trajectory (white) along a sharp pân junction in perpendicular magnetic field starting at the grey cross, where Rc is the cyclotron radius and n(x) is the electron density. (b) At lower p- and n-density Rc is reduced. In contrast to a, the trajectory results in current flow towards the right contact. (c) Curves of constant Rc as a function of n and B. The continuous lines are given by the condition that the cyclotron diameter 2Rc is commensurate with the sample width W. Snake states occur between the green dashed 2Rc=W and the red dashed Rc=lB line. The black dashed line indicates up to which field transport is dominated by bent FabryâPérot patterns (dark grey area). In the red area, scar states can occur. Below the red dashed line, Rc is smaller than the magnetic length lB and Landau levels start to dominate the transport (light grey area). (d) Conductance as a function of antisymmetric gate tuning and magnetic field is shown in the lower panel, and its derivative dG/dB (in arbitrary units) in the upper panel. Striking lines of high and low conductance with a parabola-like B-dependence can be observed. (e) In a smooth pân junction (here: linear n(x)), the cyclotron orbits become elongated along the y direction due to the additional electric field caused by the density gradient. (f)
at 120âmT. The blue dashed lines in d and f are equivalent. The inset shows lines of constant electric field Ex at the pân interface as a function of Vright,Vleft taken from electrostatic simulations.
In Fig. 3d, we show the measured conductance (bottom) and its numerical derivative dG/dB (top). The measurement exhibits strong oscillations that follow parabola-like curves. We notice that the oscillations occur on a background of strongly decreasing conductance from Gâ6e2/h to G<2e2/h. The steep decrease indicates that the transport becomes dominated by the low-density region close to the pân interface and this happens when 2Rc<W.
In a real pân interface, the density does not sharply jump from the p to the n side but evolves smoothly. An electron trajectory in such a smooth pân interface is sketched in Fig. 3e. Here the density gradient leads to an electric field Ex and its interplay with the perpendicular magnetic field results in E Ã B drift (here along y), leading to elongated cyclotron orbits. The condition Ry=const. can be studied in the measurement of Fig. 3f, where we show at B=120âmT. The measured oscillation pattern follows curves of constant Ex at the pân interface (obtained from electrostatic simulations) as shown in the inset.
Tight-binding theory of snake states
So far, we have seen that the oscillation pattern occurs in the regime where snake states are expected (that is, 2Rc<W) and that the oscillations are related to transport along the pân interface. We now present a quantitative comparison between experiment and theory. Figure 4a shows a quantum transport simulation of and
on the basis of a scalable tight-binding model19 that fully takes into account the device geometry. The simulation compares very well with the measurement shown in Fig. 3d. The parabola-like patterns are reproduced and a similarly steep decrease of conductance is obtained. In Fig. 4b, a slice following the white dashed line in Fig. 4a is shown. The visibility ÎG/G of the oscillations reaches 30% in theory and experiment and is enabled by the strong Klein collimation at the smooth pân interface.
(a) Tight-binding calculation of conductance (bottom) and its numerical derivative dG/dB (top, in arbitrary unit). The parabola-like lines seen in the experiment are well captured. (b)
along the white dashed line in a at B=90âmT. (c) Calculated x component of the local current density distribution, jx, for electrons injected from the left contact with a small d.c.-offset at B=90âmT and
âV (dashed circle in b). The complex resonance pattern in the left cavity consists of âbubblingâ trajectories that do not contribute to conductance. At the pân interface (dashed line), an alternating current pointing to the left (red) and right (blue) is observed. (dâg) Current density profiles at different
corresponding to the circles in b. As a guide to the eye, a snaking trajectory following jx is added. From one conductance maximum to the next (that is, c to g), one snake period is added. The snake state in e and f corresponds to a conductance minimum and c,d and g to a maximum.
Next we apply the KeldyshâGreenâs function method to extract local current density profiles (see Methods) at high and low conductance along this line. In Fig. 4c, we show the x component of the current density jx, taken at (dashed circle in Fig. 4b). The current is injected from the left contact using a small d.c. offset. In the left cavity, a complex resonance pattern appears, given by so-called âbubblingâ trajectories21, which are reflected before reaching the pân interface and do not contribute to current between the contacts. The pattern relevant for transport is located at the pân interface (dashed line). We observe that jx changes sign along the pân interface and that the blue and red regions penetrate the dashed line, indicating that transport is dominated by Klein-collimated snake trajectories. As a guide to the eye, we added a curve in Fig. 4dâg that follows the snake state. This is done for different
values for which G is maximal/minimal at a fixed magnetic field of 90âmT (coloured circles in Fig. 4b). In Fig. 4d, for example, the current density profile corresponds to a conductance maximum where the current points to the right at the upper edge of the pân interface. One period is added by changing
from one maximum to the next. More current density profiles evolving with
at fixed B=90âmT are shown in Supplementary Movie 1. By tracing along one of the parabola-like patterns, the current density profile of the snake state stays constant (an example is given in Supplementary Movie 2).
The conductance oscillates as a function of the ratio W/2Ry, the exact snake period corresponding to 4Ry is, however, difficult to determine using quasiclassical trajectories, since current is injected from many points under various angles resulting in a complex cusp structure similar to what was predicted in refs 22, 23, 24, 25. The excellent agreement between measurement and calculated conductance for which we could determine local current density profiles clearly indicates that the oscillations result from snake state trajectories.
There are additional parabola-like structures at lower magnetic field indicated by arrows in Fig. 4a, these structures are, however, less pronounced in the experimental data. Those resonances occur in a regime where scar states (Fig. 1b) would be expected. The parabola-like behaviour indicates that the states are commensurate with a cavity dimension. In the model, the resonances disappear for non-reflective contacts as expected for scar states.
Discussion
We investigated the magneto-conductance of a ballistic graphene pân junction in different magnetic field regimes. We have observed resonance patterns occurring in the intermediate quasiclassical regime in experiment and theory, which result from the formation of snake states at the pân interface. Among many other possibilities, these states can be used to guide electrons on arbitrary paths with a high efficiency even at very low magnetic fields. This could be used to guide electrons away from sample edges to suppress uncontrolled momentum or spin scattering. The directional scattering at the sample boundaries could be used to implement multi-terminal switches24,26. Furthermore, the similarity between Andreev reflection and Klein tunnelling is stressed in theory27 leading to a correspondence of snake states and Andreev edge states which are of theoretical28 and experimental29 interest. Our work points out that snake states are highly tunable and occur at low fields and that ballistic graphene pân junctions in a magnetic field reveal novel and intriguing phenomena.
Methods
Experimental methods
High-quality graphene is obtained by in situ current annealing30. All the measurements were performed in a variable temperature Helium cryostat with a base temperature of 1.5âK. We measured differential conductance G=dI/dV by standard lock-in technique applying an a.c. voltage of 0.1âmV at 77âHz.
For the quantum Hall data of Fig. 2b,c, we subtracted a contact resistance of 1.2âkΩ.
We extracted the cavity length L used in Fig. 2d,e from the spacing În between resonant FabryâPérot peaks in the bipolar situation and obtained Lâ0.8âμm.
Simulation methods
Real-space Greenâs function method in the tight-binding framework using a scaled graphene Hamiltonian19 is applied to simulate ballistic quantum transport in the present device, taking into account the realistic on-site energy profile obtained by three-dimensional electrostatic simulation for the self-partial capacitances of the bottom gates. All the presented conductance simulations are obtained by calculating the transmission function at zero temperature, combined with the contact resistance 1.2âkΩ. Local current densities are imaged by applying the KeldyshâGreenâs function method in the linear response regime31 on the basis of the same model Hamiltonian used for conductance simulation. At each lattice site n, the bond charge current density Jn=âmenâmâ¹Jnâm⺠is computed, where the summation runs over all the sites m nearest to n, enâm is the unit vector pointing from n to m, and â¹Jnâm⺠is the quantum statistical average of the bond charge current operator Jnâm (ref. 32). After computing for each site, the position-dependent current density profile J(x,y)=[jx(x,y),jy(x,y)] is imaged. In Fig. 4câg, the x component jx(x,y) is shown.
The low-field FabryâPérot interference contours sketched in Fig. 2e are numerically obtained from solving the resonance condition ÎΦ=2jÏ, arising from the path difference between the directly transmitted and twice reflected trajectories within the p cavity as sketched in the inset of Fig. 2e, which is found to be the major interference contribution. For such a simplified model, the phase difference is given by ÎΦ=ΦWKB+ΦAB+ΦBerry+Φ0, where is the kinetic WKB-phase, ΦAB=δA à eB/â is the AharanovâBohm phase due to the flux enclosed by the bent orbit segments,
is the Berry phase, and Φ0=Ï is a constant phase due to reflections off the two pân junction interfaces of the p cavity (smooth at the middle and sharper at the contact side). Here kF is the numerical average of
within the p cavity. The cavity size L is numerically determined and is about half of the flake length L/2=840ânm. The loop area is given by δA=Rc2(ÏâsinâÏ) with Ï=2 arcsinâ(L/2Rc) and Rc=â kF/eB. The form of the Berry phase follows from the consideration of ref. 33 with the critical field estimated by
, where the critical transmission value Tc is a parameter close to one and does not significantly influence the shape of the contours; Tc=0.95 is chosen. The contours sketched in Fig. 2e correspond to j=1,2,â¯,8.
Additional information
How to cite this article: Rickhaus, P. et al. Snake trajectories in ultraclean graphene pân junctions. Nat. Commun. 6:6470 doi: 10.1038/ncomms7470 (2015).
References
van Houten, H. et al. Coherent electron focusing with quantum point contacts in a two-dimensional electron gas. Phys. Rev. B Condens. Matter 39, 8556â8575 (1989) .
Taychatanapat, T., Watanabe, K., Taniguchi, T. & Jarillo-Herrero, P. Electrically tunable transverse magnetic focusing in graphene. Nat. Phys. 9, 225â229 (2013) .
Rickhaus, P. et al. Ballistic interferences in suspended graphene. Nat. Commun. 4, 2342 (2013) .
Grushina, A. L., Ki, D.-K. & Morpurgo., A. F. A ballistic pn junction in suspended graphene with split bottom gates. Appl. Phys. Lett. 102, 223102 (2013) .
Huard, B. et al. Transport measurements across a tunable potential barrier in graphene. Phys. Rev. Lett. 98, 236803 (2007) .
Ãzyilmaz, B. et al. Electronic transport and quantum hall effect in bipolar graphene p-n-p junctions. Phys. Rev. Lett. 99, 166804 (2007) .
Gorbachev, R. V., Mayorov, A. S., Savchenko, A. K., Horsell, D. W. & Guinea, F. Conductance of p-n-p graphene structures with "air-bridge" top gates. Nano Lett. 8, 1995â1999 (2008) .
Ye, P. D., Weiss, D. & Gerhardts, R. R. Electrons in a periodic magnetic field induced by a regular array of micromagnets. Phys. Rev. Lett. 74, 3013 (1995) .
Williams, J. R. & Marcus, C. M. Snake states along graphene p-n junctions. Phys. Rev. Lett. 107, 046602 (2011) .
Young, A. F. & Kim, P. Quantum interference and klein tunnelling in graphene heterojunctions. Nat. Phys. 5, 222â226 (2009) .
Bird, J. et al. Lead-orientation-dependent wave function scarring in open quantum dots. Phys. Rev. Lett. 82, 4691â4694 (1999) .
Novoselov, K. S. et al. Two-dimensional gas of massless dirac fermions in graphene. Nature 438, 197â200 (2005) .
Zhang, Y., Tan, Y.-W., Stormer, H. L. & Kim, P. Experimental observation of the quantum hall effect and berry's phase in graphene. Nature 438, 201 (2005) .
Williams, J. R., Dicarlo, L. & Marcus, C. M. Quantum hall effect in a gate-controlled p-n junction of graphene. Science 317, 638â641 (2007) .
Abanin, D. A. & Levitov, L. S. Quantized transport in graphene p-n junctions in a magnetic field. Science 317, 641â643 (2007) .
Tombros, N. et al. Large yield production of high mobility freely suspended graphene electronic devices on a polydimethylglutarimide based organic polymer. J. Appl. Phys. 109, 093702 (2011) .
Maurand, R. et al. Fabrication of ballistic suspended graphene with local-gating. Carbon 79, 486â492 (2014) .
Cayssol, J., Huard, B. & Goldhaber-Gordon, D. Contact resistance and shot noise in graphene transistors. Phys. Rev. B 79, 075428 (2009) .
Liu, M.-H. et al. Scalable tight-binding model for graphene. Phys. Rev. Lett. 114, 036601 (2015) .
Cheianov, V. & Falko, V. Selective transmission of Dirac electrons and ballistic magneto-resistance of n-p junctions in graphene. Phys. Rev. B 74, 041403 (2006) .
Carmier, P., Lewenkopf, C. & Ullmo, D. Semiclassical magnetotransport in graphene n-p junctions. Phys. Rev. B 84, 195428 (2011) .
Davies, N. et al. Skipping and snake orbits of electrons: singularities and catastrophes. Phys. Rev. B 85, 155433 (2012) .
Patel, A. A., Davies, N., Cheianov, V. & Fal'ko, V. I. Classical and quantum magneto-oscillations of current flow near a p-n junction in graphene. Phys. Rev. B 86, 081413 (2012) .
Milovanovic, S. P., Ramezani Masir, M. & Peeters, F. M. Spectroscopy of snake states using a graphene Hall bar. Appl. Phys. Lett. 103, 233502 (2013) .
Milovanovic, S. P., Ramezani Masir, M. & Peeters, F. M. Magnetic electron focusing and tuning of the electron current with a pn-junction. J. Appl. Phys. 115, 043719 (2014) .
Chen, J.-C., Xie, X. C. & Sun, Q.-f. Current oscillation of snake states in graphene p-n junction. Phys. Rev. B 86, 035429 (2012) .
Beenakker, C., Akhmerov, A., Recher, P. & Tworzydlo, J. Correspondence between Andreev reflection and Klein tunneling in bipolar graphene. Phys. Rev. B 77, 075409 (2008) .
Hoppe, H., Zulicke, U. & Schon, G. Andreev reflection in strong magnetic fields. Phys. Rev. Lett. 84, 1804â1807 (2000) .
Rickhaus, P., Weiss, M., Marot, L. & Schonenberger, C. Quantum hall effect in graphene with superconducting electrodes. Nano Lett. 12, 1942â1945 (2012) .
Moser, J., Barreiro, A. & Bachtold., A. Current-induced cleaning of graphene. Appl. Phys. Lett. 91, 163513 (2007) .
Cresti, A., Grosso, G. & Parravicini, G. Electronic states and magnetotransport in unipolar and bipolar graphene ribbons. Phys. Rev. B 77, 115408 (2008) .
Nikolic, B., Zarbo, L. & Souma, S. Imaging mesoscopic spin Hall flow: spatial distribution of local spin currents and spin densities in and out of multiterminal spin-orbit coupled semiconductor nanostructures. Phys. Rev. B 73, 075303 (2006) .
Shytov, A., Rudner, M. & Levitov., L. Klein backscattering and Fabry-Perot interference in graphene heterojunctions. Phys. Rev. Lett. 101, 156804 (2008) .
Acknowledgements
M.-H.L. and K.R. acknowledge financial support by the Deutsche Forschungsgemeinschaft (SFB 689 and SPP 1666). This work was further funded by the Swiss National Science Foundation, the Swiss Nanoscience Institute, the Swiss NCCR QSIT, the Graphene Flagship, the ESF programme Eurographene, the EU FP7 project SE2ND, the ERC Advanced Investigator Grant QUEST and ERC grant 258789 and OTKA grant K112918. We thank Andreas Baumgartner, Laszlo Oroszlany and Szabolcs Csonka for fruitful discussions.
Author information
Authors and Affiliations
Contributions
P.R., P.M. and E.T. fabricated the devices. The measurements were performed by P.R., P.M., R.M., E.T. and M.W. M.-H.L. performed the simulations. C.S. and K.R. guided the work. All the authors worked on the manuscript.
Corresponding author
Ethics declarations
Competing interests
The authors declare no competing financial interests.
Supplementary information
Supplementary Movie 1
In this movie we show the local current density distribution in x- (jx) and y-direction (jy) for a fixed magnetic field of 90 mT. By tuning the densities in the p- and n- cavities (via \tilde V) the number of nodes of the snake state along the interface is changing, resulting in conductance oscillations. (MOV 9652 kb)
Supplementary Movie 2
In this movie the density and magnetic field are changed such, that the conductance stays constant, i.e. we sweep along one of the parabola-like lines. Even though the pattern in the left cavity is strongly changing, the snake state along the interface remains unchanged. (MOV 8771 kb)
Rights and permissions
This work is licensed under a Creative Commons Attribution 4.0 International License. The images or other third party material in this article are included in the articleâs Creative Commons license, unless indicated otherwise in the credit line; if the material is not included under the Creative Commons license, users will need to obtain permission from the license holder to reproduce the material. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/
About this article
Cite this article
Rickhaus, P., Makk, P., Liu, MH. et al. Snake trajectories in ultraclean graphene pân junctions. Nat Commun 6, 6470 (2015). https://doi.org/10.1038/ncomms7470
Received:
Accepted:
Published:
DOI: https://doi.org/10.1038/ncomms7470
This article is cited by
-
Ballistic transport spectroscopy of spin-orbit-coupled bands in monolayer graphene on WSe2
Nature Communications (2023)
-
Electronic materials with nanoscale curved geometries
Nature Electronics (2022)
-
Upstream modes and antidots poison graphene quantum Hall effect
Nature Communications (2021)
-
Electrostatic superlattices on scaled graphene lattices
Communications Physics (2020)
-
Interplay of filling fraction and coherence in symmetry broken graphene p-n junction
Communications Physics (2020)
Comments
By submitting a comment you agree to abide by our Terms and Community Guidelines. If you find something abusive or that does not comply with our terms or guidelines please flag it as inappropriate.