Implications of low neutron star merger rates for gamma-ray bursts, r-process production and Galactic double neutron stars

This paper analyzes the latest gravitational wave data to show that the revised, lower binary neutron star merger rate creates a significant tension with the observed rates of short gamma-ray bursts, r-process element production, and Galactic double neutron star populations, suggesting potential gaps in our understanding of these astrophysical processes.

The Gist

Imagine the universe as a giant, cosmic factory. For a long time, scientists thought they knew exactly how many "products" this factory was churning out: specifically, the explosive collisions of two dead stars (neutron stars) crashing into each other.

Here is the story of a new study that suggests the factory might be producing far fewer of these collisions than we thought, and why that creates a massive puzzle for astronomers.

The Big Discovery: The "Gold" Factory

In 2017, we got our first clear look at this factory when two neutron stars smashed together (an event called GW170817). This was a "multimessenger" event, meaning we saw it in two ways:

1. Gravitational Waves: Ripples in space-time (like a shockwave in a pond).
2. Light: A burst of gamma rays and a glowing cloud of debris (a kilonova).

This single event proved two huge things:

• These crashes create Short Gamma-Ray Bursts (SGRBs), which are intense flashes of high-energy light.
• These crashes are the cosmic forges that create heavy elements (like gold, platinum, and uranium), known as "r-process" elements.

The Problem: The Factory is Slowing Down

For a few years, scientists used that one big crash to guess how often these events happen across the whole universe. Their initial guess was high: maybe thousands of crashes per year in a volume of space the size of a few billion suns. This guess fit perfectly with how many gamma-ray bursts we see and how much gold we think exists in our galaxy.

But then, the "listening devices" (gravitational wave detectors like LIGO) got much more sensitive. They listened for years, scanning the sky.

• The Result: They only found two confirmed neutron star crashes in total (GW170817 and GW190425).
• The Update: Because they found so few, the scientists had to lower their estimate. The new number is much smaller—roughly 3 to 10 times lower than the old guess.

The Puzzle: The Three Mismatched Scales

Now, the scientists (led by Maya Fishbach and her team) asked: "If the factory is producing fewer collisions, does that break our understanding of the universe?"

They compared the new, lower collision rate against three different ways we measure the universe's output. Think of it like checking a bakery's production log against three different receipts:

1. The Gamma-Ray Receipt (The "Flash" Count)

We count how many short, bright flashes of light (SGRBs) we see in the sky.

• The Conflict: Even with the new, lower collision rate, we are still seeing 3 to 18 times more gamma-ray flashes than the number of collisions should produce.
• The Analogy: Imagine the bakery says, "We only bake 10 cakes a day." But the delivery trucks are dropping off 50 boxes of cakes.
• Possible Explanations:
  • Maybe the "flash" isn't coming from a single narrow beam, but a wide, messy spray (a wider jet angle).
  • Maybe we are overcounting the flashes because some of them aren't actually from crashing stars at all.
  • Maybe the factory was much busier in the past (when the universe was younger) and is slowing down now.

2. The Gold Receipt (The "Heavy Element" Count)

We look at how much gold and platinum exist in our Milky Way galaxy.

• The Conflict: To make all the gold we see, we need a certain number of crashes. The new, lower rate is just barely enough to explain the gold, but it's on the edge. If the rate drops any further, we won't have enough collisions to explain where all our gold came from.
• The Analogy: The bakery says, "We only bake 10 cakes." But the town has 100 people who all need a cake. If the bakery really only makes 10, where did the other 90 come from?
• Possible Explanations:
  • Maybe each crash produces more gold than we thought.
  • Maybe the galaxy was "baking" much more intensely billions of years ago, so the current low rate is just the tail end of a busy era.

3. The "Waiting Room" Receipt (The "Double Neutron Star" Count)

We look at our own galaxy and count pairs of neutron stars that are currently orbiting each other, waiting to crash.

• The Conflict: Based on how many pairs we see and how fast they spiral together, our galaxy should be producing more crashes than the new gravitational wave data suggests.
• The Analogy: The bakery has a waiting list of 50 couples waiting to buy a cake. Based on how fast they usually buy, the bakery should be selling 50 cakes a day. But the sales log says they only sold 10.
• Possible Explanations:
  • Maybe we are missing some "heavy" pairs in our galaxy because they crash too fast to be seen by our radio telescopes.
  • Maybe our galaxy is just a "super-producer" compared to other galaxies.

The Bottom Line: A Cosmic Tightrope

This paper is essentially saying: "We have a problem."

The new data suggests the universe is making fewer neutron star crashes than we thought. But the evidence from light (gamma rays) and matter (gold) suggests it's making more.

What does this mean for us?
It means one of our assumptions is wrong. The scientists are now playing detective to figure out which one:

• Are the gamma-ray flashes wider than we thought?
• Is our galaxy a special case that produces more gold than others?
• Did the universe have a "golden age" of crashes billions of years ago that we aren't seeing today?

The Takeaway:
Just like a mechanic who hears a new engine sound and realizes the car might not be getting the gas mileage they thought, these astronomers are realizing their model of the universe needs a tune-up. The "low rate" of crashes is a new, tighter constraint that forces us to rethink how the universe creates its most violent explosions and its most precious metals.

The next few years of listening with gravitational wave detectors will be crucial. If they find more crashes, the tension goes away. If they find even fewer, we might have to completely rewrite the story of how the universe works.

Technical Summary

1. Problem Statement

The paper addresses a growing tension between the binary neutron star (BNS) merger rate inferred from Gravitational Wave (GW) observations and rates derived from independent astrophysical probes.

• Context: The initial BNS merger rate inferred from the single event GW170817 (110–3840 Gpc$^{-3}$ yr$^{-1}$) was consistent with the rate of short gamma-ray bursts (SGRBs), the Milky Way's (MW) r-process element mass, and the population of Galactic double neutron star (DNS) systems.
• The Issue: With the release of the latest GW catalog (GWTC-4), covering the first 8 months of the fourth observing run (O4), no new significant BNS mergers were detected. Consequently, the inferred BNS merger rate has been revised downward significantly.
• Objective: The authors investigate whether this updated, lower GW-inferred rate creates a statistical tension with the rates required to explain SGRBs, the Galactic r-process abundance, and the observed Galactic DNS population.

2. Methodology

The authors employ a multi-faceted approach combining GW data analysis with independent astrophysical constraints:

• GW Rate Inference (GWTC-4):
  • Instead of marginalizing over uncertain mass distributions (which can introduce systematics), the authors infer merger rates in specific mass bins.
  • They define three BNS mass bins based on known events: GW170817-like ($\sim$1.3 $M_\odot$), GW190425-like ($\sim$1.65 $M_\odot$), and a high-mass bin ($\sim$2.1 $M_\odot$).
  • They also estimate rates for "EM-bright" Neutron Star–Black Hole (NSBH) mergers (low-mass BHs capable of tidal disruption), specifically considering GW230529.
  • They apply conservative priors (Jeffreys prior for zero detections) and account for detection sensitivity scaling with chirp mass ($M^{5/2}$).
• SGRB Rate Comparison:
  • They compare the GW-inferred BNS rate ($R_{BNS}$) with the volumetric SGRB rate ($R_{SGRB}$) derived from cosmological observations.
  • They analyze the ratio $R_{SGRB}/R_{BNS}$ as a function of jet opening angles ($\theta_{jet}$) to determine if beaming corrections can reconcile the rates.
• r-Process Mass Calculation:
  • They recalculate the total r-process mass in the Milky Way, explicitly including the circumgalactic medium (CGM), which was often neglected in previous estimates.
  • They focus on elements with atomic number $Z \ge 56$ (above the second r-process peak) to avoid uncertainties in the first peak (s-process contamination).
  • They derive the required merger rate to produce this mass, accounting for the Milky Way's star formation history (SFH) and delay time distributions via a correction factor $C_{SFH}$.
• Galactic DNS Population:
  • They compare the GW rate to the merger rate inferred from the known population of Galactic DNS systems, applying corrections for pulsar beaming and radio survey selection effects.

3. Key Contributions & Results

A. Updated BNS Merger Rates (GWTC-4)

• Total Rate: The authors infer a total BNS merger rate of 28–300 Gpc$^{-3}$ yr$^{-1}$ (90% credibility), with a median of 110$^{+192}_{-82}$ Gpc$^{-3}$ yr$^{-1}$.
• Mass Bins:
  • GW170817-like ($\sim$1.3 $M_\odot$): 53$^{+176}_{-49}$ Gpc$^{-3}$ yr$^{-1}$.
  • GW190425-like ($\sim$1.65 $M_\odot$): 30$^{+99}_{-28}$ Gpc$^{-3}$ yr$^{-1}$.
  • High-mass ($\sim$2.1 $M_\odot$): Upper limit $< 30$ Gpc$^{-3}$ yr$^{-1}$.
• NSBH Contribution: The rate of EM-bright NSBH mergers is constrained to be low (conservatively $< 50\%$ of the BNS rate), contributing minimally to the total merger rate.

B. Tension with Short Gamma-Ray Bursts (SGRBs)

• The Discrepancy: Comparing the GW rate to the cosmological SGRB rate reveals a significant tension. The SGRB rate is 3.6 to 18 times higher than the median GW-inferred BNS rate (assuming $R_{BNS} \approx 100$ Gpc$^{-3}$ yr$^{-1}$).
• Implications:
  • Beaming: Reconciling the rates would require SGRB jets to have opening angles $\theta_{jet} \gtrsim 10^\circ$ for the majority of events. This contradicts current measurements where most jets are narrower.
  • Progenitors: It suggests that BNS mergers may only account for $\lesssim 10\%$ of observed SGRBs, implying a need for alternative progenitors (e.g., collapsars) or a high fraction of "choked" jets that fail to produce observable bursts.
  • Evolution: A steep redshift evolution (higher rates at $z \sim 0.5$ compared to $z=0$) could alleviate the tension, but this requires specific delay time distributions.

C. Tension with r-Process Production

• Revised Mass: By including the CGM, the authors estimate the MW r-process mass ($Z \ge 56$) to be 3800 $\pm$ 800 $M_\odot$, roughly double previous estimates.
• Required Rate: To produce this mass, the required NSM rate is 89–410 Gpc$^{-3}$ yr$^{-1}$.
• Status: While the lower bound of this range (89) overlaps with the GW rate, the central values suggest a tension. Reconciling them requires:
  • A higher r-process yield per merger ($>0.01 M_\odot$).
  • A steeply declining star formation history in the MW (high $C_{SFH}$).
  • Prompt merger channels (short delay times).

D. Tension with Galactic Double Neutron Stars (DNS)

• Comparison: The rate inferred from the Galactic DNS population (corrected for beaming) is 230–510 Gpc$^{-3}$ yr$^{-1}$ (or 23–51 Myr$^{-1}$ per MW).
• Tension: This is a factor of 2.3–5.1 higher than the GW-inferred rate.
• Resolution: This tension could be resolved if the MW hosts a higher merger rate than typical galaxies, or if current corrections for pulsar beaming and survey completeness are overestimated (i.e., pulsars are less beamed or surveys are more complete than assumed).

4. Significance and Physical Implications

The paper highlights a critical "rate crisis" in multimessenger astrophysics. The downward revision of the GW-inferred BNS rate challenges the standard model where BNS mergers are the sole or dominant source of SGRBs and Galactic r-process elements.

• Constraints on Physics: The tension forces a re-evaluation of:
  • Jet Physics: SGRB jets may be wider than previously thought, or a significant fraction of mergers produce "failed" jets.
  • Stellar Evolution: The delay time distribution of BNS mergers must be steep (prompt mergers) to explain r-process enrichment, which conflicts with some population synthesis models.
  • Galactic Dynamics: The Milky Way may be an outlier with a higher merger rate than the average galaxy, or our understanding of pulsar beaming in the Galactic DNS population is incomplete.
• Future Outlook: The authors note that if the O4 run concludes with zero new BNS detections, the rate will drop further (to $\sim$16–170 Gpc$^{-3}$ yr$^{-1}$), exacerbating these tensions. Future observations of kilonovae, SGRB afterglows, and deeper GW runs are essential to resolve whether the discrepancy lies in the GW rate, the astrophysical models, or the selection effects of the independent probes.

In summary, the paper argues that the low BNS merger rate inferred from GWTC-4 is increasingly difficult to reconcile with independent astrophysical evidence, suggesting that our understanding of merger physics, jet structure, or Galactic chemical evolution requires significant revision.