Eigenvalue Bounds for Multi-Particle Reduced Density Matrices of Coulombic Wavefunctions

Peter Hearnshaw

REVIEW 1 major objections 3 minor 1 cited by

Eigenvalue Bounds for Multi-Particle Reduced Density Matrices of Coulombic Wavefunctions

T0 review · 1 major / 3 minor · reviewed 2026-08-11 · deepseek-v4-flash

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

Pith's one-line read This paper proves that every K-particle reduced density matrix of a Coulombic bound state, for 2 ≤ K ≤ N−2, has eigenvalues λ_n ≤ C n^{-α_K} with α_K = 1 + 7/(3L), where L = min{K, N−K}.

desk verdict Solid new result with a small repairable gap in Theorem 1.2: the symmetry step for K>N/2 is stated but never executed. read the letter →

arxiv 2412.16073 v1 pith:BPRU2MIN submitted 2024-12-20 math-ph math.MP

Peter Hearnshaw This is my paper · ORCID

classification math-phmath.MP MSC 35J1046E3581V5581V70

keywords reduceddensitymatrixCoulombicwavefunctioneigenvaluedecaySchattenclassesBesov-Nikol'skiispacesJastrowfactorBirman-Solomyakboundsmany-electronSchrödingerequation

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

This paper proves that the K-particle reduced density matrices of Coulombic bound states, for any number K between 2 and N−2 of an N-electron atom or molecule, have eigenvalue sequences that decay faster than the trace-class rate. Concretely, if ψ is a normalized eigenfunction of the N-electron Schrödinger equation with Coulomb potential and exponential decay (as holds for discrete eigenvalues), then λ_n(Γ^(K)) ≤ C $n^{{-α_K}}$ with α_K = 1 + 7/(3L), where L = min{K, N−K}. The same bound holds for Γ^(N−K), for permuted density matrices, and therefore for the averaged K-particle operator. The interest is that it fills the previously missing range 2 ≤ K ≤ N−2, where no eigenvalue decay beyond trace class was known, and it gives quantitative rates for finite-rank approximation of the operators used in quantum chemistry.

What carries the argument

The engine is the factorization ψ(ẋ,x̂) = A(ẋ)B(x̂)C(ẋ,x̂)μ(ẋ,x̂) built from Jastrow factors exp(-(Z/2)τ) and exp(τ/4), with τ(x) = |x| - (1 + |x|^2)^{1/2}, and an exponential weight that isolates the decay. The remainder μ satisfies pointwise derivative bounds |∂^m_{x_j} μ(x)| ≤ C(1 + |x_j|^{min{2-|m|,0}} + Σ_{k≠j} |x_j - x_k|^{min{2-|m|,0}}). These bounds put μ, after localization to unit cubes, in the Besov–Nikol'skii space $N^{{7/2}}$_2 in each block of K coordinates, via a product-space lemma that builds full $R^{{3K}}$ regularity from coordinate-wise regularity. Birman–Solomyak bounds convert that Besov regularity into singular-value decay for the integral operator with kernel μ, the remaining factors A and B act as integrable weights, and C acts as a multiplier on the Schatten spaces, so it does not worsen the decay.

What would settle it

One concrete check is to compute, on unit cubes Q_n in $R^{{3K}}$, the Besov–Nikol'skii norm of μ(·,x̂); the proof requires this norm to be bounded independently of n and x̂, so any sequence of cubes where it grows without bound would falsify the central estimate. A more direct check is to approach a two-particle collision and test whether |∂^3_{x_j} φ(x)| exceeds the claimed bound C $e^{{-κ|x|'}}$(1 + |x_j|^{-1} + Σ_{k≠j}|x_j - x_k|^{-1}).

Watch

Extended reading notes

Core claim

On the paper's own terms, the central discovery is that every K-particle reduced density matrix of a Coulombic bound state lies in the weak Schatten class S_{1/α_K,∞}, equivalently its n-th eigenvalue decays as $n^{{-α_K}}$. The proof obtains this by writing the K-particle operator as (Ψ^(K))^*Ψ^(K) with Ψ^(K) a Hilbert-Schmidt operator whose kernel is ψ restricted to K of the N particle coordinates, then factorizing that kernel as A(ẋ)B(x̂)C(ẋ,x̂)μ(ẋ,x̂). The Jastrow-type factors absorb the Coulomb singularities, leaving a remainder μ with Besov regularity of order 7/2 in each coordinate; Birman–Solomyak singular-value estimates then place Ψ^(K) in S_{q,∞} with 1/q = 1/2 + 7/(6K), and the identity λ_n(Γ^(K)) = s_n(Ψ^(K))^2 converts this to the stated decay.

Load-bearing premise

The load-bearing premise is the pointwise derivative bound for the Jastrow-smoothed wavefunction φ = $e^{{-F}}$ψ stated in Proposition 2.1, which is assembled from two existing results rather than proved in full; if that bound failed, or if the constants grew differently than stated, the Besov regularity input and the final eigenvalue decay would not follow.

Editorial extensions

If this is right

For every K between 1 and N−1, including the previously open range 2 ≤ K ≤ N−2, the K-particle reduced density matrix now has eigenvalue decay at least n^{-α_K}, which for K close to N/2 gives exponent 1 + 14/(3N).
Because α_K > 1 for every allowed K, the bound is a genuine improvement over the automatic trace-class decay λ_n = O(n^{-1}).
The same rate holds for all permuted operators Γ^(K)_σ and for the averaged Γ^(K), since weak Schatten classes are vector spaces and the construction is permutation-insensitive.
The result is symmetric under K ↔ N−K, so replacing K by N−K changes nothing in the exponent.
For K = 2 the decay estimate gives quantitative control on how well the two-particle reduced density matrix, from which the ground-state energy can be recovered, is approximated by finite-rank operators.

Reading between the lines

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

If the author's cusp heuristic is right, the optimal decay for every K may in fact be the same 8/3 seen for K=1, so the exponents proved here would be conservative; testing this would require a sharper analysis of fifth-order cusps in γ^(K) along partial diagonals.
The same Jastrow-factor strategy could plausibly be adapted to the kinetic-energy density matrix or to temperature-dependent Gibbs states, where analogous factorization and Besov regularity arguments might give decay rates for higher-order correlation operators.
The product-space Besov lemma appears to be the piece that unlocks intermediate K; a natural check is whether a larger class of singular potentials with similar cusp behaviour yields the same eigenvalue rates, which would show the phenomenon is not specific to Coulomb interactions.

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

Dear [Colleague],

Quick take: this is a real result. Hearnshaw extends Sobolev's eigenvalue decay bounds to K-particle reduced density matrices for all 2≤K≤N-2, with rate n^{-(1+7/(3L))}, L=min{K,N-K}. That range was genuinely open, and the proof machinery—factorization ψ=ABC µ, Besov regularity of µ, Birman-Solomyak bounds, and a clever multiplier lemma for functions that ignore some variables—is new and coherent. There are no fitted parameters, no circular reasoning; the dependence on the author's earlier papers is for regularity estimates that are already published, not for the decay result itself.

The main soft spot is in the proof of Theorem 1.2 as printed. The argument applies Proposition 1.4 once, to the given K, and concludes λ_n(Γ^(K)) ≤ C n^{-(1+7/(3K))}. That's fine for K≤N/2. For K>N/2, the stated exponent uses N-K, and the written proof never re-runs Proposition 1.4 with the complementary K. The sentence about taking any permutation only swaps Γ^(N-K)_{σ_K} for Γ^(N-K); it does not improve the decay rate. The fix is one line: apply Proposition 1.4 with K replaced by N-K, then use (1.13) to transfer the bound. All the ingredients are there, so this is a minor but real gap in exposition—not a flaw in the underlying strategy.

I have no comparable concerns about the technical lemmas. Lemma 4.3 (tensor-product Besov inclusion) and Lemma 5.3 (multiplier for partial-variable functions) look correct, and the Birman-Solomyak application is standard. Proposition 2.1 is imported from [2] and [19] with a sketch in the appendix; that's a reasonable dependency, not a red flag.

Bottom line: this deserves a serious referee. The result is new, the proof is largely solid, and the gap is easily repairable. I'd ask the author to add the symmetry step and resubmit. Worth citing and worth discussing in a reading group if you care about spectral theory of reduced density matrices.

Best, [Your name]

Referee Report

1 major / 3 minor

Summary. The paper studies eigenvalue decay of K-particle reduced density matrices Γ^(K) for Coulombic N-electron bound-state wavefunctions. For 2 ≤ K ≤ N-2 it claims a bound λ_n(Γ^(K)) ≤ C n^{-α_K} with α_K = 1 + 7/(3 min{K,N-K}), an improvement over the trivial trace-class decay. The strategy is to factorize the wavefunction as ψ = A B C µ via Jastrow factors, use pointwise derivative bounds for µ to place its slices in Besov spaces, apply Birman-Solomyak singular-value bounds, and treat the factor C as a Schatten-class multiplier. The proof is largely self-contained except for reliance on earlier regularity results.

Significance. If the main theorem holds, it is a solid extension of Sobolev's results for K=1,N-1 to all intermediate K, with decay strictly better than trace-class. The technical apparatus (Besov spaces, Schatten-class multipliers, Birman-Solomyak bounds) is well matched to the problem. The paper contains no adjustable parameters and makes no symmetry assumptions, so the result applies to general bound states, including molecules. The author explicitly acknowledges that the exponents are likely non-optimal, which is appropriate. The main theorem's proof, however, currently does not cover K > N/2 as stated, because Proposition 1.4 is applied only with the original K. Since the repair is straightforward, the underlying approach is sound and the paper will be a useful contribution after revision.

major comments (1)

[Section 1.1, proof of Theorem 1.2] The proof as written applies Proposition 1.4 only to the given K and concludes λ_n(Γ^(K)) ≤ C n^{-(1+7/(3K))}. For K > N/2 this is strictly weaker than the claimed exponent α_K = 1+7/(3(N-K)). The sentence about permuting variables does not change the exponent in the denominator. To prove the stated bound for K > N/2, one must also apply Proposition 1.4 with the argument N-K (which is also in the range 2 ≤ · ≤ N-2) and then transfer the resulting S_{q',∞} estimate to Γ^(K) via (1.13) and Remark 1.3(1). All ingredients for this repair are present, but the theorem as printed is not established for half its stated range.

minor comments (3)

[Appendix A, proof of Proposition 2.1] The verification of condition (A.6) is delegated to [6, Corollary 4.5] without details. Since Proposition 2.1 provides the derivative bounds that drive the Besov regularity of µ, please spell out this verification or give a precise statement and proof for the exact form (2.4). The current presentation is plausible but not self-contained.
[Section 5, Lemma 5.3] The proof of Lemma 5.3 explicitly treats only q ∈ (q0,1). When q0 ≥ 1, the stated range q ∈ (q0,2) is not covered by the argument as written. The application in Lemma 5.4 has q0 = 3/4, so this does not affect the main result, but the lemma is stated in greater generality than proved.
[Section 1.1, equation (1.13)] The complex conjugation of the operator is mentioned only parenthetically in the discussion following (1.13). A precise definition of the conjugated operator in terms of the kernel of Γ would improve clarity.

Circularity Check

0 steps flagged · score 2.0 of 10

No circularity found: the derivation is self-contained; the printed proof of Theorem 1.2 has a correctness gap for K>N/2, not a circular step.

full rationale

The paper contains no fitted parameters and no claim that reduces to its own inputs. The main theorem, Theorem 1.2, is obtained from Proposition 1.4, whose proof applies Proposition 6.1 with α=2. The α=2 regularity input is Proposition 2.1, derived in Appendix A from published regularity theorems: [19, Theorem 1.5/Remark 1.6] (external to the author) and [6, Theorem 3.2] (self-authored). These cited results concern pointwise and elliptic regularity of Coulombic wavefunctions and Besov-space criteria; none of them assumes or contains the target eigenvalue-decay statement. Hence the self-citations are independent technical tools, not circular premises. The Besov lemmas and multiplier machinery are either original (Lemma 4.3) or taken from external sources [17], [23]. The factorization ψ = ABCμ in (2.10) is an identity by construction, not a fitted input. The numerical exponent is computed from the regularity order α=2 via the Birman–Solomyak formula, not selected to match the conclusion. I therefore find no circular step. A separate non-circular concern must be flagged: the proof of Theorem 1.2 as printed applies Proposition 1.4 to the single given K and concludes λ_n(Γ^(K)) ≤ Cn^{-(1+7/(3K))}. For K>N/2, the stated α_K = 1+7/(3(N-K)) is stronger than this. The sentence 'Since we could have taken any permutation of the variables of ψ in defining Ψ^(K)' only changes Γ^(N-K)_{σ_K} to Γ^(N-K); it does not apply Proposition 1.4 with K replaced by N-K. The missing step is repairable—apply Proposition 1.4 to N-K and use (1.13)—but as printed the theorem's full range is not established. This is a correctness gap, not circularity.

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

No free parameters are fitted. The result rests on standard Schatten-class and Besov-space theory, on known exponential decay of eigenfunctions, and on the Jastrow-regularity estimates of Fournais et al. and the author's prior work. No new entities are postulated.

assumptions (4)

domain assumption The eigenfunction ψ obeys exponential decay (1.5) for discrete eigenvalues of the Coulomb Hamiltonian.
Invoked in Section 1 to extract exponential decay into A, B and to define µ; standard result, see [10].
domain assumption Jastrow-transformed wavefunction φ = e^{-F}ψ satisfies the pointwise derivative bounds of Proposition 2.1, derived from [19] and [2].
The cornerstone regularity estimate; proof in Appendix A relies on [19, Theorem 1.5] and [6, Corollary 4.5].
standard math Birman-Solomyak singular value bounds (Proposition 3.2) and Schatten-class multiplier theorems (Propositions 5.1, 5.2) from [17].
Used to convert Besov regularity into S_{q,∞} inclusions and to handle the Jastrow cross-term C.
standard math Besov space embedding properties and Lemma 4.1 from [23].
Background for N^s_2 spaces and finite-difference seminorms.

how reviews work

0 comments

Cite this review

Pith. "Pith review of Eigenvalue Bounds for Multi-Particle Reduced Density Matrices of Coulombic Wavefunctions." pith.science (2026). https://pith.science/paper/BPRU2MIN

@misc{pith2026241216073,
  author       = {Pith},
  title        = {Pith review of: Eigenvalue Bounds for Multi-Particle Reduced Density Matrices of Coulombic Wavefunctions},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/BPRU2MIN}},
  note         = {Machine review of arXiv:2412.16073}
}

abstract

For bound states of atoms and molecules of $N$ electrons we consider the corresponding $K$-particle reduced density matrices, $\Gamma^{(K)}$, for $1 \le K \le N-1$. Previously, eigenvalue bounds were obtained in the case of $K=1$ and $K=N-1$ by A.V. Sobolev. The purpose of the current work is to obtain bounds in the case of $2 \le K \le N-2$. For such $K$ we label the eigenvalues of the positive, trace class operators $\Gamma^{(K)}$ by $\lambda_n(\Gamma^{(K)})$ for $n=1,2,\dots$, and obtain the bounds $\lambda_n(\Gamma^{(K)}) \le Cn^{-\alpha_K}$ for all $n$, where $\alpha_K = 1 + 7/(3L)$ and $L = \min\{K,N-K\}$.

Discussion (0). Continue with ORCID to comment.

Forward citations

Cited by 1 Pith paper

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

Eigenvalue estimates for the Coulombic one-particle density matrix and the kinetic energy density matrix
math.SP 2025-06 conditional novelty 6.0 of 10

For Coulombic many-electron eigenfunctions, the occupation-number eigenvalues decay at least as k^{-8/3} (density matrix) and k^{-2} (kinetic energy matrix), improving to k^{-10/3} and k^{-8/3} when the wavefunction v...

Reference graph

Works this paper leans on

26 extracted references · 26 canonical work pages · cited by 1 Pith paper

[2]

Hearnshaw and A.V

P. Hearnshaw and A.V. Sobolev. Analyticity of the One-Particle De nsity Matrix. Ann. Henri. Poincare, 23:707–738, 2022

work page 2022
[19]

Fournais, M

S. Fournais, M. Hoﬀmann-Ostenhof, T. Hoﬀmann-Ostenhof, and T.Ø. Sørensen. Sharp Regularity Results for Coulombic Many-Electron Wave Functions. Comm. Math. Phys. , 255:183–227, 2005

work page 2005
[1]

T. Kato. On the eigenfunctions of many-particle systems in quan tum mechanics. Comm. Pure Appl. Math., 10:151–177, 1957

work page 1957
[3]

T. Jecko. On the analyticity of electronic reduced densities for m olecules. J. Math. Phys. , 63(1),

work page
[4]

Cioslowski

J. Cioslowski. Oﬀ-diagonal derivative discontinuities in the reduce d density matrices of electronic systems. J. Chem. Phys , 153(154108), 2020

work page 2020
[5]

Cioslowski

J. Cioslowski. Reverse engineering in quantum chemistry: How to r eveal the ﬁfth-order oﬀ-diagonal cusp in the one-electron reduced density matrix without actually ca lculating it. Int. J. Quantum. Chem., 122(8), 2022. e26651

work page 2022
[6]

Hearnshaw and A.V

P. Hearnshaw and A.V. Sobolev. The diagonal behaviour of the on e-particle Coulombic density matrix. Probability and Mathematical Physics , 4(4):935–964, 2023

work page 2023
[7]

Boundedness of the Fifth Derivative for the One-Particle Coulombic Density Matrix at the Diagonal

P. Hearnshaw. Boundedness of the Fifth Derivative for the One -Particle Coulombic Density Matrix at the Diagonal. 2024. arXiv, 2301.02848. EIGENV ALUE BOUNDS FOR REDUCED DENSITY MATRICES 19

work page Pith review arXiv 2024

Show all 26 references

[8]

Jecko and C

T. Jecko and C. Noˆ us. Limited regularity of a speciﬁc electronic r educed density matrix for molecules. 2023. arXiv, 2309.06318

2023 arXiv
[9]

Jecko and C

T. Jecko and C. Noˆ us. Regularity of the (N-1)-particle electro nic reduced density matrix for molecules with ﬁxed nuclei and N electrons. 2024. arXiv, 2407.03706

2024 arXiv
[10]

B. Simon. Exponential Decay of Quantum Wave Functions. Online notes, part of B. Simon’s Online Selecta. Available at http://www.math.caltech.edu/simon/Selecta/E xponentialDecay.pdf
[11]

Friesecke

G. Friesecke. On the inﬁnitude of non-zero eigenvalues of the s ingle-electron density matrix for atoms and molecules. Proc. Math. Phys. Eng. , 459(2029):47–52, 2003

2003
[12]

Fournais and T.Ø

S. Fournais and T.Ø. Sørensen. Estimates on derivatives of Cou lombic wave functions and their electron densities. J. reine angew. Math. , 2021(775):1–38, 2021

2021
[13]

A.V. Sobolev. Eigenvalue estimates for the one-particle density matrix. J. Spectr. Theory, 12(2):857– 875, 2022

2022
[14]

A.V. Sobolev. Eigenvalue asymptotics for the one-particle dens ity matrix. Duke Math. J. , 171(17):3481–3513, 2022

2022
[15]

A.V. Sobolev. Eigenvalue asymptotics for the one-particle kinet ic energy density operator. Journal of Functional Analysis , 283(8), 2022

2022
[16]

Lieb and R

E.H. Lieb and R. Seiringer. The Stability of Matter in Quantum Mechanics . Cambridge University Press, Cambridge, 2009

2009
[17]

Birman and M.Z

M.S. Birman and M.Z. Solomyak. Estimates of singular numbers of in tegral operators. Russ. Math. Surv., 32(1):15–89, 1977

1977
[18]

Hoﬀmann-Ostenhof, T

M. Hoﬀmann-Ostenhof, T. Hoﬀmann-Ostenhof, and T.Ø. Søre nsen. Electron Wavefunctions and Densities for Atoms. Ann. Henri. Poincare , 2:77–100, 2001

2001
[20]

Birman and Solomjak M.Z

M.S. Birman and Solomjak M.Z. Spectral Theory of Self-Adjoint Operators in Hilbert Space . D. Reidel Publishing Company, 1987

1987
[21]

A. B. Aleksandrov, S. Janson, V. V. Peller, and R. Rochberg. An interesting class of operators with unusual Schatten–von Neumann behavior. In Function sp aces, interpolation theory and related topics. de Gruyter, 2002

2002
[22]

Adams and Fournier J.J.F

R.A. Adams and Fournier J.J.F. Sobolev Spaces. Academic Press, 2 edition, 2003

2003
[23]

C.S. Herz. Lipschitz Spaces and Bernstein’s Theorem on Absolut ely Convergent Fourier Transforms. Journal of Mathematics and Mechanics , 18(4):283–323, 1968

1968
[24]

Bergh and J

J. Bergh and J. L¨ ofstr¨ om.Interpolation Spaces, An Introduction . Springer-Verlag, 1976

1976
[25]

Brasseur

J. Brasseur. On restrictions of Besov functions. Nonlinear Analysis, 170:197–225, 2018

2018
[26]

Mironescu, E

P. Mironescu, E. Russ, and Y. Sire. Lifting in Besov spaces. Nonlinear Analysis, 193, 2020. 111489. Centre for the Mathematics of Quantum Theory, University of Copenhagen, Univer- sitetsparken 5, DK-2100 Copenhagen Ø, Denmark Email address : ph@math.ku.dk

2020

Pith tools

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

[1] [2]

Hearnshaw and A.V

P. Hearnshaw and A.V. Sobolev. Analyticity of the One-Particle De nsity Matrix. Ann. Henri. Poincare, 23:707–738, 2022

work page 2022

[2] [19]

Fournais, M

S. Fournais, M. Hoﬀmann-Ostenhof, T. Hoﬀmann-Ostenhof, and T.Ø. Sørensen. Sharp Regularity Results for Coulombic Many-Electron Wave Functions. Comm. Math. Phys. , 255:183–227, 2005

work page 2005

[3] [1]

T. Kato. On the eigenfunctions of many-particle systems in quan tum mechanics. Comm. Pure Appl. Math., 10:151–177, 1957

work page 1957

[4] [3]

T. Jecko. On the analyticity of electronic reduced densities for m olecules. J. Math. Phys. , 63(1),

work page

[5] [4]

Cioslowski

J. Cioslowski. Oﬀ-diagonal derivative discontinuities in the reduce d density matrices of electronic systems. J. Chem. Phys , 153(154108), 2020

work page 2020

[6] [5]

Cioslowski

J. Cioslowski. Reverse engineering in quantum chemistry: How to r eveal the ﬁfth-order oﬀ-diagonal cusp in the one-electron reduced density matrix without actually ca lculating it. Int. J. Quantum. Chem., 122(8), 2022. e26651

work page 2022

[7] [6]

Hearnshaw and A.V

P. Hearnshaw and A.V. Sobolev. The diagonal behaviour of the on e-particle Coulombic density matrix. Probability and Mathematical Physics , 4(4):935–964, 2023

work page 2023

[8] [7]

Boundedness of the Fifth Derivative for the One-Particle Coulombic Density Matrix at the Diagonal

P. Hearnshaw. Boundedness of the Fifth Derivative for the One -Particle Coulombic Density Matrix at the Diagonal. 2024. arXiv, 2301.02848. EIGENV ALUE BOUNDS FOR REDUCED DENSITY MATRICES 19

work page Pith review arXiv 2024