Abstract
The presence of periodic modulation in graphene leads to a reconstruction of the band structure and formation of minibands. In an external uniform magnetic field, a fractal energy spectrum called Hofstadter butterfly is formed. Particularly interesting in this regard are superlattices with tunable modulation strength, such as electrostatically induced ones in graphene. We perform quantum transport modeling in gate-induced square two-dimensional superlattice in graphene and investigate the relation to the details of the band structure. At low magnetic field the dynamics of carriers reflects the semi-classical orbits which depend on the mini band structure. We theoretically model transverse magnetic focusing, a ballistic transport technique by means of which we investigate the minibands, their extent and carrier type. We find a good agreement between the focusing spectra and the mini band structures obtained from the continuum model, proving usefulness of this technique. At high magnetic field the calculated four-probe resistance fit the Hofstadter butterfly spectrum obtained for our superlattice. Our quantum transport modeling provides an insight into the mini band structures, and can be applied to other superlattice geometries.
Similar content being viewed by others
Introduction
Graphene, a 2D material characterized by a linear low-energy dispersion relation, hosts charge carriers named Dirac fermions due to the resemblance of relativistic (massless) particles described by the Dirac equation. Modifying the underlying graphene lattice by a smooth periodic potential can affect the band structure through folding of the pristine graphene Dirac cone into mini bands1, formation of the secondary Dirac points, and anisotropic renormalization of velocity2,3,4,5,6. Periodic modulation has been obtained in graphene through chemical functionalization7, placing graphene on self-assembled nanostructures8, and by stacking graphene together with aligned or slightly misoriented hexagonal boron nitride (hBN), resulting in periodic moiré modulation which generates hexagonal superlattices (SLs)9,10,11,12,13,14,15. Moiré SLs were also created in low-angle twisted graphene bilayers, followed by van der Waals structures made up of few layers of graphene16,17,18,19,20,21 and other 2D materials22,23. SLs are suitable for the observation of the self-similar energy spectrum called Hofstadter butterfly13,14 and Brown-Zak oscillations24,25 that occur when the magnetic flux through the superlattice unit cell is of the order of the magnetic flux quantum, Ï0â=âh/e, and in pristine 2D crystals require unattainable magnetic fields. Also worth mentioning are the many-body phenomena present in moiré SLs26,27.
Despite their potential, artificial lattices tailored by chemical methods or moiré SLs suffer from inability to tune the strength of the periodic potential. In moiré SLs the period can be tuned to some extent via the rotation angle between the stacked layers, but they are inherent of a hexagonal symmetry. Precise control over the SL geometry, period, and strength is vital for the band structure engineering. The above limitations can be circumvented in electrostatic gate-induced SLs that allow an arbitrary design via the gates geometry, with the gate voltage being a knob for the potential strength. The experimental attempts to create gated SLs in graphene included 1D arrays of metal gates28,29,30, followed by patterned dielectric substrates31,32, and few-layer graphene patterned bottom gates for 1D33,34 and 2D SLs35,36 with down to sub-20 nm periods34 manifesting the flexibility of this approach. With the recent advance in the fabrication techniques, gated SLs with the period of a few tens of nanometers are achievable with good device quality and long electron mean free path32. This enabled observation of commensurability oscillations32,33,36, Hofstadter butterfly34,35 and Brown-Zak oscillations36 at magnetic fields of the order of a few tesla. This is more affordable compared to about 25 T required for graphene/hBN SL, where, due to small lattice constant mismatch of 1.8%, the SL period can reach up to 14 nm for aligned lattices.
While transport in quantizing magnetic field in gated SLs has been thoroughly studied, the intermediate magnetic field regime remained mostly unexplored. In the semiclassical treatment, at low magnetic field fermions undergo cyclotron motion that can be probed in transport measurements via transverse magnetic focusing (TMF). This technique has been used to experimentally probe the band structure in pristine mono-, bi-, and tri-layer graphene37, graphene/WSe2 heterostructures38, and moiré superlattices39,40, and theoretically considered in graphene pn junctions41,42 and 1D gated SLs6.
In this work, we perform a theoretical study of the TMF in 2D gate-induced square SLs, and analyze the relation between the observed TMF spectra and the miniband structure. This is complemented by the investigation of magnetotransport in the quantum Hall regime, where we observe signatures of Hofstadter butterfly, matching the numerical Hofstadter spectrum calculated for our gated SL. Our study is performed using quantum transport calculations for multiterminal structures, considering realistic experimental conditions. To model a realistic geometry, we consider SL induced by a patterned dielectric substrate31,32 with a uniform global backgate underneath the dielectric layer, and the graphene sheet sandwiched between two hBN layers lying on top (Fig. 1a). The hBN/graphene/hBN sandwich is covered by a global top gate. The voltage applied to the back gate Vbg controls the strength of the periodic modulation, while the top gate voltage Vtg is used to tune the carrier density across the SL. We follow the design of Ref. 31, with a square lattice etched in SiO2 substrate, with a lattice period λâ=â35ânm (Fig. 1b). For the ease of the calculation, we use a model function Cbg(x,ây) (see Fig. 1b) that approximates the electrostatically simulated capacitance, obtained previously by some of us43. Previous studies35,44 and modeling of previous experiment on magnetic focusing in monolayer graphene37 show good agreement between experiment and simulation for graphene superlattice devices (see Methods), and the present purely theoretical work can be regarded as a guide for further experimental magneto-transport studies.
Results
No external magnetic field
We first simulate the four-point longitudinal resistance Rxx. We consider a four terminal device shown in the right inset of Fig. 1c, where the system length Lâ=â1152ânm, width Wâ=â385 nm, and the top lead width wâ=â245 nm. With the four leads labeled in the right inset of Fig. 1c, we calculate Rxxâ=âR14,23 (for details see âMethodsâ) and show its dependence on the top and back gate voltages in the main panel of Fig. 1c. As can be seen from the map, the strength of the backgate mostly influences the superlattice modulation.
Figure 1d presents the linecuts of Rxx marked in Fig. 1c with the respective colors. Whereas at Vbgâââ0 only a single Dirac peak is visible, for increasing â£Vbg⣠second and higher order satellite Dirac peaks start to appear, as the periodic modulation gets stronger. At Vbgâ=â±â40 V (linecuts in Fig. 1d), a few higher-order Dirac points are resolved. On the other hand, changing the top gate voltage mostly tunes the carrier density in the device. The left inset of Fig. 1c shows a close-up of the boxed region with 30 Vââ¤âVbgââ¤â45 V andâââ3.5 Vââ¤âVtgââ¤âââ1.6 V, where two sharp lines are visible at both sides of the main Dirac peak, corresponding to the secondary Dirac points, and several fainter lines, corresponding to higher order Dirac points. Similar results based on the same capacitance model function have been reported in43, where two-terminal transport simulations were performed. The character of the bands can be verified in magnetotransport, as shown in the following subsections.
Low magnetic field
For a general band dispersion ϵ(k), the semiclassical equations of motion for an electron are given by45
whereâââe is the electron charge, E is the external electric field, and B the magnetic field. In the presence of constant out-of-plane magnetic field Bâ=â(0,â0,âB) only, one can obtain the relation between the shape of the Fermi contour in the momentum space Îk(t) and the carrier trajectory in real space Îr(t) 45
meaning that the cyclotron orbit is obtained by rotating the orbit in the momentum space by 90° clockwise, as illustrated in Fig. 2b, c. Carriers encircle closed orbits of electron type or hole type in the counterclockwise or clockwise direction, respectively, as determined via the group velocity, vâ=ââkϵ(k)/â, and the equation \(\hslash \dot{{{{\bf{k}}}}}=-eB{{{\bf{v}}}}\times \hat{z}\). In pristine graphene at low energy, cyclotron orbits exhibit a circular shape (Fig. 2b), but the bands formed in systems with SL modulation can be highly distorted from the original, conical shape, and thus, noncircular Fermi contours can be observed (e.g. Fig. 2c).
In a typical device designed for transverse magnetic focusing measurement, an emitter and collector [in Fig. 2a marked as 1 and 2, respectively] are located on the same edge at a center to center distance ââ=âdâ+âw, and other contacts act as absorbers. The current injected from an emitter flows along cyclotron orbits with a radius Rc depending on the magnetic field strength, and at the boundary propagates along skipping orbits. The current can end up in the collector, when the diameter or its multiples match â (Fig. 2b), or otherwise in absorber contacts. In the nonlocal resistance measurement (Fig. 2a), this results in maximum or minimum, respectively, of the resistance Rfâ=âR14,23 (see Methods). For electron-like (hole-like) orbits, the resistance maximum condition can be obtained for positive (negative) B.
In pristine graphene, a typical TMF spectrum as a function of magnetic field and voltage contains two fans of focusing peaks, one for the electron band and the other for the hole band37. In a superlattice, the emerging replicas of the Dirac cone cause a substantial modification of the TMF signal. In the following discussion, we choose Vbgâ=â40 V, such that a few higher-order Dirac peaks are present next to the main Dirac point as seen in Fig. 1c.
For transverse magnetic focusing, we choose the system geometry shown in Fig. 2a, with the distance between the bottom leads dâ=â1200ânm, their widths wâ=â100ânm, the side leads width Wâ=â1636ânm, and the top lead width Lâ=â1792ânm.
Figure 3a shows the Rf map as a function of B and Vtg. One can see several series of fans, with the focusing peaks appearing at Bâ>â0 or Bâ<â0. The map is put together with the Rxx calculated at Bâ=â0, shown in Fig. 3b, where the main and secondary Dirac peaks are seen. The sign change of the focusing peaks in Fig. 3a occurs either at the Dirac points, or van Hove singularities, and for the former, it coincides with the Rxx peaks. This confirms the sign change of the carriers when tuning Vtg, occurring as an effect of the band reconstruction due to the superlattice potential. To understand the result in detail, in Fig. 3câf we plot the miniband structures calculated at Bâ=â0, as described in Ref. 43, and the Fermi contours at Eâ=â0, at selected values of Vtg marked by arrows in Fig. 3a. Additionally, in Fig. 3gâk we show zoomed regions of the Rf map marked with the colored rectangles in Fig. 3a.
In Fig. 3c, at Vtgâ=âââ2.75 V, the Fermi level is located at the C1 subband, and the Fermi contour has a rounded shape (Figs. 2b, 3c). In the semiclassical description, electrons injected from lead 1 (Fig. 2a) with an initial velocity vâ=â(0,âv) and kâ=â(0,âk), in a moderate magnetic field travel along a rounded trajectory, and after half a period, tâ=âT/2 encircle half of the closed orbit. From Eq. (3), this corresponds to traveling a distance equal to the diameter, 2Rcâ=â2âk/eB along the x direction (Fig. 2b). For ââ=â2Rc, it leads to the first maximum of Rf. For smaller Rc, the beam is reflected at the edge, and can flow to the collector when ââ=â4Rc,â6Rc,ââ¦â, giving rise to higher order Rf peaks. In general, we can evaluate the field at which the jth maximum occurs as Bjâ=â2âjk/eâ, jâ=â1,â2,ââ¦â. We find k(Vtg) numerically and plot Bj with dashed lines in Fig. 3(g). The C1 subband cone is withinâââ2.9Vââ²âVtgââ²âââ2.5 V. We find a very good agreement with the Rf signal for up to jâââ7. Higher j are not resolved as the system enters the quantum Hall regime, and semiclassical description of the skipping orbits at the edge no longer applies. Atâââ3.1Vââ²âVtgââ²âââ2.9 V, for the V1 subband, Rf is noisy due to scattering of low-energy carriers by the periodic potential.
When the Fermi level is tuned to the van Hove singularity (at Vtgâââââ2.5 for the electron subband, and Vtgâââââ3.1 for the hole subband), the focusing signal vanishes, and smaller fans reappear. Based on the miniband structures (Fig. 3d), we interpret them as due to focusing of the secondary Dirac cones fermions. For example, at Vtgâ=âââ3.12 V (Fig. 3d) at Eâ=â0 there are tiny Dirac cones around the M and X points of the Brillouin zone.
For Vtgââ²âââ3.2 V and Vtgââ³âââ2.4 V, the miniband structures around the Fermi level get more complex, with many overlapping subbands. Nevertheless, we find ranges of Vtg where a single isolated higher-order Dirac cone is present, giving clear focusing signal [see Fig. 3e for the V2 subband, the corresponding Rf zoom in Fig. 3h, and the zoom in Fig. 3i for the C2 subband]. When there are more overlapping subbands, the signal gets very faint. Nevertheless, one can spot fans that fit well to the hole-like orbit within the V4 subband around the M point, see Fig. 3f at Vtgâ=âââ3.85 V, and zoomed Rf in Fig. 3j. A similar feature is resolved in Fig. 3k for the electron-like orbits.
To further illustrate the relation of the real-space orbits to the subbands, we calculate the current density maps at selected focusing peaks. Figure 4aâc shows the line cuts of Rf at Vtg and B marked by the respective colors in Fig. 3g, h, j. In Fig. 4a for the C1 subband, the first two focusing peaks are marked by â and â². The respective current density maps are presented in Fig. 4d, e, revealing typical current densities found for TMF calculations in pristine graphene46,47. In Fig. 4b the line cut for the V2 subband is shown. For the focusing peaks marked with â and â¿, the current density maps are shown in Fig. 4f, g. The trajectories acquire a shape close to a rectangle, consistent with the Fermi contour (Fig. 3e). For the line cut in the hole-like V4 subband (purple) (Fig. 4c), the current densities of the first two peaks are shown in Fig. 4h, i. The orbits acquire a rhombus shape, matching well the Fermi contour in the corner of the Brillouin zone (Figs. 2c, 3f). The red dashed lines show semiclassical trajectories calculated using the Fermi contours obtained from our band structures and Equation 3 with the contour starting at k for which vxââââE/âkxâ=â0. These semiclassical orbits show similarity with the current density, but the current density is obtained from quantum calculation so they are expected to be similar but not strictly identical, in particular, counter-intuitive patterns may occur at certain resonant conditions (e.g. a vertical blob on top of the rhombus-like pattern in Fig. 4h). Note that in Fig. 4fâg the current density is non-zero in the area between the vertical segments, as it contains contributions from multiple initial kx for which vx is low. The noisy background visible in Fig. 4aâc originates from scattering and resonant states due to SL which are irregular and complex due to the wave-like nature of the carriers.
Although the focusing spectrum is not symmetric with respect to the main Dirac point, it is qualitatively similar in minibands above and below the main Dirac point, except for the noisy signal for low-energy valence subbands. Let us note that the modulation induced by electrostatic gates is a complex function of Vbg and Vtg, and the band structure is significantly modified merely by changing Vtg (Fig. 3câf). This is in contrast to the SLs induced via low-angle twisted hBN substrate or twisted graphene layers, where the top gate only sweeps the carrier density without affecting the periodic modulation, and the band structure remains unchanged. This leads to the shape of the fans in Rf closer to straight lines, unlike Refs. 39,40 that observed ones close to parabolic.
High magnetic field
Now, we turn our attention to the high-field regime, where we expect the Hofstadter spectrum to emerge when the magnetic flux per SL unit cell of area A, Ïâ=âAB, is of the order of the flux quantum Ï0â=âh/e. To simulate longitudinal resistance Rxx and Hall resistance Rxy at the same time, while keeping calculations for the four-point resistances minimized, we consider a 5-terminal Hall bar sketched in Fig. 5a with the actually considered geometric dimensions shown. We compute all the 5âÃâ4â=â20 transmission functions between all pairs among the five leads, and process the data to obtain longitudinal resistance Rxxâ=âR14,23 and Hall resistance Rxyâ=âR14,52 following the Büttiker formalism48,49, according to the lead labels shown in Fig. 5a. For all the following discussions, we convert our Vtg axis to the numerically obtained average carrier density \(\bar{n}\) over the entire lattice, in order for a more transparent presentation of our high-field transport simulations.
Figure 5 b, c show the Hall conductance \({G}_{{{{\rm{H}}}}}={R}_{xy}^{-1}\) as a function of \(\bar{n}\) and magnetic field B, with the back gate voltage fixed at Vbgâ=â0 and Vbgâ=â40 V, respectively. To better highlight the Landau fans arising from the quantum Hall effect of graphene, we limit the color range toâââ40ââ¤âGH/G0ââ¤â+â40 in Fig. 5b, where G0â=âe2/h is the conductance quantum. The same limit of the color range is applied to Fig. 5c in order for a direct comparison. Line cuts at Bâ=â2 T marked by the black line on Fig. 5b and orange line on Fig. 5c are shown and compared in Fig. 5d. The lowest few quantum Hall conductance plateaus GHâ=â±â2G0,â±â6G0,â±â10G0,ââ¯âcan be clearly seen in the Vbgâ=â0 case free of superlattice potential (black line), while the combined strong magnetic field (Bâ=â2 T) and strong superlattice modulation (Vbgâ=â40 V) result in a more complex transport feature, in accordance with the predictions50,51 for a periodic 2D modulation, which can be better understood by checking Rxx, instead of Rxy (or GH), to be discussed soon below.
Without taking the inverse of Rxy, Fig. 5e shows the Hall resistance at Bâ=â0.2T with Vbgâ=â0 [red, corresponding to the \(\bar{n}\) range marked on Fig. 5b] and Vbgâ=â40 V [cyan, corresponding to the \(\bar{n}\) range marked on Fig. 5c]. The former (without superlattice) shows a single sign change at \(\bar{n}=0\), typical for graphene, while the latter (with superlattice) shows multiple sign changes at positive and negative \(\bar{n}\), in addition to the main Dirac point at \(\bar{n}=0\), consistent with our previous low-field magnetotransport simulations discussed above. Considering the same range of \(\bar{{\rm{n}}}\) and B as Fig. 5b and c, Fig. 5f shows the longitudinal resistance Rxx with the back gate voltage fixed at Vbgâ=â40 V. Since the period of our gate-controlled superlattice is λâ=â35 nm, and hence the square SL unit cell area Aâ=âλ2â=â1225 nm2, the condition Ï/Ï0â=âAB/(h/e)â=â1 is reached at Bâââ3.376 Tââ¡âB1, which is in sharp contrast with the graphene/hBN moiré superlattice that requires Bâââ24 T in order to reach Ï/Ï0â=â1, because of its periodicity limited to λââ²â14 nm and hence the area \(A=\sqrt{3}{\lambda }^{2}/2\lesssim 170\,{{{{\rm{nm}}}}}^{2}\). As is visible on Fig. 5f, the seemingly complicated Rxx map does exhibit certain features at B1 and B1/2, and vaguely at B1/3 (marked by black triangles), corresponding to Ï/Ï0â=â1,â1/2,â1/3, respectively, consistent with our calculation of the magnetic energy subbands shown in Fig. 5g, where the vertical axis of Ï/Ï0 corresponds to exactly the same B range as Fig. 5f, as well as the average density range. Good consistency between the Rxx map of Fig. 5f and the Hofstadter butterfly shown in Fig. 5g can be seen. For methods adopted to calculate Fig. 5g, see âMethodsâ.
Discussion
In summary, we theoretically investigated transport in gated superlattices based on monolayer graphene. Our zero and low magnetic field transport calculations remain in a good agreement with the continuum model band structure calculated in presence of periodic modulation. We explored the potential of TMF for probing the intricate band structure of graphene with periodic modulation. It offers possibilities to study a plethora of phenomena in superlattices, and opens the door for studies of strongly correlated systems in twisted bilayer graphene52 or in bilayer graphene superlattices53. By exploring the reconstructed band structure via magnetotransport calculations it is possible to engineer devices relying on directed electron flow due to the distortion of Fermi contour, as well as for other applications based on mini band electron optics. We also obtained high-magnetic-field response consistent with the Hofstadter spectrum calculated for a gated superlattice as a function of the gate voltage. Our modeling can be generalized to other superlattice geometries, and is promising for the investigations of future band structure engineered devices working in a broad range of magnetic fields.
Methods
Gated superlattice model
To model a realistic geometry, we consider SL induced by a patterned dielectric substrate31,32 with a uniform global backgate underneath the dielectric layer, and the graphene sheet sandwiched between two hBN layers lying on top [Fig. 1a of the main text]. The hBN/graphene/hBN sandwich is covered by a global top gate. The voltage applied to the back gate Vbg controls the strength of the periodic modulation, while the top gate voltage Vtg is used to tune the carrier density across the SL. We follow the design of Ref. 31, with a square lattice etched in SiO2 substrate, with a lattice period λâ=â35 nm. We use a model function Cbg(x,ây)
where dsmoothâ=â7.5 nm is the smoothness of the modulation, Cout/eâ=â0.77âÃâ1011âcmâ2Vâ1 and Cin/eâ=â0.22âÃâ1011âcmâ2Vâ1, and a%b means the remainder after dividing a by b. The top gate capacitance is assumed to be Ctg/eâ=â1âÃâ1012âcmâ2Vâ1. Using the parallel capacitor model, this roughly corresponds to the top hBN thickness dtâââ16.6 nm, from Ctg/eâ=âε0εhBN/edt, where ε0 is the vacuum permittivity, εhBNâ=â3 is the dielectric constant of hBN, andâââe is the electron charge. Figure 1b of the main text shows the profile of the above model function (5).
For the dual-gated graphene sample free from intrinsic doping, the carrier density is given by
Assuming that the carrier energy in graphene is given by Eâ=â±ââvFk, where â is the reduced Planck constant, vFâââ106âm/s is the Fermi velocity of graphene, and using \(\hslash {v}_{{{{\rm{F}}}}}\approx 3\sqrt{3}/8\,{{{\rm{eV}}}}\,{{{\rm{nm}}}}\), the onsite potential energy can be obtained from
in order to set the global Fermi energy at Eâ=â0 where transport occurs.
Transport calculation
For transport calculation, we use the tight-binding Hamiltonian
where ci (\({c}_{i}^{{\dagger} }\)) is an annihilation (creation) operator of an electron on site i located at riâ=â(xi,âyi). The first sum contains the nearest-neighbor hoppings with the hopping parameter tij, and the second sum describes the onsite potential energy profile. In the presence of an external magnetic field Bâ=â(0,â0,âB), the hopping integral is modified by \({t}_{ij}\to {t}_{ij}\exp (i\phi )\), where the Peierls phase \(\phi =(-e/\hslash )\int\nolimits_{{{{{\bf{r}}}}}_{i}}^{{{{{\bf{r}}}}}_{j}}{{{\bf{A}}}}\cdot d{{{\bf{r}}}}\), with A being the vector potential that satisfiesâââÃâAâ=âB, and the integral going from the site at ri to the site at rj. For a feasible simulation of realistic devices, we use the scalable tight-binding model54, where the hopping parameter becomes tijâ=ât0/sf, and the lattice spacing aâ=âa0sf, sf is the scaling factor, and we use t0â=âââ3 eV and \({a}_{0}=1/4\sqrt{3}\) nm. Transport calculations based on Hamiltonian (8) are done within the wave-function matching for the TMF, and real-space Greenâs function method in other cases, at the global Fermi energy Eâ=â0 and zero temperature. The conductance from lead i to lead j is obtained from the Landauer formula Gjiâ=â2(e2/h)Tji, where the transmission probability Tji is evaluated as a sum over the propagating modes \({T}_{ji}=\mathop{\sum}\limits_{q}{T}_{ji}^{q}\), and
Here, \({t}_{ij}^{pq}\) is the probability amplitude for the transfer from the incoming mode p in lead i to the outgoing mode q in lead j.
In the multiterminal devices, we solve the transport problem for each lead as an input, and build the conductance matrix \({{{\boldsymbol{{{{\mathcal{G}}}}}}}}\)48,49 which relates the current Ii fed to the system in lead i to the voltage at j-th lead Vj through \({I}_{i}=\mathop{\sum }\nolimits_{j = 1}^{N}{{{{\mathcal{G}}}}}_{ij}{V}_{j}\). For an N-terminal system, the matrix elements are
We set the voltage at l-th lead equal to zero, and eliminate the l-th row and column of the matrix. The reduced (Nâââ1)âÃâ(Nâââ1) matrix \(\bar{{{{\boldsymbol{{{{\mathcal{G}}}}}}}}}\) can be inverted to get \({{{\boldsymbol{{{{\mathcal{R}}}}}}}}={\bar{{{{\boldsymbol{{{{\mathcal{G}}}}}}}}}}^{-1}\), where the \({{{\boldsymbol{{{{\mathcal{R}}}}}}}}\) matrix satisfies
With the elements of matrix \({{{\boldsymbol{{{{\mathcal{R}}}}}}}}\), one can evaluate the resistance
with the current flowing from lead k to lead l, zero current in other terminals, and voltage drop measured between leads m and n.
Transverse magnetic focusing
As a numerical example of applying the above outlined Landauer-Büttiker formalism for computing the four-point resistance Eq. (13), we revisit the first TMF experiment on graphene37, considering the same probe spacing (ââ=â500ânm) and width (100ânm) but slightly different geometry of the scattering region (total length 700ânm and width 500ânm) for a 6-terminal Hall bar, made of a graphene lattice scaled by sfâ=â12. Choosing the same configuration of the leads for injector, ground, and voltage probes as the revisited experiment, the computed R61,54 as a function of the external magnetic field B perpendicular to the plane of graphene and the carrier density n is reported in Fig. 6a, showing a map rather consistent with the experiment. Due to the isotropic low-energy dispersion of graphene, the resulting cyclotron trajectory is a simple circle of radius Rcâ=ââk/eB, which is simplified from Eq. (3). By requiring the probe spacing to be equal to an integer times the cyclotron diameter, ââ=âjââ â2Rc, jâ=â1,â2,ââ¯â, together with \(k=\sqrt{\pi | n| }\) for graphene, one can solve for carrier density corresponding to the jth peak of the TMF on the B-n map of Fig. 6a:
The dashed lines on Fig. 6a show n1,ân2,ân3,ân4, matching very well with the patterns of the simulated R61,54, which requires totally 6âÃâ5â=â30 transmission functions for such a 6-terminal device, as explained above. Figure 6b, c shows two exemplary maps of transmission functions, which can look generally very different from the resulting four-point resistance.
Choosing the gauge
In the presence of an external magnetic field and semi-infinite leads, the vector potential must satisfy the translational invariance of the leads. For the magnetic field along the \(\hat{z}\) axis, the most common choice is the Landau gauge, Aâ=â(âââyB,â0,â0) or Aâ=â(0,âxB,â0) for the lead which is translationally invariant along the x or y direction, respectively. For other lead orientation, in general, the proper gauge is different. Therefore, in a multi-terminal device, the required vector potential is not uniform in the entire space. This is not a problem since adding an arbitrary curl-free component to the vector potential does not change the magnetic field. Here, we use the approach introduced in55. Assuming that the proper gauge in the 1st lead is A1(r), for another lead that is at an angle θn with respect to lead 1, the gauge can be chosen as
with
The addition of a gradient of a scalar function does not influence the requirementâââÃâAâ=âB. As an illustration of the transformation, consider A1(r)â=â(âââyB,â0,â0), and θnâ=â90°. Then, fn(r)â=âBxy,âââfn(r)â=â(By,âBx,â0), and An(r)â=â(âââyB,â0,â0)â+â(yB,âxB,â0)â=â(0,âxB,â0).
Applying the transformation (15) so that it only affects lead n is possible by defining a smooth step function ζn(r) which is only nonzero in the translationally invariant part of lead n
Then, in (15) we substitute fn(r)âââζn(r)fn(r) for lead n. In general, for the entire system we define
this completes our gauge transformation. Adopting the vector potential
we have Bâ=âââÃâA everywhere in the system, and the translation invariance in each lead is guaranteed. Importantly, curl of (19) gives exactly the desired B, regardless of the smoothness of the ζn function.
As an example, for the system used for the TMF modeling (Fig. 2a), in the vertical leads we choose the same gauge Aâ=â(0,âxB) with \(\zeta ({{{\bf{r}}}})={(\exp (-(y-{y}_{1})/{d}_{{{{\rm{step}}}}})+1)}^{-1}+{(\exp (-({y}_{2}-y)/{d}_{{{{\rm{step}}}}})+1)}^{-1}\), y1â=â1607ânm, y2â=âââ49ânm, and dstepâ=â2 nm. In the rest of the system, Aâ=â(âââyB,â0) is used. The resulting Ax and Ay profiles are shown in Fig. 7.
Hofstadter butterfly calculation
For the calculation of Hofstadter butterfly one has to consider a magnetic unit cell whose length is the least common multiple of the lattice periodicity and the periodicity introduced by the Peierls phase. For graphene, it contains more than hundreds of thousands of carbon atoms when the magnetic field strength is smaller than 1 T. However, the calculation is greatly simplified by considering a graphene ribbon. For an armchair ribbon with translational invariance along the x direction and finite width along the y direction, in the presence of a perpendicular magnetic field, dispersionless Landau levels appear near kxâ=â0, and the dispersive edge states show up at larger kx. Calculating \({E}_{{k}_{x} = 0}\) as a function of magnetic field, we get the Hofstadter butterfly of graphene. Because of the finite width of the ribbon, the spectrum can also contain edge states. With the increase of B the Landau levels elongate, and at some B the edge states are pushed to kxâ=â0, which results in the appearance of the states in the gaps. To lower the computational burden, we use the scalable tight-binding model54 with sfâ~â8 to calculate \({E}_{{k}_{x} = 0}\) as a function of magnetic field for an armchair ribbon with periodic length along the x axis equal to the superlattice period (λ). Here, in order to ensure superlattice period equal to a multiple of 3a, sf is not an integer, and the ribbon width is larger than 20λ to show the superlattice effect.
In transport, only the states at the Fermi level contribute to the conductance. Therefore, we calculate Hofstadter butterfly spectra for all Vtg values and collect the Fermi states to construct the gate-dependent Hofstadter butterfly spectrum to compare with Rxx and GH obtained from the transport calculation. Note that the spectrum in Fig. 5(g) contains energy levels that appear across the gaps, which are an artifact of the method due to the finite width of the ribbon. As mentioned above, they appear since at some value of B the edge states are pushed to kxâ=â0.
Data availability
Additional data related to this paper may be requested from the authors.
Code availability
Code may be requested from the authors upon reasonable request.
References
Wallbank, J. R., Patel, A. A., Mucha-KruczyÅski, M., Geim, A. K. & Falâko, V. I. Generic miniband structure of graphene on a hexagonal substrate. Phys. Rev. B 87, 245408 (2013).
Park, C.-H., Yang, L., Son, Y.-W., Cohen, M. L. & Louie, S. G. Anisotropic behaviours of massless Dirac fermions in graphene under periodic potentials. Nat. Phys. 4, 213â217 (2008).
Park, C.-H., Yang, L., Son, Y.-W., Cohen, M. L. & Louie, S. G. New generation of massless dirac fermions in graphene under external periodic potentials. Phys. Rev. Lett. 101, 126804 (2008).
Brey, L. & Fertig, H. A. Emerging zero modes for graphene in a periodic potential. Phys. Rev. Lett. 103, 046809 (2009).
Barbier, M., Vasilopoulos, P. & Peeters, F. M. Extra Dirac points in the energy spectrum for superlattices on single-layer graphene. Phys. Rev. B 81, 075438 (2010).
Kang, W.-H., Chen, S.-C. & Liu, M.-H. Cloning of zero modes in one-dimensional graphene superlattices. Phys. Rev. B 102, 195432 (2020).
Sun, Z. et al. Towards hybrid superlattices in graphene. Nat. Commun. 2, 559 (2011).
Zhang, Y., Kim, Y., Gilbert, M. J. & Mason, N. Electronic transport in a two-dimensional superlattice engineered via self-assembled nanostructures. npj 2D Mater. Appl. 2, 31 (2018).
Xue, J. et al. Scanning tunnelling microscopy and spectroscopy of ultra-flat graphene on hexagonal boron nitride. Nat. Mater. 10, 282â285 (2011).
Decker, R. et al. Local electronic properties of graphene on a BN substrate via scanning tunneling microscopy. Nano Lett. 11, 2291â2295 (2011).
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).
Hunt, B. et al. Massive Dirac Fermions and Hofstadter Butterfly in a van der Waals Heterostructure. Science 340, 1427â1430 (2013).
Dean, C. R. et al. Hofstadterâs butterfly and the fractal quantum Hall effect in moiré superlattices. Nature 497, 598â602 (2013).
Yu, G. L. et al. Hierarchy of Hofstadter states and replica quantum Hall ferromagnetism in graphene superlattices. Nat. Phys. 10, 525â529 (2014).
Burg, G. W. et al. Correlated insulating states in twisted double bilayer graphene. Phys. Rev. Lett. 123, 197702 (2019).
Shen, C. et al. Correlated states in twisted double bilayer graphene. Nat. Phys. 16, 520â525 (2020).
Liu, X. et al. Tunable spin-polarized correlated states in twisted double bilayer graphene. Nature 583, 221â225 (2020).
Lin, F. et al. Heteromoiré engineering on magnetic bloch transport in twisted graphene superlattices. Nano Lett. 20, 7572â7579 (2020).
de Vries, F. K. et al. Combined minivalley and layer control in twisted double bilayer graphene. Phys. Rev. Lett. 125, 176801 (2020).
Rickhaus, P. et al. Correlated electron-hole state in twisted double-bilayer graphene. Science 373, 1257â1260 (2021).
Wang, L. et al. New generation of Moiré Superlattices in doubly aligned hBN/Graphene/hBN heterostructures. Nano Lett. 19, 2371â2376 (2019).
Wang, Z. et al. Composite super-moiré lattices in double-aligned graphene heterostructures. Sci. Adv. 5, eaay8897 (2019).
Kumar, R. K. et al. High-temperature quantum oscillations caused by recurring Bloch states in graphene superlattices. Science 357, 181â184 (2017).
Barrier, J. et al. Long-range ballistic transport of Brown-Zak fermions in graphene superlattices. Nat. Commun. 11, 5756 (2020).
Wang, L. et al. Evidence for a fractional fractal quantum Hall effect in graphene superlattices. Science 350, 1231â1234 (2015).
Andrews, B. & Soluyanov, A. Fractional quantum Hall states for moiré superstructures in the Hofstadter regime. Phys. Rev. B 101, 235312 (2020).
Dubey, S. et al. Tunable superlattice in graphene to control the number of dirac points. Nano Lett. 13, 3990â3995 (2013).
Drienovsky, M. et al. Towards superlattices: Lateral bipolar multibarriers in graphene. Phys. Rev. B 89, 115421 (2014).
Kuiri, M., Gupta, G. K., Ronen, Y., Das, T. & Das, A. Large Landau-level splitting in a tunable one-dimensional graphene superlattice probed by magnetocapacitance measurements. Phys. Rev. B 98, 035418 (2018).
Forsythe, C. et al. Band structure engineering of 2D materials using patterned dielectric superlattices. Nat. Nanotechnol. 13, 566â571 (2018).
Li, Y. et al. Anisotropic band flattening in graphene with one-dimensional superlattices. Nat. Nanotechnol. 16, 525â530 (2021).
Drienovsky, M. et al. Commensurability oscillations in one-dimensional graphene superlattices. Phys. Rev. Lett. 121, 026806 (2018).
Barcons Ruiz, D. et al. Engineering high quality graphene superlattices via ion milled ultra-thin etching masks. Nat. Commun. 13, 6926 (2022).
Huber, R. et al. Gate-tunable two-dimensional superlattices in graphene. Nano Lett. 20, 8046â8052 (2020).
Huber, R. et al. Band conductivity oscillations in a gate-tunable graphene superlattice. Nat. Commun. 13, 2856 (2022).
Taychatanapat, T., Watanabe, K., Taniguchi, T. & Jarillo-Herrero, P. Electrically tunable transverse magnetic focusing in graphene. Nat. Phys. 9, 225â229 (2013).
Rao, Q. et al. Ballistic transport spectroscopy of spin-orbit-coupled bands in monolayer graphene on WSe2. Preprint at https://arxiv.org/abs/2303.01018 (2023).
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).
MilovanoviÄ, 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, S. et al. Electron optics with p-n junctions in ballistic graphene. Science 353, 1522â1525 (2016).
Chen, S.-C., Kraft, R., Danneau, R., Richter, K. & Liu, M.-H. Electrostatic superlattices on scaled graphene lattices. Commun. Phys. 3, 71 (2020).
Kraft, R. et al. Anomalous Cyclotron Motion in Graphene Superlattice Cavities. Phys. Rev. Lett. 125, 217701 (2020).
Ashcroft, N. W. & Mermin, N. D. Solid State Physics (Saunders College, Philadelphia, 1976).
Stegmann, T. & Lorke, A. Edge magnetotransport in graphene: A combined analytical and numerical study. Ann. Phys. 527, 723â736 (2015).
PetroviÄ, M. D., MilovanoviÄ, S. P. & Peeters, F. M. Scanning gate microscopy of magnetic focusing in graphene devices: quantum versus classical simulation. Nanotechnology 28, 185202 (2017).
Büttiker, M. Four-terminal phase-coherent conductance. Phys. Rev. Lett. 57, 1761â1764 (1986).
Datta, S. Electronic transport in mesoscopic systems (Cambridge University Press, Cambridge, 1995).
Thouless, D. J., Kohmoto, M., Nightingale, M. P. & den Nijs, M. Quantized hall conductance in a two-dimensional periodic potential. Phys. Rev. Lett. 49, 405â408 (1982).
Streda, P. Quantised Hall effect in a two-dimensional periodic potential. J. Phys. C: Solid State Phys. 15, L1299 (1982).
de Vries, F. K. et al. Gate-defined Josephson junctions in magic-angle twisted bilayer graphene. Nat. Nanotechnol. 16, 760â763 (2021).
Krix, Z. E. & Sushkov, O. P. Patterned bilayer graphene as a tunable strongly correlated system. Phys. Rev. B 107, 165158 (2023).
Liu, M.-H. et al. Scalable tight-binding model for graphene. Phys. Rev. Lett. 114, 036601 (2015).
Baranger, H. U. & Stone, A. D. Electrical linear-response theory in an arbitrary magnetic field: A new Fermi-surface formation. Phys. Rev. B 40, 8169â8193 (1989).
Acknowledgements
We thank 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 supports and National Center for High-performance Computing (NCHC) for providing computational and storage resources. This research was supported in part by PL-Grid Infrastructure, and by the program, Excellence Initiative â Research Universityâ for the AGH University of Science and Technology.
Author information
Authors and Affiliations
Contributions
A.M.K. and M.-H.L. performed transport calculations and wrote the manuscript with the input from all authors. S.-C.C. calculated the mini-band structures within the continuum model and the Hofstadter butterfly within the tight-binding model. M.-H.L. guided the project.
Corresponding authors
Ethics declarations
Competing interests
The authors declare no competing interests.
Additional information
Publisherâs note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
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
MreÅca-KolasiÅska, A., Chen, SC. & Liu, MH. Probing miniband structure and Hofstadter butterfly in gated graphene superlattices via magnetotransport. npj 2D Mater Appl 7, 64 (2023). https://doi.org/10.1038/s41699-023-00426-9
Received:
Accepted:
Published:
DOI: https://doi.org/10.1038/s41699-023-00426-9