---
title: Rhombohedral Hexalayer Graphene Overview
url: https://www.emergentmind.com/topics/rhombohedral-hexalayer-graphene
type: topic
---

# Rhombohedral Hexalayer Graphene Overview

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 [2401.17779] [2604.00113] [2508.15909].

## 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}\) [2512.24306] [2412.09985].

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^*\) [2401.17779]. This surface localization underlies both the strong sensitivity to perpendicular electric fields and the unusually large interaction effects reported in thin rhombohedral graphene systems [1911.04565].

## 2. Low-energy Hamiltonians and flat-band formation

Near a valley, the low-energy electronic structure is commonly written in a two-band form,
\[
H_{\rm eff}(\mathbf{k})=
\begin{pmatrix}
+\tfrac{U}{2} & (\hbar v_F \pi)^6/\gamma_1^5 \\
(\hbar v_F \pi^\dagger)^6/\gamma_1^5 & -\tfrac{U}{2}
\end{pmatrix},
\]
with \(v_F\approx10^6\) m/s, \(\gamma_1\approx0.39\)–\(0.4\) eV, and \(U\) the interlayer potential difference induced by a perpendicular displacement field. For \(U=0\), the dispersion is \(E(k)\propto k^6\); finite \(U\) opens a gap at \(k=0\) and can drive the band shape from monotonic to a very flat bottom and then to a Mexican-hat structure [2401.17779] [2604.00113] [2508.15909].

The gate-tuned flattening is particularly explicit in dual-gated devices. In the band calculations reported for hexalayer rhombohedral graphene, \(U=40\), \(90\), and \(140\) meV correspond to \(D/\varepsilon_0\approx0.65\)–\(0.85\) V/nm, and the flattening is strongest when \(U/\gamma_1\sim0.2\)–\(0.4\), producing a nearly dispersionless disk of radius \(k_{\rm flat}\) [2604.00113]. Independent theory has shown that flattening need not require a finite displacement field: for \(N=6\), self-consistent nonlinear electrostatics at \(n=0\) and \(D=0\) yields a reduced surface-band bandwidth \(W_0\approx40\) meV versus \(W_{U=0}\approx80\) meV, while hole doping toward the depletion point reduces the bandwidth to \(\approx6\) meV and enhances the low-energy density of states by \(3\)–\(5\times\) [2605.24080]. This suggests that hexalayer systems can access a low-field flat-band regime as well as the more familiar large-\(D\) 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 \(U_{\rm inv}\approx(40\text{–}45)\) meV per layer, corresponding to \(E_\perp\approx(40\text{–}45)\) mV/nm [2508.15909]. In earlier rhombohedral-graphite work specialized to the six-layer case, the gapped phase was described as carrying a valley Chern number \(C=\pm6\), and the spin-resolved quantum spin Hall order was described by \(C_s=6\) with \(\nu=1\) [1911.04565].

## 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 \(25\,\mu\)m electropolished Cu foil. Under the reported conditions, the method produced single-crystal monolayer islands \(\sim600\,\mu\)m across hosting centrally nucleated few-layer stacks up to \(9\) layers, with rhombohedral domains up to \(\sim50\) micrometers square [2006.06667]. For hexalayer material specifically, typical contiguous rhombohedral domain areas reach \(\approx20\)–\(50\,\mu{\rm m}^2\), and rhombohedral stripes are typically \(1\)–\(4\,\mu\)m wide and can exceed \(20\)–\(100\,\mu\)m in length [2006.06667].

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 \(AB\!\to\!BC\) shift, whereas displacement along zigzag preserves ABA registry. In \(40\%\) of aligned few-layer crystals, alternating ABA/ABC domains occur, while pure ABA dominates the remaining \(60\%\); no purely ABC crystals were observed under those conditions [2006.06667]. 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 [2006.06667].

Optical identification has proceeded along two complementary lines. In six-layer samples, Raman mapping under \(532\) nm excitation distinguishes ABA and ABC domains through the 2D-band full width at half maximum and the M-band profile: ABA-6L shows \({\rm FWHM}(2D)\approx50\ {\rm cm}^{-1}\), whereas ABC-6L shows \({\rm FWHM}(2D)\approx60\)–\(65\ {\rm cm}^{-1}\) with a pronounced low-frequency shoulder; the ABC-6L M-band resolves into four sub-peaks \(P_{21}\)–\(P_{24}\) at \(\approx1748\), \(1760\), \(1775\), and \(1788\ {\rm cm}^{-1}\) [2006.06667]. 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 \(N=6\), the dominant ERS peaks occur at \(4008\pm40\ {\rm cm}^{-1}\) and \(5571\pm40\ {\rm cm}^{-1}\); deviations greater than \(\sim100\ {\rm cm}^{-1}\) from these values unambiguously signal a stacking fault [2401.17779].

## 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\(_1\) and FS\(_2\), emerge and disappear through successive Lifshitz transitions; quantum oscillation frequencies trace a “leaf”-shaped region in the \(n\)–\(E_\perp\) map [2508.15909]. In Hartree–Fock calculations for \(N=6\), scanning \((n,D)\) produces a cascade of quarter metal, PIP\(_{11}\), half metal, PIP\(_{21}\), three-quarter metal, and full symmetric metal phases, with nematic regions inside the QM and PIP\(_{11}\) lobes [2508.00986]. 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 \(0.3\)–\(0.5\times10^{12}\ {\rm cm}^{-2}\) [2604.00113]. In the flat-band regime of hexalayer graphene, \(m^*\approx3\,m_e\), and at \(n\approx0.4\times10^{12}\ {\rm cm}^{-2}\) with \(\varepsilon\approx4\), the interaction parameter satisfies \(r_s>100\), well above the zero-field Wigner-crystallization threshold \(r_s\approx37\) [2604.00113]. Experimentally, the insulating Wigner region \(R2\) appears for \(n_e\approx(0.3\text{–}0.5)\times10^{12}\ {\rm cm}^{-2}\) and \(D/\varepsilon_0\approx0.78\pm0.05\) V/nm, where \(R_{xx}>1\) M\(\Omega\), the \(I\)–\(V\) curves exhibit a threshold \(V_t\sim0.5\)–\(1.5\) mV, and large hysteresis is observed [2604.00113].

Immediately adjacent to \(R2\), the \(R3\) regime is metallic in \(R_{xx}\) but hole-like in Hall response even though the net gating is electron doping. The extracted Hall density rises to \(n_h\approx0.06\times10^{12}\ {\rm cm}^{-2}\), which is \(\simeq15\%\) of the nominal electron density; two-carrier fits yield \(\mu_e\to0\) and \(\mu_h\sim1\)–\(2\times10^5\ {\rm cm}^2/{\rm Vs}\). 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 \(0.4\) T while disobeying the Středa relation [2604.00113]. 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\(_1\) at \(n\approx(-1.0\text{ to }-3.0)\times10^{12}\ {\rm cm}^{-2}\) and \(E_\perp\approx(-30\text{ to }-40)\) mV/nm, with onset \(T_c\approx170\) mK and \(B_{c\parallel}\approx6\) T, and SC\(_2\) near \(n\approx-0.5\times10^{11}\ {\rm cm}^{-2}\) and \(E_\perp\approx-36\) mV/nm, with \(T_c\approx80\) mK and \(B_{c\perp}\approx1\) mT [2508.15909]. At larger displacement fields, angle-resolved transport resolved a stripe-ordered state in which the hard transport axis is thermally activated with \(\Delta\approx4.5\) 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 \(D\approx0.9\)–\(1.1\) V nm\(^{-1}\) and \(n\approx0.6\)–\(1.2\times10^{12}\ {\rm cm}^{-2}\), with zero-resistance onset \(T_c\) in the range \(0.5\)–\(1.2\) K [2504.05129].

Magnetic-field-induced superconductivity has been reported in more than one form. In one dual-gated device, an electron-doped insulating phase at \(n\approx0.8\times10^{12}\ {\rm cm}^{-2}\) and \(D/\varepsilon_0\approx-0.36\) V nm\(^{-1}\) was driven into a superconducting state above \(\mu_0H^*_\parallel\approx0.21\) T, persisting up to \(\sim2\) T, while the insulating boundary showed a thermally activated gap \(\Delta_I\approx0.1\) meV [2512.24306]. In a separate study, Deng et al. reported that a superconducting region emerges only when a small in-plane field \(B_\parallel\) is applied, shifts toward larger \(|D|\) as \(B_\parallel\) increases, and remains robust up to \(14\) T; the maximum \(T_c\) rises from \(\sim110\) mK at \(0.4\) T to \(\sim260\) mK at \(10\) T, and quantum oscillations indicate that the parent normal state is a nematic Fermi-surface-reconstructed phase [2603.13498]. An additional report identified a superconducting dome stabilized by out-of-plane field \(1\lesssim B_\perp\lesssim4\) T, with \(T_c\approx0.5\) K and a neighboring sequence of integer and half-integer quantum Hall states sharing the same onset temperature and a Středa slope \(t=2.5\) [2507.22026]. 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 \(\Delta_{\rm CI}\approx23\) meV, a ferrovalley region near \(E_\perp\approx\pm6\) mV/nm, and a ferroelectric orbital-magnetic phase around \(E_\perp\approx-34\) mV/nm were reported; the latter shows electric-field-reversible magnetic hysteresis consistent with a \(\mathbf{P}\!\cdot\!\mathbf{M}\) multiferroic order parameter [2508.15909]. A separate study found a “bubble” phase near \(n\approx0\) and \(D/\varepsilon_0\approx\pm0.2\) V nm\(^{-1}\) with hysteretic anomalous Hall response and ferroelectric switching [2512.24306]. At moderate \(|D|\approx0.15\) V/nm and low \(|n|\lesssim3\times10^{11}\ {\rm cm}^{-2}\), 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 [2509.24672].

## 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 \(\nu=1\) with \(\sigma_{xy}=\pm e^2/h\), three insulating states at \(\nu=2\) were distinguished as spin-antiferromagnetic, spin-polarized, and valley-polarized insulators, and fractional insulating states appeared at \(\nu=1/3\) and \(2/3\) at zero magnetic field, with a stripe phase at \(\nu=1/2\) under \(B_\perp\) in a large-twist-angle device [2412.09985]. In a broader twist-angle study, the \(\nu=1\) moiré Chern insulator on the moiré-distant side survived only for \(0<\theta<1.10^\circ\), while a \(\nu=2/3\) fractional Chern insulator required still smaller \(\theta\); \(\theta=0.49^\circ\) showed a clear \(C=2/3\) state, \(\theta=0.92^\circ\) reverted to an integer \(C=1\) insulator, and by \(\theta=1.10^\circ\) no topologically non-trivial \(\nu=2/3\) state was resolved [2510.15309].

These moiré results also sharpen a theoretical controversy. Conventional Hartree–Fock predicted the opposite twist-angle dependence for the \(\nu=1\) phase boundary, whereas exact diagonalization and generalized random-phase approximation found that quantum fluctuations lower the energy of the \(C=1\) state much more strongly at small \(\theta\), thereby inverting the Hartree–Fock phase diagram [2510.15309]. In the non-moiré setting, the ideal-limit theory of rhombohedral graphene has proposed that short-range repulsion can stabilize a CP\(^5\) layer-pseudospin skyrmion lattice in the hexalayer case, reconstructing the zero-mode band into a \(C=1\) flat Chern band and producing acoustic skyrmion-lattice phonons together with gapped chiral shape modes [2511.07402]. 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.

Source: https://www.emergentmind.com/topics/rhombohedral-hexalayer-graphene