Twisted Bilayer Graphene

Twisted bilayer graphene (TBG) consists of two graphene monolayers rotated by a small relative angle. Near the first magic-angle regime, the moiré hybridization strongly suppresses the velocity and bandwidth of bands near charge neutrality. Coulomb, phonon, strain, disorder, and substrate energy scales can then compete on unusually equal terms.

That compact description contains claims of very different status:

Claim	Present status	What establishes it
a small twist produces a long-period moiré pattern	standard	calibrated real- and reciprocal-space structure
near-magic devices possess narrow active bands	established, sample dependent	tunneling, photoemission, compressibility, and transport tied to local angle
interactions reconstruct the active bands	established	chemical-potential resets, spectral-weight transfer, flavor-resolved Landau fans, and thermodynamic anomalies
every integer-filling resistance peak is a Mott insulator	not established	requires excluding one-body gaps, flavor order, topology, localization, and percolation
superconductivity occurs in TBG	established in a subset of near-magic devices	zero resistance together with critical-current, field, phase-stiffness, or spectroscopic evidence
the pairing mechanism and order-parameter symmetry are known	unresolved	no single microscopic account explains the full device-dependent record
TBG can realize Chern and quantum anomalous Hall states	established in selected regimes	hysteretic orbital magnetism and, in the strongest cases, quantized Hall response with vanishing longitudinal resistance

The durable lesson is not that one twist angle guarantees one phase diagram. It is that TBG provides a controlled setting in which geometry reshapes a Dirac Hamiltonian and makes spin, valley, sublattice, lattice relaxation, interactions, phonons, and topology experimentally comparable.

Canonical Scope

This page is the canonical home for TBG-specific theory and evidence: the Bistritzer–MacDonald continuum model, the first magic-angle regime, active-band state counting, the observed correlated and superconducting phases, valley-resolved topology, and the open questions that connect them.

Moiré Superlattices owns generic twist and mismatch geometry, moiré periods and cell areas, mini Brillouin zones, miniband construction, filling conversion, screening, and interaction-versus-kinetic scale estimates. Graphene owns the monolayer Dirac basis and perturbation ledger. van der Waals Heterostructures owns assembly, alignment metrology, gates, contacts, and interface evidence.

Flat Bands owns dispersionless-band mechanisms, projector geometry, Chern-band constraints, and flat-band ferromagnetism across platforms. Correlated Insulators in Moiré Systems owns the cross-platform distinction among Mott-like, generalized Wigner, flavor-ordered, topological, and disorder-driven insulating mechanisms. Moiré Superconductivity owns the cross-platform standards for phase coherence, pairing structure, tuning, and device probes, while Moiré Topology owns topological minibands and integer and fractional Hall classifications. Here those ideas appear only as needed to interpret TBG.

Twist Angle and Structure

Angle is a field, not only a fabrication setting

In the idealized construction, layer 1 is rotated by $+\theta/2$ and layer 2 by $-\theta/2$. Their Dirac points in a fixed valley are separated by

$$\begin{aligned} k_\theta &= 2K\sin\!\left(\frac{|\theta|}{2}\right), \\ K &= \frac{4\pi}{3a}, \end{aligned}$$

where $a\simeq0.246\,\mathrm{nm}$ is graphene’s lattice constant. The corresponding moiré period and cell area are derived at Moiré Superlattices. Around $\theta=1.1^\circ$, they are approximately

$$L_M\simeq12.8\,\mathrm{nm}, \qquad A_M\simeq142\,\mathrm{nm}^2.$$

A real device has a displacement field $\mathbf d(\mathbf r)$ rather than one global angle. Its symmetric gradient contains heterostrain, while the antisymmetric part changes the local rotation:

$$\begin{aligned} \varepsilon_{ij}^{\mathrm{rel}} &= \frac{1}{2} \left( \partial_i d_j+\partial_j d_i \right), \\ \delta\theta(\mathbf r) &= \frac{1}{2} \left( \partial_x d_y-\partial_y d_x \right). \end{aligned}$$

Raman spectroscopy, scanning probes, diffraction, and quantum oscillations weight this inhomogeneous field differently. A transport-inferred angle is therefore not automatically the local angle in the region controlling a tunneling spectrum or a weak superconducting link.

Interlayer adhesion also reconstructs the nominally rigid pattern. Near small twist, low-energy AB and BA regions expand, AA regions shrink, and strain concentrates near domain walls. In continuum models this commonly reduces the effective AA-like tunneling amplitude $w_0$ relative to the AB/BA-like amplitude $w_1$. Corrugation, pressure, dielectric environment, and microscopic parameterization all affect those numbers.

Filling conventions must be stated

Let $n$ be carrier density measured from charge neutrality. The dimensionless filling used most often in TBG is

$$\nu=nA_M.$$

With spin and valley unresolved, each isolated conduction or valence active band accommodates four carriers per moiré cell. Thus

$$\begin{aligned} \nu&=0 &&\text{at charge neutrality}, \\ \nu&=\pm4 &&\text{at the active-band edges}. \end{aligned}$$

The conduction and valence active bands together contain eight one-electron states per cell: two band labels times two spins times two valleys. Describing $\nu=\pm2$ as “half filling” can be convenient, but it is ambiguous unless the reference band and sign convention are named. Density offsets, trapped charge, quantum capacitance, and spatially varying $A_M$ must be included before an integer label is treated as a microscopic fact.

Continuum Hamiltonian

One valley couples three nearby momenta

Smooth moiré tunneling transfers momenta much smaller than the separation between graphene’s two valleys. To leading order, the valleys $\xi=\pm1$ can therefore be modeled independently, while spin is a spectator in the one-electron Hamiltonian. In the layer-sublattice basis, one common continuum convention is

$$\mathcal H_\xi(\mathbf r) = \begin{pmatrix} h_{\xi,+\theta/2} & T_\xi(\mathbf r) \\ T_\xi^\dagger(\mathbf r) & h_{\xi,-\theta/2} \end{pmatrix}.$$

For $\mathbf p=-i\hbar\nabla$, the rotated monolayer block is

$$\begin{aligned} h_{\xi,\phi} &= v_F \left[ \xi\sigma_x \left(R_{-\phi}\mathbf p\right)_x \right. \\ &\qquad\left. + \sigma_y \left(R_{-\phi}\mathbf p\right)_y \right]. \end{aligned}$$

The Pauli matrices act on sublattice. A different Bloch-phase, valley, or layer-rotation convention moves signs and phases among the terms; spectra are unchanged when the basis is transformed consistently.

The leading interlayer tunneling has three harmonics,

$$T_\xi(\mathbf r) = \sum_{j=1}^{3} T_{\xi j} e^{i\xi\mathbf q_j\cdot\mathbf r},$$

where the three $\mathbf q_j$ connect nearby layer Dirac points and have magnitude $k_\theta$. A useful matrix convention is

$$\begin{aligned} T_{\xi j} ={}& w_0\sigma_0 \\ &+ w_1 \left[ \cos\varphi_j\,\sigma_x + \xi\sin\varphi_j\,\sigma_y \right], \\ \varphi_j &= \frac{2\pi(j-1)}{3}. \end{aligned}$$

$w_0$ and $w_1$ are effective continuum parameters, not bare hopping integrals at one atomic registry. The differences $\mathbf q_i-\mathbf q_j$ generate the moiré reciprocal lattice. Expanding a Bloch state in the momentum-space network connected by those vectors turns the continuum problem into a matrix eigenvalue problem; convergence must be checked against the plane-wave cutoff.

The TBG ledger. (a) Within one valley, three transfers $\mathbf q_j$ couple the rotated layer Dirac points. (b) Relaxation expands AB/BA regions and suppresses AA-like tunneling relative to AB/BA-like tunneling. (c) The two spin-degenerate, valley-degenerate active bands span $\nu=-4$ to $+4$ around charge neutrality; narrow bandwidth $W$ and remote gaps $\Delta_\pm$ are separate quantities. (d) Structural calibration, one-electron spectroscopy, thermodynamics, and phase-coherent transport constrain different links in a many-body claim.

The dimensionless coupling organizes the magic regime

For the minimal model, the principal dimensionless parameter is

$$\alpha = \frac{w_1} {\hbar v_F k_\theta}.$$

Reducing $\theta$ increases $\alpha$. Weak-coupling perturbation theory gives

$$\frac{v_*}{v_F} = 1-3\alpha^2+O(\alpha^4),$$

so interlayer paths interfere and reduce the Dirac velocity $v_*$. A full continuum calculation, rather than this truncated series, places the first velocity zero near $\alpha\simeq0.586$ in the simplest parameterization. In the chiral limit $w_0=0$, exactly flat bands occur at discrete couplings; realistic $w_0$, relaxation, strain, remote hopping, and interactions restore dispersion and particle-hole asymmetry.

This is why the magic angle is not a material constant. A quoted value near $1.05^\circ$ or $1.1^\circ$ silently assumes values of $v_F$, $w_0$, $w_1$, structural relaxation, and a definition of “magic.” Velocity suppression, minimum bandwidth, maximum remote-band isolation, strongest correlations, and optimal superconductivity need not occur at precisely the same angle.

Magic-Angle Flat Bands

Flatness is a scale hierarchy

A useful device report separates at least five scales:

Scale	Meaning	Why it matters
$W$	active-band bandwidth	residual kinetic dispersion
$\Delta_+$, $\Delta_-$	gaps to remote conduction and valence bands	validity of an active-band projection
$E_C$	screened Coulomb scale	interaction strength before form factors
$\Gamma$	disorder and inhomogeneous broadening	whether narrow features are resolved
$k_BT$	thermal scale	which gaps and ordered phases survive

A narrow $W$ does not imply that the active bands are isolated, nor that a local-orbital Hubbard model is controlled. The TBG wavefunctions carry layer, sublattice, valley, and momentum-space texture; their projected interaction matrix elements depend on form factors and on screening by gates and remote bands.

For an active-band projection, a representative interaction is

$$H_{\mathrm{int}} = \frac{1}{2A} \sum_{\mathbf q} V(\mathbf q) :\bar\rho(-\mathbf q)\bar\rho(\mathbf q):,$$

with projected density

$$\bar\rho(\mathbf q) = \sum_{\mathbf k,m,n,f} \Lambda_{mn}^{f}(\mathbf k,\mathbf q) c^\dagger_{\mathbf k+\mathbf q,mf} c_{\mathbf k,nf}.$$

Here $m,n$ label active bands, $f$ labels spin and valley flavors, and $\Lambda$ is a Bloch-wave overlap. Discarding $\Lambda$ discards the quantum geometry and much of the sublattice structure that distinguish TBG from a featureless narrow band.

Flat bands are measured through complementary windows

Scanning tunneling spectroscopy can resolve narrow peaks, local gaps, symmetry breaking, and spatial inhomogeneity, but the tunneling matrix element and tip electrostatics affect intensity. Photoemission measures occupied dispersion and replicas but averages over a finite region and may operate at temperatures different from transport. Compressibility accesses $d\mu/dn$, while Landau fans and quantum oscillations constrain degeneracies and Fermi-surface areas. None alone reconstructs the complete interacting spectral function.

Interactions also reshape the bands. Hartree potentials depend strongly on filling and can broaden, invert, or pin portions of the active manifold; exchange can split flavors and alter topology. Agreement between a noninteracting continuum band and one spectrum at one filling is therefore not a license to reuse that band unchanged across the phase diagram.

Correlated Insulating States

Integer filling is a coordinate, not a diagnosis

The 2018 observation of an insulating state near $\nu=-2$ established that a partially filled narrow manifold could become strongly resistive. Subsequent devices found insulating or semimetallic behavior at several integer fillings, often with Landau-fan resets, compressibility anomalies, or flavor-polarized states. The detailed sequence depends on twist inhomogeneity, strain, screening, substrate alignment, and magnetic or displacement fields.

Calling every such state a Mott insulator is too specific. Competing explanations include:

• interaction-driven flavor polarization that opens a Slater-like or exchange gap;
• intervalley-coherent, valley-polarized, spin-polarized, or sublattice-polarized order;
• a Chern insulator or other topological mass;
• charge order favored by nonlocal interactions;
• a one-electron gap enhanced by interactions;
• disorder localization or percolation through an inhomogeneous gap landscape.

The Mott Insulators page gives the general diagnostic standard. In TBG, the strongest case combines activated transport with a chemical-potential jump or incompressible interval, spectroscopy of a gap and spectral-weight transfer, a consistent broken-symmetry signature, and a model that includes the measured structure.

For density $n$ and chemical potential $\mu$, the inverse electronic compressibility is conventionally written

$$\kappa^{-1} = n^2\frac{\partial\mu}{\partial n}.$$

A jump in $\mu$ across a filling interval is thermodynamic evidence for a charge gap. Negative $\partial\mu/\partial n$ can occur from exchange without implying instability of the complete gated device, because geometric capacitance and long-range electrostatics contribute to the measured response.

Flavor cascades expose interaction-driven reconstruction

In an approximate description, spin rotations and valley charge conservation produce a near-fourfold flavor space. As filling changes, exchange can favor sequential flavor polarization. Compressibility “sawteeth,” chemical-potential resets, changes in Landau-fan degeneracy, and tunneling spectra that reorganize at integer fillings support such cascade physics.

The approximate symmetry is not exact. Lattice-scale intervalley scattering, strain, substrate alignment, remote-band mixing, Zeeman coupling, and phonons select among candidate orders. A flavor count inferred from a Landau fan is also not a direct image of an order parameter: magnetic breakdown, small pockets, and reconstruction can alter the fan.

Entropy supplies another diagnostic. Re-entrant transport and thermodynamic behavior consistent with a Pomeranchuk-like effect indicate that some correlated states carry substantial spin or valley entropy. This supports local-moment-like physics in parts of the phase diagram, but it does not by itself prove a particular lattice model or zero-temperature order.

Superconductivity

What the experiments establish

Superconducting domes were first reported adjacent to a correlated state in near-magic TBG. Pressure, electrostatic tuning, and subsequent devices showed that superconductivity can also appear when a nearby correlated insulator is weak or absent. The data therefore rule out the universal claim that superconductivity must arise by doping one fixed Mott parent.

Evidence should be graded rather than collapsed into “a resistance drop”:

Observation	Supports	Important alternatives or missing step
resistance decreases on cooling	enhanced conduction or pairing fluctuations	current redistribution, metallic percolation, contact effects
resistance reaches the experimental floor	a coherent low-resistance path	an inhomogeneous filament can short the device
critical current and field suppress the state	collective superconducting response	heating and weak-link networks need calibration
nonlinear $V$–$I$ scaling near a BKT transition	two-dimensional phase unbinding	finite-size and inhomogeneity can round the scaling
tunneling or Andreev gap closes consistently with transport	pairing gap tied to the transition	pseudogaps and junction modeling remain relevant
kinetic inductance yields finite superfluid stiffness	condensate phase rigidity	extraction requires circuit and geometry calibration

For a homogeneous two-dimensional phase with stiffness $\rho_s$ in energy units, the ideal Berezinskii–Kosterlitz–Thouless jump is

$$k_BT_{\mathrm{BKT}} = \frac{\pi}{2} \rho_s\!\left(T_{\mathrm{BKT}}^{-}\right).$$

Finite area, disorder, and a distribution of local transition temperatures smear this relation. It is a stringent consistency test, not a universal fitting curve.

Pairing mechanism and gap structure remain open

The small bandwidth permits several channels to coexist:

Candidate ingredient	Attractive feature	Unresolved issue
electronic fluctuations	naturally tied to flavor, nematic, or intervalley correlations	the dominant fluctuation and controlled coupling regime are disputed
graphene phonons	provide intervalley and intravalley attraction; replica bands show strong electron–boson coupling in some samples	observed coupling does not establish that phonons dominate pairing
screened Coulomb plus phonons	allows retardation and momentum structure	screening and remote-band contributions are device dependent
quantum geometry	can contribute to superfluid weight when dispersion is very small	it constrains phase stiffness, not by itself the microscopic pairing glue

Tunneling and transport have reported anisotropic or nodal behavior in some devices, while other measurements can be fit more simply. Broken rotational symmetry may originate in the normal state, the superconducting state, strain, or their coupling. Pauli-limit comparisons are likewise not decisive without orbital depairing, $g$-factor, spin-orbit, and inhomogeneity analyses.

Micrometre-scale photoemission has observed nearly equally spaced flat-band replicas in superconducting, hBN-unaligned near-magic devices, consistent with strong coupling to an optical phonon. The same study explicitly did not establish that this coupling is the principal pairing mechanism. More recent kinetic-inductance measurements found an anisotropic temperature dependence and a superfluid stiffness larger than a conventional-dispersion estimate, strengthening the case that quantum geometry matters. These are important constraints on a theory, not closure of the mechanism debate.

Topology and Valley Physics

Approximate symmetries protect the neutral Dirac points

Neglecting intervalley scattering gives an approximate valley $U(1)$ conservation law. Physical time reversal exchanges the two valleys. Within one valley, the antiunitary combination $C_{2z}\mathcal T$ leaves momentum fixed and, in a suitable spinless basis, squares to $+1$. Together with threefold rotation, it protects the active-band Dirac crossings at charge neutrality in the ideal continuum model.

The two Dirac cones within one valley have the same chirality. Consequently, a two-band tight-binding model cannot simultaneously use exponentially localized Wannier orbitals and preserve all of the relevant valley, spatial, and antiunitary symmetries in an onsite form. This Wannier obstruction is often described as fragile topology: adding suitable remote trivial bands can remove the obstruction, even though the isolated two-band subspace remains topologically nontrivial in its symmetry representation.

This statement is not “each unbroken TBG band has a nonzero Chern number.” With $C_{2z}\mathcal T$ intact, the valley bands contain protected Dirac crossings and do not form individually isolated Chern bands. The obstruction concerns the connected two-band subspace and its symmetry-compatible localization.

Breaking the protecting symmetry can reveal Chern bands

Alignment with hexagonal boron nitride can break $C_{2z}$ by making the two graphene sublattices inequivalent. The neutral Dirac points can then gap, and valley-resolved bands may acquire nonzero Chern numbers. The values and signs depend on mass conventions, relaxation, remote-band hybridization, and which active band remains isolated.

Interactions can select one valley or flavor and convert valley Berry curvature into a net orbital magnetization and Hall response. The strongest integer quantum anomalous Hall diagnosis requires

$$\sigma_{xy} = C\frac{e^2}{h}, \qquad \sigma_{xx}\rightarrow0,$$

at zero applied magnetic field, together with reproducible magnetization reversal and a well-defined insulating gap. Hysteresis or a large anomalous Hall signal alone establishes orbital ferromagnetism more directly than exact topological quantization.

Near $\nu=3$, hBN-aligned near-magic devices have shown first a large hysteretic anomalous Hall effect and later an intrinsic quantized anomalous Hall state. Magnetic-field studies have also exposed sequences of interaction-driven Chern insulators. These results establish that topology, valley polarization, and interactions can cooperate in TBG; they do not imply that every integer correlated state is topological.

Evidence Ledger

A mature claim should report enough information to connect the following layers:

1. structure: local twist, heterostrain, relaxation, hBN alignment, gates, and disorder;
2. one-electron model: $v_F$, $w_0$, $w_1$, remote terms, displacement fields, and numerical cutoff;
3. active manifold: $W$, $\Delta_\pm$, form factors, degeneracies, and filling calibration;
4. many-body state: candidate order parameter, competing states, and approximation method;
5. observables: transport, capacitance, tunneling, photoemission, magnetism, and phase stiffness;
6. reproducibility: spatial variation, cooldown dependence, contact geometry, and device-to-device trends.

A fit at layer 5 cannot uniquely infer layer 4 if layers 1–3 are not constrained. Conversely, disagreement between devices can be scientifically useful when it tracks a measured control such as strain, screening distance, pressure, or substrate alignment.

Common Mistakes

Treating 1.1° as an exact constant

The first magic regime depends on continuum parameters, reconstruction, and the observable being optimized. State the angle uncertainty and the criterion for “magic.”

Equating a narrow density range with a narrow band

Density measures states per area. Bandwidth is an energy. Relating them requires a calibrated cell area, degeneracy, and dispersion.

Calling every resistance maximum a correlated gap

Resistance is sensitive to scattering, contacts, domains, and percolation. Add thermodynamic or spectroscopic evidence and exclude one-electron gaps.

Calling every integer state a Mott state

TBG supports flavor polarization, intervalley coherence, charge order, and Chern masses as well as local-moment-like regimes. “Correlated insulator” is the safer umbrella unless the mechanism is resolved.

Inferring pairing glue from one correlation

A phonon replica, nearby magnetic order, or large quantum-geometric stiffness constrains superconductivity but does not alone identify the attractive kernel.

Confusing valley topology with a measured charge Hall response

Opposite valleys can carry opposite Berry curvature and cancel. A nonzero charge Hall response requires valley imbalance or another time-reversal-breaking mechanism.

Open Questions

Question	What is already constrained	What would materially advance it
What is the minimal quantitative Hamiltonian?	continuum models capture the gross flat-band structure	one parameter set predicting local structure, spectroscopy, thermodynamics, and transport across filling
Which orders occur at each integer filling?	flavor reconstruction and gaps are established in many devices	order-parameter-sensitive probes tied to local strain and topology
What pairs the electrons?	superconductivity, anisotropy in some samples, phonon coupling, and geometric stiffness are observed	phase-sensitive gap measurements plus controlled isotope, screening, and strain comparisons
How universal is the phase diagram?	device-to-device variation is substantial and partly systematic	shared metrology and multi-probe measurements on the same active region
How important are remote bands?	they screen interactions and can alter topology	controlled calculations benchmarked against wide-energy spectroscopy
Can fractional topological phases be stabilized reproducibly?	integer Chern states are established	thermodynamic gaps, fractional charge or statistics diagnostics, and edge consistency
Which nonequilibrium states are intrinsic?	current, microwave, and optical driving access new regimes	calibrated heating, relaxation, and spatially resolved dynamics

Annual review is warranted because several entries in this table remain active experimental frontiers.

Exercises

Exercise 1: density scale of a near-magic device

For graphene lattice constant $a=0.246\,\mathrm{nm}$ and twist angle $\theta=1.10^\circ$, use

$$L_M = \frac{a}{2\sin(\theta/2)}$$

to estimate $L_M$, $A_M=(\sqrt3/2)L_M^2$, the density for one carrier per cell, and the density at $\nu=4$.

Solution

With $\theta/2=0.55^\circ$,

$$\begin{aligned} L_M &\simeq \frac{0.246\,\mathrm{nm}} {2\sin(0.55^\circ)} \\ &\simeq 12.8\,\mathrm{nm}. \end{aligned}$$

Therefore

$$\begin{aligned} A_M &= \frac{\sqrt3}{2}L_M^2 \simeq 142\,\mathrm{nm}^2, \\ \frac{1}{A_M} &\simeq 7.0\times10^{11}\,\mathrm{cm}^{-2}, \\ n_{\nu=4} &= \frac{4}{A_M} \simeq 2.8\times10^{12}\,\mathrm{cm}^{-2}. \end{aligned}$$

These are ideal geometric values. A filling calibration should use the locally measured moiré area and electrostatic offsets.

Exercise 2: a model-dependent magic angle

Take $w_1=110\,\mathrm{meV}$, $v_F=0.90\times10^6\,\mathrm{m\,s^{-1}}$, $a=0.246\,\mathrm{nm}$, and first-magic coupling $\alpha_1=0.586$. Using $k_\theta\simeq K\theta$ at small angle and $\hbar v_F\simeq0.592\,\mathrm{eV\,nm}$, estimate the corresponding $\theta$.

Solution

First,

$$K = \frac{4\pi}{3a} \simeq 17.0\,\mathrm{nm}^{-1}.$$

Solving $\alpha_1=w_1/(\hbar v_FK\theta)$ gives

$$\begin{aligned} \theta &\simeq \frac{0.110\,\mathrm{eV}} {(0.586)(0.592\,\mathrm{eV\,nm}) (17.0\,\mathrm{nm}^{-1})} \\ &\simeq 0.0186\,\mathrm{rad} \simeq 1.07^\circ. \end{aligned}$$

Changing $v_F$ or $w_1$ changes this estimate. That sensitivity is precisely why “the magic angle” must be accompanied by a parameter convention.

Exercise 3: count the active states

Explain why the combined valence and conduction active manifold contains eight states per moiré cell, yet the conventional filling range is $-4\leq\nu\leq4$.

Solution

There are two active band labels, one valence-like and one conduction-like. Each has two spins and two valleys, so

$$2_{\mathrm{band}} \times 2_{\mathrm{spin}} \times 2_{\mathrm{valley}} = 8$$

states per cell. Charge neutrality places the valence-like active band filled and the conduction-like active band empty. Removing all four valence states reaches $\nu=-4$; adding all four conduction states reaches $\nu=+4$. The eight-state manifold is therefore traversed over an eight-carrier interval centered at $\nu=0$.

Exercise 4: audit an insulating-state claim

A device has a resistance peak at $\nu=2$ and an activation fit over one decade in temperature. What additional evidence is needed before identifying a Mott insulator?

Solution

At minimum, one should seek:

1. local twist, strain, and density calibration;
2. a chemical-potential jump or incompressible interval;
3. spectroscopy showing a many-body gap or spectral-weight transfer;
4. Landau-fan or other evidence for the active flavor degeneracy;
5. tests against a one-electron gap, flavor-polarized Slater state, Chern state, charge order, and disorder localization;
6. consistency under field, displacement field, screening, and device repetition.

An activation fit supports a gap controlling transport. It does not identify the microscopic origin of that gap.

Exercise 5: BKT consistency

Assuming the ideal jump relation, estimate the phase stiffness immediately below a transition at $T_{\mathrm{BKT}}=1.5\,\mathrm K$. Use $k_B=0.08617\,\mathrm{meV\,K^{-1}}$.

Solution

Rearranging the jump condition,

$$\begin{aligned} \rho_s(T_{\mathrm{BKT}}^-) &= \frac{2}{\pi}k_BT_{\mathrm{BKT}} \\ k_BT_{\mathrm{BKT}} &= 0.1293\,\mathrm{meV}, \\ \rho_s(T_{\mathrm{BKT}}^-) &= \frac{2}{\pi} (0.1293\,\mathrm{meV}) \\ &\simeq 0.082\,\mathrm{meV}. \end{aligned}$$

A measured stiffness far below this value would be inconsistent with a homogeneous ideal BKT transition at $1.5\,\mathrm K$. Finite-size and inhomogeneous broadening should be modeled before drawing a stronger conclusion.

Exercise 6: valley cancellation and anomalous Hall response

Suppose the isolated active band in valley $K$ has Chern number $+1$ and its time-reversed partner in $K'$ has Chern number $-1$. Compare the Hall conductance when both valleys are equally occupied with the case in which one spin-resolved $K$ flavor is filled and its $K'$ partner is empty.

Solution

Equal occupation gives

$$C_{\mathrm{tot}} = (+1)+(-1) = 0,$$

so the charge Hall conductance cancels even though each valley has nonzero Berry curvature. If interactions select one filled spin-resolved $K$ flavor and leave the corresponding $K'$ flavor empty, that contribution has

$$\sigma_{xy} = \frac{e^2}{h}.$$

Exact quantization additionally requires a bulk charge gap, negligible dissipative conduction, and no other occupied bands contributing a compensating Chern number.

Summary

• TBG is governed by a spatially varying twist, strain, and relaxation texture, not one angle alone.
• The continuum model couples rotated Dirac Hamiltonians through three moiré tunneling harmonics; its parameters and basis convention must be stated.
• The first magic regime is a parameter-dependent region of velocity and bandwidth suppression, not an exact universal angle.
• The active conduction and valence bands contain eight states per cell and span the conventional filling interval $-4\leq\nu\leq4$.
• Correlated gaps and flavor cascades are well established, but “Mott insulator” is not a universal diagnosis.
• Superconductivity is established in selected near-magic devices; its pairing mechanism and gap symmetry remain unresolved.
• Valley symmetry and same-chirality Dirac cones produce a Wannier obstruction; symmetry breaking and flavor polarization can yield Chern and quantum anomalous Hall states.
• Trustworthy interpretation joins structure, one-electron parameters, thermodynamics, spectroscopy, transport, and reproducibility.

Connections

• Moiré Superlattices derives the generic moiré geometry, mini-zone, miniband, filling, and interaction-scale framework.
• Flat Bands supplies the platform-independent language of bandwidth, isolation, compact localization, projector geometry, Chern bands, and interaction amplification.
• Correlated Insulators in Moiré Systems supplies the filling, compressibility, order, and evidence standards needed to classify TBG insulating states.
• Moiré Superconductivity compares TBG with other moiré platforms through pairing symmetry, superfluid stiffness, BKT coherence, field response, and Josephson evidence.
• Moiré Topology places TBG valley Chern bands, flavor polarization, QAH states, and fractional candidates in the cross-platform evidence hierarchy.
• Graphene supplies the monolayer Dirac basis, valley convention, substrate masses, and strain perturbations.
• van der Waals Heterostructures develops stack assembly, local alignment, pressure, gates, contacts, and tunneling evidence.
• Hubbard Physics in Materials explains active-space projection, interaction tensors, solver validity, and model-to-material validation.
• Mott Insulators distinguishes Mott localization from band, Slater, charge-ordered, and disorder-driven insulating behavior.
• Spectral Functions gives the Green-function language behind tunneling and photoemission lineshapes.
• BCS Theory supplies the conventional pairing baseline whose assumptions are tested in narrow multicomponent bands.
• Josephson Effect develops phase-coherent and current-phase probes of superconductivity.
• Berry Curvature and Chern Numbers in Band Theory provide the geometric and response framework for valley Chern bands.
• Time Reversal explains the antiunitary symmetry that exchanges graphene valleys.

Further Reading

1. J. M. B. Lopes dos Santos, N. M. R. Peres, and A. H. Castro Neto, “Graphene Bilayer with a Twist: Electronic Structure,” Physical Review Letters 99, 256802 (2007), doi:10.1103/PhysRevLett.99.256802.
2. R. Bistritzer and A. H. MacDonald, “Moiré Bands in Twisted Double-Layer Graphene,” Proceedings of the National Academy of Sciences 108, 12233–12237 (2011), doi:10.1073/pnas.1108174108.
3. N. N. T. Nam and M. Koshino, “Lattice Relaxation and Energy Band Modulation in Twisted Bilayer Graphene,” Physical Review B 96, 075311 (2017), doi:10.1103/PhysRevB.96.075311.
4. L. Balents, “General Continuum Model for Twisted Bilayer Graphene and Arbitrary Smooth Deformations,” SciPost Physics 7, 048 (2019), doi:10.21468/SciPostPhys.7.4.048.
5. G. Tarnopolsky, A. J. Kruchkov, and A. Vishwanath, “Origin of Magic Angles in Twisted Bilayer Graphene,” Physical Review Letters 122, 106405 (2019), doi:10.1103/PhysRevLett.122.106405.
6. H. Yoo et al., “Atomic and Electronic Reconstruction at the van der Waals Interface in Twisted Bilayer Graphene,” Nature Materials 18, 448–453 (2019), doi:10.1038/s41563-019-0346-z.
7. S. Carr, S. Fang, Z. Zhu, and E. Kaxiras, “Exact Continuum Model for Low-Energy Electronic States of Twisted Bilayer Graphene,” Physical Review Research 1, 013001 (2019), doi:10.1103/PhysRevResearch.1.013001.
8. Y. Cao et al., “Correlated Insulator Behaviour at Half-Filling in Magic-Angle Graphene Superlattices,” Nature 556, 80–84 (2018), doi:10.1038/nature26154.
9. Y. Cao et al., “Unconventional Superconductivity in Magic-Angle Graphene Superlattices,” Nature 556, 43–50 (2018), doi:10.1038/nature26160.
10. M. Yankowitz et al., “Tuning Superconductivity in Twisted Bilayer Graphene,” Science 363, 1059–1064 (2019), doi:10.1126/science.aav1910.
11. X. Lu et al., “Superconductors, Orbital Magnets and Correlated States in Magic-Angle Bilayer Graphene,” Nature 574, 653–657 (2019), doi:10.1038/s41586-019-1695-0.
12. A. Kerelsky et al., “Maximized Electron Interactions at the Magic Angle in Twisted Bilayer Graphene,” Nature 572, 95–100 (2019), doi:10.1038/s41586-019-1431-9.
13. Y. Xie et al., “Spectroscopic Signatures of Many-Body Correlations in Magic-Angle Twisted Bilayer Graphene,” Nature 572, 101–105 (2019), doi:10.1038/s41586-019-1422-x.
14. Y. Choi et al., “Electronic Correlations in Twisted Bilayer Graphene near the Magic Angle,” Nature Physics 15, 1174–1180 (2019), doi:10.1038/s41567-019-0606-5.
15. D. Wong et al., “Cascade of Electronic Transitions in Magic-Angle Twisted Bilayer Graphene,” Nature 582, 198–202 (2020), doi:10.1038/s41586-020-2339-0.
16. U. Zondiner et al., “Cascade of Phase Transitions and Dirac Revivals in Magic-Angle Graphene,” Nature 582, 203–208 (2020), doi:10.1038/s41586-020-2373-y.
17. A. Rozen et al., “Entropic Evidence for a Pomeranchuk Effect in Magic-Angle Graphene,” Nature 592, 214–219 (2021), doi:10.1038/s41586-021-03319-3.
18. Y. Saito et al., “Independent Superconductors and Correlated Insulators in Twisted Bilayer Graphene,” Nature Physics 16, 926–930 (2020), doi:10.1038/s41567-020-0928-3.
19. P. Stepanov et al., “Untying the Insulating and Superconducting Orders in Magic-Angle Graphene,” Nature 583, 375–378 (2020), doi:10.1038/s41586-020-2459-6.
20. M. Oh et al., “Evidence for Unconventional Superconductivity in Twisted Bilayer Graphene,” Nature 600, 240–245 (2021), doi:10.1038/s41586-021-04121-x.
21. C. Chen et al., “Strong Electron–Phonon Coupling in Magic-Angle Twisted Bilayer Graphene,” Nature 636, 342–347 (2024), doi:10.1038/s41586-024-08227-w.
22. M. Tanaka et al., “Superfluid Stiffness of Magic-Angle Twisted Bilayer Graphene,” Nature 638, 99–105 (2025), doi:10.1038/s41586-024-08494-7.
23. H. C. Po, L. Zou, A. Vishwanath, and T. Senthil, “Origin of Mott Insulating Behavior and Superconductivity in Twisted Bilayer Graphene,” Physical Review X 8, 031089 (2018), doi:10.1103/PhysRevX.8.031089.
24. L. Zou, H. C. Po, A. Vishwanath, and T. Senthil, “Band Structure of Twisted Bilayer Graphene: Emergent Symmetries, Commensurate Approximants, and Wannier Obstructions,” Physical Review B 98, 085435 (2018), doi:10.1103/PhysRevB.98.085435.
25. J. Kang and O. Vafek, “Symmetry, Maximally Localized Wannier States, and a Low-Energy Model for Twisted Bilayer Graphene Narrow Bands,” Physical Review X 8, 031088 (2018), doi:10.1103/PhysRevX.8.031088.
26. H. C. Po, L. Zou, T. Senthil, and A. Vishwanath, “Faithful Tight-Binding Models and Fragile Topology of Magic-Angle Bilayer Graphene,” Physical Review B 99, 195455 (2019), doi:10.1103/PhysRevB.99.195455.
27. J. Ahn, S. Park, and B.-J. Yang, “Failure of Nielsen–Ninomiya Theorem and Fragile Topology in Two-Dimensional Systems with Space-Time Inversion Symmetry,” Physical Review X 9, 021013 (2019), doi:10.1103/PhysRevX.9.021013.
28. A. L. Sharpe et al., “Emergent Ferromagnetism near Three-Quarters Filling in Twisted Bilayer Graphene,” Science 365, 605–608 (2019), doi:10.1126/science.aaw3780.
29. M. Serlin et al., “Intrinsic Quantized Anomalous Hall Effect in a Moiré Heterostructure,” Science 367, 900–903 (2020), doi:10.1126/science.aay5533.
30. K. P. Nuckolls et al., “Strongly Correlated Chern Insulators in Magic-Angle Twisted Bilayer Graphene,” Nature 588, 610–615 (2020), doi:10.1038/s41586-020-3028-8.
31. E. Y. Andrei and A. H. MacDonald, “Graphene Bilayers with a Twist,” Nature Materials 19, 1265–1275 (2020), doi:10.1038/s41563-020-00840-0.
32. P. Törmä, S. Peotta, and B. A. Bernevig, “Superconductivity, Superfluidity and Quantum Geometry in Twisted Multilayer Systems,” Nature Reviews Physics 4, 528–542 (2022), doi:10.1038/s42254-022-00466-y.