Fourier expansion electromagnetic field

The electromagnetic (EM) field is of importance as a carrier of solar energy and electronic signals (radio, TV, etc.). As its name suggests, it consists of two vector fields, the electric field E and the magnetic field B. The Fourier expansion of the electromagnetic field is used in the quantization that leads to photons, light particles of well-defined energy and momentum. Further the Fourier transform of the EM field plays a role in theory of wave propagation through different media and also in light scattering.

In the absence of charges and electric currents (that is, in "free" space), both E and B can be derived from a third real vector field, the vector potential A. Below the Fourier transform of the fields E, B, and A will be discussed. It will be seen that the expansion of the vector potential A yields the expansions of the real fields E and B. Further the energy and momentum of the EM field will be expressed in the Fourier components of A.

Fourier expansion of a vector field

A real scalar function f(x), with 0 ≤ x ≤ L, has the following Fourier expansion:

${\begin{aligned}f(x)&=A_{0}+\sum _{n>0}\left[A_{n}\cos({\tfrac {2n\pi }{L}}x)+B_{n}\sin({\tfrac {2n\pi }{L}}x)\right]=A_{0}+{\frac {1}{2}}\sum _{n>0}\left[(A_{n}-iB_{n})e^{2n\pi ix/L}+(A_{n}+iB_{n})e^{-2n\pi ix/L}\right]\\&=\sum _{n\geq 0}\left[f_{n}e^{2n\pi ix/L}+{\bar {f}}_{n}e^{-2n\pi ix/L}\right],\end{aligned}}$

where the bar indicates complex conjugation and the definition of the Fourier components is obvious. For a real scalar field f(r) depending on three variables, r ≡ (x, y, z), the Fourier expansion is easily generalized; it is the following:

$f(\mathbf {r} ,t)=\sum _{\mathbf {k} }\left(f_{k}(t)e^{i\mathbf {k} \cdot \mathbf {r} }+{\bar {f}}_{k}(t)e^{-i\mathbf {k} \cdot \mathbf {r} }\right),\qquad \mathbf {k} ={\frac {2\pi }{L}}(n_{x},\;n_{y},\;n_{z}),\quad {\hbox{with}}\quad n_{x},\,n_{y},\,n_{z}=0,\;1,\;2,\ldots$

Such an expansion, labeled by a discrete (countable) set of vectors k, is always possible when f satisfies periodic boundary conditions, i.e., f(r + p,t) = f(r,t) for some finite vector p. To impose such boundary conditions, it is common to consider EM waves as if they are in a virtual cubic box of finite volume V = L³. Waves on opposite walls of the box are enforced to have the same value (usually zero). Note that the waves are not restricted to the box: the box is replicated an infinite number of times in x, y, and z direction.

The expansion above is over the first octant of the (k_x, k_y, k_z) lattice. It is often extended to run over all octants, in which case an overcomplete (linearly dependent, not completely orthogonal) basis is used. When furthermore the expansion is applied to the three components of a vector field F(r,t) separately, the result can be written concisely as follows:

$\mathbf {F} (\mathbf {r} ,t)=\sum _{\mathbf {k} }\left(\mathbf {f} _{k}(t)e^{i\mathbf {k} \cdot \mathbf {r} }+{\bar {\mathbf {f} }}_{k}(t)e^{-i\mathbf {k} \cdot \mathbf {r} }\right),\qquad \mathbf {k} ={\frac {2\pi }{L}}(n_{x},\;n_{y},\;n_{z}),\quad {\hbox{with}}\quad n_{x},\,n_{y},\,n_{z}=0,\pm \;1,\;\pm 2,\ldots$

with the following orthogonality relations,

${\begin{aligned}\iiint _{V}e^{-i\mathbf {k} \cdot \mathbf {r} }e^{i\mathbf {k'} \cdot \mathbf {r} }\mathrm {d} ^{3}\mathbf {r} &=V\delta _{\mathbf {k} ,\mathbf {k} '},\quad \iiint _{V}e^{i\mathbf {k} \cdot \mathbf {r} }e^{i\mathbf {k'} \cdot \mathbf {r} }\mathrm {d} ^{3}\mathbf {r} &=V\delta _{\mathbf {k} ,-\mathbf {k} '},\quad \iiint _{V}e^{-i\mathbf {k} \cdot \mathbf {r} }e^{-i\mathbf {k'} \cdot \mathbf {r} }\mathrm {d} ^{3}\mathbf {r} &=V\delta _{\mathbf {k} ,-\mathbf {k} '}.\end{aligned}}$

Vector potential

The magnetic field B satisfies the following Maxwell equation:

${\boldsymbol {\nabla }}\cdot \mathbf {B} (\mathbf {r} ,t)=0,$

that is, the divergence of B is zero. This equation expresses the fact that magnetic monopoles (charges) do not exist (or, rather, have never been found in nature). A divergence-free field, such as B, is a also referred to as a transverse field. By the Helmholtz decomposition, B can be written as

$\mathbf {B} (\mathbf {r} ,t)={\boldsymbol {\nabla }}\times \mathbf {A} (\mathbf {r} ,t),$

in which the vector potential A is introduced though the curl ∇×A.

The electric field obeys one of the Maxwell equations, in electromagnetic SI units it reads,

${\boldsymbol {\nabla }}\cdot \mathbf {E} (\mathbf {r} ,t)={\frac {\rho }{\epsilon _{0}}}=0$

because it is assumed that charge distributions ρ are zero. The quantity ε₀ is the electric constant. Hence, also the electric field E is transverse. Since there are no charges, the electric potential is zero and the electric field follows from A by,

$\mathbf {E} (\mathbf {r} ,t)=-{\frac {\partial \mathbf {A} (\mathbf {r} ,t)}{\partial t}}.$

The fact that E can be written this way is due to the choice of Coulomb gauge for A:

${\boldsymbol {\nabla }}\cdot \mathbf {A} (\mathbf {r} ,t)=0.$

By definition, a choice of gauge does not affect any measurable properties (the best known example of a choice of gauge is the fixing of the zero of an electric potential, for instance at infinity). The Coulomb gauge makes A transverse as well, and clearly A is parallel to E. (The time differentiation does not affect direction.) The vector fields A, B, and E are in the same plane, as they are orthogonal to the same longitudinal vector. Hence, the three fields can be written as a linear combination of two orthonormal vectors, e_x and e_y. It is often convenient to choose complex unit vectors obtained by a unitary transformation,

$\mathbf {e} ^{(1)}\equiv {\frac {-1}{\sqrt {2}}}(\mathbf {e} _{x}+i\mathbf {e} _{y})\quad {\hbox{and}}\quad \mathbf {e} ^{(-1)}\equiv {\frac {1}{\sqrt {2}}}(\mathbf {e} _{x}-i\mathbf {e} _{y})$

which are orthonormal,

$\mathbf {e} ^{(\mu )}\cdot {\bar {\mathbf {e} }}^{(\mu ')}=\delta _{\mu ,\mu '}\quad {\hbox{with}}\quad \mu ,\mu '=1,\,-1,$

and the bar indicating the complex conjugate.

Expansions

The Fourier expansion of the vector potential reads

$\mathbf {A} (\mathbf {r} ,t)=\sum _{\mathbf {k} }\sum _{\mu =-1,1}\left(\mathbf {e} ^{(\mu )}(\mathbf {k} )a_{\mathbf {k} }^{(\mu )}(t)\,e^{i\mathbf {k} \cdot \mathbf {r} }+{\bar {\mathbf {e} }}^{(\mu )}(\mathbf {k} ){\bar {a}}_{\mathbf {k} }^{(\mu )}(t)\,e^{-i\mathbf {k} \cdot \mathbf {r} }\right).$

The vector potential obeys the wave equation,

$\nabla ^{2}\mathbf {A} (\mathbf {r} ,t)={\frac {1}{c^{2}}}{\frac {\partial ^{2}\mathbf {A} (\mathbf {r} ,t)}{\partial t^{2}}}$

The substitution of the Fourier series of A into the wave equation yields for the individual terms,

$-k^{2}a_{\mathbf {k} }^{(\mu )}(t)={\frac {1}{c^{2}}}{\frac {\partial ^{2}a_{\mathbf {k} }^{(\mu )}(t)}{\partial t^{2}}}\quad \Longrightarrow \quad a_{\mathbf {k} }^{(\mu )}(t)\propto e^{-i\omega t}\quad {\hbox{with}}\quad \omega =kc\quad {\hbox{and}}\quad k\equiv |\mathbf {k} |.$

It is now an easy matter to construct the corresponding Fourier expansions for E and B from the expansion of the vector potential A. The expansion for E follows from differentiation with respect to time,

$\mathbf {E} (\mathbf {r} ,t)=i\sum _{\mathbf {k} }\sum _{\mu =-1,1}\omega \left(\mathbf {e} ^{(\mu )}(\mathbf {k} )a_{\mathbf {k} }^{(\mu )}(t)\,e^{i\mathbf {k} \cdot \mathbf {r} }-{\bar {\mathbf {e} }}^{(\mu )}(\mathbf {k} ){\bar {a}}_{\mathbf {k} }^{(\mu )}(t)\,e^{-i\mathbf {k} \cdot \mathbf {r} }\right).$

The expansion for B follows by taking the curl,

$\mathbf {B} (\mathbf {r} ,t)=i\sum _{\mathbf {k} }\sum _{\mu =-1,1}\left({\big [}\mathbf {k} \times \mathbf {e} ^{(\mu )}(\mathbf {k} ){\big ]}\;a_{\mathbf {k} }^{(\mu )}(t)\,e^{i\mathbf {k} \cdot \mathbf {r} }-{\big [}\mathbf {k} \times {\bar {\mathbf {e} }}^{(\mu )}(\mathbf {k} ){\big ]}\;{\bar {a}}_{\mathbf {k} }^{(\mu )}(t)\,e^{-i\mathbf {k} \cdot \mathbf {r} }\right).$

Fourier-expanded energy

The electromagnetic energy density is

${\mathcal {E}}_{\mathrm {Field} }={\frac {1}{2}}(\epsilon _{0}\mathbf {E} \cdot \mathbf {E} +{\frac {1}{\mu _{0}}}\mathbf {B} \cdot \mathbf {B} ),$

where μ₀ is the magnetic constant. The total energy (classical Hamiltonian) of a finite volume V is defined by

$H=\iiint _{V}{\mathcal {E}}_{\mathrm {Field} }(\mathbf {r} ,t)\mathrm {d} ^{3}\mathbf {r} .$

Use in the expansion of the electric field squared the following

${\begin{aligned}{\frac {1}{V}}\iiint _{V}&\left(\mathbf {e} ^{(\mu ')}(\mathbf {k'} )a_{\mathbf {k'} }^{(\mu ')}(t)\,e^{i\mathbf {k'} \cdot \mathbf {r} }-{\bar {\mathbf {e} }}^{(\mu ')}(\mathbf {k'} ){\bar {a}}_{\mathbf {k'} }^{(\mu ')}(t)\,e^{-i\mathbf {k'} \cdot \mathbf {r} }\right)\cdot \left(\mathbf {e} ^{(\mu )}(\mathbf {k} )a_{\mathbf {k} }^{(\mu )}(t)\,e^{i\mathbf {k} \cdot \mathbf {r} }-{\bar {\mathbf {e} }}^{(\mu )}(\mathbf {k} ){\bar {a}}_{\mathbf {k} }^{(\mu )}(t)\,e^{-i\mathbf {k} \cdot \mathbf {r} }\right)\mathrm {d} ^{3}\mathbf {r} \\&=-2\delta _{\mathbf {k'} ,\mathbf {k} }\delta _{\mu ',\mu }\,a_{\mathbf {k} }^{(\mu )}(t){\bar {a}}_{\mathbf {k} }^{(\mu )}(t)\;+\;\delta _{\mathbf {k'} ,-\mathbf {k} }\mathbf {e} ^{(\mu ')}(-\mathbf {k} )\cdot \mathbf {e} ^{(\mu )}(\mathbf {k} )\left[a_{\mathbf {k} }^{(\mu ')}(t)a_{\mathbf {-k} }^{(\mu )}(t)+{\bar {a}}_{\mathbf {k} }^{(\mu ')}(t){\bar {a}}_{\mathbf {-k} }^{(\mu )}(t)\right]\end{aligned}}$

The first term (linear in δ_k',k) will survive and will be added to the same term appearing in the expansion of B². It may be tempting to consider what the value of the following factor,

$\mathbf {e} ^{(\mu ')}(-\mathbf {k} )\cdot \mathbf {e} ^{(\mu )}(\mathbf {k} ),$

(appearing in the second term, linear in δ_k',−k) may be, but knowing its value is not necessary, because it cancels against the very same expression arising in the expansion of B².

In the expansion of B² the following appears as the factor of   $a_{\mathbf {k} }^{(\mu ')}(t)a_{\mathbf {-k} }^{(\mu )}(t)$  and  ${\bar {a}}_{\mathbf {k} }^{(\mu ')}(t){\bar {a}}_{\mathbf {-k} }^{(\mu )}(t):$

${\begin{aligned}\delta _{\mathbf {k} ,-\mathbf {k} '}&\left[\mathbf {k'} \times \mathbf {e} ^{(\mu ')}(\mathbf {k'} )\right]\cdot \left[\mathbf {k} \times \mathbf {e} ^{(\mu )}(\mathbf {k} )\right]=\delta _{\mathbf {k} ,-\mathbf {k} '}\left[-{\big (}\mathbf {k} \cdot \mathbf {k} {\big )}{\big (}\mathbf {e} ^{(\mu ')}(-\mathbf {k} )\cdot \mathbf {e} ^{(\mu )}(\mathbf {k} ){\big )}\right]-\left[-{\big (}\mathbf {k} \cdot \mathbf {e} ^{(\mu )}(\mathbf {k} ){\big )}{\big (}\mathbf {e} ^{(\mu ')}(-\mathbf {k} )\cdot \mathbf {k} ){\big )}\right]\\&=-k^{2}\;\delta _{\mathbf {k} ,-\mathbf {k} '}\;\mathbf {e} ^{(\mu ')}(-\mathbf {k} )\cdot \mathbf {e} ^{(\mu )}(\mathbf {k} )\end{aligned}}$

Recall that:

(i) In the expression of the field energy E² carries a factor ε₀ and B² the factor 1/μ₀.

(ii)  ${\frac {k^{2}}{\mu _{0}}}={\frac {\omega ^{2}}{c^{2}\mu _{0}}}={\frac {\epsilon _{0}\mu _{0}\omega ^{2}}{\mu _{0}}}=\epsilon _{0}\omega ^{2}.$

(iii) the Fourier expansion of E contains the factor ω.

All ingredients have been mentioned to arrive at the classical Hamiltonian in terms of Fourier coefficients,

$H=2V\epsilon _{0}\sum _{\mathbf {k} }\sum _{\mu =1,-1}\omega ^{2}a_{\mathbf {k} }^{(\mu )}(t){\bar {a}}_{\mathbf {k} }^{(\mu )}(t)\quad {\hbox{with}}\quad \mathbf {k} ={\frac {2\pi }{L}}(n_{x},\;n_{y},\;n_{z})\quad {\hbox{and}}\quad n_{x},\;n_{y},\;n_{z}=0,\;\pm 1,\;\pm 2,\ldots$

Fourier-expanded momentum

The electromagnetic momentum, P_EM, of EM radiation enclosed by a volume V is proportional to an integral of the Poynting vector S. In SI units:

${\begin{aligned}\mathbf {P} _{\textrm {EM}}&\equiv {\frac {1}{c^{2}}}\iiint _{V}\mathbf {S} \,{\textrm {d}}^{3}\mathbf {r} =\epsilon _{0}\iiint _{V}\mathbf {E} (\mathbf {r} ,t)\times \mathbf {B} (\mathbf {r} ,t)\,{\textrm {d}}^{3}\mathbf {r} \\&=2V\epsilon _{0}\sum _{\mathbf {k} }\sum _{\mu =1,-1}\omega \mathbf {k} a_{\mathbf {k} }^{(\mu )}(t){\bar {a}}_{\mathbf {k} }^{(\mu )}(t).\end{aligned}}$