REVIEW 3 major objections 5 minor 1 cited by
Intrinsic layer polarization and multi-flatband transport in non-centrosymmetric mixed-stacked multilayer graphene
T0 review · 3 major / 5 minor · reviewed 2026-08-15 · deepseek-v4-flash
Pith's one-line read The paper reports an intrinsic band gap in ABCBC pentalayer graphene and a $\nu=-6$ quantum Hall plateau at ~26 mT.
desk verdict First ABCBC pentalayer transport study with a solid intrinsic-gap result, but the v=-6 quantum Hall claim is internally inconsistent as written: the quoted density and field give |nu|≈16, not 6. read the letter →
The pith
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
The reading
What carries the argument
The argument runs on the chiral decomposition of multilayer graphene: ABCBC splits into a chiral trilayer (ABC, cubic dispersion) and a chiral bilayer (AB, parabolic dispersion), which hybridize at low energy. The quantitative machinery is a Slonczewski–Weiss–McClure (SWMcC) tight-binding Hamiltonian with layer-resolved onsite potentials, fitted to density-functional calculations; the non-centrosymmetric lattice supplies built-in layer potentials that survive at zero external displacement field. This Hamiltonian produces the calculated gap-versus-$D$ curve, the Lifshitz-transition map, and the Landau-level degeneracies.
What would settle it
Re-measure a fresh ABCBC device with an independent density calibration (for example from the high-field integer quantum Hall fan), and check both the zero-field activation gap and the low-field plateau: if the gap at $D = 0$ vanishes, or if the calibrated filling of the plateau is not $\nu = -6$, the central claims fail. The quoted $n = -1\times10^{10}\ \mathrm{cm}^{-2}$ and $B = 26$ mT alone give $|\nu|\approx 16$, so the calibration check is decisive.
Extended reading notes
Core claim
The central claim is that the non-centrosymmetric stacking order itself, not any external field, creates an intrinsic layer polarization in ABCBC pentalayer graphene. Because the ABC trilayer block and the AB bilayer block sit in different chemical environments, they develop opposite built-in electric fields; after inter-block hybridization a small net gap of about 0.65 meV survives at $D = 0$. The paper shows this gap rises and then falls for negative displacement fields and closes for positive ones, and that the same two-band structure—a cubic band from the ABC part and a parabolic band from the AB part—produces multiple Lifshitz transitions and Landau-level degeneracies of 4, 8, and 12. It further reports a $\nu = -6$ quantum Hall plateau at roughly 26 mT, with $R_{xy}$ reaching about 93% of $h/6e^2$, which the authors attribute to the multi-flatband structure and spontaneous symmetry breaking; they note that a conventional 12-fold degenerate zero-energy Landau level from the ABC component and a valley-polarized Chern insulator mechanism are both possible explanations.
Load-bearing premise
The load-bearing premise is that the low-field plateau is at filling factor $\nu = -6$, which rests on the density calibration $n = (D_b - D_t)/e$ and the plateau identification; the paper's own numbers for that plateau ($n = -1\times10^{10}\ \mathrm{cm}^{-2}$, $B = 26$ mT, Fig. 4b) give $|\nu|\approx 16$, so the claim depends on an unstated correction or a corrected value.
Editorial extensions
If this is right
- If the central claim is correct, ABCBC pentalayer graphene is a twist-free platform where flat bands, broken inversion symmetry, and a tunable built-in gap coexist in one crystal.
- The measured Landau-level degeneracies of 4, 8, and 12, separated by resistive ridges, directly map the Lifshitz transitions of the combined cubic and parabolic bands as the Fermi energy moves through them.
- The low-field $\nu = -6$ plateau is presented as evidence that the multi-flatband structure and spontaneous symmetry breaking can produce a quantized Hall state at very small magnetic fields, without a moiré superlattice.
- Because single-domain ABCBC shows no ferroelectric hysteresis, the paper concludes that stacking-induced polarization alone does not produce memory behavior; an extra ingredient such as domain walls is needed.
Reading between the lines
- Beyond the paper, the same SNOM-plus-nanoARPES route could be applied to ABCAC pentalayer graphene; since the paper predicts a quickly vanishing gap on both field sides there, a measurement of its gap asymmetry would directly test the chiral-decomposition picture.
- A testable extension the paper does not report is a search for correlated insulating states at partial fillings of the flat bands; if the flat bands are as active as the low-field $\nu = -6$ state suggests, compressibility or capacitance measurements should reveal interaction-driven gaps.
- We also infer that the intrinsic dipole of a single ABCBC domain should become switchable if a reversed-stacking seed or domain wall is introduced; the paper reports no ferroelectric hysteresis, so a patterned ABCBC/ABABC junction would be the natural experiment to look for stacking-order memory.
Signed reviews
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The manuscript reports a transport and NanoARPES study of ABCBC-stacked pentalayer graphene, a non-centrosymmetric mixed stacking sequence. The authors claim three central results: (i) an intrinsic band gap of about 0.65 meV at charge neutrality and zero displacement field, which responds asymmetrically and non-monotonically to an applied vertical displacement field; (ii) multiple Lifshitz transitions and Landau levels with degeneracies 4, 8, and 12 originating from coexisting cubic and parabolic bands; and (iii) a robust quantum Hall plateau at filling factor ν = -6 developing at an exceptionally low magnetic field of about 20–26 mT. The identification of ABCBC stacking is supported by SNOM imaging and by NanoARPES dispersion fitting, and four devices show consistent transport behavior. The theoretical interpretation is based on a Slonczewski–Weiss–McClure model with parameters fitted to DFT, including a 5 meV shift in the interlayer potential difference to match experiment.
Significance. If the results hold, the paper establishes mixed-stacked multilayer graphene as a platform combining broken inversion symmetry, multiple flat bands, and tunable band topology in a natural (non-moiré) crystal. The observation of an intrinsic gap at D = 0 that is reproducible across four devices and consistent with NanoARPES is a substantive advance. The explicit comparison with ABCAC stacking and the use of measured activation gaps are strengths, as is the direct imaging-based stacking identification. However, the strength of the central low-field ν = -6 claim is weakened by an internal quantitative inconsistency in the reported density versus field, and the 'quantitative agreement' of the gap calculations relies on a fitted 5 meV offset that is not flagged in the main text. These issues are fixable, but they must be resolved before the claims can be accepted as stated.
major comments (3)
- [Low-field |v| = 6 quantum Hall state, Fig. 4b] The reported density and field for the ν = -6 plateau are internally inconsistent. The text states that Rxy approaches h/6e^2 at n = -1×10^10 cm^-2 and B = 26 mT. Using ν = n h/(e B), with h/e = 4.135667×10^-11 T cm^2, n = -1×10^10 cm^-2, and B = 0.026 T gives |ν| ≈ 15.9, not 6. The measured Hall value Rxy = h/6e^2 requires either n ≈ -3.75×10^9 cm^-2 at 26 mT or B ≈ 69 mT at -1×10^10 cm^-2. This is not a rounding issue: the fan diagram in Fig. 4a should show the Streda slope dn/dB = (e/h)ν, which is not presented for the first-developed plateau. Please correct the density/field values or provide a quantitative calibration check; as written, the central 'ν = -6 at ~20 mT' claim is not internally verifiable.
- [Methods: SWMcC model] The statement 'To better satisfy the experiment, we shift Δ by 5 meV' introduces a single fitted offset into the gap calculation. The main text claims 'quantitative agreement' between the measured gap versus D and the band-structure calculations (Fig. 2b and Fig. S10b), but this agreement is not parameter-free because of the offset. Please state the 5 meV shift explicitly in the main text, provide the unshifted comparison, and quantify how the D = 0 intrinsic gap and the asymmetry depend on this offset. This is important for assessing whether the observed gap is intrinsic to ABCBC stacking or is influenced by environmental doping/substrate potentials.
- [Low-field |v| = 6 quantum Hall state, final paragraph] The manuscript explicitly states 'we cannot distinguish between these two scenarios' (conventional zero-energy Landau level versus a Chern insulator from spontaneous valley polarization). Since the abstract and introduction present the low-field ν = -6 state as a highlight and attribute it to 'the interplay between spontaneous symmetry breaking and Berry curvature,' the paper should either provide a discriminating experimental test (e.g., temperature dependence of the plateau, measurement of the Hall resistance at zero B, or a Landau fan at lower temperatures) or clearly reframe the mechanism as an unresolved open question. As written, the concluding interpretation goes beyond what the data and analysis support.
minor comments (5)
- [Abstract and Fig. 4] The abstract states '~20 mT' while the text and Fig. 4b report 26 mT; please harmonize the numbers and specify the minimum field at which the h/6e^2 value is reached, and the level of quantization at higher fields.
- [Fig. S7 caption] The caption states that the |ν| = 6 state is 'confirmed by the Streda formula,' but no Streda slope or calibration plot is shown; please add the calculated dn/dB lines to the fan diagrams or provide the slope values.
- [Fig. 3c] The band-structure and Fermi-surface panels in Fig. 3c are not labeled with the Fermi energies or the regions I–V of Fig. 3b, which makes the correspondence between theory and experiment difficult to verify.
- [Notation] The manuscript interchangeably uses 'v = -6', '|v| = 6', and 'ν = -6'; please use a single convention (e.g., ν = -6) and define ν = n h/(e B) early in the text.
- [Fig. 4b] The Hall resistance data at low field may include a longitudinal admixture due to the small Hall angle; please show error bars or multiple sweeps, and comment on the uncertainty of the '93% quantization' value.
Circularity Check
No circular derivation: gap and QH claims rest on direct transport observables; the disclosed 5 meV Δ shift and the density calibration issue are non-circular limitations.
full rationale
The central transport claims are not derived from the theory they are compared with. The 0.65 meV gap at n = D = 0 is extracted from Arrhenius fits of the measured Rxx(T) (Fig. 2b and Fig. S4), and the asymmetric D dependence is a direct experimental observation reproduced in devices M1–M4; the DFT/SWMcC calculation is an independent single-particle model with standard parameters. The only fitting step is disclosed in Methods: "To better satisfy the experiment, we shift Δ by 5 meV." This is a horizontal calibration of the interlayer-potential axis; it does not define, fit, or rename the measured gap, and the non-monotonic asymmetry of the calculated gap versus Δ is a shape property of the independent band structure, not an artifact of that shift. The Lifshitz transitions are assigned from measured Landau-level degeneracies and compared with Fermi-surface calculations, with no fitted target value. The v = -6 quantum Hall state is identified from the Hall plateau Rxy = h/6e^2, an operational definition of filling factor, and the paper explicitly states it "cannot distinguish between these two scenarios," so the theoretical scenarios are not used to generate the observation. Self-citations are methodological (SNOM identification, previous transport on related stackings) and are not load-bearing for the new claims; no uniqueness theorem or ansatz is imported from the authors' own prior work. Separately, the quoted density for the v = -6 plateau, n = -1×10^10 cm^-2 at B = 26 mT with Rxy = h/6e^2, is internally inconsistent (ν ≈ 16 from the stated calibration), which is a data-consistency and correctness issue that the authors should correct; it is not a circularity because the filling factor is read from the measured Hall resistance rather than from the theory or from a fitted parameter.
Assumptions & free parameters
free parameters (1)
- Interlayer potential shift Delta_offset =
5 meV
assumptions (5)
- domain assumption Chiral decomposition and partitioning rules describe the low-energy band structure of mixed-stacked multilayer graphene (Min and MacDonald).
- domain assumption The SWMcC tight-binding Hamiltonian with the hopping and onsite parameters in Table 1 describes the relevant bands.
- domain assumption PBE-DFT and Wannier interpolation are accurate at the sub-meV to few-meV scale relevant to the 0.65 meV intrinsic gap.
- domain assumption The measured Hall bar region is predominantly a single ABCBC stacking domain.
- domain assumption The Arrhenius activation gap at charge neutrality is a faithful measure of a single-particle band gap.
Cite this review
Pith. "Pith review of Intrinsic layer polarization and multi-flatband transport in non-centrosymmetric mixed-stacked multilayer graphene." pith.science (2026). https://pith.science/paper/2TG3EAIN
@misc{pith2026250512478,
author = {Pith},
title = {Pith review of: Intrinsic layer polarization and multi-flatband transport in non-centrosymmetric mixed-stacked multilayer graphene},
year = {2026},
howpublished = {\url{https://pith.science/paper/2TG3EAIN}},
note = {Machine review of arXiv:2505.12478}
}
read the original abstract
Graphene multilayers exhibit electronic spectra that depend sensitively on both the number of layers and their stacking order. Beyond trilayer graphene, mixed stacking sequences (alternating Bernal and rhombohedral layers) give rise to multiple coexisting low-energy bands. Here we investigate ABCBC-stacked pentalayer graphene, a less-studied non-centrosymmetric mixed sequence. This stacking can be regarded as an ABC (rhombohedral) trilayer on top of an AB (Bernal) bilayer, so its low-energy band structure contains both a cubic band and a parabolic band that hybridize. In transport measurements, we observe an intrinsic band gap at charge neutrality whose magnitude changes asymmetrically under an applied perpendicular displacement field. This behavior reflects the spontaneous layer polarization inherent to the broken inversion symmetry and mirror symmetry. By tuning the displacement field and carrier density, we drive multiple Lifshitz transitions in the Fermi surface topology and realize Landau levels with different degeneracies arising from the multi-flatband system. Remarkably, a v = -6 quantum Hall state emerges at an exceptionally low magnetic field (~20 mT), indicating the interplay between spontaneous symmetry breaking and Berry curvatures. Our results establish mixed-stacked multilayer graphene as a tunable platform with various broken symmetries and multiple flatbands, suitable for exploring emergent correlated electronic states.
Figures
Figures from the paper (1 more)
Forward citations
Cited by 1 Pith paper
-
Cascade of Even-Denominator Fractional Quantum Hall States in Mixed-Stacked Multilayer Graphene
Even-denominator fractional quantum Hall states up to filling -13/2 are observed in mixed-stacked pentalayer graphene, with numerical evidence suggesting Moore-Read type non-Abelian order.
Reference graph
Works this paper leans on
-
[43]
& Furthmüller, J
Kresse, G. & Furthmüller, J. Efficient iterative schemes for ab initio total-energy calculations using a plane-wave basis set. Phys. Rev. B 54, 11169–11186 (1996)
1996
- [44]
- [45]
-
[46]
P., Burke, K
Perdew, J. P., Burke, K. & Ernzerhof, M. Generalized Gradient Approximation Made Simple. Phys. Rev. Lett. 77, 3865–3868 (1996)
1996
-
[47]
& Vanderbilt, D
Marzari, N. & Vanderbilt, D. Maximally localized generalized Wannier functions for composite energy bands. Phys. Rev. B 56, 12847–12865 (1997)
1997
-
[48]
An updated version of wannier90: A tool for obtaining maximally -localised Wannier functions. Comput. Phys. Commun. 185, 2309–2310 (2014)
work page 2014
-
[49]
WannierTools: An open-source software package for novel topological materials. Comput. Phys. Commun. 224, 405–416 (2018)
work page 2018
-
[50]
McClure, J. W. Band Structure of Graphite and de Haas-van Alphen Effect. Phys. Rev. 108, 612–618 (1957)
work page 1957
Show all 12 references
-
[51]
McClure, J. W. Theory of Diamagnetism of Graphite. Phys. Rev. 119, 606 –613 (1960)
1960
-
[52]
Slonczewski, J. C. & Weiss, P. R. Band Structure of Graphite. Phys. Rev. 109, 272– 279 (1958)
1958
-
[53]
J., Garcia -Ruiz, A
Sarsfield, P. J., Garcia -Ruiz, A. & Fal’ko, V . I. Substrate, temperature, and magnetic field dependence of electric polarization in mixed -stacking tetralayer graphenes. Phys. Rev. Res. 6, 043324 (2024)
2024
-
[54]
Boschi, A. et al. Built-in Bernal gap in large -angle-twisted monolayer -bilayer graphene. Commun. Phys. 7, 391 (2024). Acknowledgments G.C. acknowledges support from NSF of China (grant nos. 12350005 and 12174248), National Key Research Program of China (grant nos. 2021YFA140...
2024
Reviewed August 15, 2026 · model on record in the stance chip above.
Discussion (0). Continue with ORCID to comment.