Edit: I've tried to make the question a lot simpler.
The Laplacian of the function $1/r$ is $-4\pi\delta(r)$. If I take the trace of the Hessian matrix of the function $1/r$, I find that the trace is 0 even though I expect I should find $-4\pi\delta(r)$. I understand this could be happening because the Dirac delta function isn't a function, but a measure.
Now, my questions are as follows:
- In physics, one way of defining the delta function is by appealing to Poisson's equation as is done in Jackson's E&M in the link below. Is it possible to do something similar to get a matrix analogue of the Dirac delta function, such that the trace of this matrix analogue is $-4\pi\delta(r)$?
- If the above is not possible, should this be concerning? In other words, if I can't recover the density by taking the trace of a Hessian matrix of a potential, what is that telling me about this procedure?
Link, since I can't embed pictures yet
Sorry if this question isn't well-phrased. My background is not in math.