{"id":"3fd12d93-2cbc-4fab-a81b-d993107a5ede","arxiv_id":"2505.24706","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":0,"one_line_summary":"For three-dimensional interacting fermions at zero temperature, the ground-state Wigner function approaches Thomas-Fermi theory in trace norm at rate N times a positive power of the semi-classical parameter, with Coulomb interactions allowed.","lead":"This paper proves a quantitative version of the semi-classical limit for large three-dimensional Fermi gases at zero temperature, with an explicit convergence rate of the ground-state Wigner function to the Thomas-Fermi distribution. It is the first rate stated for singular pair potentials, including Coulomb repulsion.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Uniformity of the Weyl-law input (Remark 3.5) is load-bearing and only sketched: if the O(ℏ^{-2}) constant in Proposition 3.2 grows with N, the bootstrap (5.23)–(5.25) and Theorem 1 fail.","rationale":"The reader's weakest_assumption identifies the import of Proposition 3.1 from [8] as the primary concern and mentions the uniform Weyl-law refinement in Remark 3.5 as a secondary input. I agree that both are load-bearing, but I elevate the uniformity of the Weyl law to the single most load-bearing concern. The reason is impact: the Weyl-law constant enters directly into the bootstrap (5.23)–(5.25) through (5.24). If that constant is not uniform in N, the bound (5.2) would contain an extra N-dependent factor, and the final trace-norm estimate (2.2) would lose its N^{-δ/3} smallness. By contrast, a wrong exponent in Proposition 3.1(3) would change δ but would not necessarily destroy the qualitative convergence. Furthermore, Proposition 3.1 is proved in the companion preprint [8] (modulo a scaling translation that Remark 3.4 explains), whereas the uniformity of Proposition 3.2 is only sketched in Remark 3.5 and is not fully verified in the manuscript. The paper is honest about this: it explicitly says the uniformity is not stated in [21] and gives a brief justification. That justification is plausible but is precisely the kind of omitted proof that should be checked before the theorem is accepted with the stated rate. I still recommend the same conditional verdict as the reader, because the concern is a gap in verification, not a demonstrated error. The proposed check would settle it by a careful constant-tracing through Mikkelsen's proof, or by a numerical scaling test on the actual family W_N.","tokens_in":23794,"tokens_out":21824,"duration_ms":231820,"concrete_test":"Trace the proof of [21, Theorem 1.5] as used in Proposition 3.2 and list every quantity controlling the O(ℏ^{-2}) constant. Verify that each depends only on R = ||W||_{C^{1,α}(Ω(W,ν))} and ν, and not on the diameter of Ω(W,ν), on ||W||_{C^{1,α}} outside Ω, or on any N-dependent localization scale. If any such hidden dependence appears, construct the family W_N from the Hartree minimizers and numerically evaluate Tr 1_{(−∞,0]}(−ℏ^2Δ + W_N) at several N to see whether the error relative to the Weyl term is uniformly O(ℏ^{-2}) with a fixed constant; a growing error would invalidate the bootstrap and Theorem 1.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central theorem's rate is controlled by the number estimate (5.2), whose proof bootstraps around the term dΓ(1_{|H_γ−μ|<ε}) in (5.24). This term is bounded by Tr 1_{|H_γ−μ|<ε} ≤ C N ε, which is exactly Corollary 3.1. Corollary 3.1 in turn rests on Proposition 3.2, the optimal Weyl law with error O(ℏ^{-2}) uniform over the family W_N = U + V ∗ ϱ_{γ_H} − μ. The uniformity is asserted only in Remark 3.5 and is not proved in the present paper. If the constant in (3.24) were C_N growing with N, then (5.24) would become C_N N ε, and with ε = CℏΛ^{1+a} the final bound would contain C_N Nℏ^{δ_1}, which is only small if C_N ≪ ℏ^{-δ_1}. No such control is established. The sketch in Remark 3.5 lists dependence on the C^{1,α} seminorms over Ω(W,ν) = {W < ν}, but it does not rule out hidden dependence on the diameter of Ω(W,ν), on global growth of W, or on the localization scale in the multi-scale argument. Because this input enters at the crucial bootstrap step, its failure would not just shift δ but would destroy the convergence statement. The companion-paper input Proposition 3.1 is also essential, but it has a full proof elsewhere; the uniformity of the Weyl law is an internal gap, explicitly flagged as a limitation only in a remark.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"This paper proves a quantitative semi-classical limit for the one-particle reduced density matrix of approximate ground states of N spinless fermions in R^3 with Hamiltonian H_N = Σ(-ℏ²Δ_{x_i}+U(x_i)) + N^{-1}Σ_{i<j}V(x_i-x_j), ℏ=N^{-1/3}, under Condition 1 on U and V=λ|x|^{-a}, a∈(0,1]. The main result, Theorem 1, gives ||γ_{Ψ_N}-γ_TF||_Tr ≤ C N ℏ^δ |lnℏ|^{1/4} with δ = (1/2) min((6-5a)/(16+5a),(4+15a)/(2(16+5a))) for any approximate ground state with energy error O(N^{-1/6}). The proof compares the N-body state to the Hartree minimizer via a particle-hole transformation, controls the fluctuation number operator using a regularized potential and an optimal Weyl law, and then invokes the Hartree-to-Thomas-Fermi trace estimate from a companion paper [8].","tokens_in":24134,"tokens_out":21520,"duration_ms":259795,"significance":"If correct, this is a significant advance: it provides the first explicit polynomial-in-ℏ rate for convergence of states in the combined mean-field/semiclassical limit for interacting fermions with singular potentials including Coulomb, in the trace-norm topology, going beyond the nonquantitative compactness arguments of previous works. The paper contains a detailed and mostly self-contained treatment of the new operator estimates (Section 4 and Appendix A), including a complete particle-hole conjugation computation, and the final rate is explicitly tracked through the regularization and bootstrap. The main caveats are the reliance on two external inputs: the Hartree-minimizer estimates of [8] and the uniform Weyl law of [21], the latter of which is only sketched.","major_comments":[{"comment":"The constant in the optimal Weyl law is required to be uniform over the N-dependent family W_N = U + V*ϱ_{γ_N^H} - μ, because (5.24) bounds Tr 1_{|H_{γ_N^H}-μ|<ε} ≤ C N ε and this enters the bootstrap (5.23)-(5.27). The paper only sketches this uniformity in Remark 3.5 and does not prove it or state it as an explicit assumption. If the constant in (3.24) grows with N, then (5.24) would carry a factor C_N, and the final bound (5.30) would contain C_N N ℏ^{δ_1}, which is not small without control of C_N; the proof of Theorem 1 would fail. Please either supply a complete proof of the uniformity for the family W_N or declare this uniformity as an explicit hypothesis of the theorem.","section":"§3.2 (Prop. 3.2, Remark 3.5) and §5, Eq. (5.24)"},{"comment":"The estimates (3.18)-(3.21) are imported from the companion preprint [8] and are load-bearing: (3.18) provides the commutator estimate used in Lemmas 4.3-4.4, (3.19) provides the uniform L^p bounds on ϱ_{γ_N^H} used in the regularization step, and (3.20) is the final Hartree-to-Thomas-Fermi closeness. These results are not re-derived here, and Remark 3.4 only asserts that the scaling changes are superficial. Since [8] is a preprint, the present theorem is conditional on its correctness. The manuscript should state this dependence explicitly (e.g., as a standing assumption) or include the needed statements with proofs.","section":"§3.2, Proposition 3.1 and Remark 3.4"}],"minor_comments":[{"comment":"The logarithmic exponent for the L^2 Wigner convergence is stated as |ln ℏ|^{1/4}, but the derivation in Remark 2.1 gives |ln ℏ|^{1/8}; one of the two should be corrected.","section":"§1, Eq. (1.13)"},{"comment":"The phase-space integral is written over R^{2d} although the statement is for d=3; this should be R^6 or the notation should be defined consistently with d.","section":"§3.2, Eq. (3.24)"},{"comment":"The sentence 'we include the correction due to the exchange term, i.e. which is at most of order Nℏ^{1/2}' is terse; a short estimate for (1/2N)Tr X_γ γ would improve readability.","section":"§5, Proof of Lemma 5.2"},{"comment":"Condition 1 is stated for general d while Theorem 1 is d=3; making d=3 explicit already in Condition 1 would avoid confusion.","section":"§2.1, Condition 1"}],"recommendation":"major_revision","confidential_remarks":"The main risk to the paper is the unproved uniformity in the Weyl law input. If the author can supply a verification (or the companion paper [21] contains it in the stated strength), I would be willing to accept the result. I also recommend that the editor ask the author to clarify the status of [8] and [21] (both preprints) relative to the present theorem."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Bottom line: this is the first quantitative trace-norm rate for the N-body fermionic ground state to Thomas-Fermi, and the proof is mostly sound. The central theorem is genuinely new: previous convergence results were qualitative, and no rates existed for Coulomb. The adaptation of second-quantized number estimates to the stationary problem, including the particle-hole diagonal terms, is substantial and clearly written. The paper also states its limits honestly: delta is not optimal, higher marginals are left open, and the dependence on the companion paper is flagged.\n\nThe soft spots are the two external legs. Most of the Hartree-state input comes from the companion [8]; that is not a flaw per se, since it is proved elsewhere, but it makes the theorem conditional on those estimates. The sharper issue is the uniformity of the Weyl law used in Corollary 3.1. The paper invokes Mikkelsen [21] and defends uniformity only in Remark 3.5. That uniformity is load-bearing: if the constant in (3.24) grows with N, the bootstrap in (5.23)-(5.25) collapses and Theorem 1 has no rate. I read the stress-test concern as legitimate but containable. The potentials W_N = U + V*rho_N - mu have sublevel sets that are uniformly bounded because V >= 0 and U -> infinity, and the needed C^{1,alpha} seminorms are uniform by [8]. So the missing piece is a verification, not a contradiction; it should be added or outsourced explicitly.\n\nThe rest of the proof holds together. Lemma 4.4 for diagonal terms is new and checked in the appendix. The trace-norm estimate (6.2) is standard but correct given Theorem 2. I found no circularity: the Hartree-minimizer input is parameter-free, and the final rate is not used to set constants.\n\nWho is this for: anyone working on semiclassical limits of Fermi systems or quantitative mean-field/Thomas-Fermi comparisons. It deserves serious refereeing; an editor should send it out, with the request that the Weyl uniformity be made fully precise or replaced by a self-contained statement.","headline":"First quantitative trace-norm rate from an N-body Fermi ground state to Thomas-Fermi, including Coulomb; proof is careful, but a load-bearing uniformity claim in the Weyl-law input is only sketched.","tokens_in":24643,"tokens_out":2557,"would_cite":true,"duration_ms":30686,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":["81Q20","81V70","35P20"],"pacs":[],"model":"deepseek-v4-flash","headline":"This paper proves that in three dimensions, an approximate ground state of $N$ interacting fermions has one-particle density matrix within trace distance $O(N \\hbar^\\delta |\\ln \\hbar|^{1/4})$ of Thomas-Fermi theory, where $\\hbar=N^{-1/3}$…","keywords":["semiclassical limit","Thomas-Fermi theory","fermionic ground state","trace norm convergence","Hartree minimizer","commutator estimates","Weyl law","Coulomb potential"],"falsifier":"A numerical computation of $\\|\\gamma_H^N-\\gamma_{\\mathrm{TF}}\\|_{\\mathrm{Tr}}$ for the Coulomb case $a=1$ in $d=3$ at increasing $N$ would settle the rate; if it grows faster than $N\\hbar^{1/2}|\\ln\\hbar|^{1/2}$, the power $\\delta$ in Theorem 1 is not correct.","tokens_in":23557,"feed_emoji":"⚛️","tokens_out":11500,"duration_ms":124527,"temperature":0.7,"pith_summary":"The paper aims to prove that large systems of interacting fermions at zero temperature approach Thomas-Fermi theory with a controllable error, not just in the limit. In three dimensions, and for pair potentials as singular as the Coulomb repulsion, it shows that the one-particle density matrix of any approximate ground state converges in trace norm to the Weyl quantization of the Thomas-Fermi phase-space density at an explicit power of the semiclassical parameter $\\hbar=N^{-1/3}$, up to a logarithmic correction. This matters because a quantitative rate converts a qualitative compactness result into error estimates that can be used in further asymptotic analysis, for instance in atomic and molecular energy expansions. The proof identifies the Hartree minimizer as the intermediate quantum state and controls the fluctuation number operator at the same explicit rate.","feed_headline":"Explicit rate: fermions converge to Thomas-Fermi","feed_subtitle":"Trace distance between N-body ground state and Thomas-Fermi drops at an explicit power of ħ with log correction, Coulomb included.","key_machinery":"The argument routes through the Hartree minimizer $\\gamma_H^N$, the one-particle density matrix that minimizes the Hartree functional, as an intermediate state. It uses a particle-hole transformation $R$ on the fermionic Fock space (the Hilbert space containing all particle-number sectors) that re-expresses the $N$-body state as fluctuations around $\\gamma_H^N$; the key object is the number estimate $\\langle \\Omega_N, \\mathcal{N} \\Omega_N \\rangle \\le C N \\hbar^{2\\delta}|\\ln \\hbar|^{1/2}$ for the fluctuation vector $\\Omega_N=R^*\\Psi_N$. The bound is forced through an operator inequality $\\mathcal{N} \\le \\varepsilon^{-1} d\\Gamma(|H_\\gamma-\\mu|)+d\\Gamma(1_{|H_\\gamma-\\mu|<\\varepsilon})$ on Fock space, where $H_\\gamma$ is the Hartree Hamiltonian, and it combines three imported ingredients: semiclassical commutator estimates for $\\gamma_H^N$, an ultraviolet regularization of the singular potential that preserves non-negativity of its Fourier transform, and a sharp Weyl law for Schr\\\"odinger operators with non-smooth potentials that controls the number of eigenvalues in a spectral window.","core_discovery":"The paper's central claim is that the semiclassical limit of large fermionic ground states is quantitative. For three dimensions, if $(U,V)$ satisfy Condition 1---$U$ a confining potential with bounded weighted Hessian and $V(x)=\\lambda |x|^{-a}$ with $0<a\\le 1$, including Coulomb $a=1$---then for every approximate ground state $\\Psi_N$ with energy error $\\varepsilon_N=O(N^{-1/6})$, the one-particle density matrix $\\gamma_{\\Psi_N}$ satisfies $\\|\\gamma_{\\Psi_N}-\\gamma_{\\mathrm{TF}}\\|_{\\mathrm{Tr}} \\le C N \\hbar^\\delta |\\ln \\hbar|^{1/4}$, with $\\hbar=N^{-1/3}$ and $\\delta = \\tfrac12 \\min\\left(\\frac{6-5a}{16+5a}, \\frac{4+15a}{2(16+5a)}\\right)$. Here $\\gamma_{\\mathrm{TF}}$ is the Weyl quantization of the Thomas-Fermi Wigner function $f_{\\mathrm{TF}}(x,p)=1_{|p|^2\\le C_{\\mathrm{TF}}\\rho_{\\mathrm{TF}}(x)^{2/3}}$, the generalized Fermi ball. The result upgrades previously qualitative convergence of states to an explicit rate in the semiclassical parameter.","pith_inferences":["Editorial extension: the same second-quantized strategy should yield a temperature-dependent quantitative limit at positive temperature by replacing the zero-temperature Fermi projection with a Fermi-Dirac density matrix, with the rate degrading in temperature.","Editorial extension: the exponent $\\delta$ is a byproduct of balancing regularization and spectral-gap errors, not a fundamental limit, so a direct spectral-gap argument could plausibly improve the rate for regular potentials.","Editorial extension: the trace-norm control of the ground state can serve as an input for time-dependent problems, potentially giving quantitative propagation of the semiclassical structure under Hartree-Fock dynamics."],"forward_implications":["The Wigner function of any approximate ground state converges in $L^2$ at rate $O(\\hbar^{\\delta/2}|\\ln \\hbar|^{1/8})$.","Position and momentum densities converge in $L^1$ at the trace rate $O(\\hbar^\\delta |\\ln \\hbar|^{1/4})$.","For regular pair potentials satisfying the non-negativity and integrability condition (2.6), the rate improves to $\\delta=1/2$.","For $a<1$ the logarithmic factor disappears from the bound.","The calculation formally extends to super-Coulombic potentials with $1<a<6/5$ once the required Hartree-minimizer estimates become available."],"supporting_citations":[{"why":"Supplies the Hartree-minimizer commutator estimates and the $O(N\\hbar^{1/2}|\\ln\\hbar|^{1/2})$ trace closeness to Thomas-Fermi on which the whole quantitative chain rests.","marker":"[8]"},{"why":"Provides the particle-hole transformation and the second-quantized number-estimate strategy adapted here.","marker":"[2]"},{"why":"Establishes the qualitative semiclassical limit of states that this paper makes quantitative.","marker":"[11]"},{"why":"Supplies the sharp Weyl law for non-smooth potentials used to control eigenvalue-count traces in the number estimate.","marker":"[21]"},{"why":"Guarantees the existence and fixed-point structure of the Hartree minimizers used as the intermediate state.","marker":"[22]"},{"why":"Provides the original commutator bounds for spectral projections whose Hartree-state extension is imported from [8].","marker":"[20]"}],"fun_headline_variants":["Fermions get explicit quantum-to-classical rate","Quantitative semiclassical limit for Fermi gas","Fermion ground states: explicit convergence to Thomas-Fermi","Coulomb included: explicit semiclassical rate for fermions","N-body fermions: rate of convergence to Thomas-Fermi"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The proof rests on imported estimates on the Hartree minimizer---that its commutators with position and momentum are $O(N\\hbar)$ and that its trace distance to Thomas-Fermi is $O(N\\hbar^{1/2}|\\ln\\hbar|^{1/2})$---rather than re-deriving them here.","fun_headline_variants_meta":{"raw":{"variants":["Fermions get explicit quantum-to-classical rate","Quantitative semiclassical limit for Fermi gas","Fermion ground states: explicit convergence to Thomas-Fermi","Coulomb included: explicit semiclassical rate for fermions","N-body fermions: rate of convergence to Thomas-Fermi"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000671,"raw_usage":{"total_tokens":3038,"prompt_tokens":908,"completion_tokens":2130,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":524,"completion_tokens_details":{"reasoning_tokens":2047}},"tokens_in":524,"tokens_out":2130,"duration_ms":20642,"temperature":1.0,"reasoning_tokens":2047,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-07T12:16:25.913653+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"A numerical computation of $\\|\\gamma_H^N-\\gamma_{\\mathrm{TF}}\\|_{\\mathrm{Tr}}$ for the Coulomb case $a=1$ in $d=3$ at increasing $N$ would settle the rate; if it grows faster than $N\\hbar^{1/2}|\\ln\\hbar|^{1/2}$, the power $\\delta$ in Theorem 1 is not correct.","supporting_citations":[{"cited_title":"Benedikter, M","cited_arxiv_id":null,"evidence_quote":"Provides the particle-hole transformation and the second-quantized number-estimate strategy adapted here."},{"cited_title":"Fournais, M","cited_arxiv_id":null,"evidence_quote":"Establishes the qualitative semiclassical limit of states that this paper makes quantitative."},{"cited_title":"Mikkelsen","cited_arxiv_id":null,"evidence_quote":"Supplies the sharp Weyl law for non-smooth potentials used to control eigenvalue-count traces in the number estimate."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Guarantees the existence and fixed-point structure of the Hartree minimizers used as the intermediate state."},{"cited_title":"Fournais, S","cited_arxiv_id":null,"evidence_quote":"Provides the original commutator bounds for spectral projections whose Hartree-state extension is imported from [8]."}],"review_version":1}