{"id":"c8985456-d2e7-4211-9916-66b63ec6f05b","arxiv_id":"2607.05371","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":4.0,"correctness_risk":"unknown","formal_verification":"none","parameter_count":7,"one_line_summary":"A self-consistent Hartree-Fock calculation yields the band structure and static dielectric response of a spin-polarized Wigner crystal in monolayer WSe2, showing reduced screening at low carrier density.","lead":"The paper computes the band structure and static dielectric function of a fully spin-polarized Wigner crystal in monolayer WSe2 using a self-consistent Hartree-Fock method with a Keldysh interaction. It provides a theoretical tool for screening and interaction effects in low-density 2D TMD systems, relevant to recent optical experiments.","discovery_kind":"unclear","skeptic_critique":{"model":"glm-5.2","headline":"No quantitative convergence data or solution-uniqueness checks are presented for the reported εHF(q,0), M1, and Cn values, despite acknowledged multiple numerical attractors and a free initial-state parameter β0.","rationale":"The reader correctly identified convergence data as one of several gaps and gave a CONDITIONAL verdict. I agree with that verdict. However, I would sharpen the focus: the reader's weakest_assumption centered on the HF approximation itself (missing correlations, full spin polarization, vertex corrections). While those are real and known limitations of any HF calculation — and the paper is appropriately honest about them — they are generic to the method class rather than specific to this work. The more immediately load-bearing concern is the absence of quantitative convergence evidence for the specific numbers reported. The paper acknowledges multiple numerical attractors and uses a free initial-state parameter β0, yet provides no data showing that the converged solution is unique or that the reported quantities are stable under basis-set refinement. This is the single check that would most directly settle whether the quantitative results hold. The qualitative physics (stronger modulation at larger rs, weaker screening at larger rs) follows from basic considerations and is likely robust regardless. The HF-level approximation limits the quantitative accuracy of the absolute values but does not undermine the internal consistency of the framework or the qualitative trends. The paper's analytical derivations (long-wavelength limit, selection rules, density reconstruction) are correct and provide internal consistency checks. The framework is sound; the numbers need validation. CONDITIONAL remains the appropriate verdict.","tokens_in":17604,"tokens_out":2755,"duration_ms":263799,"concrete_test":"Recompute εHF(q,0) and M1 at rs=30 with: (a) Nside = 10 vs. 20, (b) NG = 50 vs. 120, (c) Nb = 13 vs. 17, (d) two different β0 values (e.g., β0 = 0.5 and 2.0). Report the peak value of εHF and M1 for each. If either quantity shifts by more than ~10% across any pair, the reported values are not converged and the quantitative claims need revision. If they are stable to within a few percent, the numerical results are trustworthy and the CONDITIONAL verdict could be relaxed.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The paper's central quantitative output consists of specific values of M1, Cn, and εHF(q,0) at rs = 4, 10, 30, 50 (Fig. 2, Sec. 3.4–3.5). Section 3.1 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 ηΠ,' but no convergence table, plot, or quantitative statement about how much these quantities shift with parameter variation is provided anywhere in the main text or SI. Table S1 lists parameter ranges (Nside = 5–20, NG = 50–120, Nb = 13–17) but does not report which values were used for which results, nor how the outputs depend on them. This gap is load-bearing for three reasons: (1) The self-consistency procedure explicitly acknowledges multiple Hartree-preconditioned attractors (Sec. 3.1), and branch selection uses qualitative criteria (symmetry preservation, band regularity, degeneracy pattern) that introduce unquantified systematic uncertainty. (2) The initial Gaussian width β0 is a free parameter; the paper shows M1 changes from initial to converged values (e.g., rs=30: Mini=0.668 → Mfin=0.504), demonstrating the solution is not trivially the initial guess, but does not show whether different β0 values converge to the same solution. (3) The dielectric function depends on the number of empty bands Nb and broadening ηΠ in the polarizability sum (Eq. 14); without convergence data, the peak values of εHF(q,0) — the paper's main quantitative result — cannot be assessed for reliability. The qualitative trends (M1 increases with rs, ε peak decreases with rs) are physically expected and likely robust. But the specific numerical values, which are what the paper presents as results, are unverified.","agreement_with_reader":"partial"},"referee_report":{"model":"glm-5.2","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.","tokens_in":17869,"tokens_out":990,"duration_ms":239040,"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":[{"comment":"§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值","section":null}],"minor_comments":[{"comment":"§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.","section":null},{"comment":"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.","section":null},{"comment":"§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.","section":null},{"comment":"References: Refs 11 and 23 appear to be duplicate citations of the same arXiv preprint (arXiv:2512.16631). This should be consolidated.","section":null},{"comment":"The code is stated to be available 'upon reasonable request.' For reproducibility, depositing the code in a public repository would be preferable.","section":null}],"recommendation":"major_revision","confidential_remarks":"The reader's assessment and the stress-test note both correctly identify the lack of convergence data as the primary weakness. This is a genuine and load-bearing gap: the paper reports specific numerical values (M1, Cn, epsilon_HF peaks) at four rs values without any quantitative justification that these values are converged. The fix is straightforward (add a convergence table), which is why I recommend major revision rather than rejection. The theoretical framework itself is sound and standard; the issue is purely one of numerical verification."},"author_rebuttal":{"model":"glm-5.2","summary":"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.","responses":[{"response":"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_made":"yes","referee_comment":"§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."}],"tokens_in":17087,"tokens_out":264,"duration_ms":34406,"standing_objections":[]},"desk_editor":{"model":"glm-5.2","letter":"Short version: this is a legitimate application of self-consistent Hartree-Fock to a Wigner crystal in monolayer WSe2 with a Keldysh interaction, producing a static dielectric function. The framework is standard and correctly implemented. The main gap is that the paper claims convergence checks but shows no quantitative convergence data, so the specific numerical values of εHF(q,0), M1, and Cn cannot be independently assessed. Deserves a serious referee who asks for that evidence.","headline":"Sound HF Wigner-crystal framework applied to WSe2 with Keldysh interaction; numerical results need convergence evidence","tokens_in":18503,"tokens_out":549,"would_cite":false,"duration_ms":62175,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"glm-5.2","headline":"Self-consistent Hartree-Fock maps WSe2 Wigner crystal density, bands, and screening","keywords":["Wigner crystal","Hartree-Fock","dielectric function","WSe2","Keldysh interaction","monolayer TMD","Wigner crystal","static screening"],"falsifier":"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.","tokens_in":17639,"feed_emoji":"🔬","tokens_out":812,"duration_ms":142017,"temperature":0.7,"pith_summary":"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.","feed_headline":"Hartree-Fock calculation maps WSe2 Wigner crystal bands and dielectric response","feed_subtitle":"A self-consistent method links crystalline density modulation to reduced electronic screening in monolayer WSe2, with testable optical-signa","key_machinery":"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.","core_discovery":"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屏","pith_inferences":[],"forward_implications":[],"fun_headline_variants":["Hartree-Fock method maps WSe2 Wigner crystal dielectric function","Stronger Wigner localization reduces WSe2 dielectric screening","Increasing electron localization lowers WSe2 dielectric screening","Self-consistent framework links WSe2 Wigner crystal bands to screening"],"cache_read_input_tokens":0,"weakest_assumption_plain":"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.","fun_headline_variants_meta":{"raw":{"variants":["Hartree-Fock method maps WSe2 Wigner crystal dielectric function","Stronger Wigner localization reduces WSe2 dielectric screening","Increasing electron localization lowers WSe2 dielectric screening","Self-consistent framework links WSe2 Wigner crystal bands to screening"]},"model":"glm-5.2","effort":"high","cost_usd":0.0,"raw_usage":{"total_tokens":1083,"prompt_tokens":501,"completion_tokens":582,"prompt_tokens_details":null},"tokens_in":501,"tokens_out":582,"duration_ms":62267,"temperature":1.0,"reasoning_tokens":540,"cache_read_input_tokens":0,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-07-07T14:11:17.469805+00:00","model_set":{"reader":"glm-5.2"},"falsifier":"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.","supporting_citations":[],"review_version":1}