Critical Twist Angles (CAs) Overview

• Critical Twist Angles (CAs) are specific twist configurations in layered materials where changes in geometry, band topology, or transport occur, exemplified by magic-angle graphene and semiconducting bilayers.
• They are determined through conditions such as Dirac velocity renormalization, SU(2) gauge field quantization, and structural energy minima, which drive the emergence of ultra-flat bands and quantum transitions.
• CAs influence a range of applications including superconductivity, friction modulation, and optical metrology, with their effects traceable via Raman spectroscopy and tight-binding analyses.

Critical twist angles (CAs) are twist angles at which a twisted layered system undergoes a qualitative change in geometry, band topology, transport, stability, or ordering. In moiré graphene they often coincide with discrete magic angles where flat bands emerge, but the literature uses the same idea more broadly: to denote structural energy minima in finite flakes, geometric angles where a two-dimensional moiré collapses into a one-dimensional pattern, threshold angles for friction or Josephson transport, and many-body quantum critical points. Conversely, some semiconducting bilayers exhibit ultra-flat bands at large twist without a sharply defined special degree, showing that a CA is not a universal single number but a system-dependent concept (Ren et al., 2020, Tao et al., 2022, Zhu et al., 2020, An et al., 19 Jul 2025, Biedermann et al., 2024).

1. Canonical magic-angle constructions in hexagonal moiré systems

In twisted bilayer graphene (TBG), the canonical critical-angle problem is formulated in terms of the dimensionless coupling

$\alpha=\frac{w}{\hbar v_F k_\theta},$

with $k_\theta=2k_D\sin(\theta/2)$. In the symmetric Bistritzer–MacDonald description, the first magic angle is commonly around $\theta\approx 1.08^\circ$, while the chiral-limit WKB treatment shows that zero-energy flat bands occur at a discrete sequence of $\alpha_n$ determined semiclassically by

$\alpha_n \simeq \left(n+\tfrac{1}{2}\right)1.47.$

This makes the magic angles asymptotically equally spaced in inverse angle and interprets them as quantized values of an SU(2) gauge-field problem rather than as a generic property of any moiré superlattice (Eaton et al., 2022, Ren et al., 2020).

That basic condition can be renormalized rather than destroyed. In asymmetric twisted bilayer graphene, where the two layers have different Fermi velocities $v_1$ and $v_2$, the flat-band condition becomes

$$\overline\alpha \equiv \frac{w}{\hbar k_\theta \sqrt{v_1 v_2}} = \frac{1}{\sqrt{3}},$$

so the relevant velocity scale is the geometric mean $\sqrt{v_1 v_2}$. In twisted tetralayers the analogous condition is

$$\sqrt{\alpha_1^2+\alpha_4^2+z^2\alpha_{23}^2}=\frac{1}{\sqrt{3}},$$

which turns the “magic angle” into a tunable manifold in multidimensional twist-angle space rather than a unique value (Eaton et al., 2022).

Trilayer graphene generalizes the same idea further. In $k_\theta=2k_D\sin(\theta/2)$0-twisted trilayer graphene, with $k_\theta=2k_D\sin(\theta/2)$1, exactly flat bands appear at an infinite set of magic angles for every rational $k_\theta=2k_D\sin(\theta/2)$2 in the chiral limit, provided the displacement between the two moiré patterns is one of three inequivalent values $k_\theta=2k_D\sin(\theta/2)$3. The resulting “magic angle butterfly” organizes the leading magic angles over the rational axis $k_\theta=2k_D\sin(\theta/2)$4, while the associated flat bands admit a classification by rank and Chern number, including $k_\theta=2k_D\sin(\theta/2)$5 and $k_\theta=2k_D\sin(\theta/2)$6 cases (Popov et al., 2023). In multi-moiré trilayer graphene with two generally incommensurate moiré periods, the same flat-band condition becomes a magic line in $k_\theta=2k_D\sin(\theta/2)$7 space, with experimentally relevant examples including $k_\theta=2k_D\sin(\theta/2)$8 for alternating twists and $k_\theta=2k_D\sin(\theta/2)$9 for helical twists (Yang et al., 2023).

2. Flat bands without a sharply tuned critical degree

Not all flat-band mechanisms rely on fine cancellation of Dirac kinetic energy and interlayer tunneling. A distinct route is to use semiconducting bilayers whose band-edge energies depend strongly on local stacking. In bilayer $\theta\approx 1.08^\circ$0-In$\theta\approx 1.08^\circ$1Se$\theta\approx 1.08^\circ$2, the stacking-dependent valence-band-edge offset is

$\theta\approx 1.08^\circ$3

with AA regions acting as potential wells for holes and AB regions as barriers. First-principles calculations show a monotonic reduction of the highest-valence-band bandwidth from $\theta\approx 1.08^\circ$4 meV at $\theta\approx 1.08^\circ$5, to $\theta\approx 1.08^\circ$6 meV at $\theta\approx 1.08^\circ$7, to $\theta\approx 1.08^\circ$8 meV at $\theta\approx 1.08^\circ$9, and to $\alpha_n$0 meV at $\alpha_n$1. The paper explicitly states that “the appearance of ultra-flat bands here is not sensitive to the twist angle as in bilayer graphene,” and that flat bands can form at large twist angles “without specific degree” once the stacking-induced band-edge contrast is sufficiently large (Tao et al., 2022).

The same paper shows that this is not a ferroelectric-specific effect. In twisted bilayer InSe at $\alpha_n$2, the maximum VBM energy difference between AA and AB stackings is about $\alpha_n$3 eV, the topmost valence band has bandwidth $\alpha_n$4 meV, and the wavefunctions localize at AA regions. The mechanism therefore shifts the emphasis from a single magic angle to a barrier-height criterion tied to orbital character and stacking sensitivity (Tao et al., 2022).

Twisted homobilayer transition-metal dichalcogenides provide a different non-Dirac route. There, the continuum model can be mapped to a lowest-Landau-level problem with one flux quantum per moiré unit cell, and the operational magic-angle condition is

$\alpha_n$5

where $\alpha_n$6 is the first-shell Fourier coefficient of the effective periodic potential in the projected Landau-level description. The paper shows that topological flat bands occur only when the moiré potential is sufficiently strong compared with tunneling, specifically for $\alpha_n$7, and that for unstrained MoTe$\alpha_n$8 the magic angle lies around $\alpha_n$9–$\alpha_n \simeq \left(n+\tfrac{1}{2}\right)1.47.$0, shifting to $\alpha_n \simeq \left(n+\tfrac{1}{2}\right)1.47.$1 when structural relaxation is included (Morales-Durán et al., 2023).

A local tight-binding description of twisted TMD homobilayers reaches a related conclusion from the opposite direction. In a three-orbital model containing honeycomb and triangular orbitals, the topological sequence of the three highest valence bands changes from

$\alpha_n \simeq \left(n+\tfrac{1}{2}\right)1.47.$2

as twist angle increases, and the topmost band develops a flat-band “magic” region centered around $\alpha_n \simeq \left(n+\tfrac{1}{2}\right)1.47.$3 for MoTe$\alpha_n \simeq \left(n+\tfrac{1}{2}\right)1.47.$4. In the honeycomb-effective limit, the microscopic flat-band condition is

$\alpha_n \simeq \left(n+\tfrac{1}{2}\right)1.47.$5

or equivalently

$\alpha_n \simeq \left(n+\tfrac{1}{2}\right)1.47.$6

making the critical angle an orbital-interference condition rather than a Dirac-velocity zero (Crépel et al., 2024).

3. Geometric and structural critical angles

A separate class of CAs is defined structurally rather than electronically. In finite triangular or hexagonal flakes, dislocation theory identifies intrinsically preferred twist angles as local extrema of the interfacial energy, equivalently of the total dislocation-line area or the mean displacement field. The basic moiré wavelength is

$\alpha_n \simeq \left(n+\tfrac{1}{2}\right)1.47.$7

For triangular flakes with dimensionless size $\alpha_n \simeq \left(n+\tfrac{1}{2}\right)1.47.$8, the selection rules are

$\alpha_n \simeq \left(n+\tfrac{1}{2}\right)1.47.$9

while for hexagonal flakes they become

$v_1$0

These angles depend on flake size, shape, and center stacking, and the paper treats them as structural critical angles distinct from electronic magic angles (Zhu et al., 2020).

The same structural framework introduces “magic sizes” that make a target electronic angle also a structural energy minimum. For $v_1$1, the reported sequences include $v_1$2 for AA-centered triangles and $v_1$3 for AA-centered hexagons. This connects structural stability to the reproducibility of electronic flat-band experiments (Zhu et al., 2020).

Low-symmetry rectangular lattices produce a different geometric CA. In twisted rectangular bilayers, a discrete set of angles causes the ordinary two-dimensional moiré to collapse into a one-dimensional stripe pattern. The universal formula is

$v_1$4

where $v_1$5 becomes an exactly preserved common period of both layers. For twisted bilayer PdSe$v_1$6, the first critical angle is

$v_1$7

and at that angle DFT finds directionally localized flat bands, localized charge densities, and strong spin-orbit coupling along the dispersive direction (An et al., 19 Jul 2025).

This geometric perspective makes clear that a CA need not mark band flattening in a hexagonal Dirac system. It can instead mark an extremum of interfacial energy in a finite flake or a dimensional change in the moiré lattice itself (Zhu et al., 2020, An et al., 19 Jul 2025).

4. Nonequilibrium, transport, and many-body criticality

Critical angles can also be dynamical or interaction-driven. In AA-stacked twisted multilayer graphene, the static system has no equilibrium flat-band magic angles, but it does exhibit a sequence of minimum twist angles (MTAs) where the bandwidth of the central bands is locally minimized: $v_1$8 Under circularly polarized light, these MTAs evolve into non-equilibrium magic angles (NEMAs) where isolated topological flat quasienergy bands appear. The paper defines NEMAs by nonzero Chern number, positive gaps $v_1$9, and bandwidth thresholds such as $v_2$0 meV for the first range and $v_2$1 meV for a higher-order range. A displacement field further reduces the bandwidth, from about $v_2$2 meV without field to $v_2$3 meV at $v_2$4 meV and $v_2$5 meV at $v_2$6 meV (Li et al., 3 Mar 2025).

At charge neutrality in moiré bilayer graphene, the critical angle can be a many-body quantum critical point rather than a band-structure magic angle. Self-consistent Hartree–Fock on a realistic Bistritzer–MacDonald model with screened Coulomb interactions finds a continuous transition from a Kramers intervalley-coherent insulator at small twist to a symmetric Dirac semimetal above a critical value. For $v_2$7, the crossing-point analysis gives

$v_2$8

and the paper argues that the transition belongs to the relativistic Gross–Neveu–XY universality class (Biedermann et al., 2024).

In twisted three-dimensional superconductors, the critical angle can be defined by momentum-space overlap rather than flatness. A microscopic BdG treatment of a toy layered 3D superconductor finds

$v_2$9

the angle at which the projected Fermi surfaces of the two twisted superconductors cease to overlap. Below $$\overline\alpha \equiv \frac{w}{\hbar k_\theta \sqrt{v_1 v_2}} = \frac{1}{\sqrt{3}},$$0, the junction is in a high-transparency regime; above $$\overline\alpha \equiv \frac{w}{\hbar k_\theta \sqrt{v_1 v_2}} = \frac{1}{\sqrt{3}},$$1, it crosses into a low-transparency regime. The critical current remains finite even beyond $$\overline\alpha \equiv \frac{w}{\hbar k_\theta \sqrt{v_1 v_2}} = \frac{1}{\sqrt{3}},$$2, unlike in the twisted normal-metal case, and the phase drop localizes at the interface, making the twist itself act as a Josephson junction (Tani et al., 13 Aug 2025).

A non-electronic dynamical threshold appears even in confined fluids. For monolayer water inside twisted bilayer graphene, molecular dynamics shows that the friction coefficient varies with twist angle only for small angles and becomes effectively twist-independent for twists greater than $$\overline\alpha \equiv \frac{w}{\hbar k_\theta \sqrt{v_1 v_2}} = \frac{1}{\sqrt{3}},$$3. By contrast, the equilibrium number density, structure, and perpendicular dielectric constant remain essentially twist-agnostic over the studied range (Majumdar et al., 2022).

5. Metrology and spatial inhomogeneity of critical angles

Because many CA phenomena are sharply angle-dependent or rely on local angle homogeneity, metrology is part of the subject rather than an external technicality. In twisted bilayer graphene at intermediate angles, confocal Raman spectroscopy can spatially map the local twist angle between $$\overline\alpha \equiv \frac{w}{\hbar k_\theta \sqrt{v_1 v_2}} = \frac{1}{\sqrt{3}},$$4 and $$\overline\alpha \equiv \frac{w}{\hbar k_\theta \sqrt{v_1 v_2}} = \frac{1}{\sqrt{3}},$$5 using a green excitation laser. The angle is extracted directly from the moiré-superlattice-activated Raman scattering of the transverse acoustic phonon mode, and the width of the TA peak contains information on spatial twist-angle variations on length scales below the laser spot size of $$\overline\alpha \equiv \frac{w}{\hbar k_\theta \sqrt{v_1 v_2}} = \frac{1}{\sqrt{3}},$$6 nm (Schäpers et al., 2021).

A closely related need appears in cuprate twistronics. In artificially stacked BSCCO heterostructures, twist angle is treated as “a critical tunable parameter” because superconducting properties depend strongly on relative orientation, while the interface is fragile under conventional microscopy. Polarization-resolved Raman spectroscopy addresses this by using the anisotropic Bi/Sr $$\overline\alpha \equiv \frac{w}{\hbar k_\theta \sqrt{v_1 v_2}} = \frac{1}{\sqrt{3}},$$7 mode at $$\overline\alpha \equiv \frac{w}{\hbar k_\theta \sqrt{v_1 v_2}} = \frac{1}{\sqrt{3}},$$8 as an optical orientation marker. The reported heterostructure examples yield Raman-derived twists of $$\overline\alpha \equiv \frac{w}{\hbar k_\theta \sqrt{v_1 v_2}} = \frac{1}{\sqrt{3}},$$9 and $\sqrt{v_1 v_2}$0, and the method is fully non-invasive (Sardo et al., 20 Nov 2025).

These measurement works underscore a recurrent point in the CA literature: the relevant control variable is not only the nominal fabrication angle but also its local spatial distribution. This is especially important when theory predicts narrow windows, broad manifolds, or domain-dependent local configurations rather than a single ideal global angle (Schäpers et al., 2021, Sardo et al., 20 Nov 2025).

6. Conceptual synthesis

Several misconceptions are resolved once the literature is viewed comparatively. First, a critical twist angle is not synonymous with a magic angle in the TBG sense. In some systems it is a discrete flat-band angle set by $\sqrt{v_1 v_2}$1-quantization or Dirac-velocity renormalization (Ren et al., 2020, Eaton et al., 2022); in others it is a structural energy minimum of a finite flake (Zhu et al., 2020), a geometric angle where the moiré becomes one-dimensional (An et al., 19 Jul 2025), a transport crossover such as $\sqrt{v_1 v_2}$2 in twisted 3D superconductors (Tani et al., 13 Aug 2025), or a quantum critical point such as $\sqrt{v_1 v_2}$3 at charge neutrality in moiré bilayer graphene (Biedermann et al., 2024).

Second, there is no universal requirement that flat bands occur only near $\sqrt{v_1 v_2}$4. Flatness can be renormalized upward by asymmetry in Fermi velocities, producing larger magic angles in graphene multilayers (Eaton et al., 2022), or realized at large angles without a sharply defined special degree in semiconducting bilayers with strong stacking-dependent band edges (Tao et al., 2022). In TMD homobilayers and TTG, the relevant critical set can even be a line or manifold in multidimensional parameter space rather than an isolated scalar angle (Morales-Durán et al., 2023, Yang et al., 2023).

Third, angle sensitivity itself is system-dependent. TBG is sharply tuned because its magic angles rely on delicate cancellation, whereas $\sqrt{v_1 v_2}$5-In$\sqrt{v_1 v_2}$6Se$\sqrt{v_1 v_2}$7 supports ultra-flat bands over a broad angular window once $\sqrt{v_1 v_2}$8 is large (Tao et al., 2022). Structural selection rules in finite flakes accumulate toward $\sqrt{v_1 v_2}$9 with increasing flake size (Zhu et al., 2020), while rectangular lattices admit large geometric CAs such as $$\sqrt{\alpha_1^2+\alpha_4^2+z^2\alpha_{23}^2}=\frac{1}{\sqrt{3}},$$0 (An et al., 19 Jul 2025). This suggests that “criticality” in twist-angle space can arise from several non-equivalent mechanisms: interference in continuum Dirac models, stacking-controlled confinement, dislocation geometry, symmetry-protected moiré collapse, or interaction-driven order.

A comprehensive view of CAs therefore treats twist angle as a control coordinate in a larger design space that includes lattice symmetry, stacking dependence, interlayer tunneling, dielectric environment, external drive, electric field, flake geometry, and spatial homogeneity. The term remains precise only when tied to the specific instability, crossover, or geometric condition under discussion.