Detectability of kilonovae in optical surveys: post-mortem examination of the LVC O3 run follow-up

The detection of the binary neutron star (BNS) merger GW170817 and the associated electromagnetic (EM) counterpart, the 'kilonova' (kN) AT2017gfo, opened a new era in multimessenger astronomy. However, despite many efforts, it has been proven very difficult to find additional kNe, even though LIGO/Virgo has reported at least one BNS event during their latest run, O3. The focus of this work is the exploration of the sensitivity of the adopted optical surveys searching for kNe during O3. We propose ways to optimize the choices of filters and survey depth to boost the detection efficiency for these faint and fast-evolving transients in the future. In particular, we use kN models to explore the dependence on ejecta mass, geometry, viewing angle, wavelength coverage, and source distance. We find that the kN detection efficiency has a strong viewing-angle dependence, especially for filters blueward of i-band. This loss of sensitivity can be mitigated by early, deep, observations. Efficient gri counterpart searches for kNe at ∼200 Mpc would require reaching a limiting magnitude mlim = 23 mag, to ensure good sensitivity over a wide range of the model phase-space. We conclude that kN searches during O3 were generally too shallow to detect BNS optical counterparts, even under optimistic assumptions

Constraining the Observer Angle of the Kilonova AT2017gfo Associated with GW170817: Implications for the Hubble Constant

There is a strong degeneracy between the luminosity distance (D L) and the observer viewing angle (θ obs; hereafter viewing angle) of the gravitational wave (GW) source with an electromagnetic counterpart, GW170817. Here, for the first time, we present independent constraints on IMG ALIGN="MIDDLE" ALT="${θ}_{\mathrm{obs}}={32.5}-{-9.7}^{^\circ +11.7}$" SRC="apjab5799ieqn1.gif"/from broadband photometry of the kilonova (kN) AT2017gfo associated with GW170817. These constraints are consistent with independent results presented in the literature using the associated gamma-ray burst GRB 170817A. Combining the constraints on θ obs with the GW data, we find an improvement of 24% on H 0. The observer angle constraints are insensitive to other model parameters, e.g., the ejecta mass, the half-opening angle of the lanthanide-rich region and the temperature. A broad wavelength coverage extending to the near-infrared is helpful to robustly constrain θ obs. While the improvement on H 0 presented here is smaller than the one from high angular resolution imaging of the radio counterpart of GW170817, kN observations are significantly more feasible at the typical distances of such events from current and future LIGO-Virgo collaboration observing runs (D L ∼ 100 Mpc). Our results are insensitive to the assumption of the peculiar velocity of the kN host galaxy. © 2020. The American Astronomical Society

Optimizing Cadences with Realistic Light Curve Filtering for Serendipitous Kilonova Discovery with Vera Rubin Observatory

Current and future optical and near-infrared wide-field surveys have the potential of finding kilonovae, the optical and infrared counterparts to neutron star mergers, independently of gravitational-wave or high-energy gamma-ray burst triggers. The ability to discover fast and faint transients such as kilonovae largely depends on the area observed, the depth of those observations, the number of re-visits per field in a given time frame, and the filters adopted by the survey; it also depends on the ability to perform rapid follow-up observations to confirm the nature of the transients. In this work, we assess kilonova detectability in existing simulations of the LSST strategy for the Vera C. Rubin Wide Fast Deep survey, with focus on comparing rolling to baseline cadences. Although currently available cadences can enable the detection of more than 300 kilonovae out to 1400 Mpc over the ten-year survey, we can expect only 3-32 kilonovae similar to GW170817 to be recognizable as fast-evolving transients. We also explore the detectability of kilonovae over the plausible parameter space, focusing on viewing angle and ejecta masses. We find that observations in redder izy bands are crucial for identification of nearby (within 300 Mpc) kilonovae that could be spectroscopically classified more easily than more distant sources. Rubin's potential for serendipitous kilonova discovery could be increased by gain of efficiency with the employment of individual 30s exposures (as opposed to 2x15s snap pairs), with the addition of red-band observations coupled with same-night observations in g- or r-bands, and possibly with further development of a new rolling-cadence strategy

GRB 171010A/SN 2017htp: a GRB-SN at z = 0.33

The number of supernovae known to be connected with long-duration gamma-ray bursts (GRBs) is increasing and the link between these events is no longer exclusively found at low redshift (z ≲ 0.3) but is well established also at larger distances. We present a new case of such a liaison at z = 0.33 between GRB 171010A and SN 2017htp. It is the second closest GRB with an associated supernova of only three events detected by Fermi-LAT. The supernova is one of the few higher redshift cases where spectroscopic observations were possible and shows spectral similarities with the well-studied SN 1998bw, having produced a similar Ni mass (⁠ M Ni =0.33±0.02 M ⊙ MNi=0.33±0.02 M⊙ ⁠) with slightly lower ejected mass (⁠ M ej =4.1±0.7 M ⊙ Mej=4.1±0.7 M⊙ ⁠) and kinetic energy (⁠ E K =8.1±2.5× 10 51 erg EK=8.1±2.5×1051 erg ⁠). The host-galaxy is bigger in size than typical GRB host galaxies, but the analysis of the region hosting the GRB revealed spectral properties typically observed in GRB hosts and showed that the progenitor of this event was located in a very bright H ii region of its face-on host galaxy, at a projected distance of ∼ 10 kpc from its galactic centre. The star-formation rate (SFRGRB ∼ 0.2 M⊙ yr−1) and metallicity (12 + log(O/H) ∼8.15 ± 0.10) of the GRB star-forming region are consistent with those of the host galaxies of previously studied GRB–SN systems

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 Cadences with Realistic Light-curve Filtering for Serendipitous Kilonova Discovery with Vera Rubin Observatory

Current and future optical and near-infrared wide-field surveys have the potential to find kilonovae, the optical and infrared counterparts to neutron star mergers, independently of gravitational-wave or high-energy gamma-ray burst triggers. The ability to discover fast and faint transients such as kilonovae largely depends on the area observed, the depth of those observations, the number of revisits per field in a given time frame, and the filters adopted by the survey; it also depends on the ability to perform rapid follow-up observations to confirm the nature of the transients. In this work, we assess kilonova detectability in existing simulations of the Legacy Survey of Space and Time strategy for the Vera C. Rubin Wide Fast Deep survey, with focus on comparing rolling to baseline cadences. Although currently available cadences can enable the detection of &gt;300 kilonovae out to ∼1400 Mpc over the 10 year survey, we can expect only 3-32 kilonovae similar to GW170817 to be recognizable as fast-evolving transients. We also explore the detectability of kilonovae over the plausible parameter space, focusing on viewing angle and ejecta masses. We find that observations in redder izy bands are crucial for identification of nearby (within 300 Mpc) kilonovae that could be spectroscopically classified more easily than more distant sources. Rubin's potential for serendipitous kilonova discovery could be increased by gain of efficiency with the employment of individual 30 s exposures (as opposed to 2 × 15 s snap pairs), with the addition of red-band observations coupled with same-night observations in g or r bands, and possibly with further development of a new rolling-cadence strategy

Hubble constant and nuclear equation of state from kilonova spectro-photometric light curves

The merger of two compact objects of which at least one is a neutron star is signalled by transient electromagnetic emission in a kilonova (KN). This event is accompanied by gravitational waves and possibly other radiation messengers such as neutrinos or cosmic rays. The electromagnetic emission arises from the radioactive decay of heavy r-process elements synthesized in the material ejected during and after the merger. In this paper we show that the analysis of KNe light curves can provide cosmological distance measurements and constrain the properties of the ejecta. In this respect, MAAT, the new Integral Field Unit in the OSIRIS spectrograph on the 10.4 m Gran Telescopio CANARIAS (GTC), is well suited for the study of KNe by performing absolute spectro-photometry over the entire 3600-10 000 spectral range. Here, we study the most representative cases regarding the scientific interest of KNe from binary neutron stars, and we evaluate the observational prospects and performance of MAAT on the GTC to do the following: (a) study the impact of the equation of state on the KN light curve, and determine to what extent bounds on neutron star (NS) radii or compactness deriving from KN peak magnitudes can be identified and (b) measure the Hubble constant, H0, with precision improved by up to 40%, when both gravitational wave data and photometric-light curves are used. In this context we discuss how the equation of state, the viewing angle, and the distance affect the precision and estimated value of H0