We establish sharp Barron regularity for Coulombic many-electron wave functions after extraction of the universal cut-off Jastrow factors. Following the factorization of Fournais et al.~\cite[Definition~1.4]{FournaisEtAl2005}, for a Coulombic eigenfunction $ψ$ we define the successive quotients by \[ φ=e^{-F_{2,\mathrm{cut}}}ψ\quad\text{and}\quad φ_3=e^{-F_{3,\mathrm{cut}}}φ=e^{-(F_{2,\mathrm{cut}}+F_{3,\mathrm{cut}})}ψ. \] Then \[ φ,φ_3\in\mathcal{B}^s(\mathbb{R}^{3N}) \qquad\text{for every }s<2. \] This range is optimal among universal factorizations. No factor depending only on the particle number and the nuclear data, but not on the eigenfunction or its eigenvalue, can make every corresponding quotient belong to $\mathcal{B}^2$. We also determine the exact endpoint growth. Writing $\varepsilon=2-s$, we prove that, for either $u=φ$ or $u=φ_3$, there is a computable constant $M$ independent of $\varepsilon$ such that \[ \left\|u\right\|_{\mathcal{B}^{2-\varepsilon}}\leq\frac{M}{\varepsilon^2}\left\|u\right\|_{\mathcal{B}^1}. \] For the unperturbed two-electron atom we prove, with a constant independent of $\varepsilon$, \[ \left|\left\|φ_3\right\|_{\mathcal{B}^{2-\varepsilon}}-\frac{32πZ\lvertφ_3(0,0)\rvert}{\varepsilon^2}\right|\leq\frac{C}{\varepsilon}. \] Hence the quadratic rate in the upper bound is sharp whenever $\lvertφ_3(0,0)\rvert\neq0$, as is the case for the ground state.
We prove that harmonic functions with Dirichlet boundary data in Barron space, a function class tailored to wide ReLU networks with a single hidden layer and suitably bounded weights, are generally neither Lipschitz continuous nor in the Sobolev class $H^2$. A fortiori, they are not in any function class in which the norm controls the Lipschitz constant, which rules out not only Barron space regularity, but also regularity in function classes for deeper ReLU networks with bounded coefficients. They can, however, be approximated to accuracy $\sim \varepsilon$ by Barron functions of low norm $\sim |\log\varepsilon|$ in various Lebesgue and Sobolev norms (with at most two derivatives). The positive result holds on very simple domains: Half-spaces in arbitrary dimension and rectangular domains in two dimensions. As an application of this regularity theory, we obtain a priori error estimates for Deep Ritz neural PDE solvers.