vector potential

Notes for Balantis chapter 4: linear wire antennas.

February 16, 2015 ece1229 , , , , , , , , , , , , , ,

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

These are notes for the UofT course ECE1229, Advanced Antenna Theory, taught by Prof. Eleftheriades, covering ch. 4 [1] content.

Unlike most of the other classes I have taken, I am not attempting to take comprehensive notes for this class. The class is taught on slides that match the textbook so closely, there is little value to me taking notes that just replicate the text. Instead, I am annotating my copy of textbook with little details instead. My usual notes collection for the class will contain musings of details that were unclear, or in some cases, details that were provided in class, but are not in the text (and too long to pencil into my book.)

Magnetic Vector Potential.

In class and in the problem set \( \BA \) was referred to as the Magnetic Vector Potential.  I only recalled this referred to as the Vector Potential.  Prefixing this with magnetic seemed counter intuitive to me since it is generated by electric sources (charges and currents).
This terminology can be justified due to the fact that \( \BA \) generates the magnetic field by its curl. Some mention of this can be found in [4], which also points out that the Electric Potential refers to the scalar \( \phi \). Prof. Eleftheriades points out that Electric Vector Potential refers to the vector potential \( \BF \) generated by magnetic sources (because in that case the electric field is generated by the curl of \( \BF \).)

Plots of infinitesimal dipole radial dependence.

In section 4.2 of [1] are some discussions of the \( kr < 1 \), \( kr = 1 \), and \( kr > 1 \) radial dependence of the fields and power of a solution to an infinitesimal dipole system. Here are some plots of those \( k r \) dependence, along with the \( k r = 1 \) contour as a reference. All the \( \theta \) dependence and any scaling is left out.

The CDF notebook visualizeDipoleFields.cdf is available to interactively plot these, rotate the plots and change the ranges of what is plotted.

A plot of the real and imaginary parts of \( H_\phi = \frac{j k}{r} e^{-j k r} \lr{ 1-\frac{j}{k r} } \) can be found in fig. 1 and fig. 2.

infinitesimalDipoleHphiRealFig3pn

fig 1. Radial dependence of Re H_phi

infinitesimalDipoleHphiImagFig4pn

fig 2. Radial dependence of Im H_phi

 

A plot of the real and imaginary parts of \( E_r = \inv{r^2} \lr{1-\frac{j}{k r}} e^{-j k r} \) can be found in fig. 3 and fig. 4.

infinitesimalDipoleErRealFig1pn

fig 3. Radial dependence of Re E_r

infinitesimalDipoleErImagFig2pn

fig 4. Radial dependence of Im E_r

 

Finally, a plot of the real and imaginary parts of \( E_\theta = \frac{ j k }{r} \lr{1 -\frac{j}{k r} -\frac{1}{k^2 r^2} } e^{-j k r} \) can be found in fig. 5 and fig. 6.

infinitesimalDipoleEthetaRealFig5pn

fig. 5. Radial dependence of Re E_theta

infinitesimalDipoleEthetaImagFig6pn

fig. 6. Radial dependence of Im E_theta

 

Electric Far field for a spherical potential.

It is interesting to look at the far electric field associated with an arbitrary spherical magnetic vector potential, assuming all of the radial dependence is in the spherical envelope. That is

\begin{equation}\label{eqn:chapter4Notes:20}
\BA = \frac{e^{-j k r}}{r} \lr{
\rcap a_r\lr{ \theta, \phi }
+\thetacap a_\theta\lr{ \theta, \phi }
+\phicap a_\phi\lr{ \theta, \phi }
}.
\end{equation}

The electric field is

\begin{equation}\label{eqn:chapter4Notes:40}
\BE = – j \omega \BA – j \frac{1}{\omega \mu_0 \epsilon_0 } \spacegrad \lr{\spacegrad \cdot \BA }.
\end{equation}

The divergence and gradient in spherical coordinates are

\begin{equation}\label{eqn:chapter4Notes:80}
\begin{aligned}
\spacegrad \cdot \BA
&=
\inv{r^2} \PD{r}{} \lr{ r^2 A_r }
+ \inv{r \sin\theta } \PD{\theta}{} \lr{A_\theta \sin\theta}
+ \inv{r \sin\theta } \PD{\phi}{A_\phi}
\end{aligned}
\end{equation}

\begin{equation}\label{eqn:chapter4Notes:100}
\begin{aligned}
\spacegrad \psi \\
&=
\rcap \PD{r}{\psi}
+\frac{\thetacap}{r} \PD{\theta}{\psi}
+ \frac{\phicap}{r \sin\theta} \PD{\phi}{\psi}.
\end{aligned}
\end{equation}

For the assumed potential, the divergence is

\begin{equation}\label{eqn:chapter4Notes:120}
\begin{aligned}
\spacegrad \cdot \BA
&=
\frac{a_r}{r^2} \PD{r}{} \lr{ r^2 \frac{e^{-j k r}}{r} }
+ \inv{r \sin\theta } \frac{e^{-j k r}}{r} \PD{\theta}{} \lr{\sin\theta a_\theta}
+ \inv{r \sin\theta } \frac{e^{-j k r}}{r} \PD{\phi}{a_\phi} \\
&=
a_r
e^{-j k r}
\lr{
\inv{r^2}
-j k \inv{r}
}
+ \inv{r^2 \sin\theta } e^{-j k r} \PD{\theta}{} \lr{\sin\theta a_\theta}
+ \inv{r^2 \sin\theta } e^{-j k r} \PD{\phi}{a_\phi} \\
&\approx
-j k \frac{a_r}{r}
e^{-j k r}.
\end{aligned}
\end{equation}

The last approximation dropped all the \( 1/r^2 \) terms that will be small compared to \( 1/r \) contribution that dominates when \( r \rightarrow \infty \), the far field.

The gradient can now be computed

\begin{equation}\label{eqn:chapter4Notes:140}
\begin{aligned}
\spacegrad \lr{\spacegrad \cdot \BA }
&\approx
-j k
\spacegrad
\lr{
\frac{a_r}{r}
e^{-j k r}
} \\
&=
-j k \lr{
\rcap \PD{r}{}
+\frac{\thetacap}{r} \PD{\theta}{}
+ \frac{\phicap}{r \sin\theta} \PD{\phi}{}
}
\frac{a_r}{r}
e^{-j k r} \\
&=
-j k \lr{
\rcap a_r \PD{r}{} \lr{
\frac{1}{r}
e^{-j k r}
}
+\frac{\thetacap}{r^2}
e^{-j k r}
\PD{\theta}{a_r}
+
e^{-j k r}
\frac{\phicap}{r^2 \sin\theta}
\PD{\phi}{a_r}
} \\
&=
-j k \lr{
-\rcap \frac{a_r}{r^2} \lr{
1
+ j k r
}
+\frac{\thetacap}{r^2}
\PD{\theta}{a_r}
+
\frac{\phicap}{r^2 \sin\theta}
\PD{\phi}{a_r}
}
e^{-j k r} \\
&\approx
– k^2 \rcap \frac{a_r}{r}
e^{-j k r}.
\end{aligned}
\end{equation}

Again, a far field approximation has been used to kill all the \( 1/r^2 \) terms.

The far field approximation of the electric field is now possible

\begin{equation}\label{eqn:chapter4Notes:160}
\begin{aligned}
\BE
&= – j \omega \BA – j \frac{1}{\omega \mu_0 \epsilon_0 } \spacegrad \lr{\spacegrad \cdot \BA } \\
&=
– j \omega
\frac{e^{-j k r}}{r} \lr{
\rcap a_r\lr{ \theta, \phi }
+\thetacap a_\theta\lr{ \theta, \phi }
+\phicap a_\phi\lr{ \theta, \phi }
}
+ j \frac{1}{\omega \mu_0 \epsilon_0 }
k^2 \rcap \frac{a_r}{r}
e^{-j k r} \\
&=
– j \omega
\frac{e^{-j k r}}{r} \lr{
\rcap a_r\lr{ \theta, \phi }
+\thetacap a_\theta\lr{ \theta, \phi }
+\phicap a_\phi\lr{ \theta, \phi }
}
+ j \frac{c^2}{\omega }
\lr{\frac{\omega}{c}}^2 \rcap \frac{a_r}{r}
e^{-j k r}
\\
&=
– j \omega
\frac{e^{-j k r}}{r} \lr{
\thetacap a_\theta\lr{ \theta, \phi }
+\phicap a_\phi\lr{ \theta, \phi }
}.
\end{aligned}
\end{equation}

Observe the perfect, somewhat miraculous seeming, cancellation of all the radial components of the field. If \( \BA_{\textrm{T}} \) is the non-radial projection of \( \BA \), the electric far field is just

\begin{equation}\label{eqn:chapter4Notes:180}
\boxed{
\BE_{\textrm{ff}} = -j \omega \BA_{\textrm{T}}.
}
\end{equation}

Magnetic Far field for a spherical potential.

Application of the same assumed representation for the magnetic field gives
\begin{equation}\label{eqn:chapter4Notes:220}
\begin{aligned}
\BB
&=
\spacegrad \cross \BA \\
&=
\frac{\rcap}{r \sin\theta} \partial_\theta \lr{A_\phi \sin\theta}
+ \frac{\thetacap}{r} \lr{ \inv{\sin\theta} \partial_\phi A_r – \partial_r \lr{r A_\phi}}
+ \frac{\phicap}{r} \lr{ \partial_r\lr{r A_\theta} – \partial_\theta A_r} \\
&=
\frac{\rcap}{r \sin\theta} \partial_\theta \lr{
\frac{e^{-j k r}}{r} a_\phi
\sin\theta}
+ \frac{\thetacap}{r} \lr{ \inv{\sin\theta} \partial_\phi \lr{
\frac{e^{-j k r}}{r} a_r
} – \partial_r \lr{r
\frac{e^{-j k r}}{r} a_\phi
}
}
+ \frac{\phicap}{r} \lr{ \partial_r\lr{r
\frac{e^{-j k r}}{r} a_\theta
} – \partial_\theta
\lr{
\frac{e^{-j k r}}{r} a_r
}
} \\
&=
\frac{\rcap}{r \sin\theta}
\frac{e^{-j k r}}{r}
\partial_\theta \lr{
a_\phi
\sin\theta}
+ \frac{\thetacap}{r} \lr{ \inv{\sin\theta}
\frac{e^{-j k r}}{r}
\partial_\phi
a_r
– \partial_r \lr{
e^{-j k r}
}
a_\phi
}
+ \frac{\phicap}{r} \lr{
\partial_r
\lr{
e^{-j k r}
}
a_\theta

\frac{e^{-j k r}}{r}
\partial_\theta
a_r
}
\approx
j k \lr{ \thetacap a_\phi

\phicap a_\theta
}
\frac{e^{-j k r}}{r} \\
&=
-j k \rcap \cross \lr{
\thetacap a_\theta
+\phicap a_\phi
}
\frac{e^{-j k r}}{r} \\
&=
\inv{c} \BE_{\textrm{ff}}.
\end{aligned}
\end{equation}

The approximation above drops the \( 1/r^2 \) terms. Since

\begin{equation}\label{eqn:chapter4Notes:240}
\inv{\mu_0 c} = \inv{\mu_0} \sqrt{\mu_0\epsilon_0} = \sqrt{\frac{\epsilon_0}{\mu_0}} = \inv{\eta},
\end{equation}

the magnetic far field can be expressed in terms of the electric far field as
\begin{equation}\label{eqn:chapter4Notes:260}
\boxed{
\BH = \inv{\eta} \rcap \cross \BE.
}
\end{equation}

Plane wave relations between electric and magnetic fields

I recalled an identity of the form \ref{eqn:chapter4Notes:260} in [3], but didn’t think that it required a far field approximation.
The reason for this was because the Jackson identity assumed a plane wave representation of the field, something that the far field assumptions also locally require.

Assuming a plane wave representation for both fields

\begin{equation}\label{eqn:chapter4Notes:300}
\boldsymbol{\mathcal{E}}(\Bx, t) = \BE e^{j \lr{\omega t – \Bk \cdot \Bx}}
\end{equation}
\begin{equation}\label{eqn:chapter4Notes:320}
\boldsymbol{\mathcal{B}}(\Bx, t) = \BB e^{j \lr{\omega t – \Bk \cdot \Bx}}
\end{equation}

The cross product relation between the fields follows from the Maxwell-Faraday law of induction

\begin{equation}\label{eqn:chapter4Notes:340}
0 = \spacegrad \cross \boldsymbol{\mathcal{E}} + \PD{t}{\boldsymbol{\mathcal{B}}},
\end{equation}

or

\begin{equation}\label{eqn:chapter4Notes:360}
\begin{aligned}
0
&=
\Be_r \cross \BE \partial_r e^{j\lr{ \omega t – \Bk \cdot \Bx}}
+
j \omega \BB e^{j \lr{\omega t – \Bk \cdot \Bx}} \\
&=
-j \Be_r k_r \cross \BE e^{j \lr{\omega t – \Bk \cdot \Bx}}
+
j \omega \BB e^{j \lr{\omega t – \Bk \cdot \Bx}} \\
&=
\lr{ – \Bk \cross \BE + \omega \BB } j
e^{j \lr{\omega t – \Bk \cdot \Bx}},
\end{aligned}
\end{equation}

or

\begin{equation}\label{eqn:chapter4Notes:380}
\begin{aligned}
\BH
&= \frac{ k}{k c \mu_0 } \kcap \cross \BE \\
&= \inv{ \eta } \kcap \cross \BE,
\end{aligned}
\end{equation}

which also finds \ref{eqn:chapter4Notes:260}, but with much less work and less mess.

Transverse only nature of the far-field fields

Also observe that its possible to tell that the far field fields have only transverse components using the same argument that they are locally plane waves at that distance. The plane waves must satisfy the zero divergence Maxwell’s equations

\begin{equation}\label{eqn:chapter4Notes:420}
\spacegrad \cdot \boldsymbol{\mathcal{E}} = 0
\end{equation}
\begin{equation}\label{eqn:chapter4Notes:440}
\spacegrad \cdot \boldsymbol{\mathcal{B}} = 0,
\end{equation}

so by the same logic

\begin{equation}\label{eqn:chapter4Notes:480}
\Bk \cdot \BE = 0
\end{equation}
\begin{equation}\label{eqn:chapter4Notes:500}
\Bk \cdot \BB = 0.
\end{equation}

In the far field the electric field must equal its transverse projection

\begin{equation}\label{eqn:chapter4Notes:520}
\BE = \textrm{Proj}_\T \lr{-j \omega \BA
– j \frac{1}{\omega \mu_0 \epsilon_0 } \spacegrad \lr{\spacegrad \cdot \BA } }.
\end{equation}

Since by \ref{eqn:chapter4Notes:140} the scalar potential term has only a radial component, that leaves

\begin{equation}\label{eqn:chapter4Notes:540}
\BE = -j \omega \textrm{Proj}_\T \BA,
\end{equation}

which provides \ref{eqn:chapter4Notes:180} with slightly less work.

Vertical dipole reflection coefficient

In class a ground reflection scenario was covered for a horizontal dipole. Reading the text I was surprised to see what looked like the same sort of treatment section 4.7.2, but ending up with a quite different result. It turns out the difference is because the text was treating the vertical dipole configuration, whereas Prof. Eleftheriades was treating a horizontal dipole configuration, which have different reflection coefficients. These differing reflection coefficients are due to differences in the polarization of the field.

To understand these differences in reflection coefficients, consider first the field due to a vertical dipole as sketched in fig. 7, with a wave vector directed from the transmission point downwards in the z-y plane.

verticalDipoleConfigurationFig1

fig. 7. vertical dipole configuration.

 

The wave vector has direction

\begin{equation}\label{eqn:chapter4Notes:560}
\kcap = \zcap e^{\zcap \xcap \theta} = \zcap \cos\theta + \ycap \sin\theta.
\end{equation}

Suppose that the (magnetic) vector potential is that of an infinitesimal dipole

\begin{equation}\label{eqn:chapter4Notes:580}
\BA = \zcap \frac{\mu_0 I_0 l}{4 \pi r} e^{-j k r} %= \frac{A_r}{4 \pi r} e^{-j k r}
\end{equation}

The electric field, in the far field, can be computed by computing the normal projection to the wave vector direction

\begin{equation}\label{eqn:chapter4Notes:600}
\begin{aligned}
\BE
&= -j \omega \lr{\BA \wedge \kcap} \cdot \kcap \\
&= -j \omega \frac{\mu_0 I_0 l}{4 \pi r} \lr{\zcap \wedge \lr{\zcap \cos\theta
+ \ycap \sin\theta} } \lr{\zcap \cos\theta + \ycap \sin\theta} \\
&= -j \omega \frac{\mu_0 I_0 l}{4 \pi r} \lr{ \zcap \ycap \sin\theta }
\lr{\zcap \cos\theta + \ycap \sin\theta} \\
&= -j \omega \frac{\mu_0 I_0 l}{4 \pi r} \sin\theta \lr{-\ycap \cos\theta +
\zcap \sin\theta} \\
&= j \omega \frac{\mu_0 I_0 l}{4 \pi r} \sin\theta \ycap e^{\zcap \ycap \theta}.
\end{aligned}
\end{equation}

This is directed in the z-y plane rotated an additional \( \pi/2 \) past \( \kcap \). The magnetic field must then be directed into the page, along the x axis. This is sketched in fig. 8.

verticalDipoleConfigurationFig2

fig. 8. Electric and magnetic field directions

 

Referring to [2] (\eqntext 4.40) for the coefficient of reflection component

\begin{equation}\label{eqn:chapter4Notes:620}
R
=
\frac{
n_t \cos\theta_i – n_i \cos\theta_t
}
{
n_i \cos\theta_i + n_t \cos\theta_t
}
\end{equation}

This is the Fresnel equation for the case when
that corresponds to

\( \BE \) lies in the plane of incidence, and the magnetic field is completely parallel to the plane of reflection). For the no transmission case, allowing \( v_t \rightarrow 0 \), the index of refraction is \( n_t = c/v_t \rightarrow \infty \), and the reflection coefficient is \( 1 \) as claimed in section 4.7.2 of [1]. Because of the symmetry of this dipole configuration, the azimuthal angle that the wave vector is directed along does not matter.

Horizontal dipole reflection coefficient

In the class notes, a horizontal dipole coming out of the page is indicated. With the page representing the z-y plane, this is a magnetic vector potential directed along the x-axis direction

\begin{equation}\label{eqn:chapter4Notes:640}
\BA = \xcap \frac{\mu_0 I_0 l}{4 \pi r} e^{-j k r}.

\end{equation}

For a wave vector directed in the z-y plane as in \ref{eqn:chapter4Notes:560}, the electric far field is directed along

\begin{equation}\label{eqn:chapter4Notes:660}
\begin{aligned}
\lr{ \xcap \wedge \kcap } \cdot \kcap
&=
\xcap – \lr{ \xcap \cdot \kcap } \kcap \\
&=
\xcap – \lr{ \xcap \cdot \lr{
\zcap \cos\theta + \ycap \sin\theta
} } \kcap \\
&= \xcap.
\end{aligned}
\end{equation}

The electric far field lies completely in the plane of reflection. From [2] (\eqntext 4.34), the Fresnel reflection coefficients is

\begin{equation}\label{eqn:chapter4Notes:680}
R =
\frac{
n_i \cos\theta_i – n_t \cos\theta_t
}
{
n_i \cos\theta_i + n_t \cos\theta_t
},
\end{equation}

which approaches \( -1 \) when \( n_t \rightarrow \infty \). This is consistent with the image theorem summation that Prof. Eleftheriades used in class.

Azimuthal angle dependency of the reflection coefficient

Now consider a horizontal dipole directed along the y-axis. For the same wave vector direction as avove, the electric far field is now directed along

\begin{equation}\label{eqn:chapter4Notes:700}
\begin{aligned}
\lr{ \ycap \wedge \kcap } \cdot \kcap
&=
\ycap – \lr{ \ycap \cdot \kcap } \kcap \\
&=
\ycap – \lr{ \ycap \cdot \lr{
\zcap \cos\theta + \ycap \sin\theta
} } \kcap \\
&=
\ycap – \kcap \sin\theta \\
&=
\ycap – \sin\theta \lr{
\zcap \cos\theta + \ycap \sin\theta
} \\
&=
\ycap \cos^2 \theta – \sin\theta \cos\theta \zcap \\
&= \cos\theta \lr{ \ycap \cos\theta – \sin\theta \zcap } \\
&= \cos\theta \ycap e^{ \zcap \ycap \theta }.
\end{aligned}
\end{equation}

That is

\begin{equation}\label{eqn:chapter4Notes:720}
\BE =
-j \omega \frac{\mu_0 I_0 l}{4 \pi r} e^{-j k r}
\cos\theta \ycap e^{ \zcap \ycap \theta }.
\end{equation}

This far field electric field lies in the plane of incidence (a direction of \( \thetacap \) rotated by \( \pi/2 \)), not in the plane of reflection. The corresponding magnetic field should be directed along the plane of reflection, which is easily confirmed by calculation

\begin{equation}\label{eqn:chapter4Notes:740}
\begin{aligned}
\kcap \cross
\lr{ \ycap \cos\theta – \sin\theta \zcap }
&=
\lr{ \zcap \cos\theta + \ycap \sin\theta } \cross
\lr{ \ycap \cos\theta – \sin\theta \zcap } \\
&=
-\xcap \cos^2 \theta – \xcap \sin^2\theta \\
&= -\xcap.
\end{aligned}
\end{equation}

The far field magnetic field is seen to be

\begin{equation}\label{eqn:chapter4Notes:721}
\BH =
j \omega \frac{I_0 l}{4 \pi r} e^{-j k r}
\cos\theta \xcap,
\end{equation}

so a reflection coefficient of \( 1 \) is required to calculate the power loss after a single ground reflection signal bounce for this relative orientation of antenna to the target.

I fail to see how the horizontal dipole treatment in section 4.7.5 can use a single reflection coefficient without taking into account the azimuthal dependency of that reflection coefficient.

References

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

[2] E. Hecht. Optics. 1998.

[3] JD Jackson. Classical Electrodynamics. John Wiley and Sons, 2nd edition, 1975.

[4] Wikipedia. Magnetic potential — Wikipedia, The Free Encyclopedia, 2015. URL https://en.wikipedia.org/w/index.php?title=Magnetic_potential&oldid=642387563. [Online; accessed 5-February-2015].

Notes for ece1229 antenna theory

February 4, 2015 ece1229 , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , ,

I’ve now posted a first set of notes for the antenna theory course that I am taking this term at UofT.

Unlike most of the other classes I have taken, I am not attempting to take comprehensive notes for this class. The class is taught on slides that match the textbook so closely, there is little value to me taking notes that just replicate the text. Instead, I am annotating my copy of textbook with little details instead. My usual notes collection for the class will contain musings of details that were unclear, or in some cases, details that were provided in class, but are not in the text (and too long to pencil into my book.)

The notes linked above include:

  • Reading notes for chapter 2 (Fundamental Parameters of Antennas) and chapter 3 (Radiation Integrals and Auxiliary Potential Functions) of the class text.
  • Geometric Algebra musings.  How to do formulate Maxwell’s equations when magnetic sources are also included (those modeling magnetic dipoles).
  • Some problems for chapter 2 content.

Recovering the fields

February 4, 2015 ece1229 , , , , , , , , ,

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

This is a small addition to Phasor form of (extended) Maxwell’s equations in Geometric Algebra.

Relative to the observer frame implicitly specified by \( \gamma_0 \), here’s an expansion of the curl of the electric four potential

\begin{equation}\label{eqn:phasorMaxwellsWithElectricAndMagneticCharges:720}
\begin{aligned}
\grad \wedge A_{\textrm{e}}
&=
\inv{2}\lr{
\grad A_{\textrm{e}}

A_{\textrm{e}} \grad
} \\
&=
\inv{2}\lr{
\gamma_0 \lr{ \spacegrad + j k } \gamma_0 \lr{ A_{\textrm{e}}^0 – \BA_{\textrm{e}} }

\gamma_0 \lr{ A_{\textrm{e}}^0 – \BA_{\textrm{e}} } \gamma_0 \lr{ \spacegrad + j k }
} \\
&=
\inv{2}\lr{
\lr{ -\spacegrad + j k } \lr{ A_{\textrm{e}}^0 – \BA_{\textrm{e}} }

\lr{ A_{\textrm{e}}^0 + \BA_{\textrm{e}} } \lr{ \spacegrad + j k }
} \\
&=
\inv{2}\lr{
– 2 \spacegrad A_{\textrm{e}}^0 + j k A_{\textrm{e}}^0 – j k A_{\textrm{e}}^0
+ \spacegrad \BA_{\textrm{e}} – \BA_{\textrm{e}} \spacegrad
– 2 j k \BA_{\textrm{e}}
} \\
&=
– \lr{ \spacegrad A_{\textrm{e}}^0 + j k \BA_{\textrm{e}} }
+ \spacegrad \wedge \BA_{\textrm{e}}
\end{aligned}
\end{equation}

In the above expansion when the gradients appeared on the right of the field components, they are acting from the right (i.e. implicitly using the Hestenes dot convention.)

The electric and magnetic fields can be picked off directly from above, and in the units implied by this choice of four-potential are

\begin{equation}\label{eqn:phasorMaxwellsWithElectricAndMagneticCharges:760}
\BE_{\textrm{e}} = – \lr{ \spacegrad A_{\textrm{e}}^0 + j k \BA_{\textrm{e}} } = -j \lr{ \inv{k}\spacegrad \spacegrad \cdot \BA_{\textrm{e}} + k \BA_{\textrm{e}} }
\end{equation}
\begin{equation}\label{eqn:phasorMaxwellsWithElectricAndMagneticCharges:780}
c \BB_{\textrm{e}} = \spacegrad \cross \BA_{\textrm{e}}.
\end{equation}

For the fields due to the magnetic potentials

\begin{equation}\label{eqn:phasorMaxwellsWithElectricAndMagneticCharges:800}
\lr{ \grad \wedge A_{\textrm{e}} } I
=
– \lr{ \spacegrad A_{\textrm{e}}^0 + j k \BA_{\textrm{e}} } I
– \spacegrad \cross \BA_{\textrm{e}},
\end{equation}

so the fields are

\begin{equation}\label{eqn:phasorMaxwellsWithElectricAndMagneticCharges:840}
c \BB_{\textrm{m}} = – \lr{ \spacegrad A_{\textrm{m}}^0 + j k \BA_{\textrm{m}} } = -j \lr{ \inv{k}\spacegrad \spacegrad \cdot \BA_{\textrm{m}} + k \BA_{\textrm{m}} }
\end{equation}
\begin{equation}\label{eqn:phasorMaxwellsWithElectricAndMagneticCharges:860}
\BE_{\textrm{m}} = -\spacegrad \cross \BA_{\textrm{m}}.
\end{equation}

Including both electric and magnetic sources the fields are

\begin{equation}\label{eqn:phasorMaxwellsWithElectricAndMagneticCharges:900}
\BE = -\spacegrad \cross \BA_{\textrm{m}} -j \lr{ \inv{k}\spacegrad \spacegrad \cdot \BA_{\textrm{e}} + k \BA_{\textrm{e}} }
\end{equation}
\begin{equation}\label{eqn:phasorMaxwellsWithElectricAndMagneticCharges:920}
c \BB = \spacegrad \cross \BA_{\textrm{e}} -j \lr{ \inv{k}\spacegrad \spacegrad \cdot \BA_{\textrm{m}} + k \BA_{\textrm{m}} }
\end{equation}

Dual-Maxwell’s (phasor) equations in Geometric Algebra

February 3, 2015 ece1229 , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , ,

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

These notes repeat (mostly word for word) the previous notes Maxwell’s (phasor) equations in Geometric Algebra. Electric charges and currents have been replaced with magnetic charges and currents, and the appropriate relations modified accordingly.

In [1] section 3.3, treating magnetic charges and currents, and no electric charges and currents, is a demonstration of the required (curl) form for the electric field, and potential form for the electric field. Not knowing what to name this, I’ll call the associated equations the dual-Maxwell’s equations.

I was wondering how this derivation would proceed using the Geometric Algebra (GA) formalism.

Dual-Maxwell’s equation in GA phasor form.

The dual-Maxwell’s equations, omitting electric charges and currents, are

\begin{equation}\label{eqn:phasorDualMaxwellsGA:20}
\spacegrad \cross \boldsymbol{\mathcal{E}} = -\PD{t}{\boldsymbol{\mathcal{B}}} -\BM
\end{equation}
\begin{equation}\label{eqn:phasorDualMaxwellsGA:40}
\spacegrad \cross \boldsymbol{\mathcal{H}} = \PD{t}{\boldsymbol{\mathcal{D}}}
\end{equation}
\begin{equation}\label{eqn:phasorDualMaxwellsGA:60}
\spacegrad \cdot \boldsymbol{\mathcal{D}} = 0
\end{equation}
\begin{equation}\label{eqn:phasorDualMaxwellsGA:80}
\spacegrad \cdot \boldsymbol{\mathcal{B}} = \rho_m.
\end{equation}

Assuming linear media \( \boldsymbol{\mathcal{B}} = \mu_0
\boldsymbol{\mathcal{H}} \), \( \boldsymbol{\mathcal{D}} = \epsilon_0
\boldsymbol{\mathcal{E}} \), and phasor relationships of the form \(
\boldsymbol{\mathcal{E}} = \textrm{Re} \lr{ \BE(\Br) e^{j \omega t}} \) for the fields and the currents, these reduce to

\begin{equation}\label{eqn:phasorDualMaxwellsGA:100}
\spacegrad \cross \BE = – j \omega \BB – \BM
\end{equation}
\begin{equation}\label{eqn:phasorDualMaxwellsGA:120}
\spacegrad \cross \BB = j \omega \epsilon_0 \mu_0 \BE
\end{equation}
\begin{equation}\label{eqn:phasorDualMaxwellsGA:140}
\spacegrad \cdot \BE = 0
\end{equation}
\begin{equation}\label{eqn:phasorDualMaxwellsGA:160}
\spacegrad \cdot \BB = \rho_m.
\end{equation}

These four equations can be assembled into a single equation form using the GA identities

\begin{equation}\label{eqn:phasorDualMaxwellsGA:200}
\Bf \Bg
= \Bf \cdot \Bg + \Bf \wedge \Bg
= \Bf \cdot \Bg + I \Bf \cross \Bg.
\end{equation}
\begin{equation}\label{eqn:phasorDualMaxwellsGA:220}
I = \xcap \ycap \zcap.
\end{equation}

The electric and magnetic field equations, respectively, are

\begin{equation}\label{eqn:phasorDualMaxwellsGA:260}
\spacegrad \BE = – \lr{ \BM + j k c \BB} I
\end{equation}
\begin{equation}\label{eqn:phasorDualMaxwellsGA:280}
\spacegrad c \BB = c \rho_m + j k \BE I
\end{equation}

where \( \omega = k c \), and \( 1 = c^2 \epsilon_0 \mu_0 \) have also been used to eliminate some of the mess of constants.

Summing these (first scaling \ref{eqn:phasorDualMaxwellsGA:280} by \( I \)), gives Maxwell’s equation in its GA phasor form

\begin{equation}\label{eqn:phasorDualMaxwellsGA:300}
\boxed{
\lr{ \spacegrad + j k } \lr{ \BE + I c \BB } = \lr{c \rho – \BM} I.
}
\end{equation}

Preliminaries. Dual magnetic form of Maxwell’s equations.

The arguments of the text showing that a potential representation for the electric and magnetic fields is possible easily translates into GA. To perform this translation, some duality lemmas are required

First consider the cross product of two vectors \( \Bx, \By \) and the right handed dual \( -\By I \) of \( \By \), a bivector, of one of these vectors. Noting that the Euclidean pseudoscalar \( I \) commutes with all grade multivectors in a Euclidean geometric algebra space, the cross product can be written

\begin{equation}\label{eqn:phasorDualMaxwellsGA:320}
\begin{aligned}
\lr{ \Bx \cross \By }
&=
-I \lr{ \Bx \wedge \By } \\
&=
-I \inv{2} \lr{ \Bx \By – \By \Bx } \\
&=
\inv{2} \lr{ \Bx (-\By I) – (-\By I) \Bx } \\
&=
\Bx \cdot \lr{ -\By I }.
\end{aligned}
\end{equation}

The last step makes use of the fact that the wedge product of a vector and vector is antisymmetric, whereas the dot product (vector grade selection) of a vector and bivector is antisymmetric. Details on grade selection operators and how to characterize symmetric and antisymmetric products of vectors with blades as either dot or wedge products can be found in [3], [2].

Similarly, the dual of the dot product can be written as

\begin{equation}\label{eqn:phasorDualMaxwellsGA:440}
\begin{aligned}
-I \lr{ \Bx \cdot \By }
&=
-I \inv{2} \lr{ \Bx \By + \By \Bx } \\
&=
\inv{2} \lr{ \Bx (-\By I) + (-\By I) \Bx } \\
&=
\Bx \wedge \lr{ -\By I }.
\end{aligned}
\end{equation}

These duality transformations are motivated by the observation that in the GA form of Maxwell’s equation the magnetic field shows up in its dual form, a bivector. Spelled out in terms of the dual magnetic field, those equations are

\begin{equation}\label{eqn:phasorDualMaxwellsGA:360}
\spacegrad \cdot (-\BE I)= – j \omega \BB – \BM
\end{equation}
\begin{equation}\label{eqn:phasorDualMaxwellsGA:380}
\spacegrad \wedge \BH = j \omega \epsilon_0 \BE I
\end{equation}
\begin{equation}\label{eqn:phasorDualMaxwellsGA:400}
\spacegrad \wedge (-\BE I) = 0
\end{equation}
\begin{equation}\label{eqn:phasorDualMaxwellsGA:420}
\spacegrad \cdot \BB = \rho_m.
\end{equation}

Constructing a potential representation.

The starting point of the argument in the text was the observation that the triple product \( \spacegrad \cdot \lr{ \spacegrad \cross \Bx } = 0 \) for any (sufficiently continuous) vector \( \Bx \). This triple product is a completely antisymmetric sum, and the equivalent statement in GA is \( \spacegrad \wedge \spacegrad \wedge \Bx = 0 \) for any vector \( \Bx \). This follows from \( \Ba \wedge \Ba = 0 \), true for any vector \( \Ba \), including the gradient operator \( \spacegrad \), provided those gradients are acting on a sufficiently continuous blade.

In the absence of electric charges,
\ref{eqn:phasorDualMaxwellsGA:400} shows that the divergence of the dual electric field is zero. It it therefore possible to find a potential \( \BF \) such that

\begin{equation}\label{eqn:phasorDualMaxwellsGA:460}
-\epsilon_0 \BE I = \spacegrad \wedge \BF.
\end{equation}

Substituting this \ref{eqn:phasorDualMaxwellsGA:380} gives

\begin{equation}\label{eqn:phasorDualMaxwellsGA:480}
\spacegrad \wedge \lr{ \BH + j \omega \BF } = 0.
\end{equation}

This relation is a bivector identity with zero, so will be satisfied if

\begin{equation}\label{eqn:phasorDualMaxwellsGA:500}
\BH + j \omega \BF = -\spacegrad \phi_m,
\end{equation}

for some scalar \( \phi_m \). Unlike the \( -\epsilon_0 \BE I = \spacegrad \wedge \BF \) solution to \ref{eqn:phasorDualMaxwellsGA:400}, the grade of \( \phi_m \) is fixed by the requirement that \( \BE + j \omega \BF \) is unity (a vector), so
a \( \BE + j \omega \BF = \spacegrad \wedge \psi \), for a higher grade blade \( \psi \) would not work, despite satisfying the condition \( \spacegrad \wedge \spacegrad \wedge \psi = 0 \).

Substitution of \ref{eqn:phasorDualMaxwellsGA:500} and \ref{eqn:phasorDualMaxwellsGA:460} into \ref{eqn:phasorDualMaxwellsGA:380} gives

\begin{equation}\label{eqn:phasorDualMaxwellsGA:520}
\begin{aligned}
\spacegrad \cdot \lr{ \spacegrad \wedge \BF } &= -\epsilon_0 \BM – j \omega \epsilon_0 \mu_0 \lr{ -\spacegrad \phi_m -j \omega \BF } \\
\spacegrad^2 \BF – \spacegrad \lr{\spacegrad \cdot \BF} &=
\end{aligned}
\end{equation}

Rearranging gives

\begin{equation}\label{eqn:phasorDualMaxwellsGA:540}
\spacegrad^2 \BF + k^2 \BF = -\epsilon_0 \BM + \spacegrad \lr{ \spacegrad \cdot \BF + j \frac{k}{c} \phi_m }.
\end{equation}

The fields \( \BF \) and \( \phi_m \) are assumed to be phasors, say \( \boldsymbol{\mathcal{A}} = \textrm{Re} \BF e^{j k c t} \) and \( \varphi = \textrm{Re} \phi_m e^{j k c t} \). Grouping the scalar and vector potentials into the standard four vector form
\( F^\mu = \lr{\phi_m/c, \BF} \), and expanding the Lorentz gauge condition

\begin{equation}\label{eqn:phasorDualMaxwellsGA:580}
\begin{aligned}
0
&= \partial_\mu \lr{ F^\mu e^{j k c t}} \\
&= \partial_a \lr{ F^a e^{j k c t}} + \inv{c}\PD{t}{} \lr{ \frac{\phi_m}{c}
e^{j k c t}} \\
&= \spacegrad \cdot \BF e^{j k c t} + \inv{c} j k \phi_m e^{j k c t} \\
&= \lr{ \spacegrad \cdot \BF + j k \phi_m/c } e^{j k c t},
\end{aligned}
\end{equation}

shows that in
\ref{eqn:phasorDualMaxwellsGA:540}
the quantity in braces is in fact the Lorentz gauge condition, so in the Lorentz gauge, the vector potential satisfies a non-homogeneous Helmholtz equation.

\begin{equation}\label{eqn:phasorDualMaxwellsGA:550}
\boxed{
\spacegrad^2 \BF + k^2 \BF = -\epsilon_0 \BM.
}
\end{equation}

Maxwell’s equation in Four vector form

The four vector form of Maxwell’s equation follows from \ref{eqn:phasorDualMaxwellsGA:300} after pre-multiplying by \( \gamma^0 \).

With

\begin{equation}\label{eqn:phasorDualMaxwellsGA:620}
F = F^\mu \gamma_\mu = \lr{ \phi_m/c, \BF }
\end{equation}
\begin{equation}\label{eqn:phasorDualMaxwellsGA:640}
G = \grad \wedge F = – \epsilon_0 \lr{ \BE + c \BB I } I
\end{equation}
\begin{equation}\label{eqn:phasorDualMaxwellsGA:660}
\grad = \gamma^\mu \partial_\mu = \gamma^0 \lr{ \spacegrad + j k }
\end{equation}
\begin{equation}\label{eqn:phasorDualMaxwellsGA:680}
M = M^\mu \gamma_\mu = \lr{ c \rho_m, \BM },
\end{equation}

Maxwell’s equation is

\begin{equation}\label{eqn:phasorDualMaxwellsGA:720}
\boxed{
\grad G = -\epsilon_0 M.
}
\end{equation}

Here \( \setlr{ \gamma_\mu } \) is used as the basis of the four vector Minkowski space, with \( \gamma_0^2 = -\gamma_k^2 = 1 \) (i.e. \(\gamma^\mu \cdot \gamma_\nu = {\delta^\mu}_\nu \)), and \( \gamma_a \gamma_0 = \sigma_a \) where \( \setlr{ \sigma_a} \) is the Pauli basic (i.e. standard basis vectors for \R{3}).

Let’s demonstrate this, one piece at a time. Observe that the action of the spacetime gradient on a phasor, assuming that all time dependence is in the exponential, is

\begin{equation}\label{eqn:phasorDualMaxwellsGA:740}
\begin{aligned}
\gamma^\mu \partial_\mu \lr{ \psi e^{j k c t} }
&=
\lr{ \gamma^a \partial_a + \gamma_0 \partial_{c t} } \lr{ \psi e^{j k c t} }
\\
&=
\gamma_0 \lr{ \gamma_0 \gamma^a \partial_a + j k } \lr{ \psi e^{j k c t} } \\
&=
\gamma_0 \lr{ \sigma_a \partial_a + j k } \psi e^{j k c t} \\
&=
\gamma_0 \lr{ \spacegrad + j k } \psi e^{j k c t}
\end{aligned}
\end{equation}

This allows the operator identification of \ref{eqn:phasorDualMaxwellsGA:660}. The four current portion of the equation comes from

\begin{equation}\label{eqn:phasorDualMaxwellsGA:760}
\begin{aligned}
c \rho_m – \BM
&=
\gamma_0 \lr{ \gamma_0 c \rho_m – \gamma_0 \gamma_a \gamma_0 M^a } \\
&=
\gamma_0 \lr{ \gamma_0 c \rho_m + \gamma_a M^a } \\
&=
\gamma_0 \lr{ \gamma_\mu M^\mu } \\
&= \gamma_0 M.
\end{aligned}
\end{equation}

Taking the curl of the four potential gives

\begin{equation}\label{eqn:phasorDualMaxwellsGA:780}
\begin{aligned}
\grad \wedge F
&=
\lr{ \gamma^a \partial_a + \gamma_0 j k } \wedge \lr{ \gamma_0 \phi_m/c +
\gamma_b F^b } \\
&=
– \sigma_a \partial_a \phi_m/c + \gamma^a \wedge \gamma_b \partial_a F^b – j k
\sigma_b F^b \\
&=
– \sigma_a \partial_a \phi_m/c + \sigma_a \wedge \sigma_b \partial_a F^b – j k
\sigma_b F^b \\
&= \inv{c} \lr{ – \spacegrad \phi_m – j \omega \BF + c \spacegrad \wedge \BF }
\\
&= \epsilon_0 \lr{ c \BB – \BE I } \\
&= – \epsilon_0 \lr{ \BE + c \BB I } I.
\end{aligned}
\end{equation}

Substituting all of these into Maxwell’s \ref{eqn:phasorDualMaxwellsGA:300} gives

\begin{equation}\label{eqn:phasorDualMaxwellsGA:800}
-\frac{\gamma_0}{\epsilon_0}\grad G = \gamma_0 M,
\end{equation}

which recovers \ref{eqn:phasorDualMaxwellsGA:700} as desired.

Helmholtz equation directly from the GA form.

It is easier to find \ref{eqn:phasorDualMaxwellsGA:550} from the GA form of Maxwell’s \ref{eqn:phasorDualMaxwellsGA:700} than the traditional curl and divergence equations. Note that

\begin{equation}\label{eqn:phasorDualMaxwellsGA:820}
\begin{aligned}
\grad G
&=
\grad \lr{ \grad \wedge F } \\
&=
\grad \cdot \lr{ \grad \wedge F } \\
+
\grad \wedge \lr{ \grad \wedge F } \\
&=
\grad^2 F – \grad \lr{ \grad \cdot F },
\end{aligned}
\end{equation}

however, the Lorentz gauge condition \( \partial_\mu F^\mu = \grad \cdot F = 0 \) kills the latter term above. This leaves

\begin{equation}\label{eqn:phasorDualMaxwellsGA:840}
\begin{aligned}
\grad G
&=
\grad^2 F \\
&=
\gamma_0 \lr{ \spacegrad + j k }
\gamma_0 \lr{ \spacegrad + j k } F \\
&=
\gamma_0^2 \lr{ -\spacegrad + j k }
\lr{ \spacegrad + j k } F \\
&=
-\lr{ \spacegrad^2 + k^2 } F = -\epsilon_0 M.
\end{aligned}
\end{equation}

The timelike component of this gives

\begin{equation}\label{eqn:phasorDualMaxwellsGA:860}
\lr{ \spacegrad^2 + k^2 } \phi_m = -\epsilon_0 c \rho_m,
\end{equation}

and the spacelike components give

\begin{equation}\label{eqn:phasorDualMaxwellsGA:880}
\lr{ \spacegrad^2 + k^2 } \BF = -\epsilon_0 \BM,
\end{equation}

recovering \ref{eqn:phasorDualMaxwellsGA:550} as desired.

References

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

[2] C. Doran and A.N. Lasenby. Geometric algebra for physicists. Cambridge University Press New York, Cambridge, UK, 1st edition, 2003.

[3] D. Hestenes. New Foundations for Classical Mechanics. Kluwer Academic Publishers, 1999.

Maxwell’s (phasor) equations in Geometric Algebra

February 1, 2015 ece1229 , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , ,

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

In [1] section 3.2 is a demonstration of the required (curl) form for the magnetic field, and potential form for the electric field.

I was wondering how this derivation would proceed using the Geometric Algebra (GA) formalism.

Maxwell’s equation in GA phasor form.

Maxwell’s equations, omitting magnetic charges and currents, are

\begin{equation}\label{eqn:phasorMaxwellsGA:20}
\spacegrad \cross \boldsymbol{\mathcal{E}} = -\PD{t}{\boldsymbol{\mathcal{B}}}
\end{equation}
\begin{equation}\label{eqn:phasorMaxwellsGA:40}
\spacegrad \cross \boldsymbol{\mathcal{H}} = \boldsymbol{\mathcal{J}} + \PD{t}{\boldsymbol{\mathcal{D}}}
\end{equation}
\begin{equation}\label{eqn:phasorMaxwellsGA:60}
\spacegrad \cdot \boldsymbol{\mathcal{D}} = \rho
\end{equation}
\begin{equation}\label{eqn:phasorMaxwellsGA:80}
\spacegrad \cdot \boldsymbol{\mathcal{B}} = 0.
\end{equation}

Assuming linear media \( \boldsymbol{\mathcal{B}} = \mu_0 \boldsymbol{\mathcal{H}} \), \( \boldsymbol{\mathcal{D}} = \epsilon_0 \boldsymbol{\mathcal{E}} \), and phasor relationships of the form \( \boldsymbol{\mathcal{E}} = \textrm{Re} \lr{ \BE(\Br) e^{j \omega t}} \) for the fields and the currents, these reduce to

\begin{equation}\label{eqn:phasorMaxwellsGA:100}
\spacegrad \cross \BE = – j \omega \BB
\end{equation}
\begin{equation}\label{eqn:phasorMaxwellsGA:120}
\spacegrad \cross \BB = \mu_0 \BJ + j \omega \epsilon_0 \mu_0 \BE
\end{equation}
\begin{equation}\label{eqn:phasorMaxwellsGA:140}
\spacegrad \cdot \BE = \rho/\epsilon_0
\end{equation}
\begin{equation}\label{eqn:phasorMaxwellsGA:160}
\spacegrad \cdot \BB = 0.
\end{equation}

These four equations can be assembled into a single equation form using the GA identities

\begin{equation}\label{eqn:phasorMaxwellsGA:200}
\Bf \Bg
= \Bf \cdot \Bg + \Bf \wedge \Bg
= \Bf \cdot \Bg + I \Bf \cross \Bg.
\end{equation}
\begin{equation}\label{eqn:phasorMaxwellsGA:220}
I = \xcap \ycap \zcap.
\end{equation}

The electric and magnetic field equations, respectively, are

\begin{equation}\label{eqn:phasorMaxwellsGA:260}
\spacegrad \BE = \rho/\epsilon_0 -j k c \BB I
\end{equation}
\begin{equation}\label{eqn:phasorMaxwellsGA:280}
\spacegrad c \BB = \frac{I}{\epsilon_0 c} \BJ + j k \BE I
\end{equation}

where \( \omega = k c \), and \( 1 = c^2 \epsilon_0 \mu_0 \) have also been used to eliminate some of the mess of constants.

Summing these (first scaling \ref{eqn:phasorMaxwellsGA:280} by \( I \)), gives Maxwell’s equation in its GA phasor form

\begin{equation}\label{eqn:phasorMaxwellsGA:300}
\boxed{
\lr{ \spacegrad + j k } \lr{ \BE + I c \BB } = \inv{\epsilon_0 c}\lr{c \rho – \BJ}.
}
\end{equation}

Preliminaries. Dual magnetic form of Maxwell’s equations.

The arguments of the text showing that a potential representation for the electric and magnetic fields is possible easily translates into GA. To perform this translation, some duality lemmas are required

First consider the cross product of two vectors \( \Bx, \By \) and the right handed dual \( -\By I \) of \( \By \), a bivector, of one of these vectors. Noting that the Euclidean pseudoscalar \( I \) commutes with all grade multivectors in a Euclidean geometric algebra space, the cross product can be written

\begin{equation}\label{eqn:phasorMaxwellsGA:320}
\begin{aligned}
\lr{ \Bx \cross \By }
&=
-I \lr{ \Bx \wedge \By } \\
&=
-I \inv{2} \lr{ \Bx \By – \By \Bx } \\
&=
\inv{2} \lr{ \Bx (-\By I) – (-\By I) \Bx } \\
&=
\Bx \cdot \lr{ -\By I }.
\end{aligned}
\end{equation}

The last step makes use of the fact that the wedge product of a vector and vector is antisymmetric, whereas the dot product (vector grade selection) of a vector and bivector is antisymmetric. Details on grade selection operators and how to characterize symmetric and antisymmetric products of vectors with blades as either dot or wedge products can be found in [3], [2].

Similarly, the dual of the dot product can be written as

\begin{equation}\label{eqn:phasorMaxwellsGA:440}
\begin{aligned}
-I \lr{ \Bx \cdot \By }
&=
-I \inv{2} \lr{ \Bx \By + \By \Bx } \\
&=
\inv{2} \lr{ \Bx (-\By I) + (-\By I) \Bx } \\
&=
\Bx \wedge \lr{ -\By I }.
\end{aligned}
\end{equation}

These duality transformations are motivated by the observation that in the GA form of Maxwell’s equation the magnetic field shows up in its dual form, a bivector. Spelled out in terms of the dual magnetic field, those equations are

\begin{equation}\label{eqn:phasorMaxwellsGA:360}
\spacegrad \wedge \BE = – j \omega \BB I
\end{equation}
\begin{equation}\label{eqn:phasorMaxwellsGA:380}
\spacegrad \cdot \lr{ -\BB I } = \mu_0 \BJ + j \omega \epsilon_0 \mu_0 \BE
\end{equation}
\begin{equation}\label{eqn:phasorMaxwellsGA:400}
\spacegrad \cdot \BE = \rho/\epsilon_0
\end{equation}
\begin{equation}\label{eqn:phasorMaxwellsGA:420}
\spacegrad \wedge (-\BB I) = 0.
\end{equation}

Constructing a potential representation.

The starting point of the argument in the text was the observation that the triple product \( \spacegrad \cdot \lr{ \spacegrad \cross \Bx } = 0 \) for any (sufficiently continuous) vector \( \Bx \). This triple product is a completely antisymmetric sum, and the equivalent statement in GA is \( \spacegrad \wedge \spacegrad \wedge \Bx = 0 \) for any vector \( \Bx \). This follows from \( \Ba \wedge \Ba = 0 \), true for any vector \( \Ba \), including the gradient operator \( \spacegrad \), provided those gradients are acting on a sufficiently continuous blade.

In the absence of magnetic charges, \ref{eqn:phasorMaxwellsGA:420} shows that the divergence of the dual magnetic field is zero. It it therefore possible to find a potential \( \BA \) such that

\begin{equation}\label{eqn:phasorMaxwellsGA:460}
\BB I = \spacegrad \wedge \BA.
\end{equation}

Substituting this into Maxwell-Faraday \ref{eqn:phasorMaxwellsGA:360} gives

\begin{equation}\label{eqn:phasorMaxwellsGA:480}
\spacegrad \wedge \lr{ \BE + j \omega \BA } = 0.
\end{equation}

This relation is a bivector identity with zero, so will be satisfied if

\begin{equation}\label{eqn:phasorMaxwellsGA:500}
\BE + j \omega \BA = -\spacegrad \phi,
\end{equation}

for some scalar \( \phi \). Unlike the \( \BB I = \spacegrad \wedge \BA \) solution to \ref{eqn:phasorMaxwellsGA:420}, the grade of \( \phi \) is fixed by the requirement that \( \BE + j \omega \BA \) is unity (a vector), so a \( \BE + j \omega \BA = \spacegrad \wedge \psi \), for a higher grade blade \( \psi \) would not work, despite satisifying the condition \( \spacegrad \wedge \spacegrad \wedge \psi = 0 \).

Substitution of \ref{eqn:phasorMaxwellsGA:500} and \ref{eqn:phasorMaxwellsGA:460} into Ampere’s law \ref{eqn:phasorMaxwellsGA:380} gives

\begin{equation}\label{eqn:phasorMaxwellsGA:520}
\begin{aligned}
-\spacegrad \cdot \lr{ \spacegrad \wedge \BA } &= \mu_0 \BJ + j \omega \epsilon_0 \mu_0 \lr{ -\spacegrad \phi -j \omega \BA } \\
-\spacegrad^2 \BA – \spacegrad \lr{\spacegrad \cdot \BA} &=
\end{aligned}
\end{equation}

Rearranging gives

\begin{equation}\label{eqn:phasorMaxwellsGA:540}
\spacegrad^2 \BA + k^2 \BA = -\mu_0 \BJ – \spacegrad \lr{ \spacegrad \cdot \BA + j \frac{k}{c} \phi }.
\end{equation}

The fields \( \BA \) and \( \phi \) are assumed to be phasors, say \( \boldsymbol{\mathcal{A}} = \textrm{Re} \BA e^{j k c t} \) and \( \varphi = \textrm{Re} \phi e^{j k c t} \). Grouping the scalar and vector potentials into the standard four vector form \( A^\mu = \lr{\phi/c, \BA} \), and expanding the Lorentz gauge condition

\begin{equation}\label{eqn:phasorMaxwellsGA:580}
\begin{aligned}
0
&= \partial_\mu \lr{ A^\mu e^{j k c t}} \\
&= \partial_a \lr{ A^a e^{j k c t}} + \inv{c}\PD{t}{} \lr{ \frac{\phi}{c} e^{j k c t}} \\
&= \spacegrad \cdot \BA e^{j k c t} + \inv{c} j k \phi e^{j k c t} \\
&= \lr{ \spacegrad \cdot \BA + j k \phi/c } e^{j k c t},
\end{aligned}
\end{equation}

shows that in \ref{eqn:phasorMaxwellsGA:540} the quantity in braces is in fact the Lorentz gauge condition, so in the Lorentz gauge, the vector potential satisfies a non-homogeneous Helmholtz equation.

\begin{equation}\label{eqn:phasorMaxwellsGA:550}
\boxed{
\spacegrad^2 \BA + k^2 \BA = -\mu_0 \BJ.
}
\end{equation}

Maxwell’s equation in Four vector form

The four vector form of Maxwell’s equation follows from \ref{eqn:phasorMaxwellsGA:300} after pre-multiplying by \( \gamma^0 \).

With

\begin{equation}\label{eqn:phasorMaxwellsGA:620}
A = A^\mu \gamma_\mu = \lr{ \phi/c, \BA }
\end{equation}
\begin{equation}\label{eqn:phasorMaxwellsGA:640}
F = \grad \wedge A = \inv{c} \lr{ \BE + c \BB I }
\end{equation}
\begin{equation}\label{eqn:phasorMaxwellsGA:660}
\grad = \gamma^\mu \partial_\mu = \gamma^0 \lr{ \spacegrad + j k }
\end{equation}
\begin{equation}\label{eqn:phasorMaxwellsGA:680}
J = J^\mu \gamma_\mu = \lr{ c \rho, \BJ },
\end{equation}

Maxwell’s equation is

\begin{equation}\label{eqn:phasorMaxwellsGA:700}
\boxed{
\grad F = \mu_0 J.
}
\end{equation}

Here \( \setlr{ \gamma_\mu } \) is used as the basis of the four vector Minkowski space, with \( \gamma_0^2 = -\gamma_k^2 = 1 \) (i.e. \(\gamma^\mu \cdot \gamma_\nu = {\delta^\mu}_\nu \)), and \( \gamma_a \gamma_0 = \sigma_a \) where \( \setlr{ \sigma_a} \) is the Pauli basic (i.e. standard basis vectors for \R{3}).

Let’s demonstrate this, one piece at a time. Observe that the action of the spacetime gradient on a phasor, assuming that all time dependence is in the exponential, is

\begin{equation}\label{eqn:phasorMaxwellsGA:740}
\begin{aligned}
\gamma^\mu \partial_\mu \lr{ \psi e^{j k c t} }
&=
\lr{ \gamma^a \partial_a + \gamma_0 \partial_{c t} } \lr{ \psi e^{j k c t} }
\\
&=
\gamma_0 \lr{ \gamma_0 \gamma^a \partial_a + j k } \lr{ \psi e^{j k c t} } \\
&=
\gamma_0 \lr{ \sigma_a \partial_a + j k } \psi e^{j k c t} \\
&=
\gamma_0 \lr{ \spacegrad + j k } \psi e^{j k c t}
\end{aligned}
\end{equation}

This allows the operator identification of \ref{eqn:phasorMaxwellsGA:660}. The four current portion of the equation comes from

\begin{equation}\label{eqn:phasorMaxwellsGA:760}
\begin{aligned}
c \rho – \BJ
&=
\gamma_0 \lr{ \gamma_0 c \rho – \gamma_0 \gamma_a \gamma_0 J^a } \\
&=
\gamma_0 \lr{ \gamma_0 c \rho + \gamma_a J^a } \\
&=
\gamma_0 \lr{ \gamma_\mu J^\mu } \\
&= \gamma_0 J.
\end{aligned}
\end{equation}

Taking the curl of the four potential gives

\begin{equation}\label{eqn:phasorMaxwellsGA:780}
\begin{aligned}
\grad \wedge A
&=
\lr{ \gamma^a \partial_a + \gamma_0 j k } \wedge \lr{ \gamma_0 \phi/c + \gamma_b A^b } \\
&=
– \sigma_a \partial_a \phi/c + \gamma^a \wedge \gamma_b \partial_a A^b – j k
\sigma_b A^b \\
&=
– \sigma_a \partial_a \phi/c + \sigma_a \wedge \sigma_b \partial_a A^b – j k
\sigma_b A^b \\
&= \inv{c} \lr{ – \spacegrad \phi – j \omega \BA + c \spacegrad \wedge \BA }
\\
&= \inv{c} \lr{ \BE + c \BB I }.
\end{aligned}
\end{equation}

Substituting all of these into Maxwell’s \ref{eqn:phasorMaxwellsGA:300} gives

\begin{equation}\label{eqn:phasorMaxwellsGA:800}
\gamma_0 \grad c F = \inv{ \epsilon_0 c } \gamma_0 J,
\end{equation}

which recovers \ref{eqn:phasorMaxwellsGA:700} as desired.

Helmholtz equation directly from the GA form.

It is easier to find \ref{eqn:phasorMaxwellsGA:550} from the GA form of Maxwell’s \ref{eqn:phasorMaxwellsGA:700} than the traditional curl and divergence equations. Note that

\begin{equation}\label{eqn:phasorMaxwellsGA:820}
\grad F
=
\grad \lr{ \grad \wedge A }
=
\grad \cdot \lr{ \grad \wedge A }
+
\grad \wedge \lr{ \grad \wedge A }
=
\grad^2 A – \grad \lr{ \grad \cdot A },
\end{equation}

however, the Lorentz gauge condition \( \partial_\mu A^\mu = \grad \cdot A = 0 \) kills the latter term above. This leaves

\begin{equation}\label{eqn:phasorMaxwellsGA:840}
\begin{aligned}
\grad F
&=
\grad^2 A \\
&=
\gamma_0 \lr{ \spacegrad + j k }
\gamma_0 \lr{ \spacegrad + j k } A \\
&=
\gamma_0^2 \lr{ -\spacegrad + j k }
\lr{ \spacegrad + j k } A \\
&=
-\lr{ \spacegrad^2 + k^2 } A = \mu_0 J.
\end{aligned}
\end{equation}

The timelike component of this gives

\begin{equation}\label{eqn:phasorMaxwellsGA:860}
\lr{ \spacegrad^2 + k^2 } \phi = -\rho/\epsilon_0,
\end{equation}

and the spacelike components give

\begin{equation}\label{eqn:phasorMaxwellsGA:880}
\lr{ \spacegrad^2 + k^2 } \BA = -\mu_0 \BJ,
\end{equation}

recovering \ref{eqn:phasorMaxwellsGA:550} as desired.

References

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

[2] C. Doran and A.N. Lasenby. Geometric algebra for physicists. Cambridge University Press New York, Cambridge, UK, 1st edition, 2003.

[3] D. Hestenes. New Foundations for Classical Mechanics. Kluwer Academic Publishers, 1999.

%d bloggers like this: