Binary pulsar

A binary pulsar is a pulsar in a binary star: the companion's gravity changes the emitting neutron star's motion and hence its pulse arrival times. The companion need not itself be a detectable pulsar. Pulsar timing turns that motion into a measurable orbit, making the system useful for stellar mass determination and tests of general relativity.

Russell Alan Hulse and Joseph Hooton Taylor Jr. discovered this binary pulsar at Arecibo in 1974; their January 1975 paper reports a roughly $59\mathrm{ms}$ pulse period and $7.75\mathrm{h}$ orbital period. Its unseen companion was not detected as a second radio pulsar. Discovery of the orbit preceded the later binary-pulsar orbital decay test.^[ref]

The 1993 Nobel Prize in Physics recognized the discovery.^[ref]

By 1993, the system was interpreted as a pair of neutron stars.^[ref]

The surprise was a pulsar whose apparent period would not stay fixed. Russell Alan Hulse's Arecibo search used a multichannel receiver, pulsar dedispersion and computer searches over period and pulse width. Follow-up at higher time resolution retained the changing period, excluding inadequate sampling as its explanation. The observed variation exceeded the expected terrestrial Doppler shift; a predicted reversal from decreasing to increasing period appeared on 16 September 1974.^[ref]

The 1975 analysis fitted an eccentric Kepler orbit to Doppler shifts inferred from about 200 five-minute period measurements over ten days. That discovery fit was not yet the later precision arrival-time experiment.^[ref] An orbit explains both signs of the apparent period change without requiring the stellar rotation to alternate between speeding up and slowing down.

The semimajor axis $a_1$ of the pulsar's orbit about the binary barycenter, projected along the line of sight and divided by the speed of light. The orbital inclination $i$ is zero for a face-on orbit. Thus $x$ is a time, often quoted as the numerical number of light-seconds in $a_1\sin i$, not the true binary separation. The relative-orbit semimajor axis is instead $a=a_1+a_2$; conservation of momentum gives $a_1=a{m}_2/(m_1+m_2)$.

Combining Kepler's third law with that barycenter relation gives the timing form of the binary mass function, also valid for eccentric Kepler orbits:Choose $x=2\mathrm{s}$ and $P_b=8\mathrm{h}$ for an original example: $a_1\sin i=599,585\mathrm{km}$ and $f\simeq 0.0773M_{\odot }$, where $M_{\odot }$ is a solar mass. For a specified pulsar mass $m_1=1.4M_{\odot }$, the smallest allowed companion mass occurs at $\sin i=1$ and is approximately $0.698M_{\odot }$. A different assumed $m_1$ changes this bound. One measured $f$ does not determine both masses and orbital inclination.

The geometric light-travel contribution from the pulsar's position within the binary orbit. To construct an eccentric Kepler orbit, use eccentric anomaly $E$, orbital eccentricity $e$, mean anomaly $M=2\pi (t-T_0)/P_b$, and argument of periapsis $\omega$. Here $T_0$ is passage through periapsis. Coordinates in the orbital plane are $X=a_1(\cos E-e)$ and $Y=a_1\sqrt{1-e^2}\sin E$. Rotate through $\omega$ and project by $\sin i$ to obtainThe first equation is Kepler's equation: equal time intervals give equal increments of $M$, not of $E$. Since $1-e\cos E>0$ for $0\le e<1$, each $M$ has a unique $E$ on the corresponding unwrapped orbit. A constant delay is absorbed into the arbitrary arrival-time origin. The circular limit is ${\Delta }_R=x\sin (M+\omega )$, so its amplitude measures $x$. For $e>0$, the delay's offset and amplitude also depend on $\omega$; its maximum is not generally $x$.

A gradual rotation of the direction of periapsis within the orbital plane. A closed Kepler orbit has fixed argument of periapsis $\omega$; relativistic gravity can make $\omega$ advance between successive orbits. The angle change and binary-pulsar orbital decay describe different changes of the orbit.

Decreasing orbital period, distinct from spin slowdown. Joseph Hooton Taylor Jr. inferred masses using binary-pulsar periastron advance and binary-pulsar Einstein delay, then tested the independent period derivative against the general relativity quadrupole formula. Correcting differential Galactic acceleration gave indirect gravitational radiation evidence. This was not a terrestrial strain measurement.^[ref]

In particular, the raw measured derivative schematically containsHere $a_{\parallel }^{\mathrm{relative}}$ is positive when the source's recession accelerates. Its contribution follows by differentiating the leading Doppler effect. A decreasing observed orbital period alone therefore does not identify a unique energy-loss mechanism.

For fixed masses, binary orbital energy is $-G{m}_1{m}_2/(2a)$ and Kepler's third law gives $P_b\propto a^{3/2}$. Positive outgoing gravitational radiation predicted by the quadrupole formula makes the orbital energy more negative, reducing $a$ and $P_b$. This sign argument explains shrinking orbits, but a quantitative test also needs the predicted magnitude and orbital eccentricity dependence.

A tiny orbital period change accumulates into a larger timing shift. Choose the original model $P_b(t)=P_{b0}+{\dot{P}}_{b}t$, with $P_{b0}=28,800\mathrm{s}$ and ${\dot{P}}_b=-2\times 10^{-12}s/s$. Expanding the accumulated cycle count givesNegative $\delta T$ means an event of fixed orbital cycle count arrives earlier than a constant-period prediction. After twenty Julian years, the chosen period has decreased by $1.262\mathrm{ms}$, while the accumulated event shift is about $-13.83\mathrm{s}$. This is a time shift of periapsis passage; it is different from a rotation of the direction of periapsis.