🔭 astrophysics

Double Neutron Star Delay Times Across Cosmic Metallicities: The Role of Helium Star Progenitors

Using the POSYDON binary evolution code, this study demonstrates that the delay time distribution of double neutron star mergers is fundamentally shaped by metallicity-dependent helium star radii, resulting in a characteristic peak between 80–250 Myr with a complex, often double-peaked structure that explains r-process enrichment and short gamma-ray bursts across diverse cosmic environments.

Original authors: Abhishek Chattaraj, Jeff J. Andrews, Max Briel, Tassos Fragos, Seth Gossage, Vicky Kalogera, Philipp M. Srivastava, Elizabeth Teng

Published 2026-05-06

📖 5 min read🧠 Deep dive

CC BY 4.0

Original authors: Abhishek Chattaraj, Jeff J. Andrews, Max Briel, Tassos Fragos, Seth Gossage, Vicky Kalogera, Philipp M. Srivastava, Elizabeth Teng

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: Timing the Cosmic Dance

Imagine the universe as a giant dance floor. Every time two massive stars are born together as a pair, they eventually die and turn into neutron stars (super-dense city-sized balls of matter). Sometimes, these two neutron stars stay in a pair and spiral toward each other until they crash, creating a massive explosion called a "kilonova."

This paper asks a simple but tricky question: How long does it take for these pairs to crash after they are born?

Scientists call this the "delay time." If they crash immediately, they can enrich young galaxies with heavy elements. If they wait billions of years, they might explain explosions in very old galaxies. The authors wanted to know if the "metal content" of the stars (how much heavy stuff like iron they have) changes this waiting time.

The Main Character: The Helium Star

To understand the timing, the authors focused on a specific stage in the stars' lives: when one star has already died, and the other has become a Helium star (a star that has burned off its hydrogen fuel and is now just a hot, helium core).

Think of this Helium star as a balloon.

High Metallicity (Metal-rich stars): These stars are like balloons filled with thick, sticky air. Because of their composition, they puff up and get very large.
Low Metallicity (Metal-poor stars): These stars are like balloons filled with thin, light air. They stay small and compact.

The "Room Size" Rule

The authors discovered a simple rule: The size of the Helium star sets the minimum size of the dance floor.

When the Helium star is about to explode into a second neutron star, it needs room to breathe. It cannot be smaller than the star itself.

If the Helium star is big (high metallicity), the two stars must stay far apart.
If the Helium star is small (low metallicity), the two stars can get very close together.

Why does distance matter? Because gravity works like a rubber band. The closer the two neutron stars are when they are born, the faster they spiral together and crash. The farther apart they are, the longer they take to meet.

The "Kick" Factor

When the Helium star explodes, it gives the new neutron star a "kick" (like a sudden shove).

The authors checked if a lucky, perfectly aimed kick could push the stars so close together that they crash almost instantly (in just a few million years).
The Result: They found that while a kick can shrink the orbit, the physics of the explosion usually means the stars end up at least as far apart as the size of the Helium star was before it exploded. You can't cheat the "room size" rule easily.

The Findings: When Do They Crash?

Using powerful computer simulations, the authors tracked thousands of these star pairs across different metal environments. Here is what they found:

There is a "Minimum Wait Time": No matter how you tweak the physics, these pairs cannot crash sooner than about 40 million years after the stars were born. It takes that long for the stars to grow, die, and form the pair in the first place.
The "Sweet Spot": The most common time for these pairs to crash is between 80 and 250 million years. This is the peak of the distribution.
The "Long Waiters": About 20% of these pairs are slow pokes. They take more than 1 billion years to crash. This explains why we see these explosions in very old, metal-poor galaxies.
The "Fast Trackers": About 15% crash quickly (within 80 million years). This is fast enough to explain why some very young, metal-poor galaxies have heavy elements.

A Double-peak Surprise

The authors also found that the "waiting time" graph isn't a simple hill. At certain metal levels, it looks like a double-humped camel.

This happens because the stars can take two different evolutionary paths (like taking a highway vs. a backroad) to become neutron stars.
At some metal levels, both paths are popular, creating two distinct groups of crash times. At other levels, only one path is used.

The Bottom Line

The paper concludes that the "metal content" of a star acts like a thermostat for its size, which in turn sets the speed limit for how fast two neutron stars can crash.

Metal-poor stars stay small, get close, and crash relatively quickly (though still not instantly).
Metal-rich stars puff up, stay far apart, and generally take longer to crash.

This helps astronomers understand where and when to look for these cosmic collisions and explains how heavy elements (like gold and platinum) are spread across the history of the universe.

Technical Summary: Double Neutron Star Delay Times Across Cosmic Metallicities

Problem Statement
The delay time distribution (DTD) of double neutron star (DNS) mergers—the time elapsed between binary formation and merger—is a critical parameter for interpreting the contribution of DNS mergers to r-process nucleosynthesis, short gamma-ray bursts (sGRBs), and gravitational wave event rates. While previous studies have constrained the DTD using observational data from r-process enriched stars and sGRB host galaxies, the underlying theoretical shape of the DTD remains uncertain. Specifically, it is unclear how the metallicity-dependent evolution of the progenitor binary, particularly the helium (He) star phase, influences the birth orbital separation and subsequent merger timescales. Existing models often assume simplified power-law distributions for DTDs without fully accounting for the complex, metallicity-dependent physics of isolated binary evolution.

Methodology
This work employs a theory-driven approach combining detailed stellar evolution modeling with full-scale population synthesis to characterize the DNS DTD across a cosmological range of metallicities ( $Z = [10^{-4}, 2] Z_{\odot}$ ).

Stellar Evolution Modeling: The authors utilized the 1-D stellar evolution code MESA to model the evolution of single He stars and He star–NS binaries. Key physical processes examined include:
- Radiative Opacity: Analyzing how metallicity affects envelope opacity and subsequent radial expansion during the He main-sequence (HeMS).
- Mass Transfer (MT): Modeling the Roche-lobe overflow (RLO) phase (Case BB MT) and the impact of wind mass loss on orbital evolution. Wind mass loss rates were modeled using the Wolf–Rayet prescription, with a specific adjustment to address observational constraints on stripped stars (Y. Götberg et al. 2023).
- Supernova (SN) Dynamics: Computing post-SN orbital distributions assuming isotropic natal kicks of varying magnitudes to determine how kicks modify eccentricity and separation.
Population Synthesis: The authors used the POSYDON binary evolution suite to simulate the complete evolution of massive binaries from the Zero-Age Main Sequence (ZAMS) to DNS merger.
- Grids: Simulations utilized detailed grids for H-rich binaries, compact object–H-rich binaries, and compact object–He-rich binaries across eight metallicities.
- Physics Variations: Four distinct population models were generated to test robustness against binary physics uncertainties, varying the common envelope (CE) ejection efficiency ( $\alpha_{CE}$ ) and natal kick prescriptions (Asymmetric Ejecta vs. Lognormal distribution).
- Wind Correction: To address overestimates in wind mass loss for low-mass He stars in standard prescriptions, the authors artificially constrained the He star–NS phase evolution at all metallicities to follow the track of the $10^{-4} Z_{\odot}$ grid, effectively limiting unphysical orbital expansion.

Key Contributions and Results

He Star Radius as a Determinant: The study establishes that the radius of the He star during the HeMS phase sets a lower limit on the DNS birth orbital separation. Metallicity modulates this radius via opacity-driven envelope inflation; higher metallicities lead to larger radii and, consequently, wider birth orbits.
Merger Time Scaling: Since the inspiral time scales as $t_{insp} \propto a^4$ , even modest metallicity-driven differences in He star radii (e.g., a ~~5% difference between $Z_{\odot}$ and $10^{-4} Z_{\odot}$ ) translate into significant (~~20%) variations in merger times.
Robust DTD Features: Despite variations in binary physics assumptions (CE efficiency, kick models), the population synthesis reveals two robust features in the DNS DTD:
1. Minimum Delay Time: A hard lower bound exists where the bulk of DNS mergers occur no earlier than $\simeq 40$ Myr after star formation. The absolute minimum formation time (ZAMS to DNS birth) is $\simeq 20$ Myr, dictated by the nuclear timescale of the primary H-rich star.
2. Peak Delay Time: The DTD peaks strongly between 80 and 250 Myr across all metallicities.
Metallicity-Dependent Structure: The DTD shape is complex. At intermediate metallicities ( $Z \sim 10^{-2} Z_{\odot}$ ), the dominant CE formation channel bifurcates into two subchannels (He core and C/O core), imprinting a double-peaked structure on the DTD. At low metallicities ( $Z \lesssim 10^{-2} Z_{\odot}$ ), only the C/O core channel contributes, while at solar and super-solar metallicities, only the He core channel produces merging DNSs.
Natal Kicks: While ideal kick orientations can theoretically produce very short merger times, the authors find that the kick magnitudes required to produce the observed short-delay population are inconsistent with current observational constraints. Typical kick models (e.g., AsymEj) rarely produce the extreme orbital shrinkage necessary for sub-20 Myr delays.
Merger Fractions:
- $\approx 15\%$ of DNSs merge within 80 Myr, potentially explaining r-process enrichment in environments with brief star formation histories.
- $\gtrsim 20\%$ merge on delay times $> 1$ Gyr, providing a mechanism for sGRBs in old, metal-poor galaxies.
- At low metallicities ( $Z \lesssim 10^{-2} Z_{\odot}$ ), $>70\%$ of DNSs merge within a Hubble time.

Significance and Claims
The authors claim that their work provides a detailed, physically motivated link between DNS delay times and the binary evolution of their progenitor systems. By demonstrating that the DTD is not a simple power law but possesses a complex, metallicity-dependent structure with a robust peak at 80–250 Myr, the study offers a refined framework for assessing the contribution of DNS mergers to cosmic r-process enrichment and transient events.

The paper explicitly states that these models offer a foundation for future studies to convolve the derived DTDs with realistic cosmic star formation histories to predict NS-NS merger rates, sGRB rates, and kilonova transients across cosmic time. The authors caution that while their models suggest DNS mergers can plausibly contribute to enrichment in low-metallicity environments, this is not definitive proof without convolving with specific star formation histories. They also note limitations regarding the accuracy of low-mass He star wind modeling and the need for further testing of binary physics uncertainties, maintaining a modest stance on the exact magnitude of variations in the DTD shape due to these factors.