Binary systems of neutron stars are ideal laboratories for testing fundamental laws of gravity in the strong field. When one of the components is a pulsar emitting strictly periodic radio pulses, it can be used as an ultra-precise clock. By measuring shifts in pulse arrival times caused by motion in the curved spacetime of the companion, astrophysicists reconstruct the relativistic dynamics of the system. It was in this way in the 1970s that Russell Hulse and Joseph Taylor proved the existence of gravitational waves by the orbital decay rate of a pulsar. Today, thanks to sensitive next-generation radio telescopes, such tests are becoming routine, but only a few systems are as extreme as PSR J1757–1854.
The key to success was the combination of archival observations from the 64-meter Parkes telescope and the 100-meter Green Bank Telescope with new data from the South African MeerKAT array. MeerKAT provided record sensitivity in two frequency bands, partially mitigating radio-wave scattering in the turbulent interstellar medium. The total time span covered 9 years. The methodology involved precision timing using the relativistic DDH model, where the orbit is described not only by Kepler's laws but also by a set of post-Keplerian parameters—observable manifestations of Einstein's field equations. Special attention was paid to joint sampling of the parameters γ (Einstein delay, related to time dilation) and δθ. For this, Bayesian analysis with an MCMC algorithm was applied to properly account for their strong correlation.
The main result is a confident (~7σ at fixed γ) detection of the relativistic angular orbital deformation parameter δθ. This effect, a tiny correction to the classical ellipse on the order of 5×10⁻⁶, had previously been measured in only two binary pulsars: the legendary B1913+16 (requiring 40 years of observations) and the double pulsar J0737–3039. PSR J1757–1854 revealed it in just 9 years—a testament to the system’s phenomenal compactness. The value of δθ allowed a test of general relativity: subtracting the aberration contribution, dependent on the pulsar spin axis orientation, we obtained the intrinsic relativistic deformation, consistent with the GR prediction for only two of the four previously admissible geometric configurations. Thus, two scenarios requiring an exotic spin axis tilt were ruled out. Moreover, the precision measurement of periastron advance rate (¤ω = 10.364986(8) degrees per year) revealed contributions from second post-Newtonian order and the Lense-Thirring effect—the dragging of the orbit by the neutron star’s rotation. These corrections, though small (~0.0003°/yr), already noticeably shift the total mass estimate of the system.
The result turns PSR J1757–1854 into a unique tool for relativistic astrophysics. The measurement of δθ opened a new independent method for determining the pulsar spin axis orientation, which is critical for reconstructing the system’s birth scenario: the kick of the supernova that produced the second companion must have imparted a significant tilt to the orbit. The agreement of the observed gravitational-wave orbital decay rate (after subtracting kinematic corrections) with the prediction of the quadrupole formula confirms GR at a level limited only by the uncertainty in the pulsar distance. Finally, the detection of the Lense-Thirring effect contribution to periastron advance is a step toward directly measuring the neutron star’s moment of inertia.
Continued long-term monitoring of PSR J1757–1854 with MeerKAT and the upcoming SKA promises revolutionary breakthroughs. As the pulsar spin axis precesses (geodetic precession), the visible pulse shape will change, and the periastron's approach to the line of sight will improve the measurement accuracy of the Shapiro delay. By 2031, it is likely that secular change in the projected semi-major axis (¤x) caused by the Lense-Thirring effect will be detected, which, combined with geometry data, will allow the moment of inertia to be calculated. This will turn the system into a laboratory for determining the stiffness of the nuclear matter equation of state, competitive with gravitational-wave detectors.
The results will impact strong-field gravity physics, compact object astrophysics, nuclear physics, and cosmology—refining the merger rate of double neutron stars for predicting population properties.
Immediate tasks: continuing high-cadence timing with MeerKAT, refining the pulsar’s proper motion, and an independent distance measurement (e.g., via very long baseline radio interferometry) to lift the uncertainty in kinematic corrections.
The work directly connects to several key problems: the equation of state of ultra-dense matter (moment of inertia as a stiff discriminator of models); testing GR in the strong field (observing higher-order expansion terms); physics of supernova explosions (pulsar spin axis orientation as a diagnostic of the collapse kick).
🎯 The orbit of PSR J1757–1854 shrinks by 5.282 femtoseconds per second due to gravitational wave emission—that accumulates to about 0.17 milliseconds per year. In 76 million years, these two neutron stars will merge in a kilonova flash, similar to GW170817.