Strategy for identifying Vera C. Rubin Observatory kilonova candidates for targeted gravitational-wave searches

Abstract: Since the observation of the binary neutron star merger GW170817 and the associated kilonova AT2017gfo, the next joint gravitational-wave/optical kilonova has been highly anticipated. Overlapping observations between the Vera C. Rubin Observatory and the international gravitational-wave detector network are expected soon. Wide-field survey facilities, such as Rubin, can serve dual roles in gravitational-wave astronomy: conducting dedicated optical counterpart searches following gravitational-wave triggers and, through surveys such as the Legacy Survey of Space and Time (LSST), providing opportunities for fortuitous kilonova discoveries during routine operations. We use simulations to develop a strategy for identifying kilonova candidates observed by Rubin and processed by the Fink broker. These candidates can be used as astrophysical triggers for a targeted gravitational-wave search. We simulate kilonovae light-curves for the first year of Rubin with the latest observing strategy for the Wide-Fast-Deep and the Deep Drilling Fields. Assuming a kilonova rate of 250 Gpc$^{-3}$ yr$^{-1}$, we find that Rubin brokers should observe $\sim 4$ kilonovae per year with at least one alert above a signal-to-noise ratio of 5 within the gravitational-wave detector horizon ($\sim 350$ Mpc). Most of these will be faint, and detected 1-2 days following the neutron star merger. Photometric and spectroscopic follow-up will be limited to large telescopes. Using archival data from the Zwicky Transient Facility (ZTF) and our proposed selection criteria, we estimate a minimum contamination of at least 30 events per month from other transients and variables, even under our strictest selection criteria. A deep gravitational-wave search targeting Rubin kilonova candidates may lead to the next multi-messenger discovery.