Binary Pulsars and Relativistic Timing

Binary pulsars provide a rare combination: neutron stars have strong internal self-gravity, yet their orbital motion and radio pulses can be modeled phase coherently for many years. The resulting tests are not measurements of “strong gravity” in the abstract. They are comparisons between overdetermined timing observables and one theory’s predictions after propagation, clocks, astrometry, kinematics, Galactic acceleration, and stochastic noise have been controlled. This chapter builds that chain from residuals to individual and population constraints. Empirical assessments are current through 2026-08-19.

Helpful background. Weak-Field Gravity supplies barycentric propagation and external accelerations; Probability and Statistical Inference supplies likelihoods and covariance; Relativistic Two-Body Dynamics supplies conservative compact-binary motion; and Gravitational Radiation and Waveforms supplies radiative balance and multipolar emission.

Enter the timing experiment

The observable begins as a pulse time of arrival referred to an observatory clock. It becomes a barycentric timing datum only after clock transfer, station motion, a Solar-System ephemeris, relativistic propagation, and interstellar dispersion are applied. A phase-coherent spin and binary model predicts the pulse number; the residual carries the difference. Annual, chromatic, common-mode, and orbital-phase structures help identify which part of the model is wrong.

The next layer is post-Keplerian overdetermination. Keplerian parameters describe the leading orbit, while periastron advance, Einstein delay, Shapiro delay, orbital decay, and related relativistic parameters are fitted timing coordinates. In a declared theory, two independent parameters can determine the component masses. Additional parameters must agree at that same mass pair. This is why a mass–mass intersection is a theory test, whereas fitting a fresh pair for each observable is not.

Radiative tests require one more separation:

$$\begin{aligned} \dot P_b^{\rm obs}&=\dot P_b^{\rm int}+\dot P_b^{\rm Shk}+\dot P_b^{\rm Gal},\\ \dot P_b^{\rm int}&=\dot P_b^{\rm rad}+\dot P_b^{\dot m}+\dot P_b^{\rm tide}+\cdots . \end{aligned}$$

The Shklovskii and differential Galactic terms are apparent accelerations; mass loss and tides can also change the orbit. Only the corrected intrinsic derivative can be compared with quadrupole or extra-field radiation. The Double Pulsar analysis is an exemplar because it carries these corrections and higher-order timing terms through a joint model Kramer et al. 2021, §§III–VI, Table IV, and Eqs. (39)–(48).

Route through the chapter

Read the six pages in order when learning the inference chain:

1. Pulsar Timing Models and Residuals maps profiles and timestamps into structured residuals, then distinguishes spin, astrometric, chromatic, clock, ephemeris, noise, and binary signatures.
2. Post-Keplerian Parameter Tests derives the mass function and the logic by which two theory curves determine masses while further curves test consistency.
3. Orbital Decay and Radiative Damping separates the intrinsic period derivative from Shklovskii, Galactic, distance, and mass-loss corrections before comparing radiation models.
4. Dipole Radiation and Additional Fields explains why unequal compact-body charges generate a minus-one-PN channel and why equal charges can hide it.
5. Triple Systems and Strong Equivalence uses the phase-locked polarization of an inner orbit to test differential compact-body acceleration in an external field.
6. Population Timing and Theory Constraints combines systems through a selection-normalized hierarchy while retaining object-dependent charges, noise, and shared calibration.

The durable timing and strong-field formalism is reviewed by Freire and Wex 2024, §3.4.2, §3.5.1, and §§4–5. Numerical bounds should always be read with the named observing span, ephemeris, distance information, timing model, and theory map that produced them.

What each null result means

A precise GR-consistent PK intersection supports the whole declared timing-plus-GR model for that system. It does not prove that every gravity theory reduces to GR. Dipole radiation can vanish when charges are equal; a massive field may not radiate at the orbital frequency; screening can depend on environment; and scalarization can turn on only in a restricted mass range. Conversely, a residual departure does not identify new gravity until pulse-profile evolution, plasma, red noise, clock, ephemeris, astrometric, and binary-model alternatives fail discriminating checks.

The chapter therefore uses three recurring controls. First, predictions and nuisance corrections are fit with covariance, not subtracted as independent rounded numbers. Second, a proposed gravity term must keep its distinctive dependence on frequency, orbital phase, eccentricity, mass, or outer-orbit phase. Third, a population combination maps one theory consistently through every star and divides out interim priors before multiplying information. These controls make a negative result scientifically meaningful within its tested domain and make a possible anomaly difficult to manufacture through flexibility.

Review the chapter

1. Overdetermination. A timing solution measures four PK parameters but uses two to infer masses. What role do the other two play?

Solution

In a fixed theory the first two parameters select a common mass pair. The other two are then predictions evaluated at that same pair, so their residuals test the theory rather than supply permission to refit different masses. A joint likelihood must retain all parameter covariance and shared nuisance quantities. If the additional curves miss the common intersection beyond calibrated uncertainty, the complete timing-plus-theory model is inconsistent in that system.

2. Dipole blind spot. Why do asymmetric neutron-star–white-dwarf binaries and double-neutron-star binaries provide complementary information?

Solution

The leading extra-field dipole scales with the squared difference of effective charges. A neutron star and weakly bound white dwarf can have very unequal charges, amplifying that channel, whereas two neutron stars may acquire similar charges and suppress it. Double-neutron-star systems nevertheless supply exceptionally precise multi-parameter conservative and quadrupolar tests. Comparing both configurations therefore probes distinct combinations and reduces, but does not eliminate, theory-dependent blind regions.

References

• Freire, P. C. C., and Wex, N. (2024), “Gravity Experiments with Radio Pulsars,” Living Reviews in Relativity 27, article 5. DOI.
• Kramer, M., et al. (2021), “Strong-Field Gravity Tests with the Double Pulsar,” Physical Review X 11, 041050. DOI.

Further reading

• Archibald, A. M., et al. (2018), “Universality of free fall from the orbital motion of a pulsar in a stellar triple system,” Nature 559, 73–76. DOI.