---
title: Mixed-Stacked Pentalayer Graphene
url: https://www.emergentmind.com/topics/mixed-stacked-pentalayer-graphene
type: topic
---

# Mixed-Stacked Pentalayer Graphene

Searching arXiv for recent papers on mixed-stacked pentalayer graphene and related mixed-stacking graphene phenomena.
Mixed-stacked pentalayer graphene denotes five-layer graphene in which Bernal and rhombohedral stacking motifs coexist within a single crystal, so that the low-energy spectrum combines distinct chiral sectors rather than reducing to a single Bernal-like or rhombohedral-like hierarchy. In the experimentally studied non-centrosymmetric sequence ABCBC, the structure can be regarded as an ABC trilayer coupled to an AB bilayer, yielding a “3+2” chiral decomposition with coexisting cubic and parabolic bands that hybridize [2505.12478]. More generally, mixed-stacked five-layer graphene is important because broken inversion and mirror symmetries can permit out-of-plane polarization, intrinsic layer asymmetry, and transport responses unavailable in centrosymmetric pentalayers [2305.10896]. Recent work has established ABCBC as a platform for intrinsic layer polarization, multiple flatbands, Lifshitz transitions, unusual Landau-level degeneracy patterns, ultra-low-field quantum Hall phenomena, and, at higher fields, a cascade of even-denominator fractional quantum Hall states [2505.12478], [2507.20695].

## 1. Definition, stacking taxonomy, and symmetry

Mixed stacking in few-layer graphene means that a single flake contains both Bernal and rhombohedral sequences, so global inversion or mirror symmetry can be broken and an out-of-plane polarization is no longer symmetry-forbidden [2305.10896]. For pentalayer graphene, the general mixed-stacking problem can be classified either by explicit letter sequences such as ABABA, ABCAB, ABCBC, ABACA, and related permutations, or by Bernal-section decompositions separated by stacking faults [1211.5193]. Within this taxonomy, ABCBC corresponds to the decomposition \((1,3)\), while ABACA corresponds to \((2,2)\), and pure Bernal ABABA to \((4)\) [1211.5193].

Among the six pentalayer sequences discussed in the 2025 transport study, ABCBC belongs to the “none” symmetry group, with no inversion and no mirror symmetry [2505.12478]. This absence of inversion and mirror symmetry is central to its phenomenology: atomic sites in different layers experience different chemical environments, producing layer-dependent onsite energies even at zero displacement field \(D\), a built-in internal electric field, and an intrinsic band gap at charge neutrality [2505.12478]. By contrast, centrosymmetric multilayers such as ABABA and ABCBA obey \(R_{xx}(n,D)=R_{xx}(n,-D)\), do not exhibit an intrinsic \(D=0\) gap from built-in polarization, and tune gaps symmetrically under \(D\) [2505.12478].

A related but distinct symmetry classification appears in the weak-ferroelectric analysis of mixed-stacking pentalayers built as \(n{\rm ABA}m\) twins with \(n+m=2\) [2305.10896]. In that framework, asymmetric five-layer twins such as 2ABA0 and 0ABA2 are polar, while the symmetric 1ABA1 remains nonpolar because a \(z\to -z\) mirror or inversion symmetry survives [2305.10896]. This suggests that polarity in five-layer graphene is not unique to ABCBC, but rather a broader consequence of mixed stacking whenever inversion and mirror symmetry are simultaneously broken.

## 2. Electronic structure and effective descriptions

In ABCBC-stacked pentalayer graphene, the low-energy spectrum contains both an ABC-like cubic band and an AB-like parabolic band that hybridize [2505.12478]. The constituent effective Hamiltonians near \(K\) are given in two-band form. For AB bilayer graphene,
$$
H_{BLG} = -\frac{1}{2m^*}
\begin{pmatrix}
0 & (\pi^\dagger)^2 \\
\pi^2 & 0
\end{pmatrix}
+ \Delta_{BLG}\,\sigma_z + V_0^{(BLG)}\,\mathbb{I},
$$
where \(\pi = p_x + i p_y\), \(m^* = \frac{\gamma_1}{2 v_F^2}\), and \(v_F = \frac{\sqrt{3}\,a\,\gamma_0}{2\hbar}\) [2505.12478]. For ABC trilayer graphene,
$$
H_{TLG} = \alpha
\begin{pmatrix}
0 & (\pi^\dagger)^3 \\
\pi^3 & 0
\end{pmatrix}
+ \Delta_{TLG}\,\sigma_z + V_0^{(TLG)}\,\mathbb{I},
$$
with \(\alpha \approx \frac{v_F^3}{\gamma_1^2}\), modified by \(\gamma_3\), \(\gamma_4\), and \(\gamma_2\) [2505.12478].

The coupled ABCBC system is described by a \(4\times 4\) block Hamiltonian
$$
H_{ABCBC} =
\begin{pmatrix}
H_{TLG} & H_c \\
H_c^\dagger & H_{BLG}
\end{pmatrix},
$$
where \(H_c\) encodes interface hybridization across the 3|4 interface, mediated mainly by \(\gamma_1\) and \(\gamma_3\) [2505.12478]. To leading order in \(k\),
$$
H_c(\mathbf{k}) \simeq
\begin{pmatrix}
t_0 + t_2\,k^2 & \lambda_3\,\pi \\
\lambda_3\,\pi^\dagger & t_0' + t_2'\,k^2
\end{pmatrix}.
$$
In the absence of hybridization, the dispersions reduce to \(E_{parabolic}(k)=\pm \beta k^2\) and \(E_{cubic}(k)=\pm \alpha k^3\); hybridization repels crossings, produces avoided gaps, and generates multiple low-energy flatbands whose extrema underpin Lifshitz transitions and large Berry curvature [2505.12478].

The broader mixed-stacking theory developed earlier used a minimal \(\gamma_0\)–\(\gamma_1\) model and a decomposition into Bernal sections separated by faults [1211.5193]. In that scheme, low-energy eigenstates are mostly localized in each Bernal section, and the spectrum is approximated by a collection of spectra of independent sections [1211.5193]. That analytical construction classifies bands as linear, quadratic, or cubic according to the section content. For ABCBC \((1,3)\), the minimal model predicts no even sections, hence no Dirac cone, no \(F2\)-type quadratic band, and no cubic \(F3\)-\(F3\) pair; only a nearly flat \(F3\) boundary mode may remain around zero energy [1211.5193]. The 2025 ABCBC transport study instead finds a hybridized parabolic-plus-cubic low-energy manifold once realistic Slonczewski–Weiss–McClure couplings, hBN-induced onsite terms, and interface hybridization are included [2505.12478]. This contrast reflects the difference between the idealized \(\gamma_0\)–\(\gamma_1\) section model and the full SWMcC-plus-DFT description.

The modeling in the transport experiment uses \(a=2.46\) Å, \(d=3.35\) Å, an hBN-induced on-site term of 18 meV, and a field-induced interlayer potential step \(\phi=eEd\) referenced to the third graphene layer and shifted by 5 meV to match experiment [2505.12478]. Typical graphite values quoted are \(\gamma_0\approx 3.1\) eV, \(\gamma_1\approx 0.39\) eV, \(\gamma_3\approx 0.315\) eV, and \(\gamma_4\approx 0.044\) eV [2505.12478].

## 3. Intrinsic layer polarization and weak ferroelectricity

The defining consequence of non-centrosymmetry in ABCBC is intrinsic layer polarization. Because layer onsite energies differ even at \(D=0\), the trilayer and bilayer blocks experience opposite internal fields and opposite gaps [2505.12478]. Experimentally, an Arrhenius analysis of \(R_{xx}(n=0)\) yields a finite intrinsic gap \(\Delta\approx 0.65\) meV at \(D=0\) [2505.12478]. A polarization order parameter is introduced as
$$
P = \sum_{i=1}^5 s_i n_i,
$$
with alternating \(s_i\) chosen to capture the trilayer-versus-bilayer dipoles; nonzero \(P\) at \(D=0\) signals spontaneous layer polarization originating from the non-centrosymmetric stacking [2505.12478].

A closely related literature frames such behavior as elemental weak ferroelectricity in mixed-stacked few-layer graphene [2305.10896]. There the polarization is primarily electronic, arising from registry-dependent redistribution of \(\pi\)-electron charge density in the presence of symmetry-breaking interlayer couplings and an asymmetric twin boundary, with no ionic displacement [2305.10896]. The polarization is computed from layer densities as
$$
P_z = e d \sum_{L=1}^{N} L n_L
$$
with
$$
n_L = 4 \sum_{\beta,\alpha_L,p} \left[|\psi_{\beta,p}^{\alpha_L}|^2 \Theta(E_F-\varepsilon_\beta(p)) - \frac{1}{4}\right],
$$
and screening enters via self-consistent interlayer potential differences with effective out-of-plane dielectric permittivity \(\epsilon_z \approx 2.6\) [2305.10896]. For five-layer twins, the predicted polarization is “weak,” of order \(0.05\)–\(0.1\) e/\(\mu\)m, corresponding to \(\sim 1\times 10^{11}\) cm\(^{-2}\) transferred between outer surfaces [2305.10896].

The transport results on ABCBC do not report ferroelectric hysteresis, and the weak-ferroelectric study notes that no ferroelectric hysteresis is observed in the ABCBC experiment, likely due to single-domain samples [2505.12478]. Domain walls or opposite-stacking seeds such as ABABC inclusions are suggested as possible ingredients required for switchable ferroelectricity as in mixed tetralayers [2505.12478]. A plausible implication is that intrinsic layer polarization in ABCBC should be viewed as established, whereas truly switchable ferroelectric behavior in mixed-stacked pentalayers remains conditional on domain structure and sliding pathways rather than already demonstrated in the same form.

## 4. Displacement-field response, band alignment, and Fermi-surface reconstruction

The displacement field \(D\) acts on ABCBC through layer-dependent potentials satisfying \(U_{i+1}-U_i=\phi=eEd\), referenced to the middle layer [2505.12478]. Because the stack lacks inversion and mirror symmetry, external \(D\) adds to or subtracts from the built-in fields, producing a strongly asymmetric response of the gap and band alignment [2505.12478].

For \(D<0\), the gap first grows and then decreases, giving a nonmonotonic dependence. For \(D>0\), the gap rapidly closes [2505.12478]. This asymmetry is forbidden in centrosymmetric ABA and ABC trilayers, which exhibit \(R_{xx}\) symmetry under \(D\to -D\), and therefore serves as a hallmark of ABCBC [2505.12478]. The asymmetry is not merely spectroscopic: for \(D>0\), the AB bilayer conduction band overlaps the ABC trilayer valence band, producing a two-carrier regime with \(R_{xx}\propto B^2\) and \(\gtrsim 50{,}000\%\) magnetoresistance at 6 T, whereas for \(D<0\) the bands remain separated and the two-carrier signature is absent [2505.12478].

The same tuning of \(D\) and carrier density drives multiple Lifshitz transitions in the Fermi-surface topology. Landau-level maps at \(B=1\) T reveal several regions with distinct degeneracies: region II has degeneracy 4, region III has degeneracy 8, and region IV on the hole side at \(D>0\) has degeneracy 12; resistive ridges in region V mark Lifshitz transitions [2505.12478]. Band calculations show four topology changes as the chemical potential sweeps through the hybridized flatbands under \(\pm D\), consistent with the measured ridges and degeneracy changes [2505.12478]. Trigonal warping in both \(J=2\) and \(J=3\) sectors, together with their hybridization, is essential to the multi-pocket regimes [2505.12478].

More generally, trigonal warping in the ABC and AB blocks converts circular isoenergy contours into triangular or multi-pocket structures [2505.12478]. This provides a direct microscopic explanation for why LL degeneracies change between 4, 8, and 12 as the Fermi level crosses different topological regimes of the hybridized multiband system [2505.12478].

## 5. Berry curvature, Landau quantization, and quantum Hall structure

For a generic two-band Hamiltonian \(H(\mathbf{k})=d_0(\mathbf{k})\mathbb{I}+\mathbf{d}(\mathbf{k})\cdot\boldsymbol{\sigma}\), the Berry curvature of the lower band is
$$
\Omega_z(\mathbf{k}) = \frac{1}{2}\frac{\mathbf{d}\cdot\left(\partial_{k_x}\mathbf{d}\times\partial_{k_y}\mathbf{d}\right)}{|\mathbf{d}|^3}.
$$
In ABCBC, broken inversion and mirror symmetry plus inter-block hybridization produce sizable Berry curvature near band edges and avoided crossings, while trigonal warping concentrates it into valley-contrasting hot spots around \(K\) and \(K'\) [2505.12478]. This supports anomalous or valley Hall responses and facilitates weak-field Chern insulating behavior when valley degeneracy is lifted [2505.12478].

Landau quantization inherits the coexistence of \(J=2\) and \(J=3\) chiral sectors. The general chiral scaling is
$$
E_n^{(J)} \propto B^{J/2}\,\sqrt{n(n-1)\cdots(n-J+1)}.
$$
For AB bilayer-like states, \(E_n \propto B\sqrt{n(n-1)}\), with a twofold orbital \((n=0,1)\) zero-energy degeneracy per valley and spin. For ABC trilayer-like states, \(E_n \propto B^{3/2}\sqrt{n(n-1)(n-2)}\), with a threefold orbital \((n=0,1,2)\) zero-energy degeneracy per valley and spin, giving a 12-fold zero-energy Landau level after including spin and valley [2505.12478]. In ABCBC, hybridization of these sectors yields LL fans with degeneracy 4, 8, and 12 depending on density and \(D\) [2505.12478].

A particularly notable observation is the \(\nu=-6\) quantum Hall state at exceptionally low magnetic field. At \(D\approx 0\), the first developed Hall plateau is \(\nu=-6\) already by \(B\approx 26\) mT; \(R_{xy}\) reaches \(\sim 0.93\times (h/6e^2)\) and \(R_{xx}\) drops at \(n\approx -1\times 10^{10}\) cm\(^{-2}\) and \(T=14\) mK [2505.12478]. The state survives for \(D\) between about \(-0.06\) and \(+0.05\) V/nm and, for \(D\lesssim -0.2\) V/nm, remains the first LL on both electron and hole sides [2505.12478].

Two origin scenarios are considered compatible with the data. One is LL quantization from the ABC-like cubic band, whose 12-fold zero-energy LL can generate the earliest robust plateau, with partial hybridization or symmetry breaking reducing the first resolved sequence to \(\nu=\pm 6\) [2505.12478]. The other is a weak-field Chern insulator via spontaneous valley polarization, in which built-in internal fields lift valley degeneracy even at \(B\to 0\), and a small \(B\) selects a valley to yield finite Chern number and \(\nu=\pm 6\) [2505.12478]. The data are explicitly stated to be compatible with either scenario.

## 6. Correlated states, domain physics, and experimental realization

ABCBC domains were identified by scanning near-field optical microscopy on exfoliated pentalayer graphene; mixed domains appear in \(\sim 30\%\) of flakes and are smaller than pure Bernal or rhombohedral domains [2505.12478]. Selected regions were isolated by AFM cutting, encapsulated by hBN, and fabricated into dual-gated Hall bars with one-dimensional edge contacts [2505.12478]. Transport measurements were carried out at \(T=1.5\) K in a VTI system and down to 14 mK in a dilution refrigerator, with \(D=(D_b+D_t)/2\) and \(n=(D_b-D_t)/e\) controlled by top and back gates [2505.12478]. DFT calculations used VASP with PBE-GGA, \(a=2.46\) Å, and \(c=3.35\) Å, and a SWMcC model was fitted to ab initio bands near \(K\) [2505.12478].

A weak hBN–graphene moiré of period \(\sim 12.8\) nm is present in one device, but its effects are reported to be weak and negligible for the transport phenomena emphasized in the ABCBC study [2505.12478]. The later fractional quantum Hall work reinforces this point: one device shows Brown–Zak oscillations and Hofstadter features only near LL crossings, while a second device without moiré shows a near-identical high-field phase diagram, indicating that the half-filled phases are intrinsic to ABCBC-5LG [2507.20695].

At higher magnetic fields, mixed-stacked pentalayer graphene exhibits a cascade of even-denominator fractional quantum Hall states at \(\nu=-5/2\), \(-7/2\), \(-9/2\), \(-11/2\), and \(-13/2\), interwoven with conventional odd-denominator Jain states [2507.20695]. These states arise within a hybridized zeroth Landau level that inherits two zero-mode orbitals from AB and three from ABC, for up to 20-fold degeneracy before interactions and symmetry breaking [2507.20695]. Tuning \(D\) shifts the active half-filled LL between two intra-ZLL branches: a more AB-like branch hosting \(-5/2\) and \(-7/2\) near \(D\simeq 0.4\)–0.5 V/nm, and a more ABC-like branch hosting \(-9/2\), \(-11/2\), and \(-13/2\) near \(D\simeq 0.24\)–0.36 V/nm [2507.20695].

Exact-diagonalization including Coulomb interactions, LL mixing, and \(D\) finds a sixfold quasi-degenerate ground-state manifold at the relevant half fillings, described as the hallmark of Moore–Read physics on a torus [2507.20695]. Chiral-graviton spectroscopy distinguishes Pfaffian versus anti-Pfaffian tendencies, with \(-5/2\), \(-9/2\), and \(-13/2\) leaning anti-Pfaffian and \(-7/2\), \(-11/2\) leaning Pfaffian [2507.20695]. This suggests that mixed-stacked pentalayer graphene is not only a multiband quantum Hall system but also a tunable platform for displacement-field control of non-Abelian candidate phases.

Spatial inhomogeneity and mixed-domain networks provide another experimental dimension. STM on few-layer graphene on mica observed triangular networks of partial dislocations separating ABA and ABC stacked domains, with stacking-specific LDOS signatures and a pronounced peak at about \(+0.25\) eV above the Fermi level exclusively in ABA areas [1207.5427]. No pentalayer was measured explicitly in that work, but the ABA–ABC contrast and partial-dislocation network are described as generic to tri- and multilayer graphene and therefore expected to persist in five-layer systems including locally mixed ABA|ABC coexistence [1207.5427]. This suggests that mixed-stacked pentalayers may often need to be understood not only as ideal single-domain crystals but also as mesoscale networks of stacking domains, solitons, and domain-wall scattering channels.

## 7. Relation to other pentalayer stackings and open questions

The physical behavior of mixed-stacked pentalayer graphene depends sharply on stacking sequence. The most relevant contrasts stated in the source literature are summarized below.

| Stacking class | Symmetry/property | Low-energy consequence |
|---|---|---|
| ABABA, ABCBA, ABCAB, ABACA | Mirror or inversion symmetric | \(R_{xx}(n,D)=R_{xx}(n,-D)\); no intrinsic \(D=0\) gap from built-in polarization [2505.12478] |
| ABCBC | No inversion, no mirror | Intrinsic layer polarization, \(\Delta\approx 0.65\) meV at \(D=0\), asymmetric gap response, 4/8/12 LL degeneracies [2505.12478] |
| 2ABA0 / 0ABA2 vs 1ABA1 | Asymmetric twins polar; symmetric twin nonpolar | Weak electronic \(P_z\) for asymmetric twins; \(P_z=0\) for 1ABA1 [2305.10896] |

Pure rhombohedral pentalayer graphene is described as having a strongly chiral \(J=5\) band with ultra-flat dispersion and inversion symmetry that allows symmetric gap tuning under \(D\) [2505.12478]. Centrosymmetric mixed stacks such as ABABA and ABCBA retain symmetric \(D\)-dependence and lack the intrinsic \(D=0\) polarization gap of ABCBC [2505.12478]. ABCAC, although also non-centrosymmetric, would lose its gap rapidly for both \(\pm D\); the experimentally observed asymmetric persistence for \(D<0\) is therefore used to rule out ABCAC in the studied devices [2505.12478].

Several open questions are stated explicitly. Distinguishing between LL-origin and Chern-insulator-origin mechanisms for the ultra-low-field \(\nu=-6\) state requires complementary probes such as nonlocal transport, edge-state measurements, STM, or optical circular dichroism [2505.12478]. Further theory is needed to refine the effective coupling matrices \(H_c\) from SWMcC fits and quantify Berry-curvature distributions [2505.12478]. For weak ferroelectricity, the role of screening, substrate dielectric environment, doping, disorder, and domain walls remains central, and the available predictions emphasize strong sensitivity to these factors [2305.10896]. In the high-field regime, the precise spin and valley polarization of the active subspace, the detailed pairing channel of the even-denominator states, and the impact of LL mixing and intrinsic particle–hole symmetry breaking remain unresolved [2507.20695].

Taken together, the current literature presents mixed-stacked pentalayer graphene as a family of five-layer graphene crystals in which stacking faults, hybridized chiral sectors, and broken spatial symmetries reorganize the electronic structure far beyond the standard Bernal-versus-rhombohedral dichotomy. In the specific ABCBC realization, intrinsic layer polarization, asymmetric displacement-field response, multiband flatband transport, anomalously low-field integer quantum Hall behavior, and displacement-tunable even-denominator fractional quantum Hall states place the system at the intersection of band topology, spontaneous symmetry breaking, and correlated quantum Hall physics [2505.12478], [2507.20695].

Source: https://www.emergentmind.com/topics/mixed-stacked-pentalayer-graphene