In 1916, Albert Einstein deduced that accelerating mass distributions must perturb the geometry of spacetime, generating ripples in the metric that propagate across the cosmos at the speed of light: Gravitational Waves. Unlike electromagnetic waves, which represent vibrations of a vector field propagating through spacetime, gravitational waves represent ripples of spacetime itself—stretching and compressing spatial dimensions perpendicular to their direction of propagation.

1. Linearized Gravity in the Weak-Field Limit

Far from strong gravitational sources, the metric tensor can be decomposed as flat Minkowski spacetime $eta_{munu} = text{diag}(-1, 1, 1, 1)$ perturbed by an infinitesimal perturbation $h_{munu}$ ($|h_{munu}| ll 1$):

g_{munu} = eta_{munu} + h_{munu}

Defining the trace-reversed perturbation $bar{h}_{munu} = h_{munu} – frac{1}{2} eta_{munu} h$ and adopting the Lorenz gauge condition $partial^mu bar{h}_{munu} = 0$, the non-linear Einstein field equations simplify to an exact, linear wave equation:

Box bar{h}_{munu} = left( -frac{1}{c^2} frac{partial^2}{partial t^2} + nabla^2 right) bar{h}_{munu} = -frac{16pi G}{c^4} T_{munu}

In vacuum ($T_{munu} = 0$), the equation becomes $Box bar{h}_{munu} = 0$, proving that spacetime metric disturbances propagate at precisely the speed of light ($v = c$).

2. The Transverse-Traceless (TT) Gauge

Using the residual coordinate gauge freedoms of general relativity, the physical radiation field outside the source can be cast into the Transverse-Traceless (TT) Gauge ($h_{0mu} = 0$, $h^i_{;i} = 0$, and $partial^j h_{ij} = 0$). For a wave traveling along the $z$-axis ($k^mu = (omega/c, 0, 0, omega/c)$), the metric perturbation tensor possesses only two independent physical polarization states:

h_{ij}^{text{TT}}(t, z) = begin{pmatrix} h_+ & h_times & 0 \ h_times & -h_+ & 0 \ 0 & 0 & 0 end{pmatrix} e^{i(omega t – kz)}

Where $h_+$ represents the “plus” polarization and $h_times$ represents the “cross” polarization (rotated by $45^circ$). A circular ring of freely falling test particles in the $x-y$ plane deformed by an incoming $h_+$ wave stretches along the $x$-axis while compressing along the $y$-axis during the first half-cycle, reversing the distortion during the second half-cycle.

3. The Quadrupole Radiation Formula

Why are gravitational monopoles and dipoles forbidden? Conservation of total mass ($dot{M} = 0$) eliminates monopole radiation. Conservation of total linear momentum ($dot{mathbf{P}} = ddot{mathbf{d}} = 0$) eliminates mass-dipole radiation. Therefore, gravitational radiation begins at the Mass Quadrupole Moment order:

I_{ij} = int rho(t, mathbf{x}) left( x_i x_j – frac{1}{3} r^2 delta_{ij} right) d^3x

In 1918, Einstein derived the total gravitational radiation luminosity emitted by an arbitrary non-relativistic source:

L_{text{GW}} = frac{G}{5 c^5} leftlangle dddot{I}_{jk} dddot{I}^{jk} rightrangle

The coupling constant $G/c^5 approx 2.74 times 10^{-60}text{ s/erg}$ is astonishingly small, explaining why gravitational waves are imperceptible for laboratory masses and require cosmic compact binaries to be detected.

4. Binary Orbital Decay: The Hulse-Taylor Pulsar

For a circular binary system of masses $m_1$ and $m_2$ separated by orbital distance $r$ with angular velocity $Omega$, the emitted gravitational power is:

L_{text{GW}} = frac{32}{5} frac{G^4}{c^5} frac{m_1^2 m_2^2 (m_1 + m_2)}{r^5}

This steady radiation of energy drains orbital angular momentum, causing the orbital period $P_b$ to decrease over time:

dot{P}_b = -frac{192pi}{5 c^5} left( frac{2pi G}{P_b} right)^{5/3} frac{m_1 m_2}{(m_1 + m_2)^{1/3}} left( 1 + frac{73}{24} e^2 + frac{37}{96} e^4 right) (1-e^2)^{-7/2}

Russell Hulse and Joseph Taylor monitored the binary pulsar PSR B1913+16 for three decades, measuring an orbital period decay matching general relativity’s prediction to within $0.2%$, providing the first definitive indirect proof of gravitational waves (Nobel Prize 1993).