TransWikia.com

Describing a circular current loop as delta functions

Physics Asked on January 10, 2021

It would be really nice to see how Jackson got eqn 5.33 on his example problem for finding the vector potential of a circular current loop

$$
J_{phi}=Isintheta’delta(costheta’)frac{delta(r’-a)}{a}
$$

for describing this geometry:

enter image description here

I’ve found a few ‘explanations’ googling, but I’m still confused on the origin of the $sintheta$

My best guess

$$
int vec{J}cdothat{n} ; d^3x=int vec{I}cdot dvec{l}
nabla_phi left ( int vec{J}cdothat{n} ; d^3x=int vec{I}cdot dvec{l} right)
frac{1}{sintheta} frac{partial}{partialphi} int vec{J}cdothat{n} ; r^2 dr dcostheta dphi = frac{1}{sintheta} frac{partial}{partialphi} int I a dphi
frac{1}{sintheta} int J_phi r^2 dr dcostheta = frac{1}{sintheta} I a
frac{1}{sintheta} int A delta(costheta) delta(r-a) r^2 dr dcostheta =frac{1}{sintheta} I a
Aa^2=Ia
A=frac{I}{a}
therefore J_phi = I delta(costheta) frac{delta(r-a)}{a}
$$

I know the $sintheta$ doesn’t really matter.. because

$$
int sin^2theta delta(theta-pi/2)dtheta=1
$$

but I’d still like to know where it comes from.

EDIT
The explanation that “it doesn’t matter, so lets include it” bothers me. There must be some place in the derivation he picks it up. What possible reason would there be for adding it in at the end?

The change of variables [ref] doesn’t work because

$$
delta(costheta)=frac{delta(theta-pi/2)}{sinfrac{pi}{2}}
$$
so you still end up with
$$
int sin^2theta delta(theta-pi/2)dtheta
$$

I also don’t think that it comes out of the Jacobian. Because he later goes on to integrate over the solid angle

$$
int J r^2 dOmega dr = int J r^2 sintheta dtheta dphi
$$
which would, again, give something of the form $sin^2theta$

EDIT
Jackson’s next step (Eqn. 5.35)

$$
A_phi(r,theta)=frac{mu_0I}{4pi a}int r’^2dr’dOmega’ frac{sintheta’cosphi’delta(cos theta’)delta(r’-a)}{left | mathbf{x}-mathbf{x’} right |}
;
text{with } dOmega’ = dcostheta’ dphi = sintheta’ dtheta’ dphi
$$

I would agree that a $sin theta$ belongs in the formula if he only integrated over $dtheta$ not $sintheta dtheta$

5 Answers

See, you are required to find volume current density $J_phi$. Though its name is volume current density, you know it is the current flowing per unit surface area. Now the subscript $phi$ in $J_phi$ denotes it is flowing in the $hat{phi}$ direction. Now in the spherical polar co-ordinate the infinitesimal length elements along the direction $hat{r},hat{theta} rm and hat{phi}$ are dr, $rdtheta$ and $rsintheta dphi$ respectively. So $drtimes rdtheta$ is the area you are interested in.

Now looking at your figure I am interested in rewriting the current I.

$$I=Iint delta(r-a) dr= Iint delta(r-a) dr deltabig(costheta-cos(pi/2)big) d(costheta)$$ $$=Iint delta(r-a) dr deltabig(costheta-0big) d(costheta)$$ $$=Iint delta(r-a) dr deltabig(costhetabig) d(costheta)$$ $$=Iint delta(r-a) dr deltabig(costhetabig) sintheta dtheta$$ $$=Iint delta(r-a)deltabig(costhetabig) sintheta dr dtheta$$ $$=Iint frac{ delta(r-a)delta(costheta) sintheta}{r} {dr r dtheta}$$

So you see your $$J_phi = Ifrac{ delta(r-a)delta(costheta) sintheta}{r}$$

Correct answer by user22180 on January 10, 2021

EDITED ANSWER: The delta distribution $delta(x)$ is not unique. It is invariant under transformations of the form $delta(x) to f(x)delta(x)$ where $f(0) = 1$. This is because it is really a distribution and not a function. It is mathematically improper to talk about $delta(x)$ instead of $int delta(x)dx$. Derivations of the term you're interested in will not be unique either. You can show that

$$ I = Iint delta(r-a)dr = Iint delta(r-a)delta(theta - pi/2)drdtheta = Iintfrac{delta(r-a)}{r}delta(theta - pi/2)rdrdtheta. $$

From this expression, it is apparent that we can write the current as

$$ I = int J_{phi}rdrdtheta implies J_phi = I frac{delta(r-a)}{r}delta(theta - pi/2). $$

This result leaves out the $sin(theta)$ and replaces $a$ with $r$. It really makes no difference because distributions are mathematically represented by $(T, F)$ where $F$ is a set of continuous functions and $T$ is the map from $F to mathbb{R}$. While it is common for physicists to represent distributions as just the mapping (i.e., $T = (T,F)$), this is somewhat false and leads to the non-unique representation of quantities like $delta(x)$.

Answered by Ultima on January 10, 2021

If you note that $$delta(costheta)=frac{delta(theta-pi/2)}{sintheta}$$ Then you can see that the sine terms actually cancel out.

Answered by yankeefan11 on January 10, 2021

the correct answer can be easily obtained when you use the definition of begin{equation} overrightarrow{J}=rho overrightarrow{v} end{equation}

with $rho=qdelta(overrightarrow{x}-overrightarrow{x'})$, you will have:

begin{equation} overrightarrow{J}=q overrightarrow{v} delta(overrightarrow{x}-overrightarrow{x'}) end{equation}

a small moving charge $dq$ yields to a small current density $doverrightarrow{J}$, so

begin{equation} doverrightarrow{J}=dq overrightarrow{v} delta(overrightarrow{x}-overrightarrow{x'}) end{equation}

the charge will move along the distance $doverrightarrow{l}$:

begin{equation} doverrightarrow{J}=dq {{doverrightarrow{l}}over dt} delta(overrightarrow{x}-overrightarrow{x'}) end{equation}

but one defines $I={dqover dt}$

begin{equation} doverrightarrow{J}= I {doverrightarrow{l}}delta(overrightarrow{x}-overrightarrow{x'}) end{equation}

begin{equation} overrightarrow{J}= int I hspace{0.2cm}{doverrightarrow{l}}hspace{0.2cm}delta(overrightarrow{x}-overrightarrow{x'}) end{equation}

now you must write ${doverrightarrow{l}}$ according the problem as well as the Dirac's delta function. For example, if you have a circular current loop with radius a in a x-y plane and with its center and the origin, you can write ${doverrightarrow{l}}=ahspace{0.1cm}dphi hat{phi}$ and using cylindric coordinates for Dirac's delta:

begin{equation} overrightarrow{J}= int I hspace{0.2cm} ahspace{0.1cm}dphi hat{phi}hspace{0.2cm}{delta(r-a)over r} delta(phi) delta(z) end{equation}

begin{equation} overrightarrow{J}= I hspace{0.2cm} hspace{0.1cm} hat{phi}hspace{0.2cm}delta(r-a) delta(z) end{equation}

I hope that can help you!

Answered by user145629 on January 10, 2021

The following is a systematic way of obtaining the desired description for the current in the wire.

In order to represent the current, we must make use of the azimuthal symmetry to reduce the problem to one involving two coordinates $(u,z)$, check the figure below.

enter image description here

Where $u$ represents a point in the x-y plane. This is related to the spherical coordinates by $$u=rsintheta$$ $$z=rcostheta$$ Notice that rotating $u$ about the z-axis would give us the original problem. The current wrt u and z is given by $$vec{J}=idelta(u-R)delta(z)hat{phi},$$ where $i$ is the current flowing in the wire of radius $R$ and $hat{phi}$ is the unit vector corresponding tangential to the ring. Now using the transformation law for delta functions, we obtain $$vec{J}=idfrac{delta(r-R)delta(theta-frac{pi}{2})}{|J|}hat{phi},$$ where $|J|$ is the Jacobian of the transformation. This is given by $$|J|=vert detbegin{pmatrix}dfrac{partial u}{partial r} & dfrac{partial z}{partial r} dfrac{partial u}{partial theta} & dfrac{partial z}{partial theta}end{pmatrix}vert=vert detbegin{pmatrix}sintheta & costheta rcostheta & -rsinthetaend{pmatrix}vert=r.$$

Plugging in, we obtain $$vec{J}=idfrac{delta(r-R)delta(theta-frac{pi}{2})}{r}hat{phi},$$ which is the desired result.

Answered by Shoham Sen on January 10, 2021

Add your own answers!

Ask a Question

Get help from others!

© 2024 TransWikia.com. All rights reserved. Sites we Love: PCI Database, UKBizDB, Menu Kuliner, Sharing RPP