# Surface area of sphere using Dirac delta

This question is related to this one.

Suppose I want to calculate the surface area $S(R)$ of a sphere of radius $R$. I can express $S(R)$ as

$$S(R)=\int_{\mathbb{R}^3} \delta (\| \vec x \|-R) \ d \vec x$$

I would then obtain

$$S(R) = 4 \pi \int_0^\infty \delta(r-R) r^2 dr = 4 \pi R^2$$

which is the correct result.

But it seems to me that I could equivalently express $S(R)$ as

$$S(R)=\int_{\mathbb{R}^3} \delta (\| \vec x \|^2-R^2) \ d \vec x$$

which gives

$$S(R) = 4 \pi \int_0^\infty \delta(r^2-R^2) r^2 dr$$

From the property of composition of the delta with a function,

$$\delta(r^2-R^2)=\frac{\delta(r-R)+\delta(r+R)}{2R}$$

but since $r \geq 0$ I only have to consider the positive root, so that

$$S(R) = 4 \pi \int_0^\infty \frac{\delta(r-R)}{2R} r^2 dr = 2 \pi R$$

Why do I get two different results? Is something wrong with the second way of expressing $S(R)$?

• What exactly is this question asking that the other one isn't? Jul 19, 2016 at 13:14
• @DavidZ The other question was based on an error in the calculation of $\int \delta(r^2-R^2) r^2 dr$. In this question I'm asking why two expression that seem to be equivalent to me hold different results. I know that the two questions look the same, but they aren't. Jul 19, 2016 at 13:16
• Well, but that seems like the same thing the other question was asking. In the other question, you calculated the same two expressions and obtained the same two different results, and you asked which one is right (and implicitly why). I don't see what more is being asked here. Jul 19, 2016 at 13:19
• @DavidZ No, in the other question I calculated $\int \delta(r^2-R^2) r^2 dr$ and got two different results. The first was wrong, the second was correct. Now, I am calculating two different integrals which to me should hold the same result. I think I am calculating both of them correctly, and I would like to know why they hold different results because I would expect them to hold the same result. It is very similar but different. Jul 19, 2016 at 13:24
• you have to define first $\delta(r^2 - R^2)$ (and $\delta(\|x\|- R)$ and $\delta(\|x\|^2- R^2)$). change of variables can work with $\delta$ if you are careful and think to $\delta(x)$ in $\mathbb{R}^3$ as the limit of $f_a(x) = \displaystyle\frac{1_{ \|x\| < a}}{4 \pi a^3/3}$ as $a \to 0^+$ Jul 19, 2016 at 13:35

Your result is actually correct even if it may not seem intuitive at first glance. I had a very similar problem here: Compute area of a sphere through a Dirac delta The key part is that: $$\int_{R^n}{f(x)\delta(g(x))|\nabla g(x)|\,dx} = \int_{R^n}{f(x)\delta_S(x)\,dx} = \int_{S}{f(x)\,d\sigma(x)}$$ That means that your parametrisation function $$g(x)$$ is actually very important. If you want to get $$4\pi\rho^2$$ you will have to multiply your delta function with $$|\nabla g(r)| = 2r$$, $$\delta_S(r) = 2r\delta(r^2-\rho^2)$$. Then those would cancel out, and you will get: $$\int_{0}^{2\pi}{\int_{0}^{\pi}{\int_{0}^{\infty}{\delta(r^2-\rho^2)2r^3\sin\theta\,dr}\,d\theta}\,d\phi} = \int_{0}^{2\pi}{\int_{0}^{\pi}{\frac{2\rho^3\sin\theta}{2\rho}\,d\theta}\,d\phi} = 4\pi\rho^2$$
The key thing to understand is that you cannot simply substitute if your delta is of the form $$\delta(g(r))$$ where $$g(r)$$ is not the identity plus/minus a shift. Here's a nice and intuitive explanation of the problem: https://www.mathpages.com/home/kmath663/kmath663.htm
• It looks like my question is basically a duplicate of yours! I have to admit that it is not completely clear to me why the first relation should hold, even if I have (maybe too quickly) read the extensive explanations on the linked answer. However, since I am a physicist, maybe I can bring the point of view of a physicist: it's clear that there is something wrong with my attempt, because the delta function has the physical dimensions of the inverse of its argument. Therefore, it's clear that integrating $\delta(r-R)$ and $\delta(r^2-R^2)$ cannot physically give you the same result. Jan 25, 2019 at 9:54
• @valerio The first relation is a result from geometric measure theory known as the coarea formula (just applied to the Dirac delta in this case, you can check the linked threads and references in my original post to see derivations for the Dirac delta) which lets you integrate over the level sets of a function (exactly what you want since $g^{-1}(0)$ is a level set). I didn't know that the Dirac delta has the physical dimensions of its argument (if you have references do link them since I am very interested in reading more) but I guess that makes sense and explains why you get $2\pi\rho$. Jan 25, 2019 at 11:01
• It has dimensions of the inverse of its argument. It has to be, for it to make physical sense. Here is a Physics SE post about it: physics.stackexchange.com/questions/33760/…. I think the second answer is especially useful to understand why this must be the case (basically it follows from the property $\delta(\alpha x) = \delta(x)/|\alpha|$ if you think of $\alpha$ as a physical unit, for example "meters"). Jan 25, 2019 at 12:58
• @valerio I always thought of it from a mathematical point of view, so the physical perspective definitely makes sense. Note that the identity you gave is a special case to the one I presented (the $|a|$ in the denom is due to the derivative). I found a very good and intuitive explanation concerning your problem specifically: mathpages.com/home/kmath663/kmath663.htm give it a read. Jan 25, 2019 at 17:21