Pith. Sign in

REVIEW

Testing General Relativity with Pulsar Timing

Astrophysics · Submitted 30 Jul 2003 · 153 references

Pith's verdict: A competent, honest review of pulsar gravity tests; the only real soft spot is that the B1913+16 Ṗb correction budget is named but never quantified.

T0 · AI review

How this review was made: An AI referee (selfhosted/deepseek-v4.1-flash) read the paper on October 29, 2026 and answered a fixed list of questions about it. T0 is the starting level every paper gets. The levels describe how deeply Pith has checked the paper, never the authors.

  1. T0 AI review An AI read the paper and answered a fixed list of questions.
  2. T1 Audit Automatic checks that the references and files the paper relies on exist and match.
  3. T2 Compute The paper's code or calculations re-run with the same results.
  4. T3 Formal statement The claims written out in a formal math language a computer can read.
  5. T4 Proof Those claims proved, and the proof checked by computer.

How Pith checks papers →

· 1 major objection · 12 minor · Ask about this paper · The compiled record · Read the full review →
Challenge this review Comment (0) Read on arXiv ↗
Cite ▾
arXiv IDastro-ph/0307536 (version 1)
Pith ID2P7O75A2, Pith's permanent record of this paper and its review
CitationIngrid H. Stairs (UBC). "Testing General Relativity with Pulsar Timing." arXiv:astro-ph/0307536 (2003). https://pith.science/paper/2P7O75A2
Cite this reviewBibTeX for Pith's review of this paper
Claim paper

Why Pith says this

Highlight any text to comment on it.

Referee report
The paper

The paper claims

Binary pulsars turn their orbits into over-determined tests of general relativity, and the measured relativistic timing parameters agree with the theory's mass predictions to better than one percent.

The reading

A pulsar is a spinning neutron star whose radio beam sweeps past Earth once per rotation, so every arriving pulse marks a tick of an extraordinarily steady clock. Put that clock in a tight orbit around another star and the arrival times pick up the relativistic corrections to the orbit: the precession of periastron, the combined gravitational redshift and time dilation, the slow shrinking of the orbit as gravitational waves carry energy away, and the Shapiro delay suffered by pulses passing the companion. Stairs' argument is that each of these corrections, written as a post-Keplerian timing parameter, is predicted by general relativity as a function of just the two stellar masses, so measuring two parameters fixes the masses and measuring three or more over-determines the system, which makes the agreement a test of the theory rather than a fit to it. The review then collects what pulsar timing says about equivalence-principle violations: differential acceleration of neutron stars and white dwarfs in the Galactic field, preferred-frame effects, dipolar gravitational radiation, and a possible drift of Newton's constant.

What carries the argument

The load-bearing device is the set of post-Keplerian timing parameters: the periastron advance rate $\dot\omega$, the Einstein delay $\gamma$, the orbital period derivative $\dot P_b$, and the range $r$ and shape $s$ of the Shapiro delay. In general relativity each is a closed algebraic function of the two stellar masses and the Keplerian orbital elements, so plotting the mass pairs allowed by each measured parameter on a mass-mass diagram turns the test into a geometric question: do the curves intersect at a common point? Two parameters give a solution; a third or fourth makes the intersection an over-determined check, and a mismatch of the curves rules the theory out. The second device is the forced-eccentricity argument for near-circular pulsar-white-dwarf binaries, in which a violation of the strong equivalence principle or a preferred-frame effect adds a constant vector to the periastron-rotating eccentricity vector, so an upper limit on the observed eccentricity becomes a limit on the violation.

Extended reading notes

Core claim

The paper's central claim is that binary pulsars, above all the double-neutron-star systems PSR B1913+16 and PSR B1534+12, are genuine strong-field laboratories for general relativity rather than merely confirmations of a formula. For PSR B1913+16 three post-Keplerian parameters are measured, the advance of periastron $&#92;dot&#92;omega$, the redshift-plus-dilation term $&#92;gamma$ and the orbital period derivative $&#92;dot P_b$, and the measured values agree with those predicted from the masses fitted under the general-relativistic timing model to better than 0.5 percent. For PSR B1534+12 the shape $s$ of the Shapiro delay supplies a fourth parameter, and the combination of $&#92;dot&#92;omega$, $&#92;gamma$ and $s$ alone tests the theory to better than 1 percent using only quasi-static effects, which separates the quasi-static sector from the radiative one tested by $&#92;dot P_b$. In parallel, an ensemble of pulsar-white-dwarf binaries with nearly circular orbits converts upper limits on orbital eccentricity into bounds on strong-equivalence-principle violation, $&#92;Delta_{&#92;rm net} < 9&#92;times 10^{-3}$ at 95 percent confidence, and on the preferred-frame parameter $&#92;hat&#92;alpha_1 < 1.4&#92;times 10^{-4}$, while single and binary millisecond pulsars bound $&#92;hat&#92;alpha_3$ and $&#92;zeta_2$, short-period pulsar-white-dwarf pairs bound dipolar emission through $(&#92;alpha_{c1}-&#92;alpha_0)^2$, and the observed neutron-star mass distribution bounds $&#92;dot G/G$ at $(-0.6&#92;pm 4.2)&#92;times 10^{-12}$ yr$^{-1}$.

Load-bearing premise

Comparing the measured orbital period derivative with the quadrupole prediction requires subtracting the Galactic acceleration of the binary and the Shklovskii proper-motion term, and for PSR B1534+12 the distance needed for that subtraction is itself obtained by assuming general relativity, so the correction and the test lean on each other.

Editorial extensions

If this is right

  • If the over-determination holds, the quadrupole formula for gravitational radiation is confirmed in strong fields at the sub-percent level, and any alternative theory that assigns a different mass dependence to the post-Keplerian parameters is excluded by the same data.
  • Nearly circular pulsar-white-dwarf binaries bound strong-equivalence-principle violation through the forced eccentricity, and the bound improves directly as long-term timing tightens the eccentricity limits.
  • Short-period pulsar-white-dwarf pairs such as PSR B0655+64 and PSR J1012+5307 constrain scalar coupling through dipolar radiation, while the observed orbital decay of double-neutron-star systems rules out theories whose quadrupole emission predicts a positive period derivative.
  • The narrow spread of measured neutron-star masses near the Chandrasekhar mass converts into a limit on a varying Newton's constant, $&#92;dot G/G = (-0.6&#92;pm 4.2)&#92;times 10^{-12}$ yr$^{-1}$.
  • A pulsar-black-hole binary with a roughly ten-solar-mass companion and an eight-hour orbit would tighten the scalar coupling $&#92;alpha_0$ by about a factor of fifty relative to current solar-system tests.

Reading between the lines

Editorial extensions of the paper, not claims the author makes directly.

  • The over-determination is a consistency check within a fixed two-body, point-mass timing model; a classical contribution that happened to mimic the mass dependence of a post-Keplerian parameter, such as a companion mass quadrupole contributing to $&#92;dot&#92;omega$, would pass the test unnoticed, which is why the paper flags this possibility for PSR J1141$-$6545.
  • Because the distance to PSR B1534+12 is currently obtained by assuming general relativity and inverting the $&#92;dot P_b$ discrepancy, an independent VLBI parallax is the single observation that would convert that system from a consistency check into a clean test.
  • The same forced-eccentricity machinery could be inverted and applied as a stack: cross-correlating the residual eccentricity vectors of many pulsar-white-dwarf binaries against the Cosmic Microwave Background dipole direction would search for a common preferred frame rather than bounding each system separately.
  • Geodetic precession is a fifth, qualitatively different probe of the strong field, but it is currently limited by beam-shape modelling; a pulsar whose emission geometry is pinned down polarimetrically across a full precession cycle would turn the precession rate into a model-independent test.
Original abstract

Pulsars of very different types - isolated objects, and binaries with short- and long-period orbits, white-dwarf and neutron-star companions - provide the means to test both the predictions of general relativity and the viability of alternate theories of gravity. This article presents an overview of pulsars, then discusses the current status and future prospects of tests of equivalence principle violations and strong-field gravitational experiments.

The compiled record

The paper's own sentences, arranged, kept beside the review and separate from it. Each label is drawn from a closed vocabulary; the number beside it is the probability the classifier gave that label. Numbers in the tables are the authors', checked by code. No sentence in this record was written by a model.

Compiled by classifier-1.13.0 via TypeSafe API, question bank abstract-p1-20260923, 2026-09-23: 1 request, 2,329 input tokens.

What the abstract itself states, as probabilities: claims first 0.02 · states limitation 0.06 · human participants 0.02 · controlled comparison 0.04 · randomized assignment 0.01 · states quantitative result 0.03 · compares with prior results 0.06 · mentions code or data availability 0.02.

The abstract, sentence by sentence

arXiv listing abstract v1. Results and the main claim are underlined; .

Pulsars of very different types - isolated objects, and binaries with short- and long-period orbits, white-dwarf and neutron-star companions - provide the means to test both the predictions of general relativity and the viability of alternate theories of gravity. This article presents an overview of pulsars, then discusses the current status and future prospects of tests of equivalence principle violations and strong-field gravitational experiments.

background contribution method result interpretation main claim

What the classifier reports for each sentence. The number is the probability it gave the label.
#sentencerolethe authors' own workstates the main conclusionstrength
s0 Pulsars of very different types - isolated objects, and binaries with short- and long-period orbits, white-dwarf and neutron-star companions - provide the means to test both the predictions of general relativity and the viability of alternate theories of gravity. background 0.95 0.07 0.15 asserted 0.98
s1 This article presents an overview of pulsars, then discusses the current status and future prospects of tests of equivalence principle violations and strong-field gravitational experiments. contribution 1.00 0.26 0.11 asserted 0.86

What kind of paper the classifier reports this is: review 1.00 · evidence from literature 0.95.

What was examined, and what was not

sources read
  • arXiv listing abstract v1
not examined
  • full text
  • figures and tables
  • references and citing sentences
  • relations to other papers
levels run
  • P abstract
not run
  • S sections
  • T body sentences
  • C citations
  • R cross-paper relations

Every quoted sentence on this record is the authors'. Every label is from Pith's closed vocabulary, with the probability shown beside it. Every number in the measurements table is the authors', checked by code. Record as JSON

Figures
Figure 1. Top: 100 single pulses from the 253-ms pulsar B0950+08, demon￾strating pulse-to-pulse variability in shape and intensity. Bottom: Cumulative profile for this pulsar over 5 minutes (about 1200 pulses); this approaches the reproducible standard profile. Observations taken with the Green Bank Telescope [1]. Stairs, unpublished. The exact mechanism by which a pulsar radiates the energy observed as radio pulses is still a subject of vigorous debate. The basic picture of a mis￾aligned magnetic dipole, with coherent radiation from charged particles acceler￾ated along the open field lines…
Figure 2. Pulse profile shapes for PSR J1740−3052 at multiple frequencies, aligned by pulse timing. The full pulse period is displayed at each frequency. The growth of an exponential scattering tail at low frequencies is evident. All observations taken with the Green Bank Telescope [1] (Stairs, unpublished), ex￾cept for the 660-MHz profile which was acquired at the Parkes telescope [2, 123]. A short description of pulsar observing techniques is in order. As pulsars have quite steep radio spectra (e.g.,[98]), they are strongest at frequencies f0 of a few hundred MHz. At these frequencies, the…
Figure 3. Pulse profile of the fastest rotating pulsar, PSR B1937+21, observed with the 76-m Lovell telescope at Jodrell Bank Observatory [3]. The top panel shows the total-intensity profile derived from a filterbank observation (see text); the true profile shape is convolved with the response of the channel filters. The lower panel shows the full-Stokes observation with a coherent dedispersion in￾strument [125, 124]. Total intensity is indicated by black lines, and linear and circular power by red and blue lines respectively. The position angle of the linear polarization is plotted twice.…
Figure 4. “Polarization” of a nearly circular binary orbit under the influence of a forcing vector g, showing the relation between the forced eccentricity, eF , the eccentricity evolving under the general-relativistic advance of periastron, eR(t), and the angle θ. After [145]. Converting equation 15 to a limit on ∆net requires some statistical arguments to deal with the unknowns in the problem. First is the actual component of the observed eccentricity vector (or upper limit) along a given direction. Damour and Sch¨afer (1991) [49] assume the worst case of possible cancellation between the…
Figure 5. Measured neutron star masses as a function of age. The solid lines show predicted changes in the average neutron star mass corresponding to hypo￾thetical variations in G, where ζ−12 = 10 implies G/G ˙ = 10×10−12 yr−1 . From [135], used by permission. The Chandrasekhar mass, MCH, is the maximum mass possible for a white dwarf supported against gravitational collapse by electron degeneracy pressure [37]. Its value — about 1.4 M⊙ — comes directly from Newton’s constant: MCH ∼ (¯h c/G) 3/2/m2 N, where ¯h is Planck’s constant and mN is the neutron mass. All measured and constrained…
Figure 6. The parabola indicates the predicted accumulated shift in the time of periastron for PSR B1913+16, caused by the decay of the orbit. The measured values of the epoch of periastron are indicated by the data points. From [144], courtesy Joel Weisberg. is clear that the curves do intersect, at the point derived from the DDGR mass predictions. Clearly any theory of gravity which does not pass such a self-consistency test can be ruled out. However, it is possible to construct alternate theories of gravity which, while producing very different curves in the m1 − m2 plane, do pass the PSR B
Figure 7. Mass-mass diagram for the PSR B1913+16 system, using the ω˙ and γ parameters listed in
Figure 8. Mass-mass diagram for the PSR B1534+12 system. Labeled curves illustrate 68% confidence ranges of the DD parameters listed in
Figure 9. Portions of the tensor–biscalar β ′ − β ′′ plane permitted by timing observations of PSRs B1913+16, B1534+12 and B1855+09 up to 1992. Values lying above the curve labeled “a” are incompatible with the measured ω˙ and γ parameters for PSR B1913+16. The curves labeled “b” and “d” give the allowed ranges of β ′ and β ′′ for PSRs B1913+16 and B1534+12, respectively, fitting for the two neutron-star masses as well as β ′ and β ′′, using data available up to 1992. The vertical lines labeled “c” represent limits on β ′ from the SEP-violation test using PSR B1855+09 [49]. The dot at (0,0)…
Figure 10. The parameter space in the non-linear α0, β0 gravitational theory, for neutron stars described by a polytrope equation of state. The regions below the various curves are allowed by various pulsar timing limits, by solar-system tests (“1PN”) and by projected LIGO/VIRGO observations of NS-NS and NS￾BH inspiral events. The shaded region is allowed by all tests. The plane and limits are symmetric about α0 = 0. From [44]; used by permission. A different class of theories, allowing a non-linear coupling between matter and a scalar field, was later studied by Damour and Esposito-Far`ese…
Figure 11. Solid line: predicted value of the Shapiro delay in PSR J0437−4715 as a function of orbital phase, based on the observed inclination angle of 42◦±9 ◦ . For such low-eccentricity binaries, much of the Shapiro delay can be absorbed into the orbital Roemer delay; what remains is the ∼ Pb/3 periodicity shown. The points represent the timing residuals for the pulsar, binned in orbital phase, and in clear agreement with the shape predicted from the inclination angle. Reprinted by permission from Nature [139], copyright 2001, Macmillan Pub￾lishers Ltd. A different and complementary test of
Figure 12. Changes in the observed pulse profile of PSR B1913+16 throughout the 1980s, due to a changing line-of-sight cut through the emission region of the pulsar. Taken from reference [133]; used by permission. A complete discussion of GR effects in pulsar observations must mention geodetic precession, though these results are somewhat qualitative and do not (yet) provide a model-free test of GR. In standard evolutionary scenarios for double-neutron-star binaries (e.g., [29, 111]), both stellar spins are expected to be aligned with the orbital angular momentum just before the second…
Figure 13. Top: change in peak separation of the relativistic double-neutron￾star binary PSR B1913+16, as observed with the Arecibo (solid points, [141]) and Effelsberg (open circles, [85]) telescopes. Bottom: projected disappearance of PSR B1913+16 in approximately 2025. Taken from reference [85]; used by permission. magnitude is given by: Ω spin 1 = 1 2  Pb 2π −5/3 m2(4m1 + 3m2) (1 − e 2)(m1 + m2) 4/3 T 2/3 ⊙ , (38) where T⊙ ≡ GM⊙/c3 = 4.925490947 µs, as in Section 4.1. This predicted rate of precession is small; the three systems with the highest Ωspin 1 values are: PSR J11
Figure 14. Hourglass-shaped beam for PSR B1913+16 derived from the symmetric-component analysis of [143]. Taken from reference [143]; used by permission. in the 1980s [141], with a clear change in the relative heights of the two pro￾file peaks over several years (
Figure 15. Alternate proposed beam shape for PSR B1913+16, consisting of a symmetric cone plus an offset core. The red lines indicate an example cut through the emission region, as well as the predicted pulse peak ratio and sepa￾ration as functions of time. After [86], courtesy Michael Kramer. profile changes requires some kind of model of the beam shape. Kramer [86, 87] lets the rate of precession vary as another free parameter in the pulse-shape fit, and finds a value of 1.2 ◦ ± 0.2 ◦ . This is consistent with the GR prediction but still depends on the beam-shape model and is therefore not…
Figure 16. Evolution of the low-level emission surrounding the main pulse of PSR B1534+12, over a period of nearly 10 years, as measured with the Arecibo telescope [5]. Stairs et al., unpublished. during the 1990s. Clearly, careful monitoring of this pulsar’s profile is in order. 44

Figures reproduced from arXiv:astro-ph/0307536 by the authors.

References
How the literature uses it

91 arXiv papers cite this in 116 sentences

Each citing sentence is taken from the citing paper's LaTeX source and typed by Pith's classifier. A label appears at probability 0.8 or more (0.95 for support), a contest at 0.7 or more; the quoted words are the citing authors'. How each label was checked · JSON

Papers citing it · 79 sentences without a confident label

Uses its data or values 3 papers
Compares with it 1 paper
Cites it as background 33 papers

Discussion

Disagree with the review, or know something it missed? Leave a comment, highlight any passage above to comment on it, or challenge this review formally.