{"id":"315cfc98-5aa2-41f1-9e9d-f6216994e98d","arxiv_id":"2412.16073","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":0,"one_line_summary":"For bound states of N-electron atoms and molecules, the eigenvalues of the K-particle reduced density matrix decay as n^{-(1+7/(3 min{K,N-K}))}.","lead":"This paper proves new decay rates for the eigenvalues of multi-particle reduced density matrices of electrons in atoms and molecules. It fills a gap between previously known cases, giving a quantitative measure of how well these matrices can be approximated by finite-rank ones.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The proof of Theorem 1.2 only derives exponent 1+7/(3K); for K>N/2 it never applies Proposition 1.4 with K replaced by N-K, so the stated exponent 1+7/(3 min{K,N-K}) is not actually proven.","rationale":"The main theorem is plausible and the technical machinery is coherent: the Jastrow factorization in (2.10), the Besov regularity input in Lemma 4.5, the Birman-Solomyak bound in Proposition 3.2, and the multiplier result in Lemma 5.4 all line up to give Proposition 1.4. I checked that (1.13) correctly identifies singular values of Ψ^(K) with eigenvalues of both Γ^(K) and Γ^(N-K) up to permutation, and that the summability of the weights A and B in Proposition 6.1 is fine. The only place where the written proof does not match the theorem statement is the exponent for K>N/2: the proof as printed derives the bound with 1+7/(3K), not with 1+7/(3 min{K,N-K}). The reader's weakest-assumption concern about Proposition 2.1 is reasonable as a dependence on external regularity results, but I do not see an internal inconsistency there; the Appendix gives the main steps and cites the needed machinery. Since the missing symmetry step is explicit and easy to repair, the appropriate outcome is the same conditional verdict rather than a change to accept or reject.","tokens_in":16904,"tokens_out":28159,"duration_ms":238026,"concrete_test":"Take N=5, K=3 and trace the proof of Theorem 1.2 literally. The written argument gives 1/q = 1/2 + 7/18 = 8/9, hence λ_n(Γ^(3)) ≤ C n^{-16/9}, whereas the theorem claims C n^{-13/6}. To settle whether the gap is real, check whether the proof anywhere applies Proposition 1.4 with K' = N-K = 2 (so 1/q' = 1/2 + 7/12 = 13/12) and uses (1.13) to conclude λ_n(Γ^(3)) ≤ C n^{-13/6}. If it does not, the gap is confirmed; adding this explicit symmetry step repairs the proof.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The load-bearing gap is in the proof of Theorem 1.2. The written argument applies Proposition 1.4 once, to the given K, and concludes λ_n(Γ^(K)) ≤ C n^{-(1+7/(3K))}. For K > N/2 this is weaker than the stated α_K = 1+7/(3(N-K)). The following sentence about taking any permutation only replaces Γ^(N-K)_σK by Γ^(N-K); it does not change the exponent in the bound. The missing step is to apply Proposition 1.4 with K replaced by N-K (which is also between 2 and N-2) and use (1.13) to transfer the resulting S_{q',∞} bound to Γ^(K). All ingredients are already present, so the gap is easily repairable, but as printed the central theorem is not established for half of the stated range. I found no comparable flaw in the Besov, Birman-Solomyak, or multiplier lemmas; the reliance on Proposition 2.1 is external but plausible and is not the clearest textual gap.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","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.","tokens_in":17096,"tokens_out":8716,"duration_ms":68106,"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":[{"comment":"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.","section":"Section 1.1, proof of Theorem 1.2"}],"minor_comments":[{"comment":"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":"Appendix A, proof of Proposition 2.1"},{"comment":"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":"Section 5, Lemma 5.3"},{"comment":"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.","section":"Section 1.1, equation (1.13)"}],"recommendation":"major_revision","confidential_remarks":"The gap in the proof of Theorem 1.2 is local and straightforward to repair: applying Proposition 1.4 with K replaced by N-K and using (1.13) should yield the stated min{K,N-K} exponent. The revision will be quick if the author adds this symmetric argument. The reliance on Proposition 2.1 is reasonable given the published sources, but the journal may request a more detailed proof of that proposition. Overall the paper is a good fit for the journal's scope."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Dear [Colleague],\n\nQuick 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.\n\nThe 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.\n\nI 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.\n\nBottom 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.\n\nBest,\n[Your name]","headline":"Solid new result with a small repairable gap in Theorem 1.2: the symmetry step for K>N/2 is stated but never executed.","tokens_in":17691,"tokens_out":2652,"would_cite":true,"duration_ms":21843,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":["35J10","46E35","81V55","81V70"],"pacs":[],"model":"deepseek-v4-flash","headline":"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}.","keywords":["reduced density matrix","Coulombic wavefunction","eigenvalue decay","Schatten classes","Besov-Nikol'skii spaces","Jastrow factor","Birman-Solomyak bounds","many-electron Schrödinger equation"],"falsifier":"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}).","tokens_in":16644,"feed_emoji":"⚛️","tokens_out":9506,"duration_ms":76564,"temperature":0.7,"pith_summary":"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.","feed_headline":"Every K-particle density matrix now proven to decay faster than n^-1","feed_subtitle":"Coulombic bound states with N electrons now carry eigenvalue bounds λ_n ≤ C n^{-α} with α > 1 for every K between 2 and N−2.","key_machinery":"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.","core_discovery":"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.","pith_inferences":["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."],"forward_implications":["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."],"supporting_citations":[{"why":"It supplies the earlier one-particle and (N−1)-particle eigenvalue bounds that this paper extends to all intermediate K.","marker":"[13]"},{"why":"It provides the Birman–Solomyak singular-value estimates for integral operators with Besov-smooth kernels, used to place Ψ^(K) in weak Schatten classes.","marker":"[17]"},{"why":"It gives the W^{2,∞} regularity for the optimized Jastrow wavefunction from which the order-two derivative bounds on φ are obtained.","marker":"[19]"},{"why":"It supplies pointwise derivative bounds and regularity machinery combined with [19] to prove Proposition 2.1.","marker":"[2]"},{"why":"It introduces the Jastrow-factor method and derivative estimates for Coulombic wavefunctions on which the factorization and Besov regularity arguments build.","marker":"[12]"},{"why":"It provides the Besov-space criterion and the fifth-order-cusp analysis used to localize μ and to argue the exponent cannot likely beat the K=1 rate.","marker":"[6]"},{"why":"It justifies the exponential-decay assumption by showing it holds for discrete eigenvalues of the Coulomb Hamiltonian.","marker":"[10]"},{"why":"It records the formula recovering energy from the two-particle density matrix, used to motivate why the K=2 bound matters.","marker":"[16]"}],"fun_headline_variants":["Coulombic eigenvalues decay faster than n^-1 for all K","Stronger decay for Coulombic K-particle density matrices","Coulombic bound states: eigenvalue decay exponent boosted","New bounds: K-particle matrices decay faster for Coulombic","K-particle density matrices: faster eigenvalue decay proven"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"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.","fun_headline_variants_meta":{"raw":{"variants":["Coulombic eigenvalues decay faster than n^-1 for all K","Stronger decay for Coulombic K-particle density matrices","Coulombic bound states: eigenvalue decay exponent boosted","New bounds: K-particle matrices decay faster for Coulombic","K-particle density matrices: faster eigenvalue decay proven"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.001411,"raw_usage":{"total_tokens":5698,"prompt_tokens":943,"completion_tokens":4755,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":559,"completion_tokens_details":{"reasoning_tokens":4676}},"tokens_in":559,"tokens_out":4755,"duration_ms":32362,"temperature":1.0,"reasoning_tokens":4676,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-11T10:49:51.765986+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"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}).","supporting_citations":[{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"It supplies the earlier one-particle and (N−1)-particle eigenvalue bounds that this paper extends to all intermediate K."},{"cited_title":"Birman and M.Z","cited_arxiv_id":null,"evidence_quote":"It provides the Birman–Solomyak singular-value estimates for integral operators with Besov-smooth kernels, used to place Ψ^(K) in weak Schatten classes."},{"cited_title":"Fournais, M","cited_arxiv_id":null,"evidence_quote":"It gives the W^{2,∞} regularity for the optimized Jastrow wavefunction from which the order-two derivative bounds on φ are obtained."},{"cited_title":"Hearnshaw and A.V","cited_arxiv_id":null,"evidence_quote":"It supplies pointwise derivative bounds and regularity machinery combined with [19] to prove Proposition 2.1."},{"cited_title":"Fournais and T.Ø","cited_arxiv_id":null,"evidence_quote":"It introduces the Jastrow-factor method and derivative estimates for Coulombic wavefunctions on which the factorization and Besov regularity arguments build."},{"cited_title":"Hearnshaw and A.V","cited_arxiv_id":null,"evidence_quote":"It provides the Besov-space criterion and the fifth-order-cusp analysis used to localize μ and to argue the exponent cannot likely beat the K=1 rate."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"It justifies the exponential-decay assumption by showing it holds for discrete eigenvalues of the Coulomb Hamiltonian."},{"cited_title":"Lieb and R","cited_arxiv_id":null,"evidence_quote":"It records the formula recovering energy from the two-particle density matrix, used to motivate why the K=2 bound matters."}],"review_version":1}