Programmable Dirac masses in hybrid moiré--1D superlattices

Abstract

Twisted moiré Dirac systems enable powerful miniband engineering but are largely fixed once the twist angle is set, whereas unidirectional (1D) electrostatic superlattices offer continuous control of Dirac anisotropy; yet robust single-particle gaps at charge neutrality are generally difficult to obtain in either setting. Here we show that combining the two into a hybrid moiré--1D superlattice provides a gate-defined configuration space that hosts both gap-opening resonances and strongly anisotropic gapless regimes. Using full-wave continuum miniband calculations for twisted bilayer graphene, we map the charge-neutrality-point (CNP) gap versus the 1D wavevector $\mathbf G_{\rm 1D}$ and identify a Dirac--Dirac resonance condition. At resonance, a single-particle CNP gap emerges from a parity--chirality selection rule for the resonant inter-cone coupling, which can be electrically reprogrammed by layer-asymmetric modulation that switches the relative chirality and the active mass channel. The insulating phase persists within a finite near-resonant window, providing quantitative fabrication tolerances, while off-resonant settings remain gapless but enable strong suppression of the transverse Dirac velocity and continuous anisotropic band renormalization. Hybrid moiré--1D superlattices thus provide a practical route to programmable Dirac minibands and electrically selectable mass channels in coupled Dirac systems.

Figures (4)

• Figure 1: Gate-defined configuration-space phase diagram.a Schematic of a moiré--1D hybrid superlattice: twisted bilayer graphene (red and blue hexagons) under a unidirectional scalar potential (shade) with wavevector $\mathbf G_{\rm 1D}=2\pi \mathbf{L}_\mathrm{1D}/|\mathbf{L}_\mathrm{1D}|^2$. Inset illustrates $\mathbf G_{\rm 1D}$ for several resonant configurations (Eq. \ref{['eq:main_resonance']}) $n_\mathrm{1D}=1,2,3$ in a reciprocal-space. Red and blue hexagons show the BZ of each graphene layer, black hexagons show the rBZ of TBG, and $\Delta\mathbf K_\xi(\theta)$ represents the moiré Dirac-point mismatch. See Fig. S2 for the rBZ of TBG--1D. b Full-wave continuum map of the direct gap at CNP, $\Delta_\mathrm{dir}$ (in units of $\hbar v_\mathrm{F} k_\theta$, log-scale), versus $\mathbf G_{\rm 1D}$ at twist angle $\theta=2^\circ$ for a layer-symmetric (in-phase) modulation ($|u_1|=|u_2|=0.6$ and $\phi_1=\phi_2=0$; see Fig. S8 for $\theta=3^\circ$ and $5^\circ$, and numerical details in Methods and Supplementary Note S2). Resonance markers indicate $\mathbf{G}_\mathrm{1D}$ for $n_{\rm 1D}=1,2,3$, which makes Dirac--Dirac matching, while the red-dashed circle highlights the near-resonant gapped arc around the $n_{\rm 1D}=1$ resonance. This map corresponds to the default relative chirality class ($\chi_{\rm rel}=+1$) of the layer-symmetric slice; chirality switching under layer-asymmetric modulation can reassign which resonances become gapped (Eq. \ref{['eq:main_selection_rule']}, Figs. \ref{['fig:fig3']}a and \ref{['fig:fig3']}b). c DOS at a resonant configuration $n_\mathrm{1D}=1,2,3$ in b and that of a pristine TBG as a reference, illustrating the opening of a gap at CNP when $n_\mathrm{1D}$ is odd. d Band dispersion at a off-resonant configuration ($|\textbf{G}_\mathrm{1D}|=50~\mathrm{nm}$) in b, illustrating the anisotropic band flattening analogous to monolayer graphene under unidirectional scalar potentialPark2008_NatPhys_1DBreyFertig2009_PRLLi2021_NatNano_1DSL but at a weaker potential magnitude by the contribution from the moiré tunnelling (Eq. \ref{['eq:main_vstar']}, Fig. \ref{['fig:fig4']}, Supplementary Note S7)
• Figure 2: Resonant spectra: $n_\mathrm{1D}$ parity selectivity in the default chirality class. (Left) Full-wave continuum minibands near the CNP of three resonant configurations (Eq. \ref{['eq:main_resonance']}) a$n_{\rm 1D}=1$, b$n_{\rm 1D}=2$, c$n_{\rm 1D}=3$, at twist angle $\theta=2^\circ$ under a layer-symmetric 1D modulation ($|u_1|=|u_2|=u$ and $\phi_1=\phi_2=0$) with $u=0.6$. Gray polygon at $E=0$ represents the rBZ of each configuration (see Fig. S2). Gray cone shows the Dirac cone of pristine TBG in the absence of the unidirectional scalar potential for reference. (Right) DOS of each configuration with $u=0,0.2,0.4,0.6$. A robust CNP gap appears for the odd resonance $n_{\rm 1D}=1,3$, while the even resonance $n_{\rm 1D}=2$ remains gapless in this default chirality class ($\chi_{\rm rel}=+1$).
• Figure 3: Parity--chirality selection rule, chirality switching, and gap-opening channel control at Dirac--Dirac resonance. (a) Summary of the resonant gap-opening rule, organized by resonance parity (odd/even $n_{\rm 1D}$) and relative chirality class $\chi_{\rm rel}$ (same/opposite pseudospin winding of the two resonant Dirac cones). The Pauli matrix shown in each chirality column indicates the mass channel that anticommutes with the corresponding effective two-cone Dirac Hamiltonian. See Supplementary Notes S3--S4. (b) Chirality inversion switches the resonance parity for gap opening and reassigns the mass channel: band structures at the even resonance $n_{\rm 1D}=2$ for a strongly layer-asymmetric modulation in which $u_2$ is kept weak-to-moderate while $|u_1|$ is tuned from weak (left, $\chi_{\rm rel}=+1$) to strong (right, $\chi_{\rm rel}=-1$). Upon chirality switching, the gap-opening channel switches from $\sigma_z$ to $\sigma_x$, enabling an even-$n_{\rm 1D}$ CNP gap. The critical $|u_1|$ for chirality inversion in the full-wave spectrum can deviate slightly from simplified two-wave estimates due to additional dressing by remote states. (c) Pauli-channel decomposition of the effective resonant inter-cone coupling $Z\mathcal{M}^{(n_{\rm 1D})}$ for a layer-symmetric modulation (example shown for TBG--1D at $\theta=3^\circ$ with $u_1=u_2=1$ and $\phi_1=\phi_2$). Writing $Z\mathcal{M}^{(n_{\rm 1D})}=m_0 \mathbb{I}+m_x\sigma_x+m_y\sigma_y+m_z\sigma_z$, we plot the normalized magnitudes $|m_\mu|/(Z t_0)$ ($\mu=0,x,y,z$) for $n_{\rm 1D}=1$ to $4$. The resonant coupling exhibits a parity-structured Pauli content: odd $n_{\rm 1D}$ produces a prominent $\sigma_z$ component, while even $n_{\rm 1D}$ favors $\sigma_x$-type coupling, which becomes the mass channel only after chirality switching, consistent with the parity--chirality selection rule.
• Figure 4: Anisotropic band renormalization at off-resonance.a Off-resonant band structures ($|\mathbf{G}_\mathrm{1D}|=50~\mathrm{nm}$, $\theta=2^\circ$), along the direction perpendicular (left) and parallel (middle) to $\textbf{G}_\mathrm{1D}$, and DOS (right) near the CNP as a function of modulation strength $u_l$, illustrating continuously tunable anisotropy and progressive suppression of the transverse Dirac velocity in the Dirac-point neighborhood without opening a gap. (b) Map of the renormalized transverse velocity $v_\perp^*$ versus twist angle $\theta$ and modulation strength $|u_l|$, highlighting the flattening line (red dashed line) where the low-energy transverse velocity approaches zero (Eq. \ref{['eq:main_vstar']}). The $|u_l|$ required for band flattening in TBG--1D is weaker than monolayer graphene under unidirectional potential ($|u_l|\sim 2.83$) due to the contribution from the moiré tunnelling (Eq. \ref{['eq:main_vstar']}).