🔭 astrophysics

Kilonova and progenitor properties of merger-driven gamma-ray bursts

This study employs a novel Bayesian framework within the NMMA to simultaneously model kilonova and afterglow emissions for a sample of merger-driven gamma-ray bursts, successfully inferring progenitor properties (favoring binary neutron star systems for most events) and establishing correlations between ejecta mass, jet energy, and tidal deformability.

Original authors: P. Singh, G. Stratta, A. Rossi, P. T. H. Pang, M. Bulla, F. Ragosta, A. De Rosa, D. A. Kann, F. Cogato

Published 2026-07-30

📖 6 min read🧠 Deep dive

CC BY 4.0

Original authors: P. Singh, G. Stratta, A. Rossi, P. T. H. Pang, M. Bulla, F. Ragosta, A. De Rosa, D. A. Kann, F. Cogato

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

Imagine the universe as a cosmic construction site where the heaviest elements in existence—like the gold in your jewelry or the iodine in your medicine—are forged in the most violent fires imaginable. These fires aren't started by a match, but by the catastrophic crash of two incredibly dense objects: neutron stars (the ultra-compressed corpses of dead stars) or a neutron star and a black hole. When these titans collide, they don't just smash together; they fling out a cloud of super-hot, radioactive debris. As this debris cools, it glows with a brilliant, short-lived light called a "kilonova." Think of a kilonova as a cosmic firework that lasts only a few days, lighting up the darkness with the glow of newly created heavy metals.

But these crashes are often accompanied by a second, even more energetic event: a gamma-ray burst. This is like a cosmic shotgun blast, shooting a narrow, high-speed beam of energy into space. Sometimes, if we are lucky enough to be standing in the right spot (or if the beam is wide enough), we see both the firework (the kilonova) and the blast (the gamma-ray burst). For a long time, scientists could only study these events if they had a "gravitational wave" detector—a kind of cosmic ear—to tell them exactly when and where a crash happened. But what if we don't have that ear? Can we still figure out what kind of stars crashed, how much stuff they threw out, and what happened next, just by looking at the light? This is the big question that drives the research in the paper you are about to read.

The Cosmic Detective Work: Solving the Mystery of Star Crashes Without a Sound

In this paper, a team of astronomers acts like cosmic detectives, trying to solve the mystery of what happens when stars crash, even when they can't "hear" the crash with gravitational wave detectors. They focused on a specific group of six gamma-ray bursts (GRBs) that happened in the relatively nearby universe. These bursts are the "gunshots" of the cosmos, and the team suspected that each one was followed by a "kilonova" firework. The challenge? The firework's light is often mixed up with the lingering glow of the gunshot (called the "afterglow"). It's like trying to see a specific candle flame in a room that is also lit by a flickering, noisy neon sign.

To solve this, the researchers used a powerful new tool called a "Bayesian analysis" within a framework named NMMA. Imagine this as a super-smart computer simulation that tries to build a 3D model of the crash. Instead of guessing the afterglow first and then subtracting it to find the kilonova (a method that often leaves errors), this tool builds both the afterglow and the kilonova at the same time. It tests millions of different scenarios: "What if the stars were two neutron stars? What if one was a black hole? What if the explosion was tilted this way or that way?" The computer then picks the scenario that fits the observed light best.

What They Found

The team successfully modeled the light from the gamma-ray bursts in their sample, but with one important caveat: while they found evidence of a kilonova in most cases, they were unable to confirm or exclude its presence in GRB 150101B. For that specific event, the data was too ambiguous to say for sure if the light was a kilonova or just the afterglow.

Here is what their detective work revealed about the nature of these crashes:

The Cast of Characters: For four of the events (GRB 160821B, GRB 170817A, GRB 211211A, and GRB 230307A), the evidence strongly points to a crash between two neutron stars. For the other two (GRB 150101B and GRB 191019A), the data slightly favors a crash between a neutron star and a black hole, but a two-neutron-star crash is still a possible explanation.
The Debris Cloud: When these stars crash, they throw out two types of debris. One is a fast, chaotic "dynamical" splash, and the other is a slower, steady "wind" of material. The researchers found that the wind mass is actually larger than the dynamical splash. On average, the wind mass was about 0.027 M⊙ (solar masses), while the dynamical mass was about 0.012 M⊙. This confirms what computer simulations have suggested: the slower wind carries more of the heavy stuff.
The Energy Connection: The team discovered a fascinating link between the power of the jet (the gamma-ray burst) and the amount of wind material thrown out. They found a mathematical relationship: log(Mwind) = −20.23 + 0.38 log(E0,J). In plain English, this means that the more powerful the jet is, the more wind material is likely to be ejected. It's like a bigger explosion throwing out a bigger cloud of smoke.
The Shape of the Crash: By analyzing the light, they could also tell how the stars were spinning and how heavy they were. They confirmed that as the "chirp mass" (a specific measure of the binary system's weight) goes down, the "tidal deformability" (how squishy the stars are) goes up. This matches what we know about how neutron stars behave under extreme pressure.

What They Ruled Out

The paper explicitly argues against the old way of doing things. Previously, scientists would try to model the afterglow (the jet's glow) first, subtract it from the total light, and call the leftover "pure" kilonova. The authors show that this method is flawed because it doesn't account for how the two lights mix and interfere with each other. Their simultaneous modeling suggests that the old method could lead to incorrect estimates of how much mass was ejected.

How Sure Are They?

The authors are very confident in their methodology because they tested it on the most famous event of all: GRB 170817A/AT2017gfo. This was the first time we saw a kilonova with both gravitational waves and light. When they ran their new model on this event, it perfectly matched the known results, proving their tool works. For the other events, where they don't have gravitational wave data, they are confident in their conclusions based on the strength of the light data and the robustness of their statistical model, though they note that for GRB 150101B, the presence of a kilonova remains a "claim" rather than a confirmed fact.

In the end, this paper shows that even without the "sound" of gravitational waves, we can use the "light" of gamma-ray bursts and kilonovae to figure out exactly what kind of stars crashed, how much gold and heavy metal they created, and how the universe builds its most precious materials. It turns out that looking closely at the light is enough to tell a very detailed story of a cosmic collision.

Technical Summary: Kilonova and Progenitor Properties of Merger-Driven Gamma-Ray Bursts

Problem and Motivation
Gamma-ray bursts (GRBs) and kilonovae (KNs) are the expected electromagnetic (EM) counterparts of binary neutron star (BNS) and neutron star-black hole (NSBH) mergers. While the gravitational wave (GW) detection of GW170817 provided a confirmed association with EM counterparts, it remains the only such event detected via GWs to date. Without GW signals, inferring progenitor properties from EM data alone is challenging because KN emission (thermal) and afterglow emission (nonthermal) are often temporally intertwined, particularly in on-axis or near-axis events. Traditional methods often isolate the KN by extrapolating the afterglow flux and subtracting it from optical/NIR data, a procedure that introduces uncertainties by failing to model the combined thermal and nonthermal mechanisms simultaneously. This work aims to overcome these limitations by performing a simultaneous Bayesian inference of afterglow and KN emission to constrain ejecta parameters and infer binary progenitor properties for a sample of merger-driven GRBs.

Methodology
The authors analyzed a sample of six merger-driven GRBs (GRB 150101B, 160821B, 170817A, 191019A, 211211A, and 230307A) that met strict selection criteria: robust redshift measurements, host galaxy associations ( $p_{cc} \lesssim 10^{-3}$ ), redshift $z \lesssim 0.5$ , sufficient optical/NIR spectral coverage, and well-sampled X-ray data without flaring activity.

The analysis utilized the Nuclear physics and Multi-Messenger Astrophysics (NMMA) framework (version 0.2.0), which employs Bayesian inference to simultaneously model the afterglow and KN.

Afterglow Modeling: The nonthermal emission was modeled using the afterglowpy library, considering synchrotron radiation from forward shocks. The study tested both Gaussian (GS) and top-hat (TH) jet structures, varying parameters such as isotropic equivalent kinetic energy ( $E_0$ ), viewing angle ( $\iota$ ), circumburst density ( $n_0$ ), and microphysical parameters ( $\epsilon_e, \epsilon_B, p$ ).
Kilonova Modeling: The thermal KN emission was modeled using the radiative transfer code POSSIS. The ejecta were decomposed into two components: a dynamical component ( $M_{dyn}$ ) emitted at high velocities during merger, and a secular wind component ( $M_{wind}$ ) emitted from the post-merger accretion disk. The models accounted for lanthanide-rich and lanthanide-poor compositions, with specific geometric constraints (half-opening angle $\Phi$ ) differing between BNS and NSBH scenarios.
Inference: A Bayesian analysis was performed using PyMultiNest to derive posterior distributions for all parameters. Model selection was conducted via Bayes factors to determine the most plausible progenitor scenario (BNS vs. NSBH) and jet structure. In the absence of GW data, the inferred KN ejecta properties were mapped to binary progenitor properties (chirp mass $M_{Chirp}$ , mass ratio $q$ , mass-weighted tidal deformability $\tilde{\Lambda}$ ) using phenomenological relations derived from numerical simulations.

Key Results

Progenitor Identification: The analysis favored a BNS progenitor for GRB 160821B, GRB 170817A, GRB 211211A, and GRB 230307A. For GRB 150101B and GRB 191019A, the NSBH scenario was slightly preferred based on Bayes factors, though a BNS origin remained viable. GRB 150101B was the only event where the presence of a KN could not be conclusively confirmed or excluded.
Ejecta Properties: The study found that the median wind mass ( $\langle M_{wind} \rangle = 0.027^{+0.046}_{-0.019} M_\odot$ ) is higher than the median dynamical mass ( $\langle M_{dyn} \rangle = 0.012^{+0.007}_{-0.006} M_\odot$ ). The total ejected mass median was $\langle M_{Total} \rangle = 0.039^{+0.043}_{-0.033} M_\odot$ .
Correlations: A statistically significant positive correlation was identified between the collimation-corrected jet kinetic energy ( $E_{0,J}$ ) and the wind ejecta mass, expressed as $\log(M_{wind}) = -20.23 + 0.38 \log(E_{0,J})$ .
Binary Parameters: The inferred binary properties showed distinct trends: BNS mergers exhibited higher tidal deformability ( $\tilde{\Lambda} \approx 304$ ) and lower chirp masses ( $M_{Chirp} \approx 1.23 M_\odot$ ) compared to NSBH mergers ( $\tilde{\Lambda} \approx 20$ , $M_{Chirp} \approx 1.3 M_\odot$ ). The results for GRB 170817A were consistent with GW measurements, validating the EM-only methodology.
Mass Ratio and Ejecta: A trend was observed where lower mass ratios ( $q$ ) correlated with higher dynamical ejecta masses, consistent with numerical simulations suggesting stronger tidal disruption in unequal-mass binaries.

Significance and Claims
The authors claim that this work presents, for the first time, the binary progenitor properties of a sizable sample of merger-driven GRBs using only EM data. The primary significance lies in demonstrating that EM modeling is a feasible pathway to investigate binary mergers in the absence of GW observations. By simultaneously modeling afterglow and KN emission, the study avoids the systematic uncertainties associated with subtracting afterglow components. The results reinforce the connection between KN ejecta properties and binary progenitor characteristics, confirming that numerical simulation predictions (e.g., the relationship between $M_{Chirp}$ and $\tilde{\Lambda}$ ) hold true when derived from EM observations. The identified correlation between wind mass and jet energy offers new insights into the physics of accretion-to-jet energy conversion. The authors conclude that their methodology successfully reproduces known properties of GW170817 and provides a robust framework for characterizing future merger events before GW detections are available.

The Cosmic Detective Work: Solving the Mystery of Star Crashes Without a Sound

More like this