🔭 astrophysics

Neutron Star Mass across Binary Pulsar Subpopulations: Mass-Spin Correlation, Mass Distributions, and Moment of Inertia Effects

This study employs hierarchical Bayesian analysis of approximately 50 Galactic binary radio pulsars to reveal a moderate anti-correlation between neutron star mass and spin period in recycled populations, while confirming distinct formation pathways for pulsar-white dwarf and double neutron star systems and identifying a marginal mass difference based on companion white dwarf composition.

Original authors: Debatri Chattopadhyay

Published 2026-06-30

📖 5 min read🧠 Deep dive

CC BY 4.0

Original authors: Debatri Chattopadhyay

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

Imagine the universe as a giant cosmic gym. In this gym, neutron stars are the heavyweights—ultra-dense, city-sized balls of matter left over from exploded stars. Some of these heavyweights are lonely, but many are in binary pairs, dancing around a partner.

This paper is like a detective story where the author, Debatri Chattopadhyay, tries to figure out how these neutron stars got their specific "weights" (mass) and how fast they are spinning. The author looks at about 50 of these cosmic dancers and splits them into two main teams:

The "Recycled" Team (PSR–WD): A neutron star paired with a white dwarf (a dead, low-mass star). These neutron stars have been "recycled"—they spun up to incredible speeds by stealing mass and energy from their partner.
The "Double Neutron Star" Team (DNS): Two neutron stars dancing together.

Here is what the paper found, explained simply:

1. The "Heavy Spin" Mystery

The main question was: Does a heavier neutron star spin faster or slower after being recycled?

Think of it like a figure skater.

Theory A (The "Accretion" Skater): If a skater gains weight while spinning, they might spin faster because they are grabbing more energy. This would mean heavier stars spin faster (a negative correlation: high mass = low spin time).
Theory B (The "Moment of Inertia" Skater): If a skater is lighter, they are easier to spin up. A lighter skater might spin faster than a heavy one with the same push. This would mean lighter stars spin faster (a positive correlation).

The Verdict: The data shows a moderate hint that the heavier recycled stars are spinning faster (the "Accretion" theory). The author calls this "moderate evidence," meaning it's a strong clue but not a slam-dunk proof yet. The math suggests that if you have a heavier neutron star, it likely spun up more because it stole more mass from its partner.

2. The "Double Star" Twist

When the author looked only at the Double Neutron Star team (the two heavyweights dancing together), the trend flipped slightly. These systems leaned toward the idea that lighter stars spin faster. However, there were only 10 of these systems, so the author says this isn't statistically strong enough to be a rule. It's like looking at a small group of people and guessing the height of the whole population—it might be a fluke.

3. The "Partner" Effect

The paper also looked at the partners (the white dwarfs).

Helium vs. Carbon: Neutron stars paired with Helium white dwarfs seem to be slightly heavier than those paired with Carbon/Oxygen white dwarfs.
The Analogy: Imagine two people eating at a buffet. The person with the Helium partner seems to have eaten a bit more "buffet food" (mass) than the one with the Carbon partner. This makes sense because Helium partners likely had a longer, more stable feeding time, allowing the neutron star to grow bigger. However, the difference is small, so the author calls it "marginally significant"—it's a whisper, not a shout.

4. The "Bimodal" Confusion (One Hump or Two?)

Scientists have argued for years whether the weights of these recycled stars form two distinct groups (like a bimodal distribution with two humps) or just one big, messy group with a long tail.

The Paper's Take: The data looks like two humps, but the author argues it's actually one single group that is skewed. Imagine a bell curve that has been stretched out to the right. It's not two separate groups of stars; it's one group where some stars got a little extra "accretion" (mass) and ended up much heavier, creating a long tail.

5. Why Can't We Be 100% Sure?

The author explains a tricky physics problem. The reason we can't perfectly distinguish between "Heavier stars spin faster because they ate more" and "Heavier stars spin faster because of how their internal structure works" is that neutron stars are weirdly predictable.

The Analogy: Imagine trying to tell if a car is fast because it has a big engine or because it's made of light materials. If the car's weight and its engine size always go up together in a straight line, you can't tell which one is causing the speed. Similarly, for neutron stars, their "moment of inertia" (how hard they are to spin) changes almost perfectly in line with their mass. This makes the two theories look identical in our current data.

6. What's Next?

The paper concludes that we need more data to be sure.

The Forecast: To turn this "moderate hint" into a "definite discovery," we need to roughly double or triple the number of neutron stars with precise mass measurements.
The Future: New telescopes (like the SKA) will find more of these stars, but measuring their exact mass takes years of careful observation. The author predicts we might get a solid answer in the next decade or two.

Summary

The paper suggests that heavier recycled neutron stars tend to spin faster, likely because they stole more mass from their partners. However, the evidence is currently "moderate," not "conclusive," because the sample size is small and the physics of neutron stars makes two different theories look very similar. The author also confirms that neutron stars in double-star systems have a strong link between their partner's mass and the orbit's shape, a clue about how they were born.

In short: We have a strong hunch about how these cosmic heavyweights grew up, but we need to wait for more data to be absolutely certain.

Technical Summary: Neutron Star Mass across Binary Pulsar Subpopulations

Problem Statement
The formation and evolution of neutron stars (NS) in binary systems are imprinted on their masses, spin periods, and orbital parameters. While the "recycling" paradigm posits that accretion from a companion spins up NSs to millisecond periods, the physical origin of the observed mass distribution remains debated. Specifically, it is unclear whether high-mass recycled pulsars ( $\sim 1.8 M_\odot$ ) were born massive or grew via accretion. Furthermore, two competing physical mechanisms predict opposite correlations between NS mass ( $M$ ) and spin period ( $P$ ) in recycled populations:

Accretion-driven recycling: Angular momentum gain is coupled with mass gain, predicting a negative correlation ( $\rho < 0$ ; faster spins associated with higher masses).
Moment of inertia (MoI) mechanism: Lower-mass NSs have smaller moments of inertia and thus spin up more efficiently for a given accreted angular momentum, predicting a positive correlation ( $\rho > 0$ ).

Current samples are limited by measurement uncertainties and selection effects, making it difficult to distinguish these mechanisms or determine if the recycled mass distribution is intrinsically bimodal or a skewed unimodal population.

Methodology
The authors present a hierarchical Bayesian analysis of $\sim 50$ Galactic binary radio pulsars with measured NS masses, classified into Pulsar–White Dwarf (PSR–WD) and Double Neutron Star (DNS) systems.

Data Curation: The sample is compiled from established pulsar mass catalogues. Repeated measurements are combined via inverse-variance weighting, with caution applied to systems with discrepant measurements (e.g., J0348+0432, B1855+09).
Classification: Systems are classified by binary type (PSR–WD primaries, DNS primaries, DNS secondaries) rather than spin period to avoid circularity. White dwarf companions are sub-classified as Helium (He) or Carbon-Oxygen/Oxygen-Neon (CO/ONe) based on a mass threshold of $0.45 M_\odot$ .
Statistical Framework: A hierarchical bivariate Gaussian model is employed to estimate intrinsic correlations between mass and recycling proxies (spin period, orbital period, eccentricity, inclination). This approach explicitly models measurement errors as latent variables to correct for regression dilution (attenuation bias).
Selection Corrections: The analysis incorporates corrections for orbital-Doppler detectability (using the $\gamma_1$ efficiency factor) and accounts for inclination biases inherent in Shapiro-delay mass measurements.
Model Comparison: To test the MoI mechanism, the authors compare regression models of spin period against mass directly versus against the moment of inertia (derived from universal relations), using the Bayesian Information Criterion (BIC).

Key Contributions and Results

Mass–Spin Period Anti-Correlation:
The study finds moderate evidence for an anti-correlation between NS mass and spin period in the pooled recycled population ( $N=47$ primaries).
- Result: $\rho = -0.26$ with a 90% credible interval of $[-0.47, -0.01]$ and a one-sided posterior probability $P(\rho < 0) = 0.96$ (approx. $1.7\sigma$ ).
- Robustness: The result holds when including candidate DNSs and applying selection corrections.
- Subpopulation Nuance: The PSR–WD sample alone shows a suggestive trend ( $\rho = -0.19$ ), while the DNS primaries alone lean slightly positive ( $\rho = +0.13$ ), though the DNS sample is too small ( $N=10$ ) to be conclusive. The pooled significance is driven largely by the inclusion of DNS primaries, which occupy a region of lower mass and longer period, extending the negative diagonal in the mass–spin plane.
Degeneracy of Mechanisms:
While the negative correlation is consistent with accretion-driven recycling, the authors cannot statistically distinguish this from the MoI mechanism.
- Reason: Over the observed mass range, the NS moment of inertia $I(M)$ is nearly linear. Consequently, a regression of spin period on mass is statistically indistinguishable from a regression on $I(M)$ (BIC difference $\Delta \text{BIC} = -0.2$ ).
Mass Distributions and Bimodality:
- Recycled vs. Non-Recycled: Recycled primaries are significantly more massive ( $1.50 M_\odot$ ) than non-recycled DNS secondaries ( $1.30 M_\odot$ ), consistent with previous findings.
- Bimodality: Within the recycled population, a two-Gaussian mixture is preferred over a single symmetric Gaussian, but a single right-skewed (skew-normal) distribution provides the most parsimonious description (lowest BIC). This suggests the apparent bimodality may be an artifact of a single population with a high-mass tail rather than two distinct birth populations.
Companion Type and Mass:
Neutron stars with Helium white dwarf companions are marginally more massive ( $1.547 M_\odot$ ) than those with CO/ONe companions ( $1.484 M_\odot$ ), a difference of $\sim 0.06 M_\odot$ . This is of marginal statistical significance ( $p=0.11$ ) but aligns with the expectation that He-WD channels involve more extensive mass transfer.
DNS Companion Mass–Eccentricity Correlation:
In a hierarchical framework, the authors confirm a strong positive correlation between the companion (second-born) mass and orbital eccentricity in DNS systems ( $\rho = +0.82$ ). This supports the theoretical link between low-kick supernovae (producing low-mass, low-eccentricity systems) and the formation of DNSs.

Significance and Claims
The paper frames its findings as moderate or suggestive evidence rather than a definitive detection, primarily due to the modest sample size ( $N \approx 50$ ).

Interpretation: The results support a "two-channel" picture: PSR–WD systems are "accretion-grown," exhibiting the mass–spin anti-correlation and mass offsets, while DNS systems are "birth-mass-dominated," where masses reflect formation conditions rather than recycling history.
Limitations: The inability to distinguish between accretion and MoI mechanisms is a direct consequence of the near-linear relationship between mass and moment of inertia in the observed range.
Future Prospects: The authors forecast that achieving a confident $3\sigma$ detection of the mass–spin correlation will require a sample size of $N \gtrsim 125–150$ precise-mass recycled systems, likely an SKA-era achievement. They note that breaking the accretion/MoI degeneracy may require measuring the spin and mass of non-recycled NSs (e.g., second-born pulsars in DNS systems) to isolate birth properties from recycling effects.

The work reframes the debate on recycled mass distributions, suggesting that the observed high-mass tail is better described by a skewed unimodal distribution driven by accretion history than by a distinct bimodal birth population.