What is the expectation value of the 3D delta function for the Hydrogen atom ground state?

Question

I'm trying to evaluate the expectation value of some perturbation Hamiltonian $$H=\alpha \delta^3(\vec{r}),$$ where $\alpha$ is a positive constant, for the ground state wavefunction of the hydrogen atom $$\psi_{100}~\propto~ \exp[-r/a]$$ (I want to calculate the shift in the energy in first order perturbation theory). Why is it wrong to write:

$$\langle\psi_{100}|H|\psi_{100}\rangle ~\propto~ \int_0^\infty dr~ r^2 \exp[-2r/a] \alpha \delta(r) ~=~ 0$$

Does it have to do with either:

1. I can't just write $\delta^3(\vec{r}) = \delta(r)$ in the integral over $r$, or

2. I can't evaluate the delta function at zero, since it is at the endpoint of the integration limits [not inside the interval $(0,\infty)$]?

Answer 1

Hint: Rather than using spherical coordinates, which are singular where the 3D Dirac delta function has support, work instead in Cartesian coordinates $\vec{r}=(x,y,z)$ and use the defining property

$$ \iiint_{\mathbb{R}^3}\! d^3r ~f(\vec{r})~ \delta^3(\vec{r})~=~f(\vec{0}) $$

of the 3D Dirac delta function.

Answer 2

I want to post my attempt at the solution based on Qmechanic's hint (thanks!):

Rewriting in cartesian coordinates

$$ \psi_{100} \propto \exp[-\sqrt{x^2+y^2+z^2}/a] $$

$$ <\psi_{100}|H|\psi_{100}> \propto \int_{-\infty}^{+\infty} dx \int_{-\infty}^{+\infty} dy \int_{-\infty}^{+\infty} dz \exp\left[-2\sqrt{x^2+y^2+z^2}/a\right] \alpha \delta(x)\delta(y)\delta(z) =\alpha$$

Hope thats correct, and sorry if this was trivial, certainly the expected answer.

I still don't fully understand why it's not possible to directly solve this in spherical coordinates, but I understand it is related to the difficulty in evaluating $\delta(r)$ at the endpoint $r=0$ (see CuriousKev's comment).