Hulse-Taylor Binary Pulsar Test

Pulsar-timing tests have only grown sharper since Hulse-Taylor. The double pulsar PSR J0737−3039, discovered in 2003, contains two detectable pulsars in a tighter 2.45-hour orbit seen nearly edge-on; it delivers six independent post-Keplerian parameters and now tests the radiative sector to better than a part in $10^4$. The same Peters-Mathews physics, extended to high post-Newtonian order with spin and tidal terms, is the analytic backbone of the waveform templates LIGO and Virgo match against their data — the Hulse-Taylor system validated the leading-order term those templates begin from. Pulsar timing arrays now push the idea to nanohertz frequencies, watching dozens of millisecond pulsars for the correlated wobble of a gravitational-wave background from supermassive black-hole binaries, for which the first evidence appeared in 2023.

Picture the two neutron stars as a slowly tumbling mass distribution — a spinning, stretching dumbbell. A configuration whose mass quadrupole moment changes in time must radiate gravitational waves, and the radiation carries away orbital energy. Losing energy, a bound orbit does not slow down; it sinks to a smaller, faster orbit, exactly as a satellite in a thin atmosphere spirals inward and circulates more quickly as it decays. The pulsar itself is the instrument: rather than building a detector, we let the binary be both source and clock, reconstructing the orbit's inspiral from pulse arrival times and watching the period shrink by tens of microseconds over decades. The orbit's high eccentricity is what makes the effect measurable. The quadrupole luminosity scales steeply — as (separation)⁻⁵ and (velocity)⁶ — so almost all the radiation is emitted in a brief, violent burst each time the stars whip through periastron. The first figure shows that geometry: a stretched ellipse whose radiated power is concentrated at one end.

The clean logic of the test comes from over-determination. Each independently measured relativistic effect draws a curve in the plane of the two masses: periastron advance gives one, the Einstein delay another, orbital decay a third. Two curves intersect to fix both masses, leaving the third — the orbital decay — as a prediction with nothing left to adjust. Three relativistic effects, two unknowns: the system is over-constrained, so agreement is a real test rather than a tautology. The definitive way to display that agreement is a growing one — the accumulated shift in the time of periastron, which climbs quadratically and traces out the parabola general relativity draws in advance.

Everything begins with the quadrupole formula for the power radiated by a slowly moving, weakly gravitating source. It is derived in the prerequisite articles; here it is the engine that drives the orbit inward. The luminosity is set by the third time derivative of the reduced quadrupole moment, averaged over an orbit.

where: $L_{GW}$ — power radiated as gravitational waves; $Q_{ij}$ — reduced (trace-free) mass quadrupole moment; $\dddot{Q}_{ij}$ — third time derivative of $Q_{ij}$; $\langle\cdot\rangle$ — average over one orbital period; $G,c$ — Newton's constant and the speed of light. The steep power of $c$ in the denominator is why gravitational radiation is normally so feeble. Applied to a Keplerian ellipse of semi-major axis $a$ and eccentricity $e$, the orbit average of equation (1) is the Peters-Mathews luminosity. The bracketed polynomial in $e$ is the eccentricity enhancement — unity for a circle, diverging as $e\to1$.

where: $m_1,m_2$ — the two stellar masses; $a$ — semi-major axis of the relative orbit; $e$ — orbital eccentricity. For B1913+16 the enhancement factor multiplying the circular-orbit value is about 11.9. Equating that luminosity to the rate of loss of Newtonian orbital energy, and converting to a rate of change of period, yields the celebrated Peters formula — the single equation the whole test rests on.

where: $P_b$ — orbital period; $\dot P_b$ — its rate of change, dimensionless (s/s). Everything on the right is measured independently, so $\dot P_b$ is a prediction, not a fit. The masses in equation (3) come from two other relativistic effects. The first is the secular advance of periastron — the same effect that precesses Mercury's perihelion by 43″ per century, but here 35,000 times larger.

where: $\dot\omega$ — secular advance rate of the longitude of periastron. It depends only on the total mass $m_1+m_2$ and observed orbital elements, so it fixes one combination of the masses. The second is the Einstein delay $\gamma$: the combined gravitational redshift and second-order Doppler shift of the pulsar's clock, which varies around the eccentric orbit and appears as a periodic term in the timing.

where: $\gamma$ — amplitude of the combined redshift and time-dilation term in pulse arrival times. It depends on a different mass combination, so together (4) and (5) solve for both masses. Solving equations (4) and (5) simultaneously gives the two masses to six significant figures — the most precisely weighed stars outside the Solar System.

where: $m_p$ — pulsar mass; $m_c$ — companion (unseen neutron star) mass; $M_\odot$ — solar mass. These feed directly into equation (3). With the masses fixed, the predicted and observed decay rates can be compared as a single ratio — the number that constitutes the test.

where: $\dot P_b^{\rm obs}$ — intrinsic decay after kinematic corrections; $\dot P_b^{\rm GR}$ — the general-relativity quadrupole prediction. The ratio is consistent with unity to about 0.1 percent. Finally, the effect one actually plots is the running total of how early periastron arrives. Because the period drifts linearly, the accumulated timing offset grows quadratically in elapsed time — the parabola of the second figure.

where: $\Delta t_{\rm peri}$ — accumulated shift in the time of periastron; $t$ — elapsed time since a reference epoch. It reaches roughly half a minute after four decades.

The Peters prediction has clean asymptotic regimes, and two of the entries below are corrections that must be applied before the raw data can be compared to theory.

1 · Predicting the orbital decay from the measured masses (fully worked)

Using $P_b=27906.98$ s, $e=0.6171334$, $m_p=1.438\,M_\odot$, and $m_c=1.390\,M_\odot$, predict $\dot P_b$ and compare it to the corrected observed value $-2.398\times10^{-12}$ s/s.

Step 1. Evaluate the eccentricity factor. With $e^2=0.38086$ and $e^4=0.14506$, the numerator is $1+\tfrac{73}{24}(0.38086)+\tfrac{37}{96}(0.14506)=1+1.15845+0.05591=2.21436$, and the denominator is $(1-0.38086)^{7/2}=(0.61914)^{3.5}=0.18677$. So $f(e)=2.21436/0.18677=\mathbf{11.86}$ — the promised near-twelvefold boost over a circular orbit.

Step 2. Assemble the mass factor. In solar masses, $m_p m_c=1.438\times1.390=1.99882$ and $(m_p+m_c)^{1/3}=(2.828)^{1/3}=1.41438$, so $m_p m_c/(m_p+m_c)^{1/3}=1.41321\,M_\odot^{5/3}$. Converting with $M_\odot=1.989\times10^{30}$ kg gives the dimensional factor needed in equation (3).

Step 3. Assemble the period factor: $P_b/2\pi=4441.6$ s, so $(4441.6)^{-5/3}=8.05\times10^{-7}$ s⁻⁵ᐟ³. Multiplying the prefactor $192\pi/5\times(G^{5/3}/c^5)$ by the mass factor of Step 2 and by $f(e)$ from Step 1 yields $\dot P_b^{\,\rm GR}=\mathbf{-2.403\times10^{-12}}$ s/s.

Step 4. Compare. The corrected observed value is $-2.398\times10^{-12}$ s/s, so the ratio is $2.398/2.403=\mathbf{0.998}$, agreeing with the parameter-free prediction to about 0.2 percent. Every input came from other measurements, so this agreement is a genuine confirmation that the binary radiates at exactly the quadrupole rate.

2 · Periastron advance and the contrast with Mercury (fully worked)

Compute the relativistic periastron-advance rate $\dot\omega$ for B1913+16 from equation (4), express it in degrees per year, and compare it to Mercury's perihelion precession of 43″ per century.

Step 1. Gather the factors. The orbital angular frequency gives $(2\pi/P_b)^{5/3}=(2.2515\times10^{-4})^{5/3}=1.828\times10^{-6}$ s⁻⁵ᐟ³, and the total mass $2.828\,M_\odot$ gives $[G(m_1+m_2)]^{2/3}=(3.754\times10^{20})^{2/3}=5.14\times10^{13}$ (SI).

Step 2. Combine with the eccentricity prefactor $3/(1-e^2)=4.845$ and divide by $c^2=8.988\times10^{16}$: $\dot\omega=4.845\times1.828\times10^{-6}\times5.14\times10^{13}/8.988\times10^{16}=\mathbf{2.34\times10^{-9}}$ rad/s.

Step 3. Convert to degrees per year by multiplying by $180/\pi$ and by $3.156\times10^{7}$ s/yr: $\dot\omega=\mathbf{4.23^\circ\ per\ year}$ — the measured value, which closes the mass solution.

Step 4. Contrast with Mercury. Mercury's anomalous precession is $43''$ per century $=0.0119^\circ$ per century $=1.19\times10^{-4}$ degrees per year. The ratio is $4.23/(1.19\times10^{-4})\approx\mathbf{3.5\times10^{4}}$: the pulsar's periastron advances in a single day by more than Mercury's perihelion moves in a century, because $\dot\omega$ scales as $(2\pi/P_b)^{5/3}M_{\rm tot}^{2/3}$ and the binary is both far more massive and vastly more compact than the Sun-Mercury system.

The Hulse-Taylor test is a story of two radio astronomers, a giant dish, and decades of patience — the years of timing that turned one flickering source into the finest confirmation of Einstein's radiative predictions before the interferometer era.

Texts. Misner, Thorne, and Wheeler, Gravitation (1973), §36 (generation of gravitational waves) · Wald, General Relativity (1984), §4.4 · Carroll, Spacetime and Geometry (2004), §7.6 · Hartle, Gravity (2003), §23 · Schutz, A First Course in General Relativity (2009), §9.

Primary literature. Peters and Mathews, Phys. Rev. 131, 435 (1963) — orbit-averaged luminosity; Peters, Phys. Rev. 136, B1224 (1964) — the $\dot a$, $\dot e$ evolution; Hulse and Taylor, ApJ 195, L51 (1975) — discovery; Weisberg and Huang, ApJ 829, 55 (2016) — the definitive timing solution and decay test. The measured values used throughout — $P_b=27906.98$ s, $e=0.6171334$, $m_p=1.438\,M_\odot$, $m_c=1.390\,M_\odot$, corrected $\dot P_b^{\rm obs}\approx-2.398\times10^{-12}$ s/s, ratio $1.0013\pm0.0021$ — are from Weisberg and Huang (2016). The Peters formula derived via energy balance gives $f(0.617)\approx11.86$, $\dot\omega\approx4.23^\circ$/yr from equation (4), and $\dot P_b^{\rm GR}\approx-2.403\times10^{-12}$ s/s from the listed masses and period.

Prerequisites. Quadrupole Formula · Multipole Expansion of Gravitational Radiation · Linearized Einstein Equations · Isaacson Stress-Energy of Gravitational Waves. Leads to. Plus and Cross Polarizations · Geodesic Deviation and the Tidal Action of Waves.

Ask the tutor: ask me to walk through the energy-balance derivation of the Peters formula line by line, or to show exactly how ω̇ and γ pin down both neutron-star masses.