Gauge Transformations in Linearized Gravity

Linearized gravity is what Einstein's theory becomes when you look at a nearly flat patch of spacetime and keep only the first whisper of curvature. Write the metric as gμν=ημν+hμνg_{\mu\nu}=\eta_{\mu\nu}+h_{\mu\nu} with ∣hμν∣≪1|h_{\mu\nu}|\ll1, expand the field equations to first order, and the geometry dissolves into something familiar: a symmetric tensor field hμνh_{\mu\nu} on a fixed Minkowski background, obeying wave equations, radiating, carrying energy — a relativistic field theory on flat space. But a catch runs through the whole subject: not every hμνh_{\mu\nu} is physics. The full theory is invariant under arbitrary smooth relabellings of the points of spacetime (diffeomorphisms), and the shadow that invariance casts on the weak field is a gauge symmetry acting on hμνh_{\mu\nu}. Perform an infinitesimal coordinate change xμ→xμ+ξμ(x)x^\mu\to x^\mu+\xi^\mu(x), with ξμ\xi^\mu as small as hh itself, and the perturbation shifts by hμν→hμν−∂μξν−∂νξμh_{\mu\nu}\to h_{\mu\nu}-\partial_\mu\xi_\nu-\partial_\nu\xi_\mu. Two perturbations related this way are the same curved spacetime in two slightly different coordinates. This is not merely an analogy to the gauge freedom Aμ→Aμ−∂μΛA_\mu\to A_\mu-\partial_\mu\Lambda of electromagnetism — it is the same mathematical structure, with ξμ\xi^\mu playing the role of Λ\Lambda. The consequences are practical: raw metric components mean nothing until you fix a gauge or build gauge-invariant combinations; the freedom collapses the field equations to a clean wave equation; and it explains — by a bookkeeping subtraction 10−4−4=210-4-4=2 — why a gravitational wave carries exactly two polarizations. Get gauge wrong and you will, as Einstein himself once did, mistake a coordinate artifact for a physical wave.

Key takeaways
  • The weak-field metric splits as gμν=ημν+hμνg_{\mu\nu}=\eta_{\mu\nu}+h_{\mu\nu}; an infinitesimal coordinate change xμ→xμ+ξμx^\mu\to x^\mu+\xi^\mu shifts the perturbation by hμν→hμν−∂μξν−∂νξμ=hμν−2∂(μξν)h_{\mu\nu}\to h_{\mu\nu}-\partial_\mu\xi_\nu-\partial_\nu\xi_\mu=h_{\mu\nu}-2\partial_{(\mu}\xi_{\nu)}.
  • This is the linearized residue of full diffeomorphism invariance — a genuine redundancy, not a physical change — structurally identical to the gauge freedom of electromagnetism, with ξμ↔Λ\xi^\mu\leftrightarrow\Lambda.
  • The linearized Riemann tensor is invariant under any ξμ\xi^\mu, so all curvature, tidal forces, and wave strain are gauge-independent; individual components of hμνh_{\mu\nu} are not observable.
  • Counting: hμνh_{\mu\nu} has 10 components; the 4 functions ξμ\xi^\mu remove 4 and residual freedom removes 4 more, leaving 2 physical degrees of freedom — the graviton polarizations h+,h×h_+,h_\times.
  • Fixing the harmonic (Lorenz / de Donder) gauge ∂μhˉμν=0\partial^\mu\bar h_{\mu\nu}=0 collapses the field equations to decoupled wave equations □hˉμν=−16πG Tμν\Box\bar h_{\mu\nu}=-16\pi G\,T_{\mu\nu}; residual freedom with □ξμ=0\Box\xi^\mu=0 then reaches the transverse-traceless (TT) gauge.
  • Key idea: under xμ→xμ+ξμx^\mu\to x^\mu+\xi^\mu the perturbation transforms as hμν→hμν−2∂(μξν)h_{\mu\nu}\to h_{\mu\nu}-2\partial_{(\mu}\xi_{\nu)}; this is pure gauge — the linearized Riemann tensor and every curvature invariant are untouched, so physics lives only in gauge-invariant combinations.
The core insight

Under an infinitesimal coordinate change xμ→xμ+ξμx^\mu\to x^\mu+\xi^\mu the metric perturbation transforms as hμν→hμν−∂μξν−∂νξμ=hμν−2∂(μξν)=hμν−Lξημνh_{\mu\nu}\to h_{\mu\nu}-\partial_\mu\xi_\nu-\partial_\nu\xi_\mu=h_{\mu\nu}-2\partial_{(\mu}\xi_{\nu)}=h_{\mu\nu}-\mathcal{L}_\xi\eta_{\mu\nu}; this is a pure gauge redundancy, the exact linear analogue of Aμ→Aμ−∂μΛA_\mu\to A_\mu-\partial_\mu\Lambda, and the linearized Riemann tensor is invariant under it.

Where this lives today

Gauge transformations in linearized gravity are the working language of gravitational-wave astronomy. Every LIGO–Virgo–KAGRA detection is analyzed in TT gauge, where h+h_+ and h×h_\times are the templates matched against the data; matched filtering assumes the eight gauge components have been stripped away and only the invariant strain remains, and numerical-relativity simulations extract the radiated hijTTh^{\mathrm{TT}}_{ij} so it can be compared gauge-invariantly across codes. The same separation of physical from pure-gauge modes drives cosmology: the Bardeen potentials are gauge-invariant combinations of metric and matter perturbations, so cosmic-microwave-background predictions do not depend on how one slices spacetime. At the frontier, the massless spin-2 invariance δhμν=−2∂(μξν)\delta h_{\mu\nu}=-2\partial_{(\mu}\xi_{\nu)} is a constraint any quantum theory of gravity must respect, and it forbids a graviton mass without a delicate reorganization. Pulsar timing arrays now hunt a nanohertz background with the same TT-gauge strain, extending Einstein's 1916 linearization across nine decades of frequency.

The physical picture

Hold two pictures side by side. In the first, hμνh_{\mu\nu} is a tensor field on a rigid flat background, and a gauge transformation is just a smooth relabelling of the background points. Sliding every point by a tiny displacement ξμ(x)\xi^\mu(x) repaints the field's components — as rotating your axes changes the components of an ordinary vector — though no geometry has moved. That is why δhμν=−2∂(μξν)\delta h_{\mu\nu}=-2\partial_{(\mu}\xi_{\nu)} is exactly minus the Lie derivative of the background, −Lξημν-\mathcal{L}_\xi\eta_{\mu\nu}. Infinitely many hμνh_{\mu\nu} therefore describe the same spacetime; they form a "gauge orbit," and only the direction transverse to it — the gauge-invariant part — carries physics.

In the second picture we ask what a real detector feels. Because physics lives only in the gauge-invariant part, the transverse-traceless (TT) gauge throws away eight components and leaves the two — h+h_+ and h×h_\times — that stretch and squeeze proper distances between freely falling masses. A wave sweeping through the Earth distorts a ring of test masses into a pulsing ellipse, and that distortion is a real, measured, staggeringly small number.

Figure
The two polarizations h₊ and h× stretch and squeeze a ring of free test masses in orthogonal patterns; for GW150914 the peak strain was h ≈ 1.0×10⁻²¹, so across a LIGO arm of length L = 4 km the length change is ΔL = hL ≈ 4×10⁻¹⁸ m — a few thousandths of a proton diameter
What survives in TT gauge: the two polarizations a detector actually feels h₊ (plus) h× (cross) GW150914 peak strain h ≈ 1.0×10⁻²¹ · LIGO arm L = 4 km ΔL = hL ≈ 4×10⁻¹⁸ m ≈ a few thousandths of a proton diameter

The core equations

Start from the split that defines the weak field. Every equation is an expansion in hμνh_{\mu\nu} to first order.

gμν=ημν+hμν,∣hμν∣≪1(1)g_{\mu\nu}=\eta_{\mu\nu}+h_{\mu\nu},\qquad |h_{\mu\nu}|\ll1 \qquad (1)

where: ημν\eta_{\mu\nu} — Minkowski metric, signature (−,+,+,+)(-,+,+,+); hμνh_{\mu\nu} — symmetric perturbation, first order. Indices are raised and lowered with η\eta. The central object is the gauge transformation — what an infinitesimal coordinate change xμ→xμ+ξμx^\mu\to x^\mu+\xi^\mu does to the perturbation. Derivation A derives it from the tensor transformation law.

hμν→hμν′=hμν−∂μξν−∂νξμ=hμν−2∂(μξν)(2)h_{\mu\nu}\to h'_{\mu\nu}=h_{\mu\nu}-\partial_\mu\xi_\nu-\partial_\nu\xi_\mu=h_{\mu\nu}-2\partial_{(\mu}\xi_{\nu)} \qquad (2)

where: ξμ(x)\xi^\mu(x) — infinitesimal gauge vector, of order hh; ∂(μξν)≡12(∂μξν+∂νξμ)\partial_{(\mu}\xi_{\nu)}\equiv\tfrac12(\partial_\mu\xi_\nu+\partial_\nu\xi_\mu); the variation equals −Lξημν-\mathcal{L}_\xi\eta_{\mu\nu}. Gauge freedom never touches observables because the linearized Riemann tensor — which encodes tidal forces and wave strain — is invariant under (2).

Rμνρσ(1)=12(∂ν∂ρhμσ+∂μ∂σhνρ−∂μ∂ρhνσ−∂ν∂σhμρ)(3)R^{(1)}_{\mu\nu\rho\sigma}=\tfrac12\big(\partial_\nu\partial_\rho h_{\mu\sigma}+\partial_\mu\partial_\sigma h_{\nu\rho}-\partial_\mu\partial_\rho h_{\nu\sigma}-\partial_\nu\partial_\sigma h_{\mu\rho}\big) \qquad (3)

where: Rμνρσ(1)R^{(1)}_{\mu\nu\rho\sigma} — Riemann tensor to first order in hh, antisymmetric in [μν][\mu\nu] and [ρσ][\rho\sigma]. Derivation B shows δRμνρσ(1)=0\delta R^{(1)}_{\mu\nu\rho\sigma}=0 under (2). To write the field equations compactly it pays to trace-reverse the perturbation — the change of variable that makes harmonic gauge look clean.

hˉμν≡hμν−12ημνh,h≡ημνhμν,hˉ=−h(4)\bar h_{\mu\nu}\equiv h_{\mu\nu}-\tfrac12\eta_{\mu\nu}h,\qquad h\equiv\eta^{\mu\nu}h_{\mu\nu},\qquad \bar h=-h \qquad (4)

where: hˉμν\bar h_{\mu\nu} — trace-reversed perturbation; hh — the trace. Reversing twice returns the original. In these variables the linearized Einstein equations — valid in any gauge — read as one wave operator plus three divergence terms, which a gauge choice can kill.

□hˉμν−∂ρ∂νhˉρμ−∂ρ∂μhˉρν+ημν∂ρ∂σhˉρσ=−16πG Tμν(5)\Box\bar h_{\mu\nu}-\partial_\rho\partial_\nu\bar h^\rho{}_\mu-\partial_\rho\partial_\mu\bar h^\rho{}_\nu+\eta_{\mu\nu}\partial_\rho\partial_\sigma\bar h^{\rho\sigma}=-16\pi G\,T_{\mu\nu} \qquad (5)

where: □=ηρσ∂ρ∂σ\Box=\eta^{\rho\sigma}\partial_\rho\partial_\sigma — flat d'Alembertian; TμνT_{\mu\nu} — stress-energy at zeroth order in hh; GG — Newton's constant; c=1c=1. Every term after □hˉμν\Box\bar h_{\mu\nu} carries a divergence ∂ρhˉρμ\partial^\rho\bar h_{\rho\mu}. Imposing the harmonic (de Donder) gauge sets that divergence to zero — four conditions, one per free index ν\nu. Derivation C proves it is always reachable.

∂μhˉμν=0(6)\partial^\mu\bar h_{\mu\nu}=0 \qquad (6)

where: ∂μhˉμν\partial^\mu\bar h_{\mu\nu} — divergence of hˉ\bar h. "Harmonic" reflects that the coordinates satisfy □xμ=0\Box x^\mu=0; "de Donder" and "Lorenz gauge" are the gravitational and electromagnetic names. With (6) imposed, the last three terms of (5) vanish and each component of hˉμν\bar h_{\mu\nu} obeys the same sourced wave equation — predicting gravitational radiation at cc.

□hˉμν=−16πG Tμν(7)\Box\bar h_{\mu\nu}=-16\pi G\,T_{\mu\nu} \qquad (7)

where: the operator acts component-by-component, so the ten equations decouple and solve by retarded integrals. In vacuum, □hˉμν=0\Box\bar h_{\mu\nu}=0. The residual freedom left after (6) reaches the transverse-traceless gauge, where a vacuum wave shows only its two physical polarizations.

h0μTT=0,∂jhijTT=0,δijhijTT=0(8)h^{\mathrm{TT}}_{0\mu}=0,\qquad \partial^j h^{\mathrm{TT}}_{ij}=0,\qquad \delta^{ij}h^{\mathrm{TT}}_{ij}=0 \qquad (8)

where: hijTTh^{\mathrm{TT}}_{ij} — the purely spatial, transverse, traceless perturbation. In TT gauge hˉijTT=hijTT\bar h^{\mathrm{TT}}_{ij}=h^{\mathrm{TT}}_{ij}, and the survivors are exactly h+h_+ and h×h_\times.

Derivation AThe gauge law from an infinitesimal change of coordinates

The goal is to derive equation (2) from first principles — to show that hμν→hμν−∂μξν−∂νξμh_{\mu\nu}\to h_{\mu\nu}-\partial_\mu\xi_\nu-\partial_\nu\xi_\mu is not an ad-hoc rule but the unavoidable consequence of writing the same metric in shifted coordinates xμ→x′μ=xμ+ξμ(x)x^\mu\to x'^\mu=x^\mu+\xi^\mu(x), kept to first order, and to identify the variation with a Lie derivative that makes the "only a relabelling" slogan precise.

Step 1. Apply the tensor transformation law to the full metric: gμν′(x′)=∂xα∂x′μ∂xβ∂x′ν gαβ(x)g'_{\mu\nu}(x')=\dfrac{\partial x^\alpha}{\partial x'^\mu}\dfrac{\partial x^\beta}{\partial x'^\nu}\,g_{\alpha\beta}(x), exact and saying nothing yet about smallness. Why: the metric is a rank-2 tensor, so its components in a new chart are fixed entirely by the Jacobian of the coordinate change.

Step 2. Invert the coordinate change. From x′μ=xμ+ξμx'^\mu=x^\mu+\xi^\mu we get xμ=x′μ−ξμ+O(ξ2)x^\mu=x'^\mu-\xi^\mu+O(\xi^2), so the Jacobian is ∂xα∂x′μ=δμα−∂μξα+O(ξ2)\dfrac{\partial x^\alpha}{\partial x'^\mu}=\delta^\alpha_\mu-\partial_\mu\xi^\alpha+O(\xi^2). Why: to first order the shift and its inverse differ only in the sign of ξ\xi, and the derivative supplies the ∂ξ\partial\xi pieces that become the gauge term.

Step 3. Substitute and linearize. Insert the Jacobians and gαβ=ηαβ+hαβg_{\alpha\beta}=\eta_{\alpha\beta}+h_{\alpha\beta}, drop products of two small quantities, and use ∂μηαβ=0\partial_\mu\eta_{\alpha\beta}=0; the cross-terms collapse to gμν′=ημν+hμν−∂μξν−∂νξμg'_{\mu\nu}=\eta_{\mu\nu}+h_{\mu\nu}-\partial_\mu\xi_\nu-\partial_\nu\xi_\mu. Why: the constant background survives untouched, so the entire effect of the coordinate change is dumped into the perturbation — gauge acts on hh, not on η\eta.

Step 4. Read off the transformation and recognize it. Since gμν′=ημν+hμν′g'_{\mu\nu}=\eta_{\mu\nu}+h'_{\mu\nu}, matching gives hμν′=hμν−∂μξν−∂νξμh'_{\mu\nu}=h_{\mu\nu}-\partial_\mu\xi_\nu-\partial_\nu\xi_\mu, which is (2). The subtracted piece is ∂μξν+∂νξμ=Lξημν\partial_\mu\xi_\nu+\partial_\nu\xi_\mu=\mathcal{L}_\xi\eta_{\mu\nu}. Why: identifying δhμν=−Lξημν\delta h_{\mu\nu}=-\mathcal{L}_\xi\eta_{\mu\nu} shows the transformation is the flat background dragged along the flow ξ\xi — the exact linear image of a diffeomorphism.

Go deeper · Derivation B · 1 minWhy curvature does not care about gauge
Derivation BWhy curvature does not care about gauge

The goal is to prove that the linearized Riemann tensor (3) is exactly invariant under (2), so that δRμνρσ(1)=0\delta R^{(1)}_{\mu\nu\rho\sigma}=0 for every ξμ\xi^\mu. This is the most important fact in the subject: no measurement built from curvature — tidal acceleration, geodesic deviation, interferometer strain — can depend on the gauge, so the raw components of hμνh_{\mu\nu} are not observable. The strategy is direct substitution followed by a cancellation that rides on the commuting of partial derivatives.

Step 1. Feed the variation into each slot. Replace every hh in (3) by δhαβ=−∂αξβ−∂βξα\delta h_{\alpha\beta}=-\partial_\alpha\xi_\beta-\partial_\beta\xi_\alpha; the first term gives ∂ν∂ρ δhμσ=−∂ν∂ρ∂μξσ−∂ν∂ρ∂σξμ\partial_\nu\partial_\rho\,\delta h_{\mu\sigma}=-\partial_\nu\partial_\rho\partial_\mu\xi_\sigma-\partial_\nu\partial_\rho\partial_\sigma\xi_\mu, and similarly for the other three. Why: the linearized Riemann tensor is linear in hh, so its variation is the same formula with h→δhh\to\delta h — no product rule, no second-order debris.

Step 2. Lay out all eight third-derivative terms. Each of the four terms in (3) splits into two, giving eight of the schematic form ± ∂∂∂ξ\pm\,\partial\partial\partial\xi, each with three of {μ,ν,ρ,σ}\{\mu,\nu,\rho,\sigma\} on the derivatives and one on ξ\xi. Why: the cancellation is not a slogan but an explicit pairing, best seen written out.

Step 3. Pair the terms using commuting derivatives. The two terms carrying ξσ\xi_\sigma are −12∂ν∂ρ∂μξσ-\tfrac12\partial_\nu\partial_\rho\partial_\mu\xi_\sigma and +12∂ν∂ρ∂μξσ+\tfrac12\partial_\nu\partial_\rho\partial_\mu\xi_\sigma; since partial derivatives commute the ordering is irrelevant and they cancel. Why: the antisymmetries [μν][\mu\nu] and [ρσ][\rho\sigma] of Riemann generate every ∂∂∂ξ\partial\partial\partial\xi term a second time with opposite sign, and equality of mixed partials makes the two the same object.

Step 4. Sum the pairs. Every term finds a partner of opposite sign, so the total vanishes: δRμνρσ(1)=0\delta R^{(1)}_{\mu\nu\rho\sigma}=0 for arbitrary ξμ\xi^\mu. Why: gauge changes shuffle coordinates, and coordinates are exactly what the Riemann tensor is built to be independent of.

Go deeper · Derivation C · 1 minHarmonic gauge and the wave equation
Derivation CHarmonic gauge and the wave equation

The goal is threefold: to show the harmonic gauge ∂μhˉμν=0\partial^\mu\bar h_{\mu\nu}=0 can always be reached, to use it to collapse the messy field equation (5) into the clean wave equation (7), and to expose the consistency condition it demands of the source. This is where gauge freedom becomes a tool, turning Einstein's weak-field equations into the same wave equation that governs light — and where the link between gauge symmetry and energy-momentum conservation shows itself.

Step 1. Find how the divergence transforms. Under (2), δhˉμν=−∂μξν−∂νξμ+ημν∂ρξρ\delta\bar h_{\mu\nu}=-\partial_\mu\xi_\nu-\partial_\nu\xi_\mu+\eta_{\mu\nu}\partial_\rho\xi^\rho; taking its divergence and cancelling the trace piece gives ∂μhˉμν′=∂μhˉμν−□ξν\partial^\mu\bar h'_{\mu\nu}=\partial^\mu\bar h_{\mu\nu}-\Box\xi_\nu. Why: the divergence is exactly what the gauge condition constrains, so its transformation law decides whether the condition is reachable.

Step 2. Solve for the required ξ\xi. Setting ∂μhˉμν′=0\partial^\mu\bar h'_{\mu\nu}=0 demands □ξν=∂μhˉμν\Box\xi_\nu=\partial^\mu\bar h_{\mu\nu}, an inhomogeneous wave equation that always admits a retarded-Green's-function solution — so a gauge vector reaching harmonic gauge always exists. Why: reachability is guaranteed by the solvability of the wave equation, which is why every weak field can be put in harmonic form.

Step 3. Collapse the field equations. With ∂ρhˉρμ=0\partial^\rho\bar h_{\rho\mu}=0, the three divergence terms in (5) vanish and only □hˉμν=−16πG Tμν\Box\bar h_{\mu\nu}=-16\pi G\,T_{\mu\nu} survives — each component an independent sourced wave equation. Why: a single operator replaces the tangle of coupled second derivatives, and gravitational radiation falls out as its retarded solution.

Step 4. Read off the consistency condition. Take the divergence of (7): the left side gives ∂μ□hˉμν=□(∂μhˉμν)=0\partial^\mu\Box\bar h_{\mu\nu}=\Box(\partial^\mu\bar h_{\mu\nu})=0 by the gauge condition, so the right side forces ∂μTμν=0\partial^\mu T_{\mu\nu}=0. Local energy-momentum conservation is not optional — it is required for the gauge-fixed equations to be consistent. Why: this is the linearized Bianchi identity, the field-theory shadow of the fact that gauge invariance and conservation of the source are two faces of one symmetry, just as ∂μJμ=0\partial^\mu J_\mu=0 follows from Lorenz-gauge invariance in electromagnetism.

Limiting cases

Each row is the same theory viewed through a different gauge or source, showing how many physical fields survive.

CaseResultMeaning
Vacuum, harmonic gauge□hˉμν=0,  ∂μhˉμν=0\Box\bar h_{\mu\nu}=0,\ \ \partial^\mu\bar h_{\mu\nu}=0Residual freedom □ξμ=0\Box\xi^\mu=0 removes 4 further components, reaching TT gauge: exactly 2 polarizations propagate at cc.
Static source (Newtonian limit)∇2hˉ00=−16πG ρ\nabla^2\bar h_{00}=-16\pi G\,\rhoTime-independent gauge fixes the spatial coordinates; Φ=−14hˉ00\Phi=-\tfrac14\bar h_{00} obeys Poisson's equation ∇2Φ=4πGρ\nabla^2\Phi=4\pi G\rho.
Pure gaugehμν=−2∂(μξν)h_{\mu\nu}=-2\partial_{(\mu}\xi_{\nu)}For any smooth ξμ\xi^\mu the linearized Riemann tensor vanishes identically — flat spacetime in disguise, no physics at all.
Detector scale ⁣ξ¨ i=−Ri0j0(1) ξj\ddot{\!\xi}^{\,i}=-R^{(1)}_{i0j0}\,\xi^jGeodesic deviation depends only on gauge-invariant R(1)R^{(1)}; the "motion" of a single mass is pure coordinate.

Worked examples

1 · The strain of GW150914 and why it is gauge-invariant (fully worked)

A gravitational wave from a binary black-hole merger passes through a LIGO detector with orthogonal arms of length L=4L=4 km. At peak the dimensionless strain is h≈1.0×10−21h\approx1.0\times10^{-21}. Find the change in arm length, confirm it is the TT strain and not a coordinate artifact, and estimate the wave's reduced wavelength at the ∼100\sim100 Hz band.

Step 1. In TT gauge the proper distance between two freely falling masses separated by coordinate distance LL along a polarization axis is ℓ≈L(1+12hTT)\ell\approx L(1+\tfrac12 h^{\mathrm{TT}}), so δℓ/L=12hTT\delta\ell/L=\tfrac12 h^{\mathrm{TT}}. Detectors quote the strain h≡ΔL/Lh\equiv\Delta L/L as the differential effect across the two arms, restoring the factor of two: ΔL=hL\Delta L=hL. The length change is ΔL=hL=(1.0×10−21)(4×103 m)=4×10−18 m\Delta L=hL=(1.0\times10^{-21})(4\times10^3\ \mathrm{m})=4\times10^{-18}\ \mathrm{m}.

Step 2. Put that in perspective. A proton has diameter ≈1.7×10−15\approx1.7\times10^{-15} m, so ΔL/dproton≈(4×10−18)/(1.7×10−15)≈2.4×10−3\Delta L/d_{\text{proton}}\approx(4\times10^{-18})/(1.7\times10^{-15})\approx2.4\times10^{-3}. The mirrors move by a few thousandths of a proton diameter — the measurement that earned the 2017 Nobel Prize.

Step 3. Confirm gauge-invariance. A detector reads the relative proper acceleration of its test masses, governed by ξ¨ i=−Ri0j0(1)ξj\ddot{\xi}^{\,i}=-R^{(1)}_{i0j0}\xi^j. Because Rμνρσ(1)R^{(1)}_{\mu\nu\rho\sigma} is gauge-invariant (Derivation B), the same displacement is predicted in any gauge — TT gauge merely makes h+,h×h_+,h_\times manifest. A different gauge gives differently-looking hμνh_{\mu\nu} but the identical measurable ΔL\Delta L.

Step 4. Estimate the wavelength. With f≈100f\approx100 Hz near merger, λ=c/f=(3.0×108)/100=3.0×106\lambda=c/f=(3.0\times10^8)/100=3.0\times10^6 m =3000=3000 km, vastly larger than the 4 km arms — which is why the long-wavelength TT-gauge description as a nearly uniform strain across the detector is valid. λ≈3000\lambda\approx3000 km ≫L\gg L.

2 · The Newtonian limit of the Sun (fully worked)

Show that harmonic gauge reproduces Newtonian gravity, and check the weak-field assumption ∣hμν∣≪1|h_{\mu\nu}|\ll1 at the surface of the Sun (M=1.99×1030M=1.99\times10^{30} kg, R=6.96×108R=6.96\times10^8 m).

Step 1. For a static, slow source only T00=ρT_{00}=\rho matters and □→∇2\Box\to\nabla^2, so (7) gives ∇2hˉ00=−16πGρ\nabla^2\bar h_{00}=-16\pi G\rho. Define Φ≡−14hˉ00\Phi\equiv-\tfrac14\bar h_{00}; then ∇2Φ=−14(−16πGρ)=4πGρ\nabla^2\Phi=-\tfrac14(-16\pi G\rho)=4\pi G\rho. This is exactly Poisson's equation ∇2Φ=4πGρ\nabla^2\Phi=4\pi G\rho — Newtonian gravity recovered from a gauge-fixed linearized Einstein equation.

Step 2. Undo the trace reversal. With hˉ00=−4Φ\bar h_{00}=-4\Phi, the relation hμν=hˉμν−12ημνhˉh_{\mu\nu}=\bar h_{\mu\nu}-\tfrac12\eta_{\mu\nu}\bar h gives h00=−2Φh_{00}=-2\Phi and hij=−2Φ δijh_{ij}=-2\Phi\,\delta_{ij}, so ds2=−(1+2Φ) dt2+(1−2Φ) δijdxidxj\mathrm{d}s^2=-(1+2\Phi)\,\mathrm{d}t^2+(1-2\Phi)\,\delta_{ij}\mathrm{d}x^i\mathrm{d}x^j. The single potential Φ\Phi sets both the time-time and space-space parts of the metric.

Step 3. Check weak-field validity at the Sun's surface via ∣h00∣=2GM/(Rc2)=2(6.674×10−11)(1.99×1030)/[(6.96×108)(2.998×108)2]|h_{00}|=2GM/(Rc^2)=2(6.674\times10^{-11})(1.99\times10^{30})/\big[(6.96\times10^8)(2.998\times10^8)^2\big]. This evaluates to ∣h00∣≈4.2×10−6|h_{00}|\approx4.2\times10^{-6}, and Φ/c2≈2.1×10−6\Phi/c^2\approx2.1\times10^{-6}.

Step 4. Interpret. Since 4.2×10−6≪14.2\times10^{-6}\ll1, linearization is superb at the Sun's surface — corrections of order h2∼10−11h^2\sim10^{-11} are why solar-system tests probe general relativity only at parts-per-million. The weak-field expansion is valid, with time-independent gauge freedom fixing the isotropic spatial coordinates above.

Physics of the portrait

The gauge structure of linearized gravity is inseparable from one man, who made both the discovery and the classic mistake — Einstein predicted gravitational waves, then twenty years later nearly convinced himself they were a fiction, tripped up by the gauge freedom he had helped uncover.

Origins · 1 minAlbert Einstein and how this was discovered
AE
Albert Einstein
1879–1955 · weak-field gravity, gravitational waves, and the gauge that fooled him

Having completed general relativity in 1915, Albert Einstein turned at once to its weak-field limit. In 1916 he wrote the metric as a small perturbation of flat space, introduced what we now call harmonic coordinates, and found gravitational waves travelling at the speed of light; a 1918 paper corrected his own factor-of-two error and gave the quadrupole formula still used today. But the story that best illuminates gauge came later. In 1936 Einstein and Nathan Rosen submitted a paper to Physical Review titled "Are Gravitational Waves Real?" arguing — from an apparent singularity in a cylindrical-wave solution — that they were not. The referee (Howard Percy Robertson) recognized the singularity as a mere coordinate artifact, removable by a change of gauge; an indignant Einstein withdrew the paper but on reflection reached the same conclusion, and the 1937 version reversed the claim. Even Einstein could mistake a gauge effect for physics — only the invariance of curvature settles the matter. Théophile de Donder had meanwhile (1921) formalized the harmonic condition ∂μhˉμν=0\partial^\mu\bar h_{\mu\nu}=0, and in 1939 Fierz and Pauli recast the structure as the gauge theory of a massless spin-2 field. Vindication came in 2015: LIGO recorded a strain of 10−2110^{-21} and heard two black holes merge.

History
1916
— Einstein linearizes his field equations and predicts gravitational waves
1921
— de Donder formalizes the harmonic coordinate condition ∂ᵘ h̄_μν = 0
1936
— Einstein and Rosen argue waves are unphysical; a referee catches the gauge error
1939
— Fierz and Pauli recast linearized gravity as the massless spin-2 gauge theory
2016
— LIGO announces GW150914, the first direct detection of gravitational waves

Going further

Texts. Misner, Thorne & Wheeler, Gravitation, §18.1–18.3 (linearized theory and gauge) · Wald, General Relativity, §4.4 (linearized equations, degrees of freedom) · Carroll, Spacetime and Geometry, §7.1 (gauge, harmonic and TT gauge) · Schutz, A First Course in General Relativity, §8.3 · Hartle, Gravity, Ch. 21 · Maggiore, Gravitational Waves Vol. 1, §1.2. Key self-consistency checks: the gauge law δhμν=−2∂(μξν)\delta h_{\mu\nu}=-2\partial_{(\mu}\xi_{\nu)} follows from the coordinate transformation (Derivation A); substituting it into Rμνρσ(1)R^{(1)}_{\mu\nu\rho\sigma} confirms invariance (Derivation B); the wave equation □hˉμν=−16πG Tμν\Box\bar h_{\mu\nu}=-16\pi G\,T_{\mu\nu} requires ∂μTμν=0\partial^\mu T_{\mu\nu}=0 for consistency (Derivation C); the degree-of-freedom count gives 10−4−4=210-4-4=2; and the anchor numbers are ΔL=hL≈4×10−18\Delta L=hL\approx4\times10^{-18} m and ∣h00∣=2GM⊙/(R⊙c2)≈4.2×10−6|h_{00}|=2GM_\odot/(R_\odot c^2)\approx4.2\times10^{-6}. The factor 16πG16\pi G assumes c=1c=1 and signature (−,+,+,+)(-,+,+,+).

Prerequisites. The Weak-Field Approximation · The Einstein Field Equations · The Einstein Tensor · Local Conservation of Energy-Momentum. Leads to. The Newtonian Limit · Variational Derivation of the Field Equations · The Stress-Energy Tensor.

Ask the tutor: Want to see the eight ∂∂∂ξ\partial\partial\partial\xi terms cancel line by line, or walk through how the residual □ξμ=0\Box\xi^\mu=0 freedom is spent to reach TT gauge for a plane wave?

Is a gauge transformation the same thing as a coordinate transformation?

Almost. A gauge transformation in linearized gravity is an infinitesimal coordinate transformation xμ→xμ+ξμx^\mu\to x^\mu+\xi^\mu, re-read as an active change of the field hμνh_{\mu\nu} on a fixed flat background rather than a passive relabelling; the two viewpoints give identical formulas. We call it "gauge" to stress the electromagnetic analogy and to emphasize that it is a redundancy of description — two hμνh_{\mu\nu} related this way are physically indistinguishable. Finite coordinate transformations in the full theory are diffeomorphisms; these gauge transformations are their linearized shadow.

Why exactly two degrees of freedom — where do the other eight go?

Start with the 10 components of the symmetric hμνh_{\mu\nu}. The four gauge functions ξμ\xi^\mu impose the four harmonic conditions ∂μhˉμν=0\partial^\mu\bar h_{\mu\nu}=0, leaving six. But harmonic gauge does not use up all the freedom: any further ξμ\xi^\mu with □ξμ=0\Box\xi^\mu=0 preserves the condition, and this residual freedom removes four more (setting h0μ=0h_{0\mu}=0 and the trace to zero). Six minus four is two — the transverse-traceless components h+h_+ and h×h_\times, the helicity-±2\pm2 graviton states. The count 10−4−4=210-4-4=2 is the gravitational analogue of the photon's 4−2=24-2=2, with an extra subtraction because the graviton is spin-2.

Is the harmonic gauge the same as the Lorenz gauge in electromagnetism?

They are exact analogues, term for term. The Lorenz gauge ∂μAμ=0\partial^\mu A_\mu=0 turns Maxwell's equations into □Aμ=−μ0Jμ\Box A_\mu=-\mu_0 J_\mu with residual freedom □Λ=0\Box\Lambda=0; the harmonic (de Donder) gauge ∂μhˉμν=0\partial^\mu\bar h_{\mu\nu}=0 turns Einstein's into □hˉμν=−16πG Tμν\Box\bar h_{\mu\nu}=-16\pi G\,T_{\mu\nu} with residual freedom □ξμ=0\Box\xi^\mu=0. The dictionary is Aμ↔hμνA_\mu\leftrightarrow h_{\mu\nu}, Λ↔ξμ\Lambda\leftrightarrow\xi^\mu, Jμ↔TμνJ_\mu\leftrightarrow T_{\mu\nu}, and the consistency conditions match too: ∂μJμ=0\partial^\mu J_\mu=0 versus ∂μTμν=0\partial^\mu T_{\mu\nu}=0. The one difference is the extra index — gravity's source is a tensor, so there is one more layer of gauge fixing to reach the physical polarizations.

If gauge changes h_μν, how can gravitational-wave detectors measure anything?

Because detectors do not measure hμνh_{\mu\nu} — they measure the relative proper acceleration of freely falling test masses, governed by ξ¨ i=−Ri0j0(1)ξj\ddot{\xi}^{\,i}=-R^{(1)}_{i0j0}\xi^j. Since the linearized Riemann tensor is gauge-invariant, so is that acceleration and so is the interferometer fringe pattern. TT gauge is chosen not because it is "the real one" but because it is convenient: the mirrors stay at fixed coordinates and the whole physical effect sits in h+,h×h_+,h_\times. Any other gauge gives a differently-looking hμνh_{\mu\nu} but the identical measured strain — gauge invariance is what makes the measurement well-defined at all.

What actually went wrong in the 1936 Einstein–Rosen argument?

Einstein and Rosen studied cylindrical gravitational waves and found the metric became singular on an axis, reading this as evidence that no regular wave solution exists. The flaw was gauge: their coordinate system broke down on the axis, producing a coordinate singularity, not a curvature one — the Riemann tensor stayed finite, and a change to well-behaved coordinates removes the pathology. The referee (later revealed to be H. P. Robertson) spotted this, and after initial fury Einstein agreed, publishing a corrected paper in 1937. The lesson is the warning box: a feature living in the metric components but absent from the curvature invariants is an artifact of the chart, not a fact about spacetime.

Does the gauge freedom mean gravitational waves carry no energy?

No — but gauge is at the heart of the subtlety. A wave's energy is not localizable in a gauge-invariant, pointwise way, because you can always transform to a local freely-falling frame where the metric's first derivatives vanish. What is well-defined is the energy averaged over a region several wavelengths across, captured by the Isaacson stress-energy tensor built from the gauge-invariant hijTTh^{\mathrm{TT}}_{ij}; averaged, the radiated power reduces to Einstein's quadrupole formula. So waves do carry energy — it is why binary orbits decay — but gauge freedom forbids assigning that energy to a point, only to a region.

Before this

Where this leads