Long-term timing of the relativistic binary PSR J1906+0746

This study presents an 18-year timing analysis of the relativistic binary PSR J1906+0746 using data from six radio telescopes, yielding precise measurements of its orbital parameters and component masses consistent with a double neutron star system under general relativity, while also characterizing a large glitch and reporting a potential secular change in the projected semi-major axis that could indicate a massive, fast-rotating white dwarf companion.

The Gist

Imagine the universe as a giant, cosmic dance floor. On this floor, two heavy partners are spinning around each other at incredible speeds. One partner is a pulsar—a dead, super-dense star (a neutron star) that acts like a lighthouse, beaming radio waves at us with perfect regularity. The other partner is its mysterious companion.

For over 18 years, astronomers have been watching this specific dance pair, known as PSR J1906+0746, using six giant radio telescopes scattered across the globe (including the famous Arecibo, the massive FAST in China, and others). They treated the pulsar's radio beeps like a cosmic metronome, listening to see if the rhythm changed.

Here is the story of what they found, broken down into simple concepts:

1. The Perfect Clock (and the Big Glitch)

Pulsars are usually the most reliable clocks in the universe. If you set your watch by them, it would never lose a second. However, this specific pulsar had a "tantrum."

Around 2014, the pulsar suddenly sped up. Think of it like a runner who, while jogging at a steady pace, suddenly gets a burst of energy and sprints forward, then slowly settles back into a new, slightly faster rhythm. This event is called a glitch.

• The Size: It was a huge glitch, comparable to the famous "Vela" pulsar's tantrums.
• The Aftermath: The pulsar didn't just snap back; it took about 100 days to calm down, and its "spin-down" rate (how fast it slows down naturally) changed permanently.

2. The Cosmic Wobble (Geodetic Precession)

This dance pair isn't just spinning; they are wobbling. Because they are so heavy and close together, their gravity twists the fabric of space-time around them (a concept from Einstein's General Relativity).

Imagine a spinning top. As it spins, its axis wobbles in a circle. This pulsar is doing the same thing, but on a cosmic scale. Because of this wobble, the "lighthouse beam" is slowly changing its angle relative to Earth.

• The Prediction: The astronomers calculated that this wobble is so extreme that the pulsar's beam will eventually swing away from Earth entirely. By 2028, the pulsar will effectively "go dark" to us, and we won't see it again until around 2070–2090 when the wobble brings it back into view.

3. The Mystery of the Dance Partner

The biggest question was: What is the other dancer?
Since the pulsar is young and hasn't been "recycled" (sped up by eating matter from a partner), it must be the second object to form in this system. That means the first object formed, died, and became the companion we see today.

There are two main theories for who this companion is:

• Theory A: A Twin Neutron Star. Maybe the companion is another neutron star. If so, it would have to be spinning incredibly fast (like a millisecond pulsar) to explain the data.
• Theory B: A Massive White Dwarf. Or, the companion is a "White Dwarf" (a dead star, but less dense than a neutron star) that is spinning very fast.

The Clue: The astronomers noticed a tiny, subtle change in the shape of the orbit over time. It's like noticing the dance floor itself is slowly tilting. This tilt is caused by the companion's spin dragging on the orbit (a phenomenon called spin-orbit coupling).

• The amount of tilt they measured is very similar to a famous system (PSR J1141-6545) where the companion is a fast-spinning White Dwarf.
• The Twist: If the companion is a White Dwarf, it means the "dance" started with the White Dwarf forming first, spinning up, and then the pulsar was born in a supernova explosion. This flips the usual story on its head!

4. Testing Einstein's Rules

Why do we care? Because this system is a perfect laboratory for testing General Relativity (Einstein's theory of gravity).

• The astronomers measured how the orbit shrinks over time. Einstein predicted that two heavy objects orbiting each other should lose energy by emitting gravitational waves (ripples in space-time), causing them to spiral closer together.
• Their measurements of the orbit shrinking match Einstein's prediction almost perfectly. This confirms that gravity works exactly as Einstein said it does, even in these extreme conditions.

5. The "Noise" in the Signal

The data wasn't perfectly clean. There was a "hum" or "noise" in the timing data that looked like a pattern repeating every 2 years.

• The team wondered: Could there be a third dancer? Maybe a planet orbiting the pair?
• They did the math, and while a planet could fit the data, it's more likely that this 2-year pattern is just weird behavior from the pulsar's own magnetic field or atmosphere. It's like hearing a rhythm in the wind and wondering if it's a drum, when it's actually just the trees swaying.

The Big Picture

This paper is a masterclass in patience and precision. By listening to a cosmic lighthouse for nearly two decades, using telescopes from the US, China, Europe, and South Africa, the team:

1. Confirmed Einstein's gravity is correct.
2. Caught a pulsar having a massive "glitch."
3. Predicted when this pulsar will disappear from our view.
4. Found strong evidence that its partner might be a fast-spinning White Dwarf, which would make this a rare and special "triple-generation" system.

It's a reminder that even in the cold, dead vacuum of space, stars are still dancing, wobbling, and telling us secrets about how the universe works—if we just listen closely enough.

Technical Summary

Here is a detailed technical summary of the paper "Long-term timing of the relativistic binary PSR J1906+0746":

1. Problem and Motivation

PSR J1906+0746 is a young (characteristic age $\tau_c \approx 112$ kyr), relativistic binary pulsar with a 114 ms spin period orbiting a compact companion every 3.98 hours. Previous timing analyses (e.g., van Leeuwen et al. 2015) established it as a double neutron star (DNS) candidate with precise measurements of post-Keplerian (PK) parameters. However, the system's nature remains ambiguous:

• Companion Identity: While consistent with a DNS system, the companion's mass and the system's evolutionary history leave open the possibility that it is a massive, fast-rotating white dwarf (WD) formed before the pulsar (similar to PSR J1141–6545).
• Timing Noise and Glitches: The pulsar exhibits significant timing noise and a large glitch, complicating long-term timing solutions and the extraction of precise orbital parameters.
• Precision Limits: Previous measurements of the orbital period derivative ($\dot{P}_b$) and the rate of change of the projected semi-major axis ($\dot{x}$) lacked the precision required to definitively distinguish between general relativity (GR) predictions and alternative scenarios involving spin-orbit coupling.

The primary goal of this study is to conduct a coherent 18-year timing analysis to refine PK parameters, measure $\dot{x}$, and constrain the nature of the companion and the validity of GR.

2. Methodology

The authors performed a comprehensive timing analysis using data from six radio telescopes (Arecibo, FAST, Green Bank, Lovell, MeerKAT, and Nançay) spanning from May 2005 to August 2023.

• Data Reduction:
  • Data were processed using various backends (WAPP, ASP, PUPPI, GASP, BON, AFB, DFB, PTUSE).
  • Radio Frequency Interference (RFI) was removed manually or via algorithms (PSRZAP, CLFD).
  • Pulse profiles were modeled using multiple templates to account for profile evolution caused by geodetic precession.
  • Times of Arrival (ToAs) were extracted using cross-correlation against noiseless templates.
• Timing Model:
  • Coherent Connection: The authors achieved phase coherence across the entire 18.2-year dataset, bridging a large glitch that occurred around MJD 56664.
  • Red Noise Modeling: Unlike previous studies that treated timing noise as arbitrary offsets, this work utilized RUN_ENTERPRISE (Keith et al. 2022) to model red noise as a Gaussian process with a power-law Fourier basis. This allowed for simultaneous fitting of timing noise and physical parameters.
  • Binary Models: The analysis employed the DD (Damour-Deruelle) model for a theory-independent fit and the DDGR model (assuming General Relativity) to derive masses.
  • Glitch Modeling: A specific glitch model was fitted to account for the sudden frequency increase and subsequent relaxation.
• Kinematic Corrections: The observed orbital period derivative ($\dot{P}_{b,obs}$) was corrected for Galactic acceleration and the Shklovskii effect (transverse Doppler shift) to derive the intrinsic orbital decay.

3. Key Contributions

• Extended Baseline: The timing baseline was extended by a factor of ~4.5 compared to previous work, covering 18 years of data.
• Advanced Noise Modeling: The application of a Fourier-basis Gaussian process for red noise significantly improved the precision of parameter estimation, particularly for $\dot{P}_b$ and $\dot{x}$.
• First Detection of $\dot{x}$: The study reports the first $3\sigma$detection of the secular change in the projected semi-major axis ($\dot{x}$) for this system.
• Glitch Characterization: Detailed characterization of a large glitch (MJD 56664) with a fractional frequency increase comparable to the Vela pulsar.

4. Key Results

Timing Parameters and Masses

• Post-Keplerian Parameters:
  • Advance of periastron: $\dot{\omega} = 7.5841(2)$ deg yr$^{-1}$ (2.5$\times$ more precise than previous).
  • Einstein delay: $\gamma = 4.59(2) \times 10^{-4}$ s.
  • Orbital period derivative: $\dot{P}_b = -5.65(2) \times 10^{-13}$ s s$^{-1}$ (15$\times$ more precise).
• Masses (assuming GR):
  • Total mass: $M_{total} = 2.6133(1) M_\odot$.
  • Pulsar mass: $M_p = 1.316(5) M_\odot$.
  • Companion mass: $M_c = 1.297(5) M_\odot$.
  • These values are consistent with a Double Neutron Star (DNS) system and agree with GR predictions to better than 1%.

The $\dot{x}$ Anomaly

• The authors measured a secular change in the projected semi-major axis: $\dot{x} = -1.8(6) \times 10^{-13}$ lt-s s$^{-1}$.
• Implications:
  • If the companion is a Neutron Star, it would need to be a recycled pulsar spinning faster than 10 ms to generate this $\dot{x}$ via spin-orbit coupling. However, no such pulsar signal has been detected despite deep searches.
  • If the companion is a White Dwarf, a spin period of $\sim$4 minutes could explain the observed $\dot{x}$ via the Lense-Thirring effect and classical quadrupole coupling. This aligns with the system's evolutionary history where the WD formed first and was spun up by accretion from the pulsar's progenitor.
• Mass Correlation: Including $\dot{x}$ in the fit shifts the mass estimates by $\sim 3.5\sigma$ due to the correlation between $\dot{x}$ and $\gamma$. In this scenario, $M_c \approx 1.19 M_\odot$ and $M_p \approx 1.42 M_\odot$. While shifted, these masses remain consistent with a massive WD scenario.

Glitch and Timing Noise

• Glitch: A large glitch occurred at MJD 56664 with $\Delta\nu/\nu = 5.899(2) \times 10^{-6}$. It exhibited an exponential recovery with a timescale $\tau_d \approx 100$ days and a recovery fraction $Q = 0.005$.
• Quasi-periodicity: The residuals and $\dot{\nu}$ variations show a potential periodicity of $\sim$2 years. While a planetary companion was hypothesized, the authors suggest this may instead be intrinsic magnetospheric variability, noting that the 2-year period is also seen in the power spectral density of the timing noise.

5. Significance

• Testing General Relativity: The system provides a stringent test of GR. The consistency of the PK parameters with GR predictions (within 1%) reinforces the theory's validity in strong-field regimes.
• Companion Nature: The detection of $\dot{x}$ offers a unique window into the companion's nature. If confirmed as a fast-rotating WD, PSR J1906+0746 would be a critical laboratory for testing dipolar gravitational wave emission predicted by alternative theories of gravity (e.g., scalar-tensor theories).
• Future Prospects: The pulsar's emission is predicted to disappear due to geodetic precession by 2028, reappearing only in the 2070s. This makes the current 18-year dataset a crucial "snapshot" of the system's dynamics. Future observations with the Square Kilometre Array (SKA) will be required to measure proper motion precisely enough to separate kinematic effects from intrinsic orbital decay, enabling definitive tests of dipolar gravitational radiation.

In conclusion, this paper significantly refines the timing solution of PSR J1906+0746, providing the most precise mass measurements to date and presenting compelling, though not yet definitive, evidence that the companion may be a massive, fast-rotating white dwarf rather than a neutron star.