Pith. sign in

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 →

arxiv 2508.17722 v1 pith:MDCGWTPF submitted 2025-08-25 math.AP

classification math.AP MSC 35J1035B6535Q4068T07
keywords SchrödingerequationBarronspaceregularitytheoryeigenfunctionFourier-Lebesguespacesinversepowerpotentialsneuralnetworkapproximationmany-particlesystems
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 establishes that eigenfunctions of the many-particle Schrödinger operator $H = -\sum_i (1/2\mu_i)\Delta_i + V$ inherit regularity in the spectral Barron scale, provided each one-body and pairwise interaction term lies in a weighted Fourier-Lebesgue sum space rather than a plain Barron space. The potential class includes Coulomb and other inverse-power singularities, which earlier Barron-space regularity results excluded. If true, the eigenfunctions are guaranteed to be approximable by neural networks at rates independent of the high dimension $nN$, exactly the property machine-learning solvers need. The same arguments give existence, uniqueness, and Barron regularity of solutions to $(H+\rho I)u=f$ under the same assumption.

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.

Watch

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

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

  • 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.
Share X Bluesky LinkedIn Reddit HN

Signed reviews

No signed human review yet.

Editorial analysis

A structured set of objections, weighed in public.

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

Referee Report

2 major / 4 minor

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)
  1. [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.
  2. [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)
  1. [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.
  2. [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.
  3. [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.
  4. [Throughout] Several typographical errors should be corrected: 'followimg', 'cannonical', 'vatiables', 'Bassron', 'spactral', and the garbled condition in Corollary 2.7 noted above.

Circularity Check

0 steps flagged · score 0.0 of 10

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 0 free parameters · 5 assumptions · 0 invented entities

The paper introduces no fitted numbers: the constants C(V), C(V), and the regularity index s+2-n/alpha are derived from the hypothesis rather than tuned to data. No new physical entities are introduced. The results rest on standard harmonic analysis, functional analysis (Lax-Milgram, Fredholm, Kolmogorov-Riesz), and the explicit structural Assumption 2.1 on the potential.

assumptions (5)
  • standard math Fourier transform on tempered distributions S'(R^d) and standard inequalities (Hölder, Young, Parseval)
    Used throughout Sections 3 and 4 for all convolution and multiplier estimates, e.g., Lemmas 3.2, 3.4, 3.5.
  • standard math Lax-Milgram theorem
    Used in Theorem 2.9 to obtain existence and uniqueness of the weak solution in H^1.
  • standard math Fredholm alternative for compact operators on Banach spaces
    Used in Theorem 2.10 after proving compactness of R on B^{s+2} in Proposition 4.4.
  • standard math Kolmogorov-Riesz compactness criterion
    Used in Proposition 4.4 to verify total boundedness of the image set G.
  • 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)
    The entire theory is conditional on this hypothesis; Lemma 2.2 verifies it for inverse-power potentials, the intended application.

how reviews work

0 comments
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

Figures reproduced from arXiv: 2508.17722 by the authors.

Figure 1
Figure 1. The structure diagram of main results in this paper with the explicit norm estimates. The case 𝛼 = ∞ reduces to a shift estimate in Barron space; see Theorem 2.5. We recover Simon’s estimate [28, Theorem 1′ ] from Theorem 2.4 with 𝑠 = 0; see Theorem 2.6. Then, as an application of Theorem 2.6, when 𝑉𝑖 , 𝑉𝑖 𝑗 are inverse power potentials |𝑥| −𝑡 , we obtain 𝜓 ∈ B𝛾 (R 𝑁 𝑛) for 𝛾 < 2−𝑡 and naturally recover Yserentant’s… view at source ↗
Figure 2
Figure 2. Left: Embedding relations in Proposition 2.3. Each point (𝛼, 𝑠) represents a space F 𝐿 1 𝑠 (R 𝑛 ) + F 𝐿 𝛼 ′ 𝑠 (R 𝑛 ). Right: Graphical interpretation for the regularity results of eigenfunctions. The shaded area corresponds to Assumption 2.1. Potentials on each dashed line leads to the same regularity of eigenfunctions2. Example 2.8 shows the sharpness of our regularity estimates in the orange area. We postpone the … view at source ↗

Discussion (0). Continue with ORCID to comment.

Reference graph

Works this paper leans on

35 extracted references · 32 canonical work pages

  1. [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

  2. [22]

    Liao and P

    Y. Liao and P. Ming, Spectral Barron space for deep neural network approximation, SIAM J. Math. Data Sci. 7 (2025), no. 3, 1053–1076

  3. [23]

    Y. Liao, P. Ming, and H. Yu, Sharp uniform approximation for spectral Barron functions by deep neural networks , 2025. arXiv:2507.06789

  4. [1]

    A. R. Barron, Neural net approximation, Proc. 7th Yale Workshop on Adaptive and Learning Systems1 (1992), 69–72

  5. [2]

    , Universal approximation bounds for superpositions of a sigmoidal function , IEEE Trans. Inform. Theory 39 (1993), no. 3, 930–945

  6. [3]

    , Approximation and estimation bounds for artificial neural networks, Mach. Learn. 14 (1994), 115–133

  7. [4]

    Cai and J

    Z. Cai and J. Liu, Approximating quantum many-body wave functions using artificial neural networks, Phys. Rev. B 97 (2018), 035116

  8. [5]

    Carleo and M

    G. Carleo and M. Troyer, Solving the quantum many-body problem with artificial neural networks, Science 355 (2017), no. 6325, 602–606

Show all 35 references
  1. [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

  2. [7]

    K. Choo, A. Mezzacapo, and G. Carleo, Fermionic neural-network states for ab-initio electronic structure , Nat. Commun. 11 (2020), no. 1, 2368

  3. [8]

    Dus and E

    M. Dus and E. Virginie, Two-layers neural networks for Schr¨odinger eigenvalue problems, 2024. arXiv:2409.01640

  4. [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

  5. [10]

    Feng and J

    Y. Feng and J. Lu, Solution theory of Hamilton-Jacobi-Bellman equations in spectral Barron spaces, 2025. arXiv:2503.18656

  6. [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

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

  8. [13]

    , Analytic structure of many-body Coulombic wave functions, Comm. Math. Phys. 289 (2009), no. 1, 291–310

  9. [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

  10. [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

  11. [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

  12. [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

  13. [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

  14. [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

  15. [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

  16. [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

  17. [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

  18. [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

  19. [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

  20. [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

  21. [29]

    Stein and G

    E. Stein and G. Weiss, Fourier Analysis on Euclidean Spaces, Princeton University Press, 1971

  22. [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

  23. [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

  24. [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

  25. [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

  26. [34]

    , Sparse grid spaces for the numerical solution of the electronic Schr¨odinger equation, Numer. Math. 101 (2005), no. 2, 381–389

  27. [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...

Pith tools

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