Binary Neutron Star Merger Evolution and r-Process Enrichment in the Milky Way Disk

TL;DR

This study investigates whether binary neutron star mergers with evolving r-process enrichment efficiency can explain the observed heavy element abundance in the Milky Way, using multi-messenger astronomical data and statistical analysis.

Contribution

It demonstrates that models with evolving BNS merger efficiency are strongly favored and provides a quantitative framework to assess their role in galactic chemical evolution.

Findings

Evolving BNS merger scenarios are strongly preferred over non-evolving ones.

Evolved scenarios are in tension with short gamma-ray burst observations.

Evolved scenarios are consistent with gravitational-wave background constraints.

Abstract

The origin of half of the rapid neutron-capture nucleosynthesis (r-process) elements in the Universe remains an open question. Binary neutron star (BNS) mergers have been shown to face difficulties in reproducing the observed r-process enrichment in Milky Way disk stars. However, their r-process enrichment efficiency may evolve with redshift beyond the star formation history, potentially due to evolution in the merger rate or the average r-process yield in the early Universe. In this paper, we explore the possibility that BNS mergers with an evolving enrichment efficiency could serve as the sole r-process production channel. By jointly comparing gravitational-wave observations from LIGO--Virgo--KAGRA, short gamma-ray bursts, Galactic neutron star populations, and stellar abundance measurements in Milky Way disk stars, we find that scenarios with additional evolution are strongly preferred over non-evolving scenarios, with Bayes factors exceeding $10^{20}$. We quantify the required evolution in both the merger rate and yield, and directly compare them with observations and theoretical predictions. We find that the evolved scenarios are in tension with short gamma-ray burst observations and predictions from multiple population synthesis models, while remaining consistent with current stochastic gravitational-wave background constraints. Our results provide a quantitative framework for evaluating whether BNS mergers with evolving enrichment efficiency can account for the observed r-process enrichment history of Milky Way disk stars.