Rhombohedral Hexalayer Graphene Overview

• Rhombohedral hexalayer graphene is a six-layer allotrope with an ABCABC stacking sequence that localizes low-energy electronic states on its outer layers.
• Its effective low-energy Hamiltonian exhibits a k^6 dispersion, enabling electric and magnetic field tuning of insulating, metallic, and superconducting phases.
• Advanced synthesis and Raman mapping techniques reveal micron-scale domains and stacking variations critical for understanding its exotic correlated and topological behaviors.

Rhombohedral hexalayer graphene is a six-layer graphene allotrope with the stacking sequence $A\!\to\!B\!\to\!C\!\to\!A\!\to\!B\!\to\!C$ (ABCABC), in which each successive layer is displaced by one carbon–carbon bond relative to the layer below. In this stacking, the low-energy electronic states reside predominantly on the outermost layers, and the effective low-energy dispersion scales as $k^6$, producing extremely flat bands near charge neutrality. Because the band structure is highly tunable by carrier density, displacement field, and magnetic field, rhombohedral hexalayer graphene has become a platform for correlated insulators, Wigner crystallization, metallic Wigner crystal behavior, superconducting states, orbital magnetism, multiferroicity, and moiré Chern phases (Pálinkás et al., 2024, Han et al., 31 Mar 2026, Deng et al., 21 Aug 2025).

1. Crystal structure and symmetry

Rhombohedral hexalayer graphene consists of six graphene layers in the sequence ABCABC. Each layer is shifted by one carbon–carbon bond relative to the layer beneath; the in-plane lattice constant is $a_0\approx2.46$ Å and the interlayer spacing is $d\approx3.35$ Å. The structure has inversion symmetry and threefold rotational symmetry about the stacking axis; in the crystallographic description reported for rhombohedral hexalayer graphene, the stacking corresponds to space group $R\bar{3}m$, while the even-layer point-group description is $D_{3d}$ (Xie et al., 30 Dec 2025, Zheng et al., 2024).

A central structural consequence of ABC stacking is that the low-energy sector is surface dominated. In several continuum and tight-binding descriptions, the inner-layer sites are strongly hybridized, while the low-energy bands are localized on the top and bottom surface sublattices. For hexalayer rhombohedral graphene, more than $90\%$ of the flat-band amplitude resides on sublattices $A_1$ and $B_6$ within $|k|<k^*$ (Pálinkás et al., 2024). This surface localization underlies both the strong sensitivity to perpendicular electric fields and the unusually large interaction effects reported in thin rhombohedral graphene systems (Shi et al., 2019).

2. Low-energy Hamiltonians and flat-band formation

Near a valley, the low-energy electronic structure is commonly written in a two-band form,

$k^6$0

with $k^6$1 m/s, $k^6$2–$k^6$3 eV, and $k^6$4 the interlayer potential difference induced by a perpendicular displacement field. For $k^6$5, the dispersion is $k^6$6; finite $k^6$7 opens a gap at $k^6$8 and can drive the band shape from monotonic to a very flat bottom and then to a Mexican-hat structure (Pálinkás et al., 2024, Han et al., 31 Mar 2026, Deng et al., 21 Aug 2025).

The gate-tuned flattening is particularly explicit in dual-gated devices. In the band calculations reported for hexalayer rhombohedral graphene, $k^6$9, $a_0\approx2.46$0, and $a_0\approx2.46$1 meV correspond to $a_0\approx2.46$2–$a_0\approx2.46$3 V/nm, and the flattening is strongest when $a_0\approx2.46$4–$a_0\approx2.46$5, producing a nearly dispersionless disk of radius $a_0\approx2.46$6 (Han et al., 31 Mar 2026). Independent theory has shown that flattening need not require a finite displacement field: for $a_0\approx2.46$7, self-consistent nonlinear electrostatics at $a_0\approx2.46$8 and $a_0\approx2.46$9 yields a reduced surface-band bandwidth $d\approx3.35$0 meV versus $d\approx3.35$1 meV, while hole doping toward the depletion point reduces the bandwidth to $d\approx3.35$2 meV and enhances the low-energy density of states by $d\approx3.35$3–$d\approx3.35$4 (Kolář et al., 22 May 2026). This suggests that hexalayer systems can access a low-field flat-band regime as well as the more familiar large-$d\approx3.35$5 regime.

Berry curvature and topology enter naturally once inversion symmetry is broken. In the semimetallic regime studied in transport, the conduction and valence bands carry opposite Berry curvature in opposite valleys, and an electric-field-driven band inversion occurs at $d\approx3.35$6 meV per layer, corresponding to $d\approx3.35$7 mV/nm (Deng et al., 21 Aug 2025). In earlier rhombohedral-graphite work specialized to the six-layer case, the gapped phase was described as carrying a valley Chern number $d\approx3.35$8, and the spin-resolved quantum spin Hall order was described by $d\approx3.35$9 with $R\bar{3}m$0 (Shi et al., 2019).

3. Synthesis, stacking conversion, and optical identification

Large-area rhombohedral few-layer graphene was synthesized by low-pressure chemical vapor deposition on suspended Cu foils in a “back-side growth” geometry, using a 4″ cold-wall Aixtron BM-Pro reactor and a $R\bar{3}m$1m electropolished Cu foil. Under the reported conditions, the method produced single-crystal monolayer islands $R\bar{3}m$2m across hosting centrally nucleated few-layer stacks up to $R\bar{3}m$3 layers, with rhombohedral domains up to $R\bar{3}m$4 micrometers square (Bouhafs et al., 2020). For hexalayer material specifically, typical contiguous rhombohedral domain areas reach $R\bar{3}m$5–$R\bar{3}m$6, and rhombohedral stripes are typically $R\bar{3}m$7–$R\bar{3}m$8m wide and can exceed $R\bar{3}m$9–$D_{3d}$0m in length (Bouhafs et al., 2020).

The growth study showed that rhombohedral stacking was strongly correlated with Cu substrate morphology. Cu step bunching produces a wavy topography, and a directional shear-and-slip mechanism at step-bunch boundaries explains the formation of ABC regions: interlayer displacement along armchair induces an $D_{3d}$1 shift, whereas displacement along zigzag preserves ABA registry. In $D_{3d}$2 of aligned few-layer crystals, alternating ABA/ABC domains occur, while pure ABA dominates the remaining $D_{3d}$3; no purely ABC crystals were observed under those conditions (Bouhafs et al., 2020). A recurrent misconception is therefore that six-layer rhombohedral graphene is routinely obtained as a single, macroscopically pure ABC crystal by growth alone; the reported CVD route instead yielded stripe-like coexistence of Bernal and rhombohedral stacking within one crystal (Bouhafs et al., 2020).

Optical identification has proceeded along two complementary lines. In six-layer samples, Raman mapping under $D_{3d}$4 nm excitation distinguishes ABA and ABC domains through the 2D-band full width at half maximum and the M-band profile: ABA-6L shows $D_{3d}$5, whereas ABC-6L shows $D_{3d}$6–$D_{3d}$7 with a pronounced low-frequency shoulder; the ABC-6L M-band resolves into four sub-peaks $D_{3d}$8–$D_{3d}$9 at $90\%$0, $90\%$1, $90\%$2, and $90\%$3 (Bouhafs et al., 2020). For thicker rhombohedral specimens, however, the established 2D-peak method fails, and electronic Raman scattering provides a more stringent test of flawless stacking. For perfect $90\%$4, the dominant ERS peaks occur at $90\%$5 and $90\%$6; deviations greater than $90\%$7 from these values unambiguously signal a stacking fault (Pálinkás et al., 2024).

4. Electric-field tuning, fermiology, and electron crystallization

Carrier density and displacement field act as the primary control parameters in non-moiré devices. In the semimetallic regime, electron-like and hole-like pockets coexist near the band-inversion line, and two pockets, FS$90\%$8 and FS$90\%$9, emerge and disappear through successive Lifshitz transitions; quantum oscillation frequencies trace a “leaf”-shaped region in the $A_1$0–$A_1$1 map (Deng et al., 21 Aug 2025). In Hartree–Fock calculations for $A_1$2, scanning $A_1$3 produces a cascade of quarter metal, PIP$A_1$4, half metal, PIP$A_1$5, three-quarter metal, and full symmetric metal phases, with nematic regions inside the QM and PIP$A_1$6 lobes (Parra-Martínez et al., 1 Aug 2025). This suggests that the normal-state phase diagram is already highly structured before superconductivity or topological order are introduced.

Han et al. reported transport evidence for Wigner crystal and metallic Wigner crystal behavior in rhombohedral tetra-, penta-, and hexalayer graphene, including hexalayer devices, in the charge density range $A_1$7–$A_1$8 (Han et al., 31 Mar 2026). In the flat-band regime of hexalayer graphene, $A_1$9, and at $B_6$0 with $B_6$1, the interaction parameter satisfies $B_6$2, well above the zero-field Wigner-crystallization threshold $B_6$3 (Han et al., 31 Mar 2026). Experimentally, the insulating Wigner region $B_6$4 appears for $B_6$5 and $B_6$6 V/nm, where $B_6$7 M$B_6$8, the $B_6$9–$|k|<k^*$0 curves exhibit a threshold $|k|<k^*$1–$|k|<k^*$2 mV, and large hysteresis is observed (Han et al., 31 Mar 2026).

Immediately adjacent to $|k|<k^*$3, the $|k|<k^*$4 regime is metallic in $|k|<k^*$5 but hole-like in Hall response even though the net gating is electron doping. The extracted Hall density rises to $|k|<k^*$6, which is $|k|<k^*$7 of the nominal electron density; two-carrier fits yield $|k|<k^*$8 and $|k|<k^*$9–$k^6$00. The metallic Wigner crystal and the pinned Wigner crystal collapse simultaneously with increasing temperature or bias voltage, and the metallic Wigner crystal shows quantum Hall onset near $k^6$01 T while disobeying the Středa relation (Han et al., 31 Mar 2026). The coexistence of mobile holes with a transport-inert electron crystal is one of the distinctive correlated transport signatures reported so far in hexalayer rhombohedral graphene.

5. Superconductivity, orbital magnetism, and multiferroicity

Several superconducting regimes have been reported, but they occur in distinct parts of parameter space. In the semimetallic overlap region, two superconducting-like domes were reported: SC$k^6$02 at $k^6$03 and $k^6$04 mV/nm, with onset $k^6$05 mK and $k^6$06 T, and SC$k^6$07 near $k^6$08 and $k^6$09 mV/nm, with $k^6$10 mK and $k^6$11 mT (Deng et al., 21 Aug 2025). At larger displacement fields, angle-resolved transport resolved a stripe-ordered state in which the hard transport axis is thermally activated with $k^6$12 K, the easy axis remains metallic, and superconductivity appears only along the easy axis as one-dimensional-like channels. In that regime, three superconducting pockets were observed near $k^6$13–$k^6$14 V nm$k^6$15 and $k^6$16–$k^6$17, with zero-resistance onset $k^6$18 in the range $k^6$19–$k^6$20 K (Morissette et al., 7 Apr 2025).

Magnetic-field-induced superconductivity has been reported in more than one form. In one dual-gated device, an electron-doped insulating phase at $k^6$21 and $k^6$22 V nm$k^6$23 was driven into a superconducting state above $k^6$24 T, persisting up to $k^6$25 T, while the insulating boundary showed a thermally activated gap $k^6$26 meV (Xie et al., 30 Dec 2025). In a separate study, Deng et al. reported that a superconducting region emerges only when a small in-plane field $k^6$27 is applied, shifts toward larger $k^6$28 as $k^6$29 increases, and remains robust up to $k^6$30 T; the maximum $k^6$31 rises from $k^6$32 mK at $k^6$33 T to $k^6$34 mK at $k^6$35 T, and quantum oscillations indicate that the parent normal state is a nematic Fermi-surface-reconstructed phase (Deng et al., 13 Mar 2026). An additional report identified a superconducting dome stabilized by out-of-plane field $k^6$36 T, with $k^6$37 K and a neighboring sequence of integer and half-integer quantum Hall states sharing the same onset temperature and a Středa slope $k^6$38 (Nguyen et al., 29 Jul 2025). This suggests that the superconducting phenomenology is strongly regime-specific and may involve more than one pairing environment.

Orbital magnetic and multiferroic phases are equally prominent. Near charge neutrality, a correlated insulator with $k^6$39 meV, a ferrovalley region near $k^6$40 mV/nm, and a ferroelectric orbital-magnetic phase around $k^6$41 mV/nm were reported; the latter shows electric-field-reversible magnetic hysteresis consistent with a $k^6$42 multiferroic order parameter (Deng et al., 21 Aug 2025). A separate study found a “bubble” phase near $k^6$43 and $k^6$44 V nm$k^6$45 with hysteretic anomalous Hall response and ferroelectric switching (Xie et al., 30 Dec 2025). At moderate $k^6$46 V/nm and low $k^6$47, non-volatile and hysteretic anomalous Hall resistance can be toggled electrically, and small perpendicular fields produce a characteristic double sign reversal, indicating competition between distinct magnetic ground states (Krötzsch et al., 29 Sep 2025).

6. Moiré extensions, topology, and unresolved issues

When rhombohedral hexalayer graphene is aligned to hBN, the bottom graphene layer experiences a moiré potential that can further narrow and topologize the low-energy bands. In rhombohedral hexalayer graphene/hBN moiré superlattices, a switchable Chern insulator was observed at $k^6$48 with $k^6$49, three insulating states at $k^6$50 were distinguished as spin-antiferromagnetic, spin-polarized, and valley-polarized insulators, and fractional insulating states appeared at $k^6$51 and $k^6$52 at zero magnetic field, with a stripe phase at $k^6$53 under $k^6$54 in a large-twist-angle device (Zheng et al., 2024). In a broader twist-angle study, the $k^6$55 moiré Chern insulator on the moiré-distant side survived only for $k^6$56, while a $k^6$57 fractional Chern insulator required still smaller $k^6$58; $k^6$59 showed a clear $k^6$60 state, $k^6$61 reverted to an integer $k^6$62 insulator, and by $k^6$63 no topologically non-trivial $k^6$64 state was resolved (Huo et al., 17 Oct 2025).

These moiré results also sharpen a theoretical controversy. Conventional Hartree–Fock predicted the opposite twist-angle dependence for the $k^6$65 phase boundary, whereas exact diagonalization and generalized random-phase approximation found that quantum fluctuations lower the energy of the $k^6$66 state much more strongly at small $k^6$67, thereby inverting the Hartree–Fock phase diagram (Huo et al., 17 Oct 2025). In the non-moiré setting, the ideal-limit theory of rhombohedral graphene has proposed that short-range repulsion can stabilize a CP$k^6$68 layer-pseudospin skyrmion lattice in the hexalayer case, reconstructing the zero-mode band into a $k^6$69 flat Chern band and producing acoustic skyrmion-lattice phonons together with gapped chiral shape modes (Tan et al., 10 Nov 2025). The coexistence of experimentally observed Wigner crystals, switchable orbital magnetism, multiple superconducting regimes, and moiré Chern phases indicates that rhombohedral hexalayer graphene occupies a boundary region between flat-band metal, topological band system, and strongly correlated electron crystal, with the specific ground state selected by density, displacement field, magnetic-field orientation, substrate, and stacking perfection.