Peeter Joot's (OLD) Blog.

Math, physics, perl, and programming obscurity.

(A POSSIBLY WRONG) Superposition of transverse electromagnetic field solutions.

Posted by peeterjoot on August 3, 2009

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


FIXME: I’M NOT CONFIDENT THAT i GOT THIS RIGHT. An attempt to verify the final result didn’t work. I’ve either messed up along the way (perhaps right at the beginning), or my verification itself was busted. Am posting these working notes for now, and will revisit later after thinking it through again (or trying the verification again).


In ([1]), a Geometric Algebra solution for the transverse components of the electromagnetic field was found. Here we construct a superposition of these transverse fields, keeping the propagation direction fixed, and allowing for continuous variation of the wave number and angular velocity. Evaluation of this superposition integral, first utilizing a contour integral, then utilizing an inverse Fourier transform allows for the determination of the functional form of a general wave packet moving along a fixed direction in space. This wave packet will be seen to have two time dependencies, an advanced time term, and a retarded time term. The phasors will be eliminated from both the propagation and the transverse fields and will provide operators with which the transverse field can be calculated from a general propagation field wave packet.

Superposition of phasor solutions.

When the field is required to have explicit sinusoidal time and propagation direction, the field components in the transverse (to the propagation) direction were found to be

\begin{aligned}F_t &= \frac{1}{{i \left( \pm k \hat{\mathbf{z}} - \sqrt{\mu\epsilon}\frac{\omega}{c}\right) }} \boldsymbol{\nabla}_t F_z \\ &= \frac{i}{k^2 - \mu\epsilon\frac{\omega^2}{c^2}} \left( \pm k \hat{\mathbf{z}} + \sqrt{\mu\epsilon}\frac{\omega}{c}\right) \boldsymbol{\nabla}_t F_z \end{aligned} \quad\quad\quad(1)

Here F_z = F_z(x,y) = \mathbf{E}_z(x,y) + I\mathbf{B}_z(x,y)/\sqrt{\mu\epsilon} is required by construction to commute with \hat{\mathbf{z}}, but is otherwise arbitrary, except perhaps for boundary value constraints not considered here.

Removing the restriction to fixed wave number and angular velocity we can integrate over both to express a more general transverse wave propagation

\begin{aligned}F_t &= \int d\omega e^{-i\omega t} \int dk e^{ikz} \frac{i}{k^2 - \mu\epsilon\frac{\omega^2}{c^2}} \left(k \hat{\mathbf{z}} + \sqrt{\mu\epsilon}\frac{\omega}{c}\right) \boldsymbol{\nabla}_t {F_z}^{k,\omega} \end{aligned} \quad\quad\quad(3)

Inventing temporary notation for convenience we write for the frequency dependent weighting function {F_z}^{k,\omega} = F_z(x,y,k,\omega). Also note that explicit \pm k has been dropped after allowing wave number to range over both positive and negative values.

Observe that each of these integrals has the form of a Fourier integral (ignoring constant factors that can be incorporated into the weighting function). Also observe that the integral kernel has two poles on the real k-axis (or real \omega-axis). These can be utilized to evaluate one of the integrals using an upper half plane semicircular contour. FIXME: picture here.

Assuming that {F_z}^{k,\omega} is small enough at \infty (on the large contour) to be neglected, we have for the integral after integrating around the poles at k = \pm \sqrt{\mu\epsilon}{\left\lvert{\omega/c}\right\rvert}

\begin{aligned} F_t &= \pi i \int d\omega e^{-i\omega t} {\left.e^{ikz} \frac{i}{k - \sqrt{\mu\epsilon}\frac{{\left\lvert{\omega}\right\rvert}}{c}} \left(k \hat{\mathbf{z}} + \sqrt{\mu\epsilon}\frac{\omega}{c}\right) \boldsymbol{\nabla}_t {F_z}^{k,\omega}\right\vert}_{k=-\sqrt{\mu\epsilon}\frac{{\left\lvert{\omega}\right\rvert}}{c}} \\ &+ \pi i \int d\omega e^{-i\omega t} {\left.e^{ikz} \frac{i}{k + \sqrt{\mu\epsilon}\frac{{\left\lvert{\omega}\right\rvert}}{c}} \left(k \hat{\mathbf{z}} + \sqrt{\mu\epsilon}\frac{\omega}{c}\right) \boldsymbol{\nabla}_t {F_z}^{k,\omega}\right\vert}_{k=\sqrt{\mu\epsilon}\frac{{\left\lvert{\omega}\right\rvert}}{c}} \end{aligned}

Writing {F_z}^{\omega+} = F_z(x,y,\omega, k=\sqrt{\mu\epsilon}{\left\lvert{\omega}\right\rvert}/c), and {F_z}^{\omega-} = F_z(x,y,\omega, k=-\sqrt{\mu\epsilon}{\left\lvert{\omega}\right\rvert}/c), for the values of our field at the poles, we have

\begin{aligned} F_t &= \frac{\pi}{2}\int d\omega e^{-i\omega t} e^{-i \sqrt{\mu\epsilon}\frac{{\left\lvert{\omega}\right\rvert}}{c} z}\frac{\omega -\hat{\mathbf{z}}{\left\lvert{\omega}\right\rvert}}{{\left\lvert{\omega}\right\rvert}}\boldsymbol{\nabla}_t {F_z}^{\omega-} \\ &- \frac{\pi}{2}\int d\omega e^{-i\omega t} e^{i \sqrt{\mu\epsilon}\frac{{\left\lvert{\omega}\right\rvert}}{c} z} \frac{\omega +\hat{\mathbf{z}}{\left\lvert{\omega}\right\rvert}}{{\left\lvert{\omega}\right\rvert}}\boldsymbol{\nabla}_t {F_z}^{\omega+} \end{aligned}

Can one ignore the singularity at \omega = 0 since it divides out? If so, then we have

\begin{aligned} F_t &= \frac{\pi}{2}\int d\omega e^{-i\omega t} e^{-i \sqrt{\mu\epsilon}\frac{\left\lvert{\omega}\right\rvert}{c} z}(\text{sgn}(\omega) - \hat{\mathbf{z}})\boldsymbol{\nabla}_t {F_z}^{\omega-} \\ &- \frac{\pi}{2}\int d\omega e^{-i\omega t} e^{i \sqrt{\mu\epsilon}\frac{\left\lvert{\omega}\right\rvert}{c} z} (\text{sgn}(\omega) + \hat{\mathbf{z}})\boldsymbol{\nabla}_t {F_z}^{\omega+} \end{aligned}

Writing this out separately for the regions greater and lesser than zero we have

\begin{aligned} F_t &= \frac{\pi}{2}\int_{0+}^\infty d\omega e^{-i\omega t} e^{-i \sqrt{\mu\epsilon}\frac{\omega}{c} z}(1 - \hat{\mathbf{z}})\boldsymbol{\nabla}_t {F_z}^{\omega-} - \frac{\pi}{2}\int_{0+}^\infty d\omega e^{-i\omega t} e^{i \sqrt{\mu\epsilon}\frac{\omega}{c} z} (1 + \hat{\mathbf{z}})\boldsymbol{\nabla}_t {F_z}^{\omega+} \\ &- \frac{\pi}{2}\int_{-\infty}^{0-} d\omega e^{-i\omega t} e^{i \sqrt{\mu\epsilon}\frac{\omega}{c} z}(1 + \hat{\mathbf{z}})\boldsymbol{\nabla}_t {F_z}^{\omega-} - \frac{\pi}{2}\int_{-\infty}^{0-} d\omega e^{-i\omega t} e^{-i \sqrt{\mu\epsilon}\frac{\omega}{c} z} (-1 + \hat{\mathbf{z}})\boldsymbol{\nabla}_t {F_z}^{\omega+} \\  \end{aligned}

Grouping by exponentials of like sign, and integrating over a region that omits some neighborhood around the origin, we have eliminated the {\left\lvert{\omega}\right\rvert} factors

\begin{aligned}F_t &= (1 - \hat{\mathbf{z}}) \boldsymbol{\nabla}_t\int d\omega e^{-i\omega \left(t + \sqrt{\mu\epsilon}\frac{z}{c}\right)}\frac{\pi}{2}({F_z}^{\omega-}\theta(\omega) + {F_z}^{\omega+}\theta(-\omega)) \\ &-(1 + \hat{\mathbf{z}}) \boldsymbol{\nabla}_t\int d\omega e^{-i\omega \left(t - \sqrt{\mu\epsilon}\frac{z}{c}\right)}\frac{\pi}{2}({F_z}^{\omega-}\theta(\omega) + {F_z}^{\omega+}\theta(-\omega)) \end{aligned}

The Heaviside unit step function has been used to group things nicely over the doubled integration ranges, and what is left now can be observed to be Fourier transforms from the frequency to time domain. The transformed functions are to be evaluated at a shifted time in each case.

\begin{aligned}F_t &= (1 - \hat{\mathbf{z}}) \boldsymbol{\nabla}_t\frac{\pi}{2}{\left.\mathcal{F}\left( ({F_z}^{\omega-}\theta(\omega) + {F_z}^{\omega+}\theta(-\omega)) \right)\right\vert}_{t=t + \sqrt{\mu\epsilon}\frac{z}{c}} \\ &-(1 + \hat{\mathbf{z}}) \boldsymbol{\nabla}_t\frac{\pi}{2}{\left.\mathcal{F}\left( ({F_z}^{\omega-}\theta(\omega) + {F_z}^{\omega+}\theta(-\omega)) \right)\right\vert}_{t=t - \sqrt{\mu\epsilon}\frac{z}{c}} \end{aligned}

Attempting to actually evaluate this Fourier transform shouldn’t actually be necessary since we the original angular velocity and wave number domain function F(x,y,\omega,k) was arbitrary (except for its commutation with \hat{\mathbf{z}}). Instead we just suppose (once again overloading the symbol F_z) that we have a function F_z(x,y,u) that at time u takes the value

\begin{aligned}F_z(x,y,u) = \frac{\pi}{2}\mathcal{F}\left( ({F_z}^{\omega-}\theta(\omega) + {F_z}^{\omega+}\theta(-\omega)) \right) \end{aligned}

The transform will not change the grades of the propagation field F_z so this still commutes with \hat{\mathbf{z}}. We can now write

\begin{aligned}F_t &= (1 - \hat{\mathbf{z}}) \boldsymbol{\nabla}_t F_z(x,y, t + \sqrt{\mu\epsilon}{z}/{c})-(1 + \hat{\mathbf{z}}) \boldsymbol{\nabla}_t F_z(x,y, t - \sqrt{\mu\epsilon}{z}/{c}) \end{aligned} \quad\quad\quad(4)

So, after a bunch of manipulation, we find exactly how the transverse component of the field is related to the propagation direction field, and have eliminated the phasor description supplied our single frequency/wave-number field relations.

What we have left is a propagation field that has the form of an arbitrary unidirectional wave packet, traveling in either direction through the medium with speed \pm c/\sqrt{\mu\epsilon}. If all this math is right, the transverse field for an advancing wave is generated by application onto the propagation field of the operator

\begin{aligned}-(1 + \hat{\mathbf{z}}) \boldsymbol{\nabla}_t \end{aligned} \quad\quad\quad(5)

Similarly, the transverse field for a receding wave is given by application of the operator

\begin{aligned}(1 - \hat{\mathbf{z}}) \boldsymbol{\nabla}_t \end{aligned} \quad\quad\quad(6)


Given the doubt about the \omega =0 point in the integral, and the possibility for sign error and other algebraic mistakes it seems worthwhile now to go back to Maxwell’s equation and verify the correctness of these results. Then, presuming everything worked out okay, it would also be good to relate things back to the electric and magnetic field components of the field.

Maxwell’s equation in these units is

\begin{aligned}\left(\boldsymbol{\nabla}_t \pm i k \hat{\mathbf{z}} - \sqrt{\mu\epsilon}\frac{i\omega}{c}\right) F(x,y) = 0 \end{aligned} \quad\quad\quad(7)

and we want to compute the transverse direction projection from the propagation direction term

\begin{aligned}F_t &= \mathbf{E}_t + I \mathbf{B}_t = \frac{1}{{2}} (F - \hat{\mathbf{z}} F \hat{\mathbf{z}}) \\ F_z &= \mathbf{E}_z + I \mathbf{B}_z = \frac{1}{{2}} (F + \hat{\mathbf{z}} F \hat{\mathbf{z}})  \end{aligned} \quad\quad\quad(8)

Picking one of the terms of (4), we have

\begin{aligned}F_t = (1 + \pm (1 \mp \hat{\mathbf{z}}) \nabla_t) F_z(t \pm \sqrt{\mu\epsilon} z/c)  \end{aligned}

Substiting into the left hand operator of Maxwell’s equation (7) and reducing I get

\begin{aligned}((\hat{\mathbf{z}} \pm 1) {\nabla_t}^2 + \nabla_t) F_z + \frac{\mu\epsilon}{c}(1 \pm \hat{\mathbf{z}}) F_z' \end{aligned}

I don’t see how this would neccessarily equal zero?


[1] Peeter Joot. Transverse electric and magnetic fields [online].


2 Responses to “(A POSSIBLY WRONG) Superposition of transverse electromagnetic field solutions.”

  1. Eli said

    I haven’t carefully read through this yet, I think there is a physics problem with including ω=0. That corresponds to an electrostatic case, where the electric and magnetic fields are completely independent. Then, with application of boundary conditions (and in the absence of free currents and charges), both fields must go to zero inside the conducting waveguide.

    I do not really understand the notation of GA, so I’m not sure if this result is captured in your expression.

  2. peeterjoot said

    I’ll think about that. I started off with a couple of scribbled results on paper, and as I wrote them up I found mistakes. In the end I think that I corrected most of those mistakes (it was harder and more time consuming doing it without errors;) … In the end I believe that I’ve followed the math from a starting place that is no good.

    I’ve also had a peek at the phasor treatment of Doran/Lasenby’s text (which I now understand a bit better). I see that they (without any calculus) end up with

    F = (1 + \hat{\mathbf{z}})(\alpha \hat{\mathbf{x}} + \beta \hat{\mathbf{y}}) e^{i k \cdot x}

    This has a four vector dot product for their phasor term, with k having 0,3 coordinates, for angular velocity and wave number respectively. It’s similar (first factor is the same, and the second only has the directions in the transverse gradient like mine), but the (1 \pm \hat{\mathbf{z}}) terms show up even before trying to add different frequencies and in the expression for the total field not just the transverse bits.

Leave a Reply

Fill in your details below or click an icon to log in: Logo

You are commenting using your account. Log Out /  Change )

Google+ photo

You are commenting using your Google+ account. Log Out /  Change )

Twitter picture

You are commenting using your Twitter account. Log Out /  Change )

Facebook photo

You are commenting using your Facebook account. Log Out /  Change )


Connecting to %s

%d bloggers like this: