Peeter Joot's (OLD) Blog.

Math, physics, perl, and programming obscurity.

PHY450H1S. Relativistic Electrodynamics Lecture 13 (Taught by Prof. Erich Poppitz). Variational principle for the field.

Posted by peeterjoot on February 22, 2011

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


Covering chapter 4 material from the text [1].

Covering lecture notes pp.103-113: variational principle for the electromagnetic field and the relevant boundary conditions (103-105); the second set of Maxwell’s equations from the variational principle (106-108); Maxwell’s equations in vacuum and the wave equation in the non-relativistic Coulomb gauge (109-111)

Review. Our action.

\begin{aligned}S&= S_{\text{particles}} + S_{\text{interaction}} + S_{\text{EM field}}&= \sum_A \int_{x_A^i(\tau)} ds ( -m_A c )- \sum_A\frac{e_A}{c}\int dx_A^i A_i(x_A)- \frac{1}{{16 \pi c}} \int d^4 x F^{ij } F_{ij}.\end{aligned}

Our dynamics variables are

\begin{aligned}\left\{\begin{array}{l l}x_A^i(\tau) & \quad \mbox{A = 1, \cdots, N$} \\ A^i(x) & \quad \mbox{$A = 1, \cdots, N$}\end{array}\right.\end{aligned} \hspace{\stretch{1}}(2.1)$

We saw that the interaction term could also be written in terms of a delta function current, with

\begin{aligned}S_{\text{interaction}}= -\frac{1}{{c^2}} \int d^4x j^i(x) A_i(x),\end{aligned} \hspace{\stretch{1}}(2.2)


\begin{aligned}j^i(x) = \sum_A c e_A \int dx_A^i \delta^4( x - x_A(\tau)).\end{aligned} \hspace{\stretch{1}}(2.3)

Variation with respect to x_A^i(\tau) gave us

\begin{aligned}m c \frac{d{{u^i_A}}}{ds} = \frac{e}{c} u_A^j F_{ij}.\end{aligned} \hspace{\stretch{1}}(2.4)

Note that it’s easy to get the sign mixed up here. With our (+,-,-,-) metric tensor, if the second index is the summation index, we have a positive sign.

Only the S_{\text{particles}} and S_{\text{interaction}} depend on x_A^i(\tau).

The field action variation.

\paragraph{Today:} We’ll find the EOM for A^i(x). The dynamical degrees of freedom are A^i(\mathbf{x},t)

\begin{aligned}S[A^i(\mathbf{x}, t)] = -\frac{1}{{16 \pi c}} \int d^4x F_{ij}F^{ij} - \frac{1}{{c^2}} \int d^4 x A^i j_i.\end{aligned} \hspace{\stretch{1}}(3.5)

Here j^i are treated as “sources”.

We demand that

\begin{aligned}\delta S = S[ A^i(\mathbf{x}, t) + \delta A^i(\mathbf{x}, t)] - S[ A^i(\mathbf{x}, t) ] = 0 + O(\delta A)^2.\end{aligned} \hspace{\stretch{1}}(3.6)

We need to impose two conditions.
\item At spatial \infty, i.e. at {\left\lvert{\mathbf{x}}\right\rvert} \rightarrow \infty, \forall t, we’ll impose the condition

\begin{aligned}{\left.{{A^i(\mathbf{x}, t)}}\right\vert}_{{{\left\lvert{\mathbf{x}}\right\rvert} \rightarrow \infty}} \rightarrow 0.\end{aligned} \hspace{\stretch{1}}(3.7)

This is sensible, because fields are created by charges, and charges are assumed to be localized in a bounded region. The field outside charges will \rightarrow 0 at {\left\lvert{\mathbf{x}}\right\rvert} \rightarrow \infty. Later we will treat the integration range as finite, and bounded, then later allow the boundary to go to infinity.

\item at t = -T and t = T we’ll imagine that the values of A^i(\mathbf{x}, \pm T) are fixed.

This is analogous to x(t_i) = x_1 and x(t_f) = x_2 in particle mechanics.

Since A^i(\mathbf{x}, \pm T) is given, and equivalent to the initial and final field configurations, our extremes at the boundary is zero

\begin{aligned}\delta A^i(\mathbf{x}, \pm T) = 0.\end{aligned} \hspace{\stretch{1}}(3.8)


PICTURE: a cylinder in spacetime, with an attempt to depict the boundary.

Computing the variation.

\begin{aligned}\delta S[A^i(\mathbf{x}, t)]= -\frac{1}{{16 \pi c}} \int d^4 x \delta (F_{ij}F^{ij}) - \frac{1}{{c^2}} \int d^4 x \delta(A^i) j_i.\end{aligned} \hspace{\stretch{1}}(4.9)

Looking first at the variation of just the F^2 bit we have

\begin{aligned}\delta (F_{ij}F^{ij})&=\delta(F_{ij}) F^{ij} + F_{ij} \delta(F^{ij}) \\ &=2 \delta(F^{ij}) F_{ij} \\ &=2 \delta(\partial^i A^j - \partial^j A^i) F_{ij} \\ &=2 \delta(\partial^i A^j) F_{ij} - 2 \delta(\partial^j A^i) F_{ij} \\ &=2 \delta(\partial^i A^j) F_{ij} - 2 \delta(\partial^i A^j) F_{ji} \\ &=4 \delta(\partial^i A^j) F_{ij} \\ &=4 F_{ij} \partial^i \delta(A^j).\end{aligned}

Our variation is now reduced to

\begin{aligned}\delta S[A^i(\mathbf{x}, t)]&= -\frac{1}{{4 \pi c}} \int d^4 x F_{ij} \partial^i \delta(A^j) - \frac{1}{{c^2}} \int d^4 x j^i \delta(A_i) \\ &= -\frac{1}{{4 \pi c}} \int d^4 x F^{ij} \frac{\partial {}}{\partial {x^i}} \delta(A_j) - \frac{1}{{c^2}} \int d^4 x j^i \delta(A_i).\end{aligned}

We can integrate this first term by parts

\begin{aligned}\int d^4 x F^{ij} \frac{\partial {}}{\partial {x^i}} \delta(A_j)&=\int d^4 x \frac{\partial {}}{\partial {x^i}} \left( F^{ij} \delta(A_j) \right)-\int d^4 x \left( \frac{\partial {}}{\partial {x^i}} F^{ij} \right) \delta(A_j) \end{aligned}

The first term is a four dimensional divergence, with the contraction of the four gradient \partial_i with a four vector B^i = F^{ij} \delta(A_j).

Prof. Poppitz chose dx^0 d^3 \mathbf{x} split of d^4 x to illustrate that this can be viewed as regular old spatial three vector divergences. It is probably more rigorous to mandate that the four volume element is oriented d^4 x = (1/4!)\epsilon_{ijkl} dx^i dx^j dx^k dx^l, and then utilize the 4D version of the divergence theorem (or its Stokes Theorem equivalent). The completely antisymmetric tensor should do most of the work required to express the oriented boundary volume.

Because we have specified that A^i is zero on the boundary, so is F^{ij}, so these boundary terms are killed off. We are left with

\begin{aligned}\delta S[A^i(\mathbf{x}, t)]&= -\frac{1}{{4 \pi c}} \int d^4 x \delta (A_j) \partial_i F^{ij} - \frac{1}{{c^2}} \int d^4 x j^i \delta(A_i) \\ &=\int d^4 x \delta A_j(x)\left(-\frac{1}{{4 \pi c}} \partial_i F^{ij}(x) - \frac{1}{{c^2}} j^i\right)  \\ &= 0.\end{aligned}

This gives us

\begin{aligned}\boxed{\partial_i F^{ij} = \frac{4 \pi}{c} j^j}\end{aligned} \hspace{\stretch{1}}(4.10)

Unpacking these.

Recall that the Bianchi identity

\begin{aligned}\epsilon^{ijkl} \partial_j F_{kl} = 0,\end{aligned} \hspace{\stretch{1}}(5.11)

gave us

\begin{aligned}\boldsymbol{\nabla} \cdot \mathbf{B} &= 0 \\ \boldsymbol{\nabla} \times \mathbf{E} &= -\frac{1}{{c}} \frac{\partial {\mathbf{B}}}{\partial {t}}.\end{aligned} \hspace{\stretch{1}}(5.12)

How about the EOM that we have found by varying the action? One of those equations is

\begin{aligned}\partial_\alpha F^{\alpha 0} = \frac{4 \pi}{c} j^0 = 4 \pi \rho,\end{aligned} \hspace{\stretch{1}}(5.14)

since j^0 = c \rho.


\begin{aligned}F^{\alpha 0} = (\mathbf{E})^\alpha,\end{aligned} \hspace{\stretch{1}}(5.15)

we have

\begin{aligned}\boldsymbol{\nabla} \cdot \mathbf{E} = 4 \pi \rho.\end{aligned} \hspace{\stretch{1}}(5.16)

The messier one to deal with is

\begin{aligned}\partial_i F^{i\alpha} = \frac{4 \pi}{c} j^\alpha.\end{aligned} \hspace{\stretch{1}}(5.17)

Splitting out the spatial and time indexes for the four gradient we have

\begin{aligned}\partial_i F^{i\alpha}&= \partial_\beta F^{\beta \alpha} + \partial_0 F^{0 \alpha} \\ &= \partial_\beta F^{\beta \alpha} - \frac{1}{{c}} \frac{\partial {(\mathbf{E})^\alpha}}{\partial {t}} \\ \end{aligned}

The spatial index tensor element is

\begin{aligned}F^{\beta \alpha} &= \partial^\beta A^\alpha - \partial^\alpha A^\beta  \\ &= - \frac{\partial {A^\alpha}}{\partial {x^\beta}} + \frac{\partial {A^\beta}}{\partial {x^\alpha}} \\ &= \epsilon^{\alpha\beta\gamma} B^\gamma,\end{aligned}

so the sum becomes

\begin{aligned}\partial_i F^{i\alpha}&= \partial_\beta ( \epsilon^{\alpha\beta\gamma} B^\gamma) - \frac{1}{{c}} \frac{\partial {(\mathbf{E})^\alpha}}{\partial {t}} \\ &= \epsilon^{\beta\gamma\alpha} \partial_\beta B^\gamma - \frac{1}{{c}} \frac{\partial {(\mathbf{E})^\alpha}}{\partial {t}} \\ &= (\boldsymbol{\nabla} \times \mathbf{B})^\alpha - \frac{1}{{c}} \frac{\partial {(\mathbf{E})^\alpha}}{\partial {t}}.\end{aligned}

This gives us

\begin{aligned}\frac{4 \pi}{c} j^\alpha= (\boldsymbol{\nabla} \times \mathbf{B})^\alpha - \frac{1}{{c}} \frac{\partial {(\mathbf{E})^\alpha}}{\partial {t}},\end{aligned} \hspace{\stretch{1}}(5.18)

or in vector form

\begin{aligned}\boldsymbol{\nabla} \times \mathbf{B} - \frac{1}{{c}} \frac{\partial {\mathbf{E}}}{\partial {t}} = \frac{4 \pi}{c} \mathbf{j}.\end{aligned} \hspace{\stretch{1}}(5.19)

Summarizing what we know so far, we have

\begin{aligned}\boxed{\begin{aligned}\partial_i F^{ij} &= \frac{4 \pi}{c} j^j \\ \epsilon^{ijkl} \partial_j F_{kl} &= 0\end{aligned}}\end{aligned} \hspace{\stretch{1}}(5.20)

or in vector form

\begin{aligned}\boxed{\begin{aligned}\boldsymbol{\nabla} \cdot \mathbf{E} &= 4 \pi \rho \\ \boldsymbol{\nabla} \times \mathbf{B} -\frac{1}{{c}} \frac{\partial {\mathbf{E}}}{\partial {t}} &= \frac{4 \pi}{c} \mathbf{j} \\ \boldsymbol{\nabla} \cdot \mathbf{B} &= 0 \\ \boldsymbol{\nabla} \times \mathbf{E} +\frac{1}{{c}} \frac{\partial {\mathbf{B}}}{\partial {t}} &= 0\end{aligned}}\end{aligned} \hspace{\stretch{1}}(5.21)

Speed of light

\paragraph{Claim}: “c” is the speed of EM waves in vacuum.

Study equations in vacuum (no sources, so j^i = 0) for A^i = (\phi, \mathbf{A}).

\begin{aligned}\boldsymbol{\nabla} \cdot \mathbf{E} &= 0 \\ \boldsymbol{\nabla} \times \mathbf{B} &= \frac{1}{{c}} \frac{\partial {\mathbf{E}}}{\partial {t}}\end{aligned} \hspace{\stretch{1}}(6.22)


\begin{aligned}\mathbf{E} &= - \boldsymbol{\nabla} \phi - \frac{1}{{c}} \frac{\partial {\mathbf{A}}}{\partial {t}} \\ \mathbf{B} &= \boldsymbol{\nabla} \times \mathbf{A}\end{aligned} \hspace{\stretch{1}}(6.24)

In terms of potentials

\begin{aligned}0 &= \boldsymbol{\nabla} \times (\boldsymbol{\nabla} \times \mathbf{A}) \\ &= \boldsymbol{\nabla} \times \mathbf{B} \\ &= \frac{1}{{c}} \frac{\partial {\mathbf{E}}}{\partial {t}} \\ &= \frac{1}{{c}} \frac{\partial {}}{\partial {t}} \left( - \boldsymbol{\nabla} \phi - \frac{1}{{c}} \frac{\partial {\mathbf{A}}}{\partial {t}} \right) \\ &= -\frac{1}{{c}} \frac{\partial {}}{\partial {t}} \boldsymbol{\nabla} \phi - \frac{1}{{c^2}} \frac{\partial^2 \mathbf{A}}{\partial t^2} \end{aligned}

Since we also have

\begin{aligned}\boldsymbol{\nabla} \times (\boldsymbol{\nabla} \times \mathbf{A}) = \boldsymbol{\nabla} (\boldsymbol{\nabla} \cdot \mathbf{A}) - \boldsymbol{\nabla}^2 \mathbf{A},\end{aligned} \hspace{\stretch{1}}(6.26)

some rearrangement gives

\begin{aligned}\boldsymbol{\nabla} (\boldsymbol{\nabla} \cdot \mathbf{A}) = \boldsymbol{\nabla}^2 \mathbf{A}  -\frac{1}{{c}} \frac{\partial {}}{\partial {t}} \boldsymbol{\nabla} \phi - \frac{1}{{c^2}} \frac{\partial^2 \mathbf{A}}{\partial t^2}.\end{aligned} \hspace{\stretch{1}}(6.27)

The remaining equation \boldsymbol{\nabla} \cdot \mathbf{E} = 0, in terms of potentials is

\begin{aligned}\boldsymbol{\nabla} \cdot \mathbf{E} = - \boldsymbol{\nabla}^2 \phi - \frac{1}{{c}} \frac{\partial {\boldsymbol{\nabla} \cdot \mathbf{A}}}{\partial {t}} \end{aligned} \hspace{\stretch{1}}(6.28)

We can make a gauge transformation that completely eliminates 6.28, and reduces 6.27 to a wave equation.

\begin{aligned}(\phi, \mathbf{A}) \rightarrow (\phi', \mathbf{A}')\end{aligned} \hspace{\stretch{1}}(6.29)


\begin{aligned}\phi &= \phi' - \frac{1}{{c}} \frac{\partial {\chi}}{\partial {t}} \\ \mathbf{A} &= \mathbf{A}' + \boldsymbol{\nabla} \chi\end{aligned} \hspace{\stretch{1}}(6.30)

Can choose \chi(\mathbf{x}, t) to make \phi' = 0 (\forall \phi \exists \chi, \phi' = 0)

\begin{aligned}\frac{1}{{c}} \frac{\partial {}}{\partial {t}} \chi(\mathbf{x}, t) = \phi(\mathbf{x}, t)\end{aligned} \hspace{\stretch{1}}(6.32)

\begin{aligned}\chi(\mathbf{x}, t) = c \int_{-\infty}^t dt' \phi(\mathbf{x}, t')\end{aligned} \hspace{\stretch{1}}(6.33)

Can also find a transformation that also allows \boldsymbol{\nabla} \cdot \mathbf{A} = 0

\paragraph{Q:} What would that second transformation be explicitly?
\paragraph{A:} To be revisited next lecture, when this is covered in full detail.

This is the Coulomb gauge

\begin{aligned}\phi &= 0 \\ \boldsymbol{\nabla} \cdot \mathbf{A} &= 0\end{aligned} \hspace{\stretch{1}}(6.34)

From 6.27, we then have

\begin{aligned}\frac{1}{{c^2}} \frac{\partial^2 \mathbf{A}'}{\partial t^2} -\boldsymbol{\nabla}^2 \mathbf{A}' = 0\end{aligned} \hspace{\stretch{1}}(6.36)

which is the wave equation for the propagation of the vector potential \mathbf{A}'(\mathbf{x}, t) through space at velocity c, confirming that c is the speed of electromagnetic propagation (the speed of light).


[1] L.D. Landau and E.M. Lifshitz. The classical theory of fields. Butterworth-Heinemann, 1980.


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 )

Twitter picture

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

Facebook photo

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

Google+ photo

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

Connecting to %s

%d bloggers like this: