Evidence for two distinct populations of kilonova-associated Gamma Ray Bursts

Identification of Gamma Ray Burst (GRB) progenitors based on the duration of their prompt emission ($T_{90}$) has faced several roadblocks recently. Long-duration GRBs (with $T_{90} > 2s$) have traditionally been thought to be originating from the collapse of massive stars, and the short-duration ones (with $T_{90} < 2s$) from compact binary mergers. However, recent observations of a long GRB associated with a kilonova (KN) and a short GRB with supernova (SN) association demand a more detailed classification of the GRB population. In this {\it Letter}, we focus on GRBs associated with KNe, believed to be originating from mergers of binaries involving neutron stars (NS). We make use of the GRB prompt emission light curves of {\it Swift}-BAT 2022 GRB catalog and employ machine learning algorithms to study the classification of GRB progenitors. Our analysis reveals that there are five distinct clusters of GRBs, of which the KN-associated GRBs are located in two separate clusters indicating they may have been produced by different progenitors. We argue that these clusters may be due to subclasses of binary neutron star (BNS) and/or neutron star--black hole (NS-BH) mergers. We also discuss the implications of these findings for future gravitational-wave (GW) observations and how those observations may help in understanding these clusters better.