Pith. sign in

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 →

arxiv 2505.12478 v2 pith:2TG3EAIN submitted 2025-05-18 cond-mat.mes-hall

classification cond-mat.mes-hall PACS 73.43.-f73.22.Pr71.70.Di
keywords ABCBCpentalayergraphenemixedstackingnon-centrosymmetriclayerpolarizationintrinsicbandgapLifshitztransitionquantumHalleffectflatbands
verification ladder T0 review T1 audit T2 compute T3 formal

The pith

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The reading

ABCBC-stacked pentalayer graphene, a non-centrosymmetric mixed stacking of an ABC trilayer and an AB bilayer, is shown to have an intrinsic band gap of about 0.65 meV at charge neutrality even when no displacement field is applied. The paper reports that this gap responds asymmetrically and non-monotonically to a perpendicular displacement field, closing on one side and growing then shrinking on the other, in agreement with tight-binding calculations. Gate tuning drives multiple Lifshitz transitions and produces Landau levels with degeneracies 4, 8, and 12, reflecting coexisting cubic and parabolic low-energy bands. A quantum Hall plateau at filling factor $\nu = -6$ appears at about 26 mT with roughly 93% of $h/6e^2$ quantization. The interest is that a natural, twist-free graphene multilayer can combine broken inversion symmetry with multiple flat bands, offering a tunable platform for correlated and topological states.

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.

Watch

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

Editorial extensions of the paper, not claims the author makes directly.

  • 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.
Share X Bluesky LinkedIn Reddit HN

Signed reviews

No signed human review yet.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, and a circularity audit.

Referee Report

3 major / 5 minor

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)
  1. [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.
  2. [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.
  3. [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)
  1. [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.
  2. [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.
  3. [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.
  4. [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.
  5. [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

0 steps flagged · score 0.0 of 10

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 1 free parameters · 5 assumptions · 0 invented entities

The central experimental claims do not depend on invented entities; the main model input is the SWMcC Hamiltonian with parameters from DFT, plus a 5 meV shift in the interlayer potential chosen to match experiment. This shift is a free parameter for the quantitative gap comparison.

free parameters (1)
  • Interlayer potential shift Delta_offset = 5 meV
    Methods: 'To better satisfy the experiment, we shift Delta by 5 meV.' This ad hoc offset is added to the electric-field-induced on-site energies to make the calculated gap-D behavior match the transport data; it is not derived from DFT or first principles.
assumptions (5)
  • domain assumption Chiral decomposition and partitioning rules describe the low-energy band structure of mixed-stacked multilayer graphene (Min and MacDonald).
    Used in Fig. 1a to decompose ABCBC into an ABC trilayer and AB bilayer, and to interpret the low-energy bands.
  • domain assumption The SWMcC tight-binding Hamiltonian with the hopping and onsite parameters in Table 1 describes the relevant bands.
    All band structure, Fermi surface, and gap-D calculations use this model; the parameters are fitted to DFT and partly inherited from prior literature.
  • 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.
    The calculated gap versus D is compared quantitatively with the transport activation gap; no beyond-PBE corrections are applied.
  • domain assumption The measured Hall bar region is predominantly a single ABCBC stacking domain.
    Stacking is identified by SNOM contrast and NanoARPES; however, M4 shows signs of multiple doping domains (Fig. S3), so single-domain purity is assumed rather than directly proven for every device.
  • domain assumption The Arrhenius activation gap at charge neutrality is a faithful measure of a single-particle band gap.
    Thermally activated Rxx is used to extract gaps; disorder, electron-hole puddles, or interaction effects could in principle contribute.

how reviews work

0 comments
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 reproduced from arXiv: 2505.12478 by the authors.

Figure 1
Figure 1. Possible stacking orders of pentalayer graphene and th [PITH_FULL_IMAGE:figures/full_fig_p018_1.png] view at source ↗
Figure 2
Figure 2. Symmetry breaking and layer polarization of ABCBC. a, [PITH_FULL_IMAGE:figures/full_fig_p019_2.png] view at source ↗
Figure 3
Figure 3. Multiple Lifshitz transitions under D. a, n - D color plot of resistance at B = 1 T. The color bar is in the log scale. b, Corresponding phase diagram of a according to the Landau level degeneracy in each region. c, Calculated single particle band structures and the Fermi surfaces at different Fermi energy. Conduction and valence bands are drawn as red and blue respectively [PITH_FULL_IMAGE:figures/full_fig_p020_3.png] view at source ↗
Figures from the paper (1 more)
Figure 4
Figure 4. Figure 4: Robust |v| = 6 quantum hall state. a, R – n – B color plot at D = 0 when T = 1.5 K for device M1. b, Corresponding Hall resistance Rxy and Rxx as a function of magnetic field measured at D = 0, n = -1×1010 cm-2, T = 14 mK. c, R – n – D color plot at B = 1 T when T = 1.…

Discussion (0). Continue with ORCID to comment.

Forward citations

Cited by 1 Pith paper

Reviewed papers in the Pith corpus that reference this work. Sorted by Pith novelty score. Full citation record

  1. Cascade of Even-Denominator Fractional Quantum Hall States in Mixed-Stacked Multilayer Graphene

    cond-mat.mes-hall 2025-07 conditional novelty 7.0 of 10

    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

12 extracted references · 9 canonical work pages · cited by 1 Pith paper

  1. [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)

  2. [44]

    & Kohn, W

    Hohenberg, P. & Kohn, W. Inhomogeneous Electron Gas. Phys. Rev. 136, B864– B871 (1964)

  3. [45]

    & Sham, L

    Kohn, W. & Sham, L. J. Self -Consistent Equations Including Exchange and Correlation Effects. Phys. Rev. 140, A1133–A1138 (1965)

  4. [46]

    P., Burke, K

    Perdew, J. P., Burke, K. & Ernzerhof, M. Generalized Gradient Approximation Made Simple. Phys. Rev. Lett. 77, 3865–3868 (1996)

  5. [47]

    & Vanderbilt, D

    Marzari, N. & Vanderbilt, D. Maximally localized generalized Wannier functions for composite energy bands. Phys. Rev. B 56, 12847–12865 (1997)

  6. [48]

    An updated version of wannier90: A tool for obtaining maximally -localised Wannier functions. Comput. Phys. Commun. 185, 2309–2310 (2014)

  7. [49]

    WannierTools: An open-source software package for novel topological materials. Comput. Phys. Commun. 224, 405–416 (2018)

  8. [50]

    McClure, J. W. Band Structure of Graphite and de Haas-van Alphen Effect. Phys. Rev. 108, 612–618 (1957)

Show all 12 references
  1. [51]

    McClure, J. W. Theory of Diamagnetism of Graphite. Phys. Rev. 119, 606 –613 (1960)

  2. [52]

    Slonczewski, J. C. & Weiss, P. R. Band Structure of Graphite. Phys. Rev. 109, 272– 279 (1958)

  3. [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)

  4. [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...

Pith tools

Reviewed August 15, 2026 · model on record in the stance chip above.