The quantitative semi-classical limit of a large Fermi system at zero temperature

Esteban C\'ardenas

REVIEW 2 major objections 4 minor 2 cited by

The quantitative semi-classical limit of a large Fermi system at zero temperature

T0 review · 2 major / 4 minor · reviewed 2026-08-07 · deepseek-v4-flash

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

Pith's one-line read 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}$…

desk verdict 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. read the letter →

arxiv 2505.24706 v1 pith:IKFV5J46 submitted 2025-05-30 math-ph math.APmath.MP

Esteban Cárdenas This is my paper

classification math-phmath.APmath.MP MSC 81Q2081V7035P20

keywords semiclassicallimitThomas-FermitheoryfermionicgroundstatetracenormconvergenceHartreeminimizercommutatorestimatesWeyllawCoulombpotential

verification ladder T0 review T1 audit T2 compute T3 formal

The pith

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

The reading

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.

What carries the argument

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.

What would settle it

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.

Watch

Extended reading notes

Core claim

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.

Load-bearing premise

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.

Editorial extensions

If this is right

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.

Reading between the lines

Editorial extensions of the paper, not claims the author makes directly.

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.

Share X Bluesky LinkedIn Reddit HN

Editorial analysis

A structured set of objections, weighed in public.

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

Desk Editor's Note

read the 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.

The 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.

The 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.

Who 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.

Referee Report

2 major / 4 minor

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].

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 (2)

[§3.2 (Prop. 3.2, Remark 3.5) and §5, Eq. (5.24)] 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.
[§3.2, Proposition 3.1 and Remark 3.4] 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.

minor comments (4)

[§1, Eq. (1.13)] 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.
[§3.2, Eq. (3.24)] 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.
[§5, Proof of Lemma 5.2] 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.
[§2.1, Condition 1] Condition 1 is stated for general d while Theorem 1 is d=3; making d=3 explicit already in Condition 1 would avoid confusion.

Circularity Check

0 steps flagged · score 0.0 of 10

No circularity: the N-body-to-Hartree step is proved in the paper from explicit Fock-space bounds, and the Hartree-minimizer input from [8] is an independent derivation whose assumptions do not include the target trace-norm result.

full rationale

The derivation chain is not circular. Theorem 1 is obtained by the triangle inequality ||γ_ΨN − γ_TF||_Tr ≤ ||γ_ΨN − γ_HN||_Tr + ||γ_HN − γ_TF||_Tr. The first term is proved in the present paper via Theorem 2 and the trace argument in Section 6 (Eq. (6.2)); the second term is imported from [8, Theorem 1.6] as Proposition 3.1(3), which proves the Hartree-to-Thomas–Fermi closeness for the auxiliary Hartree minimizer without assuming the N-body convergence. The crucial number estimate (5.2) is derived here from the operator inequality (2.8), the commutator estimates from Proposition 3.1, and the sharp Weyl law of Mikkelsen [21] through Corollary 3.1; no quantity in Theorem 1 is fed back into these inputs. The paper explicitly flags its reliance on [8] in Remark 3.4 and Section 3.2, but that reliance is on independently proved estimates for γ_HN, not on a result equivalent to the paper's conclusion. The uniformity assertion in Remark 3.5 is a genuine internal gap—the uniformity of the O(ℏ^{-2}) constant in Proposition 3.2 is stated as a technical corollary of [21] rather than proved here—but an asserted uniformity of an external Weyl-law result is not circular; it is a completeness or robustness risk. There are no fitted parameters renamed as predictions, no equation is defined in terms of the target quantity, and the rate δ is an explicit function of a and ℏ, not a constant chosen to match the final bound. Hence the honest finding is no significant circularity.

Assumptions & free parameters 0 free parameters · 5 assumptions · 0 invented entities

No fitted parameters appear; δ is a derived exponent obtained by optimizing two auxiliary cutoffs. The central claim rests on imported quantitative estimates for Hartree minimizers and on a uniform Weyl-law corollary, both flagged in the text.

assumptions (5)

domain assumption Condition 1: U ∈ C^2_loc with exponential decay of D^2U and V(x) = λ|x|^{-a}, λ > 0, a ∈ (0,1] (Coulomb included).
The theorem is stated only under these regularity and growth conditions; the rate δ depends on a. Introduced in Condition 1 and used in Theorem 1.
domain assumption The state Ψ_N is an approximate ground state with ε_N = O(N^{-1/6}).
The central estimate is for states within this energy error; the final rate is tied to this choice. See Eq. (1.5) and Theorem 1.
standard math Hartree minimizers satisfy commutator estimates (3.18), uniform L^p density bounds (3.19), and trace bound (3.20) from [8].
Imported as black-box theorems from the companion preprint [8], summarized in Proposition 3.1 and used throughout Sections 5 and 6.
ad hoc to paper The optimal Weyl law (2.9)/(3.24) holds uniformly for the family of Hartree Hamiltonians with C^{1,α} potentials.
Mikkelsen's published theorem does not state the needed uniformity; the paper asserts it is a technical corollary in Remark 3.5. Used in Corollary 3.1 and Eq. (5.24) to control the spectral window term.
standard math Lieb-Thirring inequalities, CLR bounds, Agmon estimates, and Sobolev embeddings used in Lemma 5.1 and Proposition 3.1.
Standard tools stated in Sections 3 and 5; needed for the tail estimate (5.8) and for uniform density bounds.

how reviews work

0 comments

Cite this review

Pith. "Pith review of The quantitative semi-classical limit of a large Fermi system at zero temperature." pith.science (2026). https://pith.science/paper/IKFV5J46

@misc{pith2026250524706,
  author       = {Pith},
  title        = {Pith review of: The quantitative semi-classical limit of a large Fermi system at zero temperature},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/IKFV5J46}},
  note         = {Machine review of arXiv:2505.24706}
}

read the original abstract

In this article we consider a large system of fermions in a combined mean-field and semiclassical limit, in three dimensions. We investigate the convergence of the Wigner function of the ground state, towards the classical Thomas-Fermi theory. The main novelty of the present article is quantifying the convergence rate with respect to the semi-classical parameter. One of the main ingredients is a recent result on the validity of semi-classical commutator estimates satisfied by the Hartree theory. Singular potentials, up to the Coulomb interaction, are included.

Discussion (0). Sign in to comment.

Forward citations

Cited by 2 Pith papers

Reviewed papers in the Pith corpus that reference this work. Sorted by Pith novelty score. Full citation record

Derivation of the time-dependent Hartree equations for strongly interacting dense fermionic systems
math-ph 2025-07 conditional novelty 8.0 of 10

The rescaled Schrödinger dynamics of N dense fermions with C^2 pair interactions converges to the time-dependent Hartree dynamics with explicit rate N^{-1/24}.
Semi-classical limit of an attractive Fermi gas in one or two dimensions
math-ph 2026-02 conditional novelty 6.0 of 10

For trapped attractive Fermi gases in 1D and 2D, as N grows the ground-state energy approaches the Thomas-Fermi energy, and ground states converge via Husimi functions.

Reference graph

Works this paper leans on

22 extracted references · 21 canonical work pages · cited by 2 Pith papers

[8]

C´ ardenas, and L

E. C´ ardenas, and L. Lafleche. Commutator estimates and quantitative local Weyl’s Law for Sch¨ odinger operators with non-smooth potentials. Preprint arXiv:2501.01381 (2025)

arXiv 2025
[21]

Mikkelsen

S. Mikkelsen. Sharp semiclassical spectral asymptotics for Schr¨ odinger operators with non-smooth potentials. J. Spectr. Theory 15 (2025), no. 2, pp. 819–846

work page 2025
[1]

Benedikter

N. Benedikter. Effective Dynamics of Interacting Fermions from Semiclassical Theory to the Random Phase Approximation. Journal of Mathematical Physics, 63(8):081101, Aug. 2022

work page 2022
[2]

Benedikter, M

N. Benedikter, M. Porta, and B. Schlein. Mean–Field Evolution of Fermionic Systems . Commun. Math. Phys. 331, 1087–1131 (2014)

work page 2014
[3]

Benedikter, V

N. Benedikter, V. Jakˇ si´ c, M. Porta, C. Saffirio, B, Schlein. Mean-Field Evolution of Fermionic Mixed States. Comm. Pure Appl. Math., 69: 2250-2303 (2016)

work page 2016
[4]

Benedikter, M

N. Benedikter, M. Porta, C. Saffirio, and B. Schlein. From the Hartree Dynamics to the Vlasov Equation. Arch Rational Mech Anal 221, 273–334 (2016)

work page 2016
[5]

Benedikter, P.T

N. Benedikter, P.T. Nam, M. Porta, B. Schlein, R. Seiringer. Correlation en- ergy of a weakly interacting Fermi gas. Invent. math. 225, 885–979 (2021). https://doi.org/10.1007/s00222-021-01041-5

work page doi:10.1007/s00222-021-01041-5 2021
[6]

Derivation of Hartree-Fock Dynamics and Semiclassical Commutator Estimates for Fermions in a Magnetic Field

N. Benedikter, C. Boccato, D. Monaco, N. N. Nguyen. Derivation of Hartree-Fock Dynamics and Semiclassical Commutator Estimates for Fermions in a Magnetic Field . Preprint arXiv:2503.16001 (2025)

work page Pith review arXiv 2025

Show all 22 references

[7]

C´ ardenas.Norm convergence of confined fermionic systems at zero temperature

E. C´ ardenas.Norm convergence of confined fermionic systems at zero temperature . Lett Math Phys (114), 38 (2024)

2024
[9]

C´ ardenas, J

E. C´ ardenas, J. K. Miller, D. Mitrouskas, N. Pavlovi´ c.On the stability-instability tran- sition in large Bose-Fermi mixtures Preprint arXiv:2502.18678 (2025)

2025
[10]

M. R. Christiansen, C. Hainzl, P.T. Nam. The Random Phase Approximation for Interacting Fermi Gases in the Mean-Field Regime . Forum of Mathematics, Pi. 2023;11:e32. doi:10.1017/fmp.2023.31

2023 doi
[11]

Fournais, M

S. Fournais, M. Lewin, and J. P. Solovej. The semi-classical limit of large fermionic systems. Calc. Var. (2018) 57:105

2018
[12]

Fournais, P

S. Fournais, P. S. Madsen. Semi-classical Limit of Confined Fermionic Systems in Homogeneous Magnetic Fields. Ann. Henri Poincar´ e 21, 1401–1449 (2020). 30 ESTEBAN C ´ARDENAS

2020
[13]

Girardot, N

T. Girardot, N. Rougerie. Semiclassical Limit for Almost Fermionic Anyons. Com- mun. Math. Phys. 387, 427–480 (2021)

2021
[14]

Gottschling, P.T

N. Gottschling, P.T. Nam. Convergence of Levy–Lieb to Thomas–Fermi density func- tional. Calc. Var. 57, 146 (2018)

2018
[15]

V. Ja. Ivrii, I. M. Sigal. Asymptotics of the Ground State Energies of Large Coulomb Systems. Annals of Mathematics 138, no. 2 (1993): 243–335

1993
[16]

Lafleche

L. Lafleche. Optimal semiclassical regularity of projection operators and strong Weyl law. J. Math. Phys. 1 May 2024; 65 (5): 052104

2024
[17]

Lewin, P

M. Lewin, P. S. Madsen, A. Triay. Semi-classical limit of large fermionic systems at positive temperature. J. Math. Phys. 1 September 2019; 60 (9): 091901

2019
[18]

E. H. Lieb, B. Simon. Thomas-Fermi theory revisited. Physical Review Letters, 31 (1973), p. 681

1973
[19]

E. H. Lieb, B. Simon. The Thomas-Fermi theory of atoms, molecules and solids. Advances in mathematics, 23 (1977), pp. 22–116

1977
[20]

Fournais, S

S. Fournais, S. Mikkelsen. An optimal semiclassical bound on commutators of spectral projections with position and momentum operators. Lett Math Phys 110, 3343–3373 (2020)

2020
[22]

N. N. Nguyen. Weyl laws for interacting particles . J. Math. Phys. 1 November 2024; 65 (11): 113503 (Esteban C´ ardenas)Department of Mathematics, University of Texas at Austin, 2515 Speedway, Austin TX, 78712, USA Email address : eacardenas@utexas.edu

2024

Pith tools

Reviewed August 7, 2026 · model on record in the stance chip above.

[1] [8]

C´ ardenas, and L

E. C´ ardenas, and L. Lafleche. Commutator estimates and quantitative local Weyl’s Law for Sch¨ odinger operators with non-smooth potentials. Preprint arXiv:2501.01381 (2025)

arXiv 2025

[2] [21]

Mikkelsen

S. Mikkelsen. Sharp semiclassical spectral asymptotics for Schr¨ odinger operators with non-smooth potentials. J. Spectr. Theory 15 (2025), no. 2, pp. 819–846

work page 2025

[3] [1]

Benedikter

N. Benedikter. Effective Dynamics of Interacting Fermions from Semiclassical Theory to the Random Phase Approximation. Journal of Mathematical Physics, 63(8):081101, Aug. 2022

work page 2022

[4] [2]

Benedikter, M

N. Benedikter, M. Porta, and B. Schlein. Mean–Field Evolution of Fermionic Systems . Commun. Math. Phys. 331, 1087–1131 (2014)

work page 2014

[5] [3]

Benedikter, V

N. Benedikter, V. Jakˇ si´ c, M. Porta, C. Saffirio, B, Schlein. Mean-Field Evolution of Fermionic Mixed States. Comm. Pure Appl. Math., 69: 2250-2303 (2016)

work page 2016

[6] [4]

Benedikter, M

N. Benedikter, M. Porta, C. Saffirio, and B. Schlein. From the Hartree Dynamics to the Vlasov Equation. Arch Rational Mech Anal 221, 273–334 (2016)

work page 2016

[7] [5]

Benedikter, P.T

N. Benedikter, P.T. Nam, M. Porta, B. Schlein, R. Seiringer. Correlation en- ergy of a weakly interacting Fermi gas. Invent. math. 225, 885–979 (2021). https://doi.org/10.1007/s00222-021-01041-5

work page doi:10.1007/s00222-021-01041-5 2021

[8] [6]

Derivation of Hartree-Fock Dynamics and Semiclassical Commutator Estimates for Fermions in a Magnetic Field

N. Benedikter, C. Boccato, D. Monaco, N. N. Nguyen. Derivation of Hartree-Fock Dynamics and Semiclassical Commutator Estimates for Fermions in a Magnetic Field . Preprint arXiv:2503.16001 (2025)

work page Pith review arXiv 2025