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

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$.