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

# Magic-Angle Twisted Trilayer Graphene

Searching arXiv for the cited MATTG literature to ground the article in current arXiv records.
{"query":"magic-angle twisted trilayer graphene in-plane critical magnetic fields 2104.14026 2209.02199 2103.12083", "max_results": 10}
{"query":"magic-angle twisted trilayer graphene observation coexisting dirac bands and moire flat bands 2209.02199 higher-order van hove singularity 2106.14911", "max_results": 10}
{"query":"twisted trilayer graphene tunable platform correlated electrons 2006.00399 multi-moire trilayer graphene 2310.12961 exact flat bands 2305.16385", "max_results": 10}
Magic-angle twisted trilayer graphene (MATTG) is the alternating-twist three-layer graphene system in which the top and bottom layers are aligned and the middle layer is twisted so that adjacent interfaces have relative twists \(+\theta\) and \(-\theta\), with \(\theta\) experimentally about \(1.6^\circ\) and the relevant magic angle reported near \(1.57^\circ\) [2209.02199, 2103.12083]. In its mirror-symmetric form, its low-energy electronic structure combines moiré flat bands with a highly dispersive Dirac band, and this coexistence underlies a set of correlated and field-tunable phenomena that include correlated states near integer filling, superconductivity near \(\nu=\pm2\), Pauli-limit-violating in-plane critical fields, re-entrant superconductivity, and displacement-field-controlled phase reconstruction [2209.02199, 2103.12083, 2104.14026]. More broadly, twisted trilayer graphene is a multi-moiré platform controlled by two twist angles, relative layer displacement, and electric field, so its flat-band, topological, and even quasiperiodic regimes depend sensitively on geometry and relaxation [2006.00399, 2310.12961].

## 1. Structural definition and low-energy architecture

In the experimentally relevant mirror-symmetric geometry, the top and bottom layers are related by reflection through the middle layer. This mirror symmetry allows the trilayer problem to separate into one sector that behaves in a way mathematically related to magic-angle twisted bilayer graphene and another that supports a monolayer-like Dirac cone, so the low-energy structure contains both MATBG-like flat moiré bands and a gapless Dirac cone at the same time [2209.02199, 2106.15573]. A standard mirror-basis representation is
\[
|+\rangle \sim \frac{|1\rangle+|3\rangle}{\sqrt2},\qquad |-\rangle \sim \frac{|1\rangle-|3\rangle}{\sqrt2},
\]
with the flat-band sector coupled to the middle layer and the other sector remaining monolayer-like as long as mirror symmetry is preserved [2209.02199].

This coexistence distinguishes MATTG from magic-angle twisted bilayer graphene. In MATBG, the two twisted layers produce narrow flat bands near charge neutrality, with no symmetry-protected monolayer-like Dirac cone surviving at low energy, whereas in MATTG the trilayer structure yields a flat-band manifold plus an additional dispersive Dirac sector [2209.02199]. A complementary continuum formulation writes the ideal mirror-symmetric trilayer Hamiltonian as
\[
\tilde{\mathcal H}_\nu^\theta = \tilde{\mathcal H}^{\theta,\nu}_{\rm TBG} \oplus \tilde{\mathcal H}^{\theta,\nu}_{\rm MLG},
\]
so that the effective bilayer-like block inherits enhanced tunneling and the monolayer-like block remains metallic [2106.15573]. This decomposition is central to essentially every later discussion of superconductivity, magnetism, topology, and field response.

The middle-layer-twisted trilayer setting is also part of a larger twisted-trilayer family with two independent twist angles \(\theta_{12}\) and \(\theta_{23}\). In that broader setting, two bilayer moiré patterns coexist, their interference produces a moiré-of-moiré structure, and the extra twist degree of freedom makes twisted trilayer graphene more tunable than twisted bilayer graphene as a platform for correlation-enhanced electronic structure [2006.00399]. The mirror-symmetric magic-angle device geometry is therefore a special but experimentally central slice through a larger multi-moiré phase space.

## 2. Continuum descriptions, symmetry decomposition, and magic-angle conditions

The mirror-symmetric trilayer continuum model used for MATTG may be written, for one valley, as
\[
\mathcal{H}_{K}=\begin{bmatrix} h_{\theta/2}(\bm{k}) & T(\bm{r}) &0\\ T^{\dagger}(\bm{r}) & h_{-\theta/2}(\bm{k}) &T^{\dagger}(\bm{r}) \\ 0 & T(\bm{r})  & h_{\theta/2}(\bm{k}) \end{bmatrix},
\]
with
\[
h_{\theta}(\bm{k}) =  e^{i\frac{\theta}{2} \sigma_z} \big(v_D\bm{k} \cdot \bm{\sigma}\big) e^{-i\frac{\theta}{2} \sigma_z},
\]
and moiré-modulated interlayer tunneling \(T(\mathbf r)\) [2104.14026]. Mirror symmetry allows decomposition into even- and odd-parity sectors; one sector maps to an effective bilayer problem with tunneling enhanced by \(\sqrt{2}\), while the other is a dispersive graphene-like Dirac band [2104.14026]. This is the origin of the trilayer magic angle being shifted relative to bilayers.

More generally, arbitrary twisted trilayers with two independent twist angles are not simply “two aligned TBGs.” For generic \((\theta_{12},\theta_{23})\), the coupled momentum basis is formally infinite-dimensional, there is no global moiré Brillouin zone in the continuum limit, and the most useful diagnostics are the density of states and the merging of low-energy van Hove singularities at charge neutrality [2006.00399]. In that general setting, magic-angle behavior is organized by a continuous curve in \((\theta_{12},\theta_{23})\) space rather than a single angle, with one perturbative condition given by
\[
\alpha_{12}^2+\alpha_{23}^2=\frac{1}{3},
\]
away from the equal-angle diagonal [2006.00399]. A later multi-moiré analysis recast this structure as a “magic line” in the \((\theta_{12},\theta_{23})\) plane, with experimentally relevant examples at \((1.56^\circ,-1.56^\circ)\), \((1.85^\circ,1.85^\circ)\), \((1.41^\circ,-1.88^\circ)\), and \((1.1^\circ,-5.5^\circ)\) [2310.12961].

For commensurate rational twist-angle ratios \(p/q\), the chiral continuum limit reveals an even sharper structure. In that limit, for the three special relative displacements
\[
d=0,\qquad d=\pm d_0,\qquad d_0=\frac{1}{3p|q|}(a_1-a_2),
\]
there are exactly flat bands at an infinite set of magic angles for all coprime \(p,q\), and the exact magic-angle condition is the vanishing of a \(3\times3\) Wronskian,
\[
W(\alpha)=\det(\psi_{K_1},\psi_{K_2},\psi_{K_3})=0
\]
[2305.16385]. This exact construction shows that magic-angle trilayer physics is not limited to the experimentally dominant mirror-symmetric case; it belongs to a wider commensurate family with flat-band topologies that include \(C=\pm1\) and \(C=2\) sectors [2305.16385].

## 3. Spectroscopic observations and moiré reconstruction

Direct momentum-space and real-space spectroscopy established the defining single-particle structure of MATTG near charge neutrality. Nano-ARPES and STM/STS measurements showed the simultaneous presence of a gapless Dirac cone at \(K_{1,3}\) and a flat-band maximum near the mini-Brillouin-zone center \(\gamma\), thereby directly verifying the coexistence of a dispersive Dirac band and moiré flat bands [2209.02199]. The measured moiré reciprocal vector was
\[
0.083\pm 0.01~\text{\AA}^{-1},
\]
the momentum resolution was \(\sim 0.01~\text{\AA}^{-1}\), and the Dirac velocity was
\[
v_D \approx (1 \pm 0.04)\times 10^6~\text{m/s},
\]
essentially the monolayer value [2209.02199]. The flat-band feature remained non-dispersive over a momentum width of about \(0.2~\text{\AA}^{-1}\), while STM showed that the low-energy flat-band spectral weight was concentrated mainly in AAA regions with spatial FWHM about \(5.3\) nm [2209.02199]. The same spectroscopy also found a double-peak structure near the flat bands with an energy splitting of \(25\)–\(55\) meV, whose microscopic origin was left open and attributed only tentatively to \(C_3\)-symmetry-breaking strain, lattice relaxation, and/or electron correlations [2209.02199].

Low-temperature STM further showed that real superconducting-relevant twisted trilayers do not remain perfectly uniform. Instead, they reconstruct into near-magic-angle mirror-symmetric domains separated by localized moiré defects called twistons [2109.12631]. Representative samples displayed a small moiré wavelength
\[
\lambda \sim 9~\mathrm{nm},
\]
a larger moiré-of-moiré scale
\[
\Lambda \sim 70~\mathrm{nm},
\]
and a mismatch
\[
\delta_\theta \sim 0.2^\circ\text{--}0.3^\circ,
\]
close to values reported in superconducting transport devices [2109.12631]. For \(\Lambda \gtrsim 30~\mathrm{nm}\), the average internal twist angle within each reconstructed domain saturated near \(1.5^\circ\), so the trilayer locally locked into near-magic-angle mirror-symmetric domains, while the mismatch was absorbed by solitons and twistons [2109.12631]. The typical magic-domain lateral size was \(\sim 50~\mathrm{nm}\), comparable to the superconducting coherence length inferred from transport, and the local twist-angle standard deviation in a uniform region could be as small as \(0.03^\circ\) [2109.12631].

The reconstructed landscape has direct electronic consequences. In uniform \(1.55^\circ\) domains, the valence- and conduction-band flat-band peaks at charge neutrality were separated by \(\sim 18~\mathrm{meV}\) with average FWHM \(\sim 23~\mathrm{meV}\), narrowing to \(\sim 15~\mathrm{meV}\) near \(\nu\sim \pm2\) [2109.12631]. In magic domains the spectra resembled unstrained near-magic TTG, whereas on solitons the flat-band spectral weight was strongly suppressed, and at twistons the flat bands re-emerged but were split by roughly \(\sim 40~\mathrm{meV}\), consistent with local structure closer to \(\sim 1.8^\circ\) [2109.12631]. This spatial partition implies that superconducting and correlated states in MATTG develop not on a perfectly homogeneous background, but on a reconstructed superstructure of magic domains, solitons, and twistons.

## 4. Correlated phases and superconductivity

Transport established that MATTG supports correlated states at integer fillings and superconductivity near \(\nu=\pm2\). In dual-gated Hall-bar devices, the zero-field resistance map showed correlated resistive states at \(\nu=+1,\pm2,+3,\pm4\), while superconductivity appeared near \(\nu=\pm2\), with the strongest dome on the hole-doped side near \(\nu=-2-\delta\) and the highest critical temperature approaching
\[
T_c\approx 2.9~\mathrm{K}
\]
at
\[
|D/\varepsilon_0|\sim 0.4\text{--}0.5~\mathrm{V/nm}
\]
[2103.12083]. The most robust superconductivity occurs near \(\nu\approx -2.4\), where the midpoint transition temperature drops from about \(2.7\) K at zero field to about \(1.35\) K at \(10\) T, indicating strong but incomplete suppression by in-plane field [2103.12083].

Theoretical descriptions of pairing in MATTG split into several distinct classes. One atomistic spin-fluctuation calculation argued that local Hubbard correlations and Hartree-renormalized moiré bands generate low-energy antiferromagnetic spin fluctuations between nearby ferromagnetic instabilities, producing a spin-singlet nematic \(d\)-wave superconducting dome between \(\nu=-2\) and \(\nu=-3\) with
\[
T_c \approx 2~\mathrm{K}
\]
and strong enhancement on the electron-doped side under perpendicular electric field [2104.10176]. A different continuum-based treatment found inter-valley superconducting instabilities generated by long-wavelength charge fluctuations dressed by Coulomb interaction and longitudinal acoustic phonons, yielding degenerate spin-singlet/valley-triplet and spin-triplet/valley-singlet channels with critical temperatures of up to a few Kelvin for realistic parameters [2106.15573]. The two frameworks agree that the magic-angle flat-band structure is sufficient to support Kelvin-scale superconductivity, but they assign the dominant pairing glue to different fluctuation sectors.

Electrostatic and environmental tuning sharpened the role of Coulomb repulsion. In a double-layer device where a nearby Bernal bilayer graphene served as a screening layer, superconductivity in MATTG strengthened when the adjacent layer was compressible, and the superconducting density width \(\Delta n_{SC}\) increased as screening improved [2108.03338]. The same work reported \(T_c\sim 1.5\) K at \(n_{BLG}=0\) and \(T_c\sim 1.6\) K at large \(|n_{BLG}|\), together with Pauli-limit violation and a thermodynamic gap at \(\nu=+2\) of \(10.3\) meV [2108.03338]. Because superconductivity became stronger when Coulomb repulsion was screened, that study argued that Coulomb repulsion competes against pairing and pointed toward a pairing mechanism compatible with electron-phonon coupling and a spin-triplet, valley-singlet order parameter [2108.03338].

The displacement field is itself an active tuning parameter, not merely a band-structure perturbation. A strong-coupling slave-particle theory near \(\nu=\pm2\) argued that increasing \(D\) primarily shifts the Dirac cone and causes self-doping into the TBG-like sector, thereby driving a semimetal-to-superconductor transition rather than simply increasing a Kondo-like hybridization [2606.20849]. In that account, the parent state at even filling contains a Mott-reconstructed TBG sector plus an extra Dirac cone, and field-induced self-doping activates superconductivity by transferring charge into the lower Hubbard band [2606.20849]. A separate electrostatic study of graphene-metal contacts predicted that contact-induced charge transfer and interfacial electric fields could drive MATTG through two superconducting domes as a function of work-function difference, with a maximum over \(2\) K [2111.04451]. These analyses collectively indicate that superconductivity in MATTG is inseparable from the coupled control of filling, mirror-symmetry breaking, and the relative position of flat and Dirac sectors.

## 5. In-plane fields, symmetry protection, and anisotropy

The in-plane magnetic-field response of MATTG is one of its defining anomalies. For a conventional weak-coupling spin-singlet superconductor, the Pauli limit is
\[
B_P=(1.86\,\mathrm{T/K})\,T_c,
\]
yet transport in MATTG found superconductivity surviving to in-plane fields in excess of \(10\) T, with Pauli-violation ratios of about \(2\)–\(3\) across the superconducting dome [2103.12083]. At \(\nu=-2.28\) and \(D/\varepsilon_0=-0.41\,\mathrm{V/nm}\), the \(10\%\) criterion gave
\[
T_c^{10\%}(0)=1.56\,\mathrm{K},\qquad B_{c\parallel}(0)=9.41\,\mathrm{T},
\]
so
\[
\mathrm{PVR}\approx 3.23,
\]
while re-entrant superconductivity was observed near
\[
\nu=-2.4,\qquad D/\varepsilon_0=-0.31\,\mathrm{V/nm},
\]
where a zero-resistance phase disappeared around \(B_\parallel\sim 8\) T and reappeared above \(8\) T [2103.12083]. The low-field and high-field superconducting regions were labeled SC-I and SC-II, respectively, and the high-field superconducting phase was confined to roughly
\[
-2.52 < \nu < -2.28
\]
at intermediate displacement fields [2103.12083].

A symmetry-based explanation traced the anomalously large in-plane critical field to the combined twofold rotation and horizontal mirror symmetry \(\mathcal C_2\mathcal M_h\) of the trilayer geometry [2104.14026]. In a valley-singlet superconductor, the relevant degeneracy is
\[
\epsilon_{K'}(-\mathbf k)=\epsilon_K(\mathbf k).
\]
In MATBG an in-plane field breaks time reversal and nothing remains to enforce this relation, but in MATTG the field-coupled Hamiltonian still preserves \(\mathcal T\mathcal M_h\) and \(\mathcal C_2\mathcal M_h\), so exact intervalley degeneracy survives as long as mirror symmetry is not otherwise broken [2104.14026]. In the clean symmetric model this means the orbital mechanism does not generate an upper critical field at zero gate field; \(B_c\) is formally infinite. A perpendicular gate electric field breaks \(\mathcal M_h\), destroys \(\mathcal C_2\mathcal M_h\), hybridizes the flat bands with the dispersive Dirac band, and restores orbital pair breaking, with a representative calculation giving
\[
B_c \sim 12~\mathrm{T}
\]
at
\[
D/\epsilon = 0.1~\mathrm{V/nm}
\]
[2104.14026].

The perturbation logic may be summarized as follows.

| Perturbation | Broken symmetries | Preserved symmetry relevant to intervalley pairing |
|---|---|---|
| In-plane \(\mathbf B_\parallel\) | \(\mathcal T\), \(\mathcal C_3\), \(\mathcal M_h\) | \(\mathcal T\mathcal M_h\), \(\mathcal C_2\mathcal M_h\) |
| Gate field | \(\mathcal M_h\), \(\mathcal T\mathcal M_h\), \(\mathcal C_2\mathcal M_h\) | none of these |
| Lateral shift | \(\mathcal C_3\), \(\mathcal M_h\), \(\mathcal T\mathcal M_h\) | \(\mathcal C_2\mathcal M_h\) |

This classification implies that a gate field directly removes the symmetry protection of large \(B_c\), whereas a lateral shift mainly changes the density of states and anisotropy but does not by itself destroy the protected intervalley degeneracy [2104.14026].

Independent transport measurements probed broken rotational symmetry directly. Angle-resolved transport in a sunflower-geometry MATTG device found that at \(20\) mK the in-plane response developed a clean twofold oscillation with
\[
\alpha = 77^\circ,\qquad \eta=\frac{\Delta R}{R_0}=0.84,
\]
while a full-tensor fit gave
\[
\alpha=78^\circ,\qquad \frac{\delta \sigma}{\bar{\sigma}}=0.74
\]
[2211.01352]. The anisotropy tracked the cascade of isospin transitions, was strongest on the large-Fermi-surface side of those transitions, and was strongly suppressed in a more detuned \(\theta=1.33^\circ\) sample where correlation effects were weaker [2211.01352]. Below \(5\) K, the anisotropy was further reorganized by a low-temperature \(PT\)-breaking order detected in second-harmonic nonreciprocal transport, indicating that rotational-symmetry breaking, isospin reconstruction, and \(PT\)-breaking are intertwined rather than independent in MATTG [2211.01352].

## 6. Magnetism, quantum criticality, and heavy-fermion behavior

At charge neutrality, interaction effects in mirror-symmetric magic-angle twisted trilayer graphene can drive magnetic order even though the single-particle spectrum contains a Dirac cone in addition to flat bands. An atomistic Hubbard calculation found that turning on electron-electron interactions results in a metallic-to-antiferromagnetic transition at
\[
U_c \approx 0.6\,t,\qquad t=|\gamma_0|=2.835~\mathrm{eV},
\]
which is smaller than the corresponding Hartree-Fock critical scale in monolayer graphene, bilayer graphene, and twisted bilayer graphene [2501.18491]. In that state, the order is strongest in the middle layer and concentrated in the AAA moiré centers because the flat-band wavefunctions carry about \(50\%\) of their weight on the middle layer and about \(25\%\) on each outer layer, while the Dirac-cone wavefunctions live solely on the top and bottom layers [2501.18491]. The antiferromagnetic order opens a Mott gap in the flat-band sector but leaves the Dirac cone ungapped, so pristine charge-neutral MATTG is not driven into a fully gapped insulator by this mechanism [2501.18491].

Displacement field and encapsulation strongly modify that magnetic instability. In the same atomistic study, a perpendicular electric field hybridized the Dirac cone with the flat bands and increased the critical interaction needed for antiferromagnetism, while hBN-induced sublattice symmetry breaking strongly suppressed the ordered state and could remove it entirely in the studied range when \(\Delta_{\rm hBN}=50\) meV [2501.18491]. A separate field-theoretic analysis near charge neutrality showed that interactions close to an Ising Gross-Neveu quantum critical point strongly renormalize the two coexisting Dirac velocities of MATTG: the fast cone slows down, the slow cone speeds up, and in the infrared the velocities become equal,
\[
v_+ \to v_- \quad \text{as } E\to 0,
\]
with emergent Lorentz symmetry and strongly non-monotonic crossover behavior controlled by nearby repulsive fixed points [2112.14752]. This implies that the multivelocity Dirac structure of MATTG is itself dynamically reshaped by critical correlations.

At finite filling, especially \(v=3\), MATTG can instead realize electrically tunable heavy-fermion behavior. Transport in two dual-graphite-gated devices with twist angles \(1.52^\circ\) and \(1.45^\circ\) found a continuous displacement-field-driven transition from an antiferromagnetic semimetal to a paramagnetic heavy-fermion metal [2507.12254]. In the strong-coupling regime at
\[
D=0.9~\mathrm{V\,nm^{-1}},
\]
the resistivity at \(v=3\) showed a high-temperature
\[
R_{xx}(T)\sim \log T
\]
behavior, a coherence maximum at
\[
T^*\approx 54~\mathrm{K},
\]
and a low-temperature Fermi-liquid form
\[
R_{xx}=R_0+AT^2
\]
with
\[
A=24.7~\Omega/\mathrm{K}^2
\]
[2507.12254]. By comparison, a reference point at \(v=5.5\) had
\[
A=0.02~\Omega/\mathrm{K}^2,\qquad m^*\sim 0.06\,m_0,
\]
so Kadowaki-Woods scaling implied an effective mass in the heavy-fermion regime of roughly
\[
m^*\sim 2\,m_0
\]
[2507.12254]. Hall data and quantum oscillations further showed Fermi-surface reconstruction from a low-\(D\) frequency
\[
F=1.1~\mathrm{T}
\]
to a high-\(D\) frequency
\[
F=11.5~\mathrm{T},
\]
with critical scales
\[
D_{c1}=0.31~\mathrm{V\,nm^{-1}},\qquad D_{c2}=0.43~\mathrm{V\,nm^{-1}}
\]
marking the onset of Dirac-point alignment and a Lifshitz-like reconstruction near the quantum critical region [2507.12254]. This establishes that the coexistence of localized flat-band electrons and itinerant Dirac electrons in MATTG is sufficient to realize a field-controlled Kondo-lattice-like phase diagram in two dimensions.

## 7. Topology, higher-order singularities, and broader tunability

MATTG also hosts several forms of nontrivial band topology and singular quantum geometry. In mirror-symmetric twisted trilayer graphene, a zero-energy higher-order van Hove singularity was predicted with density-of-states divergence
\[
\rho(\omega)\sim |\omega|^{-1/3},
\]
arising from the combined merging of van Hove singularities and Dirac cones at zero energy [2106.14911]. The critical conditions are
\[
v=0,\qquad \eta=0
\]
in the symmetry-constrained \(k\cdot p\) theory near the moiré \(K\) point, and for realistic corrugation \(r=0.8\) and \(\theta=1.59^\circ\) the necessary band motion occurs near
\[
D_c/\varepsilon_0 \simeq 0.37~\mathrm{V/nm}
\]
[2106.14911]. Varying a third parameter such as corrugation drives a topological Lifshitz transition at
\[
\xi=2\gamma,
\]
where the DOS divergence becomes
\[
\rho(\omega)\sim |\omega|^{-2/5}
\]
and the local semiclassical orbits change from open to closed [2106.14911]. This singularity structure is stronger than the higher-order singularities previously discussed in twisted bilayers.

Band topology is also strongly stacking- and gate-dependent in middle-layer-twisted trilayers. A systematic study of twist angle \(\theta\), interlayer potential difference \(\Delta\), and top-bottom layer stacking \(\tau\) found that AA outer-layer stacking gives the narrowest bands near \(\theta\simeq1.5^\circ\) but remains mostly metallic, whereas AB stacking supports a wider \(1^\circ\sim1.7^\circ\) narrow-band range and, at finite \(\Delta\), opens primary and secondary gaps that isolate topological low-energy bands [2104.01570]. In the AB-like regime the isolated bands carry finite valley Chern numbers:
\[
C=\pm1\ \text{near } \theta\sim1.5^\circ,\qquad C=\pm3\ \text{for smaller angles},\qquad C=\pm2\ \text{for larger angles},
\]
with the sign reversing under \(\Delta\to-\Delta\) [2104.01570]. The same study also reported pronounced anisotropic LDOS strip patterns when \(\tau\) is the saddle-point stacking vector between AB and BA, demonstrating that stacking alone can lower the effective rotational symmetry of the double-moiré electronic structure [2104.01570].

At charge neutrality, the topological structure may be subtle even when integrated Chern numbers vanish. A mean-field Hubbard analysis found that the full four-band flat manifold and the valence and conduction two-band submanifolds each have zero multiband Chern number,
\[
C_{1+2+3+4}=0,\qquad C_{1+2}=0,\qquad C_{3+4}=0,
\]
yet the multiband Berry curvature is strongly structured near \(K\), \(K'\), and \(\Gamma\), a pattern described as hidden quantum geometry [2501.18491]. Increasing electric field reshapes that Berry-curvature texture even though the integrated Chern numbers remain zero, providing a direct tuning knob for wave-function geometry rather than only for band energies [2501.18491].

Finally, external pressure provides an additional route to magic-angle engineering. In mirror-symmetric twisted trilayer graphene of the AAA-based geometry, a full tight-binding calculation with lattice relaxation found a zero-pressure magic angle
\[
\theta_0=1.35^\circ
\]
and showed that a larger-angle device with
\[
\theta=1.89^\circ
\]
can be driven into the flat-band regime at a critical pressure of about
\[
P_c\approx 4~\mathrm{GPa}
\]
[2108.02450]. The pressure-induced flat-band state retains equal top and bottom layer weights at the band edges, preserves mirror symmetry, and supports both a high-energy long-lived plasmon near \(0.15\) eV and a low-energy flat-band collective mode whose damping is sensitive to the detailed flat-band shape [2108.02450]. Pressure, displacement field, stacking, twist-angle ratio, and local reconstruction therefore all act as independent control parameters for the topology and many-body phenomenology of magic-angle trilayers.

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