Conservation of charge with Dirac delta

Frostman
Messages
114
Reaction score
17
Homework Statement
Demonstrate conservation of charge
Relevant Equations
Dirac's property
Hello, I was reviewing a part related to electromagnetism in which the charge and current densities are defined by the Dirac delta:

##\rho(\underline{x}, t)=\sum_n e_n \delta^3(\underline{x} - \underline{x}_n(t))##
##\underline{J}(\underline{x}, t)=\sum_n e_n \delta^3(\underline{x} - \underline{x}_n(t))\frac{d\underline{x}_n}{dt}##

At this point, when I want to evaluate the current density divergence, I find a ##3## coefficient that shouldn't appear when I apply a Dirac delta property. I show you here:

##\nabla \cdot \underline{J}(\underline{x}, t)=\sum_n e_n \frac{\partial}{\partial x^i}\delta^3(\underline{x} - \underline{x}_n(t))\frac{dx_n^i}{dt}##

I apply the property for which

##\frac{\partial}{\partial x}\delta(x-y(t))\frac{dy}{dt}=-\frac{\partial}{\partial t}\delta(x-y(t))##

In the three-dimensional case, shouldn't I work like that?

##\frac{\partial}{\partial x^i}\delta^3(\underline{x} - \underline{x}_n(t))\frac{dx_n^i}{dt}=##
##\frac{\partial }{\partial x}\delta(x-x_n(t))\delta(y-y_n(t))\delta(z-z_n(t))\frac{dx_n}{dt}+
\delta(x-x_n(t))\frac{\partial }{\partial y}\delta(y-y_n(t))\delta(z-z_n(t))\frac{dy_n}{dt}+
\delta(x-x_n(t))\delta(y-y_n(t))\frac{\partial }{\partial z}\delta(z-z_n(t))\frac{dz_n}{dt}##

Applying the above property on each of the three addends should produce the same result three times, so:

##
-\frac{\partial }{\partial t}\delta(x-x_n(t))\delta(y-y_n(t))\delta(z-z_n(t))-
\delta(x-x_n(t))\frac{\partial }{\partial t}\delta(y-y_n(t))\delta(z-z_n(t))-
\delta(x-x_n(t))\delta(y-y_n(t))\frac{\partial }{\partial t}\delta(z-z_n(t))
##

Which is therefore equal to:

##-3\frac{\partial }{\partial t}\delta(x-x_n(t))\delta(y-y_n(t))\delta(z-z_n(t))=-3\frac{\partial }{\partial t}\delta^3(\underline{x} - \underline{x}_n(t))##

So when putting it into the current density divergence, I would have a ##3## factor:

##\nabla \cdot \underline{J}(\underline{x}, t)=-3\sum_n e_n \frac{\partial}{\partial t}\delta^3(\underline{x} - \underline{x}_n(t))=-3\frac{\partial \rho(\underline{x}, t)}{\partial t}##

Wrong, I would have:

##\nabla \cdot \underline{J}(\underline{x}, t) + 3\frac{\partial \rho(\underline{x}, t)}{\partial t}=0##

I can't find the point where I'm wrong, can you tell me? o_O
 
Physics news on Phys.org
Frostman said:
Homework Statement:: Demonstrate conservation of charge
Relevant Equations:: Dirac's property

Which is therefore equal to:

−3∂∂tδ(x−xn(t))δ(y−yn(t))δ(z−zn(t))=−3∂∂tδ3(x―−x―n(t))

So when putting it into the current density divergence, I would have a 3 factor:
-[\frac{\partial}{\partial t}\delta(x-x_n(t))]\delta(y-y_n(t))\delta(z-z_n(t))+similar\ for \ y\ and \ z
=-\frac{\partial}{\partial t}[\delta(x-x_n(t))\delta(y-y_n(t))\delta(z-z_n(t))]
 
Last edited:
You don't get a factor of 3. Since ##J = \sum_n \text{ something} \dfrac{d \vec{x}_n}{dt}##, then

##\nabla \cdot J = \sum_n \dfrac{\partial \text{ something}}{\partial x} \dfrac{d x_n}{dt}##
## + \sum_n \dfrac{\partial \text{ something}}{\partial y} \dfrac{d y_n}{dt}##
## + \sum_n \dfrac{\partial \text{ something}}{\partial z} \dfrac{d z_n}{dt}##

Those are three different terms, not 3 copies of the same term. There is no factor of 3.
 
So that factor ##3## really isn't there, because there's ##\sum_n## in front of everything, right?
 
Well, yes, there's a sum over ##n##, but that's not really the point.

You're computing ##\nabla \cdot J ## which involves derivatives of ##J## with respect to ##x_n^i##. There is no derivative with respect to ##t##.
 
I think your mistake was in generalizing from the one-dimensional case.

It is true that
##\frac{d}{dt} f(x-y(t)) = - \frac{d}{dx} f(x-y(t)) \frac{dy}{dt}##

But in 3 dimensions, you have:

##\frac{d}{dt} f(\vec{x}-\vec{y}(t)) ##
##\ \ \ = - \frac{d}{dx^1} f(\vec{x}-\vec{y}(t)) \frac{dy^1}{dt}##
##\ \ \ - \frac{d}{dx^2} f(\vec{x}-\vec{y}(t)) \frac{dy^2}{dt}##
## \ \ \ - \frac{d}{dx^3} f(\vec{x}-\vec{y}(t)) \frac{dy^3}{dt}##

That is not three copies of the same thing.
 
  • Like
Likes Frostman and jim mcnamara
Hi, I had an exam and I completely messed up a problem. Especially one part which was necessary for the rest of the problem. Basically, I have a wormhole metric: $$(ds)^2 = -(dt)^2 + (dr)^2 + (r^2 + b^2)( (d\theta)^2 + sin^2 \theta (d\phi)^2 )$$ Where ##b=1## with an orbit only in the equatorial plane. We also know from the question that the orbit must satisfy this relationship: $$\varepsilon = \frac{1}{2} (\frac{dr}{d\tau})^2 + V_{eff}(r)$$ Ultimately, I was tasked to find the initial...
The value of H equals ## 10^{3}## in natural units, According to : https://en.wikipedia.org/wiki/Natural_units, ## t \sim 10^{-21} sec = 10^{21} Hz ##, and since ## \text{GeV} \sim 10^{24} \text{Hz } ##, ## GeV \sim 10^{24} \times 10^{-21} = 10^3 ## in natural units. So is this conversion correct? Also in the above formula, can I convert H to that natural units , since it’s a constant, while keeping k in Hz ?
Back
Top