🔭 astrophysics

Constraining Binary Neutron Star Populations using Short Gamma-Ray Burst Observations

By analyzing 64 binary neutron star population models against Fermi-GBM short gamma-ray burst observations, this study demonstrates that only local merger rates of approximately 100 Gpc $^{-3}$ yr $^{-1}$ combined with plausible jet geometries can successfully account for the observed sGRB population, thereby ruling out lower-rate models as sole progenitors.

Original authors: Alessio Ludovico De Santis, Samuele Ronchini, Filippo Santoliquido, Marica Branchesi

Published 2026-07-22

📖 6 min read🧠 Deep dive

CC BY 4.0

Original authors: Alessio Ludovico De Santis, Samuele Ronchini, Filippo Santoliquido, Marica Branchesi

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 Cosmic Detective Story: Chasing Ghosts and Jets

Imagine the universe as a giant, dark ocean. For a long time, we only knew about the lighthouses—the bright, steady stars. But in the last decade, we've learned to listen for the splashes. When two incredibly dense, city-sized stars called neutron stars crash into each other, they create a ripple in space and time itself, known as a gravitational wave. This is like hearing a splash in the dark. But sometimes, that splash is accompanied by a blinding flash of light, a short gamma-ray burst (sGRB), which is like a massive firework exploding in the sky.

Scientists have long suspected that these fireworks are caused by the crashing stars. A famous event in 2017, where we heard the splash and saw the flash at the same time, seemed to prove it. But here's the puzzle: we can only hear a few splashes with our current listening devices, while we see thousands of fireworks. The big question is: Are the crashing stars the only thing making these fireworks? Or are there other secret firework-makers we haven't found yet? This paper is a massive detective story trying to solve that mystery by counting the stars, measuring the fireworks, and checking if the math adds up.

The Great Cosmic Firework Count

Think of binary neutron stars as a pair of cosmic dancers who eventually collide. When they crash, they might launch a super-fast jet of energy—a "relativistic jet"—that shoots out like a laser beam. If that laser points right at Earth, we see a Short Gamma-Ray Burst (sGRB), a brilliant flash of high-energy light. If the laser points away, we see nothing, even though the crash happened.

The authors of this paper, A. L. De Santis and their team, wanted to know: Do the number of crashes we predict match the number of flashes we actually see?

To find out, they didn't just guess. They built a massive digital simulation lab. They took 64 different theories about how these neutron star couples form and evolve over billions of years. Some theories say the stars are easy to pair up; others say it's a rare accident. They also had to guess how the "lasers" (jets) behave. Do they all look the same? Are they narrow beams or wide cones?

They ran their simulations through a "Monte Carlo" framework, which is basically a fancy way of saying they ran the experiment thousands of times with random variations to see what the most likely outcome is. They then compared their simulated fireworks to the real data collected by the Fermi-GBM satellite over the last 16 years.

The "Jet Fraction" Mystery

The key to the mystery is a number the authors call the jet fraction ( $f_j$ ). Imagine you have 100 neutron star crashes. If 100 of them shoot a laser at Earth, $f_j$ is 1 (or 100%). If only 10 do, $f_j$ is 0.1.

Here is the catch: $f_j$ cannot be more than 1. You can't have 150% of the crashes shooting lasers. If your math says you need 150% of the crashes to explain the number of fireworks we see, then your theory is broken. It means you are missing a lot of crashes, or there are other things making the fireworks.

The Findings: Narrow Beams and High Rates

The team tested three different shapes for the laser beams:

The Structured Jet: A beam that is super bright in the middle and fades out at the edges (like the one seen in the famous 2017 event).
The Top-Hat Jet: A simple, uniform cone of light (like a flashlight beam).
The Non-Universal Jet: A mix where every crash might have a different beam width.

The Result:
The paper found a clear split in the universe of theories.

The "Low-Rate" Theories: Some models predicted that neutron stars crash relatively rarely (fewer than 50 crashes per cubic gigaparsec per year). When the authors tried to make these rare crashes explain all the fireworks we see, the math broke. To match the observations, these models required a jet fraction greater than 1. In plain English, they would need more lasers than there are crashes. This is impossible. The paper effectively rules out these low-rate models if neutron star crashes are the only source of short gamma-ray bursts.
The "High-Rate" Theories: The models that survived the test were the ones predicting a higher crash rate, around 100 crashes per cubic gigaparsec per year. These models fit the data perfectly, requiring a jet fraction of about 0.7 to 0.8. This means that roughly 70-80% of the time, when two neutron stars crash, they successfully launch a jet. This is a very plausible, physical number.

The Beam Width Problem

The authors also played a game of "what if" regarding the shape of the laser beams. What if the beams were much wider than we think?

They found that if the crash rate is low, the only way to save the theory is to assume the laser beams are incredibly wide (opening angles of 15° to 20° or more). However, this contradicts what we actually see. Observations of the afterglows (the fading light after the flash) suggest the beams are actually quite narrow, with a median opening angle of about 6°.

So, the paper concludes: You can't have your cake and eat it too. You can't have a low number of crashes and narrow beams. If the beams are narrow (as observations suggest), there must be a high number of crashes (around 100 per unit volume) to explain the fireworks we see.

The Verdict

This study suggests that the universe is a busy place. The neutron stars are crashing more often than the "low-rate" theories predicted. The "high-rate" models, which assume the stars are efficient at merging, are the ones that line up with reality.

The paper doesn't claim to have solved every mystery, but it draws a very strong line in the sand. It tells us that if we want to explain the short gamma-ray bursts we see, we need a population of neutron star mergers that is robust and frequent. If future gravitational wave detectors find that the crash rate is actually much lower than 100, then we will have to admit that neutron stars aren't the only firework-makers in town, and we'll need to look for other cosmic culprits. But for now, the evidence points to a universe where these stellar collisions are happening often, launching powerful jets that light up the dark.

Technical Summary: Constraining Binary Neutron Star Population Synthesis Models using Short Gamma-Ray Burst Data

Problem Statement
The landmark multi-messenger detection of the binary neutron star (BNS) merger GW170817 and its associated short gamma-ray burst (sGRB) GRB 170817A confirmed that BNS mergers can produce sGRBs. However, the limited number of BNS detections by current gravitational-wave (GW) observatories leaves open the question of whether BNS mergers alone can account for the full observed sGRB population. A critical uncertainty in linking these populations is the "jet fraction" ( $f_j$ ), defined as the fraction of BNS mergers that successfully launch a relativistic jet capable of breaking out of the merger ejecta. Previous constraints on $f_j$ have varied widely, from a few tens of percent to 100%, depending on assumptions about jet geometry and the intrinsic BNS merger rate. This work addresses the tension between theoretical BNS population synthesis models and the observed sGRB population, specifically testing whether low-rate BNS models can reproduce the Fermi-GBM sGRB catalog without invoking non-physical jet fractions or unrealistic jet geometries.

Methodology
The authors employ a comprehensive, end-to-end framework connecting 64 distinct BNS population synthesis models to a 16-year archive of sGRB data from the Fermi-GBM instrument (July 2008 – May 2025). The analysis utilizes a Markov Chain Monte Carlo (MCMC) approach implemented in the MAGGPY code to simulate synthetic sGRB catalogs and infer the required jet fraction ( $f_j$ ) for each BNS model.

BNS Population Models: The study utilizes 16 base models from Iorio et al. (2023), generated using the sevn binary evolution code and cosmoRate for cosmic merger rates. These models explore uncertainties in stellar evolution, specifically varying the common-envelope efficiency ( $\alpha_{CE} \in \{0.5, 1, 3, 5\}$ ), natal kick distributions, and mass transfer stability criteria. This results in 64 unique BNS populations with varying local merger rates ( $R_{BNS}(0)$ ) and redshift evolutions.
sGRB Emission Models: Three distinct jet geometry scenarios are tested to assess the robustness of constraints:
- Universal Structured Jet: Calibrated to the GW170817/GRB 170817A afterglow, featuring a structured angular profile with a core opening angle $\theta_c \approx 3.4^\circ$ .
- Universal Top-Hat Jet: A simplified model with a uniform emission cone and a fixed core opening angle $\theta_c$ .
- Non-Universal Top-Hat Jets: Models where $\theta_c$ varies across the population according to either a Flat (uniform) or Log-Normal distribution.
Statistical Framework: The MCMC sampling adjusts free parameters (including the jet fraction $f_j$ , spectral indices, and geometric parameters) to maximize the likelihood of reproducing the observed Fermi-GBM catalog. The likelihood function combines a Poisson term for the event rate and a two-sample Cramér-von Mises (CvM) test to compare the distribution shapes of observables (fluence, peak flux, peak energy, and duration) between synthetic and real data.
Viability Metric: A BNS population is deemed "physical" if the posterior median of the inferred jet fraction satisfies $f_j \leq 1$ . Values of $f_j > 1$ indicate that the BNS population is insufficient to produce the observed sGRB rate, even if every merger launched a jet.

Key Results
The analysis yields the following constraints on BNS populations and jet physics:

Exclusion of Low-Rate Models: Under the assumption of a universal structured jet (calibrated to GW170817), BNS population models predicting local merger rates $R_{BNS}(0) \lesssim 50 \text{ Gpc}^{-3} \text{ yr}^{-1}$ are effectively ruled out as the sole progenitors of sGRBs. These models require non-physical jet fractions ( $f_j > 1$ ) to match the observed sGRB count.
Preferred Merger Rates: Models with local merger rates around $R_{BNS}(0) \approx 100 \text{ Gpc}^{-3} \text{ yr}^{-1}$ successfully reproduce the observed sGRB population with physically plausible jet fractions ( $f_j \approx 0.7 - 0.8$ ). This rate is consistent with current GWTC-4 constraints.
Jet Geometry Tension: When relaxing the structured jet assumption to a universal top-hat model, low-rate BNS populations can only be reconciled with observations if the jet core opening angles are significantly wider ( $\theta_c \gtrsim 15^\circ - 20^\circ$ ). This creates a tension with observational inferences from sGRB afterglows, which suggest a median core angle of $\langle \theta_c \rangle \approx 6^\circ$ .
Non-Universality: Relaxing the assumption of universality (allowing a distribution of jet angles) offers limited relief. While a Log-Normal distribution of angles can lower the required median $f_j$ for low-rate models by populating the tail with wide jets, it necessitates that a substantial fraction ( $>10\%$ ) of the population possesses wide jets ( $\theta_c > 15^\circ$ ), which remains in tension with the scarcity of such wide-jet observations.
Delay Times: The analysis confirms that the favored BNS populations (those with $R_{BNS}(0) \approx 100 \text{ Gpc}^{-3} \text{ yr}^{-1}$ ) exhibit short delay times ( $\langle \tau_d \rangle < 1 \text{ Gyr}$ ), consistent with recent findings correcting for selection effects in flux-limited samples.

Significance and Conclusions
The paper establishes a robust lower bound on the local BNS merger rate derived purely from electromagnetic observations of sGRBs. By demonstrating that low-rate BNS models require either non-physical jet fractions or jet geometries inconsistent with afterglow data, the study provides a critical multi-messenger constraint on binary evolution theory.

The results favor BNS population synthesis models characterized by standard or high common-envelope efficiencies ( $\alpha_{CE} \geq 1$ ) and moderate natal kicks, which naturally yield higher merger rates. Conversely, models with low ejection efficiencies or high natal kicks are disfavored. The inference of high jet production efficiencies ( $f_j \approx 0.8$ ) suggests that the conditions for launching relativistic jets are met in the majority of BNS mergers, disfavoring scenarios where prompt collapse to a black hole systematically suppresses jet formation.

Ultimately, this work highlights the power of combining GW constraints (which currently set an upper bound on the BNS rate) with sGRB population statistics (which set a lower bound) to progressively shrink the allowed parameter space for binary neutron star evolution and jet physics.