Skip to article frontmatterSkip to article content
Site not loading correctly?

This may be due to an incorrect BASE_URL configuration. See the MyST Documentation for reference.

Relativistic electrodynamics

“If [the mathematicians of the nineteenth century] had taken Maxwell’s equations to heart as Euler took Newton’s, they would have discovered...Einstein’s theory of special relativity...simply by exploring to the end the mathematical concepts to which Maxwell’s equations naturally lead.” Freeman Dyson

Special relativity on its own may seem deeply at odds with our everyday experience of the world. Yet electricity and magnetism is certainly familiar to us; could you go without a day without seeing a battery, a circuit, or some other electronic device? Surprisingly, the laws of physics tell us that these two are not unrelated theories of physics. Indeed, special relativity is implied by electromagnetic theory! So in this chapter, we will explore why special relativity and electromagnetism are deeply interrelated, and what it means when we say that electromagnetism is a relativistic field theory.

The tensor formulation of Maxwell’s equations

It is not immediately apparent that Maxwell’s equations have anything to do with special relativity, but this is more to due with a choice of notation than anything else. Typically, we express Maxwell’s equations using vector calculus notation, as follows:

E=ρ/ε0B=0×E=Bt×B=μ0J+μ0ε0Et\begin{align} \nabla \cdot \mathbf{E} &= \rho/\varepsilon_0 \\ \nabla \cdot \mathbf{B} &= 0 \\ \nabla \times \mathbf{E} &= -\frac{\partial \mathbf{B}}{\partial t} \\ \nabla \times \mathbf{B} &= \mu_0 \mathbf{J} + \mu_0 \varepsilon_0 \frac{\partial \mathbf{E}}{\partial t} \\ \end{align}

This notation is all well and good for a variety of calculations, but it also obscures a lot of information; moreover, it is quite inelegant to have a system of four partial differential equations to express the dynamics of just two fields! When using tensor calculus to express Maxwell’s equations, the advantages of tensors becomes much more clear.

To start, it is simplest to work in terms of electromagnetic potentials instead of the fields themselves. Recall that the electric (scalar) potential ϕ\phi and the magnetic (vector) potential A\mathbf{A} are defined as follows:

E=ϕAtB=×A\begin{align} \mathbf{E} &= -\nabla \phi - \dfrac{\partial \mathbf{A}}{\partial t} \\ \mathbf{B} &= \nabla \times \mathbf{A} \end{align}

We may then define the electromagnetic 4-potential AμA^\mu and the 4-current JμJ^\mu as follows:

Aμ=(1cϕ,A)Jμ=(cρ,J)\begin{align} A^\mu &= \left(\frac{1}{c} \phi, \vec A\right) \\ J^\mu &= (c\rho, \vec J) \end{align}

Then, we define the Faraday tensor (also called the electromagnetic field tensor) as:

Fij=AjxiAixjF_{ij} = \frac{\partial A_j}{\partial x^i} - \frac{\partial A_i}{\partial x^j}

Using the Faraday tensor, we can then rewrite Maxwell’s equation using just 2 tensor equations, which, like any tensor equation, is coordinate-independent[1]:

mFik+kFmi+iFkm=0iFik=μ0Jk\begin{align} \partial_{m} F_{i k}+\partial_{k} F_{m i}+\partial_{i} F_{k m}=0 \\ \partial_{i} F^{i k}=\mu_{0} J^k \end{align}

If we are willing to dispense with the fields and only consider the potentials, we can write the Maxwell equations in an even simpler form as:

Aμ=μ0Jμ\square A^\mu = \mu_0 J^\mu

Where =μμ=ημννμ\square = \partial^\mu \partial_\mu = \eta^{\mu \nu} \partial_\nu \partial_\mu is known as the d’Alembert operator (also called the d’Alembertian), and is given by:

=1c22t22\square = \frac{1}{c^2} \frac{\partial^2}{\partial t^2} - \nabla^2

Relativistic solutions to the EM wave equation

The superposition principle tells us, among other things, that the general solution to an inhomogeneous differential equation in the form L^u(r)=f(r)\hat L u(\mathbf{r}) = f(\mathbf{r}) (where L^\hat L is a linear differential operator) may be obtained in the following manner:

  1. First, solve the homogeneous version of the differential equation by setting f(r)=0f(\mathbf{r}) = 0. This yields a solution u0u_0 that satisfies L^u0=0\hat L u_0 = 0

  2. Second, solve the inhomogeneous version of the differential equation for the given source function f(r)f(\mathbf{r}). This is given by L^up=f(r)\hat L u_p = f(\mathbf{r}), where upu_p is the particular solution.

    • Frequently, the source function is chosen to be a delta function, which physically represents a point source, since a superposition of delta functions can be used to describe any arbitrarily-complex source.

  3. Then, the general solution to the differential equation L^u=f(r)\hat L u = f(\mathbf{r}) for a particular source function f(r)f(\mathbf{r}) is given by u(r)=u0+upu(\mathbf{r}) = u_0 + u_p

Once the general solution is found to the inhomogeneous PDE, we can then solve specific initial/boundary-value problems by imposing the boundary condition(s) and initial condition(s) on the general solution.

Wave solutions

To start, let us consider the homogeneous Maxwell equations in 4-potential form, given by:

Aμ=0\square A^\mu = 0

This expands out to two partial differential equations for the scalar and vector potentials:

(1c22t22)ϕ=0(1c22t22)A=0\begin{align} \left(\frac{1}{c^2} \frac{\partial^2}{\partial t^2} - \nabla^2\right) \phi = 0 \\ \left(\frac{1}{c^2} \frac{\partial^2}{\partial t^2} - \nabla^2\right) \mathbf{A} = 0 \end{align}

Both of these PDEs are wave equations, meaning that their solutions are given by electromagnetic waves. The two most common wave solutions are plane waves and spherical waves. Both are idealized solutions that (for reasons we’ve already discussed) cannot exist in the real world.

It is also possible to write out a generalized solution in terms of a Fourier transform:

Liénard–Wiechert potentials

It can be shown that the inhomogeneous Maxwell equations have the particular solution (for arbitrary JμJ^\mu) of A(h)μ=(ϕp/c,Ap)A^\mu_{(h)} = (\phi_p/c, \mathbf{A}_p), where:

ϕp(r,t)=14πε0ρ(r,tr)rrd3rAp(r,t)=μ04πJ(r,tr)rrd3r\begin{align} \phi_p(\mathbf{r}, t) &= \frac{1}{4\pi \varepsilon_0} \int \frac{\rho(\mathbf{r}', t_r)}{|\mathbf{r} - \mathbf{r}'|} d^3 \mathbf{r}' \\ \mathbf{A}_p(\mathbf{r}, t) &= \frac{\mu_0}{4\pi} \int \frac{\mathbf{J}(\mathbf{r}', t_r)}{|\mathbf{r} - \mathbf{r}'|} d^3 \mathbf{r}' \end{align}

Or alternatively, in 4-vector notation:

A(h)μ=μ04πJμ(r,tr)rrd3rA^\mu_{(h)} = \frac{\mu_0}{4\pi}\int \frac{J^\mu(\mathbf{r}', t_r)}{|\mathbf{r} - \mathbf{r}'|} d^3 \mathbf{r}'

where trt_r is called the retarded time, given by:

tr=trrct_r = t - \frac{|\mathbf{r} - \mathbf{r}'|}{c}

In the case that we have a 4-current in the form of a delta function, corresponding to a point charge, we obtain the Liénard–Wiechert potentials.

Radiation by moving charges

Derive the Larmor formula for the radiated power of an accelerated charge, and its relativistic generalization.

The electromagnetic Lagrangian

L=14μ0FμνFμνAμJμ\mathscr{L} = -\frac{1}{4\mu_0} F_{\mu \nu} F^{\mu \nu} - A_\mu J^\mu
Footnotes
  1. A useful reference, from which this derivation was based, can be found at https://profoundphysics.com/are-maxwell-equations-relativistic/