TriLayer Systematisation Overview

• TriLayer Systematisation is the systematic organization of three-layer systems that emphasizes stacking, twist, layer polarization, and symmetry constraints to guide physical responses.
• It reveals distinct phenomena in graphene, moiré materials, and correlated systems by introducing a middle-layer degree of freedom that influences spectroscopy, transport, and collective order.
• The framework employs core variables like stacking order, layer asymmetry, twist geometry, and diagnostic observables to predict electronic, phonon, and topological properties across diverse materials.

“TriLayer Systematisation” is an Editor’s term for the comparative organization of three-layer systems by stacking sequence, twist geometry, layer polarization, interlayer coupling, and the symmetry constraints that govern vibrational, electronic, magnetic, and topological responses. Across the works surveyed here, trilayers are not treated as a trivial extension of bilayers: the third layer introduces additional symmetry operations, additional interlayer channels, and a distinct middle-layer degree of freedom that can dominate spectroscopy, transport, moiré reconstruction, or collective order. In graphene this appears as the ABA–ABC contrast in phonons, Raman activity, and Hall physics; in moiré materials it appears as two independent adjacent-layer moirés and, in some cases, a moiré-of-moiré or an intrinsically trilayer modulation; in other platforms it appears as layer-count-dependent superconductivity, multipole topology, skyrmion-bag robustness, or exact three-layer entropic competition (Zan et al., 2023, Dunbrack et al., 2023, Nakatsuji et al., 29 Apr 2025, Zhang et al., 29 Jan 2025).

1. Core organizing variables

The most stable systematization of trilayers begins with five variables that recur across otherwise unrelated platforms. The first is stacking order, which fixes point-group symmetry and therefore mode activity, valley degeneracy, and selection rules. The second is layer asymmetry, imposed by gating, displacement field, disorder localization, or pressure, which can break inversion or mirror symmetry and redistribute spectral weight or current among layers. The third is twist geometry, which determines whether the trilayer is described by overlapping bilayer moirés, by a moiré-of-moiré construction, or by an intrinsically trilayer zero-sum condition in reciprocal space. The fourth is the special role of the middle layer: in several systems it is the least gate-sensitive layer, the layer on which moiré potentials sum, or the layer whose in-plane relaxation is enhanced. The fifth is the diagnostic observable, which may be Raman activity, near-field infrared response, Hall plateau sequence, Chern-number stability diagrams, shift current, edge-state counting, or bulk Meissner fraction (Dunbrack et al., 2023, Nakatsuji et al., 29 Apr 2025, Ceferino et al., 2023, Davydov et al., 13 Apr 2025).

A recurring correction to reductionist intuition is that trilayers are not always reducible either to a monolayer-plus-bilayer decomposition or to a simple moiré-of-moiré picture. ABA-stacked trilayer graphene can indeed be described at low energy as a hybrid of one linear Dirac-like band and one parabolic band, but intrinsically multilayer moiré heterostructures are defined precisely so that no proper subset of layers is singular by itself, and twisted trilayer WSe$_2$ exhibits a middle-layer potential governed by the summation rule $V_2(r) \approx -V_1(r)-V_3(r)$, which has no bilayer analogue (Xu et al., 2012, Dunbrack et al., 2023, Nakatsuji et al., 29 Apr 2025).

2. Symmetry and spectroscopic classification in trilayer graphene

Untwisted trilayer graphene provides the clearest symmetry-based taxonomy. ABA trilayer has point-group symmetry $D_{3h}$, is non-centrosymmetric and semimetallic, and its high-frequency in-plane optical phonon modes transform as $2E' + E''$. The $E'$ mode is both Raman and infrared active, while $E''$ is Raman active. ABC trilayer has point-group symmetry $D_{3d}$, is centrosymmetric and semiconducting, and its in-plane optical phonon modes transform as $2E_g + E_u$. Here $E_g$ is Raman active, $E_u$ is infrared active but Raman inactive by symmetry, and inversion-symmetry breaking can make the antisymmetric $V_2(r) \approx -V_1(r)-V_3(r)$0 mode Raman visible (Zan et al., 2023).

These symmetry distinctions generate direct spectroscopic fingerprints. At $V_2(r) \approx -V_1(r)-V_3(r)$1 gate and $V_2(r) \approx -V_1(r)-V_3(r)$2 with $V_2(r) \approx -V_1(r)-V_3(r)$3 excitation, ABA shows a single G-mode peak near $V_2(r) \approx -V_1(r)-V_3(r)$4, whereas ABC shows an additional lower-wavenumber peak near $V_2(r) \approx -V_1(r)-V_3(r)$5 besides the higher-wavenumber G-like component. Far-field infrared absorption confirms that this lower-wavenumber Raman feature coincides with the infrared-active phonon near $V_2(r) \approx -V_1(r)-V_3(r)$6–$V_2(r) \approx -V_1(r)-V_3(r)$7, so it is assigned to the $V_2(r) \approx -V_1(r)-V_3(r)$8 phonon made Raman active by inversion-symmetry breaking due to environmental doping; the device’s Dirac point is shifted to $V_2(r) \approx -V_1(r)-V_3(r)$9, indicating substantial hole doping. Under gate tuning, the symmetric Raman G mode hardens with increasing hole density in both stackings, while in ABC the low-wavenumber $D_{3h}$0-derived component softens with increasing hole density. The same study links first-order Raman intensity to electron–phonon coupling through a third-order time-dependent perturbation-theory expression and interprets the strong Raman visibility of the infrared-active $D_{3h}$1 mode in ABC as direct evidence of enhanced electron–infrared phonon coupling. Near-field infrared nanoscopy then separates off-resonant and resonant regimes: at $D_{3h}$2 and $D_{3h}$3, ABC regions are darker than ABA, whereas at the phonon resonance near $D_{3h}$4 the ABC signal is strongly amplified by oscillator-strength transfer from interband transitions to the phonon resonance, corroborating much stronger electron–infrared phonon coupling in ABC than in ABA (Zan et al., 2023).

A second, complementary spectroscopic taxonomy comes from rigid-plane interlayer shear modes. Trilayer graphene has two doubly degenerate shear branches at $D_{3h}$5: a high-frequency branch $D_{3h}$6 with normalized eigenvector $D_{3h}$7, and a low-frequency branch $D_{3h}$8 with normalized eigenvector $D_{3h}$9. Their frequencies are essentially stacking-independent, but Raman visibility is symmetry-controlled. In ABA, $2E' + E''$0 has $2E' + E''$1 symmetry and is strong in the in-plane backscattering geometry, producing a peak at $2E' + E''$2, while $2E' + E''$3 has $2E' + E''$4 symmetry and is suppressed in the same geometry. In ABC, $2E' + E''$5 has $2E' + E''$6 symmetry and is Raman-forbidden, whereas $2E' + E''$7 has $2E' + E''$8 symmetry and is Raman active, producing a peak at $2E' + E''$9. By contrast, the low layer-breathing mode is observed at $E'$0 in both stackings with similar frequency and intensity. This makes the shear modes, rather than the breathing modes, the robust symmetry-based fingerprint of trilayer stacking order (Lui et al., 2014).

3. Electrostatic control, transport, and quantum Hall taxonomy

Transport systematization in trilayer graphene is anchored by how stacking controls band overlap, gap opening, Landau-level degeneracy, and plateau sequence. In the integer quantum Hall effect, the zero-energy Landau-level degeneracy scales as $E'$1, so trilayer graphene has a 12-fold-degenerate zero-energy level for $E'$2 and $E'$3. In gated trilayers, the observed Hall plateaus $E'$4, together with the absence of the $E'$5 plateau, identify ABC stacking. Self-consistent Hartree calculations with Slonczewski–Weiss–McClure parameters reproduce the ABC sequence and exclude ABA for the reported device. More generally, the rhombohedral trilayer lacks plateaus at $E'$6 and $E'$7, whereas ABA, which lacks inversion symmetry, yields valley splitting and plateau-to-plateau steps of $E'$8 when Zeeman splitting is neglected (Kumar et al., 2011).

Double-gated transport sharpens the distinction. The external electric field $E'$9 opens a tunable band gap in ABC trilayers, for which the two-terminal square resistance at charge neutrality increases with $E''$0. In ABA trilayers, the same field increases band overlap and the charge-neutrality resistance decreases with increasing $E''$1. Under magnetic field, ABC shows plateaus at $E''$2 and $E''$3, with $E''$4 appearing under larger interlayer asymmetry; ABA develops plateaus at $E''$5 with $E''$6. This classification combines inversion symmetry, realistic next-nearest-layer couplings, and the different rearrangement of zero-energy Landau levels in the two stackings (Jhang et al., 2011).

A useful correction to a common assumption is that transistor operation in trilayers does not require a bulk band gap. In the disorder-induced field-effect proposal, the bias convention is $E''$7, $E''$8, $E''$9, and only the bottom layer is disordered. The relevant quantity is the layer-resolved local electron distribution $D_{3d}$0, which shows that the sign of $D_{3d}$1 can move current-carrying states into or away from the disordered layer. This produces large conductance ratios in both ABA and ABC, even though ABA remains gapless under bias. For strong bottom-layer disorder $D_{3d}$2 and $D_{3d}$3, the conductance-ratio peak width is $D_{3d}$4 in ABA and $D_{3d}$5 in ABC, while the “on” conductance remains typically $D_{3d}$6–$D_{3d}$7. The paper therefore identifies a bias-sensitive mobility gap, rather than a bulk spectral gap, as sufficient for switching (Xu et al., 2012).

4. Twisted and moiré trilayers

In twisted trilayers, systematization begins in reciprocal space. For an $D_{3d}$8-layer stack with Fourier components at reciprocal vectors $D_{3d}$9, a moiré-scale modulation appears when the sum $2E_g + E_u$0 is small; a singular structure is defined by the exact zero-sum condition $2E_g + E_u$1. An intrinsically $2E_g + E_u$2-layer moiré pattern is one for which the zero-sum involves all layers and no proper subset is singular by itself. This framework generalizes bilayer commensurability, explains why some trilayer patterns cannot be reduced to a “moiré of moiré,” and leads to a reduced configuration space: after modding out overall translations and “nontrivial trivial transformations,” the intrinsic trilayer configuration space becomes a 2D torus of zero-mode phases (Dunbrack et al., 2023).

Twisted trilayer WSe$2E_g + E_u$3 makes the same logic concrete in a continuum model with two adjacent-layer moirés. The geometry is set by $2E_g + E_u$4, $2E_g + E_u$5, and $2E_g + E_u$6, with moiré-of-moiré reciprocal vectors $2E_g + E_u$7. Relaxation stabilizes different commensurate domains depending on the sign of $2E_g + E_u$8: helical trilayers ($2E_g + E_u$9) form $E_g$0 and $E_g$1 domains, whereas alternating trilayers ($E_g$2) form $E_g$3 domains. The central systematizing result is the middle-layer summation rule $E_g$4. In helical $E_g$5 domains this produces a Kagome lattice of $E_g$6 maxima and a characteristic Kagome flat band; in alternating $E_g$7 domains it produces triangular quantum wells on layer 2 that are twice as deep as in a bilayer at the same twist, supporting $E_g$8-, $E_g$9-, and, at smaller angles, higher bound states. A perpendicular electric field enters as $E_u$0 and tunes layer polarization and interlayer hybridization among triangular, honeycomb, and Kagome orbital manifolds (Nakatsuji et al., 29 Apr 2025).

Symmetrically twisted trilayer graphene adds a different trilayer-specific ingredient: enhanced middle-layer relaxation. With the middle layer twisted by $E_u$1 relative to two AA-aligned outer layers, the $E_u$2 symmetry implies no net out-of-plane corrugation of the middle layer, while the adhesion forces from the outer layers add constructively in-plane. As a result, the middle-layer in-plane relaxation is enhanced by a factor of 2 relative to the bilayer case, and the continuum model carries doubled pseudogauge and in-plane pseudoscalar potentials on the central layer. In the linear regime, the pseudomagnetic-field hotspots are of order $E_u$3–$E_u$4 and are weakly dependent on twist angle for $E_u$5; the trilayer magic-angle plateau is shifted to around $E_u$6–$E_u$7, compared with $E_u$8–$E_u$9 in twisted bilayer graphene (Ceferino et al., 2023).

5. Correlated, nonlinear, magnetic, and topological trilayer phases

The most elaborate trilayer phase diagrams arise when layer polarization and interlayer Coulomb coupling are both active. In twisted trilayer graphene near $V_2(r) \approx -V_1(r)-V_3(r)$00, each insulating quantum Hall domain can be labeled by a tuple $V_2(r) \approx -V_1(r)-V_3(r)$01, the layer-specific Chern numbers counting filled fourfold-degenerate Landau levels in the top, middle, and bottom layers. Intra-layer transitions appear as vertical boundaries in the $V_2(r) \approx -V_1(r)-V_3(r)$02–$V_2(r) \approx -V_1(r)-V_3(r)$03 plane, top–bottom interlayer transitions as horizontal boundaries, and adjacent-layer interlayer transitions as slanted boundaries. The resulting stability diagram contains not only rectangular domains but also triangular and trapezoidal domains unique to three coupled layers. At certain interlayer transitions, $V_2(r) \approx -V_1(r)-V_3(r)$04 is strongly suppressed and can reach $V_2(r) \approx -V_1(r)-V_3(r)$05, behavior interpreted as evidence for strong interlayer Coulomb coupling, drag, and possible exciton condensation when outer-layer Landau levels are partially filled (Davydov et al., 13 Apr 2025).

Alternating twisted trilayer graphene at total twist $V_2(r) \approx -V_1(r)-V_3(r)$06 provides a large-angle limit in which the three layers remain Dirac-like but strongly Coulomb-coupled. Dual gates tune the total filling factor $V_2(r) \approx -V_1(r)-V_3(r)$07 and displacement field $V_2(r) \approx -V_1(r)-V_3(r)$08, while the single-particle spectrum is organized by layer-dependent potentials $V_2(r) \approx -V_1(r)-V_3(r)$09, $V_2(r) \approx -V_1(r)-V_3(r)$10, and $V_2(r) \approx -V_1(r)-V_3(r)$11, with a small $V_2(r) \approx -V_1(r)-V_3(r)$12 accounting for electron–hole asymmetry. At $V_2(r) \approx -V_1(r)-V_3(r)$13, three low-resistance states appear as $V_2(r) \approx -V_1(r)-V_3(r)$14 is swept, with representative configurations $V_2(r) \approx -V_1(r)-V_3(r)$15, $V_2(r) \approx -V_1(r)-V_3(r)$16, and $V_2(r) \approx -V_1(r)-V_3(r)$17. Hartree–Fock analysis attributes these to spin-resolved helical edge modes analogous to those of quantum spin Hall insulators, giving reduced $V_2(r) \approx -V_1(r)-V_3(r)$18 and $V_2(r) \approx -V_1(r)-V_3(r)$19. At $V_2(r) \approx -V_1(r)-V_3(r)$20, a narrow $V_2(r) \approx -V_1(r)-V_3(r)$21 dip occurs when the middle and bottom layers are each half filled while the top layer remains inert at $V_2(r) \approx -V_1(r)-V_3(r)$22, consistent with an interlayer excitonic quantum Hall state (Kim et al., 13 Sep 2025).

ABC trilayer graphene also supports a displacement-field-tunable nonlinear optical taxonomy. With inversion symmetry broken by a perpendicular displacement field $V_2(r) \approx -V_1(r)-V_3(r)$23, the DC bulk photovoltaic effect arises from shift current rather than injection current, and the allowed tensor components satisfy the $V_2(r) \approx -V_1(r)-V_3(r)$24-symmetry relations listed in the study. A low-frequency peak in $V_2(r) \approx -V_1(r)-V_3(r)$25 near $V_2(r) \approx -V_1(r)-V_3(r)$26 is shared with AB bilayer graphene, but ABC has a distinct sign change in the $V_2(r) \approx -V_1(r)-V_3(r)$27–$V_2(r) \approx -V_1(r)-V_3(r)$28 window as $V_2(r) \approx -V_1(r)-V_3(r)$29 is tuned through band closings, reopenings, and inversions between bands $V_2(r) \approx -V_1(r)-V_3(r)$30–$V_2(r) \approx -V_1(r)-V_3(r)$31 and $V_2(r) \approx -V_1(r)-V_3(r)$32–$V_2(r) \approx -V_1(r)-V_3(r)$33 near $V_2(r) \approx -V_1(r)-V_3(r)$34. The sign reversal is traced to the quantum geometric reorganization of the Berry curvature, quantum metric, and Hermitian connection, especially in the $V_2(r) \approx -V_1(r)-V_3(r)$35 and $V_2(r) \approx -V_1(r)-V_3(r)$36 transitions (Postlewaite et al., 2024).

The same trilayer logic extends well beyond graphene. In twisted trilayer magnets, rotating only the middle layer creates four interlayer-coupling patch types—FM–FM, AFM–AFM, AFM–FM, and FM–AFM—leading to a sequence of ferromagnetic, noncollinear-domain, MD-I, and MD-II phases as twist angle decreases; the MD-I to MD-II transition is first order and doubles the domain count per moiré unit cell from $V_2(r) \approx -V_1(r)-V_3(r)$37 to $V_2(r) \approx -V_1(r)-V_3(r)$38 (Kim et al., 2023). In trilayer $V_2(r) \approx -V_1(r)-V_3(r)$39-TaS$V_2(r) \approx -V_1(r)-V_3(r)$40, an odd-layer, noncentrosymmetric superconductor with point group $V_2(r) \approx -V_1(r)-V_3(r)$41 and a middle-layer mirror plane, the even-mirror, time-reversal-invariant $V_2(r) \approx -V_1(r)-V_3(r)$42 state has $V_2(r) \approx -V_1(r)-V_3(r)$43 nodes at the $V_2(r) \approx -V_1(r)-V_3(r)$44-centered pockets, while hole doping or uniaxial pressure stabilizes fully gapped time-reversal-symmetry-breaking $V_2(r) \approx -V_1(r)-V_3(r)$45 states with total Chern numbers $V_2(r) \approx -V_1(r)-V_3(r)$46 and $V_2(r) \approx -V_1(r)-V_3(r)$47, respectively (Chen et al., 2019). In vertically stacked BBH heterostructures, trilayer Topo–Tri–Topo and Tri–Topo–Tri stacks undergo interlayer-coupling-driven transitions at $V_2(r) \approx -V_1(r)-V_3(r)$48, and the appropriate invariant is not the nested Wilson loop but the multipole chiral number $V_2(r) \approx -V_1(r)-V_3(r)$49, which counts the number of zero-energy corner-state sets: $V_2(r) \approx -V_1(r)-V_3(r)$50 in Topo–Tri–Topo and $V_2(r) \approx -V_1(r)-V_3(r)$51 in Tri–Topo–Tri (Ishida et al., 2024).

6. Diagnostics, mechanics, and the breadth of the trilayer program

A mature trilayer systematization requires not only phase classification but also reliable diagnostics. In graphene, the most direct spectroscopic route to stacking identification is Raman mapping at the $V_2(r) \approx -V_1(r)-V_3(r)$52 phonon wavenumber near $V_2(r) \approx -V_1(r)-V_3(r)$53–$V_2(r) \approx -V_1(r)-V_3(r)$54: ABC is bright in this channel when inversion symmetry is broken, whereas ABA is dark. The same stacking contrast is reproduced by near-field infrared nanoscopy but Raman offers simpler instrumentation and higher throughput. A second diagnostic route is the complementary shear-mode response, with ABA displaying the $V_2(r) \approx -V_1(r)-V_3(r)$55 mode and ABC the $V_2(r) \approx -V_1(r)-V_3(r)$56 mode in the standard in-plane geometry. A third is the Hall-plateau sequence, where $V_2(r) \approx -V_1(r)-V_3(r)$57 with the missing $V_2(r) \approx -V_1(r)-V_3(r)$58 identifies ABC in the IQHE regime (Zan et al., 2023, Lui et al., 2014, Kumar et al., 2011).

Mechanical loading reveals that trilayer classification also depends on how stress is transmitted between layers. In embedded ABA trilayer graphene under uniaxial tension, both G and 2D peaks redshift and G-peak splitting into $V_2(r) \approx -V_1(r)-V_3(r)$59 and $V_2(r) \approx -V_1(r)-V_3(r)$60 is observed for the first time in trilayer graphene, with a splitting onset at $V_2(r) \approx -V_1(r)-V_3(r)$61 strain in 3LG and $V_2(r) \approx -V_1(r)-V_3(r)$62 in the adjacent 2LG segment. Under compression, the 3LG G and 2D peaks blueshift up to a critical compressive strain of about $V_2(r) \approx -V_1(r)-V_3(r)$63, after which the Raman shifts relax toward their zero-strain values, indicating out-of-plane buckling. The critical strain is roughly one fourth of the value for embedded monolayer graphene, despite the greater bending rigidity of 3LG, because the limiting weakness is interlayer van der Waals shear rather than bending stiffness alone (Tsoukleri et al., 2016).

The term “TriLayer Systematisation” also encompasses materials in which the trilayer degree of freedom is the organizing axis for bulk ordered matter rather than only for band topology or moiré kinematics. Pressurized trilayer nickelate Pr$V_2(r) \approx -V_1(r)-V_3(r)$64Ni$V_2(r) \approx -V_1(r)-V_3(r)$65O$V_2(r) \approx -V_1(r)-V_3(r)$66 provides a clean example: under pressure it undergoes a monoclinic-to-tetragonal $V_2(r) \approx -V_1(r)-V_3(r)$67 transition at about $V_2(r) \approx -V_1(r)-V_3(r)$68, the Ni–O–Ni angle increases from $V_2(r) \approx -V_1(r)-V_3(r)$69 to $V_2(r) \approx -V_1(r)-V_3(r)$70, zero resistance is achieved above $V_2(r) \approx -V_1(r)-V_3(r)$71, the maximum onset temperature reaches $V_2(r) \approx -V_1(r)-V_3(r)$72 at $V_2(r) \approx -V_1(r)-V_3(r)$73, the superconducting volume fraction is $V_2(r) \approx -V_1(r)-V_3(r)$74 at $V_2(r) \approx -V_1(r)-V_3(r)$75 and $V_2(r) \approx -V_1(r)-V_3(r)$76 at $V_2(r) \approx -V_1(r)-V_3(r)$77, and a Ginzburg–Landau fit gives $V_2(r) \approx -V_1(r)-V_3(r)$78 with $V_2(r) \approx -V_1(r)-V_3(r)$79. In this setting the trilayer is a structural and electronic control variable in its own right, not a convenience of exfoliation (Zhang et al., 29 Jan 2025).

Finally, the conceptual reach of trilayer systematization extends into model systems and engineered fields. Plasmonic trilayer moiré superlattices formed by three twisted skyrmion lattices become strictly periodic when both pairwise twist angles are bilayer-commensurate, and trilayer skyrmion bags reproduce all bilayer bag sizes under mirror-symmetric twists while also allowing nested bags. For $V_2(r) \approx -V_1(r)-V_3(r)$80 bags at disturbance strength $V_2(r) \approx -V_1(r)-V_3(r)$81, the exact Huygens trilayer robustness is $V_2(r) \approx -V_1(r)-V_3(r)$82, compared with $V_2(r) \approx -V_1(r)-V_3(r)$83 in the bilayer, and the mirror-symmetric trilayer optimum occurs at $V_2(r) \approx -V_1(r)-V_3(r)$84 (Schwab et al., 2024). In the exactly solvable three-layer model of the smectic transition, the system is “smectic” when the central layer is empty and “nematic” when all three are equally occupied; the free-energy comparison shows that the smectic is preferred only when particle tips are sufficiently wide and the density is sufficiently high, making the central-layer occupancy itself the order-parameter surrogate (King, 2023).

Taken together, these results define trilayer systematization as a general research program rather than a single material-specific classification. Its central proposition is that once three layers are present, symmetry, layer selectivity, and interlayer competition acquire a qualitatively new structure: the middle layer can become the dominant carrier of moiré potential, the preferred sink or source of charge redistribution, the layer with enhanced relaxation, the host of localized modes or quantum wells, or the channel through which collective order is stabilized or suppressed.