🔭 astrophysics

Forming Double Neutron Stars using Detailed Binary Evolution Models with POSYDON: Comparison to the Galactic Systems

Using the POSYDON binary population synthesis code, this study reveals that Galactic double neutron stars form via two distinct common envelope subchannels (Case B and Case C) requiring specific core definitions or high ejection efficiencies and low supernova kicks, thereby explaining observed orbital period splits and constraining the detailed evolutionary physics of these systems.

Original authors: Abhishek Chattaraj, Jeff J. Andrews, Simone S. Bavera, Max Briel, Debatri Chattopadhyay, Tassos Fragos, Seth Gossage, Vicky Kalogera, Konstantinos Kovlakas, Matthias U. Kruckow, Camille Liotine, Kyle

Published 2026-06-12

📖 6 min read🧠 Deep dive

CC BY 4.0

Original authors: Abhishek Chattaraj, Jeff J. Andrews, Simone S. Bavera, Max Briel, Debatri Chattopadhyay, Tassos Fragos, Seth Gossage, Vicky Kalogera, Konstantinos Kovlakas, Matthias U. Kruckow, Camille Liotine, Kyle A. Rocha, Philipp M. Srivastava, Meng Sun, Elizabeth Teng, Zepei Xing, Emmanouil Zapartas

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: A Cosmic Detective Story

Imagine the universe as a giant, chaotic dance floor. Most of the time, stars dance alone. But sometimes, two massive stars are born together, holding hands. Over millions of years, they go through a wild, dramatic dance that can end in one of two ways: they either drift apart forever, or they crash into each other and merge.

This paper is a detective story about Double Neutron Stars (DNS). These are the "heavyweights" of the cosmic dance floor—two incredibly dense, dead stars (neutron stars) orbiting each other. We know about 25 of these pairs in our own galaxy, the Milky Way. The authors used a super-computer simulation called POSYDON to figure out exactly how these pairs are formed and why they look the way they do.

The Main Discovery: Two Different Paths to the Same Destination

The biggest surprise in this paper is that there isn't just one way to make a Double Neutron Star. The authors found that these pairs form through two distinct "sub-channels," like two different highways leading to the same city.

The difference depends on the "age" and "maturity" of the second star when the drama begins:

The "Young" Path (The Helium Core Channel):
- The Scenario: Imagine a young star that hasn't had much time to grow up. It still has a core made of helium (like a teenager who hasn't finished high school yet).
- The Drama: When this young star tries to expand, it gets too close to its partner (the first neutron star). They get tangled in a giant, messy cloud of gas called a Common Envelope.
- The Outcome: To survive this tangle, the system needs a lot of energy to kick the gas cloud away. If they succeed, they end up in a very tight, fast orbit. These are the pairs that will eventually crash into each other and create gravitational waves (the "merging" ones).
The "Old" Path (The Carbon-Oxygen Core Channel):
- The Scenario: Now imagine an older, more mature star. It has had time to burn through its fuel and has a core made of carbon and oxygen (like a retiree who has finished their career).
- The Drama: This older star also gets tangled in a gas cloud with its partner.
- The Outcome: Because the older star's gas cloud is looser and easier to push away, the system survives the tangle much more easily. However, they end up in a much wider, slower orbit. These pairs are too far apart to ever crash into each other within the lifetime of the universe.

The Analogy: Think of it like two couples trying to break up a fight in a crowded room.

The Young Couple (Helium core) is fighting in a tight, cramped space. It takes a lot of effort (energy) to push the crowd away, but if they do, they end up hugging very tightly.
The Old Couple (Carbon-Oxygen core) is fighting in a spacious room. It's easier to push the crowd away, but they end up standing far apart from each other.

The "Kick" Problem: Why the Dance Floor is Slippery

When a star dies, it explodes as a supernova. This explosion often gives the new neutron star a "kick," like a cue ball being struck on a pool table.

The Old Theory: Scientists used to think these kicks were huge and random (like a cue ball being hit with a sledgehammer).
The New Finding: The authors found that if the kicks were that big, most of the pairs would fly apart and never become Double Neutron Stars.
The Solution: To match what we actually see in the sky, the "kicks" must be much gentler (like a soft tap). The paper suggests that the second star in the pair gets a very gentle kick because it has been stripped of most of its weight before it explodes. This allows the pair to stay together.

The "Recycling" of the First Star

Before the second star explodes, the first neutron star (the older one) gets a "makeover." As the second star sheds its gas, the first star eats some of it. This is like a cosmic diet plan that actually makes the star spin faster and brighter, turning it into a "recycled pulsar."

The paper shows that the "Young Path" (Helium core) provides a better diet for the first star, spinning it up faster. This explains why the merging pairs we see in the sky (Sub-population i) spin very fast, while the non-merging pairs (Sub-population ii) spin slower.

What About the Weird Ones?

The paper notes that there is a tiny group of three Double Neutron Stars that don't fit either of these two paths. They have weird orbits that are hard to explain. The authors admit their model can't easily explain these three, suggesting they might have formed in a different way entirely (perhaps by crashing into each other in a crowded star cluster, rather than being born as a pair).

Summary of Results

The Split: Double Neutron Stars form via two main paths: one for "young" donors (leading to merging pairs) and one for "old" donors (leading to non-merging pairs).
The Energy: To get the "young" path to work, the universe needs to be very efficient at throwing off gas clouds.
The Kicks: The explosions that create these stars must be gentle, not violent, to keep the pairs together.
The Mergers: The pairs that will eventually crash and create gravitational waves (like the famous GW170817 event) almost certainly come from the "Young Path."

What the Paper Does Not Say

It does not predict when the next merger will happen.
It does not claim to explain every single weird star in the universe (specifically the three "weird" ones mentioned above).
It does not include a model for how the stars "spin down" or fade out over time, which limits how perfectly they can match real-world observations right now.

In short, this paper uses a detailed computer simulation to show that the universe has two distinct recipes for making Double Neutron Stars, and only one of those recipes produces the pairs that will eventually collide and shake the fabric of space-time.

Technical Summary: Forming Double Neutron Stars using Detailed Binary Evolution Models with POSYDON

Problem Statement
Double neutron stars (DNSs) serve as critical laboratories for understanding massive binary evolution, supernova mechanisms, and gravitational wave sources. While over two dozen DNSs have been identified in the Milky Way, their formation channels remain complex and debated. Previous population synthesis studies have largely relied on rapid binary population synthesis (BPS) codes utilizing fitting formulas or libraries of single-star tracks. These approaches often lack the self-consistent treatment of binary interactions, particularly during the common envelope (CE) phase and core-collapse supernovae (SN). Furthermore, the observed Galactic DNS population exhibits a distinct bifurcation in orbital periods and eccentricities that standard models struggle to reproduce without fine-tuning. This paper addresses the need to constrain the detailed evolutionary stages leading to DNS formation using a more rigorous, self-consistent framework.

Methodology
The authors utilize POSYDON (POpulation SYnthesis with Detailed binary-evolution simulatiONs), a state-of-the-art binary population synthesis code that integrates detailed single and binary stellar evolution models based on the 1D code MESA. Unlike rapid BPS codes, POSYDON employs pre-computed stellar evolution grids and on-the-fly calculations to handle complex physics, including magnetic braking, gravitational wave radiation, CE evolution, and supernova kicks.

The study simulates $2 \times 10^6$ binary populations across 25 distinct model variations (summarized in Table 2). Key variables manipulated include:

Common Envelope (CE) Treatments: Variations in the CE ejection efficiency ( $\alpha_{CE}$ ), the definition of the core-envelope boundary (Hydrogen mass fractions of 1%, 10%, and 30%), and the treatment of post-CE mass loss (stable mass transfer vs. wind loss vs. single-phase).
Supernova (SN) Models: Different prescriptions for core-collapse outcomes (Sukhbold et al. 2016; Patton & Sukhbold 2020; Fryer et al. 2012), distinguishing between Iron-core collapse (CCSN) and Electron-Capture Supernovae (ECSN).
Natal Kick Prescriptions: Comparison of the standard Hobbs-Maxwellian distribution ( $\sigma \approx 265$ $σ \approx 265$ km/s) against "physics-based" models:
- AsymEj: Kicks derived from asymmetric ejecta mechanisms (Janka 2017), scaling with ejecta mass.
- Linear: Kicks scaling linearly with ejecta mass (Richards et al. 2023).
- Low Maxwellian: A reduced dispersion model ( $\sigma \approx 20$ km/s).

The models assume Solar metallicity and a constant star formation history. The resulting synthetic populations are qualitatively compared to the observed sample of 25 Galactic DNSs (Table 1), focusing on orbital period ( $P_{orb}$ ) and eccentricity ( $e$ ) distributions, as well as merger rates.

Key Contributions and Results

Bifurcation of Formation Channels: The study identifies a fundamental bifurcation in the CE formation channel, determined by the evolutionary state of the donor star at the onset of Roche Lobe Overflow (RLO) with the first-born neutron star:
- He Core Channel (Case B RLO): The donor initiates mass transfer while possessing a helium core. These systems typically have shorter orbital periods and higher envelope binding energies.
- C/O Core Channel (Case C RLO): The donor evolves further to develop a carbon-oxygen core before RLO. These systems possess loosely bound envelopes and result in wider orbital periods.
- Significance: This bifurcation naturally explains the observed split in the Galactic DNS population between short-period, merging systems (Sub-population i) and wide, non-merging systems (Sub-population ii).
Constraints on Common Envelope Physics: The formation of the He core channel (which produces the merging DNSs) is highly sensitive to CE efficiency. The models reveal that this channel only survives if:
- A generous core definition is used (H-fraction $\simeq 30\%$ ), or
- A high CE ejection efficiency is adopted ( $\alpha_{CE} \gtrsim 1.2$ ).
  Without these conditions, He-core systems merge within the CE. The C/O core channel is more robust and survives even with lower $\alpha_{CE}$ values.
Natal Kick Constraints: The standard Hobbs-Maxwellian kick distribution ( $\sigma \approx 265$ km/s) is ruled out as it disrupts the majority of binaries during the first supernova, resulting in merger rates orders of magnitude lower than observed. The observed Galactic DNS population is best reproduced by low-velocity kick prescriptions ( $\lesssim 50$ km/s). The "physics-based" AsymEj and Linear models, which naturally produce lower kicks for stripped progenitors, successfully reproduce both the orbital distribution and the local merger rate density ( $O(10) - O(100)$ Gpc $^{-3}$ yr $^{-1}$ ).
Pulsar Recycling and Mass Accretion: The study finds that pulsar recycling (spinning up the first-born NS) occurs predominantly during the stable Case BB mass transfer phase with the helium-star companion, rather than during the CE phase. Systems formed via the He core channel undergo more efficient recycling, consistent with the shorter spin periods observed in Sub-population (i).
Supernova Types and Progenitors:
- The first supernova is predominantly an Iron-core collapse (CCSN).
- The second supernova often involves an ultra-stripped progenitor. Approximately 60% of merging DNS progenitors are overfilling their Roche lobes at the time of the second SN.
- Merging DNSs are born from progenitors with median pre-SN envelope masses of $\sim 0.2 M_\odot$ . This suggests that a significant fraction of merging DNSs could be observed as Type Ib or Type Ic supernovae.
Merger Rates and GW Sources: The models predict that only DNSs formed through the He core channel merge within a Hubble time. Consequently, gravitational wave events like GW170817 and GW190425 are inferred to originate from this channel. The predicted local merger rates for models with physics-based kicks align with constraints from LIGO-Virgo-KAGRA and the Galactic DNS sample.

Significance and Claims
The paper claims that detailed evolutionary models are essential for resolving the formation history of DNSs. By moving beyond fitting formulas to self-consistent grids, the authors demonstrate that the observed diversity of Galactic DNSs is a natural consequence of the donor star's evolutionary state at the onset of the CE phase.

The authors modestly note that their conclusions are based on qualitative comparisons due to the current absence of a pulsar evolution model within POSYDON (which would account for selection effects and spin-down). They also acknowledge that their results are specific to Solar metallicity. However, the study establishes that:

A bifurcation in formation channels (He-core vs. C/O-core) is a robust feature of binary evolution that explains the observed orbital period split.
Efficient CE ejection (via generous core definitions or high $\alpha_{CE}$ ) is required to form the merging population.
Low natal kicks are a necessary condition for the survival of DNS binaries, supporting recent theoretical shifts away from high-velocity kick distributions.
Unstable Case BB mass transfer does not produce merging DNSs in these detailed models, challenging previous hypotheses regarding the formation of ultra-fast merging systems.

The work provides a foundation for future quantitative Bayesian analyses once pulsar evolution models are integrated into the POSYDON framework.