🔭 astrophysics

Born to be recycled: a comprehensive population synthesis of the Galactic millisecond pulsars

This paper presents a comprehensive population synthesis of Galactic millisecond pulsars using the SEVN code and self-consistent pulsar evolution models to confirm the recycling hypothesis, reproduce observed radio and gamma-ray populations, predict fewer than 220 unidentified sources in the 4FGL catalogue, and constrain the geometric and physical conditions required for their detection.

Original authors: Mattéo Sautron, Jérôme Pétri, Dipanjan Mitra, Adélie Dupuy--Junet, Marie-Eloïse Pietrin

Published 2026-03-11

📖 6 min read🧠 Deep dive

CC BY 4.0

Original authors: Mattéo Sautron, Jérôme Pétri, Dipanjan Mitra, Adélie Dupuy--Junet, Marie-Eloïse Pietrin

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 Big Picture: Cosmic Spin-Doctors

Imagine the universe is a giant dance floor. Most dancers (stars) eventually get tired, slow down, and stop dancing. But there's a special group of dancers called Pulsars. They are the "dead" hearts of massive stars that collapsed, yet they spin incredibly fast, flashing beams of light like lighthouses.

Usually, these dancers slow down over time. But then, there's a special group called Millisecond Pulsars (MSPs). These are the oldest dancers, yet they are spinning faster than anyone else—so fast they complete a rotation in just a few thousandths of a second.

The Big Question: How do old, tired stars suddenly start spinning like a top again?

The Answer: They get a "second wind" from a partner. This paper is a massive computer simulation that tries to figure out exactly how this "recycling" happens, how many of these fast spinners exist, and where they are hiding in our galaxy.

The Recipe: How the Simulation Works

The authors built a virtual universe in their computers to watch stars grow up, pair up, and evolve. Think of it like a high-tech version of The Sims, but for dying stars.

The Setup: They started with millions of binary pairs (two stars dancing together). One star is a newborn neutron star (the future pulsar), and the other is a normal main-sequence star (the partner).
The Dance (Binary Evolution): As they dance, they interact. Sometimes the partner star gets too big and spills its gas onto the neutron star.
The Spin-Up: This is the magic moment. When the neutron star eats this gas, it's like a figure skater pulling in their arms to spin faster. The gas transfers energy, spinning the neutron star up to incredible speeds. This is the "recycling" process.
The Rules: The simulation followed strict physics rules:
- Gravity: How the stars pull on each other.
- Magnetic Fields: How the stars' invisible magnetic shields decay over time (like a battery running out).
- The "Death Valley": A safety line. If a pulsar spins too slowly or its magnetic field gets too weak, it stops flashing and becomes invisible. The simulation checks if the stars survive long enough to be seen.

The Results: What Did They Find?

After running the simulation, the authors compared their virtual universe to the real one we see through telescopes. Here are the key takeaways:

1. The "Recycling" Theory Holds Up

The simulation successfully recreated the population of pulsars we actually see. This confirms that the "recycling" theory is correct: most of these fast-spinning old stars got their speed boost by stealing gas from a partner.

2. The "Spiders" are Real (and Heavy)

Some of these recycled pulsars are "Spiders."

Black Widows: The pulsar is so hungry it eats its tiny partner completely.
Redbacks: The pulsar is eating a slightly larger partner.
The simulation showed that these "spider" pulsars are surprisingly heavy, often weighing about 1.8 times the mass of our Sun, and some even up to 2.7 times. This matches real-life observations of heavy pulsars, suggesting that neutron stars might be able to be even heavier than we thought.

3. The Mystery of the "Missing" Gamma Rays

We have a telescope called Fermi that looks for high-energy light (gamma rays). We see a lot of gamma rays coming from the center of our galaxy, but we can't find enough pulsars to explain them.

The Finding: The simulation suggests that the recycled pulsars born in the open galaxy (the "spiral arms") only contribute about 5% of this mysterious gamma-ray glow.
The Conclusion: If pulsars are the only source of this glow, then the missing 95% must be coming from a different place—likely from ancient star clusters that have drifted into the center of the galaxy, or pulsars born right in the galactic center. They are there, but they are too faint or hidden for us to see yet.

4. The "Alignment" Problem

Imagine a spinning top. Usually, the axis it spins on and the direction it points might be tilted.

The Finding: The simulation shows that for most recycled pulsars, the spin axis and the orbital path align perfectly (like a top spinning straight up). About 80% of them are perfectly aligned. This makes them easier to spot if we are looking from the right angle, but harder if they are tilted away from us.

5. A Secret Backdoor: The "White Dwarf Collapse"

The paper also checked a different way to make a pulsar: crushing a White Dwarf (a dead star) until it collapses.

The Finding: This happens, but it's rare. Only a tiny fraction (less than 0.01%) of these white dwarfs turn into the "mildly recycled" pulsars we see. It's a possible backdoor, but not the main highway.

The Future: What's Next?

The paper ends with a prediction for the future.

The SKA Telescope: A new, super-powerful radio telescope (the Square Kilometre Array) is coming online. The authors predict it will find 3 times more of these recycled pulsars than we have now.
Better Gamma Ray Eyes: If we build a gamma-ray telescope 10 times more sensitive than Fermi, we might find 5 times more of these hidden fast-spinners.

The Bottom Line

This paper is a massive "stress test" of our understanding of the universe. By simulating millions of years of stellar evolution, the authors confirmed that recycling via binary partners is the main way old stars get a second wind. They also mapped out where these stars are, how heavy they are, and why we still can't find all the sources of the gamma-ray glow in the center of our galaxy.

It's like solving a cosmic jigsaw puzzle: we now have a much clearer picture of the pieces, even if a few are still hidden in the box.

1. Problem Statement

Millisecond Pulsars (MSPs) are the oldest yet fastest-spinning neutron stars ( $P < 30$ ms). The prevailing hypothesis for their formation is the "recycling" scenario, where an old neutron star in a binary system accretes matter from a companion, gaining angular momentum and spinning up. However, several aspects of this process remain poorly constrained:

Formation Channels: It is unclear if accretion-driven recycling is the sole mechanism or if alternative channels (e.g., Accretion-Induced Collapse of White Dwarfs) contribute significantly.
Population Properties: Observational biases make it difficult to derive intrinsic statistical properties of the recycled population, such as birth rates, magnetic field decay timescales, mass distributions, and spin-orbit alignment.
The GeV Excess: The origin of the gamma-ray excess in the Galactic Center is debated. It is unknown whether it stems from dark matter or a population of MSPs, and if the latter, whether they originate from the Galactic bulge, migrated globular clusters, or the Galactic field.
Spider Systems: The formation mechanisms for "spider" pulsars (Black Widows and Redbacks) and their companion evolution are not fully understood.

2. Methodology

The authors developed a comprehensive population synthesis algorithm that self-consistently models the full binary evolution alongside the pulsar's internal evolution and detection probabilities.

A. Binary Evolution (SEVN Code)

Stellar Evolution: Used the SEVN (Stellar EVolution for N-body) code to simulate binary systems starting with a Main Sequence (MS) companion and a newborn neutron star (NS).
Mass Transfer: Modeled Roche Lobe Overflow (RLO) using Eggleton's formula. Mass transfer stability was determined by critical mass ratios. Stable mass transfer was modeled with non-conservative accretion (efficiency $f_{MT} = 0.5$ ). Unstable transfer leading to Common Envelope (CE) phases was treated as a disruption channel (no spin-up during CE in this model).
Initial Conditions:
- Birth: NSs born in Galactic spiral arms with a truncated normal radial distribution ( $\mu_R = 11.3$ kpc).
- Kick Velocity: Adopted a low dispersion ( $\sigma_v = 70$ km/s) to ensure binary survival after the first supernova.
- Masses: Progenitor masses drawn from a Kroupa IMF; NS birth masses from a Gaussian distribution ( $\mu = 1.33 M_\odot$ ). Companion mass ratios drawn from a power law.
- Geometry: Isotropic distributions for spin-orbit angles ( $\alpha$ ) and viewing angles ( $\zeta$ ).

B. Pulsar Evolution (Custom Code)

Spin Evolution: Solved a system of differential equations for angular velocity ( $\Omega$ $Ω$ ), spin-orbit angle ( $\alpha$ $α$ ), and magnetic inclination ( $\chi$ $χ$ ).
- Torques: Included accretion torque ( $n_1$ ), magnetic braking ( $n_2$ ), radiation loss ( $n_3$ ), and gravitational wave (GW) braking for very fast rotators ( $\Omega \ge 960$ rad/s).
- Magnetic Field Decay: Implemented a power-law decay $B(t) = B_0(1 + t/\tau_d)^{-1/\alpha_d} + B_{min}$ . Three discrete decay timescales ( $\tau_d$ ) were used: 0.05, 1.4, and 2.3 Myr, with probabilities 0.77, 0.115, and 0.115 respectively.
- Accretion-Induced Decay: Modeled magnetic field burial due to accreted mass.
Spatial Evolution: Integrated binary trajectories in a 4-component Galactic potential (Bulge, Disk, Halo, Nucleus) using the PEFRL algorithm.
Death Valley: Defined a "death valley" in the $P-\dot{P}$ diagram based on Mitra et al. (2020) to determine when coherent radio emission ceases.

C. Detection Models

Radio: Modeled beaming fractions using polar cap geometry. Flux density calculated based on spin-down luminosity ( $\dot{E}^{0.125}$ ) and distance. Detection thresholds applied for FAST GPPS and Parkes Multibeam (PMPS) surveys, accounting for dispersion, scattering, and instrument sensitivity.
Gamma-ray: Modeled emission using the Striped Wind scenario. Luminosity derived from a fundamental plane relation ( $L_\gamma \propto B^{0.11} \dot{E}^{0.51}$ ). Detection based on Fermi/LAT sensitivity maps and geometric visibility ( $|\zeta - \pi/2| \le \chi$ ).

D. Alternative Channel (AIC)

A separate simulation of 1 million binary systems with Oxygen-Neon White Dwarfs (ONeWDs) to test the Accretion-Induced Collapse (AIC) channel.

3. Key Contributions

First Full Evolutionary Synthesis: Unlike previous studies that often started evolution after the spin-up phase or assumed constant magnetic fields, this work tracks the system from the birth of the NS through the accretion phase to the present day.
Self-Consistent Geometry: Simultaneously models the evolution of the spin-orbit angle ( $\alpha$ ) and magnetic inclination ( $\chi$ ) under the influence of accretion torques and radiation losses.
Multi-Wavelength Detection: Applies rigorous detection criteria for both radio and gamma-ray wavelengths to compare directly with specific surveys (FAST, Parkes, Fermi-LAT).
AIC Contribution: Quantifies the specific contribution of ONeWD collapse to the mildly recycled pulsar population.

4. Key Results

A. Population Statistics & P-Ċ Diagram

The simulation successfully reproduces the observed $P-\dot{P}$ distribution for radio and gamma-ray MSPs (validated by Kolmogorov-Smirnov tests).
Birth Rate: Estimated at $3.3 - 4.1 \times 10^{-6} \text{ yr}^{-1}$ .
Accretion Duration: Typical accretion phases last $10^8 - 10^9$ years.
Magnetic Field Decay: The best fit requires a dominant population with a short decay timescale ( $\tau_d = 0.05$ Myr).

B. Mass and Geometry

Masses: Recycled pulsars have a mean mass of $\sim 1.8 M_\odot$ , with a tail extending up to $2.7 M_\odot$ . This supports recent observations of massive "spider" pulsars.
Alignment: The spin-orbit alignment timescale is short ( $\sim$ kyr). Consequently, $\sim 80\%$ of the binary population has a spin-orbit angle $\alpha \approx 0^\circ$ .
Viewing Angles: High viewing angles ( $\zeta > 45^\circ$ ) are required for detection. Aligned rotation/magnetic axes ( $\chi \le 10^\circ$ ) are rarely observed, likely because accretion tends to misalign the axes or such systems become too faint.

C. Spider Systems

The simulation produces distinct Black Widow (BW) and Redback (RB) populations.
Discrepancy: The model underestimates the total number of observed spider systems (by a factor of $\sim 3$ ).
Companions: Simulated BWs often have HeWD or COWD companions, whereas observed BWs have ultra-light companions. This suggests alternative formation channels (e.g., high-metallicity disks, dynamical capture) or the need for higher accretion/evaporation efficiencies.

D. Alternative Formation (AIC)

Only $\sim 0.9\%$ of simulated ONeWD binaries collapse into NSs.
Of these, $\sim 12.2\%$ fall into the "mildly recycled" region ( $30 \le P \le 1000$ ms).
Contribution: AIC contributes $\sim 7.6 \times 10^{-3}\%$ of ONeWDs to the mildly recycled population. This channel may explain some long-period recycled pulsars not accounted for by standard core-collapse supernovae.

E. The GeV Excess

MSPs born in the Galactic field contribute only $\sim 5\%$ to the GeV excess (approx. 190 MSPs in the central 0.7 kpc).
Implication: If the GeV excess is entirely pulsar-origin, the majority must come from globular clusters that migrated inward or stars born in the bulge, which are currently undetectable.

F. Future Predictions

Unidentified Sources: Fewer than 220 unidentified pulsars are expected in the 4FGL catalog.
SKA & Future Instruments: The Square Kilometer Array (SKA) is predicted to increase radio detections by $\sim 3\times$ . A gamma-ray instrument 10x more sensitive than Fermi/LAT would increase detections by $\sim 5\times$ and reveal $\sim 10$ MSPs in the Galactic Center.

5. Significance

This study provides the most comprehensive framework to date for understanding the life cycle of recycled pulsars. By linking binary evolution, magnetic field physics, and observational selection effects, it:

Validates the Recycling Scenario: Confirms that accretion-driven formation is the dominant channel for MSPs, provided magnetic field decay is rapid.
Constrains Neutron Star Physics: Supports the existence of high-mass neutron stars ( $>2 M_\odot$ ) and suggests specific magnetic field decay timescales.
Clarifies the GeV Excess: Strongly suggests that the Galactic Center GeV excess cannot be explained solely by field-born MSPs, pointing towards a hidden population from globular clusters or the bulge.
Guides Future Surveys: Provides concrete predictions for upcoming radio (SKA) and gamma-ray missions, helping to prioritize search strategies for the remaining unidentified sources.