Edit 26/Sept/13: Fixed Typo in potential
I'm solving the following (seemingly simple) quantum-mechanical problem in four spatial dimensions. In natural units ($\hbar^2/2m=1$), the Schrödinger equation reads:
$$ \Big[-\nabla^2-\frac{24 R^2}{(\mathbf{x}^2+R^2)^2}\Big]\psi(\mathbf{x})=E\,\psi(\mathbf{x})\,, $$
where $R>0$ is a parameter simultaneously characterizing the depth and range of the potential. The potential depends only on the distance away from the origin, so I can separate variables $\psi=R_{nl}(r)\,Y_l(\vec{\theta})$ and the radial equation then reads:
$$ \Big[-\frac{\partial^2}{\partial r^2}-\frac{3}{r}\frac{\partial}{\partial r}+\frac{l(l+2)}{r^2}-\frac{24 R^2}{(r^2+R^2)^2}\Big]\,R_{nl}(r)=E_{nl}\,R_{nl}(r)\,. $$
Problem: I seem to have found an s-wave ($l=0$) non-scattering state with zero energy $E_{nl}=0$ that appears to be localized:
$$R_{n,l=0}(r)=\mathcal{N}\frac{r^2-R^2}{(r^2+R^2)^2}\qquad E_{nl}=0.$$
But, I am unable to normalize this "bound" state since the integral $\int_0^\infty dr\, r^3 |R(r)|^2$ does not converge. What is the nature of this state? Or am I just totally screwing something up?