Magnetic dipole


A physical magnetic dipole is a tiny electric current loop. If the size of the loop is infinitesimal, it is called a perfect magnetic dipole. They are characterized by a magnetic dipole moment. The electric analog is the electric dipole.

Under a magnetic field

Consider a rectangular magnetic dipole of sides aa and bb subject to a uniform magnetic field B\mathbf{B}. Any loop can be built up from an infinite number of arbitrarily small rectangles (a Riemann sum), so this discussion can be extended without loss of generality to any shape of the loop. For simplicity, center the loop at the origin and, using Cartesian coordinates, tilt the loop by an angle θ\theta from the zz axis towards the yy axis. Let B\mathbf{B} be in the zz direction.

100%

The forces on the sloping sides cancel out (they stretch the loop, but they don't rotate it). The forces horizontal to the loop F\mathbf{F} also cancel each other out linearly, but not rotationally. They apply a moment of force

N=aFsinθ x^\mathbf{N}=aF\sin \theta\ \hat{\mathbf{x}}

The magnitude of these forces is

F=IbBF=IbB

and therefore

N=IabBsinθ x^=mBsinθ x^\mathbf{N}=IabB\sin \theta\ \hat{\mathbf{x}}=mB\sin \theta\ \hat{\mathbf{x}}

or

N=m×B\boxed{\mathbf{N}=\mathbf{m}\times \mathbf{B}}

where m=Iabm=Iab is the magnetic dipole moment. This moment of force is exact on any magnetic dipole if the field is uniform. If it is not uniform, it is exact only for perfect dipoles, which have no spatial extension and therefore ignore the fact that it's not uniform.

This torque tends to align the dipole in the direction of the magnetic field and accounts for paramagnetism. It might seem then that paramagnetism is the only form of magnetism, as there's nothing here to account for the existence of a diamagnet. However, the magnetic dipole in matter is given by spinning electrons and these, due to the Pauli exclusion principle, tend to pair up with electrons of opposite spin, canceling each other out magnetically. Thus, this rotation is only observed in atoms with odd number of electrons. Even here, the alignment can be broken by thermal collisions.

If the field is uniform, the force over the current loop is zero:

F=IdI×B=I(dI)×B=0\mathbf{F}=I\oint d\mathbf{I}\times \mathbf{B}=I\left( \oint \mathbf{dI} \right)\times \mathbf{B}=\mathbf{0}

but if the field is non-uniform, for a perfect dipole of magnetic moment m\mathbf{m} we have

F=(mB)\boxed{\mathbf{F}=\nabla(\mathbf{m}\cdot \mathbf{B})}

Radiation

If the current going through the loop is alternating, the dipole moment becomes variable. The magnetic field starts to change, inducing an electric field and emitting of electromagnetic radiation.

For a simple model, let's assume the dipole is a circular loop of radius bb and that the current is alternating in a sinusoidal fashion at angular frequency ω\omega:

I(t)=I0cos(ωt)I(t)=I_{0}\cos(\omega t)

The dipole moment becomes

m(t)=πb2I(t)z^=m0cos(ωt)z^\mathbf{m}(t)=\pi b^{2}I(t)\hat{\mathbf{z}}=m_{0}\cos(\omega t)\hat{\mathbf{z}}

by setting the dipole axis on the zz axis and calling m0πb2I0m_{0}\equiv \pi b^{2}I_{0}.

Potentials

We want to find the retarded potentials. We do so in a very similar manner as Electric dipole > Potentials. The loop has no electric charge, so very simply

V=0\boxed{V=0}

Due to the current we have

A(r,t)=μ04πγI0cos(ω(trc))rdI\mathbf{A}(\mathbf{r},t)=\frac{\mu_{0}}{4\pi}\oint_{\gamma} \frac{I_{0}\cos\left( \omega\left( t- \frac{\mathfrak{r}}{c} \right) \right)}{\mathfrak{r}}d\mathbf{I}'

Diagram Magnetic dipole radiation.svg

For a point r\mathbf{r} above the xx axis, A\mathbf{A} must aim in the yy direction, since the xx component is canceled by equal contributions from both sides of the xx axis, so

A(r,t)=μ0I0b4πy^02πcos[ω(trc)]rcosϕdϕ\mathbf{A}(\mathbf{r},t)=\frac{\mu_{0}I_{0}b}{4\pi}\hat{\mathbf{y}}\int_{0}^{2\pi} \frac{\cos\left[ \omega\left( t- \frac{\mathfrak{r}}{c} \right) \right]}{\mathfrak{r}}\cos \phi'd\phi'

where cosϕ\cos \phi' picks out the yy component of dId\mathbf{I}'. Using the law of cosines

r=r2+b22rbcosψ\mathfrak{r}=\sqrt{ r^{2}+b^{2}-2rb\cos \psi }

Since ψ\psi is the angle between r\mathbf{r} and b\mathbf{b}, we have the Scalar product rbcosψ=rb=rbsinθcosϕrb\cos \psi=\mathbf{r}\cdot \mathbf{b}=rb\sin \theta \cos \phi', since

r=rsinθ x^+rcosθ y^,b=bcosϕ x^+bsinϕ y^\mathbf{r}=r\sin \theta\ \hat{\mathbf{x}}+r\cos \theta\ \hat{\mathbf{y}},\qquad \mathbf{b}=b\cos \phi'\ \hat{\mathbf{x}}+b\sin \phi'\ \hat{\mathbf{y}}

The law of cosines then becomes

r=r2+b22rbsinθcosϕ\mathfrak{r}=\sqrt{ r^{2}+b^{2}-2rb\sin \theta \cos \phi' }

We now make the perfect dipole approximation

br(Perfect dipole approximation)b\ll r\tag{Perfect dipole approximation}

which means that the loop is very small compared to the distances we're studying. We can then state

r=r212brsinθcosϕ+b2r2\mathfrak{r}=r^{2}\sqrt{ 1- 2 \frac{b}{r}\sin \theta \cos \phi'+ \frac{b^{2}}{r^{2}} }

b2/r2b^{2}/r^{2} is negligible to first order. The rest can be rewritten using the binomial expansion to first order:

rr2(1brsinθcosϕ)\mathfrak{r}\simeq r^{2}\left( 1- \frac{b}{r}\sin \theta \cos \phi' \right)

The cosine in the integral can then be rewritten as

cos[ω(trc)]cos[ω(trc)+ωbcsinθcosϕ]=cos[ω(trc)]cos(ωbcsinθcosϕ)sin[ω(trc)]sin(ωbcsinθcosϕ)=\begin{align} \cos\left[ \omega\left( t- \frac{\mathfrak{r}}{c} \right) \right]&\simeq \cos\left[ \omega\left( t- \frac{r}{c} \right)+ \frac{\omega b}{c}\sin \theta \cos \phi' \right] \\ &=\cos\left[ \omega\left( t- \frac{r}{c} \right) \right]\cos\left( \frac{\omega b}{c}\sin \theta \cos \phi' \right) \\ &\qquad-\sin\left[ \omega\left( t- \frac{r}{c} \right) \right]\sin\left( \frac{\omega b}{c}\sin \theta \cos \phi' \right) \\ &=\ldots \end{align}

using cos(α+β)=cosαcosβsinαsinβ\cos(\alpha+\beta)=\cos \alpha \cos \beta-\sin \alpha \sin \beta. We now make the second approximation:

bλcω(Large wavelength approximation)b\ll \lambda\sim \frac{c}{\omega}\tag{Large wavelength approximation}

With this we can truncate the Sine and cosine series to first order:

ωbc1cos(ωbcsinθcosϕ)1,sin(ωbcsinθcosϕ)ωbcsinθcosϕ\frac{\omega b}{c}\ll 1\quad\Rightarrow \quad \cos\left( \frac{\omega b}{c}\sin \theta \cos \phi' \right)\simeq 1,\qquad \sin\left( \frac{\omega b}{c}\sin \theta \cos \phi' \right)\simeq \frac{\omega b}{c}\sin \theta \cos \phi'

and so

=cos[ω(trc)]ωbcsinθcosϕsin[ω(trc)]\ldots=\cos\left[ \omega\left( t- \frac{r}{c} \right) \right]- \frac{\omega b}{c}\sin \theta \cos \phi'\sin\left[ \omega\left( t- \frac{r}{c} \right) \right]

Using this result and

1r±1r(1±brsinθcosϕ)\frac{1}{\mathfrak{r}_{\pm}}\simeq \frac{1}{r}\left( 1\pm \frac{b}{r}\sin \theta \cos \phi' \right)

we can compute the integral:

A(r,t)μ0I0b4πry^02π[cos[ω(trc)]++bsinθcosϕ(1rcos[ω(trc)]ωcsin[ω(trc)])]cosϕdϕ\begin{align} \mathbf{A}(\mathbf{r},t)&\simeq \frac{\mu_{0}I_{0}b}{4\pi r}\hat{\mathbf{y}}\int_{0}^{2\pi}\Big[ \cos\left[ \omega\left( t- \frac{r}{c} \right) \right]+ \\ &\qquad \qquad \quad +b\sin \theta \cos \phi'\left( \frac{1}{r}\cos\left[ \omega \left( t- \frac{r}{c} \right) \right]- \frac{\omega}{c}\sin\left[ \omega\left( t- \frac{r}{c} \right) \right] \right) \Big]\cos \phi'd\phi' \end{align}

The first term integrates to zero, since

cos[ω(trc)]02πcosϕdϕ=0\cos\left[ \omega\left( t- \frac{r}{c} \right) \right]\int_{0}^{2\pi}\cos \phi'd\phi'=0

The second and third terms both involve cos2ϕ\cos ^{2}\phi':

brsinθcos[ω(trc)]02πcos2ϕdϕπ=πbrsinθcos[ω(trc)]\frac{b}{r}\sin \theta \cos\left[ \omega\left( t- \frac{r}{c} \right) \right]\underbrace{ \int_{0}^{2\pi}\cos ^{2}\phi'd\phi' }_{ \pi }=\frac{\pi b}{r}\sin \theta \cos\left[ \omega\left( t- \frac{r}{c} \right) \right]

and

bωcsinθsin[ω(trc)]02πcos2ϕdϕπ=πbωcsinθsin[ω(trc)]\frac{b\omega}{c}\sin \theta \sin\left[ \omega\left( t- \frac{r}{c} \right) \right]\underbrace{ \int_{0}^{2\pi}\cos ^{2}\phi'd\phi' }_{ \pi }=\frac{\pi b\omega}{c}\sin \theta \sin\left[ \omega\left( t- \frac{r}{c} \right) \right]

Putting them all together, while noticing that A\mathbf{A} points in the ϕ^\hat{\boldsymbol{\phi}} direction, gives

A(r,θ,t)=μ0m04πrsinθ[1rcos[ω(trc)]ωcsin[ω(trc)]]ϕ^\mathbf{A}(r,\theta,t)=\frac{\mu_{0}m_{0}}{4\pi r}\sin \theta\left[ \frac{1}{r}\cos\left[ \omega\left( t- \frac{r}{c} \right) \right]- \frac{\omega}{c}\sin\left[ \omega\left( t- \frac{r}{c} \right) \right] \right]\hat{\boldsymbol{\phi}}

In the static limit, when ω=0\omega=0, we get the usual

A(r,t)=μ04πm0sinθr2ϕ^\mathbf{A}(r,t)=\frac{\mu_{0}}{4\pi} \frac{m_{0}\sin \theta}{r^{2}}\hat{\boldsymbol{\phi}}

In the radiation zone

rcωor equivalentlyrλ(Far field approximation)r\gg \frac{c}{\omega} \quad\text{or equivalently}\quad r\gg \lambda\tag{Far field approximation}

the first term is r2\sim r^{-2}, so it is negligible and we get

A(r,θ,t)=μ0m0ω4πcrsinθsin[ω(trc)]ϕ^\boxed{\mathbf{A}(r,\theta,t)=- \frac{\mu_{0}m_{0}\omega}{4\pi cr}\sin \theta \sin\left[ \omega\left( t- \frac{r}{c} \right) \right]\hat{\boldsymbol{\phi}}}

Fields

Since V=0\nabla V=0, the fields are quite simple:

E(r,t)=At=μ0m0ω24πcrsinθcos[ω(trc)]ϕ^\boxed{\mathbf{E}(r,t)=-\frac{ \partial \mathbf{A} }{ \partial t } =\frac{\mu_{0}m_{0}\omega ^{2}}{4\pi cr}\sin \theta \cos\left[ \omega\left( t- \frac{r}{c} \right) \right]\hat{\boldsymbol{\phi}}}

and

B(r,t)=×A=μ0m0ω24πc2rsinθcos[ω(trc)]θ^=1cr^×E\boxed{\mathbf{B}(r,t)=\nabla\times \mathbf{A}=- \frac{\mu_{0}m_{0}\omega ^{2}}{4\pi c^{2}r}\sin \theta \cos\left[ \omega\left( t- \frac{r}{c} \right) \right]\hat{\boldsymbol{\theta}}=- \frac{1}{c}\hat{\mathbf{r}}\times \mathbf{E}}

where the far field approximation is used to calculated B\mathbf{B}. Like the electric dipole case, the fields are in phase, mutually perpendicular and transverse. The ratio of their amplitudes is exactly the speed of light: E0/B0=cE_{0}/B_{0}=c.

Electromagnetic waves

These fields produce electromagnetic waves of angular frequency ω\omega. Since they are created from a point source, they are spherical waves.

We can quantify the energy emitted by the oscillating dipole by starting from the Poynting vector:

S(r,t)=1μ0(E×B)=μ0c[m0ω24πcrsinθcos[ω(trc)]]2r^\mathbf{S}(\mathbf{r},t)=\frac{1}{\mu_{0}}(\mathbf{E}\times \mathbf{B})=\frac{\mu_{0}}{c}\left[ \frac{m_{0}\omega ^{2}}{4\pi cr} \sin\theta\cos\left[ \omega\left( t- \frac{r}{c} \right) \right] \right]^{2}\hat{\mathbf{r}}

The irradiance is the magnitude of the time average of the Poynting vector over a full oscillation:

I=S=μ0m02ω432π2c3sin2θr2I=\lvert \langle \mathbf{S} \rangle \rvert =\frac{\mu_{0}m_{0}^{2}\omega^{4}}{32\pi ^{2}c^{3}} \frac{\sin^{2}\theta}{r^{2}}

Like the electric dipole, the irradiance is dependent on the angle and makes a torus shape, with no emission on the axis and maximum emission perpendicular to the axis. The average radiant power is

P=μ0m02ω412πc3\langle P \rangle =\frac{\mu_{0}m_{0}^{2}\omega^{4}}{12\pi c^{3}}

Compared to the electric dipole however, a magnetic dipole emits far less energy. We can see this by comparing powers:

PmagneticPelectric=(m0p0c)2\frac{P_\text{magnetic}}{P_\text{electric}}=\left( \frac{m_{0}}{p_{0}c} \right)^{2}

Since by definition m0=πb2I0m_{0}=\pi b^{2}I_{0} and p0=q0dp_{0}=q_{0}d this is equal to, when setting d=πbd=\pi b for comparison,

PmagneticPelectric=(ωbc)2\frac{P_\text{magnetic}}{P_\text{electric}}=\left( \frac{\omega b}{c} \right)^{2}

This is precisely the quantity that we assumed was small in the large wavelength approximation, and it is squared. This means that the electric power output is enormously greater than the magnetic one and will pretty much always drown it out. The only cases in which it appears are when the system is specifically designed to only emit magnetic dipole radiation, such as the in the pure magnetic dipole itself.