The Lagrangian for Electromagnetic Fields

There are not many ways to make a scalar Lagrangian from the field tensor. We already know that

\begin{displaymath}\bgroup\color{black} {\partial F_{\mu\nu}\over\partial x_\nu}={j_\mu\over c} \egroup\end{displaymath}

and we need to make our Lagrangian out of the fields, not just the current. Again, \bgroup\color{black}$x_\mu$\egroup cannot appear explicitly because that violates symmetries of nature. Also we want a linear equation and so higher powers of the field should not occur. A term of the form \bgroup\color{black}$mA_\mu A_\mu$\egroup is a mass term and would cause fields to fall off faster than \bgroup\color{black}${1\over r}$\egroup. So, the only reasonable choice is

\begin{displaymath}\bgroup\color{black} F_{\mu\nu}F_{\mu\nu}= 2(B^2-E^2). \egroup\end{displaymath}

One might consider

\begin{displaymath}\bgroup\color{black} e_{\mu\nu\lambda\sigma}F_{\mu\nu}F_{\lambda\sigma}=\vec{B}\cdot\vec{E} \egroup\end{displaymath}

but that is a pseudo-scalar, not a scalar. That is, it changes sign under a parity transformation. The EM interaction is known to conserve parity so this is not a real option. As with the scalar field, we need to add an interaction with a source term. Of course, we know electromagnetism well, so finding the right Lagrangian is not really guess work. The source of the field is the vector \bgroup\color{black}$j_\mu$\egroup, so the simple scalar we can write is \bgroup\color{black}$j_\mu A_\mu$\egroup.

The Lagrangian for Classical Electricity and Magnetism we will try is.

\bgroup\color{black}$\displaystyle {\cal L_{EM}}=-{1\over 4}F_{\mu\nu}F_{\mu\nu}+{1\over c}j_\mu A_\mu $\egroup
In working with this Lagrangian, we will treat each component of $A$ as an independent field.

The next step is to check what the Euler-Lagrange equation gives us.

\begin{eqnarray*}
{\partial\over\partial x_\nu}\left({\partial{\cal L}\over\part...
...\\
{\partial\over\partial x_\nu}F_{\mu\nu}&=&{j_\mu\over c} \\
\end{eqnarray*}


Note that, since we have four independent components of \bgroup\color{black}$A_\mu$\egroup as independent fields, we have four equations; or one 4-vector equation. The Euler-Lagrange equation gets us back Maxwell's equation with this choice of the Lagrangian. This clearly justifies the choice of \bgroup\color{black}${\cal L}$\egroup.

It is important to emphasize that we have a Lagrangian based, formal classical field theory for electricity and magnetism which has the four components of the 4-vector potential as the independent fields.We could not treat each component of \bgroup\color{black}$F_{\mu\nu}$\egroup as independent since they are clearly correlated. We could have tried using the six independent components of the antisymmetric tensor but it would not have given the right answer. Using the 4-vector potentials as the fields does give the right answer. Electricity and Magnetism is a theory of a 4-vector field $A_\mu$.

We can also calculate the free field Hamiltonian density, that is, the Hamiltonian density in regions with no source term. We use the standard definition of the Hamiltonian in terms of the Lagrangian.

\begin{displaymath}\bgroup\color{black} {\cal H}=\left({\partial{\cal L}\over\pa...
...x_4)}\right){\partial A_\mu\over\partial x_4}
-{\cal L} \egroup\end{displaymath}

We just calculated above that

\begin{displaymath}\bgroup\color{black} {\partial{\cal L}\over\partial(\partial A_\mu/\partial x_\nu)}=F_{\mu\nu} \egroup\end{displaymath}

which we can use to get

\begin{eqnarray*}
{\partial{\cal L}\over\partial(\partial A_\mu/\partial x_4)}&=...
...rtial A_\mu\over\partial x_4}+{1\over 4}F_{\mu\nu}F_{\mu\nu} \\
\end{eqnarray*}


\bgroup\color{black}$\displaystyle {\cal H}=F_{\mu 4}{\partial A_\mu\over\partial x_4}+{1\over 4}F_{\mu\nu}F_{\mu\nu}$\egroup
We will use this once we have written the radiation field in a convenient form. In the meantime, we can check what this gives us in general in a region with no sources.

\begin{eqnarray*}
{\cal H}&=&F_{\mu 4}\left(F_{4\mu}+{\partial A_4\over\partial ...
...
&=&{1\over 2}(E^2+B^2) + E_i{\partial\phi\over\partial x_i} \\
\end{eqnarray*}


If we integrate the last term by parts, (and the fields fall to zero at infinity), then that term contains a \bgroup\color{black}$\vec{\nabla}\cdot\vec{E}$\egroup which is zero with no sources in the region. We can therefore drop it and are left with

\begin{displaymath}\bgroup\color{black} {\cal H}={1\over 2}(E^2+B^2). \egroup\end{displaymath}

This is the result we expected, the energy density and an EM field. (Remember the fields have been decreased by a factor of \bgroup\color{black}$\sqrt{4\pi}$\egroup compared to CGS units.)

We will study the interaction between electrons and the electromagnetic field with the Dirac equation. Until then, the Hamiltonian used for non-relativistic quantum mechanics will be sufficient. We have derived the Lorentz force law from that Hamiltonian.

\begin{displaymath}\bgroup\color{black} H={1\over 2m}\left(\vec{p}+{e\over c}\vec{A}\right)^2+eA_0 \egroup\end{displaymath}

Jim Branson 2013-04-22