REVIEW 2 major objections 4 minor 35 references
Barron regularity of many particle Schr\"odinger eigenfunctions
T0 review · 2 major / 4 minor · reviewed 2026-08-15 · deepseek-v4-flash
Pith's one-line read Under a broad Fourier-Lebesgue assumption on the potential, every many-particle Schrödinger eigenfunction lies in the spectral Barron space $B^\gamma$ for all $\gamma < s+2-n/\alpha$, with explicit norm bounds.
desk verdict A genuinely new Barron regularity theorem for many-particle Schrödinger eigenfunctions under Fourier-Lebesgue potential assumptions; the proof is a solid fixed-point argument, with only minor typos to fix. 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 load-bearing machinery is a multiplier estimate in Fourier-Lebesgue spaces: when $V_i, V_{ij} \in \mathcal{F}L^1_s(\mathbb{R}^n) + \mathcal{F}L^{\alpha'}_s(\mathbb{R}^n)$ and $V_{\mathrm{ad}} \in \mathcal{F}L^1_s(\mathbb{R}^{nN})$, multiplication by $V$ is bounded from $\mathcal{F}L^p_{|s|+2\sigma\beta}$ to $\mathcal{F}L^p_{s-2(1-\sigma)\beta}$ with a constant $C(V)$ built from the norms. Combined with the resolvent bound for $(H_0+I)^{-1}$, this gives the operator $T_\lambda$ a gain of two derivative weights, and the high-frequency projector $P_K$ makes $P_KT_\lambda$ contractive. The contraction represents the high-frequency part of $\psi$ by a Neumann series in terms of its bandlimited part, placing $\psi$ in $B^{|s|}$, after which the two-weight gain lifts it to $B^{s+2-2\beta}$.
What would settle it
Take the one-particle radial potential from Example 2.8, $V(x)=\delta^2|x|^{2\delta-2}/2 - \delta(n+\delta-2)|x|^{\delta-2}/2$, whose exact eigenfunction is $e^{-|x|^\delta}$; compute its Fourier transform and check whether the decay rate is exactly $\langle\xi\rangle^{-\delta-n}$. If so, $\psi$ lies in $B^\gamma$ precisely for $\gamma<\delta$ and $\|\psi\|_{B^\gamma}$ blows up like $1/(\delta-\gamma)$ as $\gamma\to\delta$, confirming the claimed sharp index; any faster decay would invalidate the borderline regularity statement.
Extended reading notes
Core claim
On its own terms, the central discovery is Theorem 2.4: under Assumption 2.1, every $H^1$ eigenfunction $\psi$ of $H$ with eigenvalue $\lambda$ lies in the spectral Barron space $B^\gamma(\mathbb{R}^{nN})$ for every $\gamma > |s|$ with $\gamma < s+2-n/\alpha$ when $\alpha < \infty$, and in $B^{s+2}(\mathbb{R}^{nN})$ when $\alpha = \infty$. The theorem supplies two explicit norm bounds, one controlling $\|\psi\|_{B^\gamma}$ by $\|\psi\|_{B^{|s|}}$ and one by $\|\psi\|_{L^2}$. The proof proceeds from the fixed-point identity $\psi = T_\lambda \psi$, where $T_\lambda = (\lambda+1)(H_0+I)^{-1} - (H_0+I)^{-1}V$, and shows that the high-frequency part of $T_\lambda$ is a contraction on $B^{|s|}$. From this theorem the paper recovers Simon's pointwise regularity estimates and, for inverse-power potentials $|x|^{-t}$, obtains $\psi \in B^\gamma$ for $\gamma < 2-t$, which recovers Yserentant's regularity estimate for electronic wave functions with Coulomb potentials.
Load-bearing premise
The results stand on the assumption that the potential splits into one-body, pair, and residual terms whose Fourier transforms lie in a weighted $L^1$ plus $L^{\alpha'}$ space, with the index condition $2+s-|s|-n/\alpha>0$; if a physically relevant potential is more singular than this, such as an inverse power $|x|^{-t}$ with $t\ge 2$ in $n\ge 2$ dimensions, the stated Barron regularity is not obtained.
Editorial extensions
If this is right
- If the theorem is correct, every eigenfunction covered by Assumption 2.1 is a spectral Barron function with index close to $2$, so dimension-free neural network approximation rates follow on bounded domains.
- Taking $\alpha=\infty$ and $s>-1$ gives a pure shift estimate: $V\in B^s$ implies $\psi\in B^{s+2}$, filling the gap left by earlier static Schrödinger solvability results that required $s\ge 0$ and excluded pair interactions.
- The solvability theorem adds existence, uniqueness, and Barron regularity for $(H+\rho I)u=f$ with $f\in B^{\gamma-2}\cap H^{-1}$ when $\rho$ is large enough, under the same singular-potential assumption.
- For Coulomb-type pair potentials, the corollary $\psi\in B^\gamma$ for $\gamma<1$ recovers Yserentant's electronic wave function regularity through a potential-class argument rather than a Coulomb-specific Fourier multiplier.
Reading between the lines
- Beyond the paper, one might expect the same argument to extend to fermionic wave functions if antisymmetry is enforced by Slater determinants; the present statement concerns general $H^1$ eigenfunctions and does not exploit symmetry sectors.
- The borderline case $2+s-|s|-n/\alpha = 0$ is left open; if a regularity index exactly at the boundary could be proven, the divergence in the norm prefactor would disappear.
- Since the potentials are allowed to be unbounded below for $s<0$, the result suggests that sign-definiteness of $V$ is not the relevant condition for Barron regularity; a testable extension is whether more singular local singularities such as $|x|^{-t}$ with $t$ near $2$ in three dimensions can be handled by refining the decomposition.
Signed reviews
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The paper studies the regularity of many-particle Schrödinger eigenfunctions in spectral Barron spaces B^s(R^{nN}). Under Assumption 2.1, in which the one-body and pair potentials lie in F L^1_s+F L^{α'}_s and the additional many-body part lies in F L^1_s with 2+s-|s|-n/α>0, Theorem 2.4 asserts that every H^1 eigenfunction belongs to B^γ for all γ< s+2-n/α (and γ≤s+2 when α=∞), with explicit quantitative bounds. The proof is a contraction fixed-point argument on the high-frequency part of the eigenfunction, built on multiplier estimates in Lemmas 3.4–3.7. The paper also proves solvability of (H+ρI)u=f in Barron spaces under the same assumption (Theorem 2.9) and, for V∈B^s with s>-1, a Fredholm-alternative solvability result (Theorem 2.10). Applications to inverse-power potentials, recovery of Simon's estimate, and a sharpness example are included.
Significance. If the stated results hold, they constitute a substantial advance: they place eigenfunctions for singular potentials such as Coulomb interactions in the function class that guarantees dimension-free neural-network approximation rates, extending earlier Barron-space results that required bounded or nonnegative-index potentials. The assumptions are explicit and are shown to cover inverse-power potentials, the sharpness example gives a useful borderline check, and the paper recovers prior results of Simon and Yserentant as corollaries. The central contraction argument is internally coherent: the exponent bookkeeping for β, γ, s, and α is consistent, and the scope condition 2+s-|s|-n/α>0 is exactly what makes the relevant Hölder-Young estimates available.
major comments (2)
- [Section 3.3, Theorem 2.4 (second displayed estimate)] The displayed upper bound contains an extra factor of 1/(γ-|s|). In the proof, K is chosen so that \tilde μ_1(|λ+1|+C(V))⟨K⟩^{-(γ-|s|)}=1/2; combining (3.17), the B^{|s|} contraction bound, and (3.18) gives \|ψ\|_{B^γ} ≤ 2^{|s|/2+nN/4}√(ω_{nN}/(2|s|+nN)) [2\tilde μ_1(|λ+1|+C(V))]^{(γ+nN/2)/(γ-|s|)} \|ψ\|_{L^2}, with no factor 1/(γ-|s|). The proof as written does not produce the printed denominator, and for γ-|s|>1 the printed estimate is stronger than the derived one. Please correct the constant or supply the missing step.
- [Corollary 2.7] The admissible range for t is stated as t∈(0,3/2+1_{n>1/2}), but this is not the range used in the proof and appears to be a typesetting error. For n≥2 the argument requires t<2 (the decomposition uses s=0 and α< n/t together with α>n/2), and for n=1 it requires t<3/2. Consequently the case t>n in the second displayed estimate can only occur for n=1. The statement should be rephrased with the correct range and the three cases t<n, t=n, t>n should be labeled consistently with the proof.
minor comments (4)
- [Proof of Theorem 2.4] The sentence 'Since α≥1, we have 1/2≤σ≤1' is incorrect for σ=(1-α/2)_+; the correct range is 0≤σ≤1/2. The subsequent conclusion s1+2σβ<1 remains valid, but the printed inequality should be fixed.
- [Equation (3.22)] The constant in (3.22) appears to be missing a factor: from the preceding line one obtains c_{t,n}ω_n = 2π^t Γ((n-t)/2)/(Γ(t/2)Γ(n/2)), whereas (3.22) writes 2^{1-t}Γ((n-t)/2)/(Γ(t/2)Γ(n/2)). Please check and correct the displayed constant.
- [Lemma A.1, proof] In the final estimate of the proof of Lemma A.1, the choice ε=|ξ|^{-1/4} does not give the claimed bound |J_1|≤C|ξ|^{-2δ-n}; the choice should be ε=|ξ|^{-1}. The surrounding argument otherwise supports the statement of the lemma.
- [Throughout] Several typographical errors should be corrected: 'followimg', 'cannonical', 'vatiables', 'Bassron', 'spactral', and the garbled condition in Corollary 2.7 noted above.
Circularity Check
No significant circularity: eigenfunction regularity is derived from the eigenfunction equation and explicit potential hypotheses, not from the target regularity or fitted parameters.
full rationale
The central derivation chain is self-contained. Theorem 2.4 is obtained from the identity ψ=T_λψ, which follows algebraically from Hψ=λψ, combined with explicit convolution and multiplier estimates in Lemmas 3.4–3.7. Assumption 2.1 states the potential hypothesis in Fourier–Lebesgue spaces, and the constants C(V) are explicit in the potential norms; no parameter is fitted to the target regularity and no claim of Barron regularity is assumed at the outset. The fixed-point representation (3.17) gives ψ∈B^{|s|}, and Lemma 3.6 then lifts to B^{s+2−2β}, which is exactly the claimed range after choosing β=1+(s−γ)/2. The sharpness example is an independently computed eigenfunction e^{−|x|^δ} whose Fourier decay is established directly in Lemma A.1, so it does not borrow the conclusion. Self-citations [15], [22], and [23] are used for background, standard Barron-space embedding facts, and downstream neural-network approximation rates; none is load-bearing for the eigenfunction regularity proof. The noted proof typos, such as the claimed range 1/2≤σ≤1 and the exponent ε=|ξ|^{−1/4} in Lemma A.1, are correctness issues in secondary details and do not constitute circularity.
Assumptions & free parameters
assumptions (5)
- standard math Fourier transform on tempered distributions S'(R^d) and standard inequalities (Hölder, Young, Parseval)
- standard math Lax-Milgram theorem
- standard math Fredholm alternative for compact operators on Banach spaces
- standard math Kolmogorov-Riesz compactness criterion
- domain assumption Potential decomposition (1.1) with V_i,V_ij in F L^1_s+F L^{alpha'}_s and V_ad in F L^1_s (Assumption 2.1)
Cite this review
Pith. "Pith review of Barron regularity of many particle Schr\"odinger eigenfunctions." pith.science (2026). https://pith.science/paper/MDCGWTPF
@misc{pith2026250817722,
author = {Pith},
title = {Pith review of: Barron regularity of many particle Schr\"odinger eigenfunctions},
year = {2026},
howpublished = {\url{https://pith.science/paper/MDCGWTPF}},
note = {Machine review of arXiv:2508.17722}
}
abstract
This work investigates the regularity of Schr\"odinger eigenfunctions and the solvability of Schr\"odinger equations in spectral Barron space $\mathcal{B}^{s}(\mathbb{R}^{nN})$, where neural networks exhibit dimension-free approximation capabilities. Under assumptions that the potential $V$ consists of one-particle and pairwise interaction parts $V_{i},V_{ij}$ in Fourier-Lebesgue space $\mathcal{F}L_{s}^{1}(\mathbb{R}^{n})+\mathcal{F}L_{s}^{\alpha^{\prime}}(\mathbb{R}^{n})$ and an additional part $V_{\operatorname{a d}} \in \mathcal{F}L_{s}^{1}(\mathbb{R}^{nN})$, we prove that all eigenfunctions $\psi\in \bigcap_{\gamma<s+2-n/\alpha} \mathcal{B}^{\gamma}(\mathbb{R}^{nN})$ and $\psi\in \mathcal{B}^{s+2}(\mathbb{R}^{nN})$ if $\alpha=\infty$, where $1/\alpha+1/\alpha^{\prime}=1$ and $2+s-|s|-n/\alpha>0$. The assumption accommodates many prevalent singular potentials, such as inverse power potentials. Moreover, under the same assumption or a stronger assumption $V\in\mathcal{B}^{s}(\mathbb{R}^{nN})$, we establish the solvability of Schr\"odinger equations and derive compactness results for $V\in\mathcal{B}^{s}(\mathbb{R}^{nN})$ with $s>-1$.
Figures
Reference graph
Works this paper leans on
-
[15]
Y. Guo, P. Ming, and H. Yu,Generalization error estimates of machine learning methods for solving high dimensional Schr¨odinger eigenvalue problems, accepted. arXiv:2408.13511
-
[22]
Y. Liao and P. Ming, Spectral Barron space for deep neural network approximation, SIAM J. Math. Data Sci. 7 (2025), no. 3, 1053–1076
work page 2025
-
[23]
Y. Liao, P. Ming, and H. Yu, Sharp uniform approximation for spectral Barron functions by deep neural networks , 2025. arXiv:2507.06789
work page Pith review arXiv 2025
-
[1]
A. R. Barron, Neural net approximation, Proc. 7th Yale Workshop on Adaptive and Learning Systems1 (1992), 69–72
work page 1992
-
[2]
, Universal approximation bounds for superpositions of a sigmoidal function , IEEE Trans. Inform. Theory 39 (1993), no. 3, 930–945
work page 1993
-
[3]
, Approximation and estimation bounds for artificial neural networks, Mach. Learn. 14 (1994), 115–133
work page 1994
- [4]
-
[5]
G. Carleo and M. Troyer, Solving the quantum many-body problem with artificial neural networks, Science 355 (2017), no. 6325, 602–606
work page 2017
Show all 35 references
-
[6]
Z. Chen, J. Lu, Y. Lu, and S. Zhou, A regularity theory for static Schr¨odinger equations on R𝑑 in spectral Barron spaces, SIAM J. Math. Anal. 55 (2023), no. 1, 557–570
2023
-
[7]
K. Choo, A. Mezzacapo, and G. Carleo, Fermionic neural-network states for ab-initio electronic structure , Nat. Commun. 11 (2020), no. 1, 2368
2020
-
[8]
Dus and E
M. Dus and E. Virginie, Two-layers neural networks for Schr¨odinger eigenvalue problems, 2024. arXiv:2409.01640
2024 arXiv
-
[9]
W. E and S. Wojtowytsch, Some observations on high-dimensional partial differential equations with Barron data, Proceedings of the 2nd mathematical and scientific machine learning conference, 2022, pp. 253–269
2022
-
[10]
Feng and J
Y. Feng and J. Lu, Solution theory of Hamilton-Jacobi-Bellman equations in spectral Barron spaces, 2025. arXiv:2503.18656
2025
-
[11]
Folland, Introduction to Partial Differential Equations, Princeton University Press, Princeton New Jersey, 1995
G.B. Folland, Introduction to Partial Differential Equations, Princeton University Press, Princeton New Jersey, 1995. 2nd ed
1995
-
[12]
Fournais, M
S. Fournais, M. Hoffmann-Ostenhof, T. Hoffmann-Ostenhof, and T. Ø stergaard Sø rensen,Sharp regularity results for Coulombic many-electron wave functions, Comm. Math. Phys. 255 (2005), no. 1, 183–227
2005
-
[13]
, Analytic structure of many-body Coulombic wave functions, Comm. Math. Phys. 289 (2009), no. 1, 291–310
2009
-
[14]
Gao and L
X. Gao and L. Duan, Efficient representation of quantum many-body states with deep neural networks, Nat. Commun. 8 (2017), no. 1, 662
2017
-
[16]
J. Han, J. Lu, and M. Zhou, Solving high-dimensional eigenvalue problems using deep neural networks: a diffusion Monte Carlo like approach, J. Comput. Phys. 423 (2020), 13
2020
-
[17]
J. Han, L. Zhang, and W. E, Solving many-electron Schr ¨odinger equation using deep neural networks , J. Comput. Phys. 399 (2019), 8. 26 P. B. MING AND H. YU
2019
-
[18]
Hermann, Z
J. Hermann, Z. Sch ¨atzle, and F. No´e, Deep-neural-network solution of the electronic Schr¨odinger equation, Nat. Chem.12 (2020), no. 10, 891–897
2020
-
[19]
H ¨ormander, The Analysis of Linear Partial Differential Operators
L. H ¨ormander, The Analysis of Linear Partial Differential Operators. II , Vol. 257, Springer-Verlag, Berlin, 1983. Differential operators with constant coefficients
1983
-
[20]
Kato, On the eigenfunctions of many-particle systems in quantum mechanics , Comm
T. Kato, On the eigenfunctions of many-particle systems in quantum mechanics , Comm. Pure Appl. Math. 10 (1957), no. 2, 151–177
1957
-
[21]
Li and L
H. Li and L. Ying, A semigroup method for high dimensional elliptic PDEs and eigenvalue problems based on neural networks, J. Comput. Phys. 453 (2022), 15
2022
-
[24]
Lu and Y
J. Lu and Y. Lu, A priori generalization error analysis of two-layer neural networks for solving high dimensional Schr ¨odinger eigenvalue problems, Commun. Am. Math. Soc. 2 (2022), 1–21
2022
-
[25]
J. Lu, Z. Shen, H. Yang, and S. Zhang, Deep network approximation for smooth functions, SIAM J. Math. Anal. 53 (2021), no. 5, 5465–5506
2021
-
[26]
Y. Lu, H. Chen, J. Lu, L. Ying, and J. Blanchet, Machine learning for elliptic PDEs: Fast rate generalization bound, neural scaling law and minimax optimality, International conference on learning representations, 2022
2022
-
[27]
Y. Lu, J. Lu, and M. Wang, A priori generalization analysis of the deep Ritz method for solving high dimensional elliptic partial differential equations, Proceedings of thirty fourth conference on learning theory, 2021, pp. 3196–3241
2021
-
[28]
Simon, Pointwise bounds on eigenfunctions and wave packets in 𝑁-body quantum systems
B. Simon, Pointwise bounds on eigenfunctions and wave packets in 𝑁-body quantum systems. I , Proc. Amer. Math. Soc. 42 (1974), 395–401
1974
-
[29]
Stein and G
E. Stein and G. Weiss, Fourier Analysis on Euclidean Spaces, Princeton University Press, 1971
1971
-
[30]
J. Toft, K. Johansson, S. Pilipovi ´c, and N. Teofanov, Sharp convolution and multiplication estimates in weighted spaces , Anal. Appl. (Singap.) 13 (2015), no. 5, 457–480
2015
-
[31]
Wang and H
Y. Wang and H. Xie, Computing multi-eigenpairs of high-dimensional eigenvalue problems using tensor neural networks , J. Comput. Phys. 506 (2024), 112928
2024
-
[32]
Yarotsky, Error bounds for approximations with deep ReLU networks, Neural networks 94 (2017), 103–114
D. Yarotsky, Error bounds for approximations with deep ReLU networks, Neural networks 94 (2017), 103–114
2017
-
[33]
Yserentant, On the regularity of the electronic Schr¨odinger equation in Hilbert spaces of mixed derivatives, Numer
H. Yserentant, On the regularity of the electronic Schr¨odinger equation in Hilbert spaces of mixed derivatives, Numer. Math. 98 (2004), no. 4, 731–759
2004
-
[34]
, Sparse grid spaces for the numerical solution of the electronic Schr¨odinger equation, Numer. Math. 101 (2005), no. 2, 381–389
2005
-
[35]
arXiv:2307.15285
, The regularity of electronic wave functions in Barron spaces, 2025. arXiv:2307.15285. SKLMS, Institute of Computational Mathematics and Scientific/Engineering Computing, AMSS, Chinese Academy of Sciences, Beijing 100190, China School of Mathematical Sciences, University of C...
2025 arXiv
Reviewed August 15, 2026 · model on record in the stance chip above.
Discussion (0). Continue with ORCID to comment.