Maxwell's equations in curved spacetime

Electromagnetic potential

The electromagnetic potential is a covariant vector A_α, which is the undefined primitive of electromagnetism. Being a covariant vector, its components transform from one coordinate system to another according to

${\displaystyle {\overline{A}}_{\beta }(\overline{x})=\frac{\mathrm{\partial }{x}^{\gamma }}{\mathrm{\partial }{\overline{x}}^{\beta }}{A}_{\gamma }(x).}$

Electromagnetic field

The electromagnetic field is a covariant antisymmetric tensor of degree 2, which can be defined in terms of the electromagnetic potential by $${\displaystyle F_{\alpha \beta }={\mathrm{\partial }}_{\alpha }{A}_{\beta }-{\mathrm{\partial }}_{\beta }{A}_{\alpha }.}$$

To see that this equation is invariant, we transform the coordinates as described in the classical treatment of tensors: $${\displaystyle \begin{array}{cc}{\overline{F}}_{\alpha \beta }& =\frac{\mathrm{\partial }{\overline{A}}_{\beta }}{\mathrm{\partial }{\overline{x}}^{\alpha }}-\frac{\mathrm{\partial }{\overline{A}}_{\alpha }}{\mathrm{\partial }{\overline{x}}^{\beta }}\\ & =\frac{\mathrm{\partial }}{\mathrm{\partial }{\overline{x}}^{\alpha }}\left(\frac{\mathrm{\partial }{x}^{\gamma }}{\mathrm{\partial }{\overline{x}}^{\beta }}{A}_{\gamma }\right)-\frac{\mathrm{\partial }}{\mathrm{\partial }{\overline{x}}^{\beta }}\left(\frac{\mathrm{\partial }{x}^{\delta }}{\mathrm{\partial }{\overline{x}}^{\alpha }}{A}_{\delta }\right)\\ & =\frac{{\mathrm{\partial }}^2{x}^{\gamma }}{\mathrm{\partial }{\overline{x}}^{\alpha }\mathrm{\partial }{\overline{x}}^{\beta }}{A}_{\gamma }+\frac{\mathrm{\partial }{x}^{\gamma }}{\mathrm{\partial }{\overline{x}}^{\beta }}\frac{\mathrm{\partial }{A}_{\gamma }}{\mathrm{\partial }{\overline{x}}^{\alpha }}-\frac{{\mathrm{\partial }}^2{x}^{\delta }}{\mathrm{\partial }{\overline{x}}^{\beta }\mathrm{\partial }{\overline{x}}^{\alpha }}{A}_{\delta }-\frac{\mathrm{\partial }{x}^{\delta }}{\mathrm{\partial }{\overline{x}}^{\alpha }}\frac{\mathrm{\partial }{A}_{\delta }}{\mathrm{\partial }{\overline{x}}^{\beta }}\\ & =\frac{\mathrm{\partial }{x}^{\gamma }}{\mathrm{\partial }{\overline{x}}^{\beta }}\frac{\mathrm{\partial }{x}^{\delta }}{\mathrm{\partial }{\overline{x}}^{\alpha }}\frac{\mathrm{\partial }{A}_{\gamma }}{\mathrm{\partial }{x}^{\delta }}-\frac{\mathrm{\partial }{x}^{\delta }}{\mathrm{\partial }{\overline{x}}^{\alpha }}\frac{\mathrm{\partial }{x}^{\gamma }}{\mathrm{\partial }{\overline{x}}^{\beta }}\frac{\mathrm{\partial }{A}_{\delta }}{\mathrm{\partial }{x}^{\gamma }}\\ & =\frac{\mathrm{\partial }{x}^{\delta }}{\mathrm{\partial }{\overline{x}}^{\alpha }}\frac{\mathrm{\partial }{x}^{\gamma }}{\mathrm{\partial }{\overline{x}}^{\beta }}\left(\frac{\mathrm{\partial }{A}_{\gamma }}{\mathrm{\partial }{x}^{\delta }}-\frac{\mathrm{\partial }{A}_{\delta }}{\mathrm{\partial }{x}^{\gamma }}\right)\\ & =\frac{\mathrm{\partial }{x}^{\delta }}{\mathrm{\partial }{\overline{x}}^{\alpha }}\frac{\mathrm{\partial }{x}^{\gamma }}{\mathrm{\partial }{\overline{x}}^{\beta }}{F}_{\delta \gamma }.\end{array}}$$

This definition implies that the electromagnetic field satisfies $${\displaystyle {\mathrm{\partial }}_{\lambda }{F}_{\mu \nu }+{\mathrm{\partial }}_{\mu }{F}_{\nu \lambda }+{\mathrm{\partial }}_{\nu }{F}_{\lambda \mu }=0,}$$ which incorporates Faraday's law of induction and Gauss's law for magnetism. This is seen from $${\displaystyle {\mathrm{\partial }}_{\lambda }{F}_{\mu \nu }+{\mathrm{\partial }}_{\mu }{F}_{\nu \lambda }+{\mathrm{\partial }}_{\nu }{F}_{\lambda \mu }={\mathrm{\partial }}_{\lambda }{\mathrm{\partial }}_{\mu }{A}_{\nu }-{\mathrm{\partial }}_{\lambda }{\mathrm{\partial }}_{\nu }{A}_{\mu }+{\mathrm{\partial }}_{\mu }{\mathrm{\partial }}_{\nu }{A}_{\lambda }-{\mathrm{\partial }}_{\mu }{\mathrm{\partial }}_{\lambda }{A}_{\nu }+{\mathrm{\partial }}_{\nu }{\mathrm{\partial }}_{\lambda }{A}_{\mu }-{\mathrm{\partial }}_{\nu }{\mathrm{\partial }}_{\mu }{A}_{\lambda }=0.}$$

Thus, the right-hand side of that Maxwell law is zero identically, meaning that the classic EM field theory leaves no room for magnetic monopoles or currents of such to act as sources of the field.

Although there appear to be 64 equations in Faraday–Gauss, it actually reduces to just four independent equations. Using the antisymmetry of the electromagnetic field, one can either reduce to an identity (0 = 0) or render redundant all the equations except for those with {λ, μ, ν} being either {1, 2, 3}, {2, 3, 0}, {3, 0, 1}, or {0, 1, 2}.

The Faraday–Gauss equation is sometimes written $${\displaystyle F_{[\mu \nu ;\lambda ]}=F_{[\mu \nu ,\lambda ]}=\frac{1}{6}({\mathrm{\partial }}_{\lambda }{F}_{\mu \nu }+{\mathrm{\partial }}_{\mu }{F}_{\nu \lambda }+{\mathrm{\partial }}_{\nu }{F}_{\lambda \mu }-{\mathrm{\partial }}_{\lambda }{F}_{\nu \mu }-{\mathrm{\partial }}_{\mu }{F}_{\lambda \nu }-{\mathrm{\partial }}_{\nu }{F}_{\mu \lambda })=\frac{1}{3}({\mathrm{\partial }}_{\lambda }{F}_{\mu \nu }+{\mathrm{\partial }}_{\mu }{F}_{\nu \lambda }+{\mathrm{\partial }}_{\nu }{F}_{\lambda \mu })=0,}$$ where a semicolon indicates a covariant derivative, a comma indicates a partial derivative, and square brackets indicate anti-symmetrization (see Ricci calculus for the notation). The covariant derivative of the electromagnetic field is $${\displaystyle F_{\alpha \beta ;\gamma }=F_{\alpha \beta ,\gamma }-{{\mathrm{\Gamma }}^{\mu }}_{\alpha \gamma }{F}_{\mu \beta }-{{\mathrm{\Gamma }}^{\mu }}_{\beta \gamma }{F}_{\alpha \mu },}$$ where Γ^α_βγ is the Christoffel symbol, which is symmetric in its lower indices.

Electromagnetic displacement

The electric displacement field D and the auxiliary magnetic field H form an antisymmetric contravariant rank-2 tensor density of weight +1. In vacuum, this is given by

${\displaystyle {\mathcal{D}}^{\mu \nu }=\frac{1}{{\mu }_0}{g}^{\mu \alpha }{F}_{\alpha \beta }{g}^{\beta \nu }\frac{\sqrt{-g}}{c}.}$

This equation is the only place where the metric (and thus gravity) enters into the theory of electromagnetism. Furthermore, the equation is invariant under a change of scale, that is, multiplying the metric by a constant has no effect on this equation. Consequently, gravity can only affect electromagnetism by changing the speed of light relative to the global coordinate system being used. Light is only deflected by gravity because it is slower near massive bodies. So it is as if gravity increased the index of refraction of space near massive bodies.

More generally, in materials where the magnetization–polarization tensor is non-zero, we have

${\displaystyle {\mathcal{D}}^{\mu \nu }=\frac{1}{{\mu }_0}{g}^{\mu \alpha }{F}_{\alpha \beta }{g}^{\beta \nu }\frac{\sqrt{-g}}{c}-{\mathcal{M}}^{\mu \nu }.}$

The transformation law for electromagnetic displacement is

${\displaystyle {\overline{\mathcal{D}}}^{\mu \nu }=\frac{\mathrm{\partial }{\overline{x}}^{\mu }}{\mathrm{\partial }{x}^{\alpha }}\frac{\mathrm{\partial }{\overline{x}}^{\nu }}{\mathrm{\partial }{x}^{\beta }}{\mathcal{D}}^{\alpha \beta }det\left[\frac{\mathrm{\partial }{x}^{\sigma }}{\mathrm{\partial }{\overline{x}}^{\rho }}\right],}$

where the Jacobian determinant is used. If the magnetization-polarization tensor is used, it has the same transformation law as the electromagnetic displacement.

Electric current

The electric current is the divergence of the electromagnetic displacement. In vacuum,

${\displaystyle J^{\mu }={\mathrm{\partial }}_{\nu }{\mathcal{D}}^{\mu \nu }.}$

If magnetization–polarization is used, then this just gives the free portion of the current

${\displaystyle J_{\text{free}}^{\mu }={\mathrm{\partial }}_{\nu }{\mathcal{D}}^{\mu \nu }.}$

This incorporates Ampere's law and Gauss's law.

In either case, the fact that the electromagnetic displacement is antisymmetric implies that the electric current is automatically conserved:

${\displaystyle {\mathrm{\partial }}_{\mu }{J}^{\mu }={\mathrm{\partial }}_{\mu }{\mathrm{\partial }}_{\nu }{\mathcal{D}}^{\mu \nu }=0,}$

because the partial derivatives commute.

The Ampere–Gauss definition of the electric current is not sufficient to determine its value because the electromagnetic potential (from which it was ultimately derived) has not been given a value. Instead, the usual procedure is to equate the electric current to some expression in terms of other fields, mainly the electron and proton, and then solve for the electromagnetic displacement, electromagnetic field, and electromagnetic potential.

The electric current is a contravariant vector density, and as such it transforms as follows:

${\displaystyle {\overline{J}}^{\mu }=\frac{\mathrm{\partial }{\overline{x}}^{\mu }}{\mathrm{\partial }{x}^{\alpha }}{J}^{\alpha }det\left[\frac{\mathrm{\partial }{x}^{\sigma }}{\mathrm{\partial }{\overline{x}}^{\rho }}\right].}$

Verification of this transformation law: $${\displaystyle \begin{array}{cc}{\overline{J}}^{\mu }& =\frac{\mathrm{\partial }}{\mathrm{\partial }{\overline{x}}^{\nu }}\left({\overline{\mathcal{D}}}^{\mu \nu }\right)\\ & =\frac{\mathrm{\partial }}{\mathrm{\partial }{\overline{x}}^{\nu }}\left(\frac{\mathrm{\partial }{\overline{x}}^{\mu }}{\mathrm{\partial }{x}^{\alpha }}\frac{\mathrm{\partial }{\overline{x}}^{\nu }}{\mathrm{\partial }{x}^{\beta }}{\mathcal{D}}^{\alpha \beta }det\left[\frac{\mathrm{\partial }{x}^{\sigma }}{\mathrm{\partial }{\overline{x}}^{\rho }}\right]\right)\\ & =\frac{{\mathrm{\partial }}^2{\overline{x}}^{\mu }}{\mathrm{\partial }{\overline{x}}^{\nu }\mathrm{\partial }{x}^{\alpha }}\frac{\mathrm{\partial }{\overline{x}}^{\nu }}{\mathrm{\partial }{x}^{\beta }}{\mathcal{D}}^{\alpha \beta }det\left[\frac{\mathrm{\partial }{x}^{\sigma }}{\mathrm{\partial }{\overline{x}}^{\rho }}\right]+\frac{\mathrm{\partial }{\overline{x}}^{\mu }}{\mathrm{\partial }{x}^{\alpha }}\frac{{\mathrm{\partial }}^2{\overline{x}}^{\nu }}{\mathrm{\partial }{\overline{x}}^{\nu }\mathrm{\partial }{x}^{\beta }}{\mathcal{D}}^{\alpha \beta }det\left[\frac{\mathrm{\partial }{x}^{\sigma }}{\mathrm{\partial }{\overline{x}}^{\rho }}\right]+\frac{\mathrm{\partial }{\overline{x}}^{\mu }}{\mathrm{\partial }{x}^{\alpha }}\frac{\mathrm{\partial }{\overline{x}}^{\nu }}{\mathrm{\partial }{x}^{\beta }}\frac{\mathrm{\partial }{\mathcal{D}}^{\alpha \beta }}{\mathrm{\partial }{\overline{x}}^{\nu }}det\left[\frac{\mathrm{\partial }{x}^{\sigma }}{\mathrm{\partial }{\overline{x}}^{\rho }}\right]+\frac{\mathrm{\partial }{\overline{x}}^{\mu }}{\mathrm{\partial }{x}^{\alpha }}\frac{\mathrm{\partial }{\overline{x}}^{\nu }}{\mathrm{\partial }{x}^{\beta }}{\mathcal{D}}^{\alpha \beta }\frac{\mathrm{\partial }}{\mathrm{\partial }{\overline{x}}^{\nu }}det\left[\frac{\mathrm{\partial }{x}^{\sigma }}{\mathrm{\partial }{\overline{x}}^{\rho }}\right]\\ & =\frac{{\mathrm{\partial }}^2{\overline{x}}^{\mu }}{\mathrm{\partial }{x}^{\beta }\mathrm{\partial }{x}^{\alpha }}{\mathcal{D}}^{\alpha \beta }det\left[\frac{\mathrm{\partial }{x}^{\sigma }}{\mathrm{\partial }{\overline{x}}^{\rho }}\right]+\frac{\mathrm{\partial }{\overline{x}}^{\mu }}{\mathrm{\partial }{x}^{\alpha }}\frac{{\mathrm{\partial }}^2{\overline{x}}^{\nu }}{\mathrm{\partial }{\overline{x}}^{\nu }\mathrm{\partial }{x}^{\beta }}{\mathcal{D}}^{\alpha \beta }det\left[\frac{\mathrm{\partial }{x}^{\sigma }}{\mathrm{\partial }{\overline{x}}^{\rho }}\right]+\frac{\mathrm{\partial }{\overline{x}}^{\mu }}{\mathrm{\partial }{x}^{\alpha }}\frac{\mathrm{\partial }{\mathcal{D}}^{\alpha \beta }}{\mathrm{\partial }{x}^{\beta }}det\left[\frac{\mathrm{\partial }{x}^{\sigma }}{\mathrm{\partial }{\overline{x}}^{\rho }}\right]+\frac{\mathrm{\partial }{\overline{x}}^{\mu }}{\mathrm{\partial }{x}^{\alpha }}\frac{\mathrm{\partial }{\overline{x}}^{\nu }}{\mathrm{\partial }{x}^{\beta }}{\mathcal{D}}^{\alpha \beta }det\left[\frac{\mathrm{\partial }{x}^{\sigma }}{\mathrm{\partial }{\overline{x}}^{\rho }}\right]\frac{\mathrm{\partial }{\overline{x}}^{\rho }}{\mathrm{\partial }{x}^{\sigma }}\frac{{\mathrm{\partial }}^2{x}^{\sigma }}{\mathrm{\partial }{\overline{x}}^{\nu }\mathrm{\partial }{\overline{x}}^{\rho }}\\ & =0+\frac{\mathrm{\partial }{\overline{x}}^{\mu }}{\mathrm{\partial }{x}^{\alpha }}\frac{{\mathrm{\partial }}^2{\overline{x}}^{\nu }}{\mathrm{\partial }{\overline{x}}^{\nu }\mathrm{\partial }{x}^{\beta }}{\mathcal{D}}^{\alpha \beta }det\left[\frac{\mathrm{\partial }{x}^{\sigma }}{\mathrm{\partial }{\overline{x}}^{\rho }}\right]+\frac{\mathrm{\partial }{\overline{x}}^{\mu }}{\mathrm{\partial }{x}^{\alpha }}{J}^{\alpha }det\left[\frac{\mathrm{\partial }{x}^{\sigma }}{\mathrm{\partial }{\overline{x}}^{\rho }}\right]+\frac{\mathrm{\partial }{\overline{x}}^{\mu }}{\mathrm{\partial }{x}^{\alpha }}{\mathcal{D}}^{\alpha \beta }det\left[\frac{\mathrm{\partial }{x}^{\sigma }}{\mathrm{\partial }{\overline{x}}^{\rho }}\right]\frac{\mathrm{\partial }{\overline{x}}^{\rho }}{\mathrm{\partial }{x}^{\sigma }}\frac{{\mathrm{\partial }}^2{x}^{\sigma }}{\mathrm{\partial }{x}^{\beta }\mathrm{\partial }{\overline{x}}^{\rho }}\\ & =\frac{\mathrm{\partial }{\overline{x}}^{\mu }}{\mathrm{\partial }{x}^{\alpha }}{J}^{\alpha }det\left[\frac{\mathrm{\partial }{x}^{\sigma }}{\mathrm{\partial }{\overline{x}}^{\rho }}\right]+\frac{\mathrm{\partial }{\overline{x}}^{\mu }}{\mathrm{\partial }{x}^{\alpha }}{\mathcal{D}}^{\alpha \beta }det\left[\frac{\mathrm{\partial }{x}^{\sigma }}{\mathrm{\partial }{\overline{x}}^{\rho }}\right]\left(\frac{{\mathrm{\partial }}^2{\overline{x}}^{\nu }}{\mathrm{\partial }{\overline{x}}^{\nu }\mathrm{\partial }{x}^{\beta }}+\frac{\mathrm{\partial }{\overline{x}}^{\rho }}{\mathrm{\partial }{x}^{\sigma }}\frac{{\mathrm{\partial }}^2{x}^{\sigma }}{\mathrm{\partial }{x}^{\beta }\mathrm{\partial }{\overline{x}}^{\rho }}\right).\end{array}}$$

So all that remains is to show that

${\displaystyle \frac{{\mathrm{\partial }}^2{\overline{x}}^{\nu }}{\mathrm{\partial }{\overline{x}}^{\nu }\mathrm{\partial }{x}^{\beta }}+\frac{\mathrm{\partial }{\overline{x}}^{\rho }}{\mathrm{\partial }{x}^{\sigma }}\frac{{\mathrm{\partial }}^2{x}^{\sigma }}{\mathrm{\partial }{x}^{\beta }\mathrm{\partial }{\overline{x}}^{\rho }}=0,}$

which is a version of a known theorem (see Inverse functions and differentiation § Higher derivatives). $${\displaystyle \begin{array}{cc}& \frac{{\mathrm{\partial }}^2{\overline{x}}^{\nu }}{\mathrm{\partial }{\overline{x}}^{\nu }\mathrm{\partial }{x}^{\beta }}+\frac{\mathrm{\partial }{\overline{x}}^{\rho }}{\mathrm{\partial }{x}^{\sigma }}\frac{{\mathrm{\partial }}^2{x}^{\sigma }}{\mathrm{\partial }{x}^{\beta }\mathrm{\partial }{\overline{x}}^{\rho }}\\ =& \frac{\mathrm{\partial }{x}^{\sigma }}{\mathrm{\partial }{\overline{x}}^{\nu }}\frac{{\mathrm{\partial }}^2{\overline{x}}^{\nu }}{\mathrm{\partial }{x}^{\sigma }\mathrm{\partial }{x}^{\beta }}+\frac{\mathrm{\partial }{\overline{x}}^{\nu }}{\mathrm{\partial }{x}^{\sigma }}\frac{{\mathrm{\partial }}^2{x}^{\sigma }}{\mathrm{\partial }{x}^{\beta }\mathrm{\partial }{\overline{x}}^{\nu }}\\ =& \frac{\mathrm{\partial }{x}^{\sigma }}{\mathrm{\partial }{\overline{x}}^{\nu }}\frac{{\mathrm{\partial }}^2{\overline{x}}^{\nu }}{\mathrm{\partial }{x}^{\beta }\mathrm{\partial }{x}^{\sigma }}+\frac{{\mathrm{\partial }}^2{x}^{\sigma }}{\mathrm{\partial }{x}^{\beta }\mathrm{\partial }{\overline{x}}^{\nu }}\frac{\mathrm{\partial }{\overline{x}}^{\nu }}{\mathrm{\partial }{x}^{\sigma }}\\ =& \frac{\mathrm{\partial }}{\mathrm{\partial }{x}^{\beta }}\left(\frac{\mathrm{\partial }{x}^{\sigma }}{\mathrm{\partial }{\overline{x}}^{\nu }}\frac{\mathrm{\partial }{\overline{x}}^{\nu }}{\mathrm{\partial }{x}^{\sigma }}\right)\\ =& \frac{\mathrm{\partial }}{\mathrm{\partial }{x}^{\beta }}\left(\frac{\mathrm{\partial }{\overline{x}}^{\nu }}{\mathrm{\partial }{\overline{x}}^{\nu }}\right)\\ =& \frac{\mathrm{\partial }}{\mathrm{\partial }{x}^{\beta }}\left(\mathbf{4}\right)\\ =& 0.\end{array}}$$

Lorentz force density

The density of the Lorentz force is a covariant vector density given by $${\displaystyle f_{\mu }=F_{\mu \nu }{J}^{\nu }.}$$

The force on a test particle subject only to gravity and electromagnetism is $${\displaystyle \frac{d{p}_{\alpha }}{dt}={\mathrm{\Gamma }}_{\alpha \gamma }^{\beta }{p}_{\beta }\frac{d{x}^{\gamma }}{dt}+q{F}_{\alpha \gamma }\frac{d{x}^{\gamma }}{dt},}$$ where p_α is the linear 4-momentum of the particle, t is any time coordinate parameterizing the world line of the particle, Γ^β_αγ is the Christoffel symbol (gravitational force field), and q is the electric charge of the particle.

This equation is invariant under a change in the time coordinate; just multiply by ${\displaystyle dt/d\overline{t}}$ and use the chain rule. It is also invariant under a change in the x coordinate system.

Using the transformation law for the Christoffel symbol, $${\displaystyle {\overline{\mathrm{\Gamma }}}_{\alpha \gamma }^{\beta }=\frac{\mathrm{\partial }{\overline{x}}^{\beta }}{\mathrm{\partial }{x}^{ϵ}}\frac{\mathrm{\partial }{x}^{\delta }}{\mathrm{\partial }{\overline{x}}^{\alpha }}\frac{\mathrm{\partial }{x}^{\zeta }}{\mathrm{\partial }{\overline{x}}^{\gamma }}{\mathrm{\Gamma }}_{\delta \zeta }^{ϵ}+\frac{\mathrm{\partial }{\overline{x}}^{\beta }}{\mathrm{\partial }{x}^{\eta }}\frac{{\mathrm{\partial }}^2{x}^{\eta }}{\mathrm{\partial }{\overline{x}}^{\alpha }\mathrm{\partial }{\overline{x}}^{\gamma }},}$$ we get $${\displaystyle \begin{array}{cc}& \frac{d{\overline{p}}_{\alpha }}{dt}-{\overline{\mathrm{\Gamma }}}_{\alpha \gamma }^{\beta }{\overline{p}}_{\beta }\frac{d{\overline{x}}^{\gamma }}{dt}-q{\overline{F}}_{\alpha \gamma }\frac{d{\overline{x}}^{\gamma }}{dt}\\ =& \frac{d}{dt}\left(\frac{\mathrm{\partial }{x}^{\delta }}{\mathrm{\partial }{\overline{x}}^{\alpha }}{p}_{\delta }\right)-\left(\frac{\mathrm{\partial }{\overline{x}}^{\beta }}{\mathrm{\partial }{x}^{\theta }}\frac{\mathrm{\partial }{x}^{\delta }}{\mathrm{\partial }{\overline{x}}^{\alpha }}\frac{\mathrm{\partial }{x}^{\iota }}{\mathrm{\partial }{\overline{x}}^{\gamma }}{\mathrm{\Gamma }}_{\delta \iota }^{\theta }+\frac{\mathrm{\partial }{\overline{x}}^{\beta }}{\mathrm{\partial }{x}^{\eta }}\frac{{\mathrm{\partial }}^2{x}^{\eta }}{\mathrm{\partial }{\overline{x}}^{\alpha }\mathrm{\partial }{\overline{x}}^{\gamma }}\right)\frac{\mathrm{\partial }{x}^{ϵ}}{\mathrm{\partial }{\overline{x}}^{\beta }}{p}_{ϵ}\frac{\mathrm{\partial }{\overline{x}}^{\gamma }}{\mathrm{\partial }{x}^{\zeta }}\frac{d{x}^{\zeta }}{dt}-q\frac{\mathrm{\partial }{x}^{\delta }}{\mathrm{\partial }{\overline{x}}^{\alpha }}{F}_{\delta \zeta }\frac{d{x}^{\zeta }}{dt}\\ =& \frac{\mathrm{\partial }{x}^{\delta }}{\mathrm{\partial }{\overline{x}}^{\alpha }}\left(\frac{d{p}_{\delta }}{dt}-{\mathrm{\Gamma }}_{\delta \zeta }^{ϵ}{p}_{ϵ}\frac{d{x}^{\zeta }}{dt}-q{F}_{\delta \zeta }\frac{d{x}^{\zeta }}{dt}\right)+\frac{d}{dt}\left(\frac{\mathrm{\partial }{x}^{\delta }}{\mathrm{\partial }{\overline{x}}^{\alpha }}\right)p_{\delta }-\left(\frac{\mathrm{\partial }{\overline{x}}^{\beta }}{\mathrm{\partial }{x}^{\eta }}\frac{{\mathrm{\partial }}^2{x}^{\eta }}{\mathrm{\partial }{\overline{x}}^{\alpha }\mathrm{\partial }{\overline{x}}^{\gamma }}\right)\frac{\mathrm{\partial }{x}^{ϵ}}{\mathrm{\partial }{\overline{x}}^{\beta }}{p}_{ϵ}\frac{\mathrm{\partial }{\overline{x}}^{\gamma }}{\mathrm{\partial }{x}^{\zeta }}\frac{d{x}^{\zeta }}{dt}\\ =& 0+\frac{d}{dt}\left(\frac{\mathrm{\partial }{x}^{\delta }}{\mathrm{\partial }{\overline{x}}^{\alpha }}\right)p_{\delta }-\frac{{\mathrm{\partial }}^2{x}^{ϵ}}{\mathrm{\partial }{\overline{x}}^{\alpha }\mathrm{\partial }{\overline{x}}^{\gamma }}{p}_{ϵ}\frac{d{\overline{x}}^{\gamma }}{dt}\\ =& 0.\end{array}}$$

Lagrangian

In vacuum, the Lagrangian density for classical electrodynamics (in joules per cubic meter) is a scalar density $${\displaystyle \mathcal{L}=-\frac{1}{4{\mu }_0}{F}_{\alpha \beta }{F}^{\alpha \beta }\frac{\sqrt{-g}}{c}+A_{\alpha }{J}^{\alpha },}$$ where $${\displaystyle F^{\alpha \beta }=g^{\alpha \gamma }{F}_{\gamma \delta }{g}^{\delta \beta }.}$$

The 4-current should be understood as an abbreviation of many terms expressing the electric currents of other charged fields in terms of their variables.

If we separate free currents from bound currents, the Lagrangian becomes $${\displaystyle \mathcal{L}=-\frac{1}{4{\mu }_0}{F}_{\alpha \beta }{F}^{\alpha \beta }\frac{\sqrt{-g}}{c}+A_{\alpha }{J}_{\text{free}}^{\alpha }+\frac{1}{2}{F}_{\alpha \beta }{\mathcal{M}}^{\alpha \beta }.}$$

Electromagnetic stress–energy tensor

As part of the source term in the Einstein field equations, the electromagnetic stress–energy tensor is a covariant symmetric tensor $${\displaystyle T_{\mu \nu }=-\frac{1}{{\mu }_0}\left(F_{\mu \alpha }{g}^{\alpha \beta }{F}_{\beta \nu }-\frac{1}{4}{g}_{\mu \nu }{F}_{\sigma \alpha }{g}^{\alpha \beta }{F}_{\beta \rho }{g}^{\rho \sigma }\right),}$$ using a metric of signature (−, +, +, +). If using the metric with signature (+, −, −, −), the expression for ${\displaystyle T_{\mu \nu }}$ will have opposite sign. The stress–energy tensor is trace-free: $${\displaystyle T_{\mu \nu }{g}^{\mu \nu }=0}$$ because electromagnetism propagates at the local invariant speed, and is conformal-invariant.

In the expression for the conservation of energy and linear momentum, the electromagnetic stress–energy tensor is best represented as a mixed tensor density $${\displaystyle {\mathfrak{T}}_{\mu }^{\nu }=T_{\mu \gamma }{g}^{\gamma \nu }\frac{\sqrt{-g}}{c}.}$$

From the equations above, one can show that $${\displaystyle {{\mathfrak{T}}_{\mu }^{\nu }}_{;\nu }+f_{\mu }=0,}$$ where the semicolon indicates a covariant derivative.

This can be rewritten as $${\displaystyle -{{\mathfrak{T}}_{\mu }^{\nu }}_{,\nu }=-{\mathrm{\Gamma }}_{\mu \nu }^{\sigma }{\mathfrak{T}}_{\sigma }^{\nu }+f_{\mu },}$$ which says that the decrease in the electromagnetic energy is the same as the work done by the electromagnetic field on the gravitational field plus the work done on matter (via the Lorentz force), and similarly the rate of decrease in the electromagnetic linear momentum is the electromagnetic force exerted on the gravitational field plus the Lorentz force exerted on matter.

Derivation of conservation law: $${\displaystyle \begin{array}{cc}{{\mathfrak{T}}_{\mu }^{\nu }}_{;\nu }+f_{\mu }& =-\frac{1}{{\mu }_0}\left(F_{\mu \alpha ;\nu }{g}^{\alpha \beta }{F}_{\beta \gamma }{g}^{\gamma \nu }+F_{\mu \alpha }{g}^{\alpha \beta }{F}_{\beta \gamma ;\nu }{g}^{\gamma \nu }-\frac{1}{2}{\delta }_{\mu }^{\nu }{F}_{\sigma \alpha ;\nu }{g}^{\alpha \beta }{F}_{\beta \rho }{g}^{\rho \sigma }\right)\frac{\sqrt{-g}}{c}+\frac{1}{{\mu }_0}{F}_{\mu \alpha }{g}^{\alpha \beta }{F}_{\beta \gamma ;\nu }{g}^{\gamma \nu }\frac{\sqrt{-g}}{c}\\ & =-\frac{1}{{\mu }_0}\left(F_{\mu \alpha ;\nu }{F}^{\alpha \nu }-\frac{1}{2}{F}_{\sigma \alpha ;\mu }{F}^{\alpha \sigma }\right)\frac{\sqrt{-g}}{c}\\ & =-\frac{1}{{\mu }_0}\left(\left(-F_{\nu \mu ;\alpha }-F_{\alpha \nu ;\mu }\right)F^{\alpha \nu }-\frac{1}{2}{F}_{\sigma \alpha ;\mu }{F}^{\alpha \sigma }\right)\frac{\sqrt{-g}}{c}\\ & =-\frac{1}{{\mu }_0}\left(F_{\mu \nu ;\alpha }{F}^{\alpha \nu }-F_{\alpha \nu ;\mu }{F}^{\alpha \nu }+\frac{1}{2}{F}_{\sigma \alpha ;\mu }{F}^{\sigma \alpha }\right)\frac{\sqrt{-g}}{c}\\ & =-\frac{1}{{\mu }_0}\left(F_{\mu \alpha ;\nu }{F}^{\nu \alpha }-\frac{1}{2}{F}_{\alpha \nu ;\mu }{F}^{\alpha \nu }\right)\frac{\sqrt{-g}}{c}\\ & =-\frac{1}{{\mu }_0}\left(-F_{\mu \alpha ;\nu }{F}^{\alpha \nu }+\frac{1}{2}{F}_{\sigma \alpha ;\mu }{F}^{\alpha \sigma }\right)\frac{\sqrt{-g}}{c},\end{array}}$$

which is zero because it is the negative of itself (see four lines above).

Electromagnetic wave equation

The nonhomogeneous electromagnetic wave equation in terms of the field tensor is modified from the special-relativity form to

${\displaystyle ◻F_{ab}\stackrel{\text{def}}{=}{F}_{ab;}{}^d{}_d=-2R_{acbd}{F}^{cd}+R_{ae}{F}^e{}_b-R_{be}{F}^e{}_a+J_{a;b}-J_{b;a},}$

where R_acbd is the covariant form of the Riemann tensor, and ${\displaystyle ◻}$ is a generalization of the d'Alembertian operator for covariant derivatives. Using

${\displaystyle ◻A^a={{A^{a;}}^b}_b.}$

Maxwell's source equations can be written in terms of the 4-potential [ref. 2, p. 569] as

${\displaystyle ◻A^a-{A^{b;a}}_b=-{\mu }_0{J}^a}$

or, assuming the generalization of the Lorenz gauge in curved spacetime,

${\displaystyle \begin{array}{cc}{A^a}_{;a}& =0,\\ ◻A^a& =-{\mu }_0{J}^a+{R^a}_b{A}^b,\end{array}}$

where ${\displaystyle R_{ab}\stackrel{\text{def}}{=}{R^s}_{asb}}$ is the Ricci curvature tensor.

This is the same form of the wave equation as in flat spacetime, except that the derivatives are replaced by covariant derivatives and there is an additional term proportional to the curvature. The wave equation in this form also bears some resemblance to the Lorentz force in curved spacetime, where A^a plays the role of the 4-position.

For the case of a metric signature in the form (+, −, −, −), the derivation of the wave equation in curved spacetime is carried out in the article.

Nonlinearity of Maxwell's equations in a dynamic spacetime

When Maxwell's equations are treated in a background-independent manner, that is, when the spacetime metric is taken to be a dynamical variable dependent on the electromagnetic field, then the electromagnetic wave equation and Maxwell's equations are nonlinear. This can be seen by noting that the curvature tensor depends on the stress–energy tensor through the Einstein field equation

${\displaystyle G_{ab}=\frac{8\pi G}{c^4}{T}_{ab},}$

where

${\displaystyle G_{ab}\stackrel{\text{def}}{=}{R}_{ab}-\frac{1}{2}R{g}_{ab}}$

is the Einstein tensor, G is the Newtonian constant of gravitation, g_ab is the metric tensor, and R (scalar curvature) is the trace of the Ricci curvature tensor. The stress–energy tensor is composed of the stress–energy from particles, but also stress–energy from the electromagnetic field. This generates the nonlinearity.

Geometric formulation

In the differential geometric formulation of the electromagnetic field, the antisymmetric Faraday tensor can be considered as the Faraday 2-form ${\displaystyle \mathbf{F}}$. In this view, one of Maxwell's two equations is

${\displaystyle \mathrm{d}\mathbf{F}=0,}$

where ${\displaystyle \mathrm{d}}$ is the exterior derivative operator. This equation is completely coordinate- and metric-independent and says that the electromagnetic flux through a closed two-dimensional surface in space–time is topological, more precisely, depends only on its homology class (a generalization of the integral form of Gauss law and Maxwell–Faraday equation, as the homology class in Minkowski space is automatically 0). By the Poincaré lemma, this equation implies (at least locally) that there exists a 1-form ${\displaystyle \mathbf{A}}$ satisfying

${\displaystyle \mathbf{F}=\mathrm{d}\mathbf{A}.}$

The other equation is

${\displaystyle \mathrm{d}\star \mathbf{F}=\mathbf{J}.}$

In this context, ${\displaystyle \mathbf{J}}$ is the current 3-form (or even more precise, twisted 3-form), and the star ${\displaystyle \star }$ denotes the Hodge star operator. The dependence of Maxwell's equation on the metric of spacetime lies in the Hodge star operator ${\displaystyle \star }$ on 2-forms, which is conformally invariant. Written this way, Maxwell's equation is the same in any space–time, manifestly coordinate-invariant, and convenient to use (even in Minkowski space or Euclidean space and time, especially with curvilinear coordinates).

An alternative geometric interpretation is that the Faraday 2-form ${\displaystyle \mathbf{F}}$ is (up to a factor ${\displaystyle i}$) the curvature 2-form ${\displaystyle F(\mathrm{\nabla })}$ of a U(1)-connection ${\displaystyle \mathrm{\nabla }}$ on a principal U(1)-bundle whose sections represent charged fields. The connection is much like the vector potential, since every connection can be written as ${\displaystyle \mathrm{\nabla }={\mathrm{\nabla }}_0+iA}$ for a "base" connection ${\displaystyle {\mathrm{\nabla }}_0}$, and

${\displaystyle \mathbf{F}={\mathbf{F}}_0+\mathrm{d}\mathbf{A}.}$

In this view, the Maxwell "equation" ${\displaystyle \mathrm{d}\mathbf{F}=0}$ is a mathematical identity known as the Bianchi identity. The equation ${\displaystyle \mathrm{d}\star \mathbf{F}=\mathbf{J}}$ is the only equation with any physical content in this formulation. This point of view is particularly natural when considering charged fields or quantum mechanics. It can be interpreted as saying that, much like gravity can be understood as being the result of the necessity of a connection to parallel transport vectors at different points, electromagnetic phenomena, or more subtle quantum effects like the Aharonov–Bohm effect, can be understood as a result from the necessity of a connection to parallel transport charged fields or wave sections at different points. In fact, just as the Riemann tensor is the holonomy of the Levi-Civita connection along an infinitesimal closed curve, the curvature of the connection is the holonomy of the U(1)-connection.