Fast-transient Searches in Real Time with ZTFReST: Identification of Three Optically Discovered Gamma-Ray Burst Afterglows and New Constraints on the Kilonova Rate

The most common way to discover extragalactic fast transients, which fade within a few nights in the optical, is via follow-up of gamma-ray burst and gravitational-wave triggers. However, wide-field surveys have the potential to identify rapidly fading transients independently of such external triggers. The volumetric survey speed of the Zwicky Transient Facility (ZTF) makes it sensitive to objects as faint and fast fading as kilonovae, the optical counterparts to binary neutron star mergers, out to almost 200 Mpc. We introduce an open-source software infrastructure, the ZTF REaltime Search and Triggering, ZTFReST, designed to identify kilonovae and fast transients in ZTF data. Using the ZTF alert stream combined with forced point-spread-function photometry, we have implemented automated candidate ranking based on their photometric evolution and fitting to kilonova models. Automated triggering, with a human in the loop for monitoring, of follow-up systems has also been implemented. In 13 months of science validation, we found several extragalactic fast transients independently of any external trigger, including two supernovae with post-shock cooling emission, two known afterglows with an associated gamma-ray burst (ZTF20abbiixp, ZTF20abwysqy), two known afterglows without any known gamma-ray counterpart (ZTF20aajnksq, ZTF21aaeyldq), and three new fast-declining sources (ZTF20abtxwfx, ZTF20acozryr, ZTF21aagwbjr) that are likely associated with GRB200817A, GRB201103B, and GRB210204A. However, we have not found any objects that appear to be kilonovae. We constrain the rate of GW170817-like kilonovae to R < 900 Gpc-3 yr-1 (95% confidence). A framework such as ZTFReST could become a prime tool for kilonova and fast-transient discovery with the Vera Rubin Observatory

Constraining the Kilonova Rate with Zwicky Transient Facility Searches Independent of Gravitational Wave and Short Gamma-Ray Burst Triggers

The first binary neutron star merger, GW170817, was accompanied by a radioactivity-powered optical/infrared transient called a kilonova. To date, no compelling kilonova has been found in all-sky optical surveys, independently of short gamma-ray burst and gravitational-wave triggers. In this work, we searched the first 23 months of the Zwicky Transient Facility (ZTF) data stream for candidate kilonovae in the form of rapidly evolving transients. We combined ZTF alert queries with forced point-spread-function photometry and nightly flux stacking to increase our sensitivity to faint and fast transients. Automatic queries yielded >11,200 candidates, 24 of which passed quality checks and selection criteria based on a grid of kilonova models tailored for both binary neutron star and neutron star-black hole mergers. None of the candidates in our sample was deemed a possible kilonova after thorough vetting. The sources that passed our selection criteria are dominated by Galactic cataclysmic variables. We identified two fast transients at high Galactic latitude, one of which is the confirmed afterglow of long-duration GRB 190106A, the other is a possible cosmological afterglow. Using a survey simulation code, we constrained the kilonova rate for a range of models including top-hat, linearly decaying light curves, and synthetic light curves obtained with radiative transfer simulations. For prototypical GW170817-like kilonovae, we constrain the rate to be Gpc-3 yr-1 (95% confidence). By assuming a population of kilonovae with the same geometry and composition of GW170817 observed under a uniform viewing angle distribution, we obtained a constraint on the rate of R < 4029 Gpc-3 yr-1. © 2020. The American Astronomical Society

Constraining the Kilonova Rate with Zwicky Transient Facility Searches Independent of Gravitational Wave and Short Gamma-ray Burst Triggers

The first binary neutron star merger, GW170817, was accompanied by a radioactivity-powered optical/infrared transient called a kilonova. To date, no compelling kilonova has been found during optical surveys of the sky, independent of gravitational-wave triggers. In this work, we searched the first 23 months of the Zwicky Transient Facility (ZTF) data stream for candidate kilonovae in the form of rapidly evolving transients. We combined ZTF alert queries with forced point-spread-function photometry and nightly flux stacking to increase our sensitivity to faint and fast transients. Automatic queries yielded $>11,200$ candidates, 24 of which passed quality checks and strict selection criteria based on a grid of kilonova models tailored for both binary neutron star and neutron star-black hole mergers. None of the candidates in our sample was deemed a possible kilonova after thorough vetting, catalog cross-matching, and study of their color evolution. The sources that passed our selection criteria are dominated by Galactic cataclysmic variables. In addition, we identified two fast transients at high Galactic latitude, one of which is the confirmed afterglow of long-duration GRB190106A, and the other is a possible cosmological afterglow. Using a survey simulation code, we constrained the kilonova rate for a range of models including top-hat and linearly decaying light curves and synthetic light curves obtained with radiative transfer simulations. For prototypical GW170817-like kilonovae, we constrain the rate to be $R < 1775$ Gpc$^{-3}$ yr$^{-1}$ at 95% confidence level by requiring at least 2 high-significance detections. By assuming a population of kilonovae with the same geometry and composition of GW170817 observed under a uniform viewing angle distribution, we obtained a constraint on the rate of $R < 4029$ Gpc$^{-3}$ yr$^{-1}$.Comment: Submitted for publication in Ap

Optical follow-up of the neutron star-black hole mergers S200105ae and S200115j

LIGO and Virgo’s third observing run revealed the first neutron star–black hole (NSBH) merger candidates in gravitational waves. These events are predicted to synthesize r-process elements1,2 creating optical/near-infrared ‘kilonova’ emission. The joint gravitational wave and electromagnetic detection of an NSBH merger could be used to constrain the equation of state of dense nuclear matter3, and independently measure the local expansion rate of the Universe4. Here, we present the optical follow-up and analysis of two of the only three high-significance NSBH merger candidates detected to date, S200105ae and S200115j, with the Zwicky Transient Facility5. The Zwicky Transient Facility observed ~48% of S200105ae and ~22% of S200115j’s localization probabilities, with observations sensitive to kilonovae brighter than −17.5 mag fading at 0.5 mag d−1 in the g- and r-bands; extensive searches and systematic follow-up of candidates did not yield a viable counterpart. We present state-of-the-art kilonova models tailored to NSBH systems that place constraints on the ejecta properties of these NSBH mergers. We show that with observed depths of apparent magnitude ~22 mag, attainable in metre-class, wide-field-of-view survey instruments, strong constraints on ejecta mass are possible, with the potential to rule out low mass ratios, high black hole spins and large neutron star radii

Optimizing serendipitous detections of kilonovae: cadence and filter selection

The rise of multimessenger astronomy has brought with it the need to exploit all available data streams and learn more about the astrophysical objects that fall within its breadth. One possible avenue is the search for serendipitous optical/near-infrared counterparts of gamma-ray bursts (GRBs) and gravitational-wave (GW) signals, known as kilonovae. With surveys such as the Zwicky Transient Facility (ZTF), which observes the sky with a cadence of ∼3 d, the existing counterpart locations are likely to be observed; however, due to the significant amount of sky to explore, it is difficult to search for these fast-evolving candidates. Thus, it is beneficial to optimize the survey cadence for realtime kilonova identification and enable further photometric and spectroscopic observations. We explore how the cadence of wide field-of-view surveys like ZTF can be improved to facilitate such identifications. We show that with improved observational choices, e.g. the adoption of three epochs per night on a ∼nightly basis, and the prioritization of redder photometric bands, detection efficiencies improve by about a factor of two relative to the nominal cadence. We also provide realistic hypothetical constraints on the kilonova rate as a form of comparison between strategies, assuming that no kilonovae are detected throughout the long-term execution of the respective observing plan. These results demonstrate how an optimal use of ZTF increases the likelihood of kilonova discovery independent of GWs or GRBs, thereby allowing for a sensitive search with less interruption of its nominal cadence through Target of Opportunity programs

Optical follow-up of the neutron star–black hole mergers S200105ae and S200115j

International audienceThe optical follow-up and analysis of two neutron star–black hole merger candidates with the Zwicky Transient Facility did not yield viable counterparts. However, state-of-the-art kilonova models constrain the ejecta properties of these mergers

Searching for gravitational wave optical counterparts with the Zwicky Transient Facility: summary of O4a

International audienceDuring the first half of the fourth observing run (O4a) of the International Gravitational Wave Network (IGWN), the Zwicky Transient Facility (ZTF) conducted a systematic search for kilonova (KN) counterparts to binary neutron star (BNS) and neutron star-black hole (NSBH) merger candidates. Here, we present a comprehensive study of the five high-significance (FAR < 1 per year) BNS and NSBH candidates in O4a. Our follow-up campaigns relied on both target-of-opportunity observations (ToO) and re-weighting of the nominal survey schedule to maximize coverage. We describe the toolkit we have been developing, Fritz, an instance of SkyPortal, instrumental in coordinating and managing our telescope scheduling, candidate vetting, and follow-up observations through a user-friendly interface. ZTF covered a total of 2841 deg$^2$ within the skymaps of the high-significance GW events, reaching a median depth of g~20.2 mag. We circulated 15 candidates, but found no viable KN counterpart to any of the GW events. Based on the ZTF non-detections of the high-significance events in O4a, we used a Bayesian approach, nimbus, to quantify the posterior probability of KN model parameters that are consistent with our non-detections. Our analysis favors KNe with initial absolute magnitude fainter than -16 mag. The joint posterior probability of a GW170817-like KN associated with all our O4a follow-ups was 64%. Additionally, we use a survey simulation software, simsurvey, to determine that our combined filtered efficiency to detect a GW170817-like KN is 36%, when considering the 5 confirmed astrophysical events in O3 (1 BNS and 4 NSBH), along with our O4a follow-ups. Following Kasliwal et al. (2020), we derived joint constraints on the underlying KN luminosity function based on our O3 and O4a follow-ups, determining that no more than 76% of KNe fading at 1 mag/day can peak at a magnitude brighter than -17.5 mag

Correction: Discovery and confirmation of the shortest gamma-ray burst from a collapsar (Jul, 10.1038/s41550-021-01428-7, 2021)

Correction to: Nature Astronomy https://doi.org/10.1038/s41550-021-01428-7, published online 26 July 202

Kilonova Luminosity Function Constraints Based on Zwicky Transient Facility Searches for 13 Neutron Star Merger Triggers during O3

We present a systematic search for optical counterparts to 13 gravitational wave (GW) triggers involving at least one neutron star during LIGO/Virgo's third observing run (O3). We searched binary neutron star (BNS) and neutron star black hole (NSBH) merger localizations with the Zwicky Transient Facility (ZTF) and undertook follow-up with the Global Relay of Observatories Watching Transients Happen (GROWTH) collaboration. The GW triggers had a median localization area of 4480 deg2, a median distance of 267 Mpc, and false-alarm rates ranging from 1.5 to 10-25 yr-1. The ZTF coverage in the g and r bands had a median enclosed probability of 39%, median depth of 20.8 mag, and median time lag between merger and the start of observations of 1.5 hr. The O3 follow-up by the GROWTH team comprised 340 UltraViolet/Optical/InfraRed (UVOIR) photometric points, 64 OIR spectra, and three radio images using 17 different telescopes. We find no promising kilonovae (radioactivity-powered counterparts), and we show how to convert the upper limits to constrain the underlying kilonova luminosity function. Initially, we assume that all GW triggers are bona fide astrophysical events regardless of false-alarm rate and that kilonovae accompanying BNS and NSBH mergers are drawn from a common population; later, we relax these assumptions. Assuming that all kilonovae are at least as luminous as the discovery magnitude of GW170817 (-16.1 mag), we calculate that our joint probability of detecting zero kilonovae is only 4.2%. If we assume that all kilonovae are brighter than-16.6 mag (the extrapolated peak magnitude of GW170817) and fade at a rate of 1 mag day-1 (similar to GW170817), the joint probability of zero detections is 7%. If we separate the NSBH and BNS populations based on the online classifications, the joint probability of zero detections, assuming all kilonovae are brighter than-16.6 mag, is 9.7% for NSBH and 7.9% for BNS mergers. Moreover, no more than &lt;57% (&lt;89%) of putative kilonovae could be brighter than-16.6 mag assuming flat evolution (fading by 1 mag day-1) at the 90% confidence level. If we further take into account the online terrestrial probability for each GW trigger, we find that no more than &lt;68% of putative kilonovae could be brighter than-16.6 mag. Comparing to model grids, we find that some kilonovae must have M ej &lt; 0.03 M o˙, X lan &gt; 10-4, or φ &gt; 30° to be consistent with our limits. We look forward to searches in the fourth GW observing run; even 17 neutron star mergers with only 50% coverage to a depth of-16 mag would constrain the maximum fraction of bright kilonovae to &lt;25%