Correct radial operator

Consider the scalar operator H=−Δ+V(r)H=-\Delta+V(r) on ordinary Rd\mathbb R^d, d≥2d\ge2, in units ℏ2/(2m)=1\hbar^2/(2m)=1. For a normalized spherical harmonic with angular eigenvalue ℓ(ℓ+d−2)\ell(\ell+d-2), write ψ(x)=R(r)Yℓ(x^)\psi(x)=R(r)Y_\ell(\hat x) and u(r)=r(d−1)/2R(r)u(r)=r^{(d-1)/2}R(r). This transformation sends the radial norm to ∫0∞∣u∣2dr\int_0^\infty|u|^2dr. Differentiating R=r−(d−1)/2uR=r^{-(d-1)/2}u yields

Hℓu=−u′′+[V(r)+(ℓ+(d−2)/2)2−1/4r2]u. H_\ell u=-u''+ \left[V(r)+\frac{(\ell+(d-2)/2)^2-1/4}{r^2}\right]u.

The term (d−1)(d−3)/(4r2)(d-1)(d-3)/(4r^2) is contributed by the radial measure. Use the Friedrichs realization inherited from the full spatial quadratic form. In particular the regular branch has u∼rℓ+(d−1)/2u\sim r^{\ell+(d-1)/2}; the critical d=2,ℓ=0d=2,\ell=0 endpoint selects the regular, nonlogarithmic branch. An arbitrary half-line boundary condition describes a different operator.

A form bound on attraction

For d≥3d\ge3, integration by parts for smooth compactly supported ψ\psi gives

0≤∫∣∇ψ+axr2ψ∣2dx=∫∣∇ψ∣2dx+[a2−a(d−2)]∫∣ψ∣2r2dx. 0\le\int\left|\nabla\psi+a\frac{x}{r^2}\psi\right|^2dx =\int|\nabla\psi|^2dx+ [a^2-a(d-2)]\int\frac{|\psi|^2}{r^2}dx.

Choose a=(d−2)/2a=(d-2)/2 and extend by form closure:

∫∣∇ψ∣2dx≥(d−2)24∫∣ψ∣2r2dx. \boxed{\int|\nabla\psi|^2dx\ge \frac{(d-2)^2}{4}\int\frac{|\psi|^2}{r^2}dx.}

This is Hardy's inequality; see Nonnenmacher, Proposition 3.2.7.

As an application, measure lengths in a stated range LL and energies in ℏ2/(2mL2)\hbar^2/(2mL^2). For the dimensionless operator

Hsat=−Δ+λr21+r2,λ>0, H_{\rm sat}=-\Delta+\lambda\frac{r^2}{1+r^2},\qquad \lambda>0,

subtract its limiting potential λ\lambda. Since (1+r2)−1≤r−2(1+r^2)^{-1}\le r^{-2},

⟨Hsat−λ⟩≥[(d−2)24−λ]∫∣ψ∣2r2dx. \langle H_{\rm sat}-\lambda\rangle \ge\left[\frac{(d-2)^2}{4}-\lambda\right] \int\frac{|\psi|^2}{r^2}dx.

Thus λ≤(d−2)2/4\lambda\le(d-2)^2/4 excludes every spectral value below λ\lambda. The bound is sufficient; exceeding it alone does not establish a bound state. A first-order Dirac operator needs its own analysis.

resonance_profile_checks.py checks finite examples by quadrature, matrix evolution, and independent algebra.