[Click here for a PDF of this post with nicer formatting]

Maxwell’s equations with magnetic sources

The form of Maxwell’s equations to be used here are expressed in terms of \( \boldsymbol{\mathcal{E}} \) and \( \boldsymbol{\mathcal{H}} \), assume linear media, and do not assume a phasor representation

\begin{equation}\label{eqn:energyMomentumWithMagneticSources:120}
\spacegrad \cross \boldsymbol{\mathcal{E}} = – \boldsymbol{\mathcal{M}} – \mu_0 \PD{t}{\boldsymbol{\mathcal{H}}}
\end{equation}
\begin{equation}\label{eqn:energyMomentumWithMagneticSources:140}
\spacegrad \cross \boldsymbol{\mathcal{H}} = \boldsymbol{\mathcal{J}} + \epsilon_0 \PD{t}{\boldsymbol{\mathcal{E}}}
\end{equation}
\begin{equation}\label{eqn:energyMomentumWithMagneticSources:160}
\spacegrad \cdot \boldsymbol{\mathcal{E}} = \rho/\epsilon_0
\end{equation}
\begin{equation}\label{eqn:energyMomentumWithMagneticSources:180}
\spacegrad \cdot \boldsymbol{\mathcal{H}} = \rho_m/\mu_0.
\end{equation}

Energy momentum conservation

With magnetic sources the Poynting and energy conservation relationship has to be adjusted slightly. Let’s derive that result, starting with the divergence of the Poynting vector

\begin{equation}\label{eqn:energyMomentumWithMagneticSources:20}
\begin{aligned}
\spacegrad \cdot \lr{ \boldsymbol{\mathcal{E}} \cross \boldsymbol{\mathcal{H}} }
&=
\boldsymbol{\mathcal{H}} \cdot \lr{ \spacegrad \cross \boldsymbol{\mathcal{E}} }
-\boldsymbol{\mathcal{E}} \cdot \lr{ \spacegrad \cross \boldsymbol{\mathcal{H}} } \\
&=
-\boldsymbol{\mathcal{H}} \cdot \lr{ \mu_0 \partial_t \boldsymbol{\mathcal{H}} + \boldsymbol{\mathcal{M}} }
-\boldsymbol{\mathcal{E}} \cdot \lr{ \boldsymbol{\mathcal{J}} + \epsilon_0 \partial_t \boldsymbol{\mathcal{E}} } \\
&=
– \mu_0 \boldsymbol{\mathcal{H}} \cdot \partial_t \boldsymbol{\mathcal{H}} – \boldsymbol{\mathcal{H}} \cdot \boldsymbol{\mathcal{M}}
– \epsilon_0 \boldsymbol{\mathcal{E}} \cdot \partial_t \boldsymbol{\mathcal{E}} – \boldsymbol{\mathcal{E}} \cdot \boldsymbol{\mathcal{J}},
\end{aligned}
\end{equation}

or

\begin{equation}\label{eqn:energyMomentumWithMagneticSources:40}
\boxed{
\inv{2} \PD{t}{} \lr{ \epsilon_0 \boldsymbol{\mathcal{E}}^2 + \mu_0 \boldsymbol{\mathcal{H}}^2 }
+
\spacegrad \cdot \lr{ \boldsymbol{\mathcal{E}} \cross \boldsymbol{\mathcal{H}} }
=
– \boldsymbol{\mathcal{H}} \cdot \boldsymbol{\mathcal{M}}
– \boldsymbol{\mathcal{E}} \cdot \boldsymbol{\mathcal{J}}.
}
\end{equation}

The usual relationship is only modified by one additional term. Recall from electrodynamics [2] that \ref{eqn:energyMomentumWithMagneticSources:40} (when the magnetic current density \( \boldsymbol{\mathcal{M}} \) is omitted) is just one of four components of the energy momentum conservation equation

\begin{equation}\label{eqn:energyMomentumWithMagneticSources:80}
\partial_\mu T^{\mu \nu} = – \inv{c} F^{\nu \lambda} j_\lambda.
\end{equation}

Note that \ref{eqn:energyMomentumWithMagneticSources:80} was likely not in SI units. The next task is to generalize this classical relationship to incorporate the magnetic sources used in antenna theory. With an eye towards the relativistic nature of the energy momentum tensor, it is natural to assume that the remainder of the energy momentum tensor conservation relation can be found by taking the time derivatives of the Poynting vector.

\begin{equation}\label{eqn:energyMomentumWithMagneticSources:200}
\PD{t}{} \lr{ \boldsymbol{\mathcal{E}} \cross \boldsymbol{\mathcal{H}} }
=
\PD{t}{\boldsymbol{\mathcal{E}}} \cross \boldsymbol{\mathcal{H}}
+ \boldsymbol{\mathcal{E}} \cross \PD{t}{\boldsymbol{\mathcal{H}} }
=
\inv{\epsilon_0}
\lr{ \spacegrad \cross \boldsymbol{\mathcal{H}} – \boldsymbol{\mathcal{J}} } \cross \boldsymbol{\mathcal{H}}
+
\inv{\mu_0}
\boldsymbol{\mathcal{E}} \cross
\lr{
–
\spacegrad \cross \boldsymbol{\mathcal{E}} – \boldsymbol{\mathcal{M}} },
\end{equation}

or

\begin{equation}\label{eqn:energyMomentumWithMagneticSources:220}
\inv{c^2} \PD{t}{} \lr{ \boldsymbol{\mathcal{E}} \cross \boldsymbol{\mathcal{H}} }
+
\mu_0 \boldsymbol{\mathcal{J}} \cross \boldsymbol{\mathcal{H}}
+\epsilon_0
\boldsymbol{\mathcal{E}} \cross \boldsymbol{\mathcal{M}}
=
-\mu_0 \boldsymbol{\mathcal{H}} \cross \lr{ \spacegrad \cross \boldsymbol{\mathcal{H}} }
– \epsilon_0 \boldsymbol{\mathcal{E}} \cross \lr{ \spacegrad \cross \boldsymbol{\mathcal{E}} }.
\end{equation}

The \( \mu_0 \boldsymbol{\mathcal{J}} \cross \boldsymbol{\mathcal{H}} = \boldsymbol{\mathcal{J}} \cross \BB \) is a portion of the Lorentz force equation in its density form. To put \ref{eqn:energyMomentumWithMagneticSources:220} into the desired form, the remainder of the Lorentz force force equation \( \rho \boldsymbol{\mathcal{E}} = \epsilon_0 \boldsymbol{\mathcal{E}} \spacegrad \cdot \boldsymbol{\mathcal{E}} \) must be added to both sides. To extend the magnetic current term to its full dual (magnetic) Lorentz force structure, the quantity to add to both sides is \( \rho_m \boldsymbol{\mathcal{H}} = \mu_0 \boldsymbol{\mathcal{H}} \spacegrad \cdot \boldsymbol{\mathcal{H}} \). Performing these manipulations gives

\begin{equation}\label{eqn:energyMomentumWithMagneticSources:240}
\inv{c^2} \PD{t}{} \lr{ \boldsymbol{\mathcal{E}} \cross \boldsymbol{\mathcal{H}} }
+
\rho \BE + \mu_0 \boldsymbol{\mathcal{J}} \cross \boldsymbol{\mathcal{H}}
+ \rho_m \boldsymbol{\mathcal{H}}
+ \epsilon_0 \boldsymbol{\mathcal{E}} \cross \boldsymbol{\mathcal{M}}
=
\mu_0
\lr{
\boldsymbol{\mathcal{H}} \spacegrad \cdot \boldsymbol{\mathcal{H}}
-\boldsymbol{\mathcal{H}} \cross \lr{ \spacegrad \cross \boldsymbol{\mathcal{H}} }
}
+ \epsilon_0
\lr{
\boldsymbol{\mathcal{E}} \spacegrad \cdot \boldsymbol{\mathcal{E}}
–
\boldsymbol{\mathcal{E}} \cross \lr{ \spacegrad \cross \boldsymbol{\mathcal{E}} }
}.
\end{equation}

It seems slightly surprising the sign of the magnetic equivalent of the Lorentz force terms have an alternation of sign. This is, however, consistent with the duality transformations outlined in ([1] table 3.2)

\begin{equation}\label{eqn:energyMomentumWithMagneticSources:280}
\rho \rightarrow \rho_m
\end{equation}
\begin{equation}\label{eqn:energyMomentumWithMagneticSources:300}
\boldsymbol{\mathcal{J}} \rightarrow \boldsymbol{\mathcal{M}}
\end{equation}
\begin{equation}\label{eqn:energyMomentumWithMagneticSources:320}
\mu_0 \rightarrow \epsilon_0
\end{equation}
\begin{equation}\label{eqn:energyMomentumWithMagneticSources:340}
\boldsymbol{\mathcal{E}} \rightarrow \boldsymbol{\mathcal{H}}
\end{equation}
\begin{equation}\label{eqn:energyMomentumWithMagneticSources:360}
\boldsymbol{\mathcal{H}} \rightarrow -\boldsymbol{\mathcal{E}},
\end{equation}

for

\begin{equation}\label{eqn:energyMomentumWithMagneticSources:380}
\rho \BE + \mu_0 \boldsymbol{\mathcal{J}} \cross \boldsymbol{\mathcal{H}}
\rightarrow
\rho_m \BH + \epsilon_0 \boldsymbol{\mathcal{M}} \cross \lr{ -\boldsymbol{\mathcal{E}}}
=
\rho_m \BH + \epsilon_0 \boldsymbol{\mathcal{E}} \cross \boldsymbol{\mathcal{M}}.
\end{equation}

Comfortable that the LHS has the desired structure, the RHS can expressed as a divergence. Just expanding one of the differences of vector products on the RHS does not obviously show that this is possible, for example

\begin{equation}\label{eqn:energyMomentumWithMagneticSources:400}
\begin{aligned}
\Be_a \cdot
\lr{
\boldsymbol{\mathcal{E}} \spacegrad \cdot \boldsymbol{\mathcal{E}}
–
\boldsymbol{\mathcal{E}} \cross \lr{ \spacegrad \cross \boldsymbol{\mathcal{E}} }
}
&=
E_a \partial_b E_b
–
\epsilon_{a b c} E_b \epsilon_{c r s} \partial_r E_s \\
&=
E_a \partial_b E_b
–
\delta_{a b}^{[r s]} E_b \partial_r E_s \\
&=
E_a \partial_b E_b
–
E_b \lr{
\partial_a E_b
-\partial_b E_a
} \\
&=
E_a \partial_b E_b
– E_b \partial_a E_b
+ E_b \partial_b E_a.
\end{aligned}
\end{equation}

This happens to equal

\begin{equation}\label{eqn:energyMomentumWithMagneticSources:420}
\begin{aligned}
\spacegrad \cdot \lr{ \lr{E_a E_b – \inv{2} \delta_{a b} \boldsymbol{\mathcal{E}}^2 } \Be_b }
&=
\partial_b
\lr{E_a E_b – \inv{2} \delta_{a b} \boldsymbol{\mathcal{E}}^2 } \\
&=
E_b \partial_b E_a
+ E_a \partial_b E_b
–
\inv{2} \delta_{a b} 2 E_c \partial_b E_c \\
i&=
E_b \partial_b E_a
+ E_a \partial_b E_b
– E_b \partial_a E_b.
\end{aligned}
\end{equation}

This allows a final formulation of the remaining energy momentum conservation equation in its divergence form. Let

\begin{equation}\label{eqn:energyMomentumWithMagneticSources:440}
T^{a b} =
\epsilon_0 \lr{ E_a E_b – \inv{2} \delta_{a b} \boldsymbol{\mathcal{E}}^2 }
+ \mu_0 \lr{ H_a H_b – \inv{2} \delta_{a b} \boldsymbol{\mathcal{H}}^2 },
\end{equation}

so that the remaining energy momentum conservation equation, extended to both electric and magnetic sources, is

\begin{equation}\label{eqn:energyMomentumWithMagneticSources:460}
\boxed{
\inv{c^2} \PD{t}{} \lr{ \boldsymbol{\mathcal{E}} \cross \boldsymbol{\mathcal{H}} }
+
\rho \BE + \mu_0 \boldsymbol{\mathcal{J}} \cross \boldsymbol{\mathcal{H}}
+ \rho_m \boldsymbol{\mathcal{H}}
+ \epsilon_0 \boldsymbol{\mathcal{E}} \cross \boldsymbol{\mathcal{M}}
=
\Be_a \spacegrad \cdot \lr{ T^{a b} \Be_b }.
}
\end{equation}

On the LHS we have the rate of change of momentum density, the electric Lorentz force density terms, the dual (magnetic) Lorentz force density terms, and on the RHS the the momentum flux terms.

References

[1] Constantine A Balanis. Antenna theory: analysis and design. John Wiley \& Sons, 3rd edition, 2005.

[2] Peeter Joot. Relativistic Electrodynamics., chapter {Energy Momentum Tensor.} peeterjoot.com, 2011. URL https://peeterjoot.com/archives/math2011/phy450.pdf. [Online; accessed 18-February-2015].