# A question of Dirac Delta function (1 Viewer)

### Users Who Are Viewing This Thread (Users: 0, Guests: 1)

#### Psi-String

A vector function

$$V(\vec{r}) = \frac{ \hat r}{r^2}$$

If we calculate it's divergence directly:

$$\nabla \cdot \vec{V} = \frac{1}{r^2} \frac{\partial}{\partial r} \left( r^2 \frac{1}{r^2} \right) = 0$$

However, by divergence theorem, the surface integral is $$4\pi$$. This paradox can be solved by Dirac Delta function.

My problem is, it seems like we can't calculate V's divergence directly because V blow up at r=0, and we unwittingly divide the zero when $$\left( r^2 \frac{1}{r^2} \right)$$

But it seem to me that a vector function $$V'(\vec{r}) =\frac{ \hat r}{r}$$ also has the same problem, which also blows up at r=0, and we will also divide the zero when

$$\nabla \cdot \vec{V'} = \frac{1}{r^2} \frac{\partial}{\partial r} \left( r^2 \frac{1}{r} \right)$$

But we will get the accurate answer !? Why both of the function V and V' has the same problem but we only have to use Dirac Delta function for V ??

Last edited:

#### chanvincent

The potential is a scalar field... but the Electric field is a vector field...

The problem is... you cannot define a divergence on a scalar field.. in other word... the
whole thing doesn't make sense...

However, for a scalar field, you could define a gradient... the gradient of electric potential
is the electric vector field... but you still could not remove the singularity...

#### Psi-String

My function V is vector function, not electric potential.

#### arildno

Homework Helper
Gold Member
Dearly Missed
Simple thing:
The volume integral is NOT defined whenever a singularity is within it; however the surface integral about a singularity IS defined, and equals $4\pi$ IRRESPECTIVE of the shape of the surrounding surface.

It is BECAUSE of this surface shape independence that it is admissible to introduce the Dirac Delta formalism in the first place.

This has nothing to do with the divergence theorem, that has as a PRESUPPOSITION that the integrand is defined on all points in the interior.

#### Psi-String

hmm..... I have one more question....

$$\vec{E} = A e^{- \lambda r} (1+ \lambda r) \frac{\hat{r}}{r^2}$$

where E is electric field. If I want to find charge density, then

$$\rho = \epsilon _0 ( \nabla \cdot \vec{E})$$

In one way

$$\rho = \epsilon _0 ( \nabla \cdot \vec{E}) = \epsilon _0 A e^{-\lambda r} (1+ \lambda r) ( \nabla \cdot \frac{\hat{r}}{r^2}) + \epsilon _0 \frac{\hat{r}}{r^2} \left[ \nabla(A e^{-\lambda r} (1+ \lambda r)) \right] = \epsilon_0 A \left[ 4 \pi \delta ^3 (\vec{r})- \lambda ^2 e^{- \lambda r} /r \right]$$

In another way

$$\rho = \epsilon _0 ( \nabla \cdot \vec{E}) = \epsilon _0 ( \nabla \cdot \vec{E} ) = A \frac{1}{r^2} \frac{ \partial}{ \partial r} ( r^2 E_r) = A \frac{1}{r^2} \frac{ \partial}{ \partial r} ( r^2 \frac{e^{- \lambda r} (1+ \lambda r)}{r^2}) = - \epsilon_0 A \lambda ^2 e^{- \lambda r} /r$$

What is the flaw in the second way, how to avoid it??

Last edited:

#### arildno

Homework Helper
Gold Member
Dearly Missed
There is no flaw in the second method.
Again, what is mathematical "unsound" is to perform a volume integral over a point where the function is not defined.

What is so damn hard about understanding that trivial fact?

#### Meir Achuz

Homework Helper
Gold Member
There is no flaw in the second method.
Again, what is mathematical "unsound" is to perform a volume integral over a point where the function is not defined.

What is so damn hard about understanding that trivial fact?
What you say was the gospel 100 years ago. But work by Schwartz, Lighthill, and others has shown how to do just that using distribution theory.
In physics applications, the basic equation turns out to be
$$\nabla\cdot\left[\frac{\bf r}{r^3}\right]=4\pi\delta({\bf r})$$.
You can derive all the other singular quantities from that.

Last edited:

#### arildno

Homework Helper
Gold Member
Dearly Missed
Do you know what "..." means?

#### StatMechGuy

There is no flaw in the second method.
Again, what is mathematical "unsound" is to perform a volume integral over a point where the function is not defined.

What is so damn hard about understanding that trivial fact?
I don't think he actually did a volume integral at any point of that problem. He took a derivative that is now defined pretty much to leave you with the Dirac delta function.

The problem you speak of is in Jackson's E&M book, and when I worked through it I looked at a neighborhood of radius $$\epsilon$$ around the origin, took the divergence outside it, and then inside it I Taylor expanded for small r, took a divergence of the Taylor expansion, identified the delta function hiding inside, then took the $$\epsilon \rightarrow 0$$ limit and just added things up. That seemed to satisfy the TA and professor enough.

#### Psi-String

At the first place, my problem is that how should I know "what happened" at the origin?? In post 5, my second method didn't tell me any information about the origin. I think the problem is when I taking the dirivative, since the function is undifferentiable at the origin (actually it isn't defined at the origin), so the divergence I get is a function which domain is real number without origin. So I think the best way to avoid this kind of problem is to make sure that the function which I take derivative is differentiable at every point.

#### arildno

Homework Helper
Gold Member
Dearly Missed
If you take the ORDINARY derivative of your function, then indeed it doesn't make sense to take the derivative AT r=0!

The issue here, however, as Meir Achuz has pointed out (and I alluded to it myself in my first post), is that it is perfectly possible to define the Dirac delta (generalized) function by distribution theory in a mathematically rigorous manner.
However, in this rigorous manner, we have to RE-DEFINE "derivative" to what is called "weak derivative". These subtleties are generally avoided in introduction courses in, say, electro-magnetism.

Thus, in your first manner of differentiating, you are differentiating "weakly" (getting the Dirac delta term), whereas in your second method of differentiating you are differentiating ordinarily.

#### Meir Achuz

Homework Helper
Gold Member
At the first place, my problem is that how should I know "what happened" at the origin?? In post 5, my second method didn't tell me any information about the origin. I think the problem is when I taking the dirivative, since the function is undifferentiable at the origin (actually it isn't defined at the origin), so the divergence I get is a function which domain is real number without origin. So I think the best way to avoid this kind of problem is to make sure that the function which I take derivative is differentiable at every point.
You are on the right track. The best definition of the divergence is in terms of the limit of a surface integral divided by the volume as each (V and S) goes to zero. With this definition, div(rvec/r^3)=4pi\delta(r).
The form for div in any specific coordinate system depends on the derivative being defined, and so does not work for r-->0 in this case.

### The Physics Forums Way

We Value Quality
• Topics based on mainstream science
• Proper English grammar and spelling
We Value Civility
• Positive and compassionate attitudes
• Patience while debating
We Value Productivity
• Disciplined to remain on-topic
• Recognition of own weaknesses
• Solo and co-op problem solving