The landmark multi-messenger observations of the binary neutron star (BNS) merger GW170817 provided firm evidence that such mergers can produce short gamma-ray bursts (sGRBs). However, the limited number of BNS detections by current gravitational wave observatories raises the question of whether BNS mergers alone can account for the full observed sGRB population. We analysed a comprehensive set of 64 BNS population synthesis models with a Monte Carlo–based framework to reproduce the properties of sGRBs detected by GBM over the past 16 years. We considered three jet geometry scenarios: a universal structured jet calibrated to GW170817, a universal top-hat jet, and a non-universal top-hat jet with distributions of core opening angles. Our results show that models characterised by low local BNS merger rates (R_ Fermi- BNS (0) łesssim 50, Gpc ^ -3, yr ^ -1) predict too few observable sGRBs to reproduce the GBM population, effectively disfavouring them as sole progenitors. Even when relaxing the assumptions on jet geometry, low-rate models remain viable only for wide jets (þetac ≥ 15^̧irc), in tension with the narrow jet cores (þetac ≈ 6^̧irc) inferred from sGRB afterglow observations. In contrast, models with local merger rates of order R_ Fermi- BNS (0) ≈ 100, Gpc ^ -3, yr ^ -1 successfully reproduce the observed sGRB population, assuming a plausible fraction of BNS mergers launch relativistic jets and realistic jet geometries. This analysis highlights the power of combining binary neutron star population synthesis models with sGRB observations and comparing them with current gravitational wave constraints to assess the role of BNS mergers as sGRB progenitors and to identify any emerging tensions between these populations.
Santis et al. (Fri,) studied this question.