---
title: Magic-Angle Twisted Bilayer Graphene
url: https://www.emergentmind.com/topics/magic-angle-twisted-bilayer-graphene-tblg
type: topic
---

# Magic-Angle Twisted Bilayer Graphene

Magic-angle twisted bilayer graphene (tBLG) is a two-layer graphene system in which a small relative twist generates a moiré superlattice whose low-energy minibands become exceptionally narrow near the first magic-angle regime around \(1.1^\circ\). In that regime, the Dirac velocity is strongly renormalized, kinetic energy is quenched, and the electronic structure becomes a platform for correlated insulating behavior, superconductivity, nematicity, orbital ferromagnetism, and topological transport. The subject is defined as much by the moiré geometry and its symmetry structure as by the many-body phases built on top of the flat bands [1912.00587; 2008.08129].

## 1. Moiré geometry, flat-band formation, and the meaning of the “magic angle”

Twisted bilayer graphene consists of two graphene monolayers rotated by a relative angle \(\theta\). The moiré period follows
\[
L_m=\frac{a}{2\sin(\theta/2)},
\]
with \(a\approx 0.246\,\text{nm}\), and for \(\theta\ll 1\) in radians this reduces to \(L_m\approx a/\theta\). As \(\theta\) decreases, the moiré real-space cell expands into the tens-of-nanometers range and the mini Brillouin zone contracts correspondingly [1912.00587].

In the Bistritzer–MacDonald picture, twisting shifts the two monolayer Dirac cones in momentum space, while interlayer tunneling hybridizes them. Near the first magic angle, approximately \(1.1^\circ\), the low-energy moiré bands become nearly flat and the Fermi velocity is strongly suppressed. This is the canonical single-particle origin of the strong-correlation regime in tBLG [1912.00587; 2411.00854].

Later ab initio-informed modeling revised the idealized “single magic angle with vanishing velocity” narrative. In the minimal \((2+2)\)-band analysis based on the exact \(k\cdot p\) model with relaxation, pseudogauge fields, and nonlocal interlayer tunneling corrections, the low-energy bands are never perfectly flat and the Fermi velocity never vanishes. Instead, a finite “magic range” appears. Using the shifted angle \(\theta^*=\theta+\theta_{\rm DFT}\) with \(\theta_{\rm DFT}=0.1^\circ\), the lower low-energy band is maximally flat in
\[
\theta^*\in[1.01^\circ,1.14^\circ],
\]
with midpoint features near \(\theta_0^*=1.08^\circ\) [2310.12308].

This more realistic description is tied to atomic relaxation. Below about \(2^\circ\), and especially below a critical \(\theta_c\approx 1.2^\circ\), the moiré texture is better understood as shrunken AA regions, enlarged AB/BA domains, and domain walls separating them. That real-space decomposition is not secondary: it controls which low-energy orbitals dominate and why nearby dispersive domain-wall-derived bands remain important even in the magic-angle regime [2310.12308].

A practical implication is that magic-angle phenomenology is not confined to a single sharply tuned angle. Transport on a device at \(\theta=0.93^\circ\pm0.01^\circ\), about \(15\%\) below the nominal magic angle, still showed a Mott-like correlated insulator and superconductivity, motivating the language of a broader “magic range” rather than a single singular value [1902.05151].

## 2. Direct electronic structure and spectroscopic visualization of the flat bands

A major milestone was the direct momentum-resolved observation of the flat moiré miniband by nanoARPES. In an uncapped tBLG on hBN on doped Si, local ARPES on a domain with
\[
\theta=(0.96\pm0.03)^\circ,\qquad L_m=(14.7\pm0.4)\,\text{nm}
\]
revealed a sharp, weakly dispersing band near \(E_F\) around the graphene \(K\) points at room temperature. The key point was not merely a van Hove singularity in the density of states, but a direct momentum-space visualization of a flat miniband near charge neutrality [1912.00587].

NanoARPES was essential because conventional large-spot ARPES averages over twist-angle inhomogeneity, strain, reconstruction, and disorder. The experiment used a beam spot of about \(1\,\mu\text{m}\), photon energy \(95\,\text{eV}\), net energy resolution about \(34\,\text{meV}\), and UHV base pressure better than \(5\times10^{-11}\,\text{mbar}\). The flat-band signal extended over about \(\sim 0.1\,\text{\AA}^{-1}\) along \(k_x\), which was interpreted as consistent with strong real-space localization of the associated states [1912.00587].

The same measurements showed that the low-energy feature is embedded in a richer reconstructed band structure. At higher binding energy, the spectra exhibit multiple hybridized Dirac cones, moiré-zone repetition, and avoided crossings or gaps produced by the periodic moiré potential. The comparison with spectral-function simulations based on an ab initio-informed tight-binding model at \(1.12^\circ\), including relaxation and band unfolding, reproduced the flat band at \(E_F\), its broad momentum-space extent, small outgoing branches, and multiple hybridized Dirac cones [1912.00587].

Local spectroscopy provided a complementary picture. In STM/STS at the magic angle, a single flat-band DOS peak appears when the band is full or empty, while at partial filling it reconstructs into lower-band and upper-band features separated by a pseudogap. In one study at \(T=4.6\,\text{K}\), the lower-band–upper-band separation near charge neutrality was about \(38\,\text{mV}\), and the \(dI/dV\) signal at \(E_F\) displayed dips near \(\nu=0,\pm1/4,\pm1/2,\pm3/4,\pm1\), linking local spectroscopy directly to the correlated filling sequence [1904.10153].

## 3. Correlated phases: insulating states, superconductivity, and symmetry breaking

The many-body importance of the flat bands appears most clearly in transport. In a device with \(\theta=0.93^\circ\pm0.01^\circ\), the low-energy bandwidth remained \(\sim 11\,\text{meV}\), comparable to magic-angle devices, and the system exhibited a Mott-like correlated insulating state at \(n_m=+2\) together with superconductivity on the electron-doped side for
\[
2.4<n_m<3.1,
\]
with optimal doping near \(n_m\sim 2.7\). The reported superconducting transition scale was \(T_c\sim0.3\)–\(0.5\,\text{K}\), with \(H_{c\parallel}\sim0.5\,\text{T}\), showing that correlated and superconducting behavior survives well below the nominal first magic angle [1902.05151].

That same work also found additional higher-filling structure: narrow activated resistance peaks at \(n_m=\pm5\) with activation gap \(\sim 1.3\,\text{K}\approx0.1\,\text{meV}\), and resistance peaks at \(n_m=\pm12\) attributed to high-energy Dirac points. This broadened the scope of magic-angle phenomenology beyond the first miniband and suggested that partial flatness in higher-energy bands can also generate correlation effects [1902.05151].

A more explicit link between superconductivity and band structure was established by displacement-field tuning in a dual-gated near-magic-angle device at \(\theta\approx0.95^\circ\pm0.02^\circ\). There, the filling factor was defined as
\[
\nu=\frac{4n}{n_s},
\]
and the superconducting dome at zero displacement field occupied approximately \(\nu\approx 2.2\)–\(3\), with optimal doping \(\nu_c\sim2.6\)–\(2.7\). The BKT criterion \(V_{xx}\propto I^\alpha\) with \(\alpha=3\) gave \(T_{BKT}\sim0.6\,\text{K}\), and the perpendicular critical field near the strongest dome was \(B_{c\perp}\approx60\)–\(70\,\text{mT}\) [2402.11649].

The displacement-field dependence showed a pronounced competition between superconductivity and symmetry-broken order. At \(D=0\), superconductivity appeared without a half-filling resistance peak. As \(D\) increased, superconductivity weakened, while a resistance peak emerged near \(\nu\sim+2\) once \(D\gtrsim0.25\,\text{V/nm}\). Hall-density analysis, using
\[
\nu_H\equiv \frac{4n_H}{n_s},
\]
showed that superconductivity at \(D=0\) formed near a van Hove singularity around \(\nu\sim+2.7\) in an isospin-unpolarized \(g_d=4\) regime, whereas larger \(D\) shifted the van Hove singularity toward \(\nu\sim+3.0\) and drove a \(g_d\sim2\) symmetry-broken state near half-filling that suppressed superconductivity [2402.11649].

Local spectroscopy revealed that the correlated normal state is itself ordered. At partial filling of the magic-angle flat band, STM/STS found a pseudogap phase accompanied by a global stripe charge order and broken rotational symmetry. Inside a single moiré cell, the lower-band and upper-band LDOS maps became elliptical with roughly orthogonal principal axes, and the inferred local filling formed a quadrupolar charge pattern. Over larger areas, those quadrupoles aligned into stripe order, lowering the approximate moiré \(C_6\) symmetry to \(C_2\) [1904.10153].

## 4. Topology, symmetry, and low-energy effective descriptions

Flatness is not the whole story. In small-angle tBLG, the low-energy moiré Dirac bands also carry nontrivial topology protected by symmetry. Nonlocal transport in hBN-encapsulated devices with twist angles mainly between about \(1.3^\circ\) and \(1.9^\circ\) revealed pronounced nonlocal resistance peaks in both electron and hole superlattice gaps, while similar responses were absent in lower-angle samples near \(0.75^\circ\) and \(0.42^\circ\) where the relevant superlattice gaps closed. The key symmetry is
\[
(C_{2z}T)\,H(\mathbf{k})\,(C_{2z}T)^{-1}=H(\mathbf{k}),\qquad (C_{2z}T)^2=1,
\]
which trivializes Berry curvature but leaves two \(\mathbb Z_2\) invariants for the isolated two-band moiré Dirac subspace [1903.07950].

One of those invariants is the quantized Berry phase
\[
\gamma\in\{0,\pi\}\pmod{2\pi},
\]
which protects the moiré Dirac points. The other is encoded in the Wilson-loop or Wannier-center counterflow winding
\[
\mathcal W=\mathcal P\exp\!\left(i\oint \mathcal A(\mathbf k)\cdot d\mathbf k\right),
\]
and implies one pair of counter-propagating edge states per spin and valley in each superlattice gap. Experimentally, the resulting nonlocal response followed
\[
R_{NL}\sim R_L^\alpha e^{-L/\lambda},
\]
with \(\lambda_e\approx1.2\,\mu\text{m}\) and \(\lambda_h\approx1.4\,\mu\text{m}\) at \(80\,\text{K}\) [1903.07950].

At the model-building level, realistic tBLG is not captured adequately by a perfectly flat isolated two-band picture. The minimal \((2+2)\)-band model places two AA-derived orbitals and two domain-wall-derived dispersive orbitals on a honeycomb lattice and fits the low-energy bands with 13 physically motivated parameters that vary smoothly across the magic regime. Its exact Schur-complement reduction gives an effective two-band Hamiltonian,
\[
H^{\rm eff}_k(E)=H^\ast_k-H^{\rm int}_k\left(H^\triangle_k-E\right)^{-1}H^{{\rm int},\dagger}_k,
\]
providing a compact starting point for Hubbard-like many-body modeling while retaining the influence of nearby dispersive bands [2310.12308].

First-principles wavefunction calculations sharpened the real-space interpretation. Fully relaxed DFT down to \(0.99^\circ\) identified four characteristic wavefunction textures in the low-energy sector: AA-centered states forming a triangular lattice, ring-like AA\(_z\) states, domain-wall states forming a Kagome lattice, and AB/BA states forming a honeycomb lattice. By tuning interlayer coupling, the calculations tracked the emergence of the flat bands, the associated band inversion, and a further likely topological phase transition at stronger coupling, where the upper and lower flat bands exchange mirror representations and wavefunction character. The transition was associated with \(\Delta z\sim -0.08\) to \(-0.09\) Å, corresponding to roughly \(0.5\)–\(1.0\) GPa [2507.03675].

## 5. Perturbations and band-structure engineering

Magic-angle tBLG is unusually susceptible to perturbations that are modest in other Dirac materials. One route is heterostrain. STM/STS on graphene bilayers grown on \(\alpha\)-Mo\(_2\)C showed that a moderate heterostrain of only about \(\varepsilon\sim0.3\%\) can evolve a small-angle TGB from roughly \(1.5^\circ\) with \(D\approx9.5\,\text{nm}\) into a strained magic-angle TGB near \(1.1^\circ\) with \(D\approx13.5\,\text{nm}\), accompanied by the characteristic merger of two low-energy VHSs into a single flat-band peak. With still larger heterostrain, the system enters highly strained tiny-angle deformed tetragonal superlattices supporting topological helical domain-wall networks and localized domain-wall modes arranged into a hexagon-triangle-mixed frustrated lattice derived from Kagome geometry [1805.03790].

Another route is substrate-induced spin-orbit coupling. In tBLG on a TMD, proximity-induced Rashba, Ising, and sublattice terms,
\[
H_{SOC_I}=\lambda_I\tau^z s^z,\qquad
H_{SOC_R}=\lambda_R(\tau^z\sigma^x s^y-\sigma^y s^x),\qquad
H_{SL}=u\sigma^z,
\]
can reconstruct the eight near-magic-angle flat bands into a topological spin-orbit-coupled moiré manifold. For realistic values \(\lambda_R\approx16\,\text{meV}\), \(\lambda_I\approx1\,\text{meV}\), and \(u\approx1\,\text{meV}\), the theory predicts valley Chern insulating regimes and, in even or one-sided structures, a time-reversal-protected topological insulator at \(\nu=\pm2\) [2008.06528].

Charged defects provide a third form of local engineering. Tight-binding calculations with Coulomb impurities showed that in near-magic-angle tBLG the effect depends strongly on impurity position within the moiré cell. An impurity in the AA region can induce additional flattening of the low-energy bands; in the AB region it can open a gap at the moiré Dirac points by breaking the effective moiré sublattice symmetry; and in a bridge region it removes the Dirac points while leaving the system metallic through a band crossing along \(\Gamma\)-to-\(K\) [2211.01038].

Magnetic fields can also reshape the single-particle structure. In an atomistic calculation at \(\theta=1.08^\circ\), a perpendicular field produced dispersive Landau levels when the magnetic length was comparable to the AA and AB region size, while a strong in-plane field modified the low-energy bands and gap through orbital minimal coupling. In the symmetric gauge, the layer-dependent momentum shift is
\[
\Delta \mathbf k(\mathbf B)=\frac{e}{2\hbar c}\begin{pmatrix} B_y\\ -B_x\end{pmatrix}\Delta z,
\]
and the effect becomes appreciable because the intrinsic moiré momentum scale is already small near the magic angle [2307.09960].

## 6. Debates, limitations, and current directions

Several central issues remain unresolved. One is the microscopic origin of superconductivity. Screening experiments point in different directions. In a device where magic-angle tBLG was placed \(3\,\text{nm}\) from a Bernal bilayer graphene screening layer, increasing screening weakened the correlated insulating states but enhanced superconductivity, with \(T_c\) reaching about \(2.2\,\text{K}\) in the screened configuration; this was interpreted as evidence that electron-phonon coupling is the dominant pairing mechanism and Coulomb repulsion a competing interaction [2003.11072].

A later double-moiré screening experiment reached the opposite conclusion. There, a magic-angle TBG at \(\theta_{\mathrm{bTBG}}=1.15^\circ\) was placed only \(0.35\,\text{nm}\) from an electronically decoupled small-angle TBG screening layer with \(\theta_{\mathrm{tTBG}}=0.46^\circ\). Increasing the screening-layer density completely suppressed both the \(\nu=+2\) correlated insulator and superconductivity in the adjacent magic-angle layer, which was interpreted as strong support for an unconventional, electronically mediated pairing mechanism [2412.01577].

The state dependence of screening helps explain why these results are not trivially incompatible. A separate screening analysis found that interlayer coupling already enhances internal screening in magic-angle tBLG and that external dielectric engineering depends decisively on the electronic state. In a metallic RPA-like state, changing the external dielectric has little effect on the local interaction scale, whereas in a cRPA-like or insulating state it can alter the effective interaction by roughly \(40\)–\(50\%\) [1904.11765]. This suggests that screening experiments probe different interaction channels depending on geometry, layer separation, and the density of states of the screening medium.

A second unresolved issue is how literally the phrase “magic angle” should be taken. The Bistritzer–MacDonald paradigm remains the organizing principle, but realistic modeling and experiments below \(1.1^\circ\) support a broadened practical magic range rather than a single singular angle [2310.12308; 1902.05151]. A third issue concerns topology beyond the original fragile-band narrative. First-principles wavefunction work identified strong indicators of an additional topological transition under pressure or reduced angle, but did not compute topological invariants directly, so the precise classification remains open [2507.03675].

Finally, there are important experimental limitations. Room-temperature nanoARPES directly established the flat-band platform, but it did not observe the low-temperature ordered phases themselves [1912.00587]. Local probes reveal pseudogaps, charge order, and broken rotational symmetry, while transport reveals superconductivity and correlated insulators, yet a unified microscopic theory connecting flat-band topology, displacement-field band renormalization, screening environment, and pairing remains incomplete. What is firmly established is the platform: near the first magic-angle regime, tBLG realizes a moiré electronic structure in which single-particle flatness, topology, and externally tunable interaction scales become inseparable.

Source: https://www.emergentmind.com/topics/magic-angle-twisted-bilayer-graphene-tblg