Hulse-Taylor Binary

A clock in a shrinking orbit

In the summer of 1974, a graduate student named Russell Hulse was sifting through radio data from the giant Arecibo dish in Puerto Rico, hunting for pulsars — neutron stars whose lighthouse beams sweep past Earth as regular radio pulses. He found one, designated PSR B1913+16, with a spin period of about 59 milliseconds. But something was wrong. The pulse period wasn't constant: over a few hours it sped up and slowed down, drifting by hundreds of microseconds in a smooth, repeating cycle.

The explanation, worked out with his adviser Joseph Taylor, was that the pulsar is not alone. It orbits a companion — almost certainly a second neutron star — every 7.75 hours. As the pulsar swings toward us its pulses arrive a little early; as it swings away they arrive a little late. The whole system is a binary pulsar, the first ever found. And because the pulsar's pulses are as steady as an atomic clock, the binary became the most precise gravitational laboratory in the sky. Tracking the arrival times over years revealed something that no other instrument could see at the time: the orbit is slowly, inexorably shrinking.

Why the orbit loses energy

General relativity says that any mass distribution whose quadrupole moment changes in time radiates gravitational waves — ripples in spacetime that propagate at the speed of light and carry energy. A single mass moving uniformly radiates nothing; a spherically symmetric pulsation radiates nothing; but two masses whirling around each other have a quadrupole moment that oscillates, so they must radiate.

That radiation has to be paid for. The energy carried off by the waves comes out of the orbit's own energy budget. For a bound orbit, total energy is negative and scales as

E = − G m₁ m₂ / (2a)

where a is the semi-major axis. Losing energy (making E more negative) means a must shrink. By Kepler's third law the orbital period P scales as a^(3/2), so a smaller orbit means a shorter period. The binary therefore speeds up and tightens as it radiates — the orbit decays. This is the unmistakable signature: not a wave detected directly, but the recoil of a system paying its gravitational-wave bill.

The quadrupole formula and the decay rate

The rate at which a binary radiates gravitational-wave energy comes from the quadrupole formula Einstein derived in 1918. For two point masses on an eccentric orbit, Peters and Mathews (1963) integrated it to give the average rate of change of the orbital period:

dP_b/dt = − (192π / 5) (2π G / c³)^(5/3)
            × (P_b / 2π)^(−5/3)
            × [ m₁ m₂ / (m₁ + m₂)^(1/3) ]
            × f(e)

f(e) = (1 + 73/24 e² + 37/96 e⁴) / (1 − e²)^(7/2)

Two things in this formula matter enormously for the Hulse-Taylor binary. First, the decay scales steeply with how compact the orbit is: the shorter the period P_b, the faster the inspiral. A 7.75-hour orbit is tight enough for the effect to be measurable in years. Second, the eccentricity enhancement factor f(e) is large here. With e = 0.617, f(e) ≈ 11.8 — the binary radiates almost twelve times faster than a circular orbit of the same period would, because most of the radiation is emitted in the brief, fast, close passage through periastron.

Plugging in the measured masses and orbit, general relativity predicts

Ṗ_b (GR)       = − 2.40 × 10⁻¹² s/s
Ṗ_b (observed) = − 2.40 × 10⁻¹² s/s   (after a small galactic correction)

The ratio of observed to predicted decay is 0.997 ± 0.002 — agreement at the level of about 0.2 percent. That single number is the proof.

How the measurement is actually made

Astronomers never see the orbit directly; the system is an unresolved point 21,000 light-years away. Everything is reconstructed from the arrival times of the pulsar's 59-millisecond pulses, a technique called pulsar timing. The model that is fit to the data includes:

• Keplerian orbit. Five classical parameters — orbital period, eccentricity, semi-major axis projected along the line of sight, time and longitude of periastron — set the basic light-travel-time delay across the orbit (the Rømer delay).
• Periastron advance (ω̇). The orbit's ellipse slowly rotates in its plane at 4.226 degrees per year. This relativistic precession, the same effect as Mercury's perihelion shift but tens of thousands of times larger, is measured cleanly and constrains the total mass.
• Time dilation and gravitational redshift (γ). The pulsar's clock runs slow when it is deepest in the companion's gravity well and moving fastest at periastron — a combined special- and general-relativistic effect that modulates the pulse arrival times.
• Shapiro delay. Pulses that pass close to the companion are delayed by the curvature of spacetime around it, giving an independent handle on the companion's mass and orbital inclination.
• Orbital decay (Ṗ_b). Over many orbits, the steady shrinking of the period shows up as a quadratic drift in the time of periastron — the cumulative effect grows as t², making it detectable after a few years and unmistakable after decades.

Because the periastron advance and the redshift term together fix both neutron-star masses with no other assumptions, the decay rate becomes a clean, parameter-free prediction of general relativity. The observation then either confirms or refutes it. It confirmed it.

The system by the numbers

Quantity	Value	Note
Pulsar spin period	59.03 ms	The timing clock
Orbital period P_b	7.7519 hours	≈ 27,907 s
Eccentricity e	0.6171	Highly elliptical
Pulsar mass m₁	1.438 M☉	From timing
Companion mass m₂	1.390 M☉	Unseen neutron star
Semi-major axis	≈ 1.95 × 10⁶ km	≈ 2.8 R☉; comparable to the Sun's radius × a few
Separation at periastron	≈ 1.1 R☉	Stars nearly graze a solar diameter apart
Periastron advance ω̇	4.226 °/yr	vs Mercury's 43″/century
Orbital decay Ṗ_b	−2.40 × 10⁻¹² s/s	≈ −76 µs/yr
Orbital shrinkage	≈ 3.5 m/yr	The orbit tightens metres per year
Distance	≈ 21,000 ly (6.4 kpc)	In Aquila
Time to merger	≈ 300 Myr	Future neutron-star coalescence

Two of these figures deserve a second look. The orbit is shrinking by only about three and a half metres per year — out of a semi-major axis of nearly two million kilometres — yet pulsar timing measures it. And the periastron advances at 4.2° per year: at that rate the orbit's ellipse makes a complete turn in only about 85 years, an entire rotation within a human lifetime — the same relativistic effect that Mercury, in the Sun's far weaker field, takes roughly three million years to complete.

Periastron advance: Mercury, scaled up

The slow rotation of an orbit's long axis is one of the oldest tests of general relativity. Newtonian gravity from two point masses produces a perfectly closed ellipse that never precesses. Einstein's theory adds a small correction that makes the ellipse rotate. For Mercury, deep in the Sun's relatively weak field, that rotation is a mere 43 arcseconds per century — the anomaly that GR famously explained in 1915.

The Hulse-Taylor pulsar orbits inside the gravitational field of a 1.4-solar-mass neutron star, with periastron separations of barely a solar radius and orbital speeds of a few hundred kilometres per second. The relativistic correction scales roughly as the orbital compactness GM/(ac²), so the precession is enormous by comparison: 4.226° per year, about 35,000 times the Mercury rate. Measuring ω̇ to many significant figures is straightforward here, and combined with the time-dilation term it nails down both masses — the keystone that turns the decay rate into a parameter-free GR prediction.

Hulse-Taylor versus a direct LIGO detection

Both the Hulse-Taylor binary and LIGO are gravitational-wave science, but they probe opposite ends of an inspiral. It is worth being precise about what each one actually measures.

Aspect	Hulse-Taylor binary	LIGO direct detection
What is measured	Orbital decay (energy loss)	The passing wave itself (strain)
Evidence type	Indirect — recoil on the source	Direct — spacetime distortion at Earth
Stage of inspiral	Early, wide, slow (hours-long orbit)	Final seconds, merger
GW frequency	≈ 70 µHz (10⁻⁵ Hz)	≈ 35–250 Hz (audio band)
Timescale of measurement	Years to decades of timing	A fraction of a second
First result	1979–1982 (decay confirmed)	2015 (GW150914)
Instrument	Arecibo radio telescope	Laser interferometers
Recognised by	Nobel Prize 1993	Nobel Prize 2017

The Hulse-Taylor binary radiates at about 70 microhertz — far below any detector's reach, and it will not enter LIGO's audio band until its final few minutes, hundreds of millions of years from now. So the two are complementary: the binary pulsar proved gravitational waves carry energy by watching the source lose it; LIGO proved the waves are real by catching one in the act. Both got Nobel Prizes, 24 years apart, for the two halves of the same statement.

Why it matters and what came after

• First proof gravitational waves carry energy. Before 1974 it was genuinely unsettled whether gravitational waves were physical or a mathematical artefact of coordinate choices — Einstein himself wavered, and the question was debated into the 1950s. The Hulse-Taylor decay settled it: the waves remove energy at exactly the predicted rate.
• The 1993 Nobel Prize. Hulse and Taylor shared the prize for the discovery and its gravitational implications. The cumulative-periastron-shift plot — observed points falling on the GR parabola across decades — became one of the most cited figures in modern physics.
• A growing family of relativistic binaries. The Double Pulsar PSR J0737−3039, found in 2003, contains two visible pulsars and tests GR to even higher precision, confirming the decay prediction to better than 0.05 percent and adding tests like geodetic spin precession.
• The bridge to LIGO and GW170817. Every binary-neutron-star inspiral LIGO and Virgo detect is a Hulse-Taylor system caught in its final seconds. GW170817 (2017) was the merger of just such a pair, complete with an electromagnetic kilonova counterpart — the endpoint that PSR B1913+16 is slowly heading toward.
• Pulsar timing as a tool. The same timing precision now underpins pulsar timing arrays hunting for nanohertz gravitational waves from supermassive black-hole binaries across the whole sky.

Common misconceptions and edge cases

• "They detected gravitational waves in 1974." No — they detected a pulsar in a binary in 1974. The gravitational-wave evidence came from years of subsequent timing (the decay was confirmed around 1979–1982). And it is still indirect: the waves themselves were never measured, only their back-reaction on the orbit.
• "The decay rate is constant." It accelerates. As the orbit shrinks, P_b falls, and dP_b/dt scales as P_b^(−5/3), so the inspiral runs away — slowly now, catastrophically fast in the final minutes before merger.
• "You can just read the period change off directly." The raw decay is tiny — about 76 µs/yr on a 27,907-second period. What is actually fit is the cumulative shift in periastron time, which grows quadratically with elapsed time — about 40 seconds after the first thirty years, and roughly a hundred after fifty. Over short spans it is invisible; the proof lives in the long baseline.
• "A galactic correction is negligible." It is not. The Solar System and the binary accelerate differently in the Milky Way's gravitational field, producing a kinematic (Shklovskii-type) contribution to the apparent Ṗ_b. Subtracting it — using the system's distance and proper motion — is essential; getting the distance wrong is the dominant remaining uncertainty in the 0.2 percent agreement.
• "Both stars are visible pulsars." Only one is. The companion is inferred to be a neutron star from its mass and the system's history, but its radio beam (if any) does not sweep Earth, so it is detected only through its gravitational effect on the pulsar's timing.
• "The orbit is circular like most binaries you picture." It is strongly eccentric, e = 0.617. That eccentricity is what makes the system such a good laboratory — it amplifies both the periastron advance and the gravitational-wave luminosity, and it lets the redshift term be disentangled from the others.