REVIEW 2 major objections 1 minor 1 cited by
In partially isospin-polarized bilayer graphene, minority carriers form one or two Fermi pockets instead of three, producing spontaneous nematic order that breaks C3 symmetry.
Reviewed by Pith at T0; open to challenge. T0 means a machine referee read the full paper against a public rubric. the ladder, T0–T4 →
Hartree-Fock phase diagram of Bernal bilayer graphene shows nematic order with one or two minority-carrier pockets in partially polarized regimes, controlled by screened interaction strength, with an analytical criterion matching numerics.
T0 review reviewed 2026-06-29 challenge →
load-bearing objection HF calculation finds nematic pocket reduction in partially polarized bilayer graphene near the polarization boundary, but the mean-field setup likely over-stabilizes the order. the 2 major comments →
Nematic Phase Transitions in Multilayer Graphene Systems
The pith
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
The reading
Core claim
Within these partially isospin polarized states, we identify regions with three C3-symmetric pockets for the minority carriers and regions with only one or two such pockets, implying that the system develops a spontaneous nematic order that breaks C3 symmetry. We find numerically that nematicity emerges near the boundary between fully and partially polarized phases and is controlled by the strength of the screened interaction. We analytically derive a criterion for the nematic order, which agrees with our numerical results.
What carries the argument
The spontaneous reduction from three C3-symmetric minority-carrier Fermi pockets to one or two, which breaks rotational symmetry while the majority carriers remain unaffected.
Load-bearing premise
The self-consistent Hartree-Fock approximation with the chosen dual-gate screened interaction and trigonal warping term is sufficient to locate the nematic transition without fluctuation corrections.
What would settle it
Quantum-oscillation measurements that continue to resolve three distinct minority pockets throughout the partially polarized density window would rule out the predicted spontaneous symmetry breaking.
If this is right
- Nematic order appears only near the boundary separating fully polarized and partially polarized phases.
- Increasing the strength of the screened Coulomb interaction enlarges the density range where nematicity is stable.
- The derived analytical criterion for pocket merging reproduces the numerically obtained phase boundaries.
- The overall sequence of isospin-polarized states remains consistent with existing quantum-oscillation data.
Where Pith is reading between the lines
- The same pocket-reduction mechanism may operate in rhombohedral trilayer graphene once the corresponding Hartree-Fock calculation is performed.
- Nematic order could modify the pairing landscape that produces the superconducting domes reported in the same devices.
- Gate-tunable displacement fields offer a direct experimental knob to cross the nematic boundary and observe the change in oscillation frequencies.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The manuscript investigates the emergence of nematic order in multilayer graphene systems (with emphasis on Bernal bilayer graphene) using self-consistent Hartree-Fock calculations that incorporate trigonal warping in the band structure and a dual-gate screened Coulomb interaction. It reports a sequence of isospin-polarized phases as functions of carrier density and displacement field, including partially polarized states with large majority and small minority Fermi pockets. Within the partially polarized regime, regions are identified where the minority carriers occupy only one or two pockets rather than three C3-symmetric pockets, indicating spontaneous C3-breaking nematic order that appears near the boundary with fully polarized phases and is controlled by interaction strength. An analytical criterion for the onset of nematicity is derived that agrees with the numerical results.
Significance. If the central claim holds, the work provides a concrete microscopic mechanism linking observed nematicity in graphene multilayers to the boundary between fully and partially isospin-polarized states, with the screened interaction as the tuning parameter. The combination of numerical self-consistent HF results and a matching analytical criterion derived within the same framework is a positive feature, as is the use of a realistic band structure and interaction form that reproduces experimental sequences of polarized phases.
major comments (2)
- [Methods / Numerical results (near polarization boundary)] The central claim that spontaneous nematic order (manifested as reconfiguration to 1- or 2-pocket minority Fermi surfaces) emerges near the fully-to-partially polarized boundary relies entirely on self-consistent Hartree-Fock. Mean-field approximations of this type are known to overestimate the stability of broken-symmetry states, especially in the vicinity of a polarization boundary where the energy scale separating ordered and disordered configurations is small and where the paper states that nematicity is controlled by interaction strength. No fluctuation corrections, functional RG analysis, or alternative variational ansatz is presented to test whether the pocket reconfiguration survives beyond HF.
- [Analytical derivation] The analytical criterion for nematic order is obtained under the same mean-field decoupling used in the numerics and therefore inherits the same limitation; it would be necessary to show either that the criterion can be derived or tested in a controlled expansion beyond HF or that it reduces to a known result in a limit where fluctuations are demonstrably unimportant.
minor comments (1)
- [Abstract / Results] The abstract and main text should explicitly state the range of interaction strengths (or gate parameters) over which the nematic regions are found, to allow direct comparison with the experimental quantum-oscillation data cited.
Simulated Author's Rebuttal
We thank the referee for their careful reading and constructive comments on our manuscript. We respond point-by-point to the major comments below.
read point-by-point responses
-
Referee: [Methods / Numerical results (near polarization boundary)] The central claim that spontaneous nematic order (manifested as reconfiguration to 1- or 2-pocket minority Fermi surfaces) emerges near the fully-to-partially polarized boundary relies entirely on self-consistent Hartree-Fock. Mean-field approximations of this type are known to overestimate the stability of broken-symmetry states, especially in the vicinity of a polarization boundary where the energy scale separating ordered and disordered configurations is small and where the paper states that nematicity is controlled by interaction strength. No fluctuation corrections, functional RG analysis, or alternative variational ansatz is presented to test whether the pocket reconfiguration survives beyond HF.
Authors: We agree that the self-consistent Hartree-Fock approximation can overestimate the stability of broken-symmetry states near phase boundaries. All results in the manuscript, including the identification of nematic regions with 1- or 2-pocket minority Fermi surfaces, are obtained within this framework. The nematic reconfiguration is found to be controlled by screened interaction strength and appears consistently near the fully-to-partially polarized boundary. We will revise the manuscript to include an explicit discussion of the limitations of the mean-field approach. revision: partial
-
Referee: [Analytical derivation] The analytical criterion for nematic order is obtained under the same mean-field decoupling used in the numerics and therefore inherits the same limitation; it would be necessary to show either that the criterion can be derived or tested in a controlled expansion beyond HF or that it reduces to a known result in a limit where fluctuations are demonstrably unimportant.
Authors: The analytical criterion is derived using the identical Hartree-Fock decoupling employed in the self-consistent numerics, and its validity is demonstrated by direct agreement with those numerical results. We do not derive or test the criterion in an expansion beyond mean-field theory, as such an analysis lies outside the scope of the present work. revision: no
- Testing whether the nematic pocket reconfiguration survives beyond the Hartree-Fock approximation via fluctuation corrections, functional RG, or alternative variational methods.
- Deriving or testing the analytical criterion for nematic order in a controlled expansion beyond Hartree-Fock.
Circularity Check
No circularity: numerical HF phase diagram and independent analytic criterion are self-contained
full rationale
The paper computes the phase diagram via self-consistent Hartree-Fock with fixed trigonal warping and dual-gate screened Coulomb interaction, then identifies nematic regions by counting minority-carrier pockets (three C3-symmetric vs. one or two). An analytic criterion for nematicity is derived from the same mean-field equations and shown to match the numerics. Neither step reduces by construction to a fitted parameter renamed as prediction, nor relies on self-citation chains, uniqueness theorems imported from the authors, or ansatze smuggled via prior work. The derivation chain is externally falsifiable against the interaction strength and displacement field and does not collapse to its own inputs.
Axiom & Free-Parameter Ledger
axioms (1)
- domain assumption Hartree-Fock approximation suffices to capture nematic order and phase boundaries
Cite this review
Pith. "Pith review of Nematic Phase Transitions in Multilayer Graphene Systems." pith.science (2026). https://pith.science/paper/37VEIPDL
@misc{pith2026260526270,
author = {Pith},
title = {Pith review of: Nematic Phase Transitions in Multilayer Graphene Systems},
year = {2026},
howpublished = {\url{https://pith.science/paper/37VEIPDL}},
note = {Machine review of arXiv:2605.26270}
}
abstract
Recent experiments on graphene multilayers under displacement field have demonstrated a wide variety of electronically ordered phases, including valley and/or spin polarized phases as well as potentially unconventional superconducting phases. In addition, quantum oscillation measurements in Bernal bilayer graphene and rhombohedral trilayer graphene showed the presence of electronic nematic states. Here, we investigate the emergence of nematic order in these systems, with emphasis on Bernal bilayer graphene, within a self-consistent Hartree--Fock framework, using a realistic band structure with trigonal warping and a dual-gate screened Coulomb interaction. We compute the phase diagram as carrier density and displacement field are varied, and find a sequence of isospin (spin and valley) polarized states consistent with experiment, including partially polarized phases with a large Fermi pockets for isospin-majority carriers and smaller Fermi pockets for isospin-minority carriers. Within these partially isospin polarized states, we identify regions with three $C_3$-symmetric pockets for the minority carriers and regions with only one or two such pockets, implying that the system develops a spontaneous nematic order that breaks $C_3$ symmetry. We find numerically that nematicity emerges near the boundary between fully and partially polarized phases and is controlled by the strength of the screened interaction. We analytically derive a criterion for the nematic order, which agrees with our numerical results.
Figures
Forward citations
Cited by 1 Pith paper
-
Symmetry-protected cubic-touching topological surface bands with tunable singularities
Topological crystalline insulator surfaces can host symmetry-protected cubic band touchings that continuously tune between cubic, moat, and six-valley dispersions, with density-of-state singularities from power law to...
Reference graph
Works this paper leans on
-
[1]
H. Zhou, T. Xie, A. Ghazaryan, T. Holder, J. R. Ehrets, E. M. Spanton, T. Taniguchi, K. Watanabe, E. Berg, M. Serbyn, and A. F. Young, Half- and quarter-metals in rhombohedral trilayer graphene, Nature598, 429 (2021)
2021
-
[2]
H. Zhou, T. Xie, T. Taniguchi, K. Watanabe, and A. F. Young, Superconductivity in rhom- bohedral trilayer graphene, Nature598, 434 (2021)
2021
-
[3]
H. Zhou, L. Holleis, Y. Saito, L. Cohen, W. Huynh, C. L. Patterson, F. Yang, T. Taniguchi, K. Watanabe, and A. F. Young, Isospin magnetism and spin- polarized superconductivity in bernal bilayer graphene, Science375, 774 (2022), https://www.science.org/doi/pdf/10.1126/science.abm8386
-
[4]
A. M. Seiler, F. R. Geisenhof, F. Winterer, K. Watanabe, T. Taniguchi, T. Xu, F. Zhang, and R. T. Weitz, Quantum cascade of correlated phases in trigonally warped bilayer graphene, Nature608, 298 (2022); A. M. Seiler, Y. Zhumagulov, K. Zollner, C. Yoon, D. Urbaniak, F. R. 23 Geisenhof, K. Watanabe, T. Taniguchi, J. Fabian, F. Zhang, and R. T. Weitz, Layer...
2022
-
[5]
S. C. de la Barrera, S. Aronson, Z. Zheng, K. Watanabe, T. Taniguchi, Q. Ma, P. Jarillo- Herrero, and R. Ashoori, Cascade of isospin phase transitions in bernal-stacked bilayer graphene at zero magnetic field, Nature Physics18, 771 (2022)
2022
-
[6]
T. Han, Z. Lu, G. Scuri, J. Sung, J. Wang, T. Han, K. Watanabe, T. Taniguchi, L. Fu, H. Park, and L. Ju, Orbital multiferroicity in pentalayer rhombohedral graphene, Nature 623, 41 (2023)
2023
-
[7]
T. Arp, O. Sheekey, H. Zhou, C. L. Tschirhart, C. L. Patterson, H. M. Yoo, L. Holleis, E. Redekop, G. Babikyan, T. Xie, J. Xiao, Y. Vituri, T. Holder, T. Taniguchi, K. Watanabe, M. E. Huber, E. Berg, and A. F. Young, Intervalley coherence and intrinsic spin–orbit coupling in rhombohedral trilayer graphene, Nature Physics20, 1413 (2024)
2024
-
[8]
T. Han, Z. Lu, G. Scuri, J. Sung, J. Wang, T. Han, K. Watanabe, T. Taniguchi, H. Park, and L. Ju, Correlated insulator and chern insulators in pentalayer rhombohedral-stacked graphene, Nature Nanotechnology19, 181 (2024)
2024
-
[9]
Holleis, C
L. Holleis, C. L. Patterson, Y. Zhang, Y. Vituri, H. M. Yoo, H. Zhou, T. Taniguchi, K. Watan- abe, E. Berg, S. Nadj-Perge, and A. F. Young, Nematicity and orbital depairing in supercon- ducting bernal bilayer graphene, Nature Physics21, 444 (2025)
2025
-
[10]
N. Auerbach, S. Dutta, M. Uzan, Y. Vituri, Y. Zhou, A. Y. Meltzer, S. Grover, T. Holder, P. Emanuel, M. E. Huber, Y. Myasoedov, K. Watanabe, T. Taniguchi, Y. Oreg, E. Berg, and E. Zeldov, Visualizing isospin magnetic texture and intervalley exchange interaction in rhombohedral tetralayer graphene (2025), arXiv:2503.14146 [cond-mat.mes-hall]
- [11]
-
[12]
T. Han, Z. Lu, Z. Hadjri, L. Shi, Z. Wu, W. Xu, Y. Yao, A. A. Cotten, O. Sharifi Sedeh, H. Weldeyesus, J. Yang, J. Seo, S. Ye, M. Zhou, H. Liu, G. Shi, Z. Hua, K. Watanabe, T. Taniguchi, P. Xiong, D. M. Zumb¨ uhl, L. Fu, and L. Ju, Signatures of chiral superconduc- tivity in rhombohedral graphene, Nature643, 654 (2025)
2025
- [13]
-
[14]
Zhang, R
Y. Zhang, R. Polski, A. Thomson, ´E. Lantagne-Hurtubise, C. Lewandowski, H. Zhou, K. Watanabe, T. Taniguchi, J. Alicea, and S. Nadj-Perge, Enhanced superconductivity in spin–orbit proximitized bilayer graphene, Nature613, 268 (2023)
2023
-
[15]
C. Li, F. Xu, B. Li, J. Li, G. Li, K. Watanabe, T. Taniguchi, B. Tong, J. Shen, L. Lu, J. Jia, F. Wu, X. Liu, and T. Li, Tunable superconductivity in electron- and hole-doped bernal bilayer graphene, Nature631, 300 (2024)
2024
-
[16]
C. L. Patterson, O. I. Sheekey, T. B. Arp, L. F. W. Holleis, J. M. Koh, Y. Choi, T. Xie, S. Xu, Y. Guo, H. Stoyanov, E. Redekop, C. Zhang, G. Babikyan, D. Gong, H. Zhou, X. Cheng, T. Taniguchi, K. Watanabe, M. E. Huber, C. Jin, ´E. Lantagne-Hurtubise, J. Alicea, and A. F. Young, Superconductivity and spin canting in spin–orbit-coupled trilayer graphene, N...
2025
-
[17]
McCann and V
E. McCann and V. I. Fal’ko, Landau-level degeneracy and quantum hall effect in a graphite bilayer, Phys. Rev. Lett.96, 086805 (2006)
2006
-
[18]
McCann, Asymmetry gap in the electronic band structure of bilayer graphene, Phys
E. McCann, Asymmetry gap in the electronic band structure of bilayer graphene, Phys. Rev. B74, 161403 (2006)
2006
-
[19]
McCann and M
E. McCann and M. Koshino, The electronic properties of bilayer graphene, Reports on Progress in Physics76, 056503 (2013)
2013
-
[20]
Jung and A
J. Jung and A. H. MacDonald, Accurate tight-binding models for theπbands of bilayer graphene, Phys. Rev. B89, 035405 (2014)
2014
-
[21]
D. V. Chichinadze, L. Classen, Y. Wang, and A. V. Chubukov, Cascade of transitions in twisted and non-twisted graphene layers within the van hove scenario, npj Quantum Materials7, 114 (2022); Su(4) symmetry in twisted bilayer graphene: An itinerant perspective, Phys. Rev. Lett.128, 227601 (2022)
2022
-
[22]
H. Hu, B. A. Bernevig, and A. M. Tsvelik, Kondo lattice model of magic-angle twisted-bilayer graphene: Hund’s rule, local-moment fluctuations, and low-energy effective theory, Phys. Rev. Lett.131, 026502 (2023)
2023
-
[23]
Xie and S
M. Xie and S. Das Sarma, Flavor symmetry breaking in spin-orbit coupled bilayer graphene, Phys. Rev. B107, L201119 (2023)
2023
-
[24]
Y.-C. Lee, D. V. Chichinadze, and A. V. Chubukov, Crossover from ordinary to higher order 25 van hove singularity in a honeycomb system: A parquet renormalization group analysis, Phys. Rev. B109, 155118 (2024)
2024
-
[25]
J. M. Koh, J. Alicea, and E. Lantagne-Hurtubise, Correlated phases in spin-orbit-coupled rhombohedral trilayer graphene, Phys. Rev. B109, 035113 (2024)
2024
-
[26]
J. M. Koh, A. Thomson, J. Alicea, and E. Lantagne-Hurtubise, Symmetry-broken metallic orders in spin-orbit-coupled bernal bilayer graphene, Phys. Rev. B110, 245118 (2024)
2024
-
[27]
T. Wang, M. Vila, M. P. Zaletel, and S. Chatterjee, Electrical control of spin and valley in spin-orbit coupled graphene multilayers, Phys. Rev. Lett.132, 116504 (2024)
2024
-
[28]
Friedlan, H
A. Friedlan, H. Li, and H.-Y. Kee, Valley polarization, magnetization, and superconductivity in bilayer graphene near the van hove singularity, Phys. Rev. B111, 024504 (2025)
2025
-
[29]
R. D. Mayrhofer and A. V. Chubukov, Valley- and spin-polarized states in bernal bilayer graphene, Phys. Rev. B111, 245114 (2025)
2025
-
[30]
G. Parra-Mart´ ınez, A. Jimeno-Pozo, J. A. Silva-Guill´ en, and F. Guinea, Nematic and partially polarized phases in rhombohedral graphene with varying number of layers: An extensive hartree–fock study, Advanced Materialsn/a, e14008 (2025), https://advanced.onlinelibrary.wiley.com/doi/pdf/10.1002/adma.202514008
-
[31]
Ghazaryan, T
A. Ghazaryan, T. Holder, M. Serbyn, and E. Berg, Unconventional superconductivity in sys- tems with annular fermi surfaces: Application to rhombohedral trilayer graphene, Phys. Rev. Lett.127, 247001 (2021)
2021
-
[32]
Ghazaryan, T
A. Ghazaryan, T. Holder, E. Berg, and M. Serbyn, Multilayer graphenes as a platform for interaction-driven physics and topological superconductivity, Phys. Rev. B107, 104502 (2023)
2023
-
[33]
P. A. Pantale´ on, A. Jimeno-Pozo, H. Sainz-Cruz, V. Phong, T. Cea, and F. Guinea, Supercon- ductivity and correlated phases in non-twisted bilayer and trilayer graphene, Nature Reviews Physics5, 304 (2023)
2023
-
[34]
Z. Dong, A. V. Chubukov, and L. Levitov, Transformer spin-triplet superconductivity at the onset of isospin order in bilayer graphene, Phys. Rev. B107, 174512 (2023)
2023
-
[35]
Z. Dong, P. A. Lee, and L. S. Levitov, Signatures of cooper pair dynamics and quantum- critical superconductivity in tunable carrier bands, Proceedings of the National Academy of Sciences120, e2305943120 (2023), https://www.pnas.org/doi/pdf/10.1073/pnas.2305943120
-
[36]
Z. Dong, L. Levitov, and A. V. Chubukov, Superconductivity near spin and valley orders in graphene multilayers, Phys. Rev. B108, 134503 (2023). 26
2023
-
[37]
Z. M. Raines and A. V. Chubukov, Superconductivity induced by spin-orbit coupling in a two-valley ferromagnet, npj Quantum Materials11, 31 (2026)
2026
-
[38]
Z. M. Raines and A. V. Chubukov, Superconductivity via paramagnon and magnon exchange in a 2d near-ferromagnetic full metal and ferromagnetic half-metal, Phys. Rev. X16, 011059 (2026)
2026
-
[39]
X. Han, Q. Liu, Y. Wang, R. Niu, Z. Qu, Z. Wang, Z. Li, C. Han, K. Watanabe, T. Taniguchi, Z. Song, J. Liu, J. Mao, Z. Han, B. L. Chittari, J. Jung, Z. Gan, and J. Lu, Engineering the band topology in a rhombohedral trilayer graphene moir´ e superlattice, Nano Letters24, 6286 (2024), pMID: 38747346, https://doi.org/10.1021/acs.nanolett.4c00948
-
[40]
Z. M. Raines and A. V. Chubukov, Two-dimensional stoner transitions beyond mean field, Phys. Rev. B110, 235433 (2024)
2024
-
[41]
Z. M. Raines, L. I. Glazman, and A. V. Chubukov, Unconventional discontinuous transitions in isospin systems, Phys. Rev. Lett.133, 146501 (2024); Unconventional discontinuous tran- sitions in a two-dimensional system with spin and valley degrees of freedom, Phys. Rev. B 110, 155402 (2024)
2024
-
[42]
Scholle and L
R. Scholle and L. Classen, Unrestricted hartree-fock analysis of the hubbard model on single- layer graphene and bernal bilayer graphene, Phys. Rev. B113, 085132 (2026)
2026
-
[43]
Huang, T
C. Huang, T. M. R. Wolf, W. Qin, N. Wei, I. V. Blinov, and A. H. MacDonald, Spin and orbital metallic magnetism in rhombohedral trilayer graphene, Phys. Rev. B107, L121405 (2023)
2023
-
[44]
Z. Dong, M. Davydova, O. Ogunnaike, and L. Levitov, Isospin- and momentum-polarized orders in bilayer graphene, Phys. Rev. B107, 075108 (2023)
2023
-
[45]
Calvera, A
V. Calvera, A. Valenti, S. D. Huber, E. Berg, and S. A. Kivelson, Theory of coulomb driven nematicity in a multivalley two-dimensional electron gas, Phys. Rev. B111, 155135 (2025)
2025
-
[46]
In this work, we will consider the case where the sample is not in the proximity of any TMD, and the spin-orbit coupling can be safely neglected
The spin orbit coupling may be enhanced by placing the sample in proximity of a TMD like WSe2. In this work, we will consider the case where the sample is not in the proximity of any TMD, and the spin-orbit coupling can be safely neglected. 27
This paper was first reviewed by grok-4.3 on June 29, 2026.
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.