REVIEW 2 major objections 4 minor 35 references
Eigenvalue estimates for the Coulombic one-particle density matrix and the kinetic energy density matrix
T0 review · 2 major / 4 minor · reviewed 2026-08-06 · deepseek-v4-flash
Pith's one-line read For an atomic bound state, the eigenvalues of the one-particle density matrix decay at least as $k^{-8/3}$, and for antisymmetric states at least as $k^{-10/3}$.
desk verdict Solid eigenvalue bounds for Coulombic density matrices; fix the missing complex conjugate and clarify the overlap with Hearnshaw. read the letter →
The pith
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
The reading
What carries the argument
The carrying mechanism is a reduction of the eigenvalue problem to singular-value estimates for integral operators. The factorizations $\mathrm{Int}(\gamma) = \Psi^*\Psi$ and $\mathrm{Int}(\tau) = V^*V$ tie the eigenvalues to squares of singular values of $\Psi$ and $V$, whose kernels are $\psi$ and $\nabla\psi$. The main technical tool is a new singular-value bound (Theorem 4.1) for integral operators whose kernels satisfy derivative bounds with explicit factors depending on the distance to the coalescence set; the proof uses Fourier decay of such kernels (Lemma 4.3) and a sum over lattice cubes of weighted norms. Regularity input comes from derivative estimates for Coulombic eigenfunctions, with an improved version under the vanishing condition (1.10) that yields the faster decay.
What would settle it
Multiply a real bound state $\phi$ by the global phase $e^{i\pi/4}$; with the paper's definition (1.3), $\mathrm{Int}(\gamma)$ has kernel $i\,\phi(\hat{x},x)\phi(\hat{x},y)$, which is not self-adjoint, so the theorem's hypothesis that the operator is self-adjoint and non-negative fails.
Extended reading notes
Core claim
On its own terms, the paper establishes that the singular values of the integral operators with kernels $\psi(\hat{x},x)$ and $\nabla_x \psi(\hat{x},x)$ obey power-law decay with exponents determined by the regularity of $\psi$ at the coalescence points of the Coulomb potential. Concretely, $\lambda_k(\mathrm{Int}(\gamma)) \le C k^{-8/3} \rho_{3/8}$ and $\lambda_k(\mathrm{Int}(\tau)) \le C k^{-2} \rho_{1/2}$ always hold, and the vanishing condition (1.10) improves these to $\lambda_k(\mathrm{Int}(\gamma)) \le C k^{-10/3} \rho_{3/10}$ and $\lambda_k(\mathrm{Int}(\tau)) \le C k^{-8/3} \rho_{3/8}$. The constants do not depend on $\psi$, but may depend on $E$ and $N$; the norms $\rho_q$ are sums over unit cubes of powers of the one-particle density. The proof requires only polynomial decay of $\psi$ at infinity, not the exponential bound assumed in earlier work.
Load-bearing premise
The proofs assume the eigenfunction is real-valued; for complex eigenfunctions the density kernel (1.3) is not Hermitian, so the stated eigenvalues $\lambda_k(\mathrm{Int}(\gamma))$ are not defined.
Editorial extensions
If this is right
- Occupation numbers of an atomic bound state are guaranteed to decay at least as $k^{-8/3}$, and at least as $k^{-10/3}$ for antisymmetric states.
- The eigenvalues of the kinetic energy density matrix decay at least as $k^{-2}$, improving to $k^{-8/3}$ under the vanishing condition.
- Because the constants do not depend on the eigenfunction, the bounds apply uniformly across all bound states of a given atom with fixed energy and electron number.
- The statements extend to molecules with fixed nuclei, as noted in the paper.
- These bounds can be used to control errors in finite-basis approximations used in quantum chemistry.
Reading between the lines
- The faster decay for antisymmetric states suggests that Pauli exclusion itself drives the improved summability of occupation numbers, which could be tested numerically by comparing singlet and triplet pair configurations in few-electron atoms.
- The same Fourier-decay machinery likely yields analogous bounds for $n$-particle reduced density matrices with $n>1$, where the coalescence sets have more complex geometry.
- If the implicit real-valued assumption were relaxed, the density matrix would need a complex conjugate in its kernel; the stated results apply to real eigenfunctions, which always exist for this real Hamiltonian.
- One could probe sharpness by constructing model eigenfunctions with cusp-like singularities producing Fourier decay exactly at the boundary exponents.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. This paper studies the one-particle density matrix γ and the kinetic energy density matrix τ associated with an eigenfunction ψ of the N-electron atomic Schrödinger operator H = -Δ + V. The main results are the eigenvalue bounds λ_k(Int(γ)) ≤ C k^{-8/3} ρ_{3/8} and λ_k(Int(τ)) ≤ C k^{-2} ρ_{1/2} (Theorem 1.1), and improved bounds λ_k(Int(γ)) ≤ C k^{-10/3} ρ_{3/10} and λ_k(Int(τ)) ≤ C k^{-8/3} ρ_{3/8} when ψ vanishes at electron-electron coalescence points (Theorem 1.2). The constants depend only on E and N, and the right-hand sides contain explicit norms of the one-particle density ρ. The proof factorizes Int(γ) and Int(τ) as Ψ*Ψ and V*V, estimates the singular values of Ψ and V by bounding Fourier coefficients of the kernels via derivative estimates (2.3) and (2.5), and then uses Schatten/quasi-norm machinery. The exposition is clear and the proof strategy is direct and essentially self-contained, conditional on the regularity results from Fournais–Sørensen.
Significance. The proof strategy is transparent and, conditional on the cited regularity estimates, avoids the general Birman–Solomyak apparatus by directly estimating Fourier coefficients of the kernels. The explicit dependence of the bounds on the eigenfunction through lattice norms of ρ is a genuine improvement over previous results with implicit constants, and the improved decay rates under the antisymmetry/vanishing condition appear to be new. No parameters are fitted, and the earlier asymptotic results are used only to indicate sharpness, not inside the derivation. The main obstacle to accepting the paper in its present form is the incomplete definition of the density matrices for complex-valued eigenfunctions, which makes the central statements ill-posed as written; this is a fixable but load-bearing gap.
major comments (2)
- [Section 1, Eqs. (1.3)–(1.4) and Theorems 1.1–1.2] The kernels γ0 and τ0 are defined without complex conjugation. For a complex-valued eigenfunction ψ, the operator Int(γ0) with kernel γ0(x,y)=∫ψ(ˆx,x)ψ(ˆx,y)dˆx is not self-adjoint in general, so the eigenvalues λ_k(Int(γ0)) are not defined. The factorization (1.13) gives Int(γ)=Ψ*Ψ, whose kernel is ∫ \overline{ψ(ˆx,y)}ψ(ˆx,x)dˆx, not the kernel in (1.3) unless ψ is real. Similarly, ρ(x)=γ(x,x) in (1.7) equals ∫|ψ|^2 only after conjugation. A concrete failure occurs for N=1 and a hydrogenic eigenfunction with m≠0: Int(γ) is rank-one with sole eigenvalue ∫ψ^2, generally non-real. Since H is real, the results can be salvaged either by assuming ψ is real, or, preferably, by inserting the complex conjugate in (1.3)–(1.4) (and in the definition of τ). The proofs in Section 5 go through unchanged because the singular-value estimates apply identically to kernels with \overline{ψ}.
- [Section 5.1, proof of Theorem 5.1] The kernel b(ˆx)ψ(ˆx,x) is said to satisfy condition (4.1) with α=1 and z_j(x)=x_j, j=1,...,N-1. The regularity bound (2.3) used here involves λ(x)=min{|x|, 1/√2 |x-x_k|}, which includes the nuclear singularity at x=0. Thus the point x=0 must be included among the z_k in (4.1); as written, the stated list of N-1 points does not cover it. This is evidently a typographical omission—Lemma 4.3 and Theorem 4.1 allow N arbitrary singular points—but it should be corrected in the proof.
minor comments (4)
- [Theorems 1.1 and 1.2] The finiteness of the norms ρ_{3/8}, ρ_{1/2}, and ρ_{3/10} is not stated as a hypothesis. The bounds are only meaningful when these norms are finite, as noted in the text after Theorem 1.1; it would be clearer to include this assumption in the theorem statements.
- [Section 1, after (1.3)] The phrase 'a.e. x ∈ R3, a.e. y ∈ R3' would be cleaner as 'a.e. x, y ∈ R3'.
- [Proof of Lemma 4.3] The summation over 'p,s:|p|+|s|=κ' should state explicitly that p and s are multi-indices and that the sum runs over all such multi-indices.
- [Eq. (2.23)] The notation '0≤k<m' in the summation is non-standard; it should be clarified that the sum runs over multi-indices k with k_j ≤ m_j for all j and k ≠ m.
Circularity Check
No circularity: the eigenvalue bounds are derived from external regularity estimates and elementary s-value bounds, with no fitted parameters or self-citation chain carrying the claim.
full rationale
The derivation is self-contained in the relevant sense. The main theorems are proved in Section 5 by combining the exact factorization Int(γ)=Ψ*Ψ and Int(τ)=V*V (1.13) with singular-value estimates for integral operators (Theorem 4.1, proved via the Fourier-transform bounds of Lemmas 4.3 and 4.4) and pointwise derivative bounds for Coulombic eigenfunctions (Proposition 2.1 and Theorem 2.2). The derivative bounds are imported from external regularity theory, chiefly Fournais–Sørensen [15] and Hoffmann-Ostenhof et al. [13]; the cited work by Hearnshaw–Sobolev [21] is a self-citation, but it supplies a general elliptic regularity proposition (Proposition 2.6) rather than the target eigenvalue decay, and the paper actually proves Theorem 2.5 in detail rather than invoking it as a black box. The earlier asymptotic results of the author [32,33] are referenced only to assert sharpness of the exponents (1.6); they are not used inside the proof of Theorems 1.1 and 1.2. No parameter is fitted, no nontrivial quantity is renamed, and no uniqueness theorem is invoked to force a choice. The right-hand sides ρ_{3/8}, ρ_{1/2}, ρ_{3/10}, and ρ_{3/8} arise from integrating |ψ|^2 over lattice cubes during the s-value estimates, not from the eigenvalues being bounded. There is, however, a well-posedness/notation gap that is not circular: formulas (1.3)–(1.4) omit the complex conjugate, so Int(γ) and Int(τ) are self-adjoint as asserted only for real-valued ψ (or after inserting conjugates in the kernels); this affects correctness for complex eigenfunctions, not the circularity status of the argument.
Assumptions & free parameters
assumptions (5)
- domain assumption Eigenfunction ψ ∈ D(H) = W^{2,2}(R^{3N}) with Hψ = Eψ, and the one-particle density ρ has finite lattice norm ρ_q.
- standard math Derivative bound |∂_x^m ψ(x)| ≲ (1 + λ(x)^{1-|m|}) ||ψ||_{L^2(B(x,R))} (Proposition 2.1), imported from [15, Corollary 1.2].
- standard math Enhanced regularity under (1.10): e^{-F0}ψ has bounded second x-derivatives, giving |∂_x^m(e^{-F0}ψ)| ≲ (1+λ^{2-|m|}) ||ψ||_{L^2} (Theorem 2.2), built from [13] and [20].
- standard math Elliptic regularity bootstrap, Proposition 2.6, quoted from [21, Theorem 3.2] (the author's own paper with Hearnshaw).
- standard math Schatten-von Neumann s-number inequalities, Propositions 3.1 and 3.2, from [3], [1], [2].
Cite this review
Pith. "Pith review of Eigenvalue estimates for the Coulombic one-particle density matrix and the kinetic energy density matrix." pith.science (2026). https://pith.science/paper/FO3S36AR
@misc{pith2026250616178,
author = {Pith},
title = {Pith review of: Eigenvalue estimates for the Coulombic one-particle density matrix and the kinetic energy density matrix},
year = {2026},
howpublished = {\url{https://pith.science/paper/FO3S36AR}},
note = {Machine review of arXiv:2506.16178}
}
abstract
Consider a bound state (an eigenfunction) $\psi$ of an atom with $N$ electrons. We study the spectra of the one-particle density matrix $\gamma$ and of the one-particle kinetic energy density matrix $\tau$ associated with $\psi$. The paper contains two results. First, we obtain the bounds $\lambda_k(\gamma)\le C_1 k^{-8/3}$ and $\lambda_k(\tau)\le C_2 k^{-2}$ with some positive constants $C_1, C_2$ that depend explicitly on the eigenfunction $\psi$. The sharpness of these bounds is confirmed by the asymptotic results obtained by the author in earlier papers. The advantage of these bounds over the ones derived by the author previously, is their explicit dependence on the eigenfunction. Moreover, their new proofs are more elementary and direct. The second result is new and it pertains to the case where the eigenfunction $\psi$ vanishes at the particle coalescence points. In particular, this is true for totally antisymmetric $\psi$. In this case the eigenfunction $\psi$ exhibits enhanced regularity at the coalescence points which leads to the faster decay of the eigenvalues: $\lambda_k(\gamma)\le C_3 k^{-10/3}$ and $\lambda_k(\tau)\le C_4 k^{-8/3}$. The proofs rely on the estimates for the derivatives of the eigenfunction $\psi$ that depend explicitly on the distance to the coalescence points. Some of these estimates are borrowed directly from, and some are derived using the methods of a recent paper by S. Fournais and T. \O. S\o rensen.
Reference graph
Works this paper leans on
-
[20]
Eigenvalue Bounds for Multi-Particle Reduced Density Matrices of Coulombic Wavefunctions
P. Hearnshaw, Eigenvalue Bounds for Multi-Particle Reduced Density Matrices of Coulombic Wave- functions. arXiv:2412.16073 [math-ph] 2024
work page Pith review arXiv 2024
-
[1]
A.B. Aleksandrov, S. Janson, V.V. Peller and R. Rochberg, An interesting class of operators with unusual Schatten-von Neumann behavior . Function spaces, interpolation theory and related topics (Lund, 2000) , 61–149, de Gruyter, Berlin, 2002
work page 2000
-
[2]
M. S. Birman and M. Z. Solomyak, Estimates for the singular numbers of integral operators. (Rus- sian). Uspehi Mat. Nauk 32(1(193)): 17–84, 1977
work page 1977
-
[3]
M. S. Birman and M. Z. Solomyak, Spectral Theory of Selfadjoint Operators in Hilbert Space . Mathematics and its Applications (Soviet Series), D. Reidel, 1987. Translated from the 1980 Russian original by S. Khrushch¨ ev and V. Peller
work page 1987
-
[4]
J. Cioslowski, Off-diagonal derivative discontinuities in the reduced density matrices of electronic systems. The Journal of Chemical Physics 153(15): 154108, 2020
work page 2020
-
[5]
J. Cioslowski and F. Pratnicki, Universalities among natural orbitals and occupation numbers per- taining to ground states of two electrons in central potentials . The Journal of Chemical Physics 151(184107), 2019
work page 2019
-
[6]
J. Cioslowski and K. Strasburger, Angular-Momentum Extrapolations to the Complete Basis Set Limit: Why and When They Work . Journal of Chemical Theory and Computation 17(6): 3403– 3413, 2021
work page 2021
-
[7]
A. Coleman and V. Yukalov, Reduced Density Matrices , Lecture Notes in Chemistry , vol. 72. Springer-Verlag Berlin Heidelberg, 2000
work page 2000
Show all 35 references
-
[8]
J. M. Combes and L. Thomas, Asymptotic behaviour of eigenfunctions for multiparticle Schr¨ odinger operators. Comm. Math. Phys. 34: 251–270, 1973
1973
-
[9]
Davidson, Reduced Density Matrices in Quantum Chemistry
E. Davidson, Reduced Density Matrices in Quantum Chemistry . Academic Press, 1976
1976
-
[10]
Deift, W
P. Deift, W. Hunziker, B. Simon and E. Vock, Pointwise bounds on eigenfunctions and wave packets in N -body quantum systems. IV . Comm. Math. Phys. 64(1): 1–34, 1978/79
1978
-
[11]
Fournais, M
S. Fournais, M. Hoffmann-Ostenhof, T. Hoffmann-Ostenhof, T.Ø. Sørensen, The electron density is smooth away from the nuclei . Comm. Math. Phys. 228(3): 401–415, 2002
2002
-
[12]
Fournais, M
S. Fournais, M. Hoffmann-Ostenhof, T. Hoffmann-Ostenhof, T.Ø. Sørensen, Analyticity of the den- sity of electronic wavefunctions . Ark. Mat. 42(1): 87–106, 2004
2004
-
[13]
Fournais, M
S. Fournais, M. Hoffmann-Ostenhof, T. Hoffmann-Ostenhof, T.Ø. Sørensen, Sharp regularity results for Coulombic many-electron wave functions . Comm. Math. Phys. 255(1): 183–227, 2005. 22 ALEXANDER V. SOBOLEV
2005
-
[14]
Fournais, M
S. Fournais, M. Hoffmann-Ostenhof, T. Hoffmann-Ostenhof, T.Ø. Sørensen, Analytic structure of many-body Coulombic wave functions . Comm. Math. Phys. 289(1): 291–310, 2009
2009
-
[15]
Fournais and T
S. Fournais and T. Ø. Sørensen, Estimates on derivatives of Coulombic wave functions and their electron densities. Journal f¨ ur die reine und angewandte Mathematik (Crelles Journal)2021(775): 1–38, 2021
2021
-
[16]
Friesecke, The multiconfiguration equations for atoms and molecules: charge quantization and existence of solutions
G. Friesecke, The multiconfiguration equations for atoms and molecules: charge quantization and existence of solutions . Arch. Ration. Mech. Anal. 169(1): 35–71, 2003
2003
-
[17]
Friesecke, On the infinitude of non-zero eigenvalues of the single-electron density matrix for atoms and molecules
G. Friesecke, On the infinitude of non-zero eigenvalues of the single-electron density matrix for atoms and molecules . R. Soc. Lond. Proc. Ser. A Math. Phys. Eng. Sci. 459(2029): 47–52, 2003
2003
-
[18]
Froese and I
R. Froese and I. Herbst, Exponential bounds and absence of positive eigenvalues for N -body Schr¨ odinger operators. Comm. Math. Phys. 87(3): 429–447, 1982/83
1982
-
[19]
H¨ attig, W
Ch. H¨ attig, W. Klopper, A. K¨ ohn and D. P. Tew, Explicitly Correlated Electrons in Molecules . Chem. Rev 112(1): 4–74, 2012
2012
-
[21]
Hearnshaw and A
P. Hearnshaw and A. V. Sobolev, The diagonal behaviour of the one-particle Coulombic density matrix. Probab. Math. Phys. 4(4): 935–964, 2023
2023
-
[22]
Hoffmann-Ostenhof, T
M. Hoffmann-Ostenhof, T. Hoffmann-Ostenhof, and T. Ø. Sørensen, Electron wavefunctions and densities for atoms . Ann. Henri Poincar´ e2(1): 77–100, 2001
2001
-
[23]
Hoffmann-Ostenhof, T
M. Hoffmann-Ostenhof, T. Hoffmann-Ostenhof, and H. Stremnitzer, Local properties of Coulombic wave functions . Comm. Math. Phys. 163(1): 185–215, 1994
1994
-
[24]
Kato, On the eigenfunctions of many-particle systems in quantum mechanics
T. Kato, On the eigenfunctions of many-particle systems in quantum mechanics . Comm. Pure Appl. Math. 10: 151–177, 1957
1957
-
[25]
Leray, Sur les solutions de l’´ equation de Schr¨ odinger atomique et le cas particulier de deux ´ electrons
J. Leray, Sur les solutions de l’´ equation de Schr¨ odinger atomique et le cas particulier de deux ´ electrons. Trends and applications of pure mathematics to mechanics (Palaiseau, 1983) , Lecture Notes in Phys. , vol. 195, 235–247, Springer, Berlin, 1984
1983
-
[26]
Lewin, Solutions of the multiconfiguration equations in quantum chemistry
M. Lewin, Solutions of the multiconfiguration equations in quantum chemistry . Arch. Ration. Mech. Anal. 171(1): 83–114, 2004
2004
-
[27]
Lewin, Geometric methods for nonlinear many-body quantum systems
M. Lewin, Geometric methods for nonlinear many-body quantum systems. J. Funct. Anal. 260(12): 3535–3595, 2011
2011
-
[28]
Lewin, E
M. Lewin, E. H. Lieb, and R. Seiringer, Universal Functionals in Density Functional Theory . 2019. 1912.10424
2019 arXiv
-
[29]
E. H. Lieb and R. Seiringer, The stability of matter in quantum mechanics . Cambridge University Press, Cambridge, 2010
2010
-
[30]
Reed and B
M. Reed and B. Simon, Methods of modern mathematical physics. II. Fourier analysis, self- adjointness. Academic Press [Harcourt Brace Jovanovich, Publishers], New York-London, 1975
1975
-
[31]
Simon, Exponential decay of quantum wave functions,
B. Simon, Exponential decay of quantum wave functions, . http://www.math.caltech.edu/simon/Selecta/ExponentialDecay.pdf, Online notes, part of B. Simon’s Online Selecta at http://www.math.caltech.edu/simon/selecta.html
-
[32]
A. V. Sobolev, Eigenvalue asymptotics for the one-particle density matrix . Duke Math. J. 171(17): 3481–3513, 2022
2022
-
[33]
A. V. Sobolev, Eigenvalue asymptotics for the one-particle kinetic energy density operator. J. Funct. Anal. 283(8): Paper No. 109604, 35, 2022
2022
-
[34]
A. V. Sobolev, Eigenvalue estimates for the one-particle density matrix . J. Spectr. Theory 12(2): 857–875, 2022
2022
-
[35]
Szabo and N
A. Szabo and N. S. Ostlund, Modern Quantum Chemistry: Introduction to Advanced Electronic Structure Theory. Dover Publications, 1996. DENSITY MATRIX 23 Department of Mathematics, University College London, Gower Street, London, WC1E 6BT UK Email address : a.sobolev@ucl.ac.uk
1996
Reviewed August 6, 2026 · model on record in the stance chip above.
Discussion (0). Continue with ORCID to comment.