Fourier expansion electromagnetic field: Difference between revisions
imported>Paul Wormer No edit summary |
imported>Milton Beychok m (Added a wiki link) |
||
(17 intermediate revisions by one other user not shown) | |||
Line 1: | Line 1: | ||
The [[electromagnetic wave|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 | {{subpages}} | ||
The [[electromagnetic wave|electromagnetic (EM) field]] is of importance as a carrier of solar energy and of 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 of the electromagnetic field| quantization]] that leads to [[photon]]s (light particles of well-defined energy and momentum). Further the Fourier transforms of EM fields play a role in theory of wave propagation through different media and light scattering. | |||
In the absence of charges and electric currents, both '''E''' and '''B''' can be derived from a third vector field, the [[vector potential]] '''A'''. | In the absence of charges and electric currents (that is, in [[Free space (electromagnetism)|"free" space]]), both '''E''' and '''B''' can be derived from a third vector field, the [[vector potential]] '''A'''. Below, the [[Fourier transform]] of the three real fields '''E''', '''B''', and '''A''' will be introduced. It will be seen that the expansion of the [[vector potential]] '''A''' yields the expansions of the 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== | ==Fourier expansion of a vector field== | ||
For | A real scalar function ''f''(''x''), with 0 ≤ ''x'' ≤ ''L'', has the following [[Fourier series|Fourier expansion]]: | ||
:<math> | |||
\begin{align} | |||
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\ge 0}\left[ f_n e^{2n\pi ix/L} + \bar{f}_n e^{-2n\pi ix/L}\right], | |||
\end{align} | |||
</math> | |||
where the bar indicates [[complex conjugation]] and the definition of the Fourier components is obvious. Note that the summand consists of a term plus its complex conjugate and hence is real. | |||
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: | |||
:<math> | |||
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 | |||
</math> | |||
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''<sup>3</sup>. 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''<sub>''x''</sub>, ''k''<sub>''y''</sub>, ''k''<sub>''z''</sub>) 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: | |||
:<math> | :<math> | ||
\mathbf{F}(\mathbf{r}, t) = \sum_\mathbf{k} \left( \mathbf{f}_k(t) e^{i\mathbf{k}\cdot\mathbf{r}} | \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) | + \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,\;\pm2,\ldots | |||
</math> | </math> | ||
with the following orthogonality relations,<ref>In many textbooks, the second and third non-orthogonality are overlooked, but since they happen to give terms that mutually cancel, the end result is always correct. However, the necessary proofs of these cancellations complicate the derivations.</ref> | |||
:<math> | :<math> | ||
\mathbf{k} = \ | \begin{align} | ||
\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{align} | |||
</math> | </math> | ||
Line 41: | Line 63: | ||
</math> | </math> | ||
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). | 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 | 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'''<sub>''x''</sub> and '''e'''<sub>''y''</sub>. It is often convenient to choose complex unit vectors obtained by a [[unitary transformation]], | ||
:<math> | :<math> | ||
\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) | \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) | ||
Line 49: | Line 69: | ||
which are orthonormal, | which are orthonormal, | ||
:<math> | :<math> | ||
\mathbf{e}^{(\mu)}\cdot\bar{\mathbf{e}}^{(\mu')} = \delta_{\mu,\mu'}\quad\hbox{with}\quad\mu,\mu'= 1,\, -1 | \mathbf{e}^{(\mu)}\cdot\bar{\mathbf{e}}^{(\mu')} = \delta_{\mu,\mu'}\quad\hbox{with}\quad\mu,\mu'= 1,\, -1, | ||
</math> | </math> | ||
and the bar indicating the complex conjugate. | |||
==Expansions== | ==Expansions== | ||
The Fourier expansion of the vector potential reads | The Fourier expansion of the vector potential reads | ||
Line 67: | Line 89: | ||
:<math> | :<math> | ||
-k^2 a^{(\mu)}_\mathbf{k}(t) = \frac{1}{c^2} \frac{\partial^2 a^{(\mu)}_\mathbf{k}(t)}{\partial t^2} | -k^2 a^{(\mu)}_\mathbf{k}(t) = \frac{1}{c^2} \frac{\partial^2 a^{(\mu)}_\mathbf{k}(t)}{\partial t^2} | ||
\quad \Longrightarrow \quad a^{(\mu)}_\mathbf{k}(t) \propto e^{-i\omega t}\quad\hbox{with}\quad\omega = kc. | \quad \Longrightarrow \quad a^{(\mu)}_\mathbf{k}(t) \propto e^{-i\omega t}\quad\hbox{with}\quad\omega = kc\quad\hbox{and}\quad k \equiv |\mathbf{k}|. | ||
</math> | </math> | ||
It is now an easy matter to construct the corresponding Fourier expansions for '''E''' and '''B''' from the expansion of the vector potential '''A'''. | 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, | |||
:<math> | :<math> | ||
\mathbf{E}(\mathbf{r}, t) = i\sum_\mathbf{k}\sum_{\mu=-1,1} \omega \left( | \mathbf{E}(\mathbf{r}, t) = i\sum_\mathbf{k}\sum_{\mu=-1,1} \omega \left( | ||
\mathbf{e}^{(\mu)}(\mathbf{k}) a^{(\mu)}_\mathbf{k}(t) \, e^{i\mathbf{k}\cdot\mathbf{r}} - | \mathbf{e}^{(\mu)}(\mathbf{k}) a^{(\mu)}_\mathbf{k}(t) \, e^{i\mathbf{k}\cdot\mathbf{r}} - | ||
\bar{\mathbf{e}}^{(\mu)}(\mathbf{k}) \bar{a}^{(\mu)}_\mathbf{k}(t) \, e^{-i\mathbf{k}\cdot\mathbf{r}} | \bar{\mathbf{e}}^{(\mu)}(\mathbf{k}) \bar{a}^{(\mu)}_\mathbf{k}(t) \, e^{-i\mathbf{k}\cdot\mathbf{r}} | ||
\right) | \right) . | ||
</math> | |||
The expansion for '''B''' follows by taking the curl, | |||
:<math> | |||
\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^{(\mu)}_\mathbf{k}(t) \, e^{i\mathbf{k}\cdot\mathbf{r}} - | |||
\big[\mathbf{k}\times\bar{\mathbf{e}}^{(\mu)}(\mathbf{k})\big]\; \bar{a}^{(\mu)}_\mathbf{k}(t) \, e^{-i\mathbf{k}\cdot\mathbf{r}} | |||
\right) . | |||
</math> | </math> | ||
==Fourier-expanded energy== | ==Fourier-expanded energy== | ||
The electromagnetic energy density <math> | The electromagnetic energy density is | ||
:<math> | |||
\mathcal{E}_\mathrm{Field} = \frac{1}{2} ( \epsilon_0 \mathbf{E}\cdot\mathbf{E}+ \frac{1}{\mu_0}\mathbf{B}\cdot\mathbf{B}), | |||
</math> | |||
where μ<sub>0</sub> is the [[magnetic constant]]. | |||
The total energy (classical Hamiltonian) of a finite volume ''V'' is defined by | |||
:<math> | :<math> | ||
H = \iiint_V \mathcal{E}_\mathrm{Field}(\mathbf{r},t) \mathrm{d}^3\mathbf{r}. | H = \iiint_V \mathcal{E}_\mathrm{Field}(\mathbf{r},t) \mathrm{d}^3\mathbf{r}. | ||
</math> | </math> | ||
Use in the expansion of '''E'''⋅'''E''' = '''E'''<sup>2</sup> the following | |||
:<math> | |||
\begin{align} | |||
\frac{1}{V}\iiint_V | |||
&\left( | |||
\mathbf{e}^{(\mu')}(\mathbf{k'}) a^{(\mu')}_\mathbf{k'}(t) \, e^{i\mathbf{k'}\cdot\mathbf{r}} - | |||
\bar{\mathbf{e}}^{(\mu')}(\mathbf{k'}) \bar{a}^{(\mu')}_\mathbf{k'}(t) \, e^{-i\mathbf{k'}\cdot\mathbf{r}} | |||
\right) | |||
\cdot | |||
\left( | |||
\mathbf{e}^{(\mu)}(\mathbf{k}) a^{(\mu)}_\mathbf{k}(t) \, e^{i\mathbf{k}\cdot\mathbf{r}} - | |||
\bar{\mathbf{e}}^{(\mu)}(\mathbf{k}) \bar{a}^{(\mu)}_\mathbf{k}(t) \, e^{-i\mathbf{k}\cdot\mathbf{r}} | |||
\right)\mathrm{d}^3\mathbf{r} \\ | |||
&= | |||
-2\delta_{\mathbf{k'},\mathbf{k}}\delta_{\mu',\mu}\,a^{(\mu)}_\mathbf{k}(t) \bar{a}^{(\mu)}_\mathbf{k}(t)\;+\; \delta_{\mathbf{k'},-\mathbf{k}} \mathbf{e}^{(\mu')}(-\mathbf{k})\cdot\mathbf{e}^{(\mu)}(\mathbf{k}) | |||
\left[a^{(\mu')}_\mathbf{k}(t) a^{(\mu)}_\mathbf{-k}(t) | |||
+ \bar{a}^{(\mu')}_\mathbf{k}(t) \bar{a}^{(\mu)}_\mathbf{-k}(t)\right] | |||
\end{align} | |||
</math> | |||
The first term (linear in δ<sub><b>k'</b>,'''k'''</sub>) will survive and will be added to the same term appearing in the expansion of '''B'''<sup>2</sup>. It may be tempting to consider what the value will be of the following factor, | |||
:<math> | |||
\mathbf{e}^{(\mu')}(-\mathbf{k})\cdot\mathbf{e}^{(\mu)}(\mathbf{k}), | |||
</math> | |||
(appearing in the second term, linear in δ<sub><b>k'</b>,−'''k'''</sub>), but this is not needed, because it cancels against the very same expression arising in the expansion of '''B'''<sup>2</sup>. | |||
In the expansion of '''B'''<sup>2</sup> the following appears as the factor of | |||
<font style="vertical-align: text-bottom;"> | |||
<math>a^{(\mu')}_\mathbf{k}(t) a^{(\mu)}_\mathbf{-k}(t)</math></font> and <font style="vertical-align: text-bottom;"><math>\bar{a}^{(\mu')}_\mathbf{k}(t) \bar{a}^{(\mu)}_\mathbf{-k}(t): </math></font> | |||
:<math> | :<math> | ||
\begin{align} | |||
\big(\ | \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{align} | |||
</math> | </math> | ||
Recall that: | |||
(i) The expression of the field energy '''E'''<sup>2</sup> contains a factor ε<sub>0</sub> and '''B'''<sup>2</sup> the factor 1/μ<sub>0</sub>. | |||
(ii) <math> | |||
\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. | |||
</math> | |||
(iii) The Fourier expansion of '''E''' contains the factor ω. | |||
This means that all ingredients have been mentioned to arrive at the classical Hamiltonian in terms of Fourier coefficients, | |||
:<math> | :<math> | ||
\ | H = 2V\epsilon_0 \sum_\mathbf{k}\sum_{\mu=1,-1} \omega^2 | ||
a^{(\mu)}_\mathbf{k}(t)\bar{a}^{(\mu)}_\mathbf{k}(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,\;\pm1,\;\pm2,\ldots | |||
</math> | </math> | ||
==Fourier-expanded momentum== | ==Fourier-expanded momentum== | ||
The electromagnetic [[momentum]], '''P'''<sub>EM</sub>, of EM radiation enclosed by a volume ''V'' is proportional to an integral of the Poynting vector | The electromagnetic [[momentum]], '''P'''<sub>EM</sub>, of EM radiation enclosed by a volume ''V'' is proportional to an integral of the Poynting vector '''S'''. In [[SI]] units: | ||
:<math> | |||
\begin{align} | |||
\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}\\ | |||
\end{align} | |||
</math> | |||
To show that the expansion terms containing | |||
:<math> | |||
a^{(\mu')}_\mathbf{k'}(t)\;a^{(\mu)}_\mathbf{k}(t) | |||
\iiint_V e^{i(\mathbf{k'}+\mathbf{k})\cdot\mathbf{r}} \;\mathrm{d}^3\mathbf{r} | |||
\quad\hbox{and}\quad | |||
\bar{a}^{(\mu')}_\mathbf{k'}(t)\;\bar{a}^{(\mu)}_\mathbf{k}(t) | |||
\iiint_V e^{-i (\mathbf{k'}+\mathbf{k})\cdot\mathbf{r}} \;\mathrm{d}^3\mathbf{r} | |||
</math> | |||
cancel, we write | |||
:<math> | |||
\mathbf{P}_\textrm{EM} = \frac{\epsilon_0}{2} | |||
\iiint_V\left( \mathbf{E}(\mathbf{r},t)\times \mathbf{B}(\mathbf{r},t)- \mathbf{B}(\mathbf{r},t)\times \mathbf{E}(\mathbf{r},t)\right)\, \textrm{d}^3\mathbf{r} | |||
</math> | |||
Upon expansion, one finds among other terms, | |||
:<math> | |||
\delta_{\mathbf{k'},-\mathbf{k}} a^{(\mu')}_\mathbf{k'}(t)\;a^{(\mu)}_\mathbf{k}(t) | |||
\Big[ | |||
\mathbf{e}^{(\mu')}(\mathbf{k'})\times\left(\mathbf{k}\times \mathbf{e}^{(\mu)}(\mathbf{k}) \right) | |||
-\left(\mathbf{k'}\times \mathbf{e}^{(\mu')}(\mathbf{k'}) \right)\times\mathbf{e}^{(\mu)}(\mathbf{k}) | |||
\Big] | |||
</math> | |||
Twice the [[Baccab formula]] and using that "under δ<sub><b>k'</b>,−'''k'''</sub>" '''k''' and −<b>k'</b> maybe interchanged, and an obvious short-hand notation gives | |||
:<math> | |||
\begin{align} | |||
\mathbf{e}'\times(\mathbf{k}\times \mathbf{e}) &= | |||
\mathbf{k} (\mathbf{e}'\cdot \mathbf{e}) - \mathbf{e} (\mathbf{e}'\cdot \mathbf{k})= | |||
\mathbf{k} (\mathbf{e}'\cdot \mathbf{e}) + \mathbf{e} (\mathbf{e}'\cdot \mathbf{k}') | |||
=\mathbf{k} (\mathbf{e}'\cdot \mathbf{e}) \\ | |||
-(\mathbf{k}'\times \mathbf{e}')\times\mathbf{e} &= \mathbf{e}\times(\mathbf{k}'\times \mathbf{e}')= | |||
\mathbf{k}' (\mathbf{e}\cdot \mathbf{e}') - \mathbf{e}' (\mathbf{e}\cdot \mathbf{k}')= | |||
\mathbf{k}' (\mathbf{e}\cdot \mathbf{e}') + \mathbf{e}' (\mathbf{e}\cdot \mathbf{k}) | |||
=\mathbf{k}' (\mathbf{e}\cdot \mathbf{e}') . \\ | |||
\end{align} | |||
</math> | |||
And | |||
:<math> | |||
\delta_{\mathbf{k'},-\mathbf{k}}\;\left[ \mathbf{k} (\mathbf{e}'\cdot \mathbf{e}) + \mathbf{k}' (\mathbf{e}\cdot \mathbf{e}') \right] = 0 . | |||
</math> | |||
The other terms in the expansion do not offer any new problems, so that, finally, | |||
:<math> | :<math> | ||
\mathbf{P}_\textrm{EM} \ | \mathbf{P}_\textrm{EM} = 2 V \epsilon_0 \sum_\mathbf{k}\sum_{\mu=1,-1} \omega \mathbf{k} | ||
\ | a^{(\mu)}_\mathbf{k}(t)\bar{a}^{(\mu)}_\mathbf{k}(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,\;\pm1,\;\pm2,\ldots. | |||
</math> | </math> | ||
==Notes== | |||
<references /> |
Latest revision as of 23:07, 30 November 2010
The electromagnetic (EM) field is of importance as a carrier of solar energy and of 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 transforms of EM fields play a role in theory of wave propagation through different media and 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 vector field, the vector potential A. Below, the Fourier transform of the three real fields E, B, and A will be introduced. It will be seen that the expansion of the vector potential A yields the expansions of the 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:
where the bar indicates complex conjugation and the definition of the Fourier components is obvious. Note that the summand consists of a term plus its complex conjugate and hence is real. 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:
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 = L3. 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 (kx, ky, kz) 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:
with the following orthogonality relations,[1]
Vector potential
The magnetic field B satisfies the following Maxwell equation:
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
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,
because it is assumed that charge distributions ρ are zero. The quantity ε0 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,
The fact that E can be written this way is due to the choice of Coulomb gauge for A:
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, ex and ey. It is often convenient to choose complex unit vectors obtained by a unitary transformation,
which are orthonormal,
- Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle \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
- Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle \mathbf{A}(\mathbf{r}, t) = \sum_\mathbf{k}\sum_{\mu=-1,1} \left( \mathbf{e}^{(\mu)}(\mathbf{k}) a^{(\mu)}_\mathbf{k}(t) \, e^{i\mathbf{k}\cdot\mathbf{r}} + \bar{\mathbf{e}}^{(\mu)}(\mathbf{k}) \bar{a}^{(\mu)}_\mathbf{k}(t) \, e^{-i\mathbf{k}\cdot\mathbf{r}} \right). }
The vector potential obeys the wave equation,
- Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle \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,
- Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle -k^2 a^{(\mu)}_\mathbf{k}(t) = \frac{1}{c^2} \frac{\partial^2 a^{(\mu)}_\mathbf{k}(t)}{\partial t^2} \quad \Longrightarrow \quad a^{(\mu)}_\mathbf{k}(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,
- Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle \mathbf{E}(\mathbf{r}, t) = i\sum_\mathbf{k}\sum_{\mu=-1,1} \omega \left( \mathbf{e}^{(\mu)}(\mathbf{k}) a^{(\mu)}_\mathbf{k}(t) \, e^{i\mathbf{k}\cdot\mathbf{r}} - \bar{\mathbf{e}}^{(\mu)}(\mathbf{k}) \bar{a}^{(\mu)}_\mathbf{k}(t) \, e^{-i\mathbf{k}\cdot\mathbf{r}} \right) . }
The expansion for B follows by taking the curl,
- Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle \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^{(\mu)}_\mathbf{k}(t) \, e^{i\mathbf{k}\cdot\mathbf{r}} - \big[\mathbf{k}\times\bar{\mathbf{e}}^{(\mu)}(\mathbf{k})\big]\; \bar{a}^{(\mu)}_\mathbf{k}(t) \, e^{-i\mathbf{k}\cdot\mathbf{r}} \right) . }
Fourier-expanded energy
The electromagnetic energy density is
- Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle \mathcal{E}_\mathrm{Field} = \frac{1}{2} ( \epsilon_0 \mathbf{E}\cdot\mathbf{E}+ \frac{1}{\mu_0}\mathbf{B}\cdot\mathbf{B}), }
where μ0 is the magnetic constant. The total energy (classical Hamiltonian) of a finite volume V is defined by
Use in the expansion of E⋅E = E2 the following
- Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle \begin{align} \frac{1}{V}\iiint_V &\left( \mathbf{e}^{(\mu')}(\mathbf{k'}) a^{(\mu')}_\mathbf{k'}(t) \, e^{i\mathbf{k'}\cdot\mathbf{r}} - \bar{\mathbf{e}}^{(\mu')}(\mathbf{k'}) \bar{a}^{(\mu')}_\mathbf{k'}(t) \, e^{-i\mathbf{k'}\cdot\mathbf{r}} \right) \cdot \left( \mathbf{e}^{(\mu)}(\mathbf{k}) a^{(\mu)}_\mathbf{k}(t) \, e^{i\mathbf{k}\cdot\mathbf{r}} - \bar{\mathbf{e}}^{(\mu)}(\mathbf{k}) \bar{a}^{(\mu)}_\mathbf{k}(t) \, e^{-i\mathbf{k}\cdot\mathbf{r}} \right)\mathrm{d}^3\mathbf{r} \\ &= -2\delta_{\mathbf{k'},\mathbf{k}}\delta_{\mu',\mu}\,a^{(\mu)}_\mathbf{k}(t) \bar{a}^{(\mu)}_\mathbf{k}(t)\;+\; \delta_{\mathbf{k'},-\mathbf{k}} \mathbf{e}^{(\mu')}(-\mathbf{k})\cdot\mathbf{e}^{(\mu)}(\mathbf{k}) \left[a^{(\mu')}_\mathbf{k}(t) a^{(\mu)}_\mathbf{-k}(t) + \bar{a}^{(\mu')}_\mathbf{k}(t) \bar{a}^{(\mu)}_\mathbf{-k}(t)\right] \end{align} }
The first term (linear in δk',k) will survive and will be added to the same term appearing in the expansion of B2. It may be tempting to consider what the value will be of the following factor,
- Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle \mathbf{e}^{(\mu')}(-\mathbf{k})\cdot\mathbf{e}^{(\mu)}(\mathbf{k}), }
(appearing in the second term, linear in δk',−k), but this is not needed, because it cancels against the very same expression arising in the expansion of B2.
In the expansion of B2 the following appears as the factor of Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle a^{(\mu')}_\mathbf{k}(t) a^{(\mu)}_\mathbf{-k}(t)} and Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle \bar{a}^{(\mu')}_\mathbf{k}(t) \bar{a}^{(\mu)}_\mathbf{-k}(t): }
- Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle \begin{align} \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{align} }
Recall that:
(i) The expression of the field energy E2 contains a factor ε0 and B2 the factor 1/μ0.
(ii) Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle \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 ω.
This means that all ingredients have been mentioned to arrive at the classical Hamiltonian in terms of Fourier coefficients,
Fourier-expanded momentum
The electromagnetic momentum, PEM, of EM radiation enclosed by a volume V is proportional to an integral of the Poynting vector S. In SI units:
- Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle \begin{align} \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}\\ \end{align} }
To show that the expansion terms containing
- Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle a^{(\mu')}_\mathbf{k'}(t)\;a^{(\mu)}_\mathbf{k}(t) \iiint_V e^{i(\mathbf{k'}+\mathbf{k})\cdot\mathbf{r}} \;\mathrm{d}^3\mathbf{r} \quad\hbox{and}\quad \bar{a}^{(\mu')}_\mathbf{k'}(t)\;\bar{a}^{(\mu)}_\mathbf{k}(t) \iiint_V e^{-i (\mathbf{k'}+\mathbf{k})\cdot\mathbf{r}} \;\mathrm{d}^3\mathbf{r} }
cancel, we write
- Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle \mathbf{P}_\textrm{EM} = \frac{\epsilon_0}{2} \iiint_V\left( \mathbf{E}(\mathbf{r},t)\times \mathbf{B}(\mathbf{r},t)- \mathbf{B}(\mathbf{r},t)\times \mathbf{E}(\mathbf{r},t)\right)\, \textrm{d}^3\mathbf{r} }
Upon expansion, one finds among other terms,
- Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle \delta_{\mathbf{k'},-\mathbf{k}} a^{(\mu')}_\mathbf{k'}(t)\;a^{(\mu)}_\mathbf{k}(t) \Big[ \mathbf{e}^{(\mu')}(\mathbf{k'})\times\left(\mathbf{k}\times \mathbf{e}^{(\mu)}(\mathbf{k}) \right) -\left(\mathbf{k'}\times \mathbf{e}^{(\mu')}(\mathbf{k'}) \right)\times\mathbf{e}^{(\mu)}(\mathbf{k}) \Big] }
Twice the Baccab formula and using that "under δk',−k" k and −k' maybe interchanged, and an obvious short-hand notation gives
- Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle \begin{align} \mathbf{e}'\times(\mathbf{k}\times \mathbf{e}) &= \mathbf{k} (\mathbf{e}'\cdot \mathbf{e}) - \mathbf{e} (\mathbf{e}'\cdot \mathbf{k})= \mathbf{k} (\mathbf{e}'\cdot \mathbf{e}) + \mathbf{e} (\mathbf{e}'\cdot \mathbf{k}') =\mathbf{k} (\mathbf{e}'\cdot \mathbf{e}) \\ -(\mathbf{k}'\times \mathbf{e}')\times\mathbf{e} &= \mathbf{e}\times(\mathbf{k}'\times \mathbf{e}')= \mathbf{k}' (\mathbf{e}\cdot \mathbf{e}') - \mathbf{e}' (\mathbf{e}\cdot \mathbf{k}')= \mathbf{k}' (\mathbf{e}\cdot \mathbf{e}') + \mathbf{e}' (\mathbf{e}\cdot \mathbf{k}) =\mathbf{k}' (\mathbf{e}\cdot \mathbf{e}') . \\ \end{align} }
And
- Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle \delta_{\mathbf{k'},-\mathbf{k}}\;\left[ \mathbf{k} (\mathbf{e}'\cdot \mathbf{e}) + \mathbf{k}' (\mathbf{e}\cdot \mathbf{e}') \right] = 0 . }
The other terms in the expansion do not offer any new problems, so that, finally,
Notes
- ↑ In many textbooks, the second and third non-orthogonality are overlooked, but since they happen to give terms that mutually cancel, the end result is always correct. However, the necessary proofs of these cancellations complicate the derivations.