🔭 astrophysics

Shapiro Delay Measurements from Fifteen Years of PSR J1231$-$1411 Radio Observations

This paper presents a 15-year radio timing analysis of PSR J1231$-$1411 using the Green Bank and Nançay telescopes to measure its orbital inclination and constrain the masses of the pulsar and its companion via Shapiro delay, yielding results that, while providing weak mass constraints, inform future neutron star equation of state studies and recent NICER X-ray observations.

Original authors: H. Thankful Cromartie, Matthew Kerr, Scott M. Ransom, Paul S. Ray, Lucas Guillemot, Ismaël Cognard, Emmanuel Fonseca, Gilles Theureau

Published 2026-03-02

📖 6 min read🧠 Deep dive

CC BY 4.0

Original authors: H. Thankful Cromartie, Matthew Kerr, Scott M. Ransom, Paul S. Ray, Lucas Guillemot, Ismaël Cognard, Emmanuel Fonseca, Gilles Theureau

Original paper licensed under CC BY 4.0 (http://creativecommons.org/licenses/by/4.0/). ✨ This is an AI-generated explanation of the paper below. It is not written or endorsed by the authors. For technical accuracy, refer to the original paper. Read full disclaimer

The Cosmic Detective Story: Weighing a Ghost

Imagine you are trying to weigh a ghost. You can't put it on a scale, and you can't touch it. All you have is a lighthouse that flashes incredibly fast (a pulsar) and a heavy, invisible rock orbiting it (a white dwarf star).

This paper is about a team of astronomers who spent 15 years trying to figure out how heavy that ghost (the pulsar) is, using a clever trick from Einstein's theory of gravity called the Shapiro Delay.

1. The Setup: A Cosmic Dance

The star in question, PSR J1231−1411, is a "millisecond pulsar." Think of it as a cosmic lighthouse spinning 271 times every second. It's in a tight dance with a smaller, dead star (a white dwarf).

Every time the pulsar spins, it sends out a radio beam. We catch these beams on Earth with giant radio dishes (the Green Bank Telescope and the Nançay Radio Telescope). Because the pulsar spins so regularly, it acts like a perfect clock. If something messes with that clock, we know something interesting is happening.

2. The Trick: The "Gravity Traffic Jam"

Here is the magic trick: Shapiro Delay.

Imagine you are driving on a highway. Usually, you drive in a straight line at a constant speed. But imagine there is a massive, invisible hill in the middle of the road. Even if you keep your foot on the gas, the car has to slow down just a tiny bit as it climbs the hill and speeds up as it goes down.

In space, the "hill" is the gravity of the white dwarf companion.

When the pulsar is on the far side of its orbit (behind the white dwarf from our perspective), its radio signals have to pass right through the white dwarf's gravity well to get to us.
According to Einstein, gravity bends space and time. This makes the signal take a tiny bit longer to arrive than it would if the white dwarf weren't there.
This delay is called the Shapiro Delay.

The Catch: This delay is only visible if the orbit is tilted just right. If the orbit is flat (like a coin spinning on a table), the signal never passes "behind" the companion from our view. If the orbit is edge-on (like a coin spinning on its edge), the signal passes right through the gravity well.

3. The Problem: A "Troublesome" System

The astronomers wanted to measure this delay to weigh the pulsar. But PSR J1231−1411 is a bit of a troublemaker:

The Angle is "Meh": The orbit is tilted at about 80 degrees. That's close to edge-on, but not quite perfect. It's like trying to hear a whisper from someone who is slightly turned away from you. The "gravity traffic jam" is very subtle.
The Signal is Fuzzy: The pulsar is also affected by "scintillation." Imagine looking at a star through a heat haze; the light twinkles and distorts. This pulsar twinkles so much in the radio waves that sometimes it disappears completely, making it hard to get a clean reading.
The Noise: There is a lot of static (noise) in the data, like trying to hear a whisper in a crowded stadium.

4. The Solution: 15 Years of Patience and Math

To solve this, the team didn't just look at the data once. They used three different mathematical approaches (think of them as three different detectives solving the same case):

The Grid Search: They tested millions of possible combinations of weights and angles to see which one fit the data best.
The "Smart Guess" (Bayesian): They used what we know about how stars die (white dwarf evolution) to make an educated guess about how heavy the companion should be. This helped narrow down the possibilities.
The "Super-Computer" Method: They used a complex method that figured out the noise and the star's weight at the same time, rather than guessing the noise first.

5. The Results: A Heavy Ghost

After all that work, here is what they found:

The Companion: The dead star (white dwarf) is about 0.23 times the mass of our Sun.
The Pulsar: The pulsar is heavy! It weighs about 1.87 times the mass of our Sun.
- Note: Because the data was a bit fuzzy, there is a large "margin of error." The real weight could be as low as 1.2 or as high as 3.0 solar masses. However, even with the uncertainty, it confirms this is a very heavy neutron star.
The Angle: They measured the tilt of the orbit very precisely: 79.8 degrees. This is the most reliable number they got.

6. Why Does This Matter?

You might ask, "Why do we care about weighing a dead star?"

Neutron stars are the densest objects in the universe. A teaspoon of their material would weigh a billion tons. Scientists are trying to figure out what happens to matter under such extreme pressure. This is called the Equation of State.

If a neutron star is too heavy, it should collapse into a black hole.
If we find a neutron star that is very heavy (like 2 Suns), it tells us that the "stuff" inside it is incredibly strong and stiff.

This specific pulsar was also being studied by the NICER telescope (an X-ray observatory on the International Space Station). The radio measurements from this paper helped the X-ray team figure out the size and shape of the hot spots on the pulsar's surface.

The Takeaway

This paper is a story of persistence. The astronomers spent 15 years collecting data, dealing with fuzzy signals, and running complex computer simulations to measure a tiny delay in a radio signal.

Even though the result wasn't a "perfect" measurement (the error bars are still a bit wide), it provided crucial information that helped other scientists understand the physics of the universe's densest matter. It's a reminder that in science, even a "troublesome" system that refuses to give up its secrets easily can teach us a lot if we just keep listening long enough.

1. Problem Statement

The primary scientific goal is to constrain the Equation of State (EoS) of supranuclear-density matter in neutron stars. This is achieved by precisely measuring the masses of neutron stars, particularly massive ones.

The Target: PSR J1231−1411 is a millisecond pulsar (MSP) in a binary system with a white dwarf companion. It is a key target for multi-messenger astronomy, having been observed by radio telescopes (Green Bank and Nançay) and the NICER X-ray telescope.
The Challenge: The system has a high orbital inclination but not a perfectly edge-on orientation ( $i \approx 80^\circ$ ), resulting in a weak Shapiro delay signal (a general relativistic time delay caused by the companion's gravity). Furthermore, the pulsar exhibits strong diffractive scintillation, heterogeneous data from different telescope backends, and significant timing noise (both white and red). These factors make the mass measurement difficult and prone to covariance with noise parameters, leading to large uncertainties in previous attempts.

2. Methodology

The authors analyzed 15 years of radio timing data (2009–2026) from the Green Bank Telescope (GBT) and the Nançay Radio Telescope (NRT), supplemented by Fermi-LAT $\gamma$ -ray data context.

Observations:
- GBT: Included a new 22-hour campaign in Spring 2019 specifically targeting superior conjunction (where Shapiro delay is maximal) to improve constraints.
- NRT: Near-continuous observations using BON and NUPPI backends. Significant attention was paid to clock corrections (GPS vs. Rubidium) and polarization calibration changes in 2019.
Data Processing:
- Time-of-Arrival (TOA) extraction using PSRCHIVE (Matrix Template Matching).
- Noise modeling using ENTERPRISE to determine white noise parameters (EFAC, EQUAD) and search for red noise (spin noise).
- Timing analysis using PINT (Pulsar Timing Software).
Analysis Techniques: The authors employed three distinct statistical approaches to derive masses ( $m_p$ $m_{p}$ ) and companion masses ( $m_c$ $m_{c}$ ):
1. $\chi^2$ Gridding: A frequentist approach fixing noise parameters and scanning a grid of $m_c$ and $\sin(i)$ values.
2. Bayesian MCMC (Fixed Noise): Using PINT's emcee sampler with noise parameters fixed a priori. This included astrophysical priors based on the Tauris & Savonije (1999) relationship between white dwarf mass and orbital period (TS99).
3. Fully-Generalized Bayesian Fit: A novel approach using a branch of ENTERPRISE that simultaneously samples the noise and timing models, treating some parameters as non-linear to avoid hiding covariances.

3. Key Contributions

Comprehensive Noise Characterization: The paper highlights the sensitivity of mass measurements to noise modeling. It demonstrates that improper handling of clock corrections, polarization calibration, and the inclusion of the orbital period derivative ( $\dot{P}_b$ ) can drastically alter the derived white noise parameters (EQUAD), which in turn corrupts the mass measurement.
Methodological Comparison: It provides a rare comparison between fixed-noise Bayesian fits and fully-generalized simultaneous sampling, showing that for systems with weak Shapiro delays, the choice of noise model significantly impacts the posterior distributions.
New Observational Constraints: The inclusion of a targeted GBT campaign and the re-analysis of 15 years of archival data with updated calibration procedures.
Multi-messenger Context: The results serve as conservative priors for NICER X-ray light curve modeling (Salmi et al. 2024), helping to break degeneracies in the mass-to-radius ( $M/R$ ) ratio for this specific source.

4. Key Results

The study yields the following constraints (68.3% Confidence Intervals):

Orbital Inclination ( $i$ ): Measured with relatively high precision due to the geometry, despite the weak signal.
- $i = 79.80^{+3.47}_{-4.70}^\circ$ (Fixed noise MCMC with TS99 prior).
- $i = 81.06^{+3.36}_{-3.51}^\circ$ (Fully-generalized fit).
Companion Mass ( $m_c$ ):
- $0.23^{+0.09}_{-0.06} M_\odot$ (Fixed noise).
- $0.21^{+0.06}_{-0.05} M_\odot$ (Fully-generalized).
- These values are consistent with TS99 evolutionary predictions for a 1.86-day orbit.
Pulsar Mass ( $m_p$ ):
- Without mass cap: $1.87^{+1.11}_{-0.67} M_\odot$ (Fixed noise).
- With $m_p < 3 M_\odot$ prior: $1.71^{+0.70}_{-0.56} M_\odot$ .
- Fully-generalized fit: $1.62^{+0.73}_{-0.58} M_\odot$ .
Astrometry: The paper reports the first significant measurement of $\dot{P}_b$ (orbital period derivative) and parallax for this system. The derived distance ( $\sim 600$ pc) is consistent with the Dispersion Measure (DM) distance.

Critical Finding on Noise: The authors found that adding newer, higher-quality data initially worsened the mass constraint because the new data required the inclusion of $\dot{P}_b$ . Without $\dot{P}_b$ in the model, the noise parameter EQUAD absorbed the signal, leading to a high, well-constrained EQUAD value that inflated the mass uncertainty. This highlights the danger of "fixed-noise" approaches in systems with borderline-significant noise processes.

5. Significance

Equation of State (EoS) Constraints: While the mass uncertainty for PSR J1231−1411 remains large compared to "gold standard" edge-on systems (like PSR J1614−2230), the results provide a necessary lower bound and consistency check for the neutron star EoS. The companion mass confirms it is a typical low-mass white dwarf.
Methodological Cautionary Tale: The paper serves as a critical guide for analyzing "troublesome" pulsars. It demonstrates that for systems with low inclinations and weak Shapiro delays, the interplay between timing model parameters (like $\dot{P}_b$ ) and noise parameters is non-trivial. Blindly fixing noise parameters can lead to biased mass estimates.
Future Prospects: The authors argue that next-generation radio telescopes (DSA-2000, ngVLA, SKA) with higher sensitivity and flexible scheduling will be essential to overcome scintillation and achieve the precision needed to constrain the EoS for such systems.
NICER Synergy: The radio-derived inclination constraint was instrumental in the recent NICER analysis of PSR J1231−1411, proving the value of combining radio timing with X-ray pulse profile modeling.