Pith. sign in

REVIEW 4 major objections 4 minor 106 references

Wigner crystals in tetralayer rhombohedral graphene can spontaneously break threefold rotation symmetry, and the reason is kinetic: the electron solid concentrates its momentum occupation in one lobe of the trigonally warped Fermi surface.

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 →

T0 review · deepseek-v4-flash

2026-08-01 19:14 UTC pith:KID7XPBU

load-bearing objection A new zero-field nematic Wigner crystal prediction in tetralayer graphene; the physics is plausible and the numerics mostly careful, but the stability labeling has a few internal inconsistencies that need fixing. the 4 major comments →

arxiv 2607.17027 v1 pith:KID7XPBU submitted 2026-07-19 cond-mat.str-el cond-mat.mes-hall

Nematic Wigner crystals in rhombohedral multilayer graphene

classification cond-mat.str-el cond-mat.mes-hall PACS 71.10.-w73.20.Qt
keywords Wigner crystalnematicityspontaneous C3 breakingrhombohedral multilayer grapheneHartree-Focktime-dependent Hartree-Focktrigonal warpingelectron solid
verification ladder T0 review T1 audit T2 compute T3 formal T4 reserved

The pith

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

The paper asks whether a Wigner crystal — the ordered electron solid that Coulomb repulsion creates at low density — can spontaneously break the threefold rotational symmetry of the host lattice, forming a nematic Wigner crystal with no magnetic field and no external anisotropy. Working in rhombohedral tetralayer graphene, the authors claim the answer is yes. Using Hartree-Fock and time-dependent Hartree-Fock calculations projected onto the lowest conduction band, they find a stable, spin- and valley-polarized nematic Wigner crystal in a sizable region of the phase diagram at high displacement fields and low electron densities, with a second smaller region at intermediate fields and higher densities. The stated reason the nematic crystal wins is kinetic, not interaction-based: by placing its momentum occupation in one of the three lobes of the strongly trigonally warped Fermi surface, it saves roughly 0.3 meV per electron over the symmetric triangular crystal. The claim matters because a spontaneously nematic electronic crystal — without magnetic field or engineered anisotropy — has not been predicted before, and the paper supplies two concrete experimental signatures: anisotropic charge density in STM and a direction-dependent depinning voltage in angle-resolved transport.

Core claim

The central claim: in rhombohedral tetralayer graphene at low density and high displacement field, the Wigner crystal spontaneously breaks the threefold rotation C3. The Hartree-Fock ground state, e.g. at n_e=0.25×10^12 cm^-2 and V_disp=90 meV, concentrates momentum occupation in one lobe of the trigonally warped Fermi sea, beating the C3-invariant triangular crystal by a kinetic-energy gain of ~0.3 meV per electron despite higher interaction energy. Time-dependent Hartree-Fock finds all C3-invariant and half of the nematic crystals locally stable; the rest soften toward two-electron-per-cell Wigner crystals or metallic states. All are spin-valley-polarized, zero-Chern-number insulators.

What carries the argument

The engine is the three-lobed, trigonally warped Fermi surface of the lowest band of tetralayer rhombohedral graphene. Projected Hartree-Fock on a single shared momentum mesh lets the authors compare Wigner lattice shapes in energy at fixed density; time-dependent Hartree-Fock's stability matrix tests local stability at each Wigner Bloch momentum. Two diagnostics define nematicity: C3 violation of the one-body density matrix, and hyperbolic lattice distance to the perfect triangular lattice. The decisive object is momentum occupation n(k): concentrating it in one lobe lowers kinetic energy by ~0.3 meV per electron, which is what makes the nematic crystal win.

Load-bearing premise

The load-bearing premise is that projected Hartree-Fock plus TDHF captures the true ground state: only the lowest conduction band is kept (justified by an energy ratio, not exact calculation), only one-electron-per-cell crystals are converged directly, and large runs use only fully spin-valley-polarized seeds; if any of these fail — or the screening (ε_r=5, 50 nm gates) misrepresents the device — the nematic regions shift or disappear.

What would settle it

A direct falsifier: STM images in the predicted nematic region (n_e ≈ 0.25×10^12 cm^-2, V_disp ≈ 90 meV) showing a C3-symmetric triangular charge pattern, or angle-resolved transport showing identical depinning voltages in all in-plane directions. On the theory side, a converged Hartree-Fock run at N_e=144 seeded with unpolarized spin-valley states — or an exact/quantum Monte Carlo calculation in the same continuum model — that finds the C3-invariant crystal lower in energy would refute the kinetic-energy mechanism.

Watch this falsifier — get emailed when new claim-graph text bears on it.

If this is right

  • A stable nematic Wigner crystal should exist in rhombohedral tetralayer graphene at low density and high displacement field, identifiable by its anisotropic charge density in STM.
  • Angle-resolved transport in a multi-contact geometry should show a direction-dependent depinning voltage in the nematic region — a signature absent from the C3-invariant Wigner crystal.
  • All Wigner crystals in the studied phase diagram are insulating, spin- and valley-polarized, with zero Chern number and indirect gaps above 10 meV; the C3-invariant triangular crystal wins at low displacement fields.
  • Roughly half of the nematic Hartree-Fock solutions are locally unstable, softening toward two-electron-per-cell Wigner crystals or metallic, stripe, and Fermi-liquid states, which shrinks the stable nematic region.
  • The same three-lobed Fermi-surface mechanism should produce nematic Wigner crystals in other rhombohedral multilayers wherever the trigonally warped Fermi sea exists.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

  • My inference: because the mechanism is kinetic, nematic stability should track the strength of trigonal warping; layer number (trilayer versus pentalayer versus hexalayer) and strain should tune the nWC region continuously — an expectation the paper gestures at but does not quantify.
  • My inference: the ~0.3 meV per electron kinetic gain is small relative to realistic disorder and screening uncertainty, so the phase boundaries are the most fragile part of the claim; a gate-tunable transport sweep across V_disp should map the predicted nematic region and can be compared directly to the paper's phase diagram.
  • My inference: the instabilities toward two-electron-per-cell Wigner crystals suggest a paired Wigner crystal may be the true ground state in parts of the phase diagram; the paper infers this from soft modes without converging that state, so a follow-up Hartree-Fock run with two electrons per Wigner cell would either confirm the endpoint or reveal a new competing phase.

Editorial analysis

A structured set of objections, weighed in public.

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

Referee Report

4 major / 4 minor

Summary. The paper studies Wigner crystallization in rhombohedral tetralayer graphene (R4G) at low electron densities and displacement fields, using projected Hartree–Fock (HF) with one electron per Wigner unit cell and time-dependent Hartree–Fock (TDHF) to test local stability. The central claim is that the HF ground state is a fully spin- and valley-polarized Wigner crystal with zero Chern number, and that in two regions of the (n_e, V_disp) phase diagram the crystal spontaneously breaks the spinless threefold rotational symmetry, forming a nematic Wigner crystal (nWC) that is locally stable. The authors attribute the nWC's energetic advantage, in most of the nWC phase, to a kinetic-energy gain from concentrating momentum occupation in one trigonal-warping Fermi-surface lobe. They propose experimental detection via STM imaging of the anisotropic charge distribution or via direction-dependent depinning in angle-resolved transport.

Significance. If the central claim holds, this is a notable prediction: an electronic, zero-field Wigner crystal that spontaneously breaks threefold rotational symmetry, with an explicit local-stability analysis. The manuscript has real strengths: the HF comparison is performed on a common underlying momentum mesh, which makes energy differences among Wigner-lattice candidates meaningful; both the charge-neutrality and average interaction schemes give similar nematic regions; convergence residuals and energy scales are reported; and the TDHF calculation is used to identify, rather than hide, instabilities toward two-electron-per-cell Wigner crystals, metallic Wigner crystals, Fermi liquids, or stripe phases. These are appropriate tools and give the paper a falsifiable, reproducible structure. However, the stability classification is undermined by internal inconsistencies and by a threshold that appears to be applied inconsistently to the reported eigenvalues.

major comments (4)
  1. [Appendix F.2 / Fig. 10] The stated TDHF stability criterion is inconsistent with the reported data. The text defines a stable solution as one with ξ(q) > 10^-4 meV at every local minimum with q ≠ 0, yet the two showcase 'stable' nWCs in Fig. 10 have ξ_min = −2.44×10^-7 meV at (0.25, 90) and −3.37×10^-6 meV at (0.3, 100), and many other nWC panels show negative eigenvalues. By the paper's own rule these states are unstable. If the intended condition was ξ(q) > −10^-4 meV (i.e., within numerical tolerance of zero), that must be stated explicitly and justified; as written, the stable regions in Fig. 1 are not reproducible from the stated criterion.
  2. [Sec. V / Fig. 10 / Appendix D2] The text calls the nWC at n_e = 0.4×10^12 cm^-2, V_disp = 80 meV a 'stable outlier', but Fig. 10 gives ξ_min = −0.168 meV for this point, more than two orders of magnitude below the stability threshold. This is not a numerical-tolerance effect. The contradiction between Sec. V, Appendix D2, and Fig. 10 means the phase-diagram labeling is not currently trustworthy, and the exception to the kinetic-energy mechanism is misidentified. Please correct the labeling and reconcile the text with the TDHF data.
  3. [Appendix F.2 / Eq. (F13)] The TDHF stability eigenvalues of interest are in the range 10^-6 to 10^-5 meV, while the HF self-consistency residual is ‖[H_HF,P]‖_max < 5×10^-4 meV. A stability threshold of 10^-4 meV is therefore of the same order as, or larger than, the numerical uncertainty of the HF state, and the reported negative eigenvalues may be mesh or roundoff artifacts. Without convergence tests of ξ(q) with respect to k_max, N_e, and the TDHF q-grid, the assignment of 'stable' versus 'unstable' is not quantitatively established. This is load-bearing because the central claim is defined by this threshold.
  4. [Appendix B / Eq. (4)] The projection onto the lowest conduction band is justified only by a ratio estimate (interaction scale below 1/5 of the interband gap). Since the TDHF eigenvalues used to define stability are at the 10^-4 meV scale, neglected remote-band and higher-band couplings could plausibly contribute at or above this level. Please provide a concrete check of the projection for representative (n_e, V_disp) points — for example, an unprojected HF calculation at small N_e or a two-band projected model — to show that the stability classification is robust to the projection truncation.
minor comments (4)
  1. [Eq. (7)] The operator T_r is introduced in Eq. (7) as a translation operator, but the text says the minimization runs over 'displacement of the WC lattice site'. Please clarify the domain of T_r (continuous real-space translations? lattice translations?) and its action on P, since this affects the gauge-invariance statement in Appendix C2.
  2. [Reference [73]] Reference [73] is cited as 'Roll et al. (2026)' with no title, journal, or preprint number. If this is an invited talk, please cite it as such with the appropriate details, or remove it and cite a published source for the paired Wigner crystal concept.
  3. [Sec. V / Appendix D2] The phrase 'stable outlier nWC' in Sec. V directly conflicts with Appendix D2's statement that this same outlying case 'is found to be unstable according to the TDHF calculation'. Harmonize the terminology so that 'stable' and 'unstable' are used consistently across the main text and appendices.
  4. [Fig. 10] Many of the ξ_min values plotted in Fig. 10 are below 10^-4 meV in magnitude and change sign erratically between adjacent parameter points. A color scale that saturates at ±10^-4 meV, or a separate panel showing only the sign of ξ_min at local minima, would make the stability classification much easier to assess than the current heat maps.

Circularity Check

0 steps flagged

No significant circularity: the nWC is a forward HF/TDHF result, not a fit; self-cited model and interaction formalism are externally benchmarked. One internal stability-label inconsistency is a correctness risk, not a circular reduction.

full rationale

The derivation chain is forward. The inputs are a single-particle continuum model (Ref. [75], independently benchmarked against DFT in Appendix A) plus a screened Coulomb interaction and two normal-ordering schemes (Refs. [76,77]). No parameter is fitted to reproduce the nematic WC. The HF calculation scans one-electron-per-cell Wigner-lattice ansätze, and the nWC emerges as the lowest-energy solution; the TDHF stability matrix is then computed from the same Hamiltonian as an independent local-stability check. No equation makes the predicted quantity equal to an input by construction: the kinetic-energy decomposition in Sec. V is an a posteriori energy accounting, not a defining constraint, and the nematicity diagnostics merely quantify the C3 breaking they report. The self-citations to the model and interaction-scheme papers are not load-bearing in a circular sense, because the model is benchmarked externally and the CN/average scheme comparison provides cross-checks; no uniqueness theorem is invoked. I therefore find no circular reduction. I do flag an internal numerical inconsistency that is a correctness risk, not circularity: Appendix F.2 defines stability by ξ(q)>1e−4 meV at q≠0, yet Fig. 10 labels as stable nWCs with negative minima (−2.44e−7 meV at (0.25,90), −3.37e−6 meV at (0.3,100)), and Sec. V calls (0.4,80) a 'stable outlier' although Fig. 10 gives ξ_min=−0.168 meV. These labels satisfy only a looser tolerance and need clarification, but they do not amount to fitting or renaming the target result, so the circularity score remains low.

Axiom & Free-Parameter Ledger

4 free parameters · 5 axioms · 0 invented entities

The prediction rests on: (a) single-particle model parameters benchmarked against DFT (App. A); (b) the double-gated screened Coulomb interaction with ε_r=5 and gate distance 50 nm (Eq. 3); (c) method-level assumptions — single-band projection (App. B), HF over one-electron-per-cell Wigner lattices with TDHF linear-stability as the ground-state criterion (Secs. III–IV), full spin-valley polarization at N_e=144, and the hand-set TDHF threshold of 10^-4 meV. Nothing is fitted to produce the nWC; energies are computed from stated Hamiltonians. The most fragile inputs are the screening parameters and the TDHF threshold, because the nWC-vs-C3 energy competition sits at the 0.1–0.5 meV scale.

free parameters (4)
  • Continuum-model parameters (v_F, v_3, v_4, t_1, t_2, V_ISP) = v_F=542.1 meV·nm; v_3=v_4=34.0 meV·nm; t_1=355.16 meV; t_2=−7 meV; V_ISP=16.65 meV
    Single-particle Hamiltonian (1)–(2) from Ref. [75], benchmarked against DFT in Appendix A. Not fitted to the nWC result, but the three-lobed trigonal-warping Fermi surface these parameters produce is the basis of the nematic mechanism (Sec. V).
  • Screening parameters (ε_r, ξ) = ε_r=5, gate-to-gate distance ξ=50 nm
    Double-gated screened Coulomb interaction V(q) in Eq. (3). Standard choices from the rhombohedral-graphene literature; the nWC-vs-C3 energy competition (0.1–0.5 meV/electron) is potentially sensitive to these values.
  • TDHF stability thresholds = 10^-4 meV for q≠0 minima; ξ(0) > 5‖[H_HF,P]‖_max
    Appendix F.2: a solution is deemed stable only if ξ(q)>10^-4 meV at every local minimum with q≠0. Several nWC points classified unstable have minima just below zero (−6.6×10^-5 to −10^-4 meV), so the claimed stable/unstable partition is numerically sensitive.
  • Momentum cutoff k_max = ~0.04–0.09 Å^-1 depending on (n_e, V_disp)
    Selected by a heuristic separation-of-energies criterion at N_e=36 and reused at N_e=144 (App. D1, Fig. 6(n)); phase boundaries could in principle shift at larger cutoff.
axioms (5)
  • domain assumption Projection onto the lowest conduction band is valid because the interaction scale e²√(n_e)/(4πε_rε_0) is below 1/5 of the interband gap throughout the phase diagram.
    Justifies the projected model Eq. (4); stated in Sec. II and Appendix B. Remote-band mixing would alter the kinetic-energy mechanism and the TDHF kernels.
  • domain assumption The double-gate screened Coulomb potential V(q)=πξ²V_ξ tanh(ξ|q|/2)/(ξ|q|/2) with ε_r=5, ξ=50 nm is the correct effective interaction.
    Eq. (3); standard for gated devices but not independently validated in this paper at these ultralow densities.
  • domain assumption The continuum model of Ref. [75] (parameters in App. A) accurately describes the low-energy bands of R4G.
    Supported by the DFT comparison in Appendix A, but remains a fitted low-energy model; the trigonal-warping lobe structure it produces is the basis of the nWC mechanism.
  • domain assumption HF with one electron per Wigner unit cell, plus TDHF linear stability, suffices to determine the ground state.
    Sections III–IV. The global competition with 2e/cell WCs, metallic WCs, Fermi liquids, and stripes is only probed through TDHF soft modes; these states are not directly converged.
  • domain assumption At N_e=144 all HF ground states are fully spin- and valley-polarized, assumed from the N_e=36 results.
    Sec. III and App. D2: initial states at N_e=144 are restricted to fully polarized configurations; a partially polarized state lower in energy would be missed.

pith-pipeline@v1.3.0-alltime-deepseek · 31899 in / 21600 out tokens · 180248 ms · 2026-08-01T19:14:39.389381+00:00 · methodology

0 comments
read the original abstract

Recent experiments have reported evidence for Wigner crystals (WCs) in rhombohedral graphene. Here, we investigate Wigner crystallization in rhombohedral tetralayer graphene using projected Hartree-Fock (HF) calculations and time-dependent Hartree-Fock (TDHF) calculations. We first perform HF calculations with one electron per Wigner unit cell, and find nematic WCs (nWCs) that spontaneously break the threefold rotational symmetry $C_3$ and $C_3$-invariant WCs. In particular, there are two nWC regions in the phase diagram: one larger region at large displacement fields and low electron densities, and another smaller region at intermediate fields and high densities. Both the nWCs and the $C_3$-invariant WCs are valley-polarized states with zero Chern number, and have positive indirect gaps in the HF band structure. We then perform TDHF calculations to further test the local stability of the WC states. We find that all $C_3$-invariant WCs and half of the nWCs are locally stable, while the remaining nWCs are unstable towards WCs with two electrons per unit cell or metallic states. The predicted stable nWC phase can be identified experimentally by scanning tunneling microscopy through its anisotropic charge distribution or by angle-resolved transport measurements via a direction-dependent depinning voltage.

Figures

Figures reproduced from arXiv: 2607.17027 by Ashwini Malviya, Ethan Angerhofer, Jiabin Yu, Yang Ge, Zhengguang Lu.

Figure 1
Figure 1. Figure 1: FIG. 1. HF phase diagram in the CN scheme for R [PITH_FULL_IMAGE:figures/full_fig_p003_1.png] view at source ↗
Figure 2
Figure 2. Figure 2: FIG. 2. (a),(c) Band structures and (b),(d) real-space charge den [PITH_FULL_IMAGE:figures/full_fig_p003_2.png] view at source ↗
Figure 3
Figure 3. Figure 3: FIG. 3. Lowest eigenvalue branch of the TDHF stability matrix for [PITH_FULL_IMAGE:figures/full_fig_p004_3.png] view at source ↗
Figure 4
Figure 4. Figure 4: FIG. 4. Momentum-space origin of the nWC at [PITH_FULL_IMAGE:figures/full_fig_p005_4.png] view at source ↗
Figure 5
Figure 5. Figure 5: FIG. 5. Comparison of the low-energy electronic band structures calculated using density functional theory (solid blue lines) and the continuum [PITH_FULL_IMAGE:figures/full_fig_p010_5.png] view at source ↗
Figure 6
Figure 6. Figure 6: FIG. 6. Phase diagram of the CN-scheme results for [PITH_FULL_IMAGE:figures/full_fig_p014_6.png] view at source ↗
Figure 7
Figure 7. Figure 7: FIG. 7. Further details of the phase diagram of the CN-scheme results for [PITH_FULL_IMAGE:figures/full_fig_p015_7.png] view at source ↗
Figure 8
Figure 8. Figure 8: FIG. 8. Phase diagram of the average-scheme results for [PITH_FULL_IMAGE:figures/full_fig_p016_8.png] view at source ↗
Figure 9
Figure 9. Figure 9: FIG. 9. Lowest-energy TDHF neutral-mode dispersion in meV at each [PITH_FULL_IMAGE:figures/full_fig_p021_9.png] view at source ↗
Figure 10
Figure 10. Figure 10: FIG. 10. Smallest TDHF stability eigenvalue in meV at each [PITH_FULL_IMAGE:figures/full_fig_p022_10.png] view at source ↗

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Reference graph

Works this paper leans on

106 extracted references · 22 linked inside Pith

  1. [1]

    Wigner, On the interaction of electrons in metals, Phys

    E. Wigner, On the interaction of electrons in metals, Phys. Rev. 46, 1002 (1934)

  2. [2]

    Faraday Soc.34, 678 (1938)

    E.Wigner,Effectsoftheelectroninteractionontheenergylevels of electrons in metals, Trans. Faraday Soc.34, 678 (1938)

  3. [3]

    Shayegan, Wigner crystals in flat band 2d electron systems, Nat

    M. Shayegan, Wigner crystals in flat band 2d electron systems, Nat. Rev. Phys.4, 212 (2022)

  4. [4]

    C. C. Grimes and G. Adams, Evidence for a liquid-to-crystal phase transition in a classical, two-dimensional sheet of elec- trons, Phys. Rev. Lett.42, 795 (1979)

  5. [5]

    Monarkha andK

    Y. Monarkha andK. Kono,Two-Dimensional CoulombLiquids and Solids, Springer Series in Solid-State Sciences, Vol. 142 (Springer, Berlin, Heidelberg, 2004)

  6. [6]

    Pollanen, High-frequency nonlinear conductivity of a Wigner crystal, arXiv:2607.06846 [cond-mat.mes-hall] (2026)

    A.J.Schleusner,M.T.Elewa,N.R.Beysengulov,C.A.Mikolas, and J. Pollanen, High-frequency nonlinear conductivity of a Wigner crystal, arXiv:2607.06846 [cond-mat.mes-hall] (2026)

  7. [7]

    Etienne, Observation of a magnetically induced Wigner solid, Phys

    E.Y.Andrei,G.Deville,D.C.Glattli,F.I.B.Williams,E.Paris, and B. Etienne, Observation of a magnetically induced Wigner solid, Phys. Rev. Lett.60, 2765 (1988)

  8. [8]

    V.J.Goldman,M.Santos,M.Shayegan,andJ.E.Cunningham, Evidence for two-dimensional quantum Wigner crystal, Phys. Rev. Lett.65, 2189 (1990)

  9. [9]

    H. W. Jiang, R. L. Willett, H. L. Stormer, D. C. Tsui, L. N. Pfeiffer, and K. W. West, Quantum liquid versus electron solid aroundν= 1/5Landau-level filling, Phys. Rev. Lett.65, 633 (1990)

  10. [10]

    Das Sarma and A

    S. Das Sarma and A. Pinczuk, eds.,Perspectives in Quantum HallEffects: NovelQuantumLiquidsinLow-DimensionalSemi- conductor Structures(Wiley, 1997)

  11. [11]

    Y.Liu, S.Hasdemir, L.N.Pfeiffer, K.W.West, K.W.Baldwin, andM.Shayegan,ObservationofananisotropicWignercrystal, Phys. Rev. Lett.117, 106802 (2016), arXiv:1601.07135

  12. [12]

    M. S. Hossain, M. K. Ma, K. A. Villegas-Rosales, Y. J. Chung, L. N. Pfeiffer, K. W. West, K. W. Baldwin, and M. Shayegan, Anisotropic, two-dimensional, disordered Wigner solid, Phys. Rev. Lett.129, 036601 (2022), arXiv:2207.06618

  13. [13]

    B. A. Piot, Z. Jiang, C. R. Dean, L. W. Engel, G. Gervais, L. N. Pfeiffer,andK.W.West,Wignercrystallizationinaquasi-three- dimensional electronic system, Nat. Phys.4, 936 (2008)

  14. [14]

    J. Yoon, C. C. Li, D. Shahar, D. C. Tsui, and M. Shayegan, Wigner crystallization and metal-insulator transition of two- dimensional holes in GaAs atB= 0, Phys. Rev. Lett.82, 1744 (1999)

  15. [15]

    Falson, I

    J. Falson, I. Sodemann, B. Skinner, D. Tabrea, Y. Kozuka, A. Tsukazaki, M. Kawasaki, K. von Klitzing, and J. H. Smet, 6 Competing correlated states around the zero-field Wigner crys- tallizationtransitionofelectronsintwodimensions,Nat.Mater. 21, 311 (2022)

  16. [16]

    T.Smoleński,P.E.Dolgirev,C.Kuhlenkamp,A.Popert,Y.Shi- mazaki, P. Back, X. Lu, M. Kroner, K. Watanabe, T. Taniguchi, I. Esterlis, E. Demler, and A. Imamoğlu, Signatures of Wigner Crystal of electrons in a monolayer semiconductor, Nature595, 53 (2021)

  17. [17]

    Y. Zhou, J. Sung, E. Brutschea, I. Esterlis, Y. Wang, G. Scuri, R.J.Gelly,H.Heo,T.Taniguchi,K.Watanabe,G.Zaránd,M.D. Lukin,P.Kim,E.Demler,andH.Park,BilayerWignerCrystals inatransitionmetaldichalcogenideheterostructure,Nature595, 48 (2021)

  18. [18]

    Xiang, H

    Z. Xiang, H. Li, J. Xiao, M. H. Naik, Z. Ge, Z. He, S. Chen, J. Nie, S. Li, Y. Jiang, R. Sailus, R. Banerjee, T. Taniguchi, K. Watanabe, S. Tongay, S. G. Louie, M. F. Crommie, and F. Wang, Quantum melting of a disordered Wigner solid, arXiv:2402.05456 [cond-mat.str-el] (2024)

  19. [19]

    L. Wang, F. Menzel, F. Pichler, P. Knüppel, K. Watanabe, T. Taniguchi, M. Knap, and T. Smoleński, Spectroscopy of Wigner crystal polarons in an atomically thin semiconductor, arXiv:2512.16552 [cond-mat.mes-hall] (2025)

  20. [20]

    Zhang, L

    L. Zhang, L. Gu, H. S. Adlong, A. Christianen, E. Dizer, R. Ni, R. Ma, S. Park, H. Jang, T. Taniguchi, K. Watanabe, I. Esterlis, R. Schmidt, A. Imamoglu, and Y. Zhou, Wigner polarons re- veal Wigner crystal dynamics in a monolayer semiconductor, arXiv:2512.16631 [cond-mat.mes-hall] (2025)

  21. [21]

    S.Munyan,S.Ahadi,B.Guo,A.Rashidi,andS.Stemmer,Evi- dence of zero-field Wigner solids in ultrathin films of cadmium arsenide, Phys. Rev. X14, 041037 (2024), arXiv:2411.07230

  22. [22]

    Y.-C. Tsui, M. He, Y. Hu, E. Lake, T. Wang, K. Watanabe, T. Taniguchi, M. P. Zaletel, and A. Yazdani, Direct observation of a magnetic-field-induced Wigner crystal, Nature628, 287 (2024)

  23. [23]

    Xiang, J

    Z. Xiang, J. Xiao, H. Li, W. Kim, T. Wang, Z. Dong, T. Taniguchi, K. Watanabe, M. P. Zaletel, S. G. Louie, M. F. Crommie, and F. Wang, Imaging interacting two- dimensional anisotropic electrons, arXiv:2605.12761 [cond- mat.str-el] (2026)

  24. [24]

    A.M.Ettouhami,C.B.Doiron,F.D.Klironomos,R.Côté,and A. T. Dorsey, Anisotropic states of two-dimensional electrons in high magnetic fields, Phys. Rev. Lett.96, 196802 (2006)

  25. [25]

    E. Berg, M. S. Rudner, and S. A. Kivelson, Electronic liquid- crystallinephasesinaspin-orbit-coupledtwo-dimensionalelec- tron gas, Phys. Rev. B85, 035116 (2012)

  26. [26]

    G.Zhu,J.Koch,andI.Martin,Nematicquantumliquidcrystals ofbosonsinfrustratedlattices,Phys.Rev.B93,144508(2016)

  27. [27]

    Kerelsky, C

    A. Kerelsky, C. Rubio-Verdú, L. Xian, D. M. Kennes, D. Hal- bertal, N. Finney, L. Song, S. Turkel, L. Wang, K. Watanabe, T. Taniguchi, J. Hone, C. Dean, D. N. Basov, A. Rubio, and A. N. Pasupathy, Moiréless correlations in ABCA graphene, Proc. Natl. Acad. Sci. U.S.A.118, e2017366118 (2021)

  28. [28]

    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)

  29. [29]

    Taniguchi, J

    K.Liu,J.Zheng,Y.Sha,B.Lyu,F.Li,Y.Park,Y.Ren,K.Watan- abe, T. Taniguchi, J. Jia, Z. Shi, and G. Chen, Spontaneous broken-symmetry insulator and metals in tetralayer rhombohe- dral graphene, Nat. Nanotechnol.19, 188 (2024)

  30. [30]

    Taniguchi, L

    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, Nature623, 41 (2023)

  31. [31]

    H. Zhou, T. Xie, T. Taniguchi, K. Watanabe, and A. F. Young, Superconductivity in rhombohedral trilayer graphene, Nature 598, 434 (2021)

  32. [32]

    T. Han, Z. Lu, Z. Hadjri, L. Shi, Z. Wu, W. Xu, Y. Yao, A. A. Cotten,O.SharifiSedeh,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ühl, L.Fu,andL.Ju,Signaturesofchiral superconductivity in rhombohedral graphene, Nature643, 654 (2025)

  33. [33]

    R. Q. Nguyen, H.-T. Wu, E. Morissette, N. J. Zhang, P. Qin, K. Watanabe, T. Taniguchi, A. W. Hui, D. E. Feldman, and J. I. A. Li, A hierarchy of superconductivity and topo- logical charge-density-wave states in rhombohedral graphene, arXiv:2507.22026 [cond-mat.mes-hall] (2025)

  34. [34]

    J. Seo, A. A. Cotten, M. Xu, O. S. Sedeh, H. Weldeyesus, T. Han, Z. Lu, Z. Wu, S. Ye, W. Xu, J. Yang, E. Aitken, P. P. Liong, Z. Hadjri, R. Gazizulin, K. Watanabe, T. Taniguchi, M. Li, D. M. Zumbühl, and L. Ju, Family of unconventional superconductivities in crystalline Graphene, arXiv:2509.03295 [cond-mat.mes-hall] (2025)

  35. [35]

    Mak, and X

    J.Deng,J.Xie,H.Li,T.Taniguchi,K.Watanabe,J.Shan,K.F. Mak, and X. Liu, Superconductivity and ferroelectric orbital magnetism in semimetallic rhombohedral hexalayer graphene, arXiv:2508.15909 [cond-mat.supr-con] (2025)

  36. [36]

    J. Yang, O. Sharifi Sedeh, C. Yoon, S. Ye, H. Weldeye- sus, A. Cotten, T. Han, Z. Lu, Z. Hadjri, J. Seo, L. Shi, E. Aitken, P. P. Liong, Z. Wu, M. Xu, C. Scheller, M. Zheng, R. Gazizulin, K. Watanabe, T. Taniguchi, D. Laroche, M. Li, F. Zhang, D. M. Zumbühl, and L. Ju, Magnetic-field-enhanced graphene superconductivity with record Pauli-limit violation, a...

  37. [37]

    C.Zheng,C.Li,C.Zhang,K.Watanabe,T.Taniguchi,H.Yang, D. Guan, L. Liu, S. Wang, Y. Li, H. Zheng, C. Liu, J. Jia, Z. Shi, G. Chen, T. Li, and X. Liu, Multiple superconducting phasesinrhombohedralheptalayergraphene,arXiv:2605.17890 [cond-mat.mes-hall] (2026)

  38. [38]

    S. S. Kalantre, B. H. Alexander, J. May-Mann, J. Herzog- Arbeitman, M. Hocking, Q. Cao, K. Watanabe, T. Taniguchi, D. Goldhaber-Gordon, A. J. Mannix, T. Devakul, Y. H. Kwan, D. E. Parker, and A. Sharpe, Fermiology and the candidate chiral superconductor in rhombohedral tetralayer graphene, arXiv:2606.05356 [cond-mat.str-el] (2026)

  39. [39]

    Kumar, D

    M. Kumar, D. Waleffe, A. Okounkova, R. Tejani, V. T. Phong, K. Watanabe, T. Taniguchi, C. Lewandowski, J. Folk, and M. Yankowitz, Superconductivity from dual-surface carriers in rhombohedral graphene, Nat. Phys.22, 862 (2026)

  40. [40]

    R. Q. Nguyen, P. Qin, H.-T. Wu, S. Mishra, T. Wolf, J. Roll, E. Morissette, N. J. Zhang, S. Alkidim, K. Watan- abe, T. Taniguchi, A. W. Hui, D. E. Feldman, A. MacDonald, and J. I. A. Li, Coexisting charge density wave and supercon- ducting order in quantizing magnetic fields, arXiv:2607.05039 [cond-mat.mes-hall] (2026)

  41. [41]

    Nanotechnol.19, 181 (2024)

    T.Han,Z.Lu,G.Scuri,J.Sung,J.Wang,T.Han,K.Watanabe, T.Taniguchi,H.Park,andL.Ju,CorrelatedinsulatorandChern insulators in pentalayer rhombohedral-stacked graphene, Nat. Nanotechnol.19, 181 (2024)

  42. [42]

    T. Han, Z. Lu, Y. Yao, J. Yang, J. Seo, C. Yoon, K. Watan- abe, T. Taniguchi, L. Fu, F. Zhang, and L. Ju, Large quantum anomalous Hall effect in spin-orbit proximitized rhombohedral graphene, Science384, 647 (2024)

  43. [43]

    Y. Sha, J. Zheng, K. Liu, H. Du, K. Watanabe, T. Taniguchi, J. Jia, Z. Shi, R. Zhong, and G. Chen, Observation of a Chern insulator in crystalline abca-tetralayer graphene with spin-orbit coupling, Science384, 414 (2024)

  44. [44]

    H. Li, Z. Geng, J. Seo, C. Xu, Y. Jiang, S. Ye, Z. Hua, 7 J. Xie, L. Min, K. Watanabe, T. Taniguchi, K. Yasuda, X. Liu, L. Ju, J. Shan, and K. F. Mak, Competing Chern states revealed by quasiparticle charging in Moiré Rhombohedral Graphene, arXiv:2607.08710 [cond-mat.mes-hall] (2026)

  45. [45]

    J. Dong, T. Wang, T. Wang, T. Soejima, M. P. Zaletel, A. Vish- wanath, and D. E. Parker, Anomalous Hall crystals in rhombo- hedral multilayer graphene. I. interaction-driven Chern bands and fractional quantum Hall states at zero magnetic field, Phys. Rev. Lett.133, 206503 (2024)

  46. [46]

    Soejima, J

    T. Soejima, J. Dong, T. Wang, T. Wang, M. P. Zaletel, A. Vish- wanath,andD.E.Parker,AnomalousHallcrystalsinrhombohe- dral multilayer graphene. II. general mechanism and a minimal model, Phys. Rev. B110, 205124 (2024)

  47. [47]

    Z.Dong,A.S.Patri,andT.Senthil,StabilityofanomalousHall crystalsinmultilayerrhombohedralgraphene,Phys.Rev.B110, 205130 (2024)

  48. [48]

    Berrytrashcan

    B.A.BernevigandY.H.Kwan,“Berrytrashcan”modelofinter- acting electrons in rhombohedral Graphene, arXiv:2503.09692 (2025)

  49. [49]

    W.MiaoandC.Li,Variouselectroniccrystalphasesinrhombo- hedral graphene multilayers, Phys. Rev. B113, 155136 (2026)

  50. [50]

    Desrochers and A

    F. Desrochers and A. Vishwanath, Energetics of frac- tional Anomalous Hall crystals in Rhombohedral Graphene, arXiv:2607.08822 [cond-mat.str-el] (2026)

  51. [51]

    B.Zhou,H.Yang,andY.-H.Zhang,Fractionalquantumanoma- lous Hall effect in rhombohedral multilayer graphene in the moiréless limit, Phys. Rev. Lett.133, 206504 (2024)

  52. [52]

    Huang, X

    K. Huang, X. Li, S. D. Sarma, and F. Zhang, Self-consistent theory of fractional quantum anomalous Hall states in rhombo- hedral graphene, Phys. Rev. B110, 115146 (2024)

  53. [53]

    K.Huang,S.DasSarma,andX.Li,Fractionalquantumanoma- lous Hall effect in rhombohedral multilayer graphene with a strong displacement field, Phys. Rev. B111, 075130 (2025)

  54. [54]

    Desrochers, M

    F. Desrochers, M. R. Hirsbrunner, J. Huxford, A. S. Patri, T. Senthil, and Y. B. Kim, Elastic response and instabilities of Anomalous Hall crystals, Phys. Rev. Lett.136, 166503 (2026)

  55. [55]

    Zeng and A

    Y. Zeng and A. J. Millis, Berry phase dynamics of sliding elec- tron crystals, Phys. Rev. X15, 031059 (2025)

  56. [56]

    Uchida, T

    T. Uchida, T. Kawakami, and M. Koshino, Non-Abelian Chern band in rhombohedral graphene multilayers, Phys. Rev. Lett. 136, 156602 (2026)

  57. [57]

    Yoon and F

    C. Yoon and F. Zhang, Majorana crystal in rhombohedral graphene, arXiv:2603.16828 [cond-mat.mes-hall] (2026)

  58. [58]

    V. T. Phong, E. Prada, P. San-Jose, F. Guinea, and E. J. Mele, Valley valves at domain walls in symmetry-broken rhombohe- dral graphene, arXiv:2606.14878 [cond-mat.mes-hall] (2026)

  59. [59]

    T. Han, J. P. Butler, S. Ye, Z. Hua, S. Dutta, Z. Hadjri, Z. Wu, J. Yang, J. Seo, P. Pattanakanvijit, E. Aitken, K. Watanabe, T.Taniguchi, P.Xiong, E.Zeldov, Z.Lu, R.Ashoori,and L.Ju, EvidenceofmetallicWignercrystalinrhombohedralGraphene, arXiv:2604.00113 [cond-mat.mes-hall] (2026)

  60. [60]

    Z. Zhou, K. Krötzsch, R. Ayache, Y. Li, S. Joy, K. Watanabe, T.Taniguchi, M.Heiblum, P.Roulleau,andM.Banerjee,Com- peting orders driven by Wigner crystal phase in rhombohedral graphene, arXiv:2607.15014 [cond-mat.mes-hall] (2026)

  61. [61]

    J. Dong, T. Soejima, D. E. Parker, and A. Vishwanath, Crys- tals caught doping: Metallic Wigner crystals in rhombohedral Graphene, arXiv:2604.00114 [cond-mat.str-el] (2026)

  62. [62]

    J. Feng, Z. Han, M. P. Zaletel, and Z. Dong, Self-doped crystal from preempted band-inversion transitions, arXiv:2604.09820 [cond-mat.str-el] (2026)

  63. [63]

    J. Dong, O. E. Sommer, T. Soejima, D. E. Parker, and A. Vish- wanath, Phonons in electron crystals with Berry curvature, arXiv:2503.16390 [cond-mat.str-el] (2025)

  64. [64]

    MacDonald,Spinandorbitalmetallicmagnetisminrhombohe- dral trilayer graphene, Phys

    C.Huang,T.M.R.Wolf,W.Qin,N.Wei,I.V.Blinov,andA.H. MacDonald,Spinandorbitalmetallicmagnetisminrhombohe- dral trilayer graphene, Phys. Rev. B107, L121405 (2023)

  65. [65]

    Parra-Martínez, A

    G. Parra-Martínez, A. Jimeno-Pozo, J. A. Silva-Guillén, and F. Guinea, Nematic and partially polarized phases in rhombo- hedral graphene with varying number of layers: An extensive Hartree–Fock study, Adv. Mater. , e14008 (2025)

  66. [66]

    Banerjee, H

    S. Banerjee, H. D. Scammell, and M. S. Scheurer, Correlations between superconducting and resistive anisotropies, Phys. Rev. B113, 184519 (2026)

  67. [67]

    Morissette, P

    E. Morissette, P. Qin, K. Watanabe, T. Taniguchi, and J. I. A. Li, Evidence of momentum-space condensation in rhombo- hedral hexalayer graphene, arXiv:2503.09954 [cond-mat.mes- hall] (2025)

  68. [68]

    Qin, H.-T

    P. Qin, H.-T. Wu, R. Q. Nguyen, E. Morissette, N. J. Zhang, K. Watanabe, T. Taniguchi, and J. I. A. Li, Stripe order in the metallicandsuperconductingphasesofrhombohedralhexalayer graphene, arXiv:2504.05129 [cond-mat.mes-hall] (2025)

  69. [69]

    Moulopoulos and N

    K. Moulopoulos and N. W. Ashcroft, New low-density phase of interacting electrons: The paired electron crystal, Phys. Rev. Lett.69, 2555 (1992)

  70. [70]

    Moulopoulos and N

    K. Moulopoulos and N. W. Ashcroft, Many-body theory of paired electron crystals, Phys. Rev. B48, 11646 (1993)

  71. [71]

    Taut, On the paired Wigner crystal, Solid State Commun

    M. Taut, On the paired Wigner crystal, Solid State Commun. 89, 189 (1994)

  72. [72]

    Zverevich, A

    D. Zverevich, A. Levchenko, and I. Esterlis, Spin-triplet paired Wignercrystalstabilizedbyquantumgeometry,Phys.Rev.Lett. 137, 036506 (2026)

  73. [73]

    (2026), Allan H

    Roll et al. (2026), Allan H. MacDonald’s Invited talk at New Developments in 2D Correlated Materials, Fine Theoretical Physics Institute

  74. [74]

    G. A. Csáthy, D. C. Tsui, L. N. Pfeiffer, and K. W. West, Asta- bility and negative differential resistance of the Wigner solid, Phys. Rev. Lett.98, 066805 (2007)

  75. [75]

    Herzog-Arbeitman, Y

    J. Herzog-Arbeitman, Y. Wang, J. Liu, P. M. Tam, Z. Qi, Y. Jia, D. K. Efetov, O. Vafek, N. Regnault, H. Weng, Q. Wu, B. A. Bernevig, and J. Yu, Moiré fractional Chern insulators. II. first- principles calculations and continuum models of rhombohedral graphene superlattices, Phys. Rev. B109, 205122 (2024)

  76. [76]

    Y. H. Kwan, J. Yu, J. Herzog-Arbeitman, D. K. Efetov, N. Reg- nault, and B. A. Bernevig, Moiré fractional Chern insulators. III.Hartree-Fockphasediagram,magicangleregimeforChern insulator states, role of moiré potential, and Goldstone gaps in rhombohedralgraphenesuperlattices,Phys.Rev.B112,075109 (2025)

  77. [77]

    B. A. Bernevig, Z.-D. Song, N. Regnault, and B. Lian, Twisted bilayergraphene.III.interactinghamiltonianandexactsymme- tries, Phys. Rev. B103, 205413 (2021)

  78. [78]

    E.D’HokerandJ.Kaidi,Lecturesonmodularformsandstrings, arXiv:2208.07242 [hep-th] (2022)

  79. [79]

    D. J. Rowe,Nuclear collective motion: models and theory (World Scientific, 2010)

  80. [80]

    Ring and P

    P. Ring and P. Schuck,The Nuclear Many-Body Problem (Springer, Berlin, 1980)

Showing first 80 references.