Uniform Modeling of Observed Kilonovae: Implications for Diversity and the Progenitors of Merger-Driven Long Gamma-Ray Bursts

TL;DR

This study models eight kilonovae associated with short and long gamma-ray bursts to understand their diversity, ejecta properties, and implications for progenitors, revealing differences linked to GRB types and observational prospects.

Contribution

It introduces a uniform three-component kilonova model applied to multiple events, highlighting diversity and potential progenitor differences between short and long GRB-associated kilonovae.

Findings

Long GRB kilonovae have higher red ejecta masses than short GRB kilonovae.

A weak correlation exists between ejecta masses and gamma-ray energies.

The study forecasts observational capabilities of upcoming surveys for future kilonovae.

Abstract

We present uniform modeling of eight kilonovae, five following short gamma-ray bursts (GRBs; including GRB170817A) and three following long GRBs. We model their broadband afterglows to determine the relative contributions of afterglow and kilonova emission. We fit the kilonovae using a three-component model in MOSFiT that accounts for ejecta geometry, and find population median ejecta masses for the total, blue ($\kappa_{B} = 0.5$ cm^2 / g), purple ($\kappa_{P} = 3$ cm^2 / g), and red ($\kappa_{R} = 10$ cm^2 / g) components of $M_{ej, tot} = 0.085_{-0.040}^{+0.110} M_{\odot}$, $M_{ej, B} = 0.006_{-0.004}^{+0.015} M_{\odot}$, $M_{ej, P} = 0.020_{-0.010}^{+0.034} M_{\odot}$, and $M_{ej, R} = 0.051_{-0.045}^{+0.100} M_{\odot}$ (68% confidence). The kilonova of GW170817 is near the median of the sample in most derived properties, while the sample indicates great diversity. We investigate trends between the ejecta masses and the isotropic-equivalent and beaming-corrected gamma-ray energies ($E_{\gamma, iso}$, $E_{\gamma}$), as well as rest-frame durations ($T_{90, rest}$). We find long GRB kilonovae have higher median red ejecta masses ($M_{ej, R} > 0.05 M_{\odot}$) compared to on-axis short GRB kilonovae ($M_{ej, R} < 0.02 M_{\odot}$). We also observe a weak scaling between the total and red ejecta masses with $E_{\gamma, iso}$ and $E_{\gamma}$, though a larger sample is needed to establish a significant correlation. These findings imply a connection between merger-driven long GRBs and larger tidal dynamical ejecta masses, which may indicate that their progenitors are asymmetric compact object binaries. We produce representative kilonova light curves and find that the planned depths and cadences of the Rubin and Roman Observatory surveys will be sufficient for order-of-magnitude constraints on $M_{ej, B}$ (and, for Roman, $M_{ej, P}$ and $M_{ej, R}$) of future kilonovae at $z < 0.1$.