Pith. sign in

REVIEW 1 major objections 5 minor 32 references

Self-consistent Hartree-Fock maps WSe2 Wigner crystal density, bands, and screening

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 · glm-5.2

2026-07-07 14:11 UTC pith:3R4K54RB

load-bearing objection Sound HF Wigner-crystal framework applied to WSe2 with Keldysh interaction; numerical results need convergence evidence the 1 major comments →

arxiv 2607.05371 v1 pith:3R4K54RB submitted 2026-07-06 cond-mat.mes-hall

Dielectric function in WSe2

Tiberius O. Cheche , Yia-Chung Chang This is my paper
classification cond-mat.mes-hall
keywords Wigner crystalHartree-Fockdielectric functionWSe2Keldysh interactionmonolayer TMDWigner crystalstatic screening
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 develops a self-consistent Hartree-Fock method for a fully spin-polarized two-dimensional Wigner crystal on a triangular lattice, using the Keldysh interaction appropriate for monolayer transition-metal dichalcogenides. Applied to WSe2, it computes the quasiparticle band structure, the real-space crystalline charge density, and the static dielectric function ε(q,0) from a Lindhard-type interband polarizability built on the converged HF eigenvalues and eigenvectors. The central physical result is a unified microscopic picture in which increasing the density parameter rs (lower carrier density, stronger interactions) simultaneously strengthens the crystalline density modulation—quantified by the first-star order parameter M1 and the density-contrast ratio Cn—and reduces the finite-q electronic screening, pushing the peak of ε(q,0) toward unity and toward smaller momenta. The paper explicitly states these are signatures of crystalline order within the HF solution, not a thermodynamic proof of a liquid-to-crystal phase transition.

Core claim

The paper's core contribution is the construction of a single self-consistent framework that links four quantities for a spin-polarized Wigner crystal in monolayer WSe2: the real-space charge density (characterized by the order parameter M1 and contrast Cn), the Hartree-Fock band structure in the Wigner-crystal Brillouin zone, the static dielectric function εHF(q,0), and the density parameter rs. The key quantitative finding is that as rs increases from 4 to 50, M1 grows from 0.100 to 0.586, Cn grows from 2.74 to 1.44×10^3, and the peak of εHF(q,0) systematically decreases toward unity while shifting to smaller q. This means the electrons become more localized in real space and less able to屏

What carries the argument

The self-consistent Hartree-Fock cycle iterates between density harmonics ηW(Q) and plane-wave eigenvector coefficients ZnG(k) until convergence. The Keldysh interaction vK(q) = 2πe²/(εq(1+qρ0)) provides the electron-electron potential. The static dielectric function εHF(q,0) = 1 + vK(q)·ΠHF(q,0) is computed from an interband Lindhard-type polarizability using the HF band gap and density matrix elements Mnm(k,q), with momentum folding into the Wigner-crystal Brillouin zone.

Load-bearing premise

The calculation uses restricted Hartree-Fock with full spin polarization and no correlation effects, which the authors acknowledge does not constitute a thermodynamic proof of the Wigner-crystal phase and omits the self-consistent two-spin treatment needed for a genuine unpolarized crystal.

What would settle it

If a self-consistent two-spin or correlation-corrected calculation yields qualitatively different band structures, order parameters, or dielectric functions for the same rs values in WSe2, the single-spin HF predictions would be superseded.

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

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

1 major / 5 minor

Summary. The manuscript develops a self-consistent Hartree-Fock (HF) method for computing the band structure and static dielectric function of a two-dimensional Wigner crystal, applied to monolayer WSe2. The formulation uses a plane-wave basis in the Wigner-crystal Brillouin zone with a Keldysh interaction. The dielectric response is computed from a static Lindhard-type (Adler-Wiser) polarizability using the converged HF eigenvalues and eigenvectors. The authors present results for rs = 4, 10, 30, 50, showing that increasing rs strengthens the crystalline density modulation (quantified by the first-star order parameter M1 and density contrast Cn) while reducing the finite-q electronic screening. The long-wavelength limit epsilon_HF(q->0,0) -> 1 is correctly derived and numerically verified. The theoretical framework is standard and internally consistent.

Significance. The manuscript provides a clearly formulated, self-consistent HF framework connecting real-space crystalline density, quasiparticle spectrum, and static dielectric response within a single microscopic calculation adapted to the Keldysh interaction of TMD monolayers. The derivation of the long-wavelength limit (Eq. 17, SI.3) is a useful internal consistency check. The proposed experimental protocol (helicity-resolved optical spectroscopy on spin/valley-polarized WSe2) provides a falsifiable connection to ongoing experimental efforts. The work is a reasonable methodological contribution to the study of low-density 2D systems. However, the quantitative results are currently undermined by a complete absence of convergence data, which prevents independent verification of the reported numerical values.

major comments (1)
  1. §3.1 and Table S1: The manuscript states that 'convergence of the results was checked against the number of retained reciprocal vectors, the wedge mesh density, the number of empty bands retained in the polarizability, and the broadening parameter eta_Pi,' but no quantitative convergence data are provided anywhere in the main text or SI. Table S1 lists parameter ranges (N_side = 5-20, N_G = 50-120, N_b = 13-17) without specifying which values were used for the production results in Fig. 2. This is load-bearing: the central quantitative outputs (M1, Cn, and the peak values of epsilon_HF(q,0)) depend on these parameters, and without at least a convergence table or plot showing how much these quantities shift with parameter variation, the reported numbers cannot be verified or trusted. A table showing, e.g., the variation of M1 and max(epsilon_HF) with N_G and N_side at a representative rs值
minor comments (5)
  1. §3.1: The manuscript acknowledges multiple Hartree-preconditioned attractors and states that branch selection used qualitative criteria (symmetry preservation, band regularity, degeneracy pattern). A brief comment on the sensitivity of the final results to the initial Gaussian width beta_0 would strengthen the manuscript. The rs=30 example (M_ini=0.668 -> M_fin=0.504) shows the solution moves away from the initial guess, but it is not shown whether different beta_0 values converge to the same branch.
  2. Fig. 2: The dielectric function panels (d, h, l, p) have y-axis labels with overlapping tick values (e.g., '1.1' and '1.1' appearing separately). This should be cleaned up for readability.
  3. §3.3, Eq. (22): The unpolarized estimate (scaling exchange by 1/2) is described but not used for quantitative analysis. A brief comment on the expected magnitude of error from this approximation, or removal of the equation if it is not used, would improve clarity.
  4. References: Refs 11 and 23 appear to be duplicate citations of the same arXiv preprint (arXiv:2512.16631). This should be consolidated.
  5. The code is stated to be available 'upon reasonable request.' For reproducibility, depositing the code in a public repository would be preferable.

Simulated Author's Rebuttal

1 responses · 0 unresolved

We thank the referee for the careful reading and the constructive assessment. The referee's central concern—the absence of quantitative convergence data—is well taken, and we will address it in the revised manuscript.

read point-by-point responses

  1. Referee: §3.1 and Table S1: The manuscript states that convergence was checked against N_side, N_G, N_b, and eta_Pi, but no quantitative convergence data are provided. Table S1 lists parameter ranges without specifying which values were used for production results. A convergence table or plot is needed to verify the reported numbers.

    Authors: The referee is correct. The current manuscript states that convergence was checked but does not provide the quantitative data necessary for independent verification, nor does it specify which parameter values within the ranges listed in Table S1 were used for the production results in Fig. 2. This is a legitimate deficiency that we will remedy in the revised version. revision: yes

Circularity Check

0 steps flagged

No circularity found: HF equations, Keldysh interaction, and Lindhard-type polarizability are standard external results; dielectric function computed from HF bands without fitting to target quantities.

full rationale

The derivation chain is self-contained and does not exhibit circularity. The Hartree-Fock equations (Eqs. 2-6), the Keldysh interaction (Eq. 7), and the static Lindhard-type polarizability (Eq. 14) are all standard results from external literature (Lindhard [18], Stern [19], Adler [20], Wiser [21], Keldysh [13-15]). The material parameters (m*=0.4, ε=4, ρ0=1.12 nm) are taken from prior literature (Refs 13-15). The dielectric function εHF(q,0) = 1 + evK(q)ΠHF(q,0) (Eq. 16) is computed from the converged HF eigenvalues and eigenvectors via the independent-particle polarizability (Eq. 14), which involves occupied-empty interband transitions — a standard Adler-Wiser formulation for insulating periodic solids. No parameter is fitted to the dielectric function and then 'predicted' back. The order parameter M1 (Eq. 20) and density contrast Cn are computed from the self-consistent density, not fitted. The initial Gaussian width β0 is a starting guess for iteration, not a fitted parameter that constrains the output. The paper explicitly acknowledges that M1 changes from initial to converged values (e.g., rs=30: Mini=0.668 → Mfin=0.504), demonstrating the output is not trivially the input. The long-wavelength limit εHF(0,0)=1 (Eq. 17) is derived analytically from orthonormality and the insulating gap, not imposed as a fit. The convergence checks (Sec. 3.1) and branch selection criteria are numerical practice concerns, not circularity. The central claim — that increasing rs strengthens crystalline density modulation while reducing finite-q screening — follows from the self-consistent HF solution without any step where the output is defined in terms of the target result.

Axiom & Free-Parameter Ledger

7 free parameters · 5 axioms · 0 invented entities

No new physical entities, particles, forces, or dimensions are postulated. The calculation uses established theoretical frameworks (Hartree-Fock, Keldysh interaction, Adler-Wiser polarizability) applied to a known material (WSe2).

free parameters (7)
  • m* = 0.4 (in units of m0)
    Electron effective mass in WSe2, taken from prior literature, not fitted in this work.
  • epsilon = 4
    Effective dielectric constant of the medium surrounding the monolayer, taken as a fixed input parameter.
  • rho_0 = 1.12 nm
    Keldysh screening length, derived from in-plane polarizability kappa and epsilon, taken from prior literature.
  • beta_0 = Not specified numerically
    Dimensionless Gaussian width controlling the initial density profile (Eq. 10). Used only for initialization; the paper argues the converged solution is independent of this choice, though no systematic demonstration is provided.
  • eta_Pi = 5.0e-4 Ry*
    Numerical broadening parameter in the static polarizability (Eq. 14). A convergence parameter, not physical.
  • G_cut = Not specified exactly; NG=50-120 retained vectors
    Plane-wave cutoff controlling the basis size. Range given in Table S1 but exact values per calculation not stated.
  • N_side = 5-20
    Wedge-mesh resolution. Range given but exact values per figure panel not stated.
axioms (5)
  • domain assumption The ground state of the Wigner crystal is fully spin-polarized.
    Stated in Sec. 2.1: 'the ground state is assumed to be fully spin-polarized.' This is a standard assumption for high-density Wigner crystals in some regimes but is not proven here for WSe2 at the densities considered.
  • domain assumption The restricted Hartree-Fock approximation (no correlation) is adequate for the Wigner-crystal state.
    The entire calculation is at the HF level. Correlation effects, known to be significant in low-density 2D systems, are neglected. The paper does not discuss this limitation quantitatively.
  • domain assumption The static Lindhard-type polarizability (Eq. 14) without vertex corrections describes the dielectric response.
    Eq. 14 is the independent-particle (Adler-Wiser) polarizability using HF eigenvalues. Vertex corrections, which modify the response in interacting systems, are not included.
  • standard math The Keldysh interaction (Eq. 7) with isotropic dielectric environment describes the carrier-carrier interaction in monolayer WSe2.
    Standard model for 2D TMD electrostatics, widely used in the literature (Refs 13-15).
  • domain assumption The triangular lattice is the ground-state Wigner-crystal structure.
    Assumed based on Bonsall and Maradudin (Ref 3) for classical 2D Wigner crystals. The paper does not verify this for the quantum HF solution in WSe2.

pith-pipeline@v1.1.0-glm · 20504 in / 3422 out tokens · 444547 ms · 2026-07-07T14:11:17.469805+00:00 · methodology

0 comments
read the original abstract

We develop a Hartree-Fock numerical method for computing the band structure of a two-dimensional Wigner crystal in an electron gas at zero temperature. The ground state is assumed to be fully spin-polarized. Single-particle excitation spectra are evaluated in spin-conserving channel. As an application, we use the developed code to compute the static dielectric function epsilon(q,0) of a Wigner-crystal state formed in a two-dimensional transition-metal dichalcogenide, specifically monolayer WSe2. The dielectric response is obtained from the Hartree-Fock band structure and eigenfunctions through a static Lindhard-type polarizability. The method provides a theoretical tool for investigating screening, band-structure reconstruction, and interaction effects in low-density two-dimensional systems, with possible relevance for future experimental studies.

Figures

Figures reproduced from arXiv: 2607.05371 by Tiberius O. Cheche, Yia-Chung Chang.

Figure 1
Figure 1. Figure 1: Hexagonal first Brillouin zone of the triangular Wigner-crystal lattice. The irreducible [PITH_FULL_IMAGE:figures/full_fig_p005_1.png] view at source ↗
Figure 2
Figure 2. Figure 2: Evolution of the fully spin-polarized Hartree–Fock Wigner-crystal solution with the density [PITH_FULL_IMAGE:figures/full_fig_p009_2.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

32 extracted references · 32 canonical work pages · 1 internal anchor

  1. [1]

    , title =

    Wigner, E. , title =. 1934 , volume =

    work page 1934

  2. [2]

    , title =

    Wigner, E. , title =. 1938 , volume =

    work page 1938

  3. [3]

    and Maradudin, A

    Bonsall, L. and Maradudin, A. A. , title =. 1977 , volume =

    work page 1977

  4. [4]

    Trail, J. R. and Towler, M. D. and Needs, R. J. , title =. 2003 , volume =

    work page 2003

  5. [5]

    and Delyon, F

    Bernu, B. and Delyon, F. and Holzmann, M. and Baguet, L. , title =. 2011 , volume =

    work page 2011

  6. [6]

    2025 , eprint =

    Jain, Ambuj and Huang, Chunli , title =. 2025 , eprint =

    work page 2025

  7. [7]

    Zhang, Lifu and Gu, Liuxin and Adlong, Haydn S. and Christianen, Arthur and Dizer, Eugen and Ni, Ruihao and Ma, Rundong and Park, Suji and Jang, Houk and Taniguchi, Takashi and Watanabe, Kenji and Esterlis, Ilya and Schmidt, Richard and Imamoglu, Atac and Zhou, You , title =. arXiv preprint , year =. 2512.16631 , eprinttype =

    work page arXiv

  8. [8]

    and Ceperley, D

    Tanatar, B. and Ceperley, D. M. , title =. 1989 , volume =

    work page 1989

  9. [9]

    Drummond, N. D. and Radnai, Z. and Trail, J. R. and Towler, M. D. and Needs, R. J. , title =. 2004 , volume =

    work page 2004

  10. [10]

    and Heo, Hoseok and Taniguchi, Takashi and Watanabe, Kenji and Zar

    Zhou, You and Sung, Jiho and Brutschea, Elise and Esterlis, Ilya and Wang, Yao and Scuri, Giovanni and Gelly, Ryan J. and Heo, Hoseok and Taniguchi, Takashi and Watanabe, Kenji and Zar. Bilayer Wigner Crystals in a Transition Metal Dichalcogenide Heterostructure , journaltitle =. 2021 , volume =

    work page 2021

  11. [11]

    Signatures of Wigner Crystal of Electrons in a Monolayer Semiconductor , journaltitle =

    Smole. Signatures of Wigner Crystal of Electrons in a Monolayer Semiconductor , journaltitle =. 2021 , volume =

    work page 2021

  12. [12]

    and Mak, Kin Fai and Shan, Jie , title =

    Tang, Yanhao and Gu, Jie and Liu, Song and Watanabe, Kenji and Taniguchi, Takashi and Hone, James C. and Mak, Kin Fai and Shan, Jie , title =. 2022 , volume =

    work page 2022

  13. [13]

    Keldysh, L. V. , title =. 1979 , volume =

    work page 1979

  14. [14]

    and Szyniszewski, M

    Mostaani, E. and Szyniszewski, M. and Price, C. H. and Maezono, R. and Danovich, M. and Hunt, R. J. and Drummond, N. D. and Fal'ko, V. I. , title =. 2017 , volume =

    work page 2017

  15. [15]

    and Ruiz-Tijerina, D

    Danovich, M. and Ruiz-Tijerina, D. A. and Hunt, R. J. and Szyniszewski, M. and Drummond, N. D. and Fal'ko, V. I. , title =. 2018 , volume =

    work page 2018

  16. [16]

    , title =

    Lindhard, J. , title =. 1954 , volume =

    work page 1954

  17. [17]

    1967 , volume =

    Stern, Frank , title =. 1967 , volume =

    work page 1967

  18. [18]

    , title =

    Adler, Stephen L. , title =. 1962 , volume =

    work page 1962

  19. [19]

    1963 , volume =

    Wiser, Nathan , title =. 1963 , volume =

    work page 1963

  20. [20]

    Monarkha, Yu. P. and Syvokon, V. E. , title =. 2012 , volume =

    work page 2012

  21. [21]

    and Pfeiffer, L

    Deng, H. and Pfeiffer, L. N. and West, K. W. and Baldwin, K. W. and Engel, L. W. and Shayegan, M. , title =. 2019 , volume =

    work page 2019

  22. [22]

    and Pfeiffer, Loren N

    Zhao, Lili and Lin, Wenlu and Chung, Yoon Jang and Gupta, Adbhut and Baldwin, Kirk W. and Pfeiffer, Loren N. and Liu, Yang , title =. 2023 , volume =

    work page 2023

  23. [23]

    and Park, S

    Robert, C. and Park, S. and Cadiz, F. and Lombez, L. and Ren, L. and Tornatzky, H. and Rowe, A. and Paget, D. and Sirotti, F. and Yang, M. and Van Tuan, D. and Taniguchi, T. and Urbaszek, B. and Watanabe, K. and Amand, T. and Dery, H. and Marie, X. , title =. Nature Communications , year =

    work page

  24. [24]

    and Sidler, M

    Srivastava, A. and Sidler, M. and Allain, A. V. and Lembke, D. S. and Kis, A. and Imamoglu, A. , title =. Nature Physics , year =

    work page

  25. [25]

    and Yang, L

    Dey, P. and Yang, L. and Robert, C. and Wang, G. and Urbaszek, B. and Marie, X. and Crooker, S. A. , title =. Physical Review Letters , year =

    work page

  26. [26]

    and Goryca, M

    Li, J. and Goryca, M. and Yumigeta, K. and Li, H. and Tongay, S. and Crooker, S. A. , title =. Physical Review Materials , year =

    work page

  27. [27]

    and Sinitsyn, N

    Yang, L. and Sinitsyn, N. A. and Chen, W. and Yuan, J. and Zhang, J. and Lou, J. and Crooker, S. A. , title =. Nature Physics , year =

    work page

  28. [28]

    Signatures of Wigner crystal of electrons in a monolayer semiconductor , journal =

    Smole. Signatures of Wigner crystal of electrons in a monolayer semiconductor , journal =. 2021 , volume =

    work page 2021

  29. [29]

    and Menzel, F

    Wang, L. and Menzel, F. and Pichler, F. and Kn. Spectroscopy of Wigner crystal polarons in an atomically thin semiconductor , journal =. 2025 , eprint =

    work page 2025

  30. [30]

    and Dizer, Eugen and Schmidt, Richard and Imamoglu, Atac and Christianen, Arthur , title =

    Adlong, Haydn S. and Dizer, Eugen and Schmidt, Richard and Imamoglu, Atac and Christianen, Arthur , title =. arXiv preprint , year =. 2512.16651 , eprinttype =

    work page arXiv

  31. [31]

    and Chang, Yia-Chung and Lui, Chun Hung , title =

    Liu, Erfu and Wilson, Matthew and Hu, Jenny and Zimmerman, Alexandra and Mathew, Amal and Ouyang, Tianyi and Shi, Ao and Taniguchi, Takashi and Watanabe, Kenji and Heinz, Tony F. and Chang, Yia-Chung and Lui, Chun Hung , title =. arXiv preprint , year =. 2601.11914 , eprinttype =

    work page arXiv

  32. [32]

    Optical probing of Wigner crystallization in monolayer WSe$_2$ via diffraction of longitudinal excitons

    Abramov, Artem N. and Chiglintsev, Emil and Oskolkova, Tatiana and Titova, Maria and Kashchenko, Mikhail and Bandurin, Denis and Chernov, Alexander and Kravtsov, Vasily and Iorsh, Ivan V. , title =. arXiv preprint , year =. 2601.16319 , eprinttype =

    work page internal anchor Pith review Pith/arXiv arXiv