Para un hidrógeno como el átomo en 3 dimensiones espaciales, la reescritura de la parte radial
$$R(r)~=~\frac{u(r)}{r}$$
no se realiza para mantener el $u(r)$ parte regular, como OP sugiere, pero por lo general debido a que el 3D radial de la ecuación en términos de la $u$ función tiene la misma forma que una 1D ecuación de Schrödinger.
Imagino que la parte radial de la función de onda va como un poder
$$R(r) ~\sim ~ r^{p} \qquad {\rm} \qquad r~\a~ 0,
\qquad p~\~\mathbb{R}.$$
En general jardín, se puede imponer la siguiente lista de consistencia condiciones, que aparece con el más débil condición primera y la más fuerte de la condición última.
Normalizability de la función de onda
$$\infty~>~\langle\psi|\psi\rangle~=~\int d^3r~|\psi(\vec{r})|^2 ~\propto~ \int_0^{\infty} r^{2}dr~|R(r)|^2 .$$
Integrabilidad en $r=0$ de los rendimientos que el poder $p>-\frac{3}{2}$. En otras palabras, este normalizability condición por sí misma no implica que $R(r)$ o $u(r)$ debe ser regular en $r=0$, que es también la conclusión de que muchas de las otras respuestas.
La expectativa de valor de la energía potencial $V$ debe ser delimitada desde abajo,
$$-\infty~<~\langle\psi| V|\psi\rangle~=~\int d^3r~V(r)|\psi(\vec{r})|^2~\propto~-\int_0^{\infty} rdr~|R(r)|^2. $$
Integrabilidad en $r=0$ de los rendimientos que el poder $p>-1$. En otras palabras, $u(r)$ debe ser regular para $r\to 0$.
La energía cinética del operador (o, equivalentemente, el Laplaciano $\Delta$) deberían comportarse de auto-adjointly para dos funciones de onda $\psi_1(\vec{r})$$\psi_2(\vec{r})$,
$$\langle\psi_1| \Delta\psi_2\rangle~=~-\langle\vec{\nabla}\psi_1| \cdot\vec{\nabla}\psi_2\rangle,$$ without picking up pathological contributions at $r=0$. A detailed analysis shows that the powers of the radial parts of $\psi_1(\vec{r})$ and $\psi_2(\vec{r})$ should satisfy $p>-\frac{1}{2}$.
En comparación, el estado unida a soluciones no negativas $p=\ell\in \mathbb{N}_0$, y por lo tanto satisfacer estas tres condiciones.