Wave Equation for the Four-Potential

Maxwell wrote down his field equations in 1865, and buried in them was a number he had not put there. Combining the laws of electricity and magnetism produced a wave, and the wave's speed was 1/μ0ε01/\sqrt{\mu_0\varepsilon_0} — two constants measured with coils and capacitors, no optics anywhere in sight. It came out equal to the speed of light. Written in spacetime, that whole story compresses to one line:

□Aμ=μ0Jμ.\Box A^\mu = \mu_0 J^\mu .

Four coupled vector equations for E⃗\vec E and B⃗\vec B become one wave equation for one four-vector, driven by one four-current. And the operator on the left is no new invention: □=∂μ∂μ=1c2∂t2−∇2\Box = \partial_\mu\partial^\mu = \frac{1}{c^2}\partial_t^2 - \nabla^2 is the Minkowski metric written as a differential operator. The relative minus sign between the time and space terms — the only thing separating a wave equation from Laplace's — is the signature (+,−,−,−)(+,-,-,-) itself. So when □Aμ=0\Box A^\mu = 0 forces its plane waves onto the light cone, kμkμ=0k_\mu k^\mu = 0 and hence ω=c∣k⃗∣\omega = c|\vec k|, that speed is not an input. It is the metric, read back out. A null four-vector stays null under every Lorentz transformation, so every inertial observer measures the same wave speed — Einstein's second postulate, sitting inside Maxwell's equations forty years before anyone thought to postulate it.

Key takeaways
  • Signature (+,−,−,−)(+,-,-,-), p⋅p=+m2c2p\cdot p = +m^2c^2; mm is the invariant rest mass. Dictionary: Aμ=(ϕ/c,A⃗)A^\mu = (\phi/c,\vec A), Jμ=(cρ,J⃗)J^\mu = (c\rho,\vec J), ∂μ=(1c∂t,∇)\partial_\mu = (\tfrac1c\partial_t, \nabla).
  • The equation: □Aμ=μ0Jμ\boxed{\Box A^\mu = \mu_0 J^\mu}, in the Lorenz gauge ∂μAμ=0\partial_\mu A^\mu = 0. Four uncoupled wave equations, each driven by the matching component of the four-current.
  • The d'Alembertian is the metric. □≡∂μ∂μ=ημν∂μ∂ν=1c2∂t2−∇2\Box \equiv \partial_\mu\partial^\mu = \eta^{\mu\nu}\partial_\mu\partial_\nu = \frac{1}{c^2}\partial_t^2 - \nabla^2. Flip to a Euclidean signature and the minus becomes a plus: the wave equation collapses into Laplace's — no light, no retardation, no causality. Waves exist because the metric is not positive-definite.
  • The gauge is a scalar condition, and that is why this works. ∂μAμ\partial_\mu A^\mu is a four-divergence: vanish in one frame, vanish in all.
  • Counting the photon. AμA^\mu has 4 components; the Lorenz condition removes 1; the residual gauge freedom (□λ=0\Box\lambda = 0) removes 1 more. 4−1−1=24 - 1 - 1 = \mathbf{2} — the two polarisations of light, obtained by bookkeeping.
  • Key idea: cc is derived, not assumed. □Aμ=0⇒kμkμ=0⇒ω=c∣k⃗∣\Box A^\mu = 0 \Rightarrow k_\mu k^\mu = 0 \Rightarrow \omega = c|\vec k|, and kμkμk_\mu k^\mu is a Lorentz scalar, so the same cc emerges in every frame. Numerically 1/μ0ε0=2.99792458×108 m s−11/\sqrt{\mu_0\varepsilon_0} = \mathbf{2.99792458\times10^{8}\ m\,s^{-1}} — a capacitor constant and a coil constant, with no optics anywhere in it.
The core insight

The wave operator and the metric are the same object. □=ημν∂μ∂ν\Box = \eta^{\mu\nu}\partial_\mu\partial_\nu — feed two derivatives to the metric and you get the d'Alembertian, exactly as feeding it two displacements gives the interval. Everything follows: the minus sign that makes waves instead of harmonic functions, the light cone that makes causality, the null condition that makes cc frame-independent. Maxwell's equations were never "made" relativistic. They were relativistic — and it took forty years to notice that the wave operator physics had used since 1747 already had Minkowski spacetime inside it.

Where this lives today

Every antenna, fibre, radar and radio link is a boundary-value problem for □Aμ=μ0Jμ\Box A^\mu = \mu_0 J^\mu, and the engineer's version — retarded potentials, near and far zone, the 1/R1/R radiation term — is the three-dimensional shadow of this one equation. Finite-difference time-domain solvers step the hyperbolic operator forward on a grid, and the reason they can step forward at all, with a finite domain of dependence and a stability condition tied to c Δt/Δxc\,\Delta t/\Delta x, is the light cone that the minus sign in □\Box puts there. The equation's deeper career is as a template. Replace AμA^\mu by the metric perturbation hˉμν\bar h^{\mu\nu} and the harmonic gauge plays exactly the role of the Lorenz gauge: the linearised Einstein equations become □hˉμν∝Tμν\Box \bar h^{\mu\nu} \propto T^{\mu\nu} — same operator, same gauge logic, two propagating polarisations after the same subtraction, and gravitational waves travelling at cc for the same reason light does. When LIGO timed gravitational waves against a gamma-ray burst in 2017 and found their speed equal to cc to a part in 101510^{15}, what was tested was that the same □\Box — the same signature — governs both fields. In quantum field theory this page's 4−1−1=24-1-1=2 becomes the photon's two helicities.

The physical picture

Start with the covariant field equation, ∂μFμν=μ0Jν\partial_\mu F^{\mu\nu} = \mu_0 J^\nu — owned by `aq-sr-covariant-maxwell-equations`. Substituting Fμν=∂μAν−∂νAμF^{\mu\nu} = \partial^\mu A^\nu - \partial^\nu A^\mu turns it into an equation for the potential, and something awkward appears:

□Aν−∂ν(∂μAμ)=μ0Jν.\Box A^\nu - \partial^\nu\big(\partial_\mu A^\mu\big) = \mu_0 J^\nu .

The first term is what we want: a clean wave operator acting on one component at a time. The second is a menace — it couples all four components of AμA^\mu through their common divergence. But AμA^\mu is not an observable. The fields FμνF^{\mu\nu} are unchanged by Aμ→Aμ−∂μλA^\mu \to A^\mu - \partial^\mu\lambda for any scalar λ\lambda, so we may choose AμA^\mu to make our lives easy. Choose ∂μAμ=0\partial_\mu A^\mu = 0 — the Lorenz gauge, owned by `aq-sr-gauge-transformations-lorenz-gauge`. The menacing term vanishes identically, and four uncoupled wave equations remain. Two things make this more than a trick. First, ∂μAμ=0\partial_\mu A^\mu = 0 says a Lorentz scalar vanishes — zero in one frame, zero in every frame — so the gauge can be imposed globally and covariantly. (Contrast the Coulomb gauge ∇⋅A⃗=0\nabla\cdot\vec A = 0: enormously useful, but it picks out a time axis and is not preserved by boosts.) Second, once the equation reads □Aμ=μ0Jμ\Box A^\mu = \mu_0 J^\mu, every symbol is a tensor — □\Box a scalar operator, AμA^\mu and JμJ^\mu four-vectors. It has the same form in every inertial frame by inspection, which is what "manifestly covariant" means, and it is why this equation, not the four three-dimensional ones, is the honest statement of electromagnetism.

Figure
THE d'ALEMBERTIAN IS THE METRIC WEARING DERIVATIVES. Feed the metric two displacements and you get the interval, ds² = c²dt² − dx² − dy² − dz². Feed it two derivatives and you get the wave operator, □ = (1/c²)∂ₜ² − ∇². It is the SAME contraction with the SAME η^{μν} — the minus signs come from exactly one place, the signature. And the consequence is total: with the relative minus sign the equation is HYPERBOLIC, its solutions are waves travelling at c, and it has a light cone. Make the signature Euclidean (+,+,+,+) and the same operator becomes the four-dimensional LAPLACIAN — elliptic, no waves, no propagation speed, no causality, no light. The single minus sign in the metric is the entire difference between a universe with light in it and a universe without.
η^{μν} — one metric, two famous objects η^{μν} = diag(+1,−1,−1,−1) the signature. that is all. feed it two displacements feed it two derivatives THE INTERVAL ds² = η_{μν} dx^μ dx^ν = c²dt² − dx² − dy² − dz² the geometry of spacetime THE d'ALEMBERTIAN □ = η^{μν} ∂_μ ∂_ν = (1/c²)∂ₜ² − ∇² the dynamics of the field KEEP THE MINUS (+,−,−,−) □A^μ = 0 is HYPERBOLIC solutions: waves, speed c light cone · causality · LIGHT k_μ k^μ = 0 ⇒ ω = c|k|, in every frame FLIP IT (+,+,+,+) ∇₄²A^μ = 0 is ELLIPTIC solutions: harmonic functions no propagation. no cone. no radiation, no retardation, NO LIGHT

The core equations

The equation this page owns:

  □Aμ  =  μ0Jμ  in the Lorenz gauge  ∂μAμ=0(1)\boxed{\;\Box A^\mu \;=\; \mu_0 J^\mu\;}\qquad\text{in the Lorenz gauge }\ \partial_\mu A^\mu = 0 \qquad (1)

where: Aμ=(ϕ/c,A⃗)A^\mu = (\phi/c, \vec A) is the four-potential; Jμ=(cρ,J⃗)J^\mu = (c\rho, \vec J) the four-current. Physical reading: four uncoupled wave equations. The μ=0\mu = 0 component reads □ϕ=ρ/ε0\Box\phi = \rho/\varepsilon_0 (using μ0c2=1/ε0\mu_0 c^2 = 1/\varepsilon_0), the spatial ones □A⃗=μ0J⃗\Box\vec A = \mu_0\vec J. Charge density sources the scalar potential; current density sources the vector potential; nothing else couples. The operator, spelled out:

□  ≡  ∂μ∂μ  =  ημν∂μ∂ν  =  1c2∂2∂t2  −  ∇2(2)\Box \;\equiv\; \partial_\mu\partial^\mu \;=\; \eta^{\mu\nu}\partial_\mu\partial_\nu \;=\; \frac{1}{c^2}\frac{\partial^2}{\partial t^2} \;-\; \nabla^2 \qquad (2)

where: ∂μ=(1c∂t, ∇)\partial_\mu = \big(\tfrac1c\partial_t,\ \nabla\big) and ∂μ=(1c∂t, −∇)\partial^\mu = \big(\tfrac1c\partial_t,\ -\nabla\big). Physical reading: raising the index is where the minus sign is born, and there is no other source of it. Contract two four-gradients with diag(+1,−1,−1,−1)\mathrm{diag}(+1,-1,-1,-1) and you get one plus and three minuses — hyperbolic. Contract with diag(+1,+1,+1,+1)\mathrm{diag}(+1,+1,+1,+1) and you get the four-dimensional Laplacian: elliptic, with no waves in its solution set at all. The vacuum limit, and the dispersion relation it forces:

□Aμ=0,Aμ=εμe−ikνxν⟹  kμkμ=ω2c2−∣k⃗∣2=0  ⟹  ω=c ∣k⃗∣  (3)\Box A^\mu = 0,\quad A^\mu = \varepsilon^\mu e^{-i k_\nu x^\nu} \qquad\Longrightarrow\qquad \boxed{\;k_\mu k^\mu = \frac{\omega^2}{c^2} - |\vec k|^2 = 0\;}\qquad\Longrightarrow\qquad \boxed{\;\omega = c\,|\vec k|\;} \qquad (3)

where: kμ=(ω/c, k⃗)k^\mu = (\omega/c,\ \vec k) is the wave four-vector and εμ\varepsilon^\mu the polarisation four-vector. Physical reading: the wave four-vector is null — light lives on the light cone, not inside it. And because kμkμk_\mu k^\mu is a Lorentz scalar, a null kk is null in every frame: every inertial observer finds ω′=c∣k⃗′∣\omega' = c|\vec k'| with the same cc. The frame-independence of the speed of light is not an assumption bolted onto Maxwell; it is a corollary of the operator in (2). The gauge condition, and what it does not do:

∂μAμ=1c2∂ϕ∂t+∇⋅A⃗=0,with residual freedomAμ→Aμ−∂μλ  for any  □λ=0(4)\partial_\mu A^\mu = \frac{1}{c^2}\frac{\partial\phi}{\partial t} + \nabla\cdot\vec A = 0, \qquad\text{with residual freedom}\qquad A^\mu \to A^\mu - \partial^\mu\lambda\ \ \text{for any}\ \ \Box\lambda = 0 \qquad (4)

where: λ\lambda is any scalar field obeying the homogeneous wave equation. Physical reading: the Lorenz gauge does not fix the potential — only its four-divergence. A residual gauge function survives, and it is exactly enough to remove one more component of AμA^\mu. That is where the photon's polarisation count comes from.

Reducing Maxwell to a wave equation

Derivation AFrom ∂_μF^{μν} = μ₀J^ν to □A^ν = μ₀J^ν — and why the gauge choice must be a scalar one

The goal is to take the covariant inhomogeneous Maxwell equation, substitute the potential, and watch four coupled equations fall apart into four independent ones. The manoeuvre turns on a single term, ∂ν(∂μAμ)\partial^\nu(\partial_\mu A^\mu), and on our being allowed to make it vanish. What deserves attention is not that we can set ∂μAμ=0\partial_\mu A^\mu = 0, but that this condition is a Lorentz scalar equation — and that a non-scalar gauge, however convenient, would wreck the covariance we are trying to exhibit.

Step 1. Insert the potential. With Fμν=∂μAν−∂νAμF^{\mu\nu} = \partial^\mu A^\nu - \partial^\nu A^\mu, \partial_\mu F^{\mu\nu} = \partial_\mu\partial^\mu A^\nu - \partial_\mu\partial^\nu A^\mu = \mu_0 J^\nu . Why: nothing is assumed yet — an identity, true in any gauge. The construction of FμνF^{\mu\nu} from AμA^\mu belongs to `aq-sr-electromagnetic-field-tensor`; we use it, we do not rebuild it.

Step 2. Factor the second term. Partials commute, so ∂μ∂νAμ=∂ν(∂μAμ)\partial_\mu\partial^\nu A^\mu = \partial^\nu(\partial_\mu A^\mu): \Box A^\nu - \partial^\nu\big(\partial_\mu A^\mu\big) = \mu_0 J^\nu . Why: <strong>this is the coupling term, and the only obstacle.</strong> It ties the evolution of A0A^0 to the divergence of A⃗\vec A and back. Every component of AμA^\mu sits inside it, which is exactly what stops the four equations separating.

Step 3. Impose the Lorenz gauge. Gauge freedom shifts Aμ→Aμ−∂μλA^\mu \to A^\mu - \partial^\mu\lambda without touching FμνF^{\mu\nu}, and shifts the divergence by ∂μAμ→∂μAμ−□λ\partial_\mu A^\mu \to \partial_\mu A^\mu - \Box\lambda. To reach ∂μAμ=0\partial_\mu A^\mu = 0 from any starting potential we need only solve □λ=∂μAμ\Box\lambda = \partial_\mu A^\mu — an inhomogeneous wave equation, which always has a solution. The gauge is reachable from anywhere. With it imposed the coupling term dies: \boxed{\;\Box A^\nu = \mu_0 J^\nu\;} Why: reachability is not a detail. If □λ=∂μAμ\Box\lambda = \partial_\mu A^\mu had no solutions the Lorenz gauge would be an empty wish. It has them because □\Box is invertible on suitable sources — the very fact that makes retarded potentials exist.

Step 4. Check the choice was covariant. ∂μAμ\partial_\mu A^\mu contracts a lower index with an upper one: a Lorentz scalar. So ∂μ′A′μ=∂μAμ=0\partial'_\mu A'^\mu = \partial_\mu A^\mu = 0 in every frame at once, from one imposition. Why: <strong>the Lorenz gauge is the common choice that is a scalar condition, which is exactly why it is the one that yields (1).</strong> The Coulomb gauge ∇⋅A⃗=0\nabla\cdot\vec A = 0 constrains the spatial part alone: it singles out a time axis, it is not a scalar, and a boost destroys it.

Step 5. Read off the three-vector content. With μ0=1/(ε0c2)\mu_0 = 1/(\varepsilon_0 c^2) and Jμ=(cρ,J⃗)J^\mu = (c\rho, \vec J), the ν=0\nu = 0 component of (1) is □(ϕ/c)=μ0cρ\Box(\phi/c) = \mu_0 c\rho, i.e. \Box\phi = \frac{\rho}{\varepsilon_0},\qquad\text{and}\qquad \Box\vec A = \mu_0\vec J . Why: the two equations of ordinary electrodynamics — but now inside one four-vector equation, and their relative normalisation, which looks arbitrary in three-dimensional language, is fixed by nothing more than JμJ^\mu being a four-vector. *Verified in verify_aq-sr-wave-equation-four-potential.py: FμνF^{\mu\nu} is rebuilt symbolically from Aμ=(ϕ/c,A⃗)A^\mu = (\phi/c,\vec A) with arbitrary functions of (t,x,y,z)(t,x,y,z); the identity ∂μFμν=□Aν−∂ν(∂μAμ)\partial_\mu F^{\mu\nu} = \Box A^\nu - \partial^\nu(\partial_\mu A^\mu) is confirmed component by component; and the reduction to □ϕ=ρ/ε0\Box\phi = \rho/\varepsilon_0, □A⃗=μ0J⃗\Box\vec A = \mu_0\vec J is checked.*

Go deeper · Derivation B · 2 minOwned result: the speed of light is an output, and the null condition makes it frame-independent
Derivation BOwned result: the speed of light is an output, and the null condition makes it frame-independent

Set Jμ=0J^\mu = 0 and ask what the equation permits. The answer is the most consequential deduction in classical physics, and it splits into two claims usually run together. Claim A: the vacuum equation has wave solutions of speed cc. Claim B: every inertial observer measures that same speed. Claim A is nineteenth-century physics and needs only the operator. Claim B needs the operator and the fact that kμkμk_\mu k^\mu is a scalar — and it is the whole of Einstein's second postulate, arrived at without postulating anything at all.

Step 1. Substitute a plane wave. Put Aμ=εμe−ikνxνA^\mu = \varepsilon^\mu e^{-ik_\nu x^\nu} with kμ=(ω/c,k⃗)k^\mu = (\omega/c, \vec k), so kνxν=ωt−k⃗⋅x⃗k_\nu x^\nu = \omega t - \vec k\cdot\vec x. Each ∂μ\partial_\mu brings down −ikμ-ik_\mu, so □→−kμkμ\Box \to -k_\mu k^\mu: \Box A^\mu = 0 \quad\Longrightarrow\quad k_\mu k^\mu\,\varepsilon^\mu = 0 \quad\Longrightarrow\quad k_\mu k^\mu = 0 for any non-zero polarisation. Why: the wave four-vector is forced onto the <strong>light cone in momentum space</strong> — and note what did the forcing: the operator, nothing else.

Step 2. Unpack the null condition. kμkμ=ω2/c2−∣k⃗∣2=0k_\mu k^\mu = \omega^2/c^2 - |\vec k|^2 = 0, hence \omega = c\,|\vec k| ,\qquad v_{\text{phase}} = \frac{\omega}{|\vec k|} = c ,\qquad v_{\text{group}} = \frac{\mathrm{d}\omega}{\mathrm{d}|\vec k|} = c . Why: <strong>phase and group velocity coincide, both equal to cc.</strong> The vacuum is non-dispersive — so familiar it is easy to forget it is a derived property of the operator, and one that fails the moment the wave enters glass.

Step 3. Now the genuinely relativistic part. kμkμk_\mu k^\mu contracts a four-vector with itself: a Lorentz scalar. Zero in one frame, zero in all. Boost along xx: \omega' = \gamma(\omega - \beta c k_x),\qquad c k'_x = \gamma(ck_x - \beta\omega),\qquad k'_y = k_y,\ k'_z = k_z, and direct computation gives ω′2/c2−∣k⃗′∣2=ω2/c2−∣k⃗∣2\omega'^2/c^2 - |\vec k'|^2 = \omega^2/c^2 - |\vec k|^2. Why: <strong>ω\omega changes. ∣k⃗∣|\vec k| changes. Their ratio does not.</strong> That is the relativistic Doppler effect and the invariance of cc — the same statement seen from two sides.

Step 4. State it sharply. A moving observer sees a different frequency and a different wavelength — that is Doppler — but finds ω′=c∣k⃗′∣\omega' = c|\vec k'| with the same numerical cc. \boxed{\;k_\mu k^\mu = 0 \ \text{in one frame} \;\Longleftrightarrow\; k_\mu k^\mu = 0 \ \text{in all frames}\;} Why: <strong>this is Einstein&#39;s second postulate, and here it is a theorem.</strong> Einstein made the invariance of cc an axiom in 1905 because he was prepared to give up Newtonian kinematics to keep Maxwell. The logic here runs the other way: grant Maxwell and the Minkowski metric, and the invariance of cc is not an assumption — it is a line of algebra you cannot avoid. *Verified in verify_aq-sr-wave-equation-four-potential.py: the boost matrix is applied to kμk^\mu and kμ′k′μ−kμkμk'_\mu k'^\mu - k_\mu k^\mu is shown to vanish identically in β\beta; ω=c∣k⃗∣\omega = c|\vec k| is confirmed preserved for boosts up to β=0.999999\beta = 0.999999 in 50-digit arithmetic, where float64 silently fails.*

Go deeper · Derivation C · 2 minOwned result: 4 − 1 − 1 = 2. Counting the polarisations of light out of the gauge freedom
Derivation COwned result: 4 − 1 − 1 = 2. Counting the polarisations of light out of the gauge freedom

The four-potential has four components. Light has two polarisations. The gap is not swept under a rug — it is closed exactly, by two subtractions, and the second is the residual gauge freedom the Lorenz condition famously fails to remove. This is the cleanest place in classical electromagnetism to see a degree-of-freedom count done honestly, and it is the classical shadow of the statement that the photon is massless and has two helicities.

Step 1. The Lorenz condition removes one. Applying ∂μAμ=0\partial_\mu A^\mu = 0 to a plane wave Aμ=εμe−ik⋅xA^\mu = \varepsilon^\mu e^{-ik\cdot x} gives k_\mu \varepsilon^\mu = 0 , one linear constraint on the four-component polarisation vector, leaving three. Why: <strong>the polarisation must be orthogonal to the propagation direction — in the four-dimensional sense.</strong> For null kk this is weaker than it looks, and that is the crux of the next step.

Step 2. The residual gauge removes another. Inside the Lorenz gauge we may still shift Aμ→Aμ−∂μλA^\mu \to A^\mu - \partial^\mu\lambda provided □λ=0\Box\lambda = 0. Taking λ=iα e−ik⋅x\lambda = i\alpha\, e^{-ik\cdot x} — which satisfies □λ=0\Box\lambda = 0 precisely because kμkμ=0k_\mu k^\mu = 0 — shifts the polarisation by \varepsilon^\mu \;\to\; \varepsilon^\mu + \alpha\, k^\mu . Why: <strong>the change is proportional to kμk^\mu itself</strong>, so any part of εμ\varepsilon^\mu along kμk^\mu is pure gauge — set it to zero by choosing α\alpha. Note the economy: the residual gauge exists <strong>only because kk is null</strong>, the same condition that gave us cc.

Step 3. Count. Four components, minus one from k⋅ε=0k\cdot\varepsilon = 0, minus one from ε∼ε+αk\varepsilon \sim \varepsilon + \alpha k: 4 - 1 - 1 = \boxed{\;2\;} and the survivors are the two directions genuinely transverse to k⃗\vec k. Why: <strong>the two transverse polarisations of light, obtained by subtraction.</strong> Use the residual freedom to set ε0=0\varepsilon^0 = 0; then kμεμ=0k_\mu\varepsilon^\mu = 0 collapses to k⃗⋅ε⃗=0\vec k\cdot\vec\varepsilon = 0 — ordinary transversality, now visible as the leftover of a four-dimensional condition.

Step 4. Ask what a photon mass would do. A massive vector field obeys (□+κ2)Aμ=0(\Box + \kappa^2)A^\mu = 0 with κ=mc/ℏ\kappa = mc/\hbar, so kμkμ=κ2≠0k_\mu k^\mu = \kappa^2 \ne 0 — and Step 2 fails, because it required kk to be null. The second subtraction is lost: 4−1=34 - 1 = 3. Why: <strong>a massive spin-1 particle has three polarisations; a massless one has two.</strong> Light's missing longitudinal mode is not missing by accident — it is missing because kμkμ=0k_\mu k^\mu = 0, i.e. because □\Box has no mass term, i.e. because cc is the invariant speed. The polarisation count, the photon's masslessness and the constancy of cc are three names for one fact. *Verified in verify_aq-sr-wave-equation-four-potential.py: the shift εμ→εμ+αkμ\varepsilon^\mu\to\varepsilon^\mu + \alpha k^\mu is shown to preserve both □Aμ=0\Box A^\mu = 0 and k⋅ε=0k\cdot\varepsilon = 0; and the rank of the surviving polarisation space is computed to be exactly 2 for null kk and exactly 3 for timelike kk (the massive case).*

Limiting cases

CaseEquationWhat it becomes
Vacuum, Jμ=0J^\mu = 0□Aμ=0\Box A^\mu = 0plane waves, kμkμ=0k_\mu k^\mu = 0, ω=c∣k⃗∣\omega = c\lvert\vec k\rvert — light
Static, ∂t→0\partial_t \to 0−∇2Aμ=μ0Jμ-\nabla^2 A^\mu = \mu_0 J^\muPoisson: ∇2ϕ=−ρ/ε0\nabla^2\phi = -\rho/\varepsilon_0 ⇒ Coulomb; ∇2A⃗=−μ0J⃗\nabla^2\vec A = -\mu_0\vec J ⇒ Biot–Savart
Point source, retardedAμ=μ0c4πq uμuν(x−w)ν∣retA^\mu = \frac{\mu_0 c}{4\pi}\frac{q\,u^\mu}{u_\nu(x-w)^\nu}\big\vert_{\text{ret}}the Liénard–Wiechert potential (`aq-sr-lienard-wiechert-potentials`)
c→∞c\to\infty (formal Galilean)−∇2Aμ=μ0Jμ-\nabla^2 A^\mu = \mu_0 J^\mu everywhereinstantaneous action at a distance; no radiation, no causal cone
Euclidean signature∇42Aμ=0\nabla_4^2 A^\mu = 0elliptic — harmonic functions, no waves at all
Massive vector field(□+κ2)Aμ=0(\Box + \kappa^2)A^\mu = 0kμkμ=κ2≠0k_\mu k^\mu = \kappa^2 \ne 0: dispersive, three polarisations

Worked examples

1 · The speed of light, computed from a capacitor and a coil

The vacuum wave equation gives phase speed v=1/μ0ε0v = 1/\sqrt{\mu_0\varepsilon_0}. Evaluate it. Neither constant was measured with light.

Step 1 — where the speed comes from. With □=1c2∂t2−∇2\Box = \frac{1}{c^2}\partial_t^2 - \nabla^2 and c2≡1/(μ0ε0)c^2 \equiv 1/(\mu_0\varepsilon_0), the vacuum equation □A⃗=0\Box\vec A = 0 reads \mu_0\varepsilon_0\,\frac{\partial^2\vec A}{\partial t^2} = \nabla^2\vec A , a wave equation of speed v=1/μ0ε0v = 1/\sqrt{\mu_0\varepsilon_0}. μ0\mu_0 is fixed by the force between current-carrying wires; ε0\varepsilon_0 by the force between static charges. Neither has anything to do with optics.

Step 2 — the numbers. With μ0=4π×10−7 H m−1=1.25663706×10−6\mu_0 = 4\pi\times10^{-7}\ \mathrm{H\,m^{-1}} = 1.25663706\times10^{-6} and ε0=8.8541878128×10−12 F m−1\varepsilon_0 = 8.8541878128\times10^{-12}\ \mathrm{F\,m^{-1}}: \mu_0\varepsilon_0 = 1.11265006\times10^{-17}\ \mathrm{s^2\,m^{-2}},\qquad v = \frac{1}{\sqrt{1.11265006\times10^{-17}}} = \mathbf{2.99792458\times10^{8}\ m\,s^{-1}} .

Step 3 — compare with light. The measured speed of light is c=2.99792458×108 m s−1c = 2.99792458\times10^{8}\ \mathrm{m\,s^{-1}} — agreement to all nine figures. That is not a coincidence being reported; it is the reason we believe light is an electromagnetic wave. Maxwell's own version of the comparison, in 1865, was good to about one percent, and he drew the conclusion immediately.

Step 4 — the modern twist. Since the 2019 SI redefinition the logic runs backwards: cc is defined as exactly 299 792 458 m s−1299\,792\,458\ \mathrm{m\,s^{-1}}, and μ0=1.25663706212(19)×10−6 H m−1\mu_0 = 1.25663706212(19)\times10^{-6}\ \mathrm{H\,m^{-1}} is now measured. The identity c=1/μ0ε0c = 1/\sqrt{\mu_0\varepsilon_0} still holds exactly; what changed is which two of the three are handed to you and which you compute. The wave equation fixes the relation; the SI decides which end to nail down.

2 · A green laser, boosted — the frequency changes, the null condition does not

A 532 nm532\ \mathrm{nm} laser propagates along +x+x. Compute kμk^\mu, verify kμkμ=0k_\mu k^\mu = 0, boost to β=0.600\beta = 0.600 along +x+x, and verify the null condition survives.

Step 1 — the lab wave four-vector. ∣k⃗∣=2π/λ=2π/(532×10−9)=1.181050×107 m−1|\vec k| = 2\pi/\lambda = 2\pi/(532\times10^{-9}) = \mathbf{1.181050\times10^{7}\ m^{-1}}, so ω=c∣k⃗∣=3.540698×1015 rad s−1\omega = c|\vec k| = \mathbf{3.540698\times10^{15}\ rad\,s^{-1}} and ν=ω/2π=5.635197×1014 Hz\nu = \omega/2\pi = \mathbf{5.635197\times10^{14}\ Hz}. Hence k^\mu = \left(\frac{\omega}{c},\ k_x,\ 0,\ 0\right) = \big(1.181050\times10^{7},\ 1.181050\times10^{7},\ 0,\ 0\big)\ \mathrm{m^{-1}} .

Step 2 — the null check. k_\mu k^\mu = \frac{\omega^2}{c^2} - k_x^2 = (1.181050\times10^{7})^2 - (1.181050\times10^{7})^2 = \mathbf{0} \quad\text{exactly.} The time and space components are numerically equal. That is what "null" means — the arithmetic signature of a massless field.

Step 3 — boost at β=0.600\beta = 0.600 along +x+x. γ=1.2500\gamma = 1.2500, and both components transform the same way: \frac{\omega'}{c} = \gamma\left(\frac{\omega}{c} - \beta k_x\right) = k'_x = \gamma\left(k_x - \beta\frac{\omega}{c}\right) = 1.25\,(1 - 0.600)\times1.181050\times10^{7} = \mathbf{5.905249\times10^{6}\ m^{-1}} .

Step 4 — read the physics off. Both were multiplied by the same factor, \gamma(1-\beta) = \sqrt{\frac{1-\beta}{1+\beta}} = \sqrt{\frac{0.400}{1.600}} = \mathbf{0.500} , the relativistic Doppler factor. The observer running away at 0.6c0.6c sees the light redshifted by exactly two: λ′=2×532=1064 nm\lambda' = 2\times532 = \mathbf{1064\ nm}, ν′=2.817598×1014 Hz\nu' = \mathbf{2.817598\times10^{14}\ Hz}. And 1064 nm1064\ \mathrm{nm} is the Nd:YAG fundamental whose frequency-doubled output is the 532 nm532\ \mathrm{nm} line we started from — a coincidence of the numbers chosen, but a memorable one.

Step 5 — the invariant. k'_\mu k'^\mu = (5.905249\times10^{6})^2 - (5.905249\times10^{6})^2 = \mathbf{0} . Both components changed by 50%. Their difference of squares did not change at all — it was zero and it stayed zero. Since ω′=c∣k⃗′∣\omega' = c|\vec k'| follows immediately, the boosted observer measures the same cc from a wave whose colour they disagree about. That disagreement is the Doppler effect; that agreement is the second postulate; equation (3) delivers both at once.

Figure
THE WAVE FOUR-VECTOR SLIDES ALONG THE LIGHT CONE AND NEVER LEAVES IT. A 532 nm laser has k^μ = (1.181050×10⁷, 1.181050×10⁷, 0, 0) m⁻¹ — the time and space components are EQUAL, which is exactly what k_μk^μ = ω²/c² − |k|² = 0 means. Boost to β = 0.600 and BOTH components shrink by the Doppler factor √((1−β)/(1+β)) = 0.500: the wavelength stretches 532 nm → 1064 nm, the frequency drops 5.635197×10¹⁴ Hz → 2.817598×10¹⁴ Hz. The observers disagree about the colour of the light. They do NOT disagree about c, because a null vector rescaled is still null — the arrow slides ALONG the 45° line and cannot come off it. Every massive particle's four-momentum sits strictly INSIDE the cone and can be boosted to rest; light sits ON it and cannot. That is the whole difference.
k_μ k^μ = 0 is a Lorentz scalar — so it stays zero k_x ω/c the light cone: ω/c = |k| 45° — where massless fields live LAB 532 nm k^μ = (1.181050, 1.181050)×10⁷ m⁻¹ BOOSTED β = 0.600 (5.905249, 5.905249)×10⁶ m⁻¹ — 1064 nm ×0.500 WHAT CHANGES λ: 532 nm → 1064 nm ν: 5.635197×10¹⁴ → 2.817598×10¹⁴ Hz Doppler factor √((1−β)/(1+β)) = 0.500 WHAT DOES NOT k_μk^μ = ω²/c² − |k|² = 0 before: 0. after: 0. ⇒ ω′ = c|k′|, same c

Physics of the portrait

The portrait is an operator, not a person — because here the operator is the discovery.

Origins · 2 minD'Alembertian and how this was discovered
D
D'Alembertian
1747–present · the wave operator that turned out to be the metric · why the minus sign is the whole of light

The d'Alembertian □=∂μ∂μ=1c2∂t2−∇2\Box = \partial_\mu\partial^\mu = \frac{1}{c^2}\partial_t^2 - \nabla^2 is the Lorentz-invariant wave operator, and it carries the name of Jean le Rond d'Alembert, who in 1747 wrote the one-dimensional wave equation for a vibrating string and solved it as a superposition of a left-moving and a right-moving disturbance. That was a hundred and fifty years before anyone suspected the operator had anything to do with the geometry of spacetime. His equation was about strings; the "speed" in it was a property of tension and mass per unit length; and there was nothing invariant about it whatsoever — a wave on a string looks slower if you run alongside it, and that is the end of the matter. What relativity revealed is that the same operator, applied to the electromagnetic potential, is not a statement about a medium at all. It is ημν∂μ∂ν\eta^{\mu\nu}\partial_\mu\partial_\nu — the Minkowski metric with derivatives in its slots — and its "wave speed" is a property of spacetime, not of any substance, which is why no amount of running alongside changes the answer. The relative minus sign, which for d'Alembert was just what falls out of separating variables, is in Minkowski's reading the signature itself, and it is the sole reason the equation is hyperbolic — the reason its solutions have a light cone, a finite propagation speed, and somewhere for radiation to go. Everything that makes electromagnetism a theory of light rather than a theory of statics is encoded in one minus sign, and d'Alembert put it there ninety years before Maxwell and a hundred and sixty before Minkowski explained what it meant.

History
1747
— d'Alembert publishes the wave equation for a vibrating string, solved as a sum of left- and right-travelling disturbances.
1865
— Maxwell, in "A Dynamical Theory of the Electromagnetic Field", derives waves travelling at 1/√(μ₀ε₀) and identifies light as an electromagnetic wave.
1867
— Ludvig Lorenz publishes the retarded-potential formulation, introducing the gauge condition that now bears his name.
1873
— Maxwell's "Treatise on Electricity and Magnetism" presents the theory in full.
1905
— Einstein takes the invariance of c as a postulate rather than a puzzle, and rebuilds kinematics around it.
1908
— Minkowski's "Raum und Zeit" identifies the wave operator as the metric contraction ∂_μ∂^μ, making the covariance of Maxwell's theory manifest.

Going further

Texts. Griffiths, Introduction to Electrodynamics, 4th ed., §10.1 and §12.3 · Jackson, Classical Electrodynamics, 3rd ed., ch. 6 · Rindler, Introduction to Special Relativity, 2nd ed., ch. 7 · Landau & Lifshitz, The Classical Theory of Fields, 4th ed., §46–47.

Provenance. Signature (+,−,−,−)(+,-,-,-), hence p⋅p=+m2c2p\cdot p = +m^2c^2; mm is the invariant rest mass and no velocity-dependent mass appears anywhere. The dictionary Aμ=(ϕ/c,A⃗)A^\mu = (\phi/c,\vec A), ∂μ=(1c∂t,∇)\partial_\mu = (\tfrac1c\partial_t,\nabla), ∂μ=(1c∂t,−∇)\partial^\mu = (\tfrac1c\partial_t,-\nabla) is used throughout, and Fμν=∂μAν−∂νAμF^{\mu\nu} = \partial^\mu A^\nu - \partial^\nu A^\mu is rebuilt from it in the verification rather than assumed. This article's base declares the correct signature and its physics is sound — one of the two clean bases in this chapter — but its degree-of-freedom count is absent, and is supplied and machine-checked here. Every number printed on this page is recomputed in verify_aq-sr-wave-equation-four-potential.py in `mpmath` at 50 digits, against this article's own printed digits.

Scope. The four-potential belongs to `aq-sr-four-potential`; gauge freedom and the Lorenz condition to `aq-sr-gauge-transformations-lorenz-gauge`; ∂μFμν=μ0Jν\partial_\mu F^{\mu\nu} = \mu_0 J^\nu to `aq-sr-covariant-maxwell-equations`; the four-current to `aq-sr-four-current-density`. The retarded solution for a point charge is developed in `aq-sr-lienard-wiechert-potentials`, the energy those waves carry in `aq-sr-electromagnetic-stress-energy-tensor`. This page owns the reduction of Maxwell to □Aμ=μ0Jμ\Box A^\mu = \mu_0 J^\mu, the derivation of cc from the null condition, and the 4−1−1=24-1-1=2 polarisation count. A parallel treatment from the electromagnetism track lives at `aq-em-covariant-wave-equation`.

Builds on. `aq-sr-four-potential` · `aq-sr-gauge-transformations-lorenz-gauge` · `aq-sr-covariant-maxwell-equations` · `aq-sr-four-current-density`. Leads to. `aq-sr-lienard-wiechert-potentials` · `aq-sr-electromagnetic-stress-energy-tensor` · `aq-sr-covariant-lorentz-force`.

Is the appearance of c in the wave equation really a "derivation" of the speed of light, or is it circular?

It is not circular. μ0\mu_0 and ε0\varepsilon_0 are fixed by static experiments — the force between current-carrying wires, the capacitance of a parallel-plate capacitor. No light is involved in measuring either. The wave equation then predicts, with no further input, waves of speed 1/μ0ε01/\sqrt{\mu_0\varepsilon_0} — and that number equals the independently measured speed of light. A prediction confirmed, not a definition unpacked. What is conventional in the modern SI is which of {c,μ0,ε0}\{c, \mu_0, \varepsilon_0\} you declare exact; the relation among them is physics, and was physics in 1865.

Why is the Lorenz gauge "better" than the Coulomb gauge? Are the physical fields different?

The physical fields are identical — that is what gauge invariance means. The Lorenz gauge is better for this purpose because it is a scalar condition: impose it once, it holds in every frame, and □Aμ=μ0Jμ\Box A^\mu = \mu_0 J^\mu is manifestly covariant. In the Coulomb gauge ∇⋅A⃗=0\nabla\cdot\vec A = 0 the scalar potential obeys ∇2ϕ=−ρ/ε0\nabla^2\phi = -\rho/\varepsilon_0, which looks like instantaneous action at a distance. It is not — the apparent instantaneity cancels against a compensating term in ∂tA⃗\partial_t\vec A — but the cancellation is invisible, whereas in the Lorenz gauge causality is manifest term by term. Use Coulomb gauge to count photons; use Lorenz gauge to see relativity.

If A^μ has four components and light has two polarisations, where do the other two go?

One is eaten by the Lorenz condition, kμεμ=0k_\mu\varepsilon^\mu = 0. The other by the residual gauge freedom: because kk is null, λ∝e−ik⋅x\lambda \propto e^{-ik\cdot x} satisfies □λ=0\Box\lambda = 0, so εμ→εμ+αkμ\varepsilon^\mu \to \varepsilon^\mu + \alpha k^\mu changes nothing physical — the part of εμ\varepsilon^\mu along kμk^\mu is pure gauge. Two genuinely transverse directions remain. The count is exact and it is not a coincidence: it is the classical version of a massless spin-1 particle having two helicities, and it fails for a massive vector field precisely because kμkμ≠0k_\mu k^\mu \ne 0 breaks the residual-gauge construction.

Does the Lorenz gauge condition constrain the physics, or only the bookkeeping?

Only the bookkeeping — with one caveat. Imposing ∂μAμ=0\partial_\mu A^\mu = 0 is always possible (solve □λ=∂μAμ\Box\lambda = \partial_\mu A^\mu), so it excludes no physical configuration. But it is not automatic: consistency requires ∂μJμ=0\partial_\mu J^\mu = 0. Take ∂μ\partial_\mu of equation (1): □(∂μAμ)=μ0 ∂μJμ\Box(\partial_\mu A^\mu) = \mu_0\,\partial_\mu J^\mu. If charge were not conserved the right-hand side would not vanish, the Lorenz condition would not propagate, and the equation would be inconsistent. The wave equation and charge conservation are not independent statements — see `aq-sr-charge-invariance-continuity`.

Why "null" rather than just "the wave travels at c"? What does the language buy?

It buys frame-independence, which "travels at cc" does not obviously have. "The wave travels at cc" is a statement about a speed, and speeds are frame-dependent in every other context — which is exactly why the constancy of cc was a crisis rather than a triviality. "kμkμ=0k_\mu k^\mu = 0" says a scalar vanishes, and scalars are the same for everyone by construction. The null language does not merely restate the fact; it explains why the fact is frame-independent — and it transfers immediately to any field whose quanta live on the light cone.

The base article says the retarded solution "expresses causality". Is causality an input or an output here?

Both, in different senses. The light cone is an output: forced by the hyperbolic character of □\Box, which is forced by the signature. But choosing the retarded Green function over the advanced one is an input — □G=δ4\Box G = \delta^4 has both solutions, and the equation prefers neither. We impose retardation because effects follow causes. The wave equation gives you a cone; it does not tell you which way to go round it. That choice is developed in `aq-sr-lienard-wiechert-potentials`, whose invariant Green function Gret∝θ(Δt) δ((Δx)2)G_{\text{ret}} \propto \theta(\Delta t)\,\delta\big((\Delta x)^2\big) carries the step function that does the choosing.

Before this

Where this leads