Pith. sign in

REVIEW 4 major objections 5 minor 64 references

A Wigner crystal with a hole at its zone center can beat the filled crystal near the melting transition.

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 →

Challenge this review Re-run · record.json Download PDF Read on arXiv ↗

T0 review · deepseek-v4-flash

2026-08-04 01:07 UTC pith:T7CQKLKB

load-bearing objection Serious VMC study with a genuinely new dressed-hole ansatz; the self-doping window is a real variational result but rests on a 0.5–1% energy difference with no finite-size scaling, so treat it as promising rather than established. the 4 major comments →

arxiv 2608.00170 v1 pith:T7CQKLKB submitted 2026-07-31 cond-mat.str-el cond-mat.mes-hall

Shape of Wigner Crystals and Hole Self-Doping in a Mexican-Hat Dispersion

Minho Luke Kim , Xiao-Gang Wen This is my paper
classification cond-mat.str-el cond-mat.mes-hall
keywords Wigner crystalMexican-hat dispersionhole self-dopingvariational Monte CarlopolaronJastrow factorrhombohedral multilayer graphenemetallic Wigner crystal
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 studies how a Wigner crystal forms from the ring-like Fermi surface of a Mexican-hat dispersion, where kinetic energy rises as c₂k² + c₄k⁴ with a negative c₂. It shows that the lowest-energy crystal orbitals must deplete electrons near k=0 to lower kinetic energy. Even after such orbital shaping, a doped crystal with a single hole at the Brillouin-zone center can be lower in energy than the undoped crystal near the crystallization transition, provided the hole is dressed by a polaron factor that lets neighboring electrons relax toward the vacancy. This supports the possibility of a metallic Wigner crystal whose conductance is carried by holes, and the paper estimates the hole dispersion with hopping t₁ ≈ 0.13–0.18 E_c, t₂ ≈ 0, and an effective mass of about 0.07 m_e.

Core claim

The central claim is that, for a Mexican-hat dispersion with parameters calibrated to rhombohedral multilayer graphene, the ground state near the Wigner-crystal transition is a self-doped charge-density wave with a hole at Γ, not the commensurate crystal. The paper constructs this state by taking Bloch orbitals of an auxiliary periodic potential (the CDW ansatz), removing the Γ orbital, and multiplying by a vacancy-conditioned Gaussian dressing factor that enhances electron density around the missing site (the pCDW ansatz). Variational Monte Carlo energies show pCDWΓ lies below the commensurate CDW for electron densities around n_e ≈ 0.3–0.6×10¹² cm⁻², with a gain of order 0.5–1% of the Coul

What carries the argument

The load-bearing object is the polaron-dressed CDW wavefunction Ψ_pCDW_k = J Σ_R e^{-ik·R} F_R D_R, where D_R is the Slater determinant with the Wannier orbital at site R removed (a pinned vacancy) and F_R = exp[A_p Σ_i exp(-|r_i-R|²/2ξ_p²)] is a Gaussian dressing that attracts electrons toward the vacancy, screening it. This electron–vacancy correlation, added on top of the standard Jastrow factor J, is what tips the energy balance toward self-doping. The CDW orbitals themselves come from an auxiliary periodic single-particle Hamiltonian H_aux = c₂k² + c₄k⁴ − 2V₀ Σ cos(G_j·r), whose lowest-band Bloch states form the crystal.

Load-bearing premise

The polaron dressing that makes the doped crystal win contributes only about half a percent to one percent of the Coulomb energy, while quantum geometry, trigonal warping, and band-topological effects — which the paper omits — can shift competing-state energies by the same few-percent scale, so a few-percent correction could close the self-doping window.

What would settle it

Compute the CDW versus pCDWΓ energy difference with the same variational Monte Carlo setup but with a trigonal-warping term of magnitude ~1% of E_c added; if the ordering reverses, the self-doping claim fails for real rhombohedral graphene. Alternatively, a transport or scanning-tunneling experiment at n_e ≈ 0.4×10¹² cm⁻² that resolves the Fermi surface and finds no hole pocket at Γ would falsify the prediction.

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

If this is right

  • Near the crystallization transition in rhombohedral multilayer graphene, a metallic Wigner crystal with hole carriers at Γ is a plausible ground state, not merely an excited state.
  • The hole's effective mass, a few times lighter than the bare band mass, explains hole conduction with lighter mass than band calculations.
  • Energy differences between holes at Γ, K, and M yield a tight-binding model with t₁ ≈ 0.13–0.18 E_c and t₂ ≈ 0, implying nearly nearest-neighbor-only hopping.
  • The polaron dressing mechanism — density relaxation around a vacancy — is generic and could apply to other Wigner-crystal systems with non-quadratic dispersions.
  • To test the self-doping window, experiments should look for a hole pocket at Γ coexisting with crystalline order near n_e ≈ 0.4×10¹² cm⁻².

Where Pith is reading between the lines

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

  • The paper tests only a single vacancy; if the finite-density generalization preserves the polaron gain, a finite concentration of holes may form, making the metallic WC a true phase rather than a near-degenerate point.
  • Quantum geometry and trigonal warping, set aside here, could either close or widen the self-doping window; the paper's own estimate puts their scale at a few percent of E_c, the same order as the polaron gain, so a systematic treatment of those terms is the natural next step.
  • A direct diagnostic would be to compute the static structure factor or electron density correlation of the trial pCDW state to confirm metallicity, a direction the paper notes but does not pursue.
  • The auxiliary-periodic-potential method for shaping Wannier orbitals could be exported to other flat or ring-dispersion systems where Gaussian orbitals fail.

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 / 5 minor

Summary. Using variational Monte Carlo with Slater-Jastrow wavefunctions, this paper studies Wigner crystallization in a Mexican-hat dispersion epsilon_k = c_2 k^2 + c_4 k^4. The authors introduce annular (aWC) orbitals with a high-pass momentum filter and CDW orbitals obtained from the lowest band of an auxiliary periodic potential, all supplemented by a short- and long-range Jastrow factor. They calibrate c_2 = -190 meV nm^2, c_4 = 1000 meV nm^4, and epsilon = 10 so that the optimized QFL and commensurate crystal energies cross at densities consistent with the experimental Lifshitz and Wigner-crystal transitions. The central claim is that a one-hole doped crystal, hCDW_Gamma, and especially a mobile polaron-dressed variant pCDW_Gamma (Eq. 10), lies below the commensurate CDW for n_e ~ 0.3-0.6 x 10^12 cm^-2, with the polaron dressing contributing 0.5-1% of the Coulomb energy per particle. From energy differences at Gamma, K, and M, the authors estimate hole hoppings t_1 ~ 0.13-0.18 E_c, t_2 ~ 0, and a renormalized hole mass ~ 0.07 m_e.

Significance. The variational machinery is a genuine step beyond Gaussian Wigner-crystal ansatze: the aWC/CDW comparison supplies an internal optimization check, and the paper validates the Ewald/Coulomb and kinetic-energy samplers against the Madelung energy and Hartree-Fock limits (App. A.4-A.5). The pCDW construction is an interesting and novel way to include vacancy-conditioned density relaxation, and the extracted hole mass is a falsifiable, order-of-magnitude prediction. If the energy ordering survives finite-size and parameter-robustness checks, the result would provide a concrete self-doping mechanism for metallic Wigner crystals in rhombohedral graphene. However, the central effect is a sub-percent energy difference at a single system size, and the paper itself identifies omitted band effects at the few-percent scale; these caveats currently make the headline claim suggestive rather than established.

major comments (4)
  1. [App. A.3, Fig. 3, App. B.2] The central crossing pCDW_Gamma vs CDW rests on an energy difference of order 0.5-1% of E_c per particle (App. B.2, Fig. 5). All VMC runs use 64 electrons (CDW) and 63 electrons (hCDW/pCDW) on the torus; no N=127/256 check or finite-size extrapolation is reported. The Ewald potential and the long-range Jastrow u_LR(q) explicitly depend on the torus dimensions (Eqs. A9-A13), so the single-hole chemical potential, obtained as an O(N) total-energy difference, can be contaminated by O(1/N) box effects of the same magnitude as the effect. Please report Delta(n_e) = E_CDW - E_pCDW_Gamma as a function of N (e.g., 64/63, 128/127, 256/255) with statistical error bars, and provide a finite-size extrapolation at n_e = 0.4 x 10^12 cm^-2 before drawing the phase conclusion.
  2. [App. A.5] The band parameters and dielectric constant are not independent inputs: c_2, c_4, and epsilon are chosen semi-iteratively so that the optimized variational energies of QFL and CDW/aWC reproduce the experimental Lifshitz density (~0.7 x 10^12 cm^-2) and crystallization density (~0.4-0.5 x 10^12 cm^-2). Thus the phrase 'reproduce the two transitions' in the abstract is a calibration statement, and the location of the hole-self-doping window is not an independent prediction. This is acceptable if framed as constrained modeling, but the paper should explicitly state the calibrated nature in the abstract and conclusions, and should test robustness of the pCDW-CDW ordering under variations of c_2, c_4, epsilon within experimental uncertainty (e.g., +/-20% in c_2 and c_4).
  3. [Eq. (10), Discussion] The doped state pCDW_Gamma contains exactly one vacancy (63 electrons on a 64-site torus; Eq. 10, Eq. B1). A negative single-hole chemical potential is necessary for an infinitesimal-doping instability, but the paper concludes that 'a hole-self-doping instability' occurs over 'a substantial density window' without evaluating the energy as a function of hole density or including hole-hole interactions (the Discussion explicitly sets them aside). The finite-density conclusion requires the energy curvature with respect to the number of holes, which is not provided. Please add a multi-hole calculation at a few densities (or a well-justified estimate of the hole-hole repulsion) to support the claimed instability window.
  4. [Discussion] The Discussion states that quantum geometry, trigonal warping, and band-topological effects 'will play a role in the fine energetics of the various competing states as they are within few percent of Coulomb energy scale.' Since the polaron gain that makes pCDW_Gamma beat CDW is only 0.5-1% of E_c (App. B.2), the omitted terms are of the same order as the effect driving the headline result. I am not arguing that these effects are necessarily destructive, but the claim as it stands has no quantified error budget. Please estimate the leading geometric/warping correction in the n_e ~ 0.4 x 10^12 cm^-2 window (e.g., by adding a Berry-curvature or band-anisotropy correction to the auxiliary Hamiltonian) or clearly label the result as a statement valid only within the Mexican-hat model.
minor comments (5)
  1. [Figs. 2 and 3] The y-axis shows (E - E_offset)/(sqrt(n_e) e^2/epsilon), where E_offset is a density-dependent polynomial. State explicitly that E_offset is a common shift applied to all states at a given density, and give the fit coefficients or uncertainties; otherwise the reader cannot reconstruct absolute energy differences.
  2. [Fig. 5] Error bars are displayed only for the hCDW_K curve; the caption says other states have similar errors. A table with numerical values of the energy differences and their statistical errors at n_e = 0.4 x 10^12 cm^-2 would be much more informative for assessing the 0.5-1% E_c polaron gain.
  3. [Ref. [45]] Reference [45] currently appears as a footnote reading 'References [2, 10, 37, 53, 54] are relevant in the supplemental material.' This looks like a placeholder and should be replaced by a proper citation to the supplemental material.
  4. [Eq. (2)] Define psi_gWC(k) explicitly and specify the normalization of psi_aWC. The long algebraic tail is controlled by xi, but the normalization of the orbital is not given.
  5. [Inset of Fig. 2] Z_Gamma = |u_{G=0}(Gamma)|^2 is used as a diagnostic for reciprocal-vector mixing. Define u_G(k) in the main text or point explicitly to App. A.1, since the main text otherwise leaves this notation undefined.

Circularity Check

0 steps flagged

No significant circularity: the hole-self-doping comparison is a variational output, not a fitted target; acknowledged calibration and finite-size/model limitations are correctness risks, not tautologies.

full rationale

The paper's derivation chain is self-contained for the claimed new result. The band parameters c2, c4, and the dielectric constant are explicitly calibrated to reproduce the Lifshitz and WC transition densities (Terminology section and Appendix A.5: the parameters are chosen so that the crystal state wins over the QFL at approximately 0.5e12 cm^-2), and the abstract correctly labels this as calibration rather than prediction. The central claim—that the pCDW(Gamma) one-vacancy state falls below the commensurate CDW near n_e ~0.4e12 cm^-2—is an independent variational Monte Carlo comparison: no parameter was fitted to make pCDW win, and the energy difference is not a restatement of the inputs. The polaron dressing is defined explicitly in Eqs. (10)-(12) and computed in Appendix B; citing Refs. [43,44] for similar VMC techniques is not load-bearing because the needed formulas and sampling procedures are given in the paper. The finite-size limitation (64 vs 63 electrons, no scaling) and the acknowledged neglect of quantum geometry/trigonal warping at the few-percent-of-Coulomb scale are real scientific risks and may affect the robustness of the crossing, but they are not equivalence-by-construction or tautological reduction. No circular step satisfying the quoted-evidence requirement is present.

Axiom & Free-Parameter Ledger

8 free parameters · 7 axioms · 0 invented entities

The model has two genuinely load-bearing fitted inputs — c_2 and c_4 (with ε) — placed to reproduce the experimental transition densities, plus the usual VMC variational parameters. The self-doping conclusion is an output of the variational comparison, but it is computed in a restricted single-vacancy subspace with a partial optimization and at an energy scale (≈1% E_c) comparable to the neglected band effects listed in the axioms. No new physical entities (particles, forces, dimensions) are introduced: the polaron is a variational correlation factor F_R({r_i}) (Eq. 12), not an independent quasiparticle.

free parameters (8)
  • c_2 (quadratic band coefficient) = −190 meV nm²
    Chosen, jointly with c_4 and ε, so that the optimized QFL–CDW energy comparison reproduces the experimental Lifshitz transition at ≈0.7×10¹² cm⁻² (App. A.5).
  • c_4 (quartic band coefficient) = 1000 meV nm⁴
    Chosen so that the CDW wins over the QFL at ≈0.4–0.5×10¹² cm⁻², matching the experimental crystallization density (App. A.5).
  • dielectric constant ε = 10
    Set by hand (Introduction) to fix the Coulomb scale e²√n_e/ε = 10.2 meV at n_e = 0.5×10¹² cm⁻².
  • Jastrow parameters A, B, A_LR = grid-scanned (A: 0.5–12, B: 0.1–0.4, A_LR: 0–15 in n_e=1 units)
    Variational VMC parameters optimized per state and per density (App. A.5); required by the Slater–Jastrow ansatz.
  • auxiliary potential V_0 and diagonalization density n_diag = grid-scanned (V_0: 0.2–20 meV; n_diag: 0.2–1.2×10¹² cm⁻²)
    Define the CDW orbital subspace via H_aux (Eq. 4) and control how orbitals are rescaled across densities.
  • polaron amplitude A_p and width ξ_p = A_p up to 1.8; ξ_p 0.5–2.0 in √n_e units (optimized 1.2 and 0.8 at n_e=0.4×10¹² cm⁻²)
    Fitted per density with the other five parameters frozen (partial optimization, App. A.5); the polaron energy gain (0.5–1% E_c) is the load-bearing effect for the self-doping claim.
  • aWC variational set C, k_F, ξ = e.g., C=3, ξ=1, k_F=4.0 at n_e=0.4×10¹² cm⁻²
    Gaussian width, momentum cutoff, and smearing for the annular-WC orbital (Eq. 2); variational.
  • E_offset polynomial coefficients = a = −3.23, b = 2.71, c = −1.60
    Density-dependent offset applied to all energies in Figs. 2–3 for display; not physical, and it does not change relative energies.
axioms (7)
  • standard math Slater–Jastrow variational principle: sampled local-energy expectation is an upper bound on the ground-state energy of H = Σ(c_2 k² + c_4 k⁴) + Coulomb.
    Used for every energy comparison in Figs. 2–3.
  • domain assumption The universal Jastrow form (Eqs. 7–8) with u_LR(q) ∝ (q² + q_0²)^{-5/4} adequately captures correlations in the quartic-dispersion Wigner crystal.
    App. A.3 derives the 1/q^{5/2} Gaskell form from RPA-type reasoning for a quartic dispersion, but the adequacy of the form is not independently verified.
  • domain assumption The cusp condition u'(0) = 0 holds for the quartic dispersion.
    App. A.3; used to argue universality of the Jastrow factor, stated without an exact derivation.
  • domain assumption The lowest band of H_aux (Eq. 4) with momentum cutoff |G| ≤ 3|b_1| spans the relevant crystalline orbital subspace.
    App. A.1 gives a 4th-shell suppression estimate (δE ~ |V_0|²/(256 c_4 |b_1|⁴)) but no convergence check at 61 bands.
  • domain assumption Neglect of quantum geometry, trigonal warping, and band-topological effects in the Mexican-hat model.
    Introduction; the paper acknowledges these affect energetics at the few-percent-of-E_c scale (Discussion), the same scale as the self-doping energy gain.
  • domain assumption Single-vacancy subspace: exactly one hole on the torus; hole–hole interactions and multi-polaron coherence are ignored.
    63-electron simulations and Eq. 13 fit a single-hole tight-binding model; Discussion states that hole–hole interaction is ignored.
  • standard math The Wannier Bloch phase convention (ψ_{1k} real and positive at r=0) is admissible for the polaron determinant.
    App. B footnote 47; a gauge choice for Eq. 9 that removes relative sign ambiguities.

pith-pipeline@v1.3.0-alltime-deepseek · 5445 in / 6091 out tokens · 247804 ms · 2026-08-04T01:07:46.778575+00:00 · methodology

0 comments
read the original abstract

We study Wigner crystals (WCs) induced by a strong Coulomb interaction from the ring-like Fermi surface of a Mexican-hat dispersion $\epsilon_k= c_2k^2+c_4k^4$. We design orbital shape in order to minimize the energy of the WC, and find that a low ground-state energy requires an orbital shape with a depletion of electrons near $k=0$. To capture the Coulomb-induced correlations, we include a Jastrow factor as well as a factor describing the correlation between electrons and doped vacancies. Using variational Monte Carlo calculations, we calibrate the effective band parameters $c_2$ and $c_4$ to reproduce the two transitions observed experimentally as the electron density is lowered: from a spin-valley-polarized Fermi liquid with a disk-like Fermi surface, to one with a ring-like Fermi surface, and finally to a WC. We find that, even with an optimized orbital shape that depletes electrons near $k=0$, a WC with hole self-doping near $k=0$ can still be energetically favorable near the WC transition, provided that the electron-vacancy correlation is included. We also estimate the dispersion of the doped hole.

Figures

Figures reproduced from arXiv: 2608.00170 by Minho Luke Kim, Xiao-Gang Wen.

Figure 1
Figure 1. Figure 1: (Left) In the weak-interaction limit, an induced CDW may favor a lattice constant that matches the outer edge of the ring-like Fermi surface associated with the Mexican-hat dispersion, thereby gapping out the outer Fermi surface. The inner Fermi surface remains, giving rise to an mWC state. The shaded region in momentum space denotes the occupied electronic states, and kWC is the momentum transfer generate… view at source ↗
Figure 2
Figure 2. Figure 2: Energies of the commensurate crystal ans¨atze [PITH_FULL_IMAGE:figures/full_fig_p002_2.png] view at source ↗
Figure 3
Figure 3. Figure 3: (Top) Energies of the commensurate CDW and aWC compared with those of one-vacancy hCDW and dressed pCDW states with a hole at Γ, where the hole creation is most energetically competitive. We find that the hCDWΓ is energetically competitive from around 0.4 × 1012 cm−2 . The dressed pCDWΓ state performs even better with noticeable improvements compared to CDW and hCDWΓ. Range of errors are given in the ribbo… view at source ↗
Figure 4
Figure 4. Figure 4: Difference of the total energy of the system in units of Coulomb energy for hCDWK and hCDWM states against hCDWΓ (in lines) and pCDWK and pCDWM states against pCDWΓ (in dashed lines). The energy difference can be used to fit coefficients of the effective hole tight-binding model for the nearest and next-nearest neighbor terms. Errors for hCDWK is displayed; other states have similar size of errors. Therefo… view at source ↗
Figure 5
Figure 5. Figure 5: Energies of the doped CDW ans¨atze per electron compared against hCDW [PITH_FULL_IMAGE:figures/full_fig_p015_5.png] view at source ↗
Figure 6
Figure 6. Figure 6: Visualization of the single-determinant density for the gWC (a), aWC (b), and commensurate CDW (c) [PITH_FULL_IMAGE:figures/full_fig_p016_6.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

64 extracted references · 5 linked inside Pith

  1. [1]

    On the interaction of electrons in metals,

    E. Wigner, “On the interaction of electrons in metals,” Phys. Rev.46(Dec, 1934) 1002–1011

    1934

  2. [2]

    Ground state of the two-dimensional electron gas,

    B. Tanatar and D. M. Ceperley, “Ground state of the two-dimensional electron gas,”Phys. Rev. B39(Mar,

  3. [3]

    Correlation energy and spin polarization in the 2d electron gas,

    C. Attaccalite, S. Moroni, P. Gori-Giorgi, and G. B. Bachelet, “Correlation energy and spin polarization in the 2d electron gas,”Phys. Rev. Lett.88(Jun, 2002) 256601

    2002

  4. [4]

    Phase diagram of the low-density two-dimensional homogeneous electron gas,

    N. D. Drummond and R. J. Needs, “Phase diagram of the low-density two-dimensional homogeneous electron gas,” Phys. Rev. Lett.102(Mar, 2009) 126402

    2009

  5. [5]

    Chiral wigner crystal phases induced by berry curvature,

    S. Joy, L. Levitov, and B. Skinner, “Chiral wigner crystal phases induced by berry curvature,”Phys. Rev. Lett.135 (Dec, 2025) 256502

    2025

  6. [6]

    Wigner crystal of a two-dimensional electron gas with a strong spin-orbit interaction,

    P. G. Silvestrov and O. Entin-Wohlman, “Wigner crystal of a two-dimensional electron gas with a strong spin-orbit interaction,”Phys. Rev. B89(Apr, 2014) 155103

    2014

  7. [7]

    Phases intermediate be- tween a two-dimensional electron liquid and wigner crys- tal,

    B. Spivak and S. A. Kivelson, “Phases intermediate be- tween a two-dimensional electron liquid and wigner crys- tal,”Phys. Rev. B70(Oct, 2004) 155114

    2004

  8. [8]

    Self-doping instabil- ity of the wigner-mott insulator,

    S. Pankov and V. Dobrosavljevi´ c, “Self-doping instabil- ity of the wigner-mott insulator,”Phys. Rev. B77(Feb,

  9. [9]

    Charge density waves in a quantum plasma,

    Z. Han, S. Zhang, and X. Dai, “Charge density waves in a quantum plasma,”Phys. Rev. B100(Oct, 2019) 155132

    2019

  10. [10]

    Dynamical defects in a two-dimensional wigner crystal: Self-doping and kinetic magnetism,

    K.-S. Kim, I. Esterlis, C. Murthy, and S. A. Kivelson, “Dynamical defects in a two-dimensional wigner crystal: Self-doping and kinetic magnetism,”Phys. Rev. B109 (Jun, 2024) 235130

    2024

  11. [11]

    Single and paired point defects in a 2d wigner crystal,

    L. Cˆ andido, P. Phillips, and D. M. Ceperley, “Single and paired point defects in a 2d wigner crystal,”Phys. Rev. Lett.86(Jan, 2001) 492–495

    2001

  12. [12]

    Andreev- lifshitz supersolid revisited for a few electrons on a square lattice. i,

    G. Katomeris, F. Selva, and J.-L. Pichard, “Andreev- lifshitz supersolid revisited for a few electrons on a square lattice. i,”The European Physical Journal B - Condensed Matter and Complex Systems31no. 3, (Feb, 2003) 401– 412

    2003

  13. [13]

    Hybrid phase at the quantum melting of the wigner crystal,

    H. Falakshahi and X. Waintal, “Hybrid phase at the quantum melting of the wigner crystal,”Phys. Rev. Lett. 94(Jan, 2005) 046801

    2005

  14. [14]

    Persistent currents in two dimensions: New regimes induced by the interplay between electronic correlations and disorder,

    Z. A. N´ emeth and J.-L. Pichard, “Persistent currents in two dimensions: New regimes induced by the interplay between electronic correlations and disorder,”The Euro- pean Physical Journal B - Condensed Matter and Com- plex Systems45no. 1, (May, 2005) 111–128

    2005

  15. [15]

    Imaging two-dimensional generalized wigner crystals,

    H. Li, S. Li, E. C. Regan, D. Wang, W. Zhao, S. Kahn, K. Yumigeta, M. Blei, T. Taniguchi, K. Watanabe, S. Tongay, A. Zettl, M. F. Crommie, and F. Wang, “Imaging two-dimensional generalized wigner crystals,” Nature597no. 7878, (Sep, 2021) 650–654

    2021

  16. [16]

    Fractional quantum anomalous hall effect in multilayer graphene,

    Z. Lu, T. Han, Y. Yao, A. P. Reddy, J. Yang, J. Seo, K. Watanabe, T. Taniguchi, L. Fu, and L. Ju, “Fractional quantum anomalous hall effect in multilayer graphene,” Nature626no. 8000, (Feb, 2024) 759–764

    2024

  17. [17]

    Direct observation of a magnetic-field-induced wigner crystal,

    Y.-C. Tsui, M. He, Y. Hu, E. Lake, T. Wang, K. Watan- abe, T. Taniguchi, M. P. Zaletel, and A. Yazdani, “Direct observation of a magnetic-field-induced wigner crystal,” Nature628no. 8007, (2024) 287–292

    2024

  18. [18]

    Imaging quantum melting in a disordered 2d wigner solid,

    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. Ton- gay, S. G. Louie, M. F. Crommie, and F. Wang, “Imaging quantum melting in a disordered 2d wigner solid,”Science388no. 6748, (2025) 736–740, https://www.science.org/doi/pdf/10.1126/science.ado7136

    doi:10.1126/science.ado7136 2025

  19. [19]

    Visualizing the im- pact of quenched disorder on 2d electron wigner solids,

    Z. Ge, C. Smith, Z. He, Y. Yang, Q. Li, Z. Xiang, J. Xiao, W. Zhou, S. Kahn, M. Erdi,et al., “Visualizing the im- pact of quenched disorder on 2d electron wigner solids,” arXiv:2510.12009

    Pith/arXiv arXiv

  20. [20]

    Compatibility of crystalline order and the quantized hall effect,

    B. I. Halperin, Z. Teˇ sanovi´ c, and F. Axel, “Compatibility of crystalline order and the quantized hall effect,”Phys. Rev. Lett.57(Aug, 1986) 922

    1986

  21. [21]

    Cooperative ring exchange theory of the fractional quantized hall effect,

    S. Kivelson, C. Kallin, D. P. Arovas, and J. R. Schrief- fer, “Cooperative ring exchange theory of the fractional quantized hall effect,”Phys. Rev. Lett.56(Feb, 1986) 873–876

    1986

  22. [22]

    Theory of quan- tum anomalous hall phases in pentalayer rhombohedral graphene moir´ e structures,

    Z. Dong, A. S. Patri, and T. Senthil, “Theory of quan- tum anomalous hall phases in pentalayer rhombohedral graphene moir´ e structures,”Phys. Rev. Lett.133(Nov,

  23. [23]

    Anomalous hall crys- tals in rhombohedral multilayer graphene. i. interaction- driven chern bands and fractional quantum hall states at zero magnetic field,

    J. Dong, T. Wang, T. Wang, T. Soejima, M. P. Zaletel, A. Vishwanath, and D. E. Parker, “Anomalous hall crys- tals in rhombohedral multilayer graphene. i. interaction- driven chern bands and fractional quantum hall states at zero magnetic field,”Phys. Rev. Lett.133(Nov, 2024) 206503

    2024

  24. [24]

    Fractional quan- tum anomalous hall effect in rhombohedral multilayer graphene in the moir´ eless limit,

    B. Zhou, H. Yang, and Y.-H. Zhang, “Fractional quan- tum anomalous hall effect in rhombohedral multilayer graphene in the moir´ eless limit,”Phys. Rev. Lett.133 (Nov, 2024) 206504

    2024

  25. [25]

    Anomalous hall crystals in rhombohedral multilayer graphene. ii. general mechanism and a minimal model,

    T. Soejima, J. Dong, T. Wang, T. Wang, M. P. Zale- tel, A. Vishwanath, and D. E. Parker, “Anomalous hall crystals in rhombohedral multilayer graphene. ii. general mechanism and a minimal model,”Phys. Rev. B110 (Nov, 2024) 205124

    2024

  26. [26]

    Stability of anoma- lous hall crystals in multilayer rhombohedral graphene,

    Z. Dong, A. S. Patri, and T. Senthil, “Stability of anoma- lous hall crystals in multilayer rhombohedral graphene,” Phys. Rev. B110(Nov, 2024) 205130

    2024

  27. [27]

    Y. H. Kwan, J. Yu, J. Herzog-Arbeitman, D. K. Efetov, N. Regnault, and B. A. Bernevig, “Moir´ e fractional chern 7 insulators. iii. hartree-fock phase diagram, magic angle regime for chern insulator states, role of moir´ e potential, and goldstone gaps in rhombohedral graphene superlat- tices,”Phys. Rev. B112(Aug, 2025) 075109

    2025

  28. [28]

    Moir´ e fractional chern insulators. iv. fluctuation-driven collapse in multiband exact di- agonalization calculations on rhombohedral graphene,

    J. Yu, J. Herzog-Arbeitman, Y. H. Kwan, N. Regnault, and B. A. Bernevig, “Moir´ e fractional chern insulators. iv. fluctuation-driven collapse in multiband exact di- agonalization calculations on rhombohedral graphene,” Phys. Rev. B112(Aug, 2025) 075110

    2025

  29. [29]

    Parent berry curvature and the ideal anomalous hall crystal,

    T. Tan and T. Devakul, “Parent berry curvature and the ideal anomalous hall crystal,”Phys. Rev. X14(Nov,

  30. [30]

    Quantum melting a wigner crystal into hall liquids,

    A. P. Reddy and L. Fu, “Quantum melting a wigner crystal into hall liquids,”Phys. Rev. B113(Apr, 2026) L161403

    2026

  31. [31]

    Variational wave-function analysis of the fractional anomalous hall crystal,

    T. Tan, J. May-Mann, and T. Devakul, “Variational wave-function analysis of the fractional anomalous hall crystal,”Phys. Rev. Lett.135(Jul, 2025) 036604

    2025

  32. [32]

    New classes of quantum anomalous hall crystals in multilayer graphene,

    B. Zhou and Y.-H. Zhang, “New classes of quantum anomalous hall crystals in multilayer graphene,”Phys. Rev. Lett.135(Jul, 2025) 036501

    2025

  33. [33]

    Extended quantum anomalous hall states in graphene/hbn moir´ e superlattices,

    Z. Lu, T. Han, Y. Yao, Z. Hadjri, J. Yang, J. Seo, L. Shi, S. Ye, K. Watanabe, T. Taniguchi, and L. Ju, “Extended quantum anomalous hall states in graphene/hbn moir´ e superlattices,”Nature637no. 8048, (Jan, 2025) 1090– 1095

    2025

  34. [34]

    Chern insulators at integer and frac- tional filling in moir´ e pentalayer graphene,

    D. Waters, A. Okounkova, R. Su, B. Zhou, J. Yao, K. Watanabe, T. Taniguchi, X. Xu, Y.-H. Zhang, J. Folk, and M. Yankowitz, “Chern insulators at integer and frac- tional filling in moir´ e pentalayer graphene,”Phys. Rev. X15(Feb, 2025) 011045

    2025

  35. [35]

    Signatures of Chiral Superconductivity in Rhombohedral Graphene,

    T. Han, Z. Lu, Y. Yao, L. Shi, J. Yang, J. Seo, S. Ye, Z. Wu, M. Zhou, H. Liu, G. Shi, Z. Hua, K. Watanabe, T. Taniguchi, P. Xiong, L. Fu, and L. Ju, “Signatures of Chiral Superconductivity in Rhombohedral Graphene,” (Aug., 2024) ,arXiv:2408.15233

    Pith/arXiv arXiv 2024

  36. [36]

    Extreme anisotropy in the metallic and superconducting phases of rhombohedral hexalayer graphene,

    P. Qin, H.-T. Wu, R. Q. Nguyen, E. Morissette, N. J. Zhang, K. Watanabe, T. Taniguchi, and J. Li, “Extreme anisotropy in the metallic and superconducting phases of rhombohedral hexalayer graphene,”arXiv:2504.05129

    Pith/arXiv arXiv

  37. [37]

    Evidence of metallic wigner crystal in rhombohedral graphene,

    T. Han, J. P. Butler, S. Ye, Z. Hua, S. Dutta, Z. Had- jri, Z. Wu, J. Yang, J. Seo, P. Pattanakanvijit,et al., “Evidence of metallic wigner crystal in rhombohedral graphene,”arXiv:2604.00113

    arXiv

  38. [38]

    Sublattice structure and topology in spontaneously crystallized electronic states,

    Y. Zeng, D. Guerci, V. Cr´ epel, A. J. Millis, and J. Cano, “Sublattice structure and topology in spontaneously crystallized electronic states,”Phys. Rev. Lett.132(Jun,

  39. [39]

    Various electronic crystal phases in rhombohedral graphene multilayers,

    W. Miao and C. Li, “Various electronic crystal phases in rhombohedral graphene multilayers,”Phys. Rev. B113 (Apr, 2026) 155136

    2026

  40. [40]

    λ-jellium model for the anomalous hall crystal,

    T. Soejima, J. Dong, A. Vishwanath, and D. E. Parker, “λ-jellium model for the anomalous hall crystal,”Phys. Rev. Lett.135(Oct, 2025) 186505

    2025

  41. [41]

    Crystals caught doping: Metallic wigner crystals in rhombohedral graphene,

    J. Dong, T. Soejima, D. E. Parker, and A. Vishwanath, “Crystals caught doping: Metallic wigner crystals in rhombohedral graphene,”arXiv:2604.00114

    arXiv

  42. [42]

    Self- doped crystal from preempted band-inversion transi- tions,

    J. Feng, Z. Han, M. P. Zaletel, and Z. Dong, “Self- doped crystal from preempted band-inversion transi- tions,”arXiv:2604.09820

    Pith/arXiv arXiv

  43. [43]

    Topolog- ical chiral superconductivity beyond pairing in a fermi liquid,

    M. Kim, A. Timmel, L. Ju, and X.-G. Wen, “Topolog- ical chiral superconductivity beyond pairing in a fermi liquid,”Phys. Rev. B111(Jan, 2025) 014508

    2025

  44. [44]

    Variational monte carlo optimization of topological chiral supercon- ductors,

    M. L. Kim, A. Timmel, and X.-G. Wen, “Variational monte carlo optimization of topological chiral supercon- ductors,”Phys. Rev. B114(Jul, 2026) 014506

    2026

  45. [45]

    References [2, 10, 37, 53, 54] are relevant in the supple- mental material

  46. [46]

    While it does not matter due to the symmetry, it is con- venient to defineG 1 andG 2 to be 120 degrees apart and G3 =−G 1 −G 2

  47. [47]

    We choose the convention that for each Bloch eigenvector, the wavefunction atr= 0 is real and positive

    We note that the Bloch functions have arbitrary phases: ψ1k →e iθ(k)ψ1k. We choose the convention that for each Bloch eigenvector, the wavefunction atr= 0 is real and positive

  48. [48]

    Here, since there areKvec- tors which translateK/Mto their respective degenerate points, there is guaranteed degenerate mixing between these orbitals

    The physics is different forKandMsectors as these are at the Brillouin zone edges. Here, since there areKvec- tors which translateK/Mto their respective degenerate points, there is guaranteed degenerate mixing between these orbitals

  49. [49]

    Wigner crystallization in bernal bilayer graphene,

    S. Joy and B. Skinner, “Wigner crystallization in bernal bilayer graphene,”arXiv:2310.07751

    Pith/arXiv arXiv

  50. [50]

    Rhombohedral graphene: A tale of many crystals,

    A. Abouelkomsan, F. Gaggioli, D. Guerci, and L. Fu, “Rhombohedral graphene: A tale of many crystals,”to appear

  51. [51]

    Second, technically, just enforcing bothk and−kare occupied is enough, butL x =L y gives the additionalC 4 symmetry

    First, the centralk= 0 point is degenerate under 4- fold rotation. Second, technically, just enforcing bothk and−kare occupied is enough, butL x =L y gives the additionalC 4 symmetry

  52. [52]

    In reality, this is a smoothly decaying function, and for this Jastrow to decay by 1/e, the effective length isr≈ 3.21B

  53. [53]

    Self- attention neural network for solving correlated electron problems in solids,

    M. Geier, K. Nazaryan, T. Zaklama, and L. Fu, “Self- attention neural network for solving correlated electron problems in solids,”Phys. Rev. B112(Jul, 2025) 045119

    2025

  54. [54]

    Some static and dy- namical properties of a two-dimensional wigner crystal,

    L. Bonsall and A. A. Maradudin, “Some static and dy- namical properties of a two-dimensional wigner crystal,” Phys. Rev. B15(Feb, 1977) 1959–1973. 8 Appendix A: Methodology In this section, we present the details of the construction of the trial wavefunctions and their properties. In the end, we also discuss the variational Monte Carlo routine and the sam...

    1977

  55. [55]

    To ensure that this is well-defined on a torus, we first define a set of lattice points

    The CDW auxiliary periodic potential model We employed a triangular lattice configuration for the Wigner crystal. To ensure that this is well-defined on a torus, we first define a set of lattice points. Defining the lattice vectors asa 1 = (1,0) anda 2 = 1 2 , √ 3 2 , the lattice points for the Wigner crystal are positioned atR m,n =ma 1 +na 2 for 0≤m, n ...

  56. [56]

    All numerics are run atn e = 1 units, and therefore the whole torus is then scaled by the factor q 2√ 3 so that the total area is ensured to be 64. As outlined in the main letter, the individual orbitals of the commensurate Wigner crystal is obtained by introducing a periodic potential as an auxiliary potential, which will act as creating orbitals which a...

  57. [57]

    The individual single particle orbitals are now given as plane wave statese ikj ·ri

    Quarter F ermi Liquid Defining the torus directions as ˆx= (1,0) and ˆy= (0,1), this time we set the torus lengths to be equal on both directions. The individual single particle orbitals are now given as plane wave statese ikj ·ri . For the ordinary Fermi liquid, we fill fromk= 0 until 65 states are occupied, ensuring that the 90-degree rotation symmetry ...

  58. [58]

    We define the short-ranged Jastrowu SR as uSR =−A 1 + r B + r2 2B2 e− r B ,(A7) which satisfies the cusp conditionu ′ SR = 0.Btherefore effectively sets the range of theu SR[52]

    Electron correlation Jastrow F actor The Jastrow factor used for all our trial states involve a short-range correlator and a long-range correlator. We define the short-ranged Jastrowu SR as uSR =−A 1 + r B + r2 2B2 e− r B ,(A7) which satisfies the cusp conditionu ′ SR = 0.Btherefore effectively sets the range of theu SR[52]. For all purposes,B would be th...

  59. [59]

    Kinetic Energy for Slater-Jastrow wavefunction As per VMC protocol, kinetic energy is sampled by taking the local energy

    Sampling the Energy a. Kinetic Energy for Slater-Jastrow wavefunction As per VMC protocol, kinetic energy is sampled by taking the local energy. Given an operator, the expectation value is given as ⟨O⟩= R ψ∗OψdRR ψ∗ψdR = R ψ∗ψ Oψ ψ dRR ψ∗ψdR = Oψ ψ M C .(A19) 11 whereRare the collective coordinates for all particles involved. Therefore, for the kinetic en...

  60. [60]

    First, the QFL, gWC, and aWC result samples⟨k 2⟩/ne,⟨k 4⟩/n2 e, and⟨ 1 r ⟩/√ne

    Parameter selection and optimization For the selection of the Mexican hat parameter, we take a semi-iterative approach. First, the QFL, gWC, and aWC result samples⟨k 2⟩/ne,⟨k 4⟩/n2 e, and⟨ 1 r ⟩/√ne. Therefore, we can scale these results to produce the Fermi liquid energy at arbitary densities: E(ne) N =c 2ne⟨k2⟩ne=1 +c 4ne⟨k4⟩ne=1 + e2√ne ϵ ⟨V⟩ ne=1 (A26...

  61. [61]

    For each run, we equilibriate for 10 5 steps, with each step moving one particle for a distance of 0.03Lin case of the Fermi liquid and 0.03L 2 in case of the Wigner crystal

    VMC Routine and Sampling As per standard VMC protocol, we sample the probability distribution|Ψ({r i})|2 with the Metropolis-Hastings algorithm. For each run, we equilibriate for 10 5 steps, with each step moving one particle for a distance of 0.03Lin case of the Fermi liquid and 0.03L 2 in case of the Wigner crystal. To avoid autocorrelation effects, we ...

  62. [62]

    Kinetic Energy for Polaron dressed states The methodology presented in the previous section only applies for ordinary Slater-Jastrow wavefunction has to be modified for the polaronic pCDW states. Starting from the pCDW wavefunction stated: ΨpCDW k ({ri}) =J({r i}) X R e−ik·RFR({ri})DR({ri}),(B5) a useful trick is that this wavefunction can still be writte...

  63. [63]

    It is instructive to see how much energy gain per particle the polaron dressing induces, including for pCDW K and pCDW M states against their hCDW counterparts

    Improvement achieved by pCDW compared to hCDW states In the main letter, we presented the energetics result for pCDW and hCDW with a hole at Γ. It is instructive to see how much energy gain per particle the polaron dressing induces, including for pCDW K and pCDW M states against their hCDW counterparts. Figure 5 shows the difference of all doped states at...

  64. [64]

    Single particle orbitals with removed contributions near k= 0 A classic Wigner crystal ansatz in real space is in Slater-Jastrow form with Gaussian single particle orbitals, which can be written as ΨgW C({r}) = det[ψi(rj)] exp  X i<j u(|ri −r j|)   .(C1) where{r}is the collective electron coordinates,r i is the electron coordinate for electroni, andu(...