Mergers Fall Short: Non-merger Channels Required for Galactic Heavy Element Production

By integrating gravitational-wave, gamma-ray burst, and stellar abundance data within a unified likelihood framework, this study demonstrates that neither neutron star-black hole nor fast-merging binary neutron star mergers can serve as the dominant additional source of Galactic r-process elements, thereby confirming the necessity of non-merger production channels.

The Gist

The Big Mystery: Where Did the Heavy Stuff Come From?

Imagine the Universe as a giant kitchen. For a long time, we knew how to make the "light" ingredients like hydrogen and helium (the flour and sugar of the cosmos). But the "heavy" ingredients—like gold, platinum, and the rare earth element Europium (which we'll call Eu for short)—were a mystery.

Scientists knew these heavy elements are made in violent, neutron-rich explosions. For a while, the leading theory was that merging neutron stars (two dead stars crashing into each other) were the main chefs in this kitchen. This was confirmed when we heard the "crash" of two neutron stars (GW170817) and saw the resulting "kitchen firework" (a kilonova) that created heavy elements.

But there's a problem.

If neutron star mergers were the only chefs, the "recipe" for our galaxy (the Milky Way) wouldn't taste right. When we look at the oldest stars in our galaxy's disk, they have a lot of heavy elements (Europium) even when they have very few "iron" ingredients.

Think of it like baking a cake:

• Iron is the flour.
• Europium is the chocolate chips.
• Time is the baking process.

If you only add chocolate chips after you've baked a huge amount of flour (which is what the standard "delayed merger" theory predicts), your early cakes would have no chocolate. But the old stars (the early cakes) are full of chocolate. This means something else must have been adding chocolate chips immediately after the kitchen opened, long before the standard mergers could get to work.

The Suspects: Who Else Could Be Cooking?

Scientists proposed two new suspects to explain the early chocolate chips:

1. Fast-Merging Neutron Stars: Neutron stars that crash into each other almost instantly after being born, rather than waiting millions of years.
2. Neutron Star-Black Hole Mergers: A neutron star crashing into a black hole, which might be a more efficient "chocolate maker."

The authors of this paper asked: "Can these two suspects be the main reason our galaxy has so much heavy stuff?"

The Investigation: A Statistical Taste Test

The researchers built a super-advanced computer model (a "taste test") that simulates the history of the Milky Way. They fed it data from:

• Gravitational Waves: The "sound" of stars crashing (from LIGO/Virgo).
• Gamma-Ray Bursts: The "flash" of the explosion.
• Star Spectra: The actual chemical fingerprints of stars in our galaxy.

They didn't just guess; they used a rigorous statistical method to see if the data could possibly fit if these suspects were the main chefs.

The Verdict: The Suspects Are Innocent (of the Main Crime)

The results were surprising and definitive. The paper concludes that neither of these merger types can be the primary source of heavy elements in our galaxy. Here is why, broken down by scenario:

Scenario 1: The "Detectable" Fast Mergers

• The Theory: Maybe there are tons of neutron stars crashing instantly, and our current detectors can hear them all.
• The Result: To explain the heavy elements in old stars, the model says 98% of all neutron star crashes would have to be "fast" ones.
• The Analogy: It's like saying that 98% of all car accidents in the US happen within 5 minutes of the car being built. While possible, it contradicts everything we know about how cars (and stars) are built and driven. It's just too extreme to be true.

Scenario 2: The "Undetectable" Fast Mergers

• The Theory: Maybe there are so many fast crashes happening that our detectors are missing them (perhaps because they are too weird or quiet).
• The Result: To explain the heavy elements, the model says there would need to be 17 times more of these invisible crashes than the ones we actually see.
• The Analogy: Imagine you see 10 people sneezing in a room. To explain the smell of a rotting fish, you have to assume there are 170 other people sneezing that you can't see or hear. It's a huge stretch. It requires a "ghost population" of stars that is theoretically very unlikely to exist in such massive numbers.

Scenario 3: The Black Hole Partners

• The Theory: Maybe neutron stars crashing into black holes are the secret sauce.
• The Result: Even if we give black holes the benefit of the doubt (assuming they crash instantly and make lots of heavy stuff), they still can't do the job alone. They might contribute a little bit, but they can't be the main source.

The Real Conclusion: We Need a New Chef

Since the "merger" suspects (both fast and slow, both star-star and star-black hole) can't explain the data, the paper concludes that we are missing a major ingredient.

There must be a non-merger channel.

• The Analogy: We've been looking for the chocolate chips in the "crash" bin. But the evidence suggests the chocolate chips are actually coming from a completely different machine—perhaps a specific type of massive star explosion (like a magnetorotational supernova) or a "collapsar" (a star collapsing into a black hole) that we haven't fully accounted for yet.

Why This Matters

This paper is a "reality check" for astrophysics. It tells us that while neutron star mergers are real and important, they aren't the whole story. If we want to understand how the heavy elements in our jewelry and electronics were made, we need to look for a different kind of cosmic event that happens quickly and frequently in the early history of our galaxy.

In short: The "crash" theory is part of the recipe, but it's not the whole cake. We need to find the other chef.

Technical Summary

1. Problem Statement

The origin of the heaviest elements in the Universe, synthesized via the rapid neutron-capture process (r-process), remains a central question in astrophysics. While the detection of the binary neutron star (BNS) merger GW170817 confirmed that such mergers produce r-process elements, a significant discrepancy exists between theoretical models of BNS mergers and observational data from the Milky Way.

• The Discrepancy: Observations of metal-poor stars in the Milky Way disk reveal a high abundance of Europium (Eu, a proxy for r-process elements) even at low metallicity ([Fe/H]). This implies that r-process enrichment began very early in the Galaxy's history, "promptly" following star formation.
• The Conflict: Standard BNS merger models predict a "delayed" population, where mergers occur significantly after star formation due to gravitational wave inspiral timescales. Chemical evolution models using these delayed rates fail to reproduce the observed early enrichment and the specific negative slope of the [Eu/Fe] versus [Fe/H] trend seen in disk stars.
• The Hypothesis: To resolve this, it has been proposed that additional channels exist, specifically:
  1. Fast-merging BNS systems: Systems that merge almost immediately after formation (e.g., via dynamical interactions in dense clusters).
  2. Neutron Star–Black Hole (NSBH) mergers: Which may be more efficient ejectors or have different delay times.
• The Goal: This paper rigorously tests whether these merger-based channels can quantitatively account for the required r-process enrichment in the Milky Way without contradicting existing gravitational-wave (GW) and electromagnetic observations.

2. Methodology

The authors employ a unified, likelihood-based Bayesian inference framework that integrates multiple data sources to constrain r-process production channels.

• Chemical Evolution Model:

  • Utilizes a one-zone model of the Milky Way disk, treating the interstellar medium (ISM) as a well-mixed reservoir.
  • Tracks the evolution of [Fe/H] and [Eu/Fe] over cosmic time.
  • Incorporates standard sources for Iron (Type Ia and Core-Collapse Supernovae) and tests r-process contributions from various merger populations.
  • Key Assumption: The r-process yield per event is derived from component masses, spins, and the Equation of State (EOS), using semi-analytical fits calibrated to numerical relativity simulations.

• Data Integration:

  • Stellar Abundances: Uses [Eu/Fe] vs. [Fe/H] data from Milky Way disk stars (Ref. [52]).
  • Gravitational Waves: Incorporates population inference results from GWTC-4.0 (LIGO-Virgo-KAGRA) for local merger rates and mass distributions.
  • Electromagnetic Constraints: Uses short Gamma-Ray Burst (sGRB) host galaxy observations to constrain delay-time distributions (DTD) for BNS mergers.
  • Pulsar Data: Uses constraints on neutron star radii and compactness to refine ejecta mass estimates.

• Inference Framework:

  • Instead of using point estimates, the authors propagate uncertainties in merger rates, DTDs, ejecta yields, and stellar measurements through a full Bayesian likelihood analysis.
  • They marginalize over population parameters to determine the probability space where merger channels can satisfy the chemical evolution constraints.

3. Key Contributions

• Unified Likelihood Framework: The study moves beyond "point-estimate" comparisons by rigorously propagating uncertainties across the entire parameter space (rates, yields, delays, and observational errors).
• Testing "Hidden" Populations: The paper explicitly models scenarios where fast-merging BNS systems might be undetectable by current GW pipelines (e.g., due to high eccentricity or short in-band signals), testing if a "hidden" population could solve the enrichment problem.
• Quantitative Constraints on NSBH: It provides the first rigorous, likelihood-based assessment of NSBH mergers as a secondary channel alongside BNS mergers, rather than treating them in isolation.

4. Results

The analysis evaluates three specific scenarios, all of which fail to explain the observed r-process enrichment using merger channels alone:

A. Scenario 1: GW-Detectable Fast-Merging BNS

• Assumption: Fast-merging BNS systems are detectable by current GW searches and follow the same detection efficiency as delayed systems.
• Result: To reproduce the observed [Eu/Fe] trend, the model requires that >98% of all BNS mergers be "fast-merging" (merging with zero delay relative to star formation).
• Implication: This implies a local volumetric rate for fast mergers of $\sim 93$ Gpc$^{-3}$ yr$^{-1}$, vastly dominating the delayed population. This contradicts theoretical expectations for the formation of such systems and creates tension with the observed sGRB delay-time distribution.

B. Scenario 2: GW-Undetectable Fast-Merging BNS

• Assumption: Fast-merging BNS systems exist but are missed by current GW searches (e.g., due to eccentricity).
• Result: The model requires the rate of these "hidden" fast mergers to be $\sim 17$ times higher than the rate of delayed BNS mergers in the local Universe (rising to $\sim 27$ times at $z=2$).
• Implication: Even if invisible to GW detectors, the sheer number of events required to explain the chemical enrichment is physically implausible given current understanding of stellar dynamics in globular clusters and galactic nuclei.

C. Scenario 3: NSBH Mergers as a Secondary Channel

• Assumption: A combination of standard delayed BNS mergers and NSBH mergers (with prompt delays) contributes to enrichment.
• Result: Even under optimistic assumptions (NSBH merging promptly), the combined contribution of BNS and NSBH mergers is limited to <10% of the total r-process budget.
• Implication: More than 90% of the r-process enrichment must come from a third, unidentified source that tracks the Star Formation Rate (SFR) directly.

5. Significance and Conclusion

• Exclusion of Merger-Dominance: The paper definitively concludes that neither fast-merging BNS systems nor NSBH mergers can serve as the dominant source of r-process elements in the Milky Way disk without generating severe tension with existing GW, sGRB, and stellar abundance data.
• Necessity of Non-Merger Channels: The results strongly imply that the "missing" prompt enrichment source must be a non-merger channel. Candidates include magnetorotational supernovae, collapsars, or giant flares from magnetars, which can produce r-process elements promptly after star formation.
• Robustness: The conclusions hold even when varying the Star Formation Rate (SFR) model (e.g., assuming a constant SFR) or assuming significantly higher yields for fast-merging systems (e.g., 2x the yield of delayed systems).
• Theoretical Challenge: The findings pose a major challenge to theories suggesting that dynamical formation in dense clusters is the primary driver of early r-process enrichment. It suggests that the "prompt" r-process source is likely a distinct astrophysical phenomenon from compact binary mergers.

In summary, while mergers contribute to the r-process, they are insufficient to explain the chemical history of the Milky Way disk. The search for the primary source of early heavy elements must look beyond binary mergers.