---
abstract: 'The detection of the binary neutron star merger GW170817 and the associated electromagnetic (EM) counterpart, the “kilonova" (kN) AT2017gfo, opened a new era in multi-messenger astronomy. However, despite many efforts, it has been proven very difficult to find additional EM counterparts to the many LIGO/Virgo compact binary merger candidates since then. We explore the sensitivity of ongoing optical surveys searching for kNe and propose ways to optimize the choices of filters and survey depth to boost the detection efficiency for these faint and fast-evolving transients. In particular, we use kN models to explore the dependence on ejecta mass, geometry, viewing angle, wavelength coverage and source distance. Considering that the O3 LIGO/Virgo binary neutron star merger candidates were at distances about 5 times further than GW170817, we find that efficient $gri$ counterpart search would require reaching a limiting magnitude $m_{\rm lim}=$ 23 mag, given the model phase-space. We conclude that kN searches during O3 were generally too shallow to observe optical counterparts, even under optimistic assumptions, thus providing a plausible explanation for the lack of detections.'
author:
- |
    A. Sagués Carracedo,$^{1}$[^1] M. Bulla,$^{1,2}$ U. Feindt$^{1}$ and A. Goobar$^{1}$\
    $^{1}$ The Oskar Klein Centre for Cosmoparticle,Department of Physics, Stockholm University, AlbaNova, Stockholm SE106 91, Sweden\
    $^{2}$ Nordita, KTH Royal Institute of Technology and Stockholm University, Roslagstullsbacken 23, SE-106 91 Stockholm, Sweden\
bibliography:
- 'mnras.bib'
date: 'Accepted XXX. Received YYY; in original form ZZZ'
title: 'Detectability of kilonovae in optical surveys: [*post-mortem*]{} examination of the LVC O3 run follow-up'
---

\[firstpage\]

neutron star mergers – gravitational waves – surveys

Introduction
============

On August 17th, 2017, during the second run of Advanced LIGO [@AdvanceLigo2015] and Advanced Virgo [@2015CQGra..32b4001A] collaborations, the first gravitational-wave (GW) signal from the merger of two neutron stars (BNS) was detected, GW170817. Multiple independent detections associated with this event were made throughout the entire electromagnetic (EM) spectrum, including a short gamma-ray burst (@Abbott2017; @2017ApJ...848L..14G; @2017ApJ...848L..15S; @2017ApJ...848L..21A; @2017ApJ...848L..25H; @2017Sci...358.1579H; @2017ApJ...848L..20M; @2017Sci...358.1556C;@Troja17) and the kilonova (kN) AT2017gfo found at optical and IR wavelengths (@Cowperthwaite2017; @Kasliwal2017; @Tanvir2017; @SoaresSantos2017; @Andreoni17; @spitzer19). A record number of studies were carried out on this single event, and many open questions remain with regards to the rates of neutron mergers, the nucleosynthesis of heavy elements, the mass range and equation of state of neutron stars, the diversity in the kN explosion population and their use as BNS mergers as distance yardsticks in cosmology. At only $\sim 40$ Mpc distance, in the near-by galaxy NGC 4993, and close to face-on towards us [@Hotokezaka19], GW170817/AT2017gfo was extremely favourable for detection.

On April 1st, 2019, LIGO/Virgo started the third observing run (O3) on the search of GWs. Since then, astronomers work closely with the LIGO and Virgo collaboration (LVC), who provide search maps with their public alerts used by the astronomy community to carry out target of opportunity (ToO) follow-up programs. One year later, by March 27 2020, the O3 run was prematurely suspended due to the COVID-19 pandemic. In spite of 14 merger alerts been issued by the LVC including at least one neutron star, i.e., either BNS or neutron star - black hole (NS-BH), no kN was securely identified. In this work we explore possible explanations for the lack of optical counterparts, in spite of rapid responses by the optical astronomy community to the LVC alerts. We quantify the kN detection probability for the typical specifications of the follow-up surveys adopted, and suggest potential improvements. Furthermore, we investigate the feasibility of serendipitous kN detection, independently of LVC triggers.

Significant volumes of space are explored with high cadence thanks to the advent of optical sky surveys such as PanSTARRS [@panstarrs], the AllSky Automated Survey for Supernovae [ASASSN, @asasn], the Dark Energy Survey [DES, @des], the Asteroid Terrestrial impact Last Alert System [ATLAS, @atlas], the Zwicky Transient Facility [ZTF, @ztf], the Gravitational-Wave Optical Transient Observer [GOTO, @goto], the Mobile Astronomical System of TElescope Robots [MASTER, @master], the deca-degree optical transient imager [DDOTI, @ddoti], MeerLICHT [@meerlicht], KMT-Net [@kmtnet], the Rapid Action Telescope for Transient Objects [TAROT, @tarot] or the Global Rapid Advanced Network Devoted to the Multi-messenger Addicts [GRANDMA, @antier2020; @GRANDMA]. These surveys have also been carrying out regular targeted observations searching for potential EM counterparts of GW alerts in the regions posted by the LVC during O3, whenever the weather and the observable region of the sky made it possible.

In this work, we combine theoretical kN models with survey simulation tools to explore the detection probability. We focus on optical wavelengths, making use of the ZTF $gri$ filter system, at effective wavelengths 4814, 6422 and 7883 Å, as a representative set of photometric observations. Viewing-angle dependent models created with the radiative transfer code <span style="font-variant:small-caps;">possis</span> [@possis] were generated, and we choose a [*fiducial*]{} model that best fits the photometric data of AT2017gfo as our bench mark, also when comparing other models.

In Section \[sec:models\], we describe the properties of synthetic kN lightcurves from <span style="font-variant:small-caps;">possis</span>. In Section \[sec:surveysimulations\], we present the details of the observing strategies explored, along with the distribution of transients generated and the detection criterion considered. Section \[sec:LVCtriggers\] focuses on the detection probabilities for the LVC alerts in O3. In Section \[sec:synthetic\_population\], we describe the population of synthetic kNe selected under various scenarios. Section \[sec:rates\] considers attainable constraints on kN rates from serendipitous searches, unrelated to LVC alerts. Finally in Section \[sec:discussion\], we discuss and summarize the results.

Kilonova models {#sec:models}
===============

In this work, we consider kN models synthesized using the Monte Carlo radiative transfer code <span style="font-variant:small-caps;">possis</span> [@Bulla2015; @possis]. Assuming homologous expansion and analytic functions for the wavelength- and time-dependence of the opacity, <span style="font-variant:small-caps;">possis</span> performs time-dependent radiative transfers simulations for multi-dimensional models. Time-dependent spectral energy distributions (SEDs) are computed for different viewing angles and used to construct multi-band lightcurves.

Compared to the version of <span style="font-variant:small-caps;">possis</span> in [@possis], here we use an improved version, where the temperature is no longer parameterized and uniform throughout the ejecta, but estimated self-consistently from the mean intensity of the radiation field at each time and in each zone. Furthermore, we adopt thermalization efficiencies $\epsilon_\mathrm{th}$ from [@Barnes2016] rather than assuming $\epsilon_\mathrm{th}=0.5$ as in [@possis]. For more information about the code, we refer the reader to [@possis].

A two-dimensional geometry tailored to BNS mergers is used, comprising two ejecta components as in [@possis]: a lanthanide-rich component distributed within an angle $\pm\phi$ around the merger/equatorial plane, and a lanthanide-free component at higher latitudes (see fig. 1 in @possis). A density profile $\propto \varv^{-3}$, where $\varv$ is the velocity, is adopted for both components between $\varv_\mathrm{min}=0.025$c (instead of $\varv_\mathrm{min}=0$ as in @possis) and $\varv_\mathrm{max}=0.3$c (where c is the speed of light). The three free model parameters explored in this work are:

-   the total ejecta mass, $m_\mathrm{\rm ej}$

-   the half-opening angle, $\phi$

-   the viewing angle, $\Theta$

For $m_\mathrm{\rm ej}$, we consider $0.01$ to $0.10$ $M_{\odot}$ with a mass step of $0.01$ $M_{\odot}$, and compute half-opening angles, $\phi$, from 15 to 75$^\circ$ with steps of 15$^\circ$. A total number of 50 simulations are generated, each comprising predictions for viewing angles $\Theta$ from pole ($\cos\Theta=1$) to equator ($\cos\Theta=0$). SEDs used in this work are made available at .

![image](figures/lcs_models.pdf){width="\textwidth"}

Synthetic lightcurves {#sec:lcs}
---------------------

The ugrizyJH photometric data of AT2017gfo (@Andreoni17; @Arcavi17; @Chornock17; @Cowperthwaite2017; @Drout17; @Evans17; @Kasliwal2017; @Pian17; @Smartt17; @Tanvir2017; @Troja17; @Utsumi17; @Valenti17) can be reproduced by a [*fiducial*]{} model with $m_{\rm ej}=0.06\,M_{\odot}$ and $\phi=60^{\circ}$ for viewing angles close to polar, as shown in Fig. \[fig:lcs\_phot\] for the three filters of interest in this work, $gri$. The figure also shows that kNe are fast-evolving transients with absolute peak magnitudes around -16 mag, i.e., relatively faint compared to many types of supernovae.

The absolute peak magnitudes as a function of ejecta mass ranges from -12.6 to -17.4 mag, with a decay rate for 3 days after peak ($\Delta m_3$) varying from 0.4 to 2.2 mag, depending on the inclination angle and the broadband filter used, as shown in Fig. \[fig:peak\_dmag\]. For the fiducial model ($m_{\rm ej}=0.06\,\mathrm{M_{\sun}}$ and $\phi=60^{\circ}$), the optical peak magnitudes vary from -13.8 to , and $\Delta m_3$ from 0.6 to 1.5 mag for different viewing angles.

Finally, Fig. \[fig:mag\_z\] shows the observed r-band magnitude as a function of redshift, corresponding to the range of absolute magnitudes at peak, emphasizing the faintness of kNe, and thus the limited redshift range that can be probed with small telescopes. Brighter lightcurves result from larger $m_{\rm ej}$, smaller $\phi$, and $\Theta$ close to viewing angles towards the pole. The difference in brightness and shape of the synthetic kN lightcurves is less pronounced at redder wavelengths. $\Delta m_3$ is higher for shorter wavelengths and smaller $m_{\rm ej}$, being the difference small when varying $\phi$.

![Peak absolute magnitude (top) and decline rate at day 3 after peak (bottom) as a function of the ejecta mass in the models. Synthetic lightcurves at polar (left) and equatorial (right) viewing angles are shown, where the shaded areas represent variations in half-opening angles $\phi$ from 15 to 75$^{\circ}$. The stars indicate the values for the fiducial model.[]{data-label="fig:peak_dmag"}](figures/peak_dmag.pdf){width="\columnwidth"}

![Apparent $r$-band magnitude at peak vs redshift/distance for the entire grid of models. The orange region represents the fiducial model ($m_{\rm ej} = 0.06\,\mathrm{M_{\sun}}
    $ and $\phi=60^\circ$), but for different viewing angles, while the blue region the rest of the models for different $m_{\rm ej}$, $\phi$ and $\Theta$.[]{data-label="fig:mag_z"}](figures/mag_z_r.pdf){width="\columnwidth"}

Survey simulations {#sec:surveysimulations}
==================

We perform survey simulations using `simsurvey`[^2] [@simsurvey], a python based simulation tool for astronomical surveys. `simsurvey` uses Monte Carlo methods to generate synthetic lightcurves of transients, e.g., supernovae or kNe. It uses a timeseries of SEDs for the transient model, from which the photometry is calculated for a given filter, including cosmological redshift. The input for the simulation is a survey schedule containing the information about cadence, field of view, pointings, filter choice, sky brightness and zero points of each exposure. `simsurvey` is a tool developed within the ZTF collaboration, and it has been used successfully for forecasting rates of supernovae [@simsurvey].

Besides including the range of kN models with random orientations described above, we also account for dimming by dust in the interstellar medium of the host galaxy, as well as Galactic extinction based on the Schlegel, Finkbeiner, Davis (SFD98) reddening maps. In both cases, we use the standard Milky-Way extinction law for the wavelength dependence [@ccm89dust]. For this work, we have picked a region of the sky with a mean Milky Way extinction of $E(B-V)=0.08$ mag.

The output of `simsurvey` is a python dictionary that contains all the parameters of lightcurves for both detected and undetected transients.

Adopted survey plan {#subsec:survey_plan}
-------------------

We simulate ToO triggers under optimal conditions, i.e., shortly after the merger of two neutron stars. The synthetic observations start between 2.4 to 72 hours after the merger, and continue for seven days, with observations twice per night. Although we use ZTF $gri$ filters for the simulations, we expect little dependence on the specifics of these filters, i.e., our findings should be applicable to most optical follow-up campaigns during O3. We simulate lightcurves down to a hypothetical $m_{\rm lim} =$ 25 mag in all filters. Magnitudes are in AB system throughout this work.

Transient generation {#subsec:transient_generator}
--------------------

We generate a fixed number of synthetic kN lightcurves randomly distributed in the field of view of the instrument. Transients are placed at distances sampled from a distribution based on the volumetric rate as a function of redshift up to $z=0.15$.

The kN viewing angle is randomized, i.e., has constant probability for $\cos{\Theta}$. Thus, because of the assumed axial symmetry in the models, viewing angles closer to the merger plane are more probable to happen than those closer to the pole.

Detection criterion {#sec:detection_criteria}
-------------------

Next, we describe the adopted detection criterion. In dedicated searches of kNe following a GW alert, a basic requirement is non-detection [*prior*]{} to the merger time, as well as non-detection after sufficient time has elapsed, beyond which the lightcurves are expected to have faded. To have a non-detection we require a [*prior*]{} visit with no signal passing the threshold. However, in this study we simplify the situation in that we do not simulate data before the merger, nor after the lightcurves fade. Furthermore, we separate the synthetic observations by at least one hour, a common strategy to discard potential false positives from asteroids. We require at least two photometric detections with a SNR $\ge 5$.

LIGO/Virgo alerts in O3 {#sec:LVCtriggers}
=======================

LVC alerts involving neutron stars from LVC in O3, 6 BNS and 8 NS-BH mergers, were followed-up by ground-based optical/NIR telescopes around the globe, following the ability of each telescope to cover the merger location and weather conditions permitting. The searches were performed trying to optimize coverage of the candidate skymap provided by the LVC. Unfortunately, no identification of a kN was made, nor any other EM detection that could be securely associated to the GW signal. In this section, we explore the detection probability for kNe at the distances reported in the LVC alerts, with the purpose to explore possible explanations for the non-detections. Since the modelled grid are tailored to BNS mergers (Section \[sec:models\]), we will focus on the six BNS triggers listed in Table \[tab:listBNS\]. The efficiency for detection of NS-BH counterparts have to be treated separately, as their lightcurves are predicted to be brighter and longer lasting compared to those from BNS mergers [@Rosswog2017; @Barbieri2019; @Kawaguchi2020].

[*S190425z*]{} was the first BNS candidate in O3. The LALInference localization pipeline [@2015PhRvD..91d2003V] provided 7461 deg$^2$ in the 90% credible region and an estimated luminosity distance of 156 $\pm$ 41 Mpc. The day after, a new alert for a possible BNS was sent ([*S190426c*]{}), with a better localization (1131 deg$^2$) but more than twice further away (377 $\pm$ 100 Mpc). These first two candidates were followed-up within the GROWTH collaboration (@2019ApJ...885L..19C,@Goldstein19), probing the fast response to the low latency skymap and coverage capacity. No alerts with BNS candidates were released in the following few months until the beginning of September with [*S190901ap*]{} and [*S190910h*]{}. The skymap for both of these events was poorly localized,14753 deg$^2$ and 24264 deg$^2$, and the luminosity distances inferred to be 241 $\pm$ 79 Mpc and 233 $\pm$ 88 Mpc, respectively. These were the last events detected before the break in October and the end of O3a. During O3b, two non-retracted BNS candidates were sent, [*S191213g*]{} and [*S200213t*]{}, both at a mean distance of 201 Mpc ($\pm$ 81 and 80 Mpc, respectively). The closest BNS candidate in O3, at 156 Mpc, is four times more distant than GW170817 (40 Mpc), and the rest of the events happened more than five times further away. Another remarkable factor is the poor localization comparing with GW170817, 28 deg$^2$ 90% confident area [@2017PhRvL.119p1101A].

We seek to understand which is the limiting magnitude required to detect a kN. We study the detectability at different luminosity distances. We focus on the results at 200 Mpc as a representative distance for the BNS candidates in O3.

      BNS      $D_{\rm L}$ (Mpc)   $m_{{\rm lim}, gr}$   $m_{{\rm lim}, gri}$
  ----------- ------------------- --------------------- ----------------------
   S190425z           156                 24.5                    23
   S190426c           377                &gt;25                   25
   S190901ap          241                &gt;25                   24
   S190910h           233                &gt;25                   24
   S191213g           201                  25                    23.5
   S200213t           201                  25                    23.5

  : List of BNS merger candidates during LVC O3 run, together with their mean estimated distance and the minimum $m_{\rm lim}$ that gives detection probability more than 90%, integrating over all the viewing angles for the fiducial model. We show results for observations with $gr$ ($m_{{\rm lim}, gr}$) and $gri$ ($m_{{\rm lim}, gri}$).[]{data-label="tab:listBNS"}

Detection efficiencies for fiducial kN model
--------------------------------------------

The efficiency is calculated from the fraction of the synthetic lightcurve population that fulfill the adopted detection criterion. As expected, the kN detection probability improves for short distances and viewing angles closer to the pole, i.e., where the line of sight points to the lanthanide-free region. Fig. \[fig:hist2deff\] shows the detection efficiency for limiting magnitudes $21$, $22$ and $23$ mag, assuming the fiducial model matching AT2017gfo, where the only free parameters are distance and viewing angle. We also indicate the estimated mean distances from the GW signals of the six BNS trigger by LIGO/Virgo during O3. The shape of the $90 \%$ probability contours shows the dependency with the viewing angle for the fiducial model, the probability drops significantly towards equatorial viewing angles in which the photons are heavily absorbed by the lanthanides. As can be appreciated, unlike GW170817, none of the much more distant six LVC BNS events during O3 had a sizeable probability of detection at a limiting magnitude of $\lesssim$21 mag, which was very typical for the O3 follow-up campaigns. For $22$ mag, only favourable viewing angles can be probed, whereas surveys reaching $23$ mag would have had a good chance to detect a counterpart. We also find that redder wavelength coverage, like $i$-band reaching the same limiting magnitude, significantly increases the volume in which a kN is likely observable.

Table \[tab:listBNS\], shows the minimum limiting magnitude at which the integrated detection probability over all viewing angles is more than 90% at distance of the BNS merger candidates from LVC in O3 assuming the fiducial model.

![image](figures/cad0_6e-02_60_eff_hist21.pdf){width=".32\textwidth"} ![image](figures/cad0_6e-02_60_eff_hist22.pdf){width=".32\textwidth"} ![image](figures/cad0_6e-02_60_eff_hist23.pdf){width=".32\textwidth"}

Generalized kilonova models {#sec:all_grid}
---------------------------

The probability of kN detection depends on the assumed physical parameters of the model. Fig. \[fig:p\_22\] indicates the estimated probabilities at 200 Mpc for observations with $gr$ and $gri$, all at an assumed limiting magnitude of 22 mag. Results for other limiting magnitudes can be found in the appendix \[app:det\_prob\]. The probabilities are calculated for two subsets of detections: events with viewing angles pointing to the lanthanide-free component $\Theta_{\rm LF}= \Theta < 90^{\circ}º - \phi$ or to the lanthanide-rich component $\Theta_{\rm LR}=\Theta > 90^{\circ}º - \phi$. Table \[tab:det\_prob\] shows probabilities at 40 Mpc (GW170817), 100 Mpc and 200 Mpc, for the fiducial model and two extremes cases, i.e., models that provide the brightest and the faintest kN.

The detection probability strongly depends on the physical parameters of the kN and the orientation, i.e., large $m_{\rm ej}$, small $\phi$, and $\Theta_{\rm LF}$ are associated with brighter kNe and thus provide higher probability values.

The observations are constraining on the parameter space of the model that provides high detection probability. We aim to be sensitive to a wide range of parameters with our observations. The inclination of the merger is unknown when the ToO is triggered. Viewing angles close to face-on provide the brighter lightcurves, but they are less likely to happen (see section \[subsec:transient\_generator\]). Therefore, it is preferable to plan for less desirable circumstances. 23 mag depth in $gri$ is needed to constraint most of the parameter space at 200 Mpc.

  ----- ----- ----------- ----------- ------------ ------------ ----------- ----------- ------------ ------------ ----------- ----------- ------------ ------------
               fiducial                                          brightest                                         faintest                            
   Mpc   lim   LF$_{gr}$   LR$_{gr}$   LF$_{gri}$   LR$_{gri}$   LF$_{gr}$   LR$_{gr}$   LF$_{gri}$   LR$_{gri}$   LF$_{gr}$   LR$_{gr}$   LF$_{gri}$   LR$_{gri}$
   40    21      1.00        0.73         1.00         0.95        1.00        1.00         1.00         1.00        0.74        0.29         0.83         0.57
   40    22      1.00        0.90         1.00         1.00        1.00        1.00         1.00         1.00        0.84        0.50         0.90         0.80
   40    23      1.00        0.99         1.00         1.00        1.00        1.00         1.00         1.00        0.87        0.71         0.98         0.94
   100   21      0.79        0.27         0.96         0.59        0.98        0.90         1.00         1.00        0.17        0.00         0.30         0.11
   100   22      0.97        0.56         1.00         0.89        1.00        1.00         1.00         1.00        0.33        0.06         0.51         0.40
   100   23      1.00        0.81         1.00         1.00        1.00        1.00         1.00         1.00        0.46        0.31         0.73         0.68
   200   21      0.12        0.00         0.36         0.04        0.61        0.37         0.91         0.75        0.00        0.00         0.00         0.00
   200   22      0.54        0.10         0.87         0.36        0.95        0.81         1.00         1.00        0.00        0.00         0.00         0.00
   200   23      0.92        0.38         1.00         0.76        1.00        0.92         1.00         1.00        0.06        0.03         0.31         0.24
  ----- ----- ----------- ----------- ------------ ------------ ----------- ----------- ------------ ------------ ----------- ----------- ------------ ------------

![image](figures/p_200_22.pdf){width="\textwidth"}

Synthetic kilonova population detected in simulations {#sec:synthetic_population}
=====================================================

Next, we explore the model properties for the population of kNe that pass the detection criterion for a given magnitude limit. We group the viewing angles of the detected sample in two sets: one with orientation such that the line of sight points to the lanthanide-free component ($\Theta_{\rm LF}$) and the other pointing to the lanthanide-rich component ($\Theta_{\rm LR}$), defined in Section \[sec:all\_grid\]. We also separate the detections for synthetic observations into those that include $i$-band and those that do not. To better assess the selected samples under different considerations, we scale the number of detections to a BNS rate of 1000 Gpc$^{-1}$yr$^{-1}$ [@Scolnic2017] for 20000 sq. degrees. The summary of our findings is shown in Table \[tab:output\_sim\].

   $m_{\rm lim}$   $n_{{\rm KNe}, gr}$   $n_{{\rm KNe}, gri}$   $\Theta_{{\rm LF}, gr}$   $\Theta_{{\rm LF}, gri}$   $\bar{\Theta}_{gr}$   $\bar{\Theta}_{gri}$   $\bar{D}_{{\rm L}, gr}$   $\bar{D}_{{\rm L}, gri}$
  --------------- --------------------- ---------------------- ------------------------- -------------------------- --------------------- ---------------------- ------------------------- --------------------------
       20.5                 1                     3                      0.45                       0.33                     36                     43                      65                         78
        21                  2                     5                      0.38                       0.32                     39                     43                      100                       125
        22                  8                     20                     0.36                       0.33                     40                     43                      156                       196
        23                 28                     71                     0.38                       0.33                     39                     43                      245                       302
        24                 96                    213                     0.38                       0.28                     40                     45                      369                       432
        25                 242                   411                     0.26                       0.18                     46                     53                      458                       502

Given the red nature of kNe, including $i$-band to the same depth as $g$ and $r$ in the survey leads to more than 50% percent improvement on the number of detections. kN lightcurves show a slower decay rate in redder bands (smaller $\Delta m_3$, see Fig. \[fig:peak\_dmag\]), they are less affected by the flux suppression of the lanthanide-rich component and are less affected by dust extinction.

The increase in detections with limiting magnitude is exponential, as more volume becomes available, . In our simulations, we find that $g$-band observations, although useful for background rejection and characterization of the kN properties, do not add to the detection efficiency. All the kNe detected with SNR $\ge 5$ do also pass the detection criterion in $r$ and $i$-band.

Selection effects {#sec:selection_effects}
-----------------

The distribution of detections in viewing angle and distance is shown in Fig. \[fig:hist2d\]. The detected population is affected by Malmquist bias [@Malmquist], i.e., given a magnitude limit, only the brightest objects in the population are detected. For the synthetic population, as volume increases with distance, most kNe are generated at farther distances. Furthermore, equatorial viewing angles are more common (see Section \[subsec:transient\_generator\]). Thus, the brightest events correspond to lower distances and polar viewing angles. The combination of the two effects gives a median viewing angle of about 40$^{\circ}$ for the population of selected kNe, assuming the fiducial model.

![image](figures/det_fid.pdf){width="\textwidth"}

When searching for transients in the survey alerts, astronomers look for at least two data points with a sufficient signal to noise and time separation to consider a reliable discovery. The chances to observe a lightcurve that can be identified as a kN increase when the transient is observable for a more extended period. The duration of the observable lightcurve depends on the properties of the kN and vary with distance. Fig. \[fig:last5sigma\] shows the range from the minimum to the maximum length of the detected lightcurves per distance for the fiducial model, being length defined as the post-merger epoch of the last 5$\sigma$ data point.

![Range from the minimum to the maximum last detection (in units of days) from the detected lightcurves as a function of distance, observing with $gr$ and $gri$ at different limiting magnitudes for the fiducial model.[]{data-label="fig:last5sigma"}](figures/last5sigma_redshift_range.pdf){width="\columnwidth"}

Feasibility for serendipitous detections {#sec:rates}
========================================

Next, we explore the feasibility for kN detections in time-domain surveys, irrespective of external GW triggers, for a volumetric BNS rate of 1000 Gpc$^{-3}$yr$^{-1}$, consistent with [@Scolnic2017]. We make several simplifying assumptions, including losses due to candidate selection inefficiencies, bad weather, or other observing conditions. Every kNe observable within the covered volume with lightcurves of at least two 5$\sigma$ data points with a minimum separation of one hour is detected. Furthermore, a constant magnitude limit is assumed for the entire duration of an all-sky survey, i.e., 20000 sq.degrees. We assume that the distribution of BNS is homogeneous on the sky and rescale the number of detections using efficiencies calculated in Section \[sec:LVCtriggers\]. With the fiducial model, we find that for $m_{\rm lim}$ = 20.5 mag, 3 kNe are observable per year, around 20 for limiting magnitude of 22 mag and 71 for 25 mag, assuming constant depth with $gri$ wavelength coverage (see Fig. \[fig:nkne\] and Table \[tab:output\_sim\]).

![Expected number of KNe passing our detection criterion for a survey covering 20000 square degrees of the sky, assuming constant limiting magnitudes from 20.5 to 23 mag, observing with $gr$ (green) and $gri$ (red) filters. A BNS rate of $1000$ Gpc$^{-3}$yr$^{-1}$ is assumed for this calculation [@Scolnic2017]. The solid lines correspond to the fiducial model, while the shaded areas show the values for the entire grid of parameters.[]{data-label="fig:nkne"}](figures/n_kne_yr.pdf){width="\columnwidth"}

![Probability for KN detection (horizontal axis) and lightcurve duration (vertical axis) for different limiting magnitudes (different lines) and assuming a source at 200 Mpc. The post-merger epoch in units of days indicates the last data point (with SNR$\ge5$) in the lightcurves. The green area on the right hand side of the plot highlights the probability values above 0.9, i.e., likely detections. The curves are computed from the average over the entire population of detected synthetic lightcurves for the fiducial model observed with $gri$.[]{data-label="fig:last5p"}](figures/p_last.pdf){width="\columnwidth"}

Summary and Conclusion {#sec:discussion}
======================

The detection of GW-EM counterparts has proven to be very challenging. During the LIGO/Virgo O3 run, six BNS triggers were issued, and in spite of large efforts from the entire time-domain astronomy community, no transient resembling AT2017gfo was found. Furthermore, fast-evolving and faint kN transients are too rare to facilitate serendipitous discoveries. The fast evolution in bolometric luminosity and color typically leads to lightcurves with few data points, hard to identify photometrically, and even more to get a good quality spectrum. In some cases, the lightcurve would fade below threshold within the first day from the merger. Fig. \[fig:last5p\] shows the detection probability for various lengths of lightcurve coverage demanding SNR $\ge 5$, for limiting magnitudes ranging from 21 to 25 mag. For a given limiting magnitude, the efficiency is reduced when requiring longer lightcurves. E.g., for $m_{\rm lim} = 23$ mag, there is more than 90% probability for detecting a kilonova at 200 Mpc over 3 days since the merger. Requiring a duration of 5 days of the observable lightcurve, the efficiency drops down to about 10%. For viewing angles towards the lanthanide-free component, reaching 22 mag is required in $gri$, and 23 mag for the viewing angles towards the lanthanide-rich component, to reach 90% probability with a 3 days lightcurve.

In addition to the intrinsic observational challenges posed by the rapid evolution of kN lightcurves, the distances to the BNS triggers during the third run of LIGO and Virgo have been around 200 Mpc, about five times larger than for GW170817. Some of the events were detected with a single interferometer, which translates into poor localization. Large search areas, in combination with long distances, poses great challenges for the detection of EM counterparts. Survey scheduling under such conditions is even harder, considering the lack of knowledge of the underlying kN population. During O3, the effort to cover these large search areas on the sky was remarkable (@2020arXiv200400025G; @GRANDMA; @2020ApJ...890..131A; @Ackley2020; @2020arXiv200309437V; @2020MNRAS.492.5916W; @2020arXiv200400025G; @Goldstein19), but the adopted depth was generally insufficient to identify the counterpart, according to the models we have studied. Several collaborations such as DECam-GROWTH[^3] [@GROWTH_DECam-GROWTH], ENGRAVE [@Ackley2020] or GRANDMA [@antier2020; @GRANDMA] had triggered ToO observations with a global multi-telescope network when it was feasible. It has been shown that current capabilities can afford to cover large areas on the sky, but the reached depth was insufficient to detect the EM.

Clearly, GW170817 appears to have been an exceptionally lucky occurrence, both with respect to distance, localisation and orientation, well outside the bulk of the BNS mergers in O3 (see Fig. \[fig:hist2d\]). Due to their low rate, e.g., compared with other types of transients like supernovae, the probability of another nearby kN with viewing angle close to the pole, is very low. The asymmetry of the kN emission and the higher likelihood of viewing angles close to the merger plane means that most of the kNe are going to be produced under conditions much less favourable for detection than GW170817.

In agreement with previous studies of follow up campaigns during O3 (@2020MNRAS.492..863C,@2019ApJ...880L...4H), we conclude that deeper observations are required to ensure high detection probability. At a distance of 200 Mpc, the required magnitude to observe a kN similar to AT2017gfo at lanthanide-rich viewing angles with three days lasting lightcurve is at least 23 mag. For an AT2017gfo-like event at that distance, pointing towards us with a viewing angle within the lanthanide-free component, $m_{\rm lim} = 22$ mag could be sufficient. This values change for different physical parameters of the kN model (see Section \[sec:lcs\]). The optical lightcurves of kNe become very red in just a few days. The nucleosynthesis of heavy nuclei translates into longer-lived emission at longer wavelengths. The blue emission fades rapidly, and it suffers more from absorption at viewing angles close to the merger plane, where the dynamical ejecta leaves most of the lanthanides. Thus, observations in red filters tend to be much more efficient. In this work we show that adding $i$-band observations in an optical survey yields, on average, more than two times more detections than e.g., observing with only $g$- and $r$-band. The $g$-band lightcurves fade very quickly and, in some cases, would not even provide any detection. However, observations in bluer bands may play important roles in candidate vetting and characterization of kN events. Detection of the very fast decay in bluer bands, besides the longer lasting redder wavelengths is an excellent feature for the identification.

On going ground-based telescopes can reach 23 mag under good atmospheric conditions and enough integration time on the source. However, due to the very fast evolving luminosities, continued deep observations over multiple days is required, something that is hard to achieve with current instruments due to unpredictable weather patterns and moon phase among other factors. The Vera C. Rubin Observatory (VRO/LSST) with the Legacy Survey of Space and Time [LSST, @lsst] with its 8-meter-class telescope, planned to be online by 2023, will discover transients that are too far or too faint for the current available surveys, including serendipitous discoveries of kNe [@setzer19].

In this work, we have used models tailored to BNS. However, we expect that results for NS-BH lightcurves would change quantitatively but not qualitatively. In which case, the lack of EM counterpart detections for the distant NS-BH LVC candidates during O3 would also be justified. In conclusion, detection and identification of GW electromagnetic counterparts from BNS mergers is very challenging. The survey magnitude depths during follow-up campaigns in O3 have been generally insufficient, considering the asymmetry of the kN emission, the sometimes poor localization and distances, typically 5 times further than GW170817. Despite the lack of kN detections during O3, the multimessenger astronomy community have demonstrated a remarkable ability to respond promptly to LVC alerts, even those with wide search areas. Thus, the future is promising, also as more detectors join the gravitational waves observatory network, leading to better localization. Thus, we expect that with improved search strategies described here, future campaigns will be successful.

Acknowledgements {#acknowledgements .unnumbered}
================

The authors acknowledge support from the G.R.E.A.T research environment, funded by [*Vetenskapsrådet*]{}, the Swedish Research Council.

detection probability {#app:det_prob}
=====================

In section \[sec:all\_grid\], we discuss how likely is a kilonova at 200 Mpc to be detected assuming constant depth in $gr$ and $gri$. We explore the entire parameter space of the models for ejecta masses ($m_{\rm ej}$) varying between 0.01 to 0.10 $M_{\sun}$, half-opening angles ($\phi$) from 15 to 75 degrees and orientations towards the lanthanide-free ($\Theta_{\rm LF}$) or the lanthanide-rich ($\Theta_{\rm LR}$) components. Fig.  \[fig:p\_22\] presents the results for $m_{\rm lim} = 22$. Fig. \[fig:p\_21\_23\] shows results for $m_{\rm lim} = 21$ and $m_{\rm lim} = 23$ mag.

![image](figures/p_200_21.pdf){width="\textwidth"} ![image](figures/p_200_23.pdf){width="\textwidth"}

\[lastpage\]

[^1]: E-mail: ana.sagues-carracedo@fysik.su.se

[^2]: <https://github.com/ZwickyTransientFacility/simsurvey>

[^3]: <http://growth.caltech.edu/>
