# Observational constraints on the optical and near-infrared emission from the neutron star-black hole binary merger S190814bv

F. d'Ammando, P. d'Avanzo, V. d'Elia, A. Nicuesa Guelbenzu, P. T. O'Brien, A. Sagués Carracedo, L. Wyrzykowski, K. Ackley, L. Amati, C. Barbieri, et al.

► To cite this version:

F. d'Ammando, P. d'Avanzo, V. d'Elia, A. Nicuesa Guelbenzu, P. T. O'Brien, et al. Observational constraints on the optical and near-infrared emission from the neutron star-black hole binary merger S190814bv. Astronomy and Astrophysics - A&A, 2020, 643, pp.A113. 10.1051/0004-6361/202037669. hal-02904245

HAL Id: hal-02904245

https://hal.science/hal-02904245

Submitted on 5 May 2023

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
# Observational constraints on the optical and near-infrared emission from the neutron star–black hole binary merger candidate S190814bv*

K. Ackley¹, L. Amati², C. Barbieri³,⁴,⁵, F. E. Bauer⁶,⁷,⁸, S. Benetti⁹, M. G. Bernardini³, K. Bhirombhakdi¹⁰,
M. T. Botticella¹¹, M. Branchesi¹²,¹³, E. Brocato¹⁴,¹⁵, S. H. Bruun¹⁶, M. Bulla¹⁷, S. Campana³, E. Cappellaro⁹,
A. J. Castro-Tirado¹⁸, K. C. Chambers¹⁹, S. Chaty²⁰,²¹, T.-W. Chen²²,²³, R. Ciolfi⁹,²⁴, A. Coleiro²⁰,
C. M. Copperwheat²⁵, S. Covino³, R. Cutter²⁶, F. D'Ammando²⁷, P. D'Avanzo³, G. De Cesare², V. D'Elia¹⁴,²⁸,
M. Della Valle¹¹, L. Denneau¹⁹, M. De Pasquale²⁹, V. S. Dhillon³⁰,³¹, M. J. Dyer³⁰, N. Elias-Rosa⁹,³², P. A. Evans³³,
R. A. J. Eyles-Ferris³³, A. Fiore⁹,³⁴, M. Fraser³⁵, A. S. Fruchter¹⁰, J. P. U. Fynbo³⁶,³⁷, L. Galbany³⁸, C. Gall¹⁶,
D. K. Galloway¹, F. I. Getman¹¹, G. Ghirlanda³, J. H. Gillanders³⁹, A. Gomboc⁴⁰, B. P. Gompertz²⁶,
C. González-Fernández⁴¹, S. González-Gaitán⁴², A. Grado¹¹, G. Greco⁴³,⁴⁴, M. Gromadzki⁴⁵, P. J. Groot⁴⁶,⁴⁷,⁴⁸,⁴⁹,
C. P. Gutiérrez⁵⁰, T. Heikkilä⁵¹, K. E. Heintz⁵², J. Hjorth¹⁶, Y.-D. Hu¹⁸,³⁸, M. E. Huber¹⁹, C. Inserra⁵³, L. Izzo¹⁶,
J. Japelj⁵⁴, A. Jerkstrand²², Z. P. Jin⁵⁵, P. G. Jonker⁴⁶,⁵⁶, E. Kankare⁵¹, D. A. Kann¹⁸, M. Kennedy⁵⁷, S. Kim⁶,⁵⁸,
S. Klose⁵⁹, E. C. Kool²², R. Kotak⁵¹, H. Kuncarayakti⁵¹,⁶⁰, G. P. Lamb³³, G. Leloudas⁶¹, A. J. Levan²⁶,⁴⁶,★★,
F. Longo⁶², T. B. Lowe¹⁹, J. D. Lyman²⁶, E. Magnier¹⁹, K. Maguire⁶³, E. Maiorano², I. Mandel¹,⁶⁴, M. Mapelli³⁴,
S. Mattila⁵¹, O. R. McBrien³⁹, A. Melandri³, M. J. Michałowski⁵⁶, B. Milvang-Jensen³⁶,³⁷, S. Moran⁵¹, L. Nicastro²,
M. Nicholl⁶⁴, A. Nicuesa Guelbenzu⁵⁹, L. Nuttal⁶⁶, S. R. Oates²⁶,⁶⁴, P. T. O'Brien³³, F. Onori⁶⁷, E. Palazzi²,
B. Patricelli¹⁴,⁶⁸,⁶⁹, A. Perego⁷⁰, M. A. P. Torres³¹,⁵⁶,⁷¹, D. A. Perley²⁵, E. Pian², G. Pignata⁷,⁷², S. Piranomonte¹⁴,
S. Poshyachinda⁷³, A. Possenti⁷⁴,⁷⁵, M. L. Pumo⁹,⁷⁶,⁷⁷, J. Quirola-Vásquez⁶,⁷, F. Ragosta¹¹,⁷⁸,⁷⁹, G. Ramsay⁸⁰,
A. Rau²³, A. Rest¹⁰,⁸¹, T. M. Reynolds⁵¹, S. S. Rosetti³³, A. Rossi², S. Rosswog²², N. B. Sabha⁸²,
A. Sagués Carracedo⁸³, O. S. Salafia³,⁴, L. Salmon³⁵, R. Salvaterra⁸⁴, S. Savaglio⁸⁵, L. Sbordone⁸⁶, P. Schady⁸⁷,
P. Schipani¹¹, A. S. B. Schultz¹⁹, T. Schweyer²²,²³, S. J. Smartt³⁹, K. W. Smith³⁹, M. Smith⁵⁰, J. Sollerman²²,
S. Srivastav³⁹, E. R. Stanway²⁶, R. L. C. Starling³³, D. Steeghs²⁶, G. Stratta²,⁴⁴, C. W. Stubbs⁸⁸, N. R. Tanvir³³,
V. Testa¹⁴, E. Thrane¹, J. L. Tonry¹⁹, M. Turatto⁹, K. Ulaczyk²⁶,⁴⁵, A. J. van der Horst⁸⁹,⁹⁰, S. D. Vergani⁹¹,
N. A. Walton⁴¹, D. Watson³⁶,³⁷, K. Wiersema²⁶,³³, K. Wiik⁵¹, Ł. Wyrzykowski⁴⁵, S. Yang²²,
S.-X. Yi⁴⁶, and D. R. Young³⁹

*(Affiliations can be found after the references)*

Received 5 February 2020 / Accepted 18 June 2020

## ABSTRACT

**Context.** Gravitational wave (GW) astronomy has rapidly reached maturity, becoming a fundamental observing window for modern astrophysics. The coalescences of a few tens of black hole (BH) binaries have been detected, while the number of events possibly including a neutron star (NS) is still limited to a few. On 2019 August 14, the LIGO and Virgo interferometers detected a high-significance event labelled S190814bv. A preliminary analysis of the GW data suggests that the event was likely due to the merger of a compact binary system formed by a BH and a NS.
**Aims.** In this paper, we present our extensive search campaign aimed at uncovering the potential optical and near infrared electromagnetic counterpart of S190814bv. We found no convincing electromagnetic counterpart in our data. We therefore use our non-detection to place limits on the properties of the putative outflows that could have been produced by the binary during and after the merger.

**Methods.** Thanks to the three-detector observation of S190814bv, and given the characteristics of the signal, the LIGO and Virgo Collaborations delivered a relatively narrow localisation in low latency – a 50% (90%) credible area of 5 deg² (23 deg²) – despite the relatively large distance of 267 ± 52 Mpc. Electromagnetic counterparts of GRAVITATIONAL WAVE SOURCES AT THE V-Ery Large Telescope collaboration members carried out an intensive multi-epoch, multi-instrument observational campaign to identify the possible optical and near infrared counterpart of the event. In addition, the ATLAS, GOTO, GRAWITA-VST, PAN-STARRS, and VINROUGE projects also carried out a search on this event. In this paper, we describe the combined observational campaign of these groups.

**Results.** Our observations allow us to place limits on the presence of any counterpart and discuss the implications for the kilonova (KN), which was possibly generated by this NS–BH merger, and for the strategy of future searches. The typical depth of our wide-field observations, which cover most of the projected sky localisation probability (up to 99.8%, depending on the night and filter considered), is r ~ 22 (resp. K ~ 21) in the optical (resp. near infrared). We reach deeper limits in a subset of our galaxy-targeted observations, which cover a total ~50% of the

<footnote>* Full Table 3 is only available at the CDS via anonymous ftp to cdsarc.u-strasbg.fr (130.79.128.5) or via http://cdsarc.u-strasbg.fr/viz-bin/cat/J/A+A/643/A113</footnote>

<footnote>** Corresponding author: A. J. Levan, e-mail: a.levan@astro.ru.nl</footnote>
galaxy-mass-weighted localisation probability. Altogether, our observations allow us to exclude a KN with large ejecta mass $M \ge 0.1 M_\odot$ to a high (>90%) confidence, and we can exclude much smaller masses in a sub-sample of our observations. This disfavours the tidal disruption of the neutron star during the merger.

**Conclusions.** Despite the sensitive instruments involved in the campaign, given the distance of S190814bv, we could not reach sufficiently deep limits to constrain a KN comparable in luminosity to AT 2017gfo on a large fraction of the localisation probability. This suggests that future (likely common) events at a few hundred megaparsecs will be detected only by large facilities with both a high sensitivity and large field of view. Galaxy-targeted observations can reach the needed depth over a relevant portion of the localisation probability with a smaller investment of resources, but the number of galaxies to be targeted in order to get a fairly complete coverage is large, even in the case of a localisation as good as that of this event.

**Key words.** gravitational waves – stars: neutron – supernovae: general

# 1. Introduction

The discovery of the binary black hole (BH) merger event GW150914 (Abbott et al. 2016) was a major landmark in the history of physics. It was the first detection of gravitational waves (GWs) and the beginning of GW astronomy. The detection of the first confirmed binary neutron star (NS) merger, GW170817 (Abbott et al. 2017a), and the subsequent discovery of its electromagnetic (EM) counterparts – the short GRB 170817A (Abbott et al. 2017b; Goldstein et al. 2017; Savchenko et al. 2017) and the UV, optical, and IR transient AT2017gfo (Coulter et al. 2017a; Lipunov et al. 2017; Tanvir et al. 2017; Soares-Santos et al. 2017; Valenti et al. 2017) – was a second major breakthrough, and marked the beginning of multi-messenger astrophysics with GWs (Abbott et al. 2017c).

The subsequent investigation of GW170817 convincingly linked NS-NS mergers with short duration gamma-ray bursts (e.g. Lyman et al. 2018; Dobie et al. 2018; Mooley et al. 2018; Lazzati et al. 2018; Resmi et al. 2018; Ghirlanda et al. 2019; Lamb et al. 2019a; Margutti et al. 2018; Nynka et al. 2018; Troja et al. 2018a; D'Avanzo et al. 2018) – a link for which the evidence had been accumulating for some time (Fong & Berger 2013; Tanvir et al. 2013; Berger et al. 2013; Jin et al. 2016). In addition, the identification of its host galaxy, NGC 4993 (see Levan et al. 2017), and an assessment of its cosmological recession velocity (Hjorth et al. 2017) permitted the first measurement of a cosmological parameter (the Hubble constant) using the GW distance measurement, thanks to the “standard siren” nature of compact object binaries (Abbott et al. 2017d). The optical and near-infrared (NIR) monitoring campaigns of the transient also unveiled, for the first time, the developing kilonova (KN) emission (Arcavi et al. 2017a; Chornock et al. 2017; Covino et al. 2017; Cowperthwaite et al. 2017; Drout et al. 2017; Evans et al. 2017; Kasliwal et al. 2017; McCully et al. 2017; Nicholl et al. 2017; Pian et al. 2017; Shappee et al. 2017; Smartt et al. 2017; Tanvir et al. 2017) due to the production and decay of r-process elements (e.g. Metzger et al. 2010; Kasen et al. 2017), demonstrating that NS-NS mergers are indeed a major source of these elements (Gall et al. 2017; Watson et al. 2019), as previously suggested (Lattimer et al. 1977; Eichler et al. 1989; Li & Paczyński 1998; Freiburghaus et al. 1999).

Following the success of the GW170817 follow-up campaign, considerable effort has been expended in mounting similar campaigns with the aim of discovering and characterising the counterparts of new GW events. To optimise the science return of the demanding observations of GW counterparts, a large fraction of the GW and EM community in member states of the European Southern Observatory (ESO) has gathered together to form the ElectromagNetic counterparts of GRAvitational wave sources at the VEry Large Telescope (ENGRAVE) consortium<fnref n="1" />. This paper introduces the collaboration and our first major campaign to search for an EM counterpart to a GW source,

the NS-BH event merger candidate S190814bv, reported during the O3 run of the LIGO Scientific Collaboration and the Virgo Collaboration (LVC)<fnref n="2" />.

S190814bv was detected by the LVC on 2019 Aug 14 21:10:39.01 UT (MJD 58709.88240; The LIGO Scientific Collaboration and the Virgo Collaboration 2019a), and an alert was issued on 2019 Aug 14 21:31:40 UT (MJD 58709.89699), approximately 21 min after the merger. The source was localised to a 50% credible region of 133 deg<fnref n="2" /> in the initial report (90% credible region of 772 deg<fnref n="2" />), which was reduced to 5 deg<fnref n="2" /> at 50% and 23 deg<fnref n="2" /> at 90% half a day later (The LIGO Scientific Collaboration and the Virgo Collaboration 2019b) making this the best-localised candidate GW event so far. The source distance (as inferred directly from the GW observations) is 267 ± 52 Mpc. The estimated False-Alarm Rate (FAR) is extremely low, at $2.033 \times 10^{-33}$ Hz (1 per 1.559 × 10²⁵ years).

Preliminary parameter estimation (Veitch et al. 2015) indicated that the lighter object had a mass $M_2 < 3 M_\odot$, while the heavier object had a mass $M_1 > 5 M_\odot$, making this a NS-BH candidate according to the LVC classification criteria. The same preliminary analysis pointed to a negligible probability of any disrupted material remaining outside the final compact object (given by the parameter HasRemnant <1%), implying that in this case an EM counterpart was unlikely. We note that the classification of one of the components as a NS is based solely on the mass being <3 $M_\odot$ and that a low mass BH is not ruled out by such low-latency classification.

No γ-ray, X-ray, or neutrino signal could be connected to the event (Molkov et al. 2019; Kocevksi 2019; Ohno et al. 2019; Pilia et al. 2019a,b; Sugizaki et al. 2019; Palmer et al. 2019; Evans et al. 2019; Cai et al. 2019; Svinkin et al. 2019; Ice-Cube Collaboration 2019; Ageron et al. 2019; Alvarez-Muniz et al. 2019). The relatively small localisation region led to a world-wide follow-up effort with optical and NIR telescopes (e.g. Gomez et al. 2019a; Andreoni et al. 2020; Dobie et al. 2019; Watson et al. 2020; Antier et al. 2020; Vieira et al. 2020). The ENGRAVE collaboration activated its search programmes to try to discover, or set limits on, an EM counterpart to S190814bv. ENGRAVE members ran wide-field searches for EM counterparts and the ATLAS, GOTO, GRAWITA-VST, Pan-STARRS, and VINROUGE projects also triggered their searches on this event. No promising EM counterpart was detected. In this paper we combine our ENGRAVE ESO/Very Large Telescope (VLT) data, a number of other narrow-field facilities, and the wide-field programmes to report the combined results of our search for a counterpart (Sects. 2 and 3). We place limits on the presence of a counterpart (Sect. 3.6), and discuss the implications of these limits for NS-BH mergers and future searches (Sect. 4). Unless otherwise specified, errors are given at 68% confidence level ($1\sigma$), upper limits are given at $3\sigma$, and magnitudes are in the AB system. When needed, we assume a flat FLRW cosmology with $H_0 = 70 \text{ km s}^{-1} \text{ Mpc}^{-1}$ and $\Omega_m = 0.3$.

<footnote n="1">¹ http://www.engrave-eso.org/</footnote>

<footnote n="2">² https://gracedb.ligo.org/superevents/public/03/</footnote>
**2. Wide-field survey observations and results**

We employed two different approaches to search for an optical or NIR counterpart to S190814bv. A number of wide-field facilities (with Field of View, FoV, of 1 deg² or more) were used to tile the LVC GW sky localization probability maps (skymap) with the aim of covering as much of the 2D probability of S190814bv localisation as possible. These telescopes (with apertures of 0.4 m–4 m) were situated in La Palma, Chile and Hawaii, giving a spread of latitude and temporal coverage. The second approach was to target specific galaxies in the 3D sky region with larger aperture (2 m–8 m) telescopes and smaller FoV cameras. In the next subsections we summarise the search for transients with the different facilities.

**2.1. The search for transients with GOTO**

The Gravitational wave Optical Transient Observer (GOTO<fnref n="3" />) is a robotic array of wide-field optical telescopes sited at the Roque de los Muchachos Observatory in La Palma. It is operated by the University of Warwick on behalf of an international collaboration. The hardware is modular in design and optimised to autonomously respond to GW events, being able to cover large areas of sky quickly. At the time of S190814bv, GOTO was equipped with 4 active unit telescopes, each having an aperture of 40 cm at *f*/2.5 and featuring a 50-megapixel CCD detector. This corresponds to a plate-scale of 1′.25 pix⁻¹ and a FoV of 5.9 deg² per camera. Observations of the S190814bv error box were automatically scheduled on the basis of a ranked tiling pattern derived from the available skymaps (Dyer et al. 2018; Gompertz et al. 2020a). Each tile was observed using sequences of 60 s or 90 s exposures with the GOTO-*L* filter, which is a wide filter covering 400–700 nm (slightly wider than Sloan *g*+*r* combined). The bulk of the final localisation probability, as given by the *LALInference* skymap (*LALInference.v1.fits*), was covered in 8 observable tiles (Fig. 1), with 89.6% probability covered over the timespan MJD = 58710.091–58710.230 (5.09–8.34 h after the GW event). Additional observations were obtained the following night, which covered 94.1%. The Moon was closer to the relevant tiles during this second night, affecting the zeropoints achieved in the exposures sets. Most tiles were observed multiple times (see Table B.2), though observing conditions were not optimal, given the presence of the Moon nearby and poor weather. The probability regions were also close to GOTO’s lower declination limit, meaning some of the GW probability region could not be observed. Individual exposures are median combined in groups of 3–6 subsequent images and reached 5σ limiting magnitudes covering *V* = 17.4–19.1 mag. These are derived by photometrically calibrating our photometry to AAVSO Photometric All-Sky Survey (APASS<fnref n="4" />) stars using the *V* band.

Raw GOTO images are transferred from La Palma to the Warwick data centre in real-time and processing begins minutes after acquisition using the GOTOphoto pipeline (GOTO collaboration, in prep.). Image-level processing includes detector corrections, astrometry tied to Gaia and photometric zero-points using a large number of field stars. Difference imaging was performed on the median exposures using recent survey observations as reference. Source candidates were initially filtered using a trained classifier and cross-matched against a variety of catalogues, including the Minor Planet Center (MPC)<fnref n="5" /> and

Pan-STARRS1 3π survey (Chambers et al. 2016). The classifier employs a random forest algorithm based on image features, largely following the procedure of Bloom et al. (2012). It was trained using an injected source data-set. Human vetting was performed on the resulting candidates using a web-based marshall interface. No viable optical counterpart candidates could be associated with S190814bv (Ackley et al. 2019).

**2.2. The search for transients with the VST**

The VLT Survey Telescope (VST; Capaccioli & Schipani 2011) is a 2.6 m facility located at Cerro Paranal, Chile, and managed by ESO. The telescope is equipped with OmegaCam (Kuijken 2011), a 268-megapixel camera with a FoV of 1 × 1 deg² and a resolution of 0′.21 pix⁻¹. The filter set includes Sloan *ugriz* filters. The telescope is operated in service mode. The GW trigger follow-up is performed using Guaranteed Time allocated to the Italian VST and OmegaCam Consortium.

The monitoring of the S190814bv sky area started on MJD 58710.36 (Grado et al. 2019a,b; Yang et al. 2019). For each pointing and epoch we obtained three dithered exposures for a total exposure time of 135 s. The pointings were visited up to five times during a period of two weeks (see Table B.2). All exposures were obtained using the *r* filter.

On the first night we imaged 15 deg², covering 53.6% of the localisation probability of the preliminary BAYESTAR skymap (which was the only map available at that time), and 60.7% of the final LALInference skymap probability. Starting from the second epoch we revised the pointing list to optimise the sky coverage for the updated LALInference skymap. The survey area increased to 23 deg² covering a maximum of 87.7% of the localisation probability, as shown in Fig. 1.

Details of the image processing and candidate detection are given in Brocato et al. (2018) and Grado et al. (2020). The area identified by the skymap is not fully covered by VST archive observations, hence we used both Pan-STARRS1 (Chambers et al. 2016) and DECam (Abbott et al. 2018a) archive images as templates for the comparison. To select the candidates, we applied a random-forest machine learning algorithm trained on previous search instances (Yang 2018) and then visually inspected the candidates with the highest score. We detected a number of transients from which we removed objects detected only at one epoch and/or associated with stellar sources in the template images. The final list includes 27 transients (reported in Table 2). Out of these, 21 objects were already discovered by other surveys and had been registered on the on the Transient Name Server (TNS)<fnref n="6" />. In Table 2 we include also three transients reported in TNS by other groups that are detected on our images but are below the detection threshold of our search and were therefore not independently detected in our search.

All of these candidates show a slow evolution in the two weeks of the observing campaign (Δ*m* < 1 mag between the first and last detection). Therefore we tentatively exclude that any of them are associated with the GW event (see Sect. 4.4 and Appendix C for more details on the rejected candidates). The limiting magnitudes of the stacked images were estimated by means of artificial star experiments. The limiting magnitude for each pointing, defined as the magnitude at which 50% of the artificial stars are recovered, is reported in Table B.2. The limits are shallower in the first epoch because of a high background due to full Moon.

<footnote n="3">³ https://goto-observatory.org</footnote>

<footnote n="4">⁴ https://www.aavso.org/apass</footnote>

<footnote n="5">⁵ https://minorplanetcenter.net/cgi-bin/checkmp.cgi</footnote>

<footnote n="6">⁶ https://wis-tns.weizmann.ac.il/</footnote>
S190814bv - Sky Localization and Coverage

<figure><img src="image_2.png" /><figcaption>**Fig. 1.** Coverage maps from the wide-field surveys as listed in Table 1 with the probability contours of the initial skymap (BAYESTAR) and the refined skymap LALInference.</figcaption></figure>

## 2.3. The search for transients with VISTA

The 4.1 m Visible and Infrared Survey Telescope for Astronomy (VISTA) is sited at Cerro Paranal and operated by the European Southern Observatory (ESO). The VISTA InfraRed CAMera (VIRCAM) has 16 detectors arranged in a sparse array, and conventionally six pointings are combined with offsets to form a contiguous "tile" of ~1.6 deg² (Sutherland et al. 2015).

Observations were made with VISTA under the VINROUGE programme at three epochs, the first over several nights post-merger (beginning at MJD 58711.17), the second around a week later, and the final epoch roughly seven weeks post-merger which was used as our primary reference template. We only observed in the Ks band (2.15 µm), to optimise our search for a

red KN component. A large majority (>90%, see Table 1) of the LALInference localisation area (referred to here as "VISTA-wide") was covered at all three epochs, as shown in Fig. 1. The single tile covering the highest likelihood region was re-imaged six times to provide deeper limits in that area (referred to as "VISTA-deep", enclosing ~21% of the sky localisation probability). Full details of the area covered, timing and representative depth reached are given in Tables 1 and B.2.

Initial processing of the data was performed using a pipeline based on the VISTA Data Flow System (VDFS; González-Fernández et al. 2018) modified for on-the-fly processing. Subsequently, the VINROUGE in-house pipeline for the Identification of GW counterparts through NIR Image
<table><caption>**Table 1.** Summary of wide-field survey coverage and typical limiting magnitudes.</caption><thead><tr><th>Telescope</th><th>Start MJD</th><th>Time after GW</th><th>Probability coverage</th><th>Limiting mag</th><th>Filter</th></tr></thead><tbody><tr><td>ATLAS</td><td>58709.52</td><td>-8.7 h</td><td>99.8%</td><td>18.0</td><td>c</td></tr><tr><td>GOTO</td><td>58710.09</td><td>+5.0 h</td><td>89.6%</td><td>18.7</td><td>L</td></tr><tr><td>VST</td><td>58710.37</td><td>+11.5 h</td><td>60.7%</td><td>20.9</td><td>r</td></tr><tr><td>Pan-STARRS1</td><td>58710.528</td><td>+15.50 h</td><td>89.4%</td><td>20.6, 20.3</td><td>i<sub>P1</sub>, z<sub>P1</sub></td></tr><tr><td>ATLAS</td><td>58710.60</td><td>+17.23 h</td><td>99.8%</td><td>18.0</td><td>o</td></tr><tr><td>GOTO</td><td>58711.09</td><td>+1.2 d</td><td>94.1%</td><td>18.1</td><td>L</td></tr><tr><td>VISTA-wide</td><td>58711.17</td><td>+1.3–3.4 d</td><td>94%</td><td>21.0</td><td>K<sub>s</sub></td></tr><tr><td>VISTA-deep</td><td>58711.24</td><td>+1.4 d</td><td>21%</td><td>22.0</td><td>K<sub>s</sub></td></tr><tr><td>VST</td><td>58711.2</td><td>+1.5 d</td><td>71.5%</td><td>21.9</td><td>r</td></tr><tr><td>ATLAS</td><td>58711.5</td><td>+1.6 d</td><td>99.8%</td><td>17.6</td><td>o</td></tr><tr><td>Pan-STARRS1</td><td>58713.5</td><td>+3.6 d</td><td>70.4%</td><td>21.9</td><td>z<sub>P1</sub></td></tr><tr><td>VST</td><td>58714.2</td><td>+4.3 d</td><td>87.7%</td><td>21.7</td><td>r</td></tr><tr><td>Pan-STARRS1</td><td>58716.5</td><td>+6.6 d</td><td>70.7%</td><td>23.0</td><td>z<sub>P1</sub></td></tr><tr><td>VST</td><td>58717.1</td><td>+7.2 d</td><td>87.7%</td><td>21.8</td><td>r</td></tr><tr><td>VISTA-wide</td><td>58719.05</td><td>+9.2–10.5 d</td><td>94%</td><td>21.2</td><td>K<sub>s</sub></td></tr><tr><td>VISTA-deep</td><td>58720.15</td><td>+10.3 d</td><td>21%</td><td>22.0</td><td>K<sub>s</sub></td></tr><tr><td>VST</td><td>58724.4</td><td>+14.5 d</td><td>87.7%</td><td>22.0</td><td>r</td></tr><tr><td>VISTA-wide</td><td>58750.1</td><td>+40–41 d</td><td>94%</td><td>21.0</td><td>K<sub>s</sub></td></tr><tr><td>VISTA-deep</td><td>58751.1</td><td>+41 d</td><td>21%</td><td>22.0</td><td>K<sub>s</sub></td></tr></tbody></table>

**Notes.** The start MJD refers to the start of observations on that night (for reference, the GW trigger occurred at MJD 58709.882). The given limiting magnitude is the median magnitude of the individual tiles that covered the probability listed. All times are in the observer frame.

Subtraction (IGNIS) was used to aid the search for transient sources. Using object lists generated from the VDFS pipeline for both the science and template images, positions were cross-checked to create a list in which the majority of objects visible across multiple epochs were removed, along with objects associated with error flags.

Template and science images were paired based on area coverage, with templates resized to encompass the entire science area per image. Coordinates were aligned through a combination of the astrometry.net software solve-field, and SExtractor (Bertin & Arnouts 1996) and swarp to match template image positions directly to their corresponding science images. Template images were then subtracted from science images using the hotpants tool<fnref n="7" /> (Becker 2015).

We searched for potential counterparts in the subtracted images using two approaches: first by eyeballing the regions around obvious galaxies, particularly those thought to be in the distance range of interest, and secondly through an automated search for sources.

Candidate transients from the automated search were culled based on various criteria, in particular, objects within low confidence regions (e.g. tile edges), with <5 σ detection significance, coincident with foreground stars or the bright cores of galaxies (for which the subtractions often left scars), or unusually sharp images suggestive of hot pixels rather than stellar sources. Moving sources were identified by reference to the MPC. A final check involved human vetting of remaining candidates (typically 1–10 per science image).

This process was repeated for all available science data (including the VISTA-deep field) across the three epochs, cross-matching them over with as many template files as could be attributed to each. Of the sources found in the automated procedure, all were deemed to be image artefacts. Similarly, no convincing new sources were found in the eyeball search, with the

exception of the known transient AT2019noq, which was found to have AB magnitudes $K = 20.12 \pm 0.07$ (at MJD 58711.23) and $K = 20.06 \pm 0.07$ (at MJD 58719.25). This source was marginally below the adopted significance threshold in the subtracted image.

## 2.4. The search for transients with ATLAS

The Asteroid Terrestrial-impact Last Alert System (ATLAS) is a high cadence, near-Earth asteroid (NEA) survey with two telescopes located on two separate sites in Hawaii (Mauna Loa and Haleakala). The f/2 telescopes are 0.5 m aperture with 10.56k × 10.56k pixel CCD cameras (Tonry et al. 2018a). The plate scale is 1′.86 pix⁻¹, giving each camera a FoV of 29.2 deg². Both units scan the sky between -40° < δ < +80° with a cadence of approximately two days, weather permitting. ATLAS survey mode uses two composite filters – “cyan” and “orange” (c and o, respectively). Cyan covers the Sloan g and r filters and orange covers the Sloan r and i filters.

A typical NEA survey observing cycle is comprised of a sequence of 4 slightly dithered exposures (which we call quads), each lasting 30 s, with overheads and processing requiring an additional 10 s. The 4 exposures are typically separated by 15 minutes within a 1 h period to allow for detection and linking of fast-moving objects. ATLAS frequently adjusts this NEA optimised schedule to carry out similar sequences of quads over the sky area of a GW sky map (e.g. Stalder et al. 2017). Observations are processed by an automatic pipeline to produce de-trended, sky-flattened images. These are corrected astrometrically with respect to the ICRS using Gaia stellar positions, and corrected photometrically with respect to a custom built reference catalogue (Refcat2; Tonry et al. 2018b). Difference imaging is employed to identify transients in the survey data and source extraction and measurement are carried out as described in Tonry et al. (2018a). All detections with $S/N \ge 5$ are read into a database at Queen's University Belfast and we require 3 or more detections at $S/N \ge 5$ to form an object detection. After

<footnote n="7">⁷ http://www.astro.washington.edu/users/becker/v2.0/hotpants.htm</footnote>
<table><caption>**Table 2. Transients detected with the VST (r-band).**</caption><thead><tr><th>Name</th><th>TNS name</th><th>RA (hms)</th><th>Dec (dms)</th><th>MJD</th><th>Mag (err)</th><th>Note</th><th>Tiling</th></tr></thead><tbody><tr><td>VST J005109.17-221740.7</td><td>AT2019qbu</td><td>00:51:09.173</td><td>-22:17:40.69</td><td>58715.16</td><td>21.19 ± 0.09</td><td>1</td><td>T13</td></tr><tr><td>VST J004414.33-250744.3</td><td>AT2019qby</td><td>00:44:14.334</td><td>-25:07:44.32</td><td>58711.27</td><td>21.35 ± 0.02</td><td>1</td><td>T1</td></tr><tr><td>VST J005653.99-275921.4</td><td>AT2019qbz</td><td>00:56:53.987</td><td>-27:59:21.37</td><td>58711.30</td><td>20.89 ± 0.07</td><td>1</td><td>T24</td></tr><tr><td>VST J004548.54-264939.0</td><td>AT2019qca</td><td>00:45:48.540</td><td>-26:49:39.01</td><td>58725.33</td><td>21.51 ± 0.05</td><td>1</td><td>T9</td></tr><tr><td>VST J004619.06-260843.2</td><td>AT2019qcb</td><td>00:46:19.062</td><td>-26:08:43.19</td><td>58714.24</td><td>21.55 ± 0.04</td><td>1</td><td>T8</td></tr><tr><td>VST J005349.82-244549.6</td><td>AT2019qcc</td><td>00:53:49.820</td><td>-24:45:49.58</td><td>58725.35</td><td>22.14 ± 0.09</td><td>1</td><td>T21</td></tr><tr><td>VST J004656.70-252236.7</td><td>AT2019npd</td><td>00:46:56.711</td><td>-25:22:36.43</td><td>58717.22</td><td>19.87 ± 0.06</td><td>2</td><td>T8</td></tr><tr><td>VST J004847.88-251823.5</td><td>AT2019noq</td><td>00:48:47.882</td><td>-25:18:23.46</td><td>58711.23</td><td>19.96 ± 0.03</td><td>2</td><td>T16</td></tr><tr><td>VST J005605.55-243826.4</td><td>AT2019nve</td><td>00:56:05.510</td><td>-24:38:26.40</td><td>58710.38</td><td>20.53 ± 0.08</td><td>2</td><td>T21</td></tr><tr><td>VST J004659.45-230559.5</td><td>AT2019nyv</td><td>00:46:59.451</td><td>-23:05:59.50</td><td>58711.27</td><td>21.20 ± 0.09</td><td>2</td><td>T12</td></tr><tr><td>VST J005542.30-244149.9</td><td>AT2019nvd</td><td>00:55:42.301</td><td>-24:41:49.93</td><td>58725.35</td><td>21.09 ± 0.05</td><td>2</td><td>T21</td></tr><tr><td>VST J005002.82-224118.8</td><td>AT2019mwp</td><td>00:50:02.820</td><td>-22:41:18.78</td><td>58711.26</td><td>20.57 ± 0.13</td><td>2</td><td>T13</td></tr><tr><td>VST J004804.40-234750.9</td><td>AT2019ntm</td><td>00:48:04.398</td><td>-23:47:50.94</td><td>58711.25</td><td>21.07 ± 0.11</td><td>2</td><td>T10</td></tr><tr><td>VST J010001.84-264251.3</td><td>AT2019ntr</td><td>01:00:01.843</td><td>-26:42:51.32</td><td>58711.29</td><td>21.21 ± 0.15</td><td>2</td><td>T27</td></tr><tr><td>VST J005012.07-261152.6</td><td>AT2019ntp</td><td>00:50:12.072</td><td>-26:11:52.56</td><td>58711.23</td><td>21.04 ± 0.03</td><td>2</td><td>T16</td></tr><tr><td>VST J005305.56-242138.7</td><td>AT2019npz</td><td>00:53:05.560</td><td>-24:21:38.71</td><td>58714.25</td><td>20.87 ± 0.08</td><td>2</td><td>T15</td></tr><tr><td>VST J004616.81-242221.2</td><td>AT2019nxe</td><td>00:46:16.814</td><td>-24:22:21.19</td><td>58714.25</td><td>20.87 ± 0.07</td><td>2</td><td>T7</td></tr><tr><td>VST J004320.49-255302.1</td><td>SN2019mbq</td><td>00:43:20.493</td><td>-25:53:02.07</td><td>58715.17</td><td>18.83 ± 0.04</td><td>2</td><td>T2</td></tr><tr><td>VST J004901.74-231404.9</td><td>AT2019nuj</td><td>00:49:01.738</td><td>-23:14:04.93</td><td>58715.16</td><td>21.72 ± 0.18</td><td>2</td><td>T13</td></tr><tr><td>VST J004133.33-234432.0</td><td>AT2019npe</td><td>00:41:33.330</td><td>-23:44:31.95</td><td>58717.15</td><td>21.51 ± 0.07</td><td>2</td><td>T0</td></tr><tr><td>VST J005552.40-254659.8</td><td>AT2019npw</td><td>00:55:52.399</td><td>-25:46:59.81</td><td>58711.24</td><td>21.34 ± 0.05</td><td>2</td><td>T23</td></tr><tr><td>VST J004330.16-224329.4</td><td>AT2019nsm</td><td>00:43:30.160</td><td>-22:43:29.35</td><td>58717.15</td><td>21.39 ± 0.09</td><td>2</td><td>T5</td></tr><tr><td>VST J005531.60-225808.5</td><td>AT2019num</td><td>00:55:31.602</td><td>-22:58:08.48</td><td>58717.14</td><td>21.39 ± 0.09</td><td>2</td><td>T19</td></tr><tr><td>VST J005243.34-233753.6</td><td>AT2019nva</td><td>00:52:43.339</td><td>-23:37:53.64</td><td>58717.13</td><td>21.36 ± 0.07</td><td>2</td><td>T14</td></tr><tr><td>VST J005646.69-250933.3</td><td>AT2019nqw</td><td>00:56:46.693</td><td>-25:09:33.29</td><td>58725.35</td><td>20.77 ± 0.03</td><td>2</td><td>T21</td></tr><tr><td>VST J005806.46-245014.3</td><td>AT2019nzd</td><td>00:58:06.456</td><td>-24:50:14.28</td><td>58714.23</td><td>21.18 ± 0.10</td><td>3</td><td>T25</td></tr><tr><td>VST J005756.90-243400.5</td><td>AT2019nys</td><td>00:57:56.904</td><td>-24:34:00.48</td><td>58714.23</td><td>21.31 ± 0.10</td><td>3</td><td>T25</td></tr><tr><td>VST J005332.32-234958.5</td><td>SN2019npv</td><td>00:53:32.316</td><td>-23:49:58.50</td><td>58717.14</td><td>21.62 ± 0.09</td><td>3</td><td>T20</td></tr></tbody></table>

**Notes.** 1 – New transient candidates which we first reported in TNS. We verified that AT2019qcb and AT2019qcc are also visible in Pan-STARRS1 subtractions. AT2019qcb (VST J004619.06-260843.2) is detected at $z_{P1} = 21.1 \pm 0.1$ mag, on MJD 58716.53. It is a nuclear transient coincident with a compact galaxy (Kron mag) $r_{P1} = 18.68$ mag. AT2019qcc (VST J005349.82-244549.6) is detected in two Pan-STARRS1 subtractions at $z_{P1} = 21.6 \pm 0.2$ mag on both MJD 58713.54 and MJD 58716.54 and has a flat light curve in $z_{P1}$. It is also coincident with a probable compact galaxy $r_{P1} = 20.7$ (Kron mag). 2 – Independent discoveries already reported in TNS by other groups. 3 – Candidates that are detected on our images but are below the threshold of our search criteria.

such objects are defined, they are subject to various quality filters, machine learning algorithms and cross-matching to known minor planet, star and galaxy catalogues.

ATLAS was serendipitously observing the S190814bv skymap region several hours before the GW detection during its normal survey mode. Hence any recent, young and bright transients would have been identified. Seven pointings of ATLAS covered the entirety of the LALInference skymap, and the first pre-discovery observation of the map started at MJD 58709.52 (8.7 h before S190814bv). The coverage continued until MJD 58709.635, or 2.8 h later (Fig. 1). Only the Haleakala telescope observed, in the cyan (c) filter. Some of the earliest exposures were affected by cloud cover and moonlight. No new transient objects which are not cross-matched with stars or known AGN were found in our images. ATLAS re-observed the field on the next two subsequent nights, in the o band (with the Mauna Loa unit). The second night of observations began at MJD 58710.602 (17.26 h after S190814bv) covering 99.8% of the localisation area within a 1 h period. The third night of observations began at MJD 58711.6, again covering 99.8% of the LALInference skymap probability. In none of the three post-event epochs did we find any new transients within the GW localisation area of ATLAS.

## 2.5. The search for transients with Pan-STARRS1

The Pan-STARRS system (Chambers et al. 2016) comprises 2 × 1.8 m telescopes on Haleakala, each with a 1.4-Gigapixel camera mounted at the Cassegrain f/4.4 focus of each unit. Here we describe observations with the Pan-STARRS1 telescope (PS1) and the camera GPC1. The GPC1 is composed of sixty Orthogonal Transfer Array devices (OTAs), each of which has a detector area of 4846 × 4868 pixels. The 10 micron pixels (0.26) give a focal plane of 418.88 mm in diameter or 3.0 degrees. This provides a FoV area of 7.06 deg², and an active region of about 5 deg² (see Chambers et al. 2016, for a description of the focal plane gaps). The five filter system (generally denoted grizyP1) is described in Tonry et al. (2012) and Chambers et al. (2016). For filters in common, the PS1 filters have similar transmission profiles as those from the Sloan Digital Sky Survey (SDSS; Abazajian et al. 2009). Images from Pan-STARRS1 are processed immediately with the Image Processing Pipeline (Magnier et al. 2016; Waters et al. 2016). The existence of the PS1 3π Survey data (Chambers et al. 2016) provides a ready-made template image of the whole sky north of $\delta = -30^\circ$, and we furthermore have proprietary $i_{P1}$ data in a band between $-40^\circ < \delta < -30^\circ$, giving reference sky images in the $i_{P1}$ band down to this lower declination limit. All individual PS1 images have a reference sky subtracted from them
<table><caption>**Table 3.** All PS1 objects discovered for S190814bv.</caption><thead><tr><th>PS1 ID</th><th>RA J2000</th><th>Dec J2000</th><th>Classification type</th><th>Disc. epoch</th><th>Disc. mag</th><th>Host</th><th>Redshift</th><th>IAU ID</th><th>Prob. contour</th></tr></thead><tbody><tr><td>PS19ekf</td><td>00:46:57.39</td><td>-24:21:42.6</td><td>Probable SN</td><td>58710.547</td><td>19.66 (i)</td><td>WISEA J004657.40-242142.6</td><td>-</td><td>AT2019nbp</td><td>40</td></tr><tr><td>PS19epf</td><td>00:48:48.77</td><td>-25:18:23.4</td><td>Probable SN</td><td>58710.585</td><td>19.93 (i)</td><td>WISEA J004847.51-251823.0</td><td>-</td><td>AT2019noq</td><td>20</td></tr><tr><td>PS19eph</td><td>00:49:51.99</td><td>-24:16:17.7</td><td>Probable SN</td><td>58710.545</td><td>19.46 (i)</td><td>6dFJ0049520-241618</td><td>0.435622</td><td>AT2019nor</td><td>10</td></tr><tr><td>PS19epw<sup>(1)</sup></td><td>00:46:56.71</td><td>-25:22:36.6</td><td>VS in NGC 253</td><td>58710.586</td><td>20.28 (i)</td><td>NGC 253</td><td>0.0008</td><td>AT2019npd</td><td>50</td></tr><tr><td>PS19epx</td><td>00:56:50.42</td><td>-24:20:50.0</td><td>Probable SN</td><td>58710.587</td><td>20.66 (z)</td><td>WISEA J005650.42-242050.3</td><td>-</td><td>AT2019nqp</td><td>80</td></tr><tr><td>PS19epz</td><td>00:50:26.34</td><td>-25:52:57.8</td><td>Probable SN</td><td>58713.541</td><td>21.85 (z)</td><td>faint, uncatalogued host</td><td>-</td><td>AT2019nuw</td><td>20</td></tr><tr><td>PS19eqa</td><td>00:50:21.01</td><td>-23:42:46.7</td><td>Probable SN</td><td>58713.541</td><td>21.75 (z)</td><td>WISEA J005021.03-234246.0</td><td>-</td><td>AT2019nux</td><td>50</td></tr><tr><td>PS19eqb</td><td>00:50:50.39</td><td>-25:29:29.5</td><td>Probable SN</td><td>58713.541</td><td>21.01 (z)</td><td>PSO J012.7099-25.4915</td><td>-</td><td>AT2019nuy</td><td>20</td></tr><tr><td>PS19eqc</td><td>00:49:52.26</td><td>-25:31:25.6</td><td>Probable SN</td><td>58713.541</td><td>21.89 (z)</td><td>PSO J012.4678-25.5238</td><td>-</td><td>AT2019nuz</td><td>20</td></tr><tr><td>PS19eqd</td><td>00:52:43.39</td><td>-23:37:54.0</td><td>Probable SN</td><td>58713.541</td><td>21.49 (z)</td><td>PSO J013.1807-23.6317</td><td>-</td><td>AT2019nva</td><td>70</td></tr><tr><td>PS19eqe</td><td>00:46:51.16</td><td>-25:25:39.3</td><td>VS in NGC 253</td><td>58713.541</td><td>21.72 (z)</td><td>NGC 253</td><td>-</td><td>AT2019nvb</td><td>50</td></tr><tr><td>PS19eqf<sup>(2)</sup></td><td>00:52:18.32</td><td>-26:19:42.0</td><td>SN II</td><td>58713.543</td><td>21.31 (z)</td><td>WISEA J005218.36-261942.5</td><td>0.070</td><td>AT2019nvc</td><td>50</td></tr><tr><td>PS19eqg</td><td>00:55:42.39</td><td>-24:41:50.2</td><td>Probable SN</td><td>58713.541</td><td>21.47 (z)</td><td>PSO J013.9262-24.6973</td><td>-</td><td>AT2019nvd</td><td>60</td></tr><tr><td>PS19eqh<sup>(3)</sup></td><td>00:56:05.51</td><td>-24:38:26.3</td><td>Probable SN</td><td>58713.541</td><td>21.30 (z)</td><td>PSO J014.0230-24.6407</td><td>-</td><td>AT2019nve</td><td>60</td></tr><tr><td>PS19eqi<sup>(4)</sup></td><td>00:53:32.30</td><td>-23:49:58.6</td><td>SN Ib</td><td>58713.541</td><td>21.26 (z)</td><td>WISEA J005332.35-234955.8</td><td>0.056</td><td>SN2019npv</td><td>70</td></tr><tr><td>PS19eqj<sup>(5)</sup></td><td>00:55:52.39</td><td>-25:46:59.7</td><td>SN IIb</td><td>58713.544</td><td>21.35 (z)</td><td>PSO J013.9687-25.7831</td><td>0.163</td><td>AT2019npw</td><td>70</td></tr><tr><td>PS19eqk<sup>(6)</sup></td><td>00:56:46.71</td><td>-25:09:33.4</td><td>Probable SN</td><td>58713.542</td><td>21.23 (z)</td><td>PSO J014.1947-25.1593</td><td>-</td><td>AT2019nqw</td><td>80</td></tr><tr><td>PS19eqo</td><td>00:48:16.08</td><td>-25:28:14.9</td><td>Probable SN</td><td>58713.539</td><td>20.89 (z)</td><td>WISEA J004816.11-252814.8</td><td>-</td><td>AT2019nvr</td><td>40</td></tr><tr><td>PS19eqp</td><td>00:52:37.75</td><td>-26:11:41.4</td><td>Probable SN</td><td>58713.549</td><td>21.44 (z)</td><td>WISEA J005237.72-261142.4</td><td>-</td><td>AT2019nvs</td><td>50</td></tr><tr><td>PS19eqq<sup>(7)</sup></td><td>00:50:12.06</td><td>-26:11:52.8</td><td>SN Ic-BL</td><td>58713.541</td><td>21.31 (z)</td><td>WISEA J005012.11-261154.7</td><td>-</td><td>AT2019ntp</td><td>50</td></tr><tr><td>PS19erd</td><td>00:55:19.23</td><td>-26:11:50.7</td><td>Probable SN</td><td>58716.542</td><td>21.41 (z)</td><td>WISEA J005519.14-261150.9</td><td>-</td><td>AT2019ofb</td><td>70</td></tr></tbody></table>

**Notes.** The host galaxies are the primary names as catalogued now in NED. Spectroscopically classified events are noted, and “Probable SN” means the light curve points we have are consistent with it being an unrelated supernova. There are two objects in the nearby galaxy NGC 253, which are certain variable stars in the outskirts of the disc (labelled VS in the table). The final column gives the 2D skymap probability contour within which the transient position lies. Some of the sources are not in the redshift range of the GW event (i.e. 0.046–0.068). A machine-readable file with all photometry for these candidates is available at the CDS. (¹)Discovered by DECam-Growth: DG19hqpgc. (²)Classified by Rodriguez et al. (2019). (³)Offset by 3′′8 from galaxy WISEA J005605.37-243830.5, but coincident with faint uncatalogued stellar source. (⁴)Discovered by DECam-Growth: DG19wxnjc. Classified by Gomez et al. (2019b), De et al. (2019a), and Jonker et al. (2019). (⁵)Discovered by DECam-Growth: DG19wgmjc. Classified by Tucker et al. (2019a). (⁶)Discovered by DECam-Growth: DG19xczjc. (⁷)Discovered by DECam-Growth: DG19gcwjc. Classification reported in Wiesner et al. (2019a), but no redshift given.

and sources with at least two detections with $S/N \ge 5$ significance and spatially coincident to within 0.′5 are detected and measured. PS1 typically observes in a quad sequence similar to ATLAS, with a set of 4 × 45 s exposures taken across a time span of 1 h to identify and link moving sources. The PS1 team can intervene at any moment and direct the telescope to observe a LVC GW sky map with a flexible choice of filter, exposure time, coverage, and dither and stack strategy. The difference images can be combined into deeper stacks or processed individually and the sources resulting from these are read into a large database at Queen’s University Belfast. A series of quality control filters, machine learning algorithms and cross-matches against minor planet, stellar and galaxy catalogues are automatically run and human scanning occurs for all objects not coincident with known solar system objects, stars or catalogued AGN (see Smartt et al. 2016a,b, for more details). At the detection time of S190814bv (2019 Aug 14 21:10:39.01), Hawaii was in day time and PS1 began observing the field at 2019 Aug 15 12:40:37 UT, or 15.50 h after the LVC discovery time.

On the first night of observation the individual 45 s exposures (called “warps”) were processed individually to search for any fading transient over the 2 h 33 min period of observation. Image sensitivities are estimated by injecting 500 point sources per sky-cell across a range of magnitudes and the limiting magnitude is defined when 50% of the sources are recovered (described in the content of database table *DiffDetEffMeta* in Flewelling et al. 2016). Each chip exposure is warped onto a pre-defined tessellation (called skycells, see Chambers et al. 2016), and the limits refer to these skycells. The 45 s exposures were combined into a nightly stack on the first three nights of observing. The stacks are made by median combining the warps of each skycell. On the first observing night 25 exposures were combined in $i_{P1}$ and 31 exposures in $z_{P1}$ in each skycell stack, giving a typical exposure

time of 1125 s and 1395 s in $i_{P1}$ and $z_{P1}$, respectively. We did not find any fading transient, but the true constraints are weak due to the dither strategy and fill factor.

For the two subsequent nights, we did not process the individual images, rather we combined all the $z_{P1}$ band warps into a nightly stack. The effective exposure times were 12480 s and 13440 s on these respective nights. These were deeper than the $3\pi$ reference stacks in $z_{P1}$ in this sky region, so over the next four weeks PS1 observed the region in $z_{P1}$ to make a deeper and more uniform reference stack. The limiting magnitudes of the skycells on the three nights observing of S190814bv were calculated using the new, custom-made deeper reference-stack for template subtraction. The final sky coverage is plotted in Fig. 1.

All images were processed through the Image Processing Pipeline described above. Detections coincident with known stellar objects from the *Gaia* DR1 (Gaia Collaboration 2016), Guide Star Catalogue<fnref n="8" />, Two Micron All Sky Survey (2MASS, Skrutskie et al. 2006), SDSS DR12 (Alam et al. 2015) and PS1 (Chambers et al. 2016) catalogues were rejected. Additionally any objects coincident with known AGN were identified and excised from the transient search list. As discussed in Smartt et al. (2016a,b), the AGN identification is based mostly on the Véron-Cetty & Véron (2001) and MILLIQUAS<fnref n="9" /> catalogues (Flesch 2015). The resultant objects are spatially cross-matched against known galaxies (mostly through the NASA Extragalactic Database, NED<fnref n="10" />) and all are visually inspected. The objects discovered are listed in Table 3, along with their likely classification (Huber et al. 2019; Srivastas et al. 2019; Smartt et al. 2019a). None of these objects is a viable counterpart of S190814bv (see Sect. 4.4 for details on the candidate rejection).

<footnote n="8">⁸ https://archive.stsci.edu/gsc/</footnote>

<footnote n="9">⁹ http://quasars.org/milliquas.htm</footnote>

<footnote n="10">¹⁰ https://ned.ipac.caltech.edu/</footnote>
<figure><img src="image_3.png" /></figure>

**Fig. 2.** LVC skymap of S190814bv (LALInference.v1.fits) in greyscale, with galaxies selected through HOGWARTS over-plotted (black dots). Only the higher probability Northern region of the skymap is included, since the Southern region was not followed up by ENGRAVE. The size of the symbol of each galaxy (i.e. the black dots) is proportional to the probability of hosting the GW event (see Table B.1) given the skymap and a weighing scheme following Arcavi et al. (2017a). The various instruments are illustrated with different colours as in the figure legend, and the typical limiting magnitudes and filters used are given in Table B.1.

### 3. Galaxy targeted searches

In addition to the wide-field survey coverage, the unusually tight localisation map of S190814bv (5 deg² at 50%), and the distance estimate available from the GW signal (267 ± 52 Mpc), allowed us to define a coordinated programme of multi-wavelength observations of galaxies within the localisation region (e.g. Nissanke et al. 2013; Evans et al. 2016a; Gehrels et al. 2016). While these images cannot cover the whole 2D skymap, they can (often) significantly improve upon the depth of the wide-field surveys for a select number of high-luminosity galaxies (Fig. 2). To identify galaxies with the highest probability of hosting the GW event we utilised the HOGWARTS code<fnref n="11" /> (Salmon et al. 2020), which ranks galaxies in the Galaxy List for the Advanced Detector Era (GLADE; Dálya et al. 2018) catalogue according to their probability of containing the corresponding merger given the 3D localization probability density (Singer et al. 2016), and based on the expectation that NS–BH merger rates follow the galaxy mass distribution (Arcavi et al. 2017a).

Since the expected colours and luminosities of the counterparts of NS–BH mergers (see Sect. 4.1) still have significant uncertainties (largely due to the lack of observational constraints), our goal was to obtain multi-colour (optical and NIR) imaging, which we prioritised over observing a greater number of galaxies. This strategy enabled our observations to be sensitive to counterparts that were either blue (e.g. disc-wind driven), or very red due to high lanthanide opacities in dynamical ejecta. While our observations targeted the most luminous galaxies weighted for

<figure><img src="image_4.png" /></figure>

**Fig. 3.** VLT/HAWK-I image of a galaxy targeted field. A number of catalogued galaxies (and at least one uncatalogued galaxy likely at the same redshift) are visible in the field. The insets show each individual galaxy as well as the resulting subtraction, demonstrating the absence of variable sources to the limits of the data in any of these possible host galaxies. Each galaxy inset is labelled by the HyperLEDA identifier from the GLADE catalogue, and the corresponding limiting magnitudes are listed in Table B.1.

the localization probability, individual telescope pointings were refined in order to capture additional (lower luminosity) galaxies within the localisation volume of the LALInference skymap. We obtained a series of coordinated observations using the Gamma-Ray burst Optical and Near-IR Detector (GROND), the Liverpool Telescope (LT), the Nordic Optical Telescope (NOT), the Telescopio Nazionale Galileo (TNG), the Very Large Telescope (VLT) and the *William Herschel* Telescope (WHT). In total, over 400 multi-wavelength (grizJHK) images of the 67 most probable galaxies within the 3D volume were obtained in the ten days following the merger. When generating target lists for each telescope, we attempted to avoid unnecessary duplication of observations, while ensuring that the highest probability galaxies were observed to the greatest possible depth. In practice weather, seeing and other scheduling constraints meant that some duplication was unavoidable. Our global coverage is shown in Fig. 2 and a list of observed galaxies in order of decreasing probability (as defined in Sect. 3.7) is given in Table B.1. An example set of observations with the VLT High Acuity Wide field Ks band Imager (HAWK-I; Pirard et al. 2004; Casali et al. 2006; Kissler-Patig et al. 2008; Siebenmorgen et al. 2011) is shown in Fig. 3.

Our techniques for searching for transient objects depended on the nature of the data available. All images obtained were manually, and rapidly, compared against existing optical survey data, in particular the PS1 3π survey (e.g. as was done in Coulter et al. 2017b, when AT2017gfo was first discovered). Given the brightness and proximity of the Moon at the time of the observations, only the VLT data exceeded the depth of the PS1 3π survey. For some observations these data remained the best comparison.

In most other cases, when reference images were subsequently obtained, we performed PSF-matched image subtraction using the *hotpants* code. The residual images were then manually inspected to identify any possible transient sources. We

<footnote n="11">¹¹ https://gwtool.watcher telescope.ie/</footnote>
limited our search to the circle centred on the nucleus of the
investigated galaxy, with radius 1.5×R<sub>25</sub>, where R<sub>25</sub> is the galaxy
isophotal radius at B = 25 mag arcsec<sup>-212</sup>. In general, all galaxies
are well subtracted except for some of the brighter nuclei
that leave notable residuals and prevent the search for transients.
We confirmed that no transients are identified in the difference
images. To quantify the depth of these images we inserted artificial
stars into the images in different positions within the galaxy
search radius with a range of magnitudes and estimated their
recovery in our difference images. The limiting magnitude is
defined as the average magnitude of the faintest artificial stars
that can be visually identified, where optical and NIR photometry
is calibrated against the PS1 and 2MASS catalogues, respectively.
In all cases these stars had S/N ~ 3. We found that the
limiting magnitude is fairly constant at different positions within
the galaxy search radius with the exception of the nuclei.

**3.1. VLT observations**

The Very Large Telescope (VLT) is a facility operated by ESO
on Cerro Paranal in the Atacama Desert of northern Chile
(Arsenault et al. 2006), which consists of four individual 8.2 m
telescopes (UT1–4). We obtained observations of 16 high-
priority galaxies with the VLT using *i* and *z* imaging with
the FOcal Reducer and low dispersion Spectrograph (FORS –
Appenzeller et al. 1998, mounted at UT1), and NIR imaging
for further 17 galaxies with HAWK-I (Kissler-Patig et al. 2008,
mounted at UT4; see Fig. 3 for an example image) in the *K*<sub>s</sub>
band (see Table B.1). All VLT data were reduced using the stan-
dard EsoReflex graphical environment (v2.9.1; Freudling et al.
2013). The observations were performed in one epoch of FORS
imaging on 2019 Aug 16 and three epochs of HAWK-I imaging
on 2019 Aug 16, 22–23 and 2019 Sept 23–24. Given the likely
slow rise time of KNe in the NIR bands, the first two HAWK-I
*K*<sub>s</sub>-band epochs were intended to be sensitive to the peak of the
KN a few days after merger time. This complemented the FORS
optical imaging within the first 24–48 h, which is more sensitive
to early emission.

FORS observations consisted of 3 × 100 s observations in the *i* band, although one field was erroneously observed in the *z* band for the same exposure time. These images reached significantly deeper limiting magnitudes (*i* ~ 23–24.5 mag) across the field than those obtained by smaller aperture telescope searches. The cores of some galaxies were, on occasion, saturated, removing our ability to detect transients close to the nucleus.

**3.2. WHT observations**

A series of optical and NIR observations of 17 galaxies in the
sample were taken with the *William Herschel* Telescope (WHT,
Boksenberg 1985) from 2019 Aug 14–22. Optical observations
of 17 galaxies were obtained in the *r* band using the Auxiliary-
port CAMera (ACAM, Benn et al. 2008) instrument on 2019
Aug 15, while NIR observations of 12 galaxies were taken in
the *K*<sub>s</sub> band using the Long-slit Intermediate Resolution Infrared
Spectrograph (LIRIS, Acosta-Pulido et al. 2002) over the fol-
lowing nights. Both the LIRIS and ACAM images were reduced
using standard IRAF procedures<sup>13</sup> and the custom LIRIS pack-
age for LIRIS<sup>14</sup>.

**3.3. TNG observations**

Optical and NIR images of a subset of 19 galaxies were carried
out with the Italian 3.6 m Telescopio Nazionale Galileo (TNG,
Poretti 2018), situated on La Palma, using the optical DOLORes
(Molinari et al. 1999) and near-infrared NICS (Oliva et al. 2001)
instruments. Ten galaxies were observed in the *r* band with rela-
tively short (120 s) exposures obtained on 2019 Aug 15 between
02:13 and 02:54 UT (D'Avanzo et al. 2019a). Image reduction
was carried out using standard IRAF procedures. Astrometry was
performed using the USNO-B1.0<sup>15</sup> catalogue. The typical upper
limit is *r* ~ 22.8 (3σ detection limit).

NIR observations of five galaxies were undertaken on 2019 Aug 16, Aug 20 and Sept 5 (usually between 03:00 and 05:00 UT). In addition, the galaxy ESO474-026 was observed on 2019 Aug 17. Each galaxy was observed for 1200 s total exposure time with the *K*<sub>s</sub> filter (D'Avanzo et al. 2019b). Image reduction was carried out using the *jitter* task of the ESO-eclipse package<sup>16</sup>. Astrometry and photometric zeropoints were calculated using the 2MASS<sup>17</sup> catalogue. The typical upper limit is *K*<sub>s</sub> ~ 19.7 – 20.9 (3σ detection limit).

**3.4. GROND observations**

We observed 36 galaxies simultaneously in *g', r', i', z', J, H, Ks* with GROND (Greiner et al. 2008), mounted at the 2.2 m MPG telescope at the ESO La Silla Observatory. For each galaxy we obtained an average exposure of 2.1 min in the optical bands and 3.9 min in the NIR bands, and the data were reduced using the GROND pipeline (Krühler et al. 2008), which applies bias and flat-field corrections, stacks images and provides an astrometric calibration. The observations reached typical 3σ detection limits of 20–22 mag in the *r'* band and 17.5–19.5 mag in *Ks*.

**3.5. LT observations**

The Liverpool Telescope (LT; Steele et al. 2004) is a 2 m fully
robotic telescope on the Canary island of La Palma, Spain.
A total of 19 galaxies were observed using the IO:O imaging
camera. IO:O has a 10 × 10 arcmin² FoV and was operated with
a 2×2 binning, providing a pixel scale of 0.′/3 pix⁻¹. The observa-
tions were made between 01:38 and 05:34 UT on 2019 Aug 15.
For all fields, 2 × 150 s exposures in *r* band were obtained, and
for some of the highest-probability candidates we also obtained
2 × 150 s exposures in *i* band. Reduced images were provided by
the IO:O pipeline and stacked with SWarp<fnref n="18" />. Image subtraction
of these data were performed using our own subtraction tools
rather than hotpants, and detection limits were measured per-
forming PSF photometry at fixed positions over a grid around
the centre of each image to determine the median and standard
deviation of the sky background as measured in each aperture.
Using these measurements, with the calibration tied to the PS1
photometric standards in each field, we derived 3σ limiting mag-
nitudes of 20.3 in both bands.

**3.6. Galaxy catalogue incompleteness**

While searches for EM counterparts targeting known galaxies
within the localisation region of a GW are eminently feasible for

<footnote n="12">¹² http://leda.univ-lyon1.fr/</footnote>

<footnote n="13">¹³ https://iraf-community.github.io</footnote>

<footnote n="14">¹⁴ http://www.ing.iac.es/Astronomy/instruments/liris/liris_ql.html</footnote>

<footnote n="15">¹⁵ https://www.usno.navy.mil/USNO/astrometry/optical-IR-prod/icas/fchpix</footnote>

<footnote n="16">¹⁶ https://www.eso.org/sci/software/eclipse/</footnote>

<footnote n="17">¹⁷ https://irsa.ipac.caltech.edu/Missions/2mass.html</footnote>

<footnote n="18">¹⁸ https://www.astromatic.net/software/swarp</footnote>
<figure><img src="image_5.png" /></figure>

**Fig. 4.** 50, 95, and 99% probability regions for S190814bv. Boxes indicate the northern and southern regions of the map as discussed in the text. Galaxies which have a spectroscopic redshift in NED, and lie within the 95% contour at a distance ±3σ that of S190814bv are marked in black; galaxies with an inconsistent spectroscopic redshift are plotted in yellow. The inhomogeneous coverage of NED in the northern contour is clearly visible.

nearby events, such an approach becomes less effective as dis-
tance increases. There are two reasons for this: firstly, the density
of galaxies per unit area on the sky increases such that tiling the
GW map becomes more efficient; secondly, the completeness of
galaxy catalogues drops off precipitously beyond 200–300 Mpc.
Nonetheless, for S190814bv we used targeted deep optical and
NIR observations of some of the most likely host galaxies in
the probability map as reported by the HOGWARTS ranking tool.
We thus needed to determine the completeness of the galaxy
catalogues that HOGWARTS used for the position and distance of
S190814bv.

To assess the completeness of GLADE, we queried NED for any galaxies within the 95% probability region of the LALInference skymap, and with a listed spectroscopic redshift. This resulted in 5209 galaxies, of which 1376 have a spectroscopic redshift within 3σ of the S190814bv distance luminosity marginalised over the whole sky (267 ± 52 Mpc). We plot the positions of these galaxies in Fig. 4. What is apparent from Fig. 4 is that the completeness of the NED database varies across the map, with a sharp drop in the number of galaxies above Dec = −25°. This is almost certainly due to the lack of coverage of the 2dF galaxy redshift survey above Dec = −25° (Colless et al. 2001).

We attempt to quantify this varying incompleteness in NED
in order to determine what fraction of stellar luminosity and
mass we have covered in the galaxy targeted search. To this end,
we selected galaxy candidates from the PS1 3π catalogue, as
this is the deepest, most homogeneous public imaging catalogue
available over the whole skymap. PS1 reaches a limiting magni-
tude of 98% completeness for point sources of 22.5–23 mag in
each of g<sub>P1</sub>, r<sub>P1</sub>, i<sub>P1</sub>, with extended source completeness being
about 0.5 mag brighter (Chambers et al. 2016). For reference, an
apparent magnitude limit of ~22 mag corresponds to an absolute
magnitude ~−15.1 mag at a distance of 267 Mpc. This absolute
magnitude is comparable to that of the Small Magellanic Cloud,
and so at the distance of S190814bv PS1 is essentially complete
to all galaxies of relevance. While we may miss some very low
surface brightness dwarf galaxies, these contain so little stellar
mass that they can be ignored for our purposes (see Sect. 3.7 for
a discussion of this).

To create our galaxy candidate catalogue for S190814bv, we
queried the PS1 database (Flewelling et al. 2016) for all sources

<figure><img src="image_6.png" /></figure>

**Fig. 5.** Example 15" × 15" Pan-STARRS gri cutouts around extended sources identified by our cuts.

within the northern 95% localisation region. In order to select
only extended objects, we require that gPSF − gKron > 0.1 mag,
and rPSF − rKron > 0.1 mag, and in addition that the source has
n detections > 10 within the PS1 catalogue. Finally, we limit our-
selves to the brightest galaxies in the field, setting a threshold
of rKron < 20 mag, which is equivalent to an absolute magnitude
of ~−17 mag at the distance of S190814bv. We also mask out
regions in our catalogue around the Sculptor Galaxy NGC 253
and globular cluster NGC 288, which both contain a large num-
ber of spurious detections in the PS1 catalogue.

Visual inspection of a random sample of sources from our
extended source catalogue confirms that the majority (≥90%) are
indeed galaxies (Fig. 5). The small number of sources brighter
than *r* = 14 mag in the catalogue all appear to be saturated,
bright stars rather than galaxies, so we impose a brightness cut-
off at *r* = 14 mag. We are finally left with 23 466 candidate
galaxies within the northern 95% localisation probability region
of S190814bv. In order to better assess the issue of complete-
ness, we cross-matched our PS1 galaxy catalogue against the
NED and GLADE galaxy lists, requiring a matching radius of
<1′.5. We show the fraction of galaxies that have an associated
NED or GLADE counterpart in Fig. 6. The GLADE complete-
ness reaches a maximum of ~80% for galaxies brighter than
*r* ~ 16.5 mag, but drops rapidly at fainter magnitudes, with
a completeness of only ~50% at *r* = 17.5 mag, and ≤20%
for *r* < 18 mag; the NED completeness is substantially lower
between *r* ~ 15.5–17.5 mag.

The completeness of galaxy catalogues in the context of
gravitational wave searches was also recently considered by
Kulkarni et al. (2018). These authors employed a different
approach to this work, by using nearby supernovae with known
distances to generate a random sample of galaxies. From this,
Kulkarni et al. then determine the fraction with extant spectro-
scopic redshifts. While Kulkarni et al. look at a somewhat closer
distance (<200 Mpc), it is nonetheless encouraging that their
“Relative Completeness Fraction” of 75% (which is implicitly
weighted by host mass), is comparable to our completeness for
the most luminous galaxies.

3.7. Probability covered by our targeted search

As demonstrated in the previous section, while the GLADE cat-
alogue is somewhat incomplete in terms of galaxy number in
<figure><img src="image_7.png" /><figcaption>Fig. 6. Orange and green lines, left axis: percentage of sources in the PS1 catalogue that have an associated NED and GLADE, respectively, cross-matched galaxy with spectroscopic redshift as a function of magnitude. Blue line, right axis: histogram of galaxy counts as a function of magnitude, from our Pan-STARRS-derived extended source catalogue.</figcaption></figure>

the localisation volume, it contains the majority of the most luminous (and hence most massive) galaxies. In order to quantitatively estimate the efficiency of our galaxy-targeted search, in terms of covering the GW localisation probability, we proceed here to assign a definite probability of being the actual host of S190814bv to all catalogued galaxies in the volume, accounting for the mentioned incompleteness. The full list of targeted galaxies and the corresponding observations are reported in Table B.1. Let us consider galaxies as point-like objects, and let $i$ be an index that runs on all galaxies that are located within the volume $V_α$ that contains a given fraction $α$ of the GW 3D posterior localisation probability density $P_{3D}(RA, Dec, d_L)$ (i.e. the 3D skymap). We assume the probability $P_{gal,i}$ that the merger has taken place within galaxy $i$ to be proportional to the product between $P_{3D}(\mathbf{x}_i)$, namely the GW localisation probability density per unit volume at the galaxy position $\mathbf{x}_i = (RA_i, Dec_i, d_{L,i})$, and $R_{NS-BH,i}$, that is the NS–BH merger rate in galaxy $i$ (which in principle depends on its present properties and on its history). Artale et al. (2019) have shown, combining state-of-the-art compact binary population synthesis models and cosmological simulations, that the NS–BH rate in galaxies at low redshift correlates almost linearly with the galaxy total stellar mass ($R_{NS-BH} \propto M^{0.8}$), with some scatter driven by differences in galaxy merger histories, specific star formation rate and metallicity evolution. Based on these results, for simplicity we assume $R_{NS-BH,i} \propto M_i$ and we use the galaxy $K_s$-band luminosity $L_K$ as a proxy for galaxy mass, so that $R_{NS-BH,i} \propto L_{K,i}$. This leads to

$$P_{gal,i} = A P_{3D}(\mathbf{x}_i) L_{K,i}, \qquad (1)$$

which is similar to the galaxy ranking score used by Arcavi et al. (2017b), but with B-band replaced by $K_s$-band luminosity (which is a better tracer of galaxy mass). In order to compute the normalisation constant $A$, we impose the condition

$$\alpha = \sum_{i=1}^{N} P_{gal,i} = A \sum_{i=1}^{N} P_{3D}(\mathbf{x}_i) L_{K,i}, \qquad (2)$$

where $N$ is the total number of galaxies in the volume $V_α$ (equality (2) is strictly valid only when $α = 1$, but in practice it remains correct to an excellent approximation for $α$ close to one). Since our catalogue only contains a fraction of the actual galaxies in

the volume, we need to split the sum on the RHS of Eq. (2) into two terms

$$\sum_{i=1}^{N} P_{3D}(\mathbf{x}_i) L_{K,i} = \sum_{i=1}^{N_{\text{cat}}} P_{3D}(\mathbf{x}_i) L_{K,i} + \sum_{i=N_{\text{cat}}+1}^{N} P_{3D}(\mathbf{x}_i) L_{K,i}, \qquad (3)$$

where $N_{\text{cat}}$ is the number of GLADE galaxies within $V_α$. Assuming the remaining $K_s$-band luminosity (present in the volume, but missing from the catalogue) to be uniformly distributed within the volume, we can approximate the last term as

$$\sum_{i=N_{\text{cat}}+1}^{N} P_{3D}(\mathbf{x}_i) L_{K,i} \sim \langle P_{3D} \rangle_{V_\alpha} (L_{\text{TOT}} - L_{\text{cat}}), \qquad (4)$$

where $L_{\text{cat}}$ is the total $K_s$-band luminosity in GLADE galaxies, $L_{\text{TOT}}$ is the total $K_s$-band luminosity in the localisation volume, and $\langle P_{3D} \rangle_{V_α}$ is the GW 3D localisation probability density averaged over the volume, that is $\langle P_{3D} \rangle_{V_α} = α/|V_α|$, where $|V_α|$ is the extent of the localisation volume (e.g. in comoving Mpc$^3$).

This finally gives the normalisation constant as

$$A = \frac{\alpha}{\sum_{i=1}^{N_{\text{cat}}} P_{3D}(\mathbf{x}_i) L_{K,i} + \alpha(L_{\text{TOT}} - L_{\text{cat}})/|V_\alpha|} \qquad (5)$$

The total $K_s$-band luminosity in the localisation volume can be estimated as $L_{\text{TOT}} \sim |V_α| \times j$, where $j \sim (7 \pm 1.5) \times 10^8 L_\odot h \text{Mpc}^{-3}$ is the $K_s$-band local luminosity density (Hill et al. 2010).

In order to construct the volume that contains 95% of the localisation probability from the latest public 3D GW localisation probability density based on GW parameter estimation (the LALInference skymap; The LIGO Scientific Collaboration and the Virgo Collaboration 2019b), we employ a 3D greedy binning approach. The 3D bins are defined by dividing the sky into tiles using a healpix grid (with $N_{side} = 1024$), and further dividing the distance coordinate into 3000 linearly-spaced bins between 0 and 700 Mpc. The probability contained in each bin is assigned based on the 3D skymap (following Singer et al. 2016). The bins are then summed in order of decreasing probability density until the enclosed probability equals 95%, which defines the desired localisation volume. The extent of the obtained volume is $|V_{95\%}| \approx 1.4 \times 10^5 \text{Mpc}^3$, which gives $L_{\text{TOT}} \sim (6.9 \pm 1.5) \times 10^{13} L_\odot$. 1061 GLADE galaxies fall within this volume. Only 45% of these have a $K_s$-band measurement reported in the catalogue, due to the 2MASS magnitude limit. To circumvent this problem we utilise our own VINROUGE observations to obtain $K_s$-band magnitudes for a large fraction of the galaxies. These data were processed through the VISTA Data Flow System (González-Fernández et al. 2018), which provides outputs in the same form as for other VISTA public surveys, including catalogue counts and photometric calibration (per tile) for each observation. We therefore determine and extract $K_s$-band magnitudes for sources on each tile and cross match the resulting catalogues with our GLADE output. This results in 876 matches, providing an 82% completeness. For the brightest galaxies we use the $K_s-J$ colour from 2MASS and the redshifts reported in GLADE to k-correct the VISTA magnitudes, in order to compute the corresponding luminosities. For the remaining galaxies, we use the median k-correction (which amounts to 0.10 mag). We finally compute the $K_s$-band luminosity as

$$\log(L_K/L_\odot) = 0.4(3.27 - \tilde{K}) - \log(1+z) + 2\log(d_L/d_0), \qquad (6)$$

where 3.27 is the absolute $K_s$-band magnitude of the Sun (Willmer 2018) (in the Vega system), $\tilde{K}$ is the k-corrected
<figure><img src="image_8.png" /><figcaption>**Fig. 7.** Cumulative $K_s$-band luminosity distribution of GLADE galaxies in the S190814bv localisation region, and of our targeted galaxies, compared to the expected distribution in the 95% localisation volume assuming the $K_s$-band luminosity density and distribution from Hill et al. (2010). The grey area in the latter distribution represents the $1\sigma$ uncertainty.</figcaption></figure>

$K_s$-band magnitude of the galaxy, and $d_0 = 10$ pc. Summing over all galaxies with a $K_s$-band magnitude measurement, we obtain $L_{\text{cat}} \sim 7.8 \times 10^{13} L_{\odot}$, that is, the GLADE catalogue for S190814bv is essentially complete in terms of $K_s$-band luminosity. This can be seen in Fig. 7, which compares the cumulative $K_s$-band luminosity distribution of GLADE galaxies in the S190814bv 95% localisation volume (red line) and that of our targeted galaxies (blue line) with the expected distribution in the same volume, based on a Schechter fit to the local galaxy luminosity distribution (Hill et al. 2010). This comparison indicates that, despite the incompleteness of current catalogues, galaxy-targeting-based searches are still viable out to these distances, as already suggested, for example, by Hanna et al. (2014), Evans et al. (2016b) and Gehrels et al. (2016).

All the galaxies in our targeted search (see Table B.1) apart from five have a measured $K_s$-band magnitude reported in GLADE (from 2MASS). In other words, even though the HOGWARTS code selects the galaxies based on their *B*-band luminosity, the resulting sample is generally bright in $K_s$-band as well. We compute $L_K$ (following Eq. (6)) and therefore $P_{\text{gal},i}$ (Eq. (1)) using our VISTA magnitudes, as explained above. The resulting distribution of covered probability as a function of limiting magnitude in different bands is shown in Fig. 8. The sum of the probabilities over the targeted galaxies in our search amounts to ~50%. This does not enable us to place stringent limits on the properties of the putative EM counterpart of S190814bv using the galaxy-targeted search alone, but it nevertheless shows that targeted searches still have a reasonable chance of detecting a counterpart at ~250 Mpc.

# 4. Discussion

NS-BHs are hybrid merger events that offer insights into a range of behaviours that are not accessible through other mergers. They have both a larger total mass and a larger chirp mass than NS-NS systems. Thus they should produce a stronger GW signal that can be observed out to greater distances. No extant NS-BH systems are known and their range of masses, and astrophysical rates still therefore have very few observational constraints. Population synthesis models show NS-BH systems may be some-

<figure><img src="image_9.png" /><figcaption>**Fig. 8.** Probability (as defined in Eq. (1)) covered by our galaxy-targeted search as a function of limiting absolute magnitude in different bands.</figcaption></figure>

what rarer than NS-NS (e.g., Abadie et al. 2010; Dominik et al. 2015; Eldridge & Stanway 2016; Kruckow et al. 2018; Giacobbo & Mapelli 2018; Neijssel et al. 2019), but with significant uncertainties (Belczynski et al. 2016). There is also tentative evidence that they may contribute to the known population of cosmological short GRBs (Gompertz et al. 2020b).

During the last phase of the NS-BH coalescence, the NS can be partially or totally disrupted by the BH tidal field or swallowed directly by the BH without any significant mass left outside the merger remnant (Lattimer & Schramm 1976; Shibata & Taniguchi 2011). EM emission is expected when the disruption occurs before the NS reaches the innermost stable circular orbit (ISCO) of the BH. Tidal disruption depends on the mass ratio of the two compact objects, on the BH spin, and on the NS Equation of State (EoS, Shibata & Taniguchi 2011; Kyutoku et al. 2011; Foucart 2012; Foucart et al. 2018). Simulations in Newtonian gravity show that the NS can also be disrupted over several orbits (Rosswog 2005; Davies et al. 2005). The properties of the progenitors dictate the mass ejected in tidal tails, and the potential formation of a disc wind. In turn, the properties of the ejecta (mass, electron fraction, entropy, and expansion velocity) determine the nucleosynthetic outcome, and hence the contribution that such binaries may make to the heavy element budget of the Universe (see e.g. Rosswog et al. 2017; Just et al. 2015; Roberts et al. 2017). The presence and properties of ejecta and disc determine also the possible formation of a relativistic jet and hence electromagnetic emission as a short GRB. The EM counterparts could also be much more varied than in the NS-NS case (Rosswog 2017).

Finally, the combination of a GW-detected NS-BH binary with an EM counterpart would enable a standard siren measurement of the Hubble constant and other cosmological parameters (e.g., Schutz 1986; Nissanke et al. 2010) out to larger distances than attainable via NS-NS binaries. In what follows, using our limits, we place constraints on the properties of the putative KN and GRB jet that might have been associated with S190814bv.

## 4.1. Constraints on kilonova emission

### 4.1.1. Comparison to AT2017gfo-like kilonovae

Currently, the only KN detected alongside a GW trigger is AT2017gfo, the KN that accompanied GW170817. While we have to be cautious since that source was classified as a NS-NS
<figure><img src="image_10.png" /></figure>

**Fig. 9.** $3\sigma$ or 50% completeness upper limits from the wide-field instrument follow-up campaign. The data used are referenced in Sect. 4.1. The kilonova models, representing an AT2017gfo-like evolution, are shifted to the luminosity distance measure from the LVC skymap (267 Mpc; LALInference.v1.fits), and the shaded regions represent the $1\sigma$ confidence interval (±52 Mpc). Absolute magnitudes assume a distance of 267 Mpc. Foreground extinction is not included.

merger (Abbott et al. 2017a), it is nonetheless prudent to compare it to our limits for S190814bv because it is the only high-confidence KN to date. Foreground<fnref n="19" /> and host galaxy extinction is assumed to be negligible in this analysis.

Figure 9 presents our wide-field follow-up limits (ATLAS, GOTO, PS1, VISTA and VST), plotted against phenomenological fits<fnref n="20" /> to the AT2017gfo light curve based on data from Andreoni et al. (2017), Arcavi et al. (2017a), Chornock et al. (2017), Cowperthwaite et al. (2017), Drout et al. (2017), Evans et al. (2017), Kasliwal et al. (2017), Pian et al. (2017), Smartt et al. (2017), Tanvir et al. (2017), Troja et al. (2017), Utsumi et al. (2017) and Valenti et al. (2017). We find that some of the early VST observations were deep enough to detect a KN of similar brightness to AT2017gfo if one occurred within the $1\sigma$ distance confidence interval. The first VISTA-deep observation also constrains an AT2017gfo-like KN down to the S190814bv distance, and several PS1 frames constrain the near end of the distance distribution. However the large distance to this event precludes a strong statement on whether an AT2017gfo-like event would have been detected by PS1, VISTA or VST. Our deepest limits do exclude KNe (within the relevant frames) similar to those which have been claimed to accompany GRB 130603B (Tanvir et al. 2013; Berger et al. 2013), GRB 050709 (Jin et al. 2016), GRB 060614 (Yang et al. 2015), and GRB 150101B (Gompertz et al. 2018; Troja et al. 2018b). These were all brighter than AT2017gfo at similar epochs to our sampling (Gompertz et al. 2018). Note, though, that some claimed KNe are fainter than AT2017gfo, such as those accompanying GRB 160821B (Lamb et al. 2019b; Troja et al. 2019) and GRB 070809 (Jin et al. 2020).

Our galaxy-targeted observations are able to place significantly deeper limits over a fraction of the error box. In partic-

<figure><img src="image_11.png" /></figure>

**Fig. 10.** Our galaxy targeted limits for S190814bv, alongside the equivalent AT2017gfo KN models. Due to the different distances of the observed galaxies, the data and models are presented in absolute magnitudes. For limits below the model lines, our observations would have uncovered a transient comparable to AT2017gfo, had it been present. We also show the apparent magnitude of the data and models when shifted to the luminosity distance of S190814bv (267 Mpc). Foreground extinction is not included.

ular, early observations, which exceeded $r > 22$ mag, are well below the expected brightness of an AT2017gfo-like KN. Data taken in the IR on timescales of 5–10 days reaching $K > 21$ mag are also competitive. This is illustrated in Fig. 10, which shows the comparison (in absolute magnitudes) between the galaxy-targeted limits and AT2017gfo. Solid lines in the figure show our phenomenological fits to the AT2017gfo light curve in the listed bands, converted to absolute magnitudes at the distances of the targeted galaxies. Downward pointing triangles show the absolute limiting magnitudes of our galaxy-targeted observations. As can be seen, early observations from the WHT and TNG in the $r$-band (Sects. 3.2 and 3.3), and later observations in the IR from HAWK-I are the most constraining (Sect. 3.1), and we are confident in these cases that our observations would have uncovered a KN similar to AT2017gfo if it had been present in the targeted galaxies.

In order to assess more quantitatively the ability of our observations to uncover a putative AT2017gfo-like transient, we can combine our wide-field and galaxy-targeted observations as described in Appendix A. Assuming as our EM counterpart model an AT2017gfo-like event whose flux is scaled by a constant factor, we can derive the covered probability as a function of the ratio between our limiting flux in the most constraining observations and that of AT2017gfo, which is shown in Fig. 11, where we show the covered probability in the $r, i, z$ and $K$ bands (red, purple, blue and cyan lines, respectively), and the combined probability (black line) that corresponds to having a constraining observation in at least one band. For the few, highest probability galaxies the most constraining limit arises from our galaxy targeted programme, but for the majority of the localisation volume the most constraining limits are through wide field observations, in particular from VST ($r$ band), PS1 ($z$-band) and VISTA ($K_s$-band). Our search is therefore sensitive to an AT2017gfo-like KN over ~40% of the localisation probability (as defined in Sect. 3.7), and over ~80% to a transient with the same temporal behaviour, but brighter by a factor of 2. The most constraining observations (due to both depth and coverage) are those in the $r$ and $K_s$ bands.

<footnote n="19">¹⁹ The typical value of $E(B-V)$ over the skymap is <0.1.</footnote>

<footnote n="20">²⁰ As described on the ENGRAVE webpage, http://www.engrave-eso.org, these Bazin et al. (2011) model fits are purely phenomenological (cf. Gompertz et al. 2018), and describe the temporal evolution of AT2017gfo when shifted to the luminosity distance (267 Mpc) of S190814bv.</footnote>
<figure><img src="image_12.png" /><figcaption>**Fig. 11.** Covered probability (defined as described in Appendix A) at which we are sensitive to KN of a given brightness relative to AT2017gfo, based on both our galaxy targeted observations and wide-field limits. Coloured lines show the covered probabilities in four different bands, listed in the legend. The black line shows the combined covered probability.</figcaption></figure>

### 4.1.2. Constraints on the ejecta and on the binary properties

There are fundamental differences between the merger of two NSs (e.g. Ciolfi et al. 2017; Radice et al. 2018; Wollaeger et al. 2018; Shibata & Hotokezaka 2019) and that of a NS and a BH (e.g. Foucart et al. 2013; Kyutoku et al. 2013; Fernández et al. 2017; Tanaka et al. 2014). In the latter case less extreme values for the NS–BH mass ratio, larger BH spin and stiffer NS EoS favour the disruption of the NS before the ISCO, enabling the formation of the accretion disc, tidal tails, and unbound ejecta. This material forms different components from which EM signals can originate (Rantsiou et al. 2008; Pannarale & Ohme 2014; Foucart 2012; Hinderer et al. 2019; Foucart et al. 2019; Barbieri et al. 2019). The KN emission for a given merger is a function of the mass deposited in the various components of the KN, including low-electron-fraction tidal tails (Foucart et al. 2014; Kiuchi et al. 2015; Roberts et al. 2017; Kyutoku et al. 2018) and the neutrino- and viscosity-driven less neutron-rich winds (Fernández & Metzger 2013; Just et al. 2015). It is therefore relevant to compare the observational limits on any KN emission from S190814bv with the expectations of NS–BH models.

We use, for that purpose, the multi-component, anisotropic, NS–BH-specific KN model presented in Barbieri et al. (2019), which builds on the NS–NS KN model of Perego et al. (2017). In this model, three outflow components produce KN emission: (1) the tidal ejecta (which are concentrated close to the orbital plane and have the shape of a crescent); two disc-related winds, namely (2) the neutrino-driven wind from the inner part of the accretion disc and (3) the viscosity-driven wind that results from small scale turbulence of magnetic origin inside the disc. For simplicity, since the neutrino-driven wind (2) is expected to unbind only a small fraction of the disc mass in NS–BH remnants, we neglect that component. As a further simplification, we fix the average (root mean square) velocity $v$ and the (grey) opacity $\kappa$ of the remaining two components to plausible values, namely $v_t = 0.3c$ and $\kappa_t = 15 \text{ cm}^2\text{g}^{-1}$ for the tidal ejecta and $v_w = 0.1c$ and $\kappa_w = 5 \text{ cm}^2\text{g}^{-1}$ for the viscous disc wind, based on their expected velocity and composition: the tidal ejecta are typically expected to retain a very low electron frac-

<figure><img src="image_13.png" /><figcaption>**Fig. 12.** Limits on the tidal ejecta mass ($M_t$) and secular disc wind mass ($M_w$) that we obtain by comparing the NS–BH KN model from Barbieri et al. (2019) to the limits derived from our search (both galaxy-targeted and wide-field). The colour map shows the confidence level at which we can exclude each pair ($M_w, M_t$).</figcaption></figure>

tion $Y_e < 0.2$ (e.g. Fernández et al. 2017) leading to efficient r-process nucleosynthesis of Lanthanides and hence a high opacity $\kappa_t > 10 \text{ cm}^2\text{g}^{-1}$ (Tanaka et al. 2020); disc wind outflows feature a wider range of $Y_e$, due to viscous heating and neutrino irradiation from the inner part of the disc. Differently from NS–NS mergers, though, the absence of shocks and of intense neutrino production by a meta-stable neutron star remnant are likely to cause the disc wind to remain significantly neutron-rich (e.g. Fernández & Metzger 2013; Just et al. 2015). This justifies our choice of $\kappa_w = 5 \text{ cm}^2\text{g}^{-1}$, which is appropriate for outflows with intermediate $Y_e \sim 0.25–0.35$ (Tanaka et al. 2020). We assume $\theta_v = 30^\circ$ as the viewing angle$^{21}$ – measured with respect to the total angular momentum axis – which is the most likely value for a GW-detected inspiral (see Schutz 2011). We are left with the total masses of the two components, $M_t$ and $M_w$, as free parameters. By requiring the resulting light curves to be compatible with our upper limits (both from the galaxy-targeted and from the wide-field searches), following the method outlined in Appendix A, we obtain the constraints on $M_t$ and $M_w$ shown in Fig. 12. The colour map in the figure shows the confidence level at which we can exclude each pair ($M_w, M_t$), denoted as $1 - P(M_w, M_t)$. The region to the lower left of the white dashed line is constrained only by galaxy-targeted observations, while the outer region is constrained mostly by wide-field observations. Large tidal ejecta masses $M_t > 0.05 M_\odot$ are excluded with high confidence $>95\%$, and we can exclude the region $M_t > 0.01 M_\odot$ and $M_w > 0.1 M_\odot$ at approximately one sigma confidence.

By employing numerical-relativity-based fitting formulae that link the properties of the outflows to those of the progenitor binary (Foucart et al. 2018; Kawaguchi et al. 2016), the limits can also be translated into constraints on the NS–BH binary intrinsic properties, again following Barbieri et al. (2019, see also Barbieri et al. 2020). By assuming the disc wind mass to be 30% of the total disc mass (e.g. Fernández & Metzger 2013; Just et al. 2015; Fernández et al. 2019), we take our representative limits on the disc and tidal ejecta masses (corresponding to the region excluded at $1\sigma$ confidence in Fig. 12) to be $M_{\text{disc}} < 0.3 M_\odot$ and $M_t < 10^{-2} M_\odot$, respectively. Figure 13 shows the NS–BH parameter space allowed by our limits, for three

<footnote>$^{21}$ This parameter has only a minor influence on the light curve, so this assumption does not affect our results significantly.</footnote>
<figure><img src="image_14.png" /><figcaption>**Fig. 13.** Constraints on the BH spin ($a_{BH}$) and the BH mass ($M_{BH}$) of the NS–BH binary, assuming remnant disc and tidal ejecta mass limits of $M_{disc} < 0.3 M_{\odot}$ and $M_{t} < 10^{-2} M_{\odot}$, which correspond to approximately 1 sigma exclusion confidence (see text). The SFHo EoS has been adopted to compute the NS tidal deformability.</figcaption></figure>

<figure><img src="image_15.png" /><figcaption>**Fig. 14.** Same as Fig. 13, but for the DD2 equation of state.</figcaption></figure>

possible NS masses, assuming the SFHo NS EoS (Steiner et al. 2013). Different colour shades show the allowed parameter region of the binary for a fixed NS mass (reported near the edge of the region). Figure 14 shows the corresponding limits assuming the DD2 EoS (Typel et al. 2010; Hempel et al. 2012), which is stiffer than SFHo. These two EoSs are representative of the uncertainties in the NS EoS obtained from present nuclear and astrophysical constraints (e.g. Oertel et al. 2017), as well as from constraints derived from GW170817 (Abbott et al. 2018b).

Based on these results, we can therefore exclude that the progenitor NS–BH binary produced a large amount of ejecta. This is consistent with the negligible probability for remnant material left after the merger as reported by the LIGO/Virgo Collaboration in low latency (The LIGO Scientific Collaboration and the Virgo Collaboration 2019a). This indicates that most likely the NS was not disrupted by tidal forces during the final part of the inspiral towards its BH companion, which disfavours high (aligned) BH spins and small mass ratio (or both), as shown quantitatively in Figs. 12–14. At very low mass ratios $q = M_{BH}/M_{NS}$, the maximum allowed BH spin actually increases with decreasing BH mass. Although the NS tidal disruption is more likely for lower BH masses, the frac-

tion of unbound mass (on which our constraint is tighter) is much smaller (Foucart et al. 2019). As a caveat, we note that the dynamical ejecta fitting formula by Kawaguchi et al. (2016) was calibrated only for mass ratios $3 \le M_{BH}/M_{NS} \le 7$ (but see Foucart et al. 2019). We note that a small or absent amount of mass left outside the remnant BH was found to be the most likely outcome of BH–NS mergers in the population synthesis simulations described in Zappa et al. (2019).

## 4.2. Constraints of GRB afterglow-like emission

By assuming that S190814bv launched a short GRB jet, and that all short GRBs have a similar jet structure to that seen in GW170817, we can use the upper limits on any prompt γ-ray emission to constrain the inclination of the system (e.g. Salafia et al. 2019; Saleem et al. 2020; Song et al. 2019). We employ the two jet structures of Lamb et al. (2019a) (see also Resmi et al. 2018; Salafia et al. 2019), both of which are compatible with the afterglow of GRB 170817A, namely a Gaussian and a two-component structure. For both structures, the central-core isotropic-equivalent kinetic energy is $E_{K,iso} = 10^{52}$ erg and the Lorentz factor is $\Gamma = 100$. The two-component structure has “wings” with 10% of the core kinetic isotropic equivalent energy and $\Gamma = 5$. The core half-opening angles are $\theta_c = 0.09$ rad for the Gaussian structure, and $\theta_c = 0.07$ rad for the two-component structure. Figure 15 shows the reported Fermi/GBM upper-limit (Kocevksi 2019 – pink line), assuming a soft, ~1 s duration burst (see Goldstein et al. 2016) at 267 Mpc. Using the two jet structure models described above, the observed isotropic equivalent γ-ray energy for an off-axis observer can be found using the method in Ioka & Nakamura 2019 (which is equivalent to that described in Salafia et al. 2015). We assume a 10% efficiency for energy dissipated as γ-rays by the jet and include opacity due to pair-production where the Lorentz factor is $\Gamma \lesssim 20–30$ following the method in Lamb & Kobayashi (2016, 2017, see also Matsumoto et al. 2019). The top panel of Fig. 15 shows the resulting isotropic-equivalent emitted gamma-ray energy, as a function of the viewing angle, for the Gaussian model (orange dotted line) and the two-component model (blue dash-dotted line) respectively. If there had been a successful GRB 170817A-like jet, the figure shows that the system should be inclined at $>10^{\circ}$, or with a ($\theta_v - \theta_c$) $\ge 5^{\circ}$ where $\theta_c$ is the jet’s core opening angle and $\theta_v$ is the viewing angle from the central axis. The cosmological population of short GRBs typically have an isotropic γ-ray energy in the range $10^{49} \lesssim E_{\gamma,iso} \lesssim 10^{52}$ erg (Fong et al. 2015), thus a successful-GRB producing jet (if any) may have had a lower efficiency or core energy than those assumed here and the lower-limit on the off-axis angle could be smaller.

The afterglow for each structure is shown in the bottom panel of Fig. 15 at an inclination of $10^{\circ}$ and $20^{\circ}$ (thick and thin lines, respectively). The *r* band VST upper limits are shown as red triangles and the *i*- and *z*-band PS1 upper limits as purple and blue triangles, respectively. For the afterglow light curves we assume an ambient density $n = 10^{-3}$ cm$^{-3}$, microphysical parameters $\epsilon_B = \epsilon_c$² = 0.01, and an electron distribution index of *p* = 2.15. In Fig. 15 we only show the *r* band light-curve, noting that the difference in magnitude for *r*–*z* is $\delta m_{r-z} \sim 0.2$ for our model parameters. Our model afterglow light curve is too faint to be constrained by the upper limits. However, for these parameters we can rule out an environment with an ambient density $n \gtrsim 1$ cm$^{-3}$ for a system inclined at ($\theta_v - \theta_c$) $\sim 5^{\circ}$, where we have assumed our energy and microphysical parameters are typical (e.g. Fong et al. 2015; Gompertz et al. 2015).
<figure><img src="image_16.png" /><figcaption>**Fig. 15.** Top panel: the upper limits from Fermi/GBM (Kocevksi 2019) provide an inclination limit for a GRB with a GRB 170817A-like jet structure (pink line with downward arrows). We assume two possible structures following the phenomenology of those in Lamb et al. (2019a), namely a Gaussian (orange) and a two-component structure (blue). The parameters are reported in the text. Bottom panel: r band afterglow for each structure (Gaussian as an orange dotted line, two-component as a blue dash-dotted line) at an inclination of 10° and 20° (thinner lines) assuming an ambient number density $n = 10^{-3}$ cm$^{-3}$, microphysical parameters $\epsilon_B = \epsilon_e^2 = 0.01$, and $p = 2.15$. Upper limits in the r band from VST, and i and z band limits from PS1 are shown as triangles. Models are calculated for a luminosity distance 267 Mpc, the shaded region on the 10° afterglow indicates a ±52 Mpc uncertainty in the luminosity distance.</figcaption></figure>

## 4.3. Comparison to other studies

S190814bv has also been the target of further follow-up reported by other groups. Gomez et al. (2019a) present a study of 96 GLADE galaxies in the 50% error region, which represent 70% of the integrated luminosity in the covered region. They estimate that they cover 25% of all galaxies in the overall localisation area, and are complete down to 0.75 $L^*$ (Schechter 1976) in the region they cover. The typical limiting magnitude they obtain after image subtraction is $i = 22.2$ mag, which is equivalent to $M_i = -14.9$ mag and therefore fainter than AT2017gfo at the distance of S190814bv, at an observing time of ≈ 36 hrs after the GRB. These limits are comparable to the ones we reach with VST in r at a similar time. They rule out KNe with $M_{ej} > 0.01 M_\odot$ in their observations, but the incomplete coverage prevents their observations from being constraining at a high confidence.

Andreoni et al. (2020) present the results of the EM counterpart search by the GROWTH collaboration, using both DECam wide-field tiling observations and targeted spectroscopic and photometric observations of detected candidate transients. All candidates are found to be unrelated SNe. The wide-field coverage is very complete (>98%), with the most constraining limit reached at 3.4 days, $z > 22.3$ mag. This limit is comparable to our PS1 z band limits at a similar time. Comparing their limits to different models, they constrain, depending on model and distance, the ejecta mass to $M_{ej} < 0.03...0.1 M_\odot$, which is consistent with our results.

Watson et al. (2020) present shallow wide-field observations with the DDOTI imager, covering the entire main probability

region (90% of the total probability) down to an unfiltered limit of $w > 17.9$ mag half a day after trigger. They find no candidate transients. They are able to rule out typical on-axis sGRBs but would not have detected a KN similar to AT2017gfo at the distance of S190814bv.

Dobie et al. (2019) present their search for radio transients using ASKAP. They find a single significant transient, AT2019osy, which they suggest is likely associated with a low luminosity AGN. ENGRAVE observations of AT2019osy will be presented in a companion paper in preparation.

Antier et al. (2020) report on rapid early follow-up by the Global Rapid Advanced Network Devoted to the Multi-messenger Addicts (GRANDMA), reaching limits of 17–18.5 mag within the first two hours, earlier but shallower than the GOTO limits presented in this work. These place no constraints on the potential KN emission associated with S190814bv.

Finally, during the revision of this manuscript a preprint was circulated by Vieira et al. (2020), describing the wide-field optical search by the Canada-France-Hawaii Telescope (CFHT). The search reaches limits comparable to ours, despite covering a lower total localisation probability. The corresponding constraints on the putative KN ejecta masses are similar to ours, even though we caution that they are obtained using non BHNS-specific KN models.

## 4.4. Ruling out identified transients as counterparts

The worldwide intensive efforts to identify a counterpart to S190814bv had led to the identification of multiple transients within the error localisation, even though this remains one of the smallest regions available. In part this is due to the deep observations that were capable of identifying transients sources fainter than 22 mag. In addition to the transients identified here through our searches, additional counterparts have been found by other groups (Andreoni et al. 2020; Gomez et al. 2019a; Dobie et al. 2019; Vieira et al. 2020). In total approximately 75 unique optical transients were identified. In principle, each of these should be considered a potential counterpart unless it can be ruled out through follow-up observations. There are various routes that such an approach can take. Firm reasons for rejection include:

1. The identification of the transient source in imaging taken prior to the detection of the GW event (Pre.Det)$^{22}$.

2. A spectrum of the source or host galaxy that places the source outside of the plausible 3D-GW volume (i.e. too distant, or too close, (Spec Host.z)).

3. A spectrum which identifies the source as a different kind of transient event, for which the progenitors are known (SN).

4. The source is actually moving, normally because it is an asteroid, but in one case a high proper motion star (Ast, HPM).

In addition there are further indications which can be used to disfavour sources, but offer a less secure rejection of their association with S190814bv such as:

5. A photometric redshift which is inconsistent with the 3D-GW volume (Phot. Host.z).

6. A lightcurve which does not match the expectations for the counterparts of NS-NS or NS-BH mergers, but is in keeping with a supernova (SN)?

7. No obvious underlying host galaxy (No Host).

8. A source which is nuclear in its host galaxy and therefore likely to be related to AGN activity (AGN?).

<footnote>$^{22}$ The definition in parenthesis is that used in Table C.1.</footnote>
These latter scenarios (5–8) are not as robust as 1–4 since each has potential pitfalls. Photometric redshifts often have significant associated uncertainty and are prone to catastrophic failure. It is possible that a photometric redshift which formally places the event outside the GW horizon could be in error. The use of lightcurves requires some assumptions as to the nature of the electromagnetic emission from the GW event. Since we have only a single well sampled kilonova, and a handful of events identified superposed to short-GRB afterglows, this provides a limited observational comparison. Furthermore, we have no kilonova clearly associated to black-hole neutron star mergers. Nonetheless there are strong reasons to expect low ejecta masses and hence strong limits on the associated luminosities and timescales, hence the photometric evolution can provide a constraint. The lack of an obvious host galaxy would at first sight suggest a distant object, or a large kick to the progenitor. In most cases we would expect to be able to identify a host within ~100 kpc of the transient location. The absence of such a host would disfavour an association with the GW event. However, it should also be noted that some short GRBs arise from an apparently “hostless” population (e.g. Berger 2010; Tunnicliffe et al. 2014), and such an event could be missed. Finally, while most nuclear activity is due to either AGN activity or nuclear starbursts, there are suggestions that mergers could be driven at much higher rates within accretion discs around supermassive black holes (Bartos et al. 2017; Stone et al. 2017). Hence, while nuclear events would apparently be disfavoured as counterparts this should not rule them out. Indeed, the transient AT2019osy, identified as a nuclear radio transient, was an object of interest in the error region of S190814bv. For this reason we separate events which are ruled out for one of the reasons 5–8 from those firmly ruled out via 1–4.

In Appendix C we present a summary of all transient sources identified in the error region of S190814bv by our searches and those of others. We also indicate the reasons that each of these can be rejected (or not). Of the 73 sources presented there 36 are ruled out robustly, 29 are unlikely based on weaker constraints, the remaining 8 have little information to make such distinctions. There is some value in ascertaining if any of these eight events could be plausible counterparts through future observations, for example to obtain host redshifts. However, we also note that these events are not ruled out due to a paucity of observational constraints, rather than any particular diagnostics which would indicate they are likely related to S190814bv. Indeed, given the small error localisation of S190814bv and its relatively high distance, this list of transients provides some indication of the challenge that will remain in identifying robust EM counterparts to GW sources even in the 4-detector era.

## 5. Conclusions

S190814bv was unique amongst the GW detections to date in having an exceptionally small error box. This in turn made it plausible to search for EM emission via targeting of known galaxies. However, such searches were hampered by the large distance to the event (267 ± 52 Mpc compared to ≈40 Mpc for GW170817). Although unprecedented in previous observations, systems like S190814bv may well be more common in the future thanks to both the increasing sensitivity (hence range) of the detectors, and the addition of further GW observatories such as KAGRA in the Kamiokande underground site, Japan, and LIGO-India (Abbott et al. 2018c). Hence, it is relevant to consider what may be the most effective route to the identification of counterparts in this era (corresponding to O4 and beyond). It is striking

that at these distances relatively sensitive wide-field searches such as PS1, VST and VISTA do not, in general, reach sufficiently deep limits to constrain a KN comparable in luminosity to AT2017gfo. This suggests that the current generation of wide-field facilities may not be especially well suited to the majority of candidates in the future, where observations may need to reach $r > 23$ mag to probe a reasonable fraction of KN parameter space (see Sagués Carracedo et al. 2020 and Coughlin et al. 2020 who reach similar conclusions). While some wide-field facilities may be able to attain sufficient depth over a significant fraction of future events (e.g. DECam, BlackGEM, Rubin Observatory LSST) it may well be the case that events at ~300 Mpc may only be detectable by 8 m class observatories. The requirement to observe such events will depend sensitively on where the true event rate of NS-NS and NS-BH lies. At the higher end, optical/IR observers can focus on more nearby events. However, should the event rate lie at the lower end it is necessary to consider if 8 m telescope resources may be needed to identify counterparts, and ELT-like resources required for their follow-up. It therefore continues to be case that the effort should be expended on extending the GW detector network such that the 3D probability volumes for the GW events can become tractable for such observations to be plausible.

**Acknowledgements.** Based on observations collected at the European Southern Observatory under ESO programmes 1102.D-0353(E), 1102.D-0353(F), 1102.D-0353(Q), 1102.D-0353(G), 0103.D-0070(A), 0103.D-0070(B), 0103.D-0703(A), 0103.D-0722(A), 0103.A-9099(A), 198.D-2010(D) and 60.A-9285(A). ATLAS is primarily funded through NEO NASA grants NN12AR55G, 80NSSC18K0284, and 80NSSC18K1575. The ATLAS science products have been made possible through the contributions of the University of Hawaii IfA, the Queen's University Belfast, the Space Telescope Science Institute, and the South African Astronomical Observatory. PanSTARRS is primarily funded through NEO NASA grants NASA Grants NNX08AR22G, NNX14AM74G. The Pan-STARRS science products for LIGO-Virgo follow-up are made possible through the contributions of the University of Hawaii IfA and the Queen's University Belfast. The Gravitational-wave Optical Transient Observer (GOTO) project acknowledges the support of the Monash-Warwick Alliance; Warwick University; Monash University; Sheffield University; Leicester University; Armagh Observatory & Planetarium; the National Astronomical Research Institute of Thailand (NARIT); University of Portsmouth; Turku University and the Instituto de Astrofísica de Canarias (IAC). Part of the funding for GROND was generously granted from the Leibniz-Prize to Prof. G. Hasinger (DFG grant HA 1850/28-1). The Liverpool Telescope is operated on the island of La Palma by Liverpool John Moores University in the Spanish Observatorio del Roque de los Muchachos of the Instituto de Astrofísica de Canarias with financial support from the UK Science and Technology Facilities Council. The WHT and its over-ride programme are operated on the island of La Palma by the Isaac Newton Group of Telescopes in the Spanish Observatorio del Roque de los Muchachos of the Instituto de Astrofísica de Canarias; part of these data were taken under program (19A)N3. FEB thanks CONICYT Basal AFB-170002 and Chile's Ministry of Economy fund IC120009. MGB, PDA and AM acknowledge support from ASI grant I/004/11/3. MBr, EC, AP and SpI acknowledge support from MIUR (PRIN 2017 grant 20179ZF5KS). EB, EM and MT acknowledge funding from GRAWITA. SHB is indebted to the Danish National Research Foundation (DNRF132) for support. SCA acknowledges support from grant MAE0065741. EC acknowledges the support of the H2020 OPTICON programme 730890. TWC acknowledges the Humboldt Foundation and Marie Skłodowska-Curie grant 842471. MDP thanks Istanbul University for support. PAE acknowledges UKSA support. RAJEF is supported by an STFC studentship. MF is supported by a Royal Society - SFI University Research Fellowship. LG was funded by the EU H2020 programme under MSCA grant no. 839090. CG, JH and LI were supported by a research grant from VILLUM FONDEN (project 16599). CG and LI were supported by a research grant from VILLUM FONDEN (25501). GGh acknowledges the PRIN MIUR "Figaro" for financial support. AGo acknowledges financial support from the Slovenian Research Agency (grants P1-0031, I0-0033, and J1-8136). BPG, AJL and JDL acknowledge support from ERC grant 725246 (TEDE, PI Levan). SGG acknowledges support by FCT Fundação para a Ciência e Tecnologia and by Project PTDC/FIS-AST-31546. GGr acknowledges the ESCAPE H2020 project no. 824064. MG is supported by the Polish NCN MAESTRO grant 2014/14/A/ST9/00121. PJG acknowledges support from NOVA and from the South African NRF SARChI grant 111692.
CPG and MS acknowledge support from EU/FP7-ERC grant no. 615929. KEH acknowledges support by a Project Grant from The Icelandic Research Fund. YDH acknowledges support from the China Scholarships Council. JJ acknowledges support from NOVA and NWO-FAPESP grant for instrumentation. AJ acknowledges funding from the European Research Council (ERC). ZPJ was supported by the Foundation for Distinguished Young Scholars of Jiangsu Province (no. BK20180050). PGJ acknowledges funding from the ERC under Consolidator Grant agreement no. 647208. DAK acknowledges Spanish research project RTI2018-098104-J100 (GRBPhot). SKI acknowledges support by DFG grant K1 766/16-3. ECK acknowledges support from the GREAT research environment. GPL acknowledges support from STFC via grant ST/N000757/1. GL was supported by a research grant (19054) from VILLUM FONDEN. KM acknowledges support from the ERC (grant no. 758638). IM is partially supported by OzGrav (ARC project CE17010000). MM acknowledges support from ERC through ERC-2017-CoG no. 770017. MJM acknowledges the National Science Centre, Poland, grant 2018/30/E/ST9/00208. BMJ and DW are supported by Independent Research Fund Denmark grant DFF-7014-00017. MN is supported by a Royal Astronomical Society Research Fellowship. ANG acknowledges support by grant DFG K1 766/16-3. PTOB acknowledges funding from STFC. SRO gratefully acknowledges the support of the Leverhulme Trust. FO acknowledges the support of the H2020 Hemera program, grant no. 730970. MAPT was supported by grants RYC-2015-17854 and AYA2017-83216-P. EP acknowledges financial support from INAF. GP is supported by the Millennium Science Initiative through grant IC120009. MLP is partially supported by a "Linea 2" project of the Catania University. JQV acknowledges support from CONICYT folio 21180886. TMR acknowledges the support of the Vilho, Yrjo and Kalle Vaisala Foundation. ARO acknowledges support from Premiale LBT 2013. SR is supported by VR grants 2016-03657_3 and the research environment grant GREAT, Dnr. 2016-06012, and the Swedish National Space board, Dnr. 107/16. OSS acknowledges the Italian Ministry of Research (MIUR) grant 1.05.06.13. LSa acknowledges the Irish Research Council Scholarship no. GOIPG/2017/1525. SJS acknowledges support from STFC Grant ST/P000312/1. ERS and DS acknowledge funding from UK STFC CG ST/P000495/1. RLCS acknowledges funding from STFC. DS acknowledges support from STFC via grant ST/T007184/1. SDV acknowledges the support of the CNES. LW supported by Polish NCN DAINA 2017/27/L/ST9/03221. The Cosmic DAWN center is funded by the Danish National Research Foundation.

## References

Abadie, J., Abbott, B. P., Abbott, R., et al. 2010, Classical Quantum Gravity, 27, 173001

Abazajian, K. N., Adelman-McCarthy, J. K., Agüeros, M. A., et al. 2009, ApJS, 182, 543

Abbott, B. P., Abbott, R., Abbott, T. D., et al. 2016, Phys. Rev. Lett., 116, 061102

Abbott, B. P., Abbott, R., Abbott, T. D., et al. 2017a, Phys. Rev. Lett., 119, 161101

Abbott, B. P., Abbott, R., Abbott, T. D., et al. 2017b, ApJ, 848, L13

Abbott, B. P., Abbott, R., Abbott, T. D., et al. 2017c, ApJ, 848, L12

Abbott, B. P., Abbott, R., Abbott, T. D., et al. 2017d, Nature, 551, 85

Abbott, T. M. C., Abdalla, F. B., Allam, S., et al. 2018a, ApJS, 239, 18

Abbott, B. P., Abbott, R., Abbott, T. D., et al. 2018b, Phys. Rev. Lett., 121, 161101

Abbott, B. P., Abbott, R., Abbott, T. D., et al. 2018c, Liv. Rev. Relativ., 21, 3

Ackley, K., Dyer, M., Eyles, R., et al. 2019, GRB Coordinates Network, 25337, 1

Acosta-Pulido, J., Ballesteros, E., Barreto, M., et al. 2002, Newsletter Isaac Newton Group Telescopes, 6, 22

Ageron, M., Baret, B., Coleiro, A., et al. 2019, GRB Coordinates Network, 25330, 1

Alam, S., Albareti, F. D., Allende Prieto, C., et al. 2015, ApJS, 219, 12

Alvarez-Muniz, J., Pedreira, F., Zas, E., Hampert, K. H., & Schimp, M. 2019, GRB Coordinates Network, 25409, 1

Andreoni, I., Ackley, K., Cooke, J., et al. 2017, PASA, 34, e069

Andreoni, I., Goldstein, D. A., Ahumada, T., et al. 2019a, GRB Coordinates Network, 25362, 1

Andreoni, I., Goldstein, D. A., Dobie, D., & Kasliwal, M. M. 2019b, GRB Coordinates Network, 25488, 1

Andreoni, I., Goldstein, D. A., Kasliwal, M. M., et al. 2020, ApJ, 890, 131

Antier, S., Agayeva, S., Aivazyan, V., et al. 2020, MNRAS, 492, 3904

Appenzeller, I., Fricke, K., Fürtig, W., et al. 1998, The Messenger, 94, 1

Arcavi, I., Hosseinzadeh, G., Howell, D. A., et al. 2017a, Nature, 551, 64

Arcavi, I., McCully, C., Hosseinzadeh, G., et al. 2017b, ApJ, 848, L33

Arnaboldi, M., Petr-Gotzens, M., Rejkuba, M., et al. 2010, The Messenger, 139, 6

Arsenault, R., Hubin, N., Stroebele, S., et al. 2006, The Messenger, 123, 6

Artale, M. C., Mapelli, M., Giacobbo, N., et al. 2019, MNRAS, 487, 1675

Barbieri, C., Salafia, O. S., Perego, A., Colpi, M., & Ghirlanda, G. 2019, A&A, 625, A152

Barbieri, C., Salafia, O. S., Perego, A., Colpi, M., & Ghirlanda, G. 2020, Eur. Phys. J. A, 56, 8

Bartos, I., Kocsis, B., Haiman, Z., & Márka, S. 2017, ApJ, 835, 165

Bauer, F. E., Fruchter, A. S., Gonzalez Lopez, J., et al. 2019, GRB Coordinates Network, 25801, 1

Bazin, G., Ruhlmann-Kleider, V., Palanque-Delabrouille, N., et al. 2011, A&A, 534, A43

Becker, A. 2015, Astrophysics Source Code Library [record ascl:1504.004]

Belczynski, K., Holz, D. E., Bulik, T., & O'Shaughnessy, R. 2016, Nature, 534, 512

Benn, C., Dee, K., & Agócs, T. 2008, in CAMAC: A New Imager/Spectrograph for the William Herschel Telescope, SPIE Conf. Ser., 7014, 70146X

Berger, E. 2010, ApJ, 722, 1946

Berger, E., Fong, W., & Chornock, R. 2013, ApJ, 774, L23

Bertin, E., & Arnouts, S. 1996, A&AS, 117, 393

Bloom, J. S., Richards, J. W., Nugent, P. E., et al. 2012, PASP, 124, 1175

Boksenberg, A. 1985, Vistas Astron, 28, 531

Brocato, E., Branchesi, M., Cappellaro, E., et al. 2018, MNRAS, 474, 411

Brunn, S. H., Sagues Carracedo, A., Chen, T.-W., et al. 2019, GRB Coordinates Network, 25384, 1

Cai, C., Yi, Q. B., Xiao, S., et al. 2019, GRB Coordinates Network, 25365, 1

Capaccioli, M., & Schipani, P. 2011, The Messenger, 146, 2

Cartier, R., Briceno, C., Olivares, F., et al. 2019, GRB Coordinates Network, 25784, 1

Casali, M., Pirard, J.-F., Kissler-Patig, M., et al. 2006, in SPIE Conf. Ser., Proc. SPIE, 6269, 62690W

Castro-Tirado, A. J., Valeev, A. F., Hu, Y.-D., et al. 2019, GRB Coordinates Network, 25543, 1

Chambers, K. C., Magnier, E. A., Metcalfe, N., et al. 2016, ArXiv e-prints [arXiv:1612.05560]

Chen, T.-W., Nicuesa Guelbenzu, A., Fraser, M., et al. 2019, GRB Coordinates Network, 25372, 1

Chornock, R., Berger, E., Kasen, D., et al. 2017, ApJ, 848, L19

Ciofi, R., Kastaun, W., Giacomazzo, B., et al. 2017, Phys. Rev. D, 95, 063016

Colless, M., Dalton, G., Maddox, S., et al. 2001, MNRAS, 328, 1039

Coughlin, M. W., Dietrich, T., Antier, S., et al. 2020, MNRAS, 492, 863

Coulter, D. A., Foley, R. J., Kilpatrick, C. D., et al. 2017a, Science, 358, 1556

Coulter, D., Kilpatrick, C., Siebert, M., et al. 2017b, GRB Coordinates Network, 21529, 1

Covino, S., Wiersema, K., Fan, Y. Z., et al. 2017, Nat. Astron., 1, 791

Cowperthwaite, P. S., Berger, E., Villar, V. A., et al. 2017, ApJ, 848, L17

Dálya, G., Galgóczi, G., Dobos, L., et al. 2018, MNRAS, 479, 2374

D'Avanzo, P., Campana, S., Salafia, O. S., et al. 2018, A&A, 613, L1

D'Avanzo, P., Melandri, A., Izzo, L., et al. 2019a, GRB Coordinates Network, 25331, 1

D'Avanzo, P., Rossi, A., Greco, G., et al. 2019b, GRB Coordinates Network, 25361, 1

Davies, M. B., Levan, A. J., & King, A. R. 2005, MNRAS, 356, 54

De, K., Jencson, J., Kasliwal, M. M., Goldstein, D. A., & Andreoni, I. 2019a, GRB Coordinates Network, 25478, 1

De, K., Goldstein, D., Andreoni, I., et al. 2019b, GRB Coordinates Network, 25348, 1

Dichiarì, S., Troja, E., Cenko, S. B., et al. 2019, GRB Coordinates Network, 25374, 1

Dimitriadis, G., Brown, J. S., Seibert, M. R., et al. 2019, GRB Coordinates Network, 25395, 1

Dobie, D., Kaplan, D. L., Murphy, T., et al. 2018, ApJ, 858, L15

Dobie, D., Stewart, A., Murphy, T., et al. 2019, ApJ, 887, L13

Dominik, M., Berti, E., O'Shaughnessy, R., et al. 2015, ApJ, 806, 263

Drout, M. R., Piro, A. L., Shappee, B. J., et al. 2017, Science, 358, 1570

Dyer, M. J., Dhillon, V. S., Littlefair, S., et al. 2018, in Observatory Operations: Strategies, Processes, and Systems VII, SPIE Conf. Ser., 10704, 107040C

Eichler, D., Livio, M., Piran, T., & Schramm, D. N. 1989, Nature, 340, 126

Eldridge, J. J., & Stanway, E. R. 2016, MNRAS, 462, 3302

Evans, P. A., Osborne, J. P., Kennea, J. A., et al. 2016a, MNRAS, 455, 1522

Evans, P. A., Kennea, J. A., Barthelmy, S. D., et al. 2016b, MNRAS, 460, L40

Evans, P. A., Cenko, S. B., Kennea, J. A., et al. 2017, Science, 358, 1565

Evans, P. A., Kennea, J. A., Tohuvavohu, A., et al. 2019, GRB Coordinates Network, 25400, 1

Fernández, R., & Metzger, B. D. 2013, MNRAS, 435, 502

Fernández, R., Foucart, F., Kasen, D., et al. 2017, Classical Quantum Gravity, 34, 154001

Fernández, R., Tchekhovskoy, A., Quataert, E., Foucart, F., & Kasen, D. 2019, MNRAS, 482, 3373

Flesch, E. W. 2015, PASA, 32, e010
Flewelling, H. A., Magnier, E. A., Chambers, K. C., et al. 2016, ArXiv e-prints
[arXiv:1612.05243]

Fong, W., & Berger, E. 2013, ApJ, 776, 18

Fong, W., Berger, E., Margutti, R., & Zauderer, B. A. 2015, ApJ, 815, 102

Foucart, F. 2012, Phys. Rev. D, 86, 124007

Foucart, F., Deaton, M. B., Duez, M. D., et al. 2013, Phys. Rev. D, 87, 084006

Foucart, F., Deaton, M. B., Duez, M. D., et al. 2019, Phys. Rev. D, 90, 024026

Foucart, F., Hinderer, T., & Nissanke, S. 2018, Phys. Rev. D, 98, 081501

Foucart, F., Duez, M. D., Kidder, L. E., et al. 2019, Phys. Rev. D, 99, 103025

Freiburghaus, C., Rosswog, S., & Thielemann, F.-K. 1999, ApJ, 525, L121

Freudling, W., Romaniello, M., Bramich, D. M., et al. 2013, A&A, 559, A96

Gaia Collaboration (Brown, A. G. A., et al.) 2016, A&A, 595, A2

Gall, C., Hjorth, J., Rosswog, S., Tanvir, N. R., & Levan, A. J. 2017, ApJ, 849, L19

Gehrels, N., Cannizzo, J. K., Kanner, J., et al. 2016, ApJ, 820, 136

Ghirlanda, G., Salafia, O. S., Paragi, Z., et al. 2019, Science, 363, 968

Giacobbo, N., & Mapelli, M. 2018, MNRAS, 480, 2011

Goldstein, D. A., & Anand, S. 2019, GRB Coordinates Network, 25394, 1

Goldstein, A., Burns, E., Hamburg, R., et al. 2016, ArXiv e-prints
[arXiv:1612.02395]

Goldstein, A., Veres, P., Burns, E., et al. 2017, ApJ, 848, L14

Goldstein, D., Perley, D. A., Andreoni, I., & Kasliwal, M. M. 2019a, GRB Coordinates Network, 25355, 1

Goldstein, D. A., Andreoni, I., Zhou, R., et al. 2019b, GRB Coordinates Network, 25391, 1

Goldstein, D. A., Andreoni, I., Hankins, M., et al. 2019c, GRB Coordinates Network, 25393, 1

Gomez, S., Hosseinzadeh, G., Cowperthwaite, P. S., et al. 2019a, ApJ, 884, L55

Gomez, S., Hosseinzadeh, G., Berger, E., et al. 2019b, GRB Coordinates Network, 25483, 1

Gompertz, B. P., van der Horst, A. J., O'Brien, P. T., Wynn, G. A., & Wiersema, K. 2015, MNRAS, 448, 629

Gompertz, B. P., Levan, A. J., Tanvir, N. R., et al. 2018, ApJ, 860, 62

Gompertz, B. P., Cutter, R., Steeghs, D., et al. 2020a, MNRAS, 497, 726

Gompertz, B. P., Levan, A. J., & Tanvir, N. R. 2020b, ApJ, 895, 58

González-Fernández, C., Hodgkin, S. T., Irwin, M. J., et al. 2018, MNRAS, 474, 5459

Grado, A., Cappellaro, E., Getman, F., et al. 2019a, GRB Coordinates Network, 25371, 1

Grado, A., Cappellaro, E., Brocato, E., et al. 2019b, GRB Coordinates Network, 25669, 1

Grado, A., Cappellaro, E., Covino, S., et al. 2020, MNRAS, 492, 1731

Greiner, J., Bornemann, W., Clemens, C., et al. 2008, PASP, 120, 405

Hanna, C., Mandel, I., & Vousden, W. 2014, ApJ, 784, 8

Hempel, M., Fischer, T., Schaffner-Bielich, J., & Liebendörfer, M. 2012, ApJ, 748, 70

Herner, K., Palmese, A., Soares-Santos, M., et al. 2019a, GRB Coordinates Network, 25373, 1

Herner, K., Palmese, A., Soares-Santes, M., et al. 2019b, GRB Coordinates Network, 25398, 1

Hill, D. T., Driver, S. P., Cameron, E., et al. 2010, MNRAS, 404, 1215

Hinderer, T., Nissanke, S., Foucart, F., et al. 2019, Phys. Rev. D, 100, 063021

Hjorth, J., Levan, A. J., Tanvir, N. R., et al. 2017, ApJ, 848, L31

Hu, Y.-D., Castro-Tirado, A. J., Valeev, A. F., et al. 2019, GRB Coordinates Network, 25588, 1

Huber, M., Smith, K. W., Chambers, K., et al. 2019, GRB Coordinates Network, 25356, 1

IceCube Collaboration 2019, GRB Coordinates Network, 25321, 1

Ioka, K., & Nakamura, T. 2019, MNRAS, 487, 4884

Japelj, J., Jin, Z., Kankare, E., et al. 2019a, GRB Coordinates Network, 25447, 1

Japelj, J., Kankare, E., Kool, E., et al. 2019b, GRB Coordinates Network, 25526, 1

Jin, Z.-P., Hotokezaka, K., Li, X., et al. 2016, Nat. Commun., 7, 12898

Jin, Z.-P., Covino, S., Liao, N.-H., et al. 2020, Nat. Astron., 4, 77

Jonker, P., Maguire, K., Fraser, M., et al. 2019, GRB Coordinates Network, 25454, 1

Just, O., Bauswein, A., Ardevol Pulpillo, R., Goriely, S., & Janka, H.-T. 2015, MNRAS, 448, 541

Kasen, D., Metzger, B., Barnes, J., Quataert, E., & Ramirez-Ruiz, E. 2017, Nature, 551, 80

Kasliwal, M. M., Nakar, E., Singer, L. P., et al. 2017, Science, 358, 1559

Kawaguchi, K., Kyutoku, K., Shibata, M., & Tanaka, M. 2016, ApJ, 825, 52

Kissler-Patig, M., Pirard, J.-F., Casali, M., et al. 2008, A&A, 491, 941

Kiuchi, K., Sekiguchi, Y., Kyutoku, K., et al. 2015, Phys. Rev. D, 92, 064034

Kocevksi, D. 2019, GRB Coordinates Network, 25326, 1

Kruckow, M. U., Tauris, T. M., Langer, N., Kramer, M., & Izzard, R. G. 2018, MNRAS, 481, 1908

Krühler, T., Küpcü Yoldas, A., Greiner, J., et al. 2008, ApJ, 685, 376

Kuijken, K. 2011, The Messenger, 146, 8

Kulkarni, S. R., Perley, D. A., & Miller, A. A. 2018, ApJ, 860, 22

Kyutoku, K., Okawa, H., Shibata, M., & Taniguchi, K. 2011, Phys. Rev. D, 84, 064018

Kyutoku, K., Ioka, K., & Shibata, M. 2013, Phys. Rev. D, 88, 041503

Kyutoku, K., Kiuchi, K., Sekiguchi, Y., Shibata, M., & Taniguchi, K. 2018, Phys. Rev. D, 97, 023009

Lamb, G. P., & Kobayashi, S. 2016, ApJ, 829, 112

Lamb, G. P., & Kobayashi, S. 2017, MNRAS, 472, 4953

Lamb, G. P., Lyman, J. D., Levan, A. J., et al. 2019a, ApJ, 870, L15

Lamb, G. P., Tanvir, N. R., Levan, A. J., et al. 2019b, ApJ, 883, 48

Lattimer, J. M., & Schramm, D. N. 1976, ApJ, 210, 549

Lattimer, J. M., Mackie, F., Ravenhall, D. G., & Schramm, D. N. 1977, ApJ, 213, 225

Lazzati, D., Perna, R., Morsony, B. J., et al. 2018, Phys. Rev. Lett., 120, 241103

Levan, A. J., Lyman, J. D., Tanvir, N. R., et al. 2017, ApJ, 848, L28

Li, L.-X., & Paczyński, B. 1998, ApJ, 507, L59

Lipunov, V. M., Gorbovskoy, E., Kornilov, V. G., et al. 2017, ApJ, 850, L1

Lopez-Cruz, O., Castro-Tirado, A. J., Macri, L., et al. 2019a, GRB Coordinates Network, 25419, 1

Lopez-Cruz, O., Castro-Tirado, A. J., Macri, L., et al. 2019b, GRB Coordinates Network, 25571, 1

Lyman, J. D., Lamb, G. P., Levan, A. J., et al. 2018, Nat. Astron., 2, 751

Magnier, E.A., Chambers, K.C., Flewelling, H.A., et al. 2016, ArXiv e-prints
[arXiv:1612.05240]

Margutti, R., Alexander, K. D., Xie, X., et al. 2018, ApJ, 856, L18

Matsumoto, T., Nakar, E., & Piran, T. 2019, MNRAS, 486, 1563

McCully, C., Hiramatsu, D., Howell, D. A., et al. 2017, ApJ, 848, L32

Metzger, B. D., Martínez-Pinedo, G., Darbha, S., et al. 2010, MNRAS, 406, 2650

Molinari, E., Conconi, P., Pucillo, M., & Monai, S. 1999, in Looking Deep in the Southern Sky, eds. R. Morganti, & W. J. Couch, 157

Molkov, S., Mereghetti, S., Savchenko, V., et al. 2019, GRB Coordinates Network, 25323, 1

Mooley, K. P., Frail, D. A., Dobie, D., et al. 2018, ApJ, 868, L11

Neijssel, C. J., Vigna-Gómez, A., Stevenson, S., et al. 2019, MNRAS, 490, 3740

Nicholl, M., Berger, E., Kasen, D., et al. 2017, ApJ, 848, L18

Nissanke, S., Holz, D. E., Hughes, S. A., Dalal, N., & Sievers, J. L. 2010, ApJ, 725, 496

Nissanke, S., Kasliwal, M., & Georgieva, A. 2013, ApJ, 767, 124

Nynka, M., Ruan, J. J., Haggard, D., & Evans, P. A. 2018, ApJ, 862, L19

Oertel, M., Hempel, M., Klähn, T., & Typel, S. 2017, Rev. Mod. Phys., 89, 015007

Ohno, M., Axelsson, M., Longo, F., et al. 2019, GRB Coordinates Network, 25386, 1

Oliva, E., Marconi, A., Maiolino, R., et al. 2001, A&A, 369, L5

Palmer, D. M., Barthelmy, S. D., Lien, A. Y., et al. 2019, GRB Coordinates Network, 25341, 1

Pannarale, F., & Ohme, F. 2014, ApJ, 791, L7

Perego, A., Radice, D., & Bernuzzi, S. 2017, ApJ, 850, L37

Pian, E., D'Avanzo, P., Benetti, S., et al. 2017, Nature, 551, 67

Pilia, M., Pittori, C., Tavani, M., et al. 2019a, GRB Coordinates Network, 25327, 1

Pilia, M., Pittori, C., Tavani, M., et al. 2019b, GRB Coordinates Network, 25335, 1

Pirard, J.-F., Kissler-Patig, M., Moorwood, A., et al. 2004, in Ground-based Instrumentation for Astronomy, eds. A. F. M. Moorwood, & M. Iye, Proc. SPIE, 5492, 1763

Poretti, E. 2018, https://www.doi.org/10.5281/zenodo.1892959

Radice, D., Perego, A., Hotokezaka, K., et al. 2018, ApJ, 869, 130

Rantsiou, E., Kobayashi, S., Laguna, P., & Rasio, F. A. 2008, ApJ, 680, 1326

Resmi, L., Schulze, S., Ishwara-Chandra, C. H., et al. 2018, ApJ, 867, 57

Roberts, L. F., Lippuner, J., Duez, M. D., et al. 2017, MNRAS, 464, 3907

Rodriguez, A., Meza-Retamal, N., Quirola, J., et al. 2019, GRB Coordinates Network, 25423, 1

Rossi, A., D'Avanzo, P., Cappellaro, E., et al. 2019, GRB Coordinates Network, 25383, 1

Rosswog, S. 2005, ApJ, 634, 1202

Rosswog, S. 2017, Phys. Online J., 10, 131

Rosswog, S., Feindt, U., Korobkin, O., et al. 2017, Classical Quantum Gravity, 34, 104001

Ruan, J., Vieira, N., Haggard, D., et al. 2019, GRB Coordinates Network, 25443, 1

Sagués Carracedo, A., Bulla, M., Feindt, U., & Goobar, A. 2020, MNRAS, submitted [arXiv:2004.06137]

Salafia, O. S., Ghisellini, G., Pescalli, A., Ghirlanda, G., & Nappo, F. 2015, MNRAS, 450, 3549

Salafia, O. S., Ghirlanda, G., Ascenzi, S., & Ghisellini, G. 2019, A&A, 628, A18
Saleem, M., Resmi, L., Arun, K. G., & Mohan, S. 2020, ApJ, 891, 130

Salmon, L., Hanlon, L., Jeffrey, R. M., & Martin-Carrillo, A. 2020, A&A, 634, A32

Savchenko, V., Ferrigno, C., Kuulkers, E., et al. 2017, ApJ, 848, L15

Schechter, P. 1976, ApJ, 203, 297

Schutz, B. F. 1986, Nature, 323, 310

Schutz, B. F. 2011, Classical Quantum Gravity, 28, 125023

Shappee, B. J., Simon, J. D., Drout, M. R., et al. 2017, Science, 358, 1574

Shibata, M., & Hotokezaka, K. 2019, Ann. Rev. Nucl. Part. Sci., 69, 41

Shibata, M., & Taniguchi, K. 2011, Liv. Rev. Relativ., 14, 6

Siebenmorgen, R., Carraro, G., Valenti, E., et al. 2011, The Messenger, 144, 9

Singer, L. P., Chen, H.-Y., Holz, D. E., et al. 2016, ApJ, 829, L15

Skrutskie, M. F., Cutri, R. M., Stiening, R., et al. 2006, AJ, 131, 1163

Smartt, S. J., Chambers, K. C., Smith, K. W., et al. 2016a, MNRAS, 462, 4094

Smartt, S. J., Chambers, K. C., Smith, K. W., et al. 2016b, ApJ, 827, L40

Smartt, S. J., Chen, T.-W., Jerkstrand, A., et al. 2017, Nature, 551, 75

Smartt, S. J., Smith, K. W., Srivastav, S., et al. 2019a, GRB Coordinates Network, 25455, 1

Smartt, S. J., Malesani, D. B., Smith, K. W., et al. 2019b, GRB Coordinates Network, 25386, 1

Soares-Santos, M., Holz, D. E., Annis, J., et al. 2017, ApJ, 848, L16

Soares-Santos, M., Tucker, D., Allam, S., et al. 2019a, GRB Coordinates Network, 25336, 1

Soares-Santos, M., Annis, J., Garcia, A., et al. 2019b, GRB Coordinates Network, 25425, 1

Soares-Santos, M., Annis, J., Garcia, A., et al. 2019c, GRB Coordinates Network, 25438, 1

Soares-Santos, M., Annis, J., Herner, K., et al. 2019d, GRB Coordinates Network, 25486, 1

Song, H.-R., Ai, S.-K., Wang, M.-H., et al. 2019, ApJ, 881, L40

Srivastas, S., Huber, M., Smartt, S. J., et al. 2019, GRB Coordinates Network, 25417, 1

Stalder, B., Tonry, J., Smartt, S. J., et al. 2017, ApJ, 850, 149

Steele, I. A., Smith, R. J., Rees, P. C., et al. 2004, in The Liverpool Telescope: Performance and First Results, SPIE Conf. Ser., 5489, 679

Steiner, A. W., Hempel, M., & Fischer, T. 2013, ApJ, 774, 17

Stewart, A., Dobie, D., Murphy, T., et al. 2019, GRB Coordinates Network, 25487, 1

Stone, N. C., Metzger, B. D., & Haiman, Z. 2017, MNRAS, 464, 946

Sugizaki, M., Kawai, N., Negoro, H., et al. 2019, GRB Coordinates Network, 25329, 1

Sutherland, W., Emerson, J., Dalton, G., et al. 2015, A&A, 575, A25

Svinkin, D., Golenetskii, S., Aptekar, R., et al. 2019, GRB Coordinates Network, 25369, 1

Tanaka, M., Hotokezaka, K., Kyutoku, K., et al. 2014, ApJ, 780, 31

Tanaka, M., Kato, D., Gaigalas, G., & Kawaguchi, K. 2020, MNRAS, 496, 1369

Tanvir, N. R., Levan, A. J., Fruchter, A. S., et al. 2013, Nature, 500, 547

Tanvir, N. R., Levan, A. J., González-Fernández, C., et al. 2017, ApJ, 848, L27

The LIGO Scientific Collaboration and the Virgo Collaboration 2019a, GRB Coordinates Network, 25324, 1

The LIGO Scientific Collaboration and the Virgo Collaboration 2019b, GRB Coordinates Network, 25333, 1

Tonry, J. L., Stubbs, C. W., Lykke, K. R., et al. 2012, ApJ, 750, 99

Tonry, J. L., Denneau, L., Heinze, A. N., et al. 2018a, PASP, 130, 064505

Tonry, J. L., Denneau, L., Flewelling, H., et al. 2018b, ApJ, 867, 105

Troja, E., Piro, L., van Eerten, H., et al. 2017, Nature, 551, 71

Troja, E., Piro, L., Ryan, G., et al. 2018a, MNRAS, 478, L18

Troja, E., Ryan, G., Piro, L., et al. 2018b, Nat. Commun., 9, 4089

Troja, E., Castro-Tirado, A. J., Becerra González, J., et al. 2019, MNRAS, 489, 2104

Tucker, D., Butner, M., Wiesner, M., et al. 2019a, GRB Coordinates Network, 25484, 1

Tucker, D., Allam, S., Wiesner, M., et al. 2019b, GRB Coordinates Network, 25379, 1

Tunnicliffe, R. L., Levan, A. J., Tanvir, N. R., et al. 2014, MNRAS, 437, 1495

Typel, S., Röpke, G., Klähn, T., Blaschke, D., & Wolter, H. H. 2010, Phys. Rev. C, 81, 015803

Utsumi, Y., Tanaka, M., Tominaga, N., et al. 2017, PASJ, 69, 101

Valenti, S., David, J., Yang, S., et al. 2017, ApJ, 848, L24

Veitch, J., Raymond, V., Farr, B., et al. 2015, Phys. Rev. D, 91, 042003

Véron-Cetty, M. P., & Véron, P. 2001, A&A, 374, 92

Vieira, N., Ruan, J. J., Haggard, D., et al. 2020, ApJ, 895, 96

Waters, C. Z., Magnier, E. A., Price, P. A., et al. 2016, ArXiv e-prints [arXiv:1612.05245]

Watson, D., Hansen, C. J., Selsing, J., et al. 2019, Nature, 574, 497

Watson, A. M., Butler, N. R., Lee, W. H., et al. 2020, MNRAS, 492, 5916

Wiesner, M., Butner, M., Allam, S., et al. 2019a, GRB Coordinates Network, 25596, 1

Wiesner, M., Butner, M., Tucker, D., et al. 2019b, GRB Coordinates Network, 25540, 1

Willmer, C. N. A. 2018, ApJS, 236, 47

Wollaeger, R. T., Korobkin, O., Fontes, C. J., et al. 2018, MNRAS, 478, 3298

Yang, S. 2018, Ph.D. thesis, http://paduaresearch.cab.unipd.it/11854/

Yang, B., Jin, Z.-P., Li, X., et al. 2015, Nat. Commun., 6, 7323

Yang, S., Cappellaro, E., Grado, A., et al. 2019, GRB Coordinates Network, 25748, 1

Zappa, F., Bernuzzi, S., Pannarale, F., Mapelli, M., & Giacobbo, N. 2019, Phys. Rev. Lett, 123, 041102

1 School of Physics and Astronomy, Monash University, Clayton, Victoria 3800, Australia

2 INAF – Osservatorio di Astrofisica e Scienza dello Spazio di Bologna, Via Piero Gobetti 93/3, 40129 Bologna, Italy

3 INAF / Brera Astronomical Observatory, Via Bianchi 46, 23807 Merate, LC, Italy

4 INFN – Sezione di Milano-Bicocca, Piazza della Scienza 3, 20126 Milano, MI, Italy

5 University of Milano-Bicocca, Department of Physics “G. Occhiali”, Piazza della Scienza 3, 20126 Milano, Italy

6 Instituto de Astrofísica and Centro de Astroingeniería, Facultad de Física, Pontificia Uni versidad Católica de Chile, Casilla 306, Santiago 22, Chile

7 Millennium Institute of Astrophysics (MAS), Nuncio Monseñor Sótero Sanz 100, Providencia, Santiago, Chile

8 Space Science Institute, 4750 Walnut Street, Suite 205, Boulder, CO 80301, USA

9 INAF, Osservatorio Astronomico di Padova, 35122 Padova, Italy

10 Space Telescope Science Institute, 3700 San Martin Drive, Baltimore, MD 21218, USA

11 INAF Osservatorio Astronomico di Capodimonte, Via Moiariello 16, 80131 Napoli, Italy

12 Gran Sasso Science Institute, Viale F. Crispi 7, 67100 L’Aquila, AQ, Italy

13 INFN – Laboratori Nazionali del Gran Sasso, 67100 L’Aquila, AQ, Italy

14 INAF, Osservatorio Astronomico di Roma, Via di Frascati 33, 00078 Monteporzio Catone, RM, Italy

15 INAF – Osservatorio Astronomico d’Abruzzo, Via M. Maggini snc, 64100 Teramo, Italy

16 DARK, Niels Bohr Institute, University of Copenhagen, Lyngbyvej 2, 2100 Copenhagen Ø, Denmark

17 Nordita, KTH Royal Institute of Technology and Stockholm University, Roslagstullsbacken 23, 106 91 Stockholm, Sweden

18 Instituto de Astrofísica de Andalucía (IAA-CSIC), Glorieta de la Astronomía s/n, 18008 Granada, Spain

19 Institute for Astronomy, University of Hawai’i, 2680 Woodlawn Drive, Honolulu, HI 96822, USA

20 APC, Univ Paris Diderot, CNRS/IN2P3, CEA/Irfu, Obs de Paris, Sorbonne Paris Cité, France

21 AIM, CEA, CNRS, Université Paris-Saclay, Université Paris Diderot, Sorbonne Paris Cité, 91191 Gif-sur-Yvette, France

22 The Oskar Klein Centre, Department of Astronomy, Stockholm University, AlbaNova 10691, Stockholm, Sweden

23 Max-Planck-Institut für Extraterrestrische Physik, Giessenbacherstraße 1, 85748 Garching, Germany

24 INFN, Sezione di Padova, 35131 Padova, Italy

25 Astrophysics Research Institute, IC2 building, Liverpool Science Park, 146 Brownlow Hill, Liverpool L3 5RF, UK

26 Department of Physics, University of Warwick, Coventry CV4 7AL, UK

27 INAF – Istituto di radioastronomia Bologna, Bologna, Italy

28 ASI Science Data Centre, Via del Politecnico snc, 00133 Rome, Italy

29 Faculty of Science, Department of Astronomy and Space Sciences, Istanbul University, Beyazıt, 34119 Istanbul, Turkey

30 Department of Physics and Astronomy, University of Sheffield, Sheffield S3 7RH, UK
31 Instituto de Astrofísica de Canarias, 38205 La Laguna, Tenerife, Spain

32 Institute of Space Sciences (ICE, CSIC), Campus UAB, Carrer de Can Magrans s/n, 08193 Barcelona, Spain

33 School of Physics and Astronomy, University of Leicester, University Road, Leicester LE1 7RH, UK

34 Physics and Astronomy Department Galileo Galilei, University of Padova, Padova, Italy

35 School of Physics, O’Brien Centre for Science North, University College Dublin, Belfield, Dublin 4, Ireland

36 Cosmic Dawn Center (DAWN), Copenhagen Ø, Denmark

37 Niels Bohr Institute, University of Copenhagen, Lyngbyvej 2, 2100 Copenhagen Ø, Denmark

38 Departamento de Física Teórica y del Cosmos, Universidad de Granada, 18071 Granada, Spain

39 Astrophysics Research Centre, School of Mathematics and Physics, Queen's University Belfast, Belfast BT7 1NN, UK

40 University of Nova Gorica, Center for Astrophysics and Cosmology, Vipavska 13, 5000 Nova Gorica, Slovenia

41 Institute of Astronomy, University of Cambridge, Madingley Road, Cambridge CB3 0HA, UK

42 CENTRA-Centro de Astrofísica e Gravitação and Departamento de Física, Instituto Superior Técnico, Universidade de Lisboa, Avenida Rovisco Pais, 1049-001 Lisboa, Portugal

43 Università degli Studi di Urbino “Carlo Bo”, 61029 Urbino, Italy

44 INFN, Sezione di Firenze, 50019 Sesto Fiorentino, Firenze, Italy

45 Astronomical Observatory, University of Warsaw, Al. Ujazdowskie 4, 00-478 Warszawa, Poland

46 Department of Astrophysics/IMAPP, Radboud University, PO Box 9010, 6500 GL Nijmegen, The Netherlands

47 Department of Astronomy, University of Cape Town, Private Bag X3, Rondebosch 7701, South Africa

48 South African Astronomical Observatory, PO Box 9, Observatory 7935, South Africa

49 The Inter-University Institute for Data Intensive Astronomy, University of Cape Town, Private Bag X3, Rondebosch 7701, South Africa

50 Department of Physics and Astronomy, University of Southampton, Southampton SO17 1BJ, UK

51 Department of Physics and Astronomy, University of Turku, Vesilinantie 5, Turku 20014, Finland

52 Centre for Astrophysics and Cosmology, Science Institute, University of Iceland, Dunhagi 5, 107 Reykjavík, Iceland

53 School of Physics & Astronomy, Cardiff University, Queens Buildings, The Parade, Cardiff CF24 3AA, UK

54 Anton Pannekoek Institute for Astronomy, University of Amsterdam, Science Park 904, 1098 XH Amsterdam, The Netherlands

55 Key Laboratory of Dark Matter and Space Astronomy, Purple Mountain Observatory, Chinese Academy of Sciences, Nanjing 210008, PR China

56 SRON, Netherlands Institute for Space Research, Sorbonnelaan 2, 3584 CA Utrecht, The Netherlands

57 School of Physics and Astronomy, The University of Manchester, Manchester M13 9PL, UK

58 Max-Planck-Institut für Astronomie, Königstuhl 17, 69117 Heidelberg, Germany

59 Thüringer Landessternwarte Tautenburg, Sternwarte 5, 07778 Tautenburg, Germany

60 Finnish Centre for Astronomy with ESO (FINCA), 20014 University of Turku, Finland

61 DTU Space, National Space Institute, Technical University of Denmark, Elektrovej 327, 2800 Kongens Lyngby, Denmark

62 Università degli Studi di Trieste and INFN, Sezione di Trieste, 34127 Trieste, Italy

63 School of Physics, Trinity College Dublin, University of Dublin, College Green, Dublin 2, Ireland

64 Institute for Gravitational Wave Astronomy and School of Physics and Astronomy, University of Birmingham, Edgbaston, Birmingham B15 2TT, UK

65 Astronomical Observatory Institute, Faculty of Physics, Adam Mickiewicz University, ul. Słoneczna 36, 60-286 Poznań, Poland

66 Institute of Cosmology and Gravitation, University of Portsmouth, Portsmouth PO1 3FX, UK

67 Istituto di Astrofisica e Planetologia Spaziali (INAF), Via del Fosso del Cavaliere 100, Roma 00133, Italy

68 Università di Pisa, Largo B. Pontecorvo 3, 56127 Pisa, Italy

69 INFN - Sezione di Pisa, Largo B. Pontecorvo 3, 56127 Pisa, Italy

70 Department of Physics, Trento University, Via Sommarive 14, 38123 Povo, Trento, Italy

71 Departamento de Astrofísica, Universidad de La Laguna, 38206 La Laguna, Tenerife, Spain

72 Departamento de Ciencias Físicas, Universidad Andres Bello, Avda. Republica 252, Santiago, Chile

73 National Astronomical Research Institute of Thailand, Ministry of Science and Technology, Chiang Mai 50180, Thailand

74 INAF - Osservatorio Astronomico di Cagliari, Via della Scienza 5, 09047 Selargius, CA, Italy

75 University of Cagliari, Dept of Physics, S.P. Monserrato-Sestu Km 0, 700, 09042 Monserrato, CA, Italy

76 Dipartimento di Fisica e Astronomia “E. Majorana”, Università degli studi di Catania, Via Santa Sofia 64, 95123 Catania, Italy

77 INFN - Laboratori Nazionali del Sud, Via Santa Sofia 62, 95123 Catania, Italy

78 INFN, Sezione di Napoli, Complesso Universitario di Monte S. Angelo, Via Cintia Edificio 6, 80126 Napoli, Italy

79 Department of Physics, Via Cinthia, 80126 Fuorigrotta, Naples, Italy

80 Armagh Observatory and Planetarium, Armagh BT61 9DG, UK

81 Department of Physics and Astronomy, Johns Hopkins University, Baltimore, MD 21218, USA

82 Institut für Astro- und Teilchenphysik, Universität Innsbruck, Technikerstrasse 25/8, 6020 Innsbruck, Austria

83 The Oskar Klein Centre, Department of Physics, Stockholm University, AlbaNova, 10691 Stockholm, Sweden

84 INAF - Istituto di Astrofisica Spaziale e Fisica Cosmica di Milano, Via A. Corti 12, 20133 Milano, Italy

85 Physics Department, University of Calabria, Via P. Bucci, 87036 Rende, Italy

86 European Southern Observatory, Alonso de Córdova, 3107, Vitacura, Santiago 763-0355, Chile

87 Department of Physics, University of Bath, Bath BA2 7AY, UK

88 Department of Physics, Harvard University, Cambridge, MA 02138, USA

89 Department of Physics, The George Washington University, 725 21st Street NW, Washington, DC 20052, USA

90 Astronomy, Physics, and Statistics Institute of Sciences (APSIS), The George Washington University, Washington, DC 20052, USA

91 GEPI, Observatoire de Paris, PSL University, CNRS, 5 Place Jules Janssen, 92190 Meudon, France
## Appendix A: From limiting magnitudes to limits in the model parameter space

For an EM counterpart model defined by an intrinsic light curve $dL/dv(v, t)$ (specific luminosity at a given rest-frame frequency $v$, as a function of rest-frame post-merger time $t$), one can in principle use the galaxy-targeted and wide-field observations to exclude the presence of such emission in a given galaxy. We define here a framework that allows us to combine the results of different searches with a heterogeneous range of telescopes. We work under the simplifying assumption that each observation has a well-defined limiting flux, above which we can exclude a detection with high confidence.

A galaxy-targeted search consists of a set of observations of $N_{gal}$ galaxies, each observed $N_{obs,i}$ times ($i$ here runs on the galaxies). Each observation takes place at a time $t_{i,j}$ post-merger, and reaches a limiting flux $F_{lim,i,j}$ in a band whose central frequency is $v_{i,j}$. We can exclude that galaxy $i$ hosted the putative EM counterpart as long as

$$F_{\text{lim},i,j} < F_{\text{model},i,j} = \frac{(1+z_i)}{4\pi d_{L,i}^2} \frac{dL}{dv} \left( (1+z_i)v_{i,j} \frac{t_{i,j}}{(1+z_i)} \right), \quad (A.1)$$

for any index $j$ running over the $N_{obs,i}$ observations of that galaxy. In that case, the observations contribute a total of $P_{gal,i}$ to the confidence at which the particular EM counterpart model can be excluded. Formally, calling $\xi$ the set of parameters and assumptions that define a particular EM counterpart model, we can exclude $\xi$ with a confidence defined by

$$1 - P(\xi) = \sum_{i=1}^{N_{gal}} P_{gal,i} \mathcal{E}_i(\xi), \quad (A.2)$$

with

$$\mathcal{E}_i(\xi) = 1 - \prod_{j=1}^{N_{obs,i}} H(F_{\text{lim},i,j} - F_{\text{model}}(\xi, v_{i,j}, t_{i,j})), \quad (A.3)$$

where $H$ is the Heaviside step function. The quantity $\mathcal{E}_{i,j}(\xi)$ is 1 if at least one observation of galaxy $i$ is constraining (i.e. it satisfies inequality reported in Eq. (A.1)), and 0 otherwise. In the absence of a detection, the quantity defined by Eq. (A.2) is most commonly referred to as the “covered probability” with respect to a particular source model.

In the case of a wide-field search, one can define the corresponding exclusion confidence as<fnref n="23" />

$$1 - P(\xi) = \sum_{i=1}^{N_{\text{tiles}}} P_{\text{tile},i} \int_{0}^{\infty} dr \frac{dP}{dr} \mathcal{E}_i(r, \xi), \quad (A.4)$$

with

$$\mathcal{E}_i(r, \xi) = 1 - \prod_{j=1}^{N_{\text{obs},i}} H(F_{\text{lim},i,j} - F_{\text{model}}(r, \xi, v_{i,j}, t_{i,j})). \quad (A.5)$$

Here the sum runs over a number $N_{tiles}$ of non-overlapping sky tiles (e.g. a healpix tessellation), each observed $N_{obs,i}$ times, at epochs $t_{i,j}$, with limiting fluxes $F_{lim,i,j}$ in bands whose central frequencies are $v_{i,j}$. $P_{tile,i}$ is the (2D) skymap probability density integrated over the tile, while $dP/dr$ defines how this probability density is distributed over luminosity distance $r \equiv d_L$ (i.e. sky-position-conditional distance probability density of the tile in the 3D skymap – Singer et al. 2016). Here $\mathcal{E}_i(r, \xi)$ equals 1 up to the distance beyond which the putative EM counterpart becomes too faint to be detected by the wide-field observations of tile $i$, and 0 for longer distances. This general framework allows for combining constraints from different wide-field searches.

The results from wide-field and galaxy-targeted searches can be combined conservatively by taking the most constraining between the two for each particular EM counterpart model.

## Appendix B: Data tables

<figure><img src="image_17.png" /><figcaption>Fig. B.1. GW probability versus galaxy-weighted probability in wide-field observations. Comparison between the LALInference GW sky localisation probability density ($P_{2D}$) integrated over the tiles of our wide-field observations and the sum of individual galaxy probabilities $P_{gal}$ (Sect. 3.7) over galaxies that fall in the same tile. The data are reported in Table B.2. Due to the wide field of view of these facilities, the galaxy density averages out, and the probabilities computed in the two ways are very similar. We note that the larger scatter for Pan-STARRS observations is due to the fact that observations are divided into smaller “sky cells”. Grouping sky cells in the same observation would reduce the scatter.</figcaption></figure>

In this appendix, we provide tables which list exhaustively both our galaxy-targeted (Table B.1) and wide-field (Table B.2) observations. Figure B.1 compares two different possible definitions of the probability contained in a wide-field observation, namely the LALInference sky localisation probability integrated over the observed tile, and the sum of the individual galaxy probabilities contained in the same tile, showing that the two are essentially equivalent in the case of wide-field observations.

<footnote n="23">²³ For wide-field observations we assume the localisation probability density to follow the GW 3D skymap, that is, we do not distribute such probability to galaxies as in the galaxy-targeted search case. This relies on the assumption that the galaxy density averages out on scales as large as those probed by wide-field observations. Figure B.1 shows that this is a good approximation in our case, and this removes the uncertainty on the catalogue completeness.</footnote>
<table><caption>**Table B.1.** Galaxies targeted in our search, in order of decreasing, combined, probability ($P_{3D} \times P_{gal}$).</caption><thead><tr><th rowspan="2">Galaxy</th><th colspan="5">Targeted galaxies</th><th colspan="4">Observations of each galaxy</th></tr><tr><th>RA<br>(Deg)</th><th>Dec<br>(Deg)</th><th>d<sub>L</sub><br>(Mpc)</th><th>K<sub>s</sub><br>(Vega)</th><th>log(L<sub>K</sub>)<br>(L<sub>☉</sub>)</th><th>P<sub>3D</sub><br>(10<sup>-5</sup>Mpc<sup>-3</sup>)</th><th>P<sub>gal</sub></th><th>Telescope/<br>Instrument</th><th>ΔT<br>(days)</th><th>Band<br>(AB)</th></tr></thead><tbody><tr><td>198197</td><td>12.091079</td><td>-25.126814</td><td>297.2</td><td>12.4</td><td>11.30</td><td>24.341</td><td>0.02607</td><td>LT/I/O/O<br>LT/I/O/O</td><td>0.20<br>0.21</td><td>r<br>i</td><td>20.3<br>20.3</td></tr><tr><td>ESO474-035</td><td>13.173257</td><td>-25.733852</td><td>271.0</td><td>12.5</td><td>11.15</td><td>31.691</td><td>0.02402</td><td>VLT/HAWKI<br>VLT/HAWKI</td><td>1.37<br>4.51</td><td>K<br>K</td><td>22.3<br>22.0</td></tr><tr><td>797191</td><td>12.178233</td><td>-23.773075</td><td>253.9</td><td>12.3</td><td>11.22</td><td>24.334</td><td>0.02150</td><td>VLT/HAWKI<br>VLT/HAWKI</td><td>1.39<br>4.54</td><td>K<br>K</td><td>22.5<br>21.5</td></tr><tr><td>773232</td><td>12.793847</td><td>-25.954172</td><td>279.4</td><td>12.6</td><td>11.17</td><td>26.130</td><td>0.02098</td><td>VLT/HAWKI<br>VLT/HAWKI</td><td>1.35<br>4.52</td><td>K<br>K</td><td>22.5<br>22.4</td></tr><tr><td>IC1587</td><td>12.180364</td><td>-23.561686</td><td>195.8</td><td>11.4</td><td>11.35</td><td>13.990</td><td>0.01673</td><td>WHT/ACAM<br>VLT/FORS</td><td>0.29<br>1.40</td><td>r<br>i</td><td>22.2<br>24.0</td></tr><tr><td>2798981</td><td>13.803607</td><td>-26.330975</td><td>237.3</td><td>11.7</td><td>11.37</td><td>12.437</td><td>0.01585</td><td>LT/I/O/O<br>VLT/HAWKI</td><td>0.25<br>4.49</td><td>r<br>K</td><td>20.3<br>22.5</td></tr><tr><td>787700</td><td>12.874431</td><td>-24.642494</td><td>274.2</td><td>13.2</td><td>10.91</td><td>34.979</td><td>0.01535</td><td>VLT/HAWKI<br>VLT/HAWKI</td><td>1.42<br>6.32</td><td>K<br>K</td><td>21.4<br>21.0</td></tr><tr><td>2864</td><td>12.256170</td><td>-23.811317</td><td>233.7</td><td>12.6</td><td>11.02</td><td>25.547</td><td>0.01453</td><td>LT/I/O/O<br>LT/I/O/O</td><td>0.22<br>0.23</td><td>r<br>i</td><td>20.3<br>20.3</td></tr><tr><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td>WHT/ACAM<br>VLT/FORS</td><td>0.30<br>1.43</td><td>r<br>i</td><td>21.7<br>23.5</td></tr><tr><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td>WHT/ACAM<br>VLT/HAWKI</td><td>0.28<br>0.28</td><td>r<br>i</td><td>20.7<br>20.7</td></tr><tr><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td>ESO2.2/GROND<br>ESO2.2/GROND</td><td>0.44<br>0.44</td><td>K<br>r</td><td>18.5<br>21.0</td></tr><tr><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td>ESO2.2/GROND<br>VLT/HAWKI</td><td>1.48<br>4.52</td><td>K<br>K</td><td>22.4<br>22.0</td></tr></tbody></table>

**Notes.** For each galaxy we report: HyperLEDA name (except for those followed by 2MX in parentheses, which are from 2MASSX), right ascension, declination and luminosity distance as listed in the GLADE catalogue*, K<sub>s</sub>-band magnitude as measured in our own VISTA/VINROUGE data. K<sub>s</sub>-band luminosity, 3D skymap probability density per unit volume at the galaxy position, and localisation probability (Eq. (1)), followed by a list of observations of that galaxy, for which we specify the telescope, the post-merger time at which the observation was performed, the band, and the limiting magnitude. *The distances in GLADE are computed assuming a flat cosmology with $H_0 = 70\,\text{km s}^{-1}\,\text{Mpc}^{-1}$ and $\Omega_m = 0.27$ (Dálya et al. 2018). The difference in luminosity distance when compared to our assumed cosmology ($H_0 = 70\,\text{km s}^{-1}\,\text{Mpc}^{-1}$ and $\Omega_m = 0.30$) is <1 Mpc within 400 Mpc, but we anyway report the distance computed with our assumed cosmology.
<table><caption>Table B.1. continued.</caption><thead><tr><th rowspan="2">Galaxy</th><th colspan="4">Targeted galaxies</th><th rowspan="2">K<sub>s</sub><br>(Vega)</th><th rowspan="2">log(L<sub>K</sub>)<br>(L<sub>☉</sub>)</th><th rowspan="2">P<sub>3D</sub><br>(10<sup>-5</sup>Mpc<sup>-3</sup>)</th><th rowspan="2">P<sub>gal</sub></th><th colspan="4">Observations of each galaxy</th></tr><tr><th>RA<br>(Deg)</th><th>Dec<br>(Deg)</th><th>d<sub>L</sub><br>(Mpc)</th><th>P</th><th>Telescope/<br>Instrument</th><th>ΔT<br>(days)</th><th>Band</th><th>m<sub>lim</sub><br>(AB)</th></tr></thead><tbody><tr><td>2694</td><td>11.543399</td><td>-24.650192</td><td>219.8</td><td>12.1</td><td>11.17</td><td>18.183</td><td>0.01448</td><td></td><td>LT/IO:O</td><td>0.24</td><td>r</td><td>20.3</td></tr><tr><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td>ESO2.2/GROND</td><td>0.50</td><td>K</td><td>18.0</td></tr><tr><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td>ESO2.2/GROND</td><td>0.50</td><td>r</td><td>20.0</td></tr><tr><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td>VLT/HAWKI</td><td>1.46</td><td>K</td><td>21.8</td></tr><tr><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td>VLT/HAWKI</td><td>6.34</td><td>K</td><td>21.0</td></tr><tr><td>769203</td><td>12.800914</td><td>-26.313063</td><td>262.8</td><td>12.6</td><td>11.09</td><td>20.970</td><td>0.01374</td><td></td><td>WHT/ACAM</td><td>0.29</td><td>r</td><td>22.5</td></tr><tr><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td>VLT/FORS</td><td>1.41</td><td>i</td><td>23.0</td></tr><tr><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td>WHT/LIRIS</td><td>3.27</td><td>K</td><td>20.0</td></tr><tr><td>133717</td><td>12.320091</td><td>-26.219179</td><td>330.6</td><td>12.1</td><td>11.52</td><td>7.504</td><td>0.01329</td><td></td><td>LT/IO:O</td><td>0.33</td><td>r</td><td>20.3</td></tr><tr><td>198196</td><td>11.870618</td><td>-25.440655</td><td>265.8</td><td>12.9</td><td>11.01</td><td>23.869</td><td>0.01301</td><td></td><td>ESO2.2/GROND</td><td>0.47</td><td>K</td><td>17.5</td></tr><tr><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td>ESO2.2/GROND</td><td>0.47</td><td>r</td><td>21.0</td></tr><tr><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td>VLT/HAWKI</td><td>1.41</td><td>K</td><td>22.5</td></tr><tr><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td>VLT/HAWKI</td><td>6.31</td><td>K</td><td>22.5</td></tr><tr><td>787067</td><td>13.064037</td><td>-24.698730</td><td>224.5</td><td>12.8</td><td>10.91</td><td>29.229</td><td>0.01289</td><td></td><td>WHT/ACAM</td><td>0.30</td><td>r</td><td>21.7</td></tr><tr><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td>WHT/LIRIS</td><td>2.32</td><td>K</td><td>20.0</td></tr><tr><td>3235517</td><td>12.103154</td><td>-25.595707</td><td>272.1</td><td>13.1</td><td>10.98</td><td>24.497</td><td>0.01249</td><td></td><td>WHT/LIRIS</td><td>2.29</td><td>K</td><td>20.0</td></tr><tr><td>798968</td><td>12.643789</td><td>-23.618547</td><td>227.8</td><td>12.4</td><td>11.07</td><td>19.656</td><td>0.01234</td><td></td><td>WHT/LIRIS</td><td>1.33</td><td>K</td><td>20.0</td></tr><tr><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td>VLT/FORS</td><td>1.34</td><td>i</td><td>23.5</td></tr><tr><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td>WHT/LIRIS</td><td>3.17</td><td>K</td><td>20.0</td></tr><tr><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td>WHT/LIRIS</td><td>3.17</td><td>K</td><td>20.0</td></tr><tr><td>2998</td><td>12.828166</td><td>-26.168060</td><td>285.2</td><td>13.0</td><td>11.02</td><td>21.414</td><td>0.01195</td><td></td><td>VLT/FORS</td><td>1.42</td><td>r</td><td>23.0</td></tr><tr><td>IC1588</td><td>12.740481</td><td>-23.557995</td><td>240.6</td><td>12.4</td><td>11.11</td><td>17.184</td><td>0.01187</td><td></td><td>VLT/FORS</td><td>1.34</td><td>i</td><td>24.0</td></tr><tr><td>00485495-2504100 (2MX)</td><td>12.228970</td><td>-25.069471</td><td>240.6</td><td>13.1</td><td>10.84</td><td>31.409</td><td>0.01165</td><td></td><td>WHT/LIRIS</td><td>3.20</td><td>K</td><td>20.0</td></tr><tr><td>3235434</td><td>13.788131</td><td>-25.455704</td><td>229.5</td><td>12.6</td><td>11.01</td><td>20.751</td><td>0.01151</td><td></td><td>TNG/NICS</td><td>1.22</td><td>K</td><td>19.0</td></tr><tr><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td>VLT/FORS</td><td>1.46</td><td>i</td><td>24.0</td></tr><tr><td>783013</td><td>13.656960</td><td>-25.067099</td><td>221.7</td><td>12.5</td><td>11.00</td><td>20.588</td><td>0.01117</td><td></td><td>TNG/NICS</td><td>5.23</td><td>K</td><td>19.0</td></tr><tr><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td>WHT/ACAM</td><td>0.31</td><td>r</td><td>22.2</td></tr><tr><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td>TNG/NICS</td><td>1.25</td><td>K</td><td>19.0</td></tr><tr><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td>TNG/NICS</td><td>1.25</td><td>K</td><td>19.0</td></tr><tr><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td>WHT/LIRIS</td><td>1.26</td><td>K</td><td>20.0</td></tr><tr><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td>VLT/FORS</td><td>1.47</td><td>i</td><td>24.0</td></tr><tr><td>768450.0</td><td>13.722800</td><td>-26.381190</td><td>259.0</td><td>12.3</td><td>11.22</td><td>12.416</td><td>0.01100</td><td></td><td>VLT/HAWKI</td><td>1.42</td><td>K</td><td>21.4</td></tr><tr><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td>VLT/HAWKI</td><td>1.42</td><td>K</td><td>21.8</td></tr><tr><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td>VLT/HAWKI</td><td>6.32</td><td>K</td><td>21.0</td></tr></tbody></table>
<table><thead><tr><th rowspan="2">Galaxy</th><th colspan="5">Targeted galaxies</th><th colspan="5">Observations of each galaxy</th></tr><tr><th>RA<br>(Deg)</th><th>Dec<br>(Deg)</th><th>d<sub>L</sub><br>(Mpc)</th><th>K<sub>s</sub><br>(Vega)</th><th>log(L<sub>K</sub>)<br>(L<sub>☉</sub>)</th><th>P<sup>3D</sup><br>(10<sup>-5</sup>Mpc<sup>-3</sup>)</th><th>P<sub>gal</sub></th><th>Telescope/<br>Instrument</th><th>ΔT<br>(days)</th><th>Band</th><th>m<sub>lim</sub><br>(AB)</th></tr></thead><tbody><tr><td>133715</td><td>12.441799</td><td>-26.443037</td><td>242.2</td><td>12.6</td><td>11.04</td><td>16.618</td><td>0.00975</td><td>LT/IO-O</td><td>0.34</td><td>r</td><td>20.3</td></tr><tr><td>3231</td><td>13.704722</td><td>-26.371256</td><td>236.5</td><td>12.3</td><td>11.14</td><td>13.226</td><td>0.00974</td><td>ESO2.2/GROND</td><td>0.45</td><td>K</td><td>18.5</td></tr><tr><td>00530427-2610148 (2MX)</td><td>13.267828</td><td>-26.170792</td><td>243.0</td><td>13.1</td><td>10.87</td><td>23.859</td><td>0.00955</td><td>VLT/HAWKI</td><td>1.44</td><td>K</td><td>22.5</td></tr><tr><td>78830</td><td>13.490571</td><td>-24.543142</td><td>295.8</td><td>13.1</td><td>11.04</td><td>15.322</td><td>0.00905</td><td>VLT/HAWKI</td><td>6.36</td><td>K</td><td>22.5</td></tr><tr><td>786964</td><td>12.473102</td><td>-24.707197</td><td>233.1</td><td>13.4</td><td>10.70</td><td>30.588</td><td>0.00820</td><td>WHT/ACAM</td><td>0.31</td><td>r</td><td>21.4</td></tr><tr><td>786999</td><td>13.269273</td><td>-24.704401</td><td>231.0</td><td>13.4</td><td>10.71</td><td>28.883</td><td>0.00800</td><td>WHT/ACAM</td><td>0.31</td><td>r</td><td>22.0</td></tr><tr><td>792107</td><td>11.771893</td><td>-24.238703</td><td>291.8</td><td>13.3</td><td>10.94</td><td>17.133</td><td>0.00797</td><td>VLT/FORS</td><td>1.39</td><td>i</td><td>24.0</td></tr><tr><td>133716</td><td>12.364808</td><td>-26.538301</td><td>223.9</td><td>12.3</td><td>11.11</td><td>11.189</td><td>0.00778</td><td>WHT/ACAM</td><td>0.28</td><td>r</td><td>21.7</td></tr><tr><td>133702</td><td>12.242800</td><td>-25.693450</td><td>296.0</td><td>13.7</td><td>10.78</td><td>20.375</td><td>0.00658</td><td>WHT/LIRIS</td><td>1.22</td><td>K</td><td>20.0</td></tr><tr><td>100480</td><td>23.594837</td><td>-32.835316</td><td>286.2</td><td>11.3</td><td>11.72</td><td>2.269</td><td>0.00648</td><td>TNG/NICS</td><td>1.27</td><td>K</td><td>20.0</td></tr><tr><td>133703</td><td>11.472051</td><td>-23.772461</td><td>234.2</td><td>12.6</td><td>11.00</td><td>11.794</td><td>0.00640</td><td>VLT/FORS</td><td>1.35</td><td>i</td><td>24.0</td></tr><tr><td>3264</td><td>13.806392</td><td>-26.321253</td><td>184.1</td><td>11.7</td><td>11.16</td><td>8.011</td><td>0.00623</td><td>TNG/NICS</td><td>5.30</td><td>K</td><td>20.0</td></tr><tr><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td>LT/IO-O</td><td>0.29</td><td>r</td><td>20.3</td></tr><tr><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td>VLT/FORS</td><td>1.49</td><td>i</td><td>24.5</td></tr><tr><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td>ESO2.2/GROND</td><td>0.37</td><td>K</td><td>17.5</td></tr><tr><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td>ESO2.2/GROND</td><td>0.37</td><td>r</td><td>21.0</td></tr><tr><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td>VLT/HAWKI</td><td>1.44</td><td>K</td><td>22.8</td></tr><tr><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td>VLT/HAWKI</td><td>6.36</td><td>K</td><td>22.0</td></tr><tr><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td>TNG/NICS</td><td>1.23</td><td>K</td><td>20.0</td></tr><tr><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td>VLT/FORS</td><td>1.48</td><td>i</td><td>24.5</td></tr><tr><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td>TNG/NICS</td><td>5.26</td><td>K</td><td>20.0</td></tr><tr><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td>ESO2.2/GROND</td><td>0.34</td><td>K</td><td>19.0</td></tr><tr><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td>ESO2.2/GROND</td><td>0.34</td><td>r</td><td>20.5</td></tr><tr><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td>LT/IO-O</td><td>0.21</td><td>r</td><td>20.3</td></tr><tr><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td>LT/IO-O</td><td>0.22</td><td>i</td><td>20.3</td></tr><tr><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td>ESO2.2/GROND</td><td>0.49</td><td>K</td><td>19.0</td></tr><tr><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td>ESO2.2/GROND</td><td>0.49</td><td>r</td><td>21.0</td></tr><tr><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td>WHT/LIRIS</td><td>2.19</td><td>K</td><td>20.0</td></tr><tr><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td>ESO2.2/GROND</td><td>0.36</td><td>K</td><td>17.5</td></tr><tr><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td>ESO2.2/GROND</td><td>0.36</td><td>r</td><td>21.0</td></tr><tr><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td>VLT/HAWKI</td><td>1.42</td><td>K</td><td>21.4</td></tr><tr><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td>VLT/HAWKI</td><td>6.32</td><td>K</td><td>21.0</td></tr></tbody></table>
<table><caption>Table B.1. continued.</caption><thead><tr><th rowspan="2">Galaxy</th><th colspan="5">Targeted galaxies</th><th colspan="5">Observations of each galaxy</th></tr><tr><th>RA<br>(Deg)</th><th>Dec<br>(Deg)</th><th>d<sub>L</sub><br>(Mpc)</th><th>K<sub>s</sub><br>(Vega)</th><th>log(L<sub>K</sub>)<br>(L<sub>☉</sub>)</th><th>P<sub>3D</sub><br>(10<sup>-5</sup>Mpc<sup>-3</sup>)</th><th>P<sub>gal</sub></th><th>Telescope/<br>Instrument</th><th>ΔT<br>(days)</th><th>Band</th><th>m<sub>lim</sub><br>(AB)</th></tr></thead><tbody><tr><td>101138</td><td>13.804421</td><td>-24.044033</td><td>203.1</td><td>11.8</td><td>11.21</td><td>7.085</td><td>0.00613</td><td>ESO2.2/GROND</td><td>0.52</td><td>K</td><td>19.5</td></tr><tr><td>771842</td><td>12.881038</td><td>-26.077213</td><td>294.6</td><td>14.0</td><td>10.69</td><td>19.981</td><td>0.00530</td><td>ESO2.2/GROND</td><td>0.52</td><td>r</td><td>22.0</td></tr><tr><td>198221</td><td>12.385612</td><td>-26.538588</td><td>301.2</td><td>13.3</td><td>10.95</td><td>9.862</td><td>0.00473</td><td>VLT/FORS</td><td>1.42</td><td>z</td><td>24.5</td></tr><tr><td>ESO474-036</td><td>13.192388</td><td>-22.975018</td><td>213.1</td><td>11.3</td><td>11.45</td><td>3.038</td><td>0.00459</td><td>VLT/HAWKI</td><td>1.44</td><td>K</td><td>22.3</td></tr><tr><td>2947</td><td>12.666141</td><td>-26.813395</td><td>316.8</td><td>12.7</td><td>11.22</td><td>5.106</td><td>0.00455</td><td>VLT/HAWKI</td><td>6.36</td><td>K</td><td>22.5</td></tr><tr><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td>ESO2.2/GROND</td><td>0.33</td><td>K</td><td>19.0</td></tr><tr><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td>ESO2.2/GROND</td><td>0.33</td><td>r</td><td>22.0</td></tr><tr><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td>WHT/ACAM</td><td>0.28</td><td>r</td><td>21.2</td></tr><tr><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td>ESO2.2/GROND</td><td>0.51</td><td>K</td><td>17.5</td></tr><tr><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td>ESO2.2/GROND</td><td>0.51</td><td>r</td><td>20.0</td></tr><tr><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td>TNG/NICS</td><td>1.29</td><td>K</td><td>20.0</td></tr><tr><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td>VLT/FORS</td><td>1.37</td><td>i</td><td>24.5</td></tr><tr><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td>TNG/NICS</td><td>5.33</td><td>K</td><td>20.0</td></tr><tr><td>198225</td><td>12.174922</td><td>-23.368631</td><td>329.4</td><td>12.9</td><td>11.17</td><td>5.247</td><td>0.00418</td><td>WHT/ACAM</td><td>0.27</td><td>r</td><td>22.0</td></tr><tr><td>802484</td><td>12.810972</td><td>-23.291632</td><td>268.2</td><td>13.2</td><td>10.91</td><td>8.806</td><td>0.00388</td><td>WHT/ACAM</td><td>0.30</td><td>r</td><td>22.0</td></tr><tr><td>804101</td><td>11.210597</td><td>-23.159029</td><td>222.7</td><td>12.5</td><td>11.00</td><td>6.527</td><td>0.00353</td><td>LT/IO:O</td><td>0.34</td><td>r</td><td>20.3</td></tr><tr><td>198205</td><td>12.540436</td><td>-23.280048</td><td>256.8</td><td>13.6</td><td>10.71</td><td>12.082</td><td>0.00337</td><td>WHT/LIRIS</td><td>3.23</td><td>K</td><td>20.0</td></tr><tr><td>133698</td><td>14.254662</td><td>-23.837297</td><td>221.6</td><td>12.2</td><td>11.11</td><td>4.664</td><td>0.00324</td><td>WHT/ACAM</td><td>0.32</td><td>r</td><td>22.5</td></tr><tr><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td>ESO2.2/GROND</td><td>0.38</td><td>K</td><td>19.0</td></tr><tr><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td>ESO2.2/GROND</td><td>0.38</td><td>r</td><td>21.0</td></tr><tr><td>198217</td><td>12.346334</td><td>-26.507498</td><td>298.9</td><td>13.8</td><td>10.74</td><td>10.667</td><td>0.00314</td><td>VLT/HAWKI</td><td>1.44</td><td>K</td><td>22.3</td></tr><tr><td>142558</td><td>12.332034</td><td>-26.476397</td><td>253.3</td><td>13.8</td><td>10.59</td><td>14.724</td><td>0.00305</td><td>VLT/HAWKI</td><td>6.36</td><td>K</td><td>22.0</td></tr><tr><td>101147</td><td>22.938530</td><td>-32.415623</td><td>240.8</td><td>12.0</td><td>11.27</td><td>2.868</td><td>0.00289</td><td>VLT/HAWKI</td><td>1.44</td><td>K</td><td>22.8</td></tr><tr><td>798818</td><td>12.726862</td><td>-23.631887</td><td>316.2</td><td>13.9</td><td>10.76</td><td>9.227</td><td>0.00289</td><td>ESO2.2/GROND</td><td>0.44</td><td>K</td><td>19.0</td></tr><tr><td>773323</td><td>12.467500</td><td>-25.946250</td><td>305.6</td><td>14.4</td><td>10.52</td><td>15.905</td><td>0.00287</td><td>ESO2.2/GROND</td><td>0.44</td><td>r</td><td>21.0</td></tr><tr><td>805757</td><td>11.720657</td><td>-23.033381</td><td>249.3</td><td>13.6</td><td>10.66</td><td>10.878</td><td>0.00270</td><td>VLT/LIRIS</td><td>3.30</td><td>K</td><td>23.0</td></tr><tr><td>198201</td><td>12.637869</td><td>-23.295488</td><td>237.0</td><td>13.6</td><td>10.64</td><td>11.357</td><td>0.00266</td><td>LT/IO:O</td><td>0.20</td><td>r</td><td>20.3</td></tr><tr><td>795185</td><td>10.888593</td><td>-23.952913</td><td>307.3</td><td>13.2</td><td>11.02</td><td>4.174</td><td>0.00235</td><td>LT/IO:O</td><td>0.32</td><td>r</td><td>20.3</td></tr><tr><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td>WHT/ACAM</td><td>0.27</td><td>r</td><td>22.5</td></tr><tr><td>131961</td><td>23.018082</td><td>-33.802876</td><td>307.5</td><td>11.9</td><td>11.55</td><td>1.089</td><td>0.00208</td><td>WHT/LIRIS</td><td>2.22</td><td>K</td><td>20.0</td></tr><tr><td>133690</td><td>10.639115</td><td>-22.866308</td><td>292.3</td><td>11.8</td><td>11.54</td><td>1.055</td><td>0.00199</td><td>ESO2.2/GROND</td><td>0.35</td><td>K</td><td>19.0</td></tr><tr><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td>ESO2.2/GROND</td><td>0.35</td><td>r</td><td>22.0</td></tr><tr><td>773149</td><td>12.815700</td><td>-25.960900</td><td>300.1</td><td>15.1</td><td>10.24</td><td>20.028</td><td>0.00189</td><td>WHT/ACAM</td><td>0.26</td><td>r</td><td>21.7</td></tr><tr><td>PGC002953</td><td>12.679955</td><td>-22.888277</td><td>221.5</td><td>12.9</td><td>10.84</td><td>4.890</td><td>0.00183</td><td>ESO2.2/GROND</td><td>0.43</td><td>K</td><td>19.0</td></tr><tr><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td>ESO2.2/GROND</td><td>0.43</td><td>r</td><td>20.0</td></tr><tr><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td>VLT/FORS</td><td>1.40</td><td>i</td><td>24.5</td></tr><tr><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td>WHT/ACAM</td><td>0.29</td><td>r</td><td>21.0</td></tr></tbody></table>
<table><caption>Table B.1. continued.</caption><thead><tr><th rowspan="2">Galaxy</th><th colspan="6">Targeted galaxies</th><th colspan="5">Observations of each galaxy</th></tr><tr><th>RA<br>(Deg)</th><th>Dec<br>(Deg)</th><th>d<sub>L</sub><br>(Mpc)</th><th>K<sub>s</sub><br>(Vega)</th><th>log(L<sub>k</sub>)<br>(L<sub>☉</sub>)</th><th>P<sub>3d</sub><br>(10<sup>-5</sup>Mpc<sup>-3</sup>)</th><th>P<sub>gal</sub></th><th>Telescope/<br>Instrument</th><th>ΔT<br>(days)</th><th>Band</th><th>m<sub>lim</sub><br>(AB)</th></tr></thead><tbody><tr><td>131965</td><td>22.342892</td><td>-32.849354</td><td>228.0</td><td>12.1</td><td>11.19</td><td>1.922</td><td>0.00159</td><td>ESO2.2/GROND</td><td>0.41</td><td>K</td><td>19.0</td></tr><tr><td>IC1581</td><td>11.442945</td><td>-25.920193</td><td>222.0</td><td>13.0</td><td>10.81</td><td>4.485</td><td>0.00155</td><td>ESO2.2/GROND</td><td>0.41</td><td>r</td><td>21.0</td></tr><tr><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td>LT/IO:O</td><td>0.29</td><td>r</td><td>20.3</td></tr><tr><td>684102</td><td>22.549915</td><td>-32.903667</td><td>292.0</td><td>12.8</td><td>11.11</td><td>1.651</td><td>0.00115</td><td>ESO2.2/GROND</td><td>0.42</td><td>K</td><td>18.0</td></tr><tr><td>131958</td><td>24.334146</td><td>-33.330853</td><td>232.2</td><td>12.6</td><td>11.01</td><td>2.003</td><td>0.00109</td><td>ESO2.2/GROND</td><td>0.48</td><td>r</td><td>21.0</td></tr><tr><td>777373.0</td><td>12.718400</td><td>-25.577060</td><td>225.6</td><td>16.0</td><td>9.61</td><td>27.183</td><td>0.00060</td><td>ESO2.2/GROND</td><td>0.38</td><td>K</td><td>19.0</td></tr><tr><td>ESO540-029</td><td>12.120861</td><td>-22.147444</td><td>220.1</td><td>13.4</td><td>10.67</td><td>1.513</td><td>0.00038</td><td>ESO2.2/GROND</td><td>0.38</td><td>r</td><td>21.0</td></tr><tr><td>ESO540-013</td><td>10.242757</td><td>-22.083128</td><td>286.5</td><td>12.5</td><td>11.24</td><td>0.398</td><td>0.00037</td><td>ESO2.2/GROND</td><td>0.35</td><td>r</td><td>20.3</td></tr><tr><td>95210</td><td>22.145348</td><td>-32.296467</td><td>288.3</td><td>14.0</td><td>10.62</td><td>1.587</td><td>0.00036</td><td>ESO2.2/GROND</td><td>0.55</td><td>K</td><td>19.0</td></tr><tr><td>131959</td><td>23.877426</td><td>-32.970612</td><td>288.2</td><td>14.4</td><td>10.49</td><td>2.018</td><td>0.00034</td><td>ESO2.2/GROND</td><td>0.55</td><td>r</td><td>21.0</td></tr><tr><td>ESO474-026</td><td>11.781363</td><td>-24.370647</td><td>114.9</td><td>11.1</td><td>10.98</td><td>0.516</td><td>0.00027</td><td>ESO2.2/GROND</td><td>0.43</td><td>K</td><td>19.5</td></tr><tr><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td>LT/IO:O</td><td>0.24</td><td>r</td><td>20.0</td></tr><tr><td>172818</td><td>11.202063</td><td>-22.116533</td><td>288.6</td><td>13.8</td><td>10.72</td><td>0.713</td><td>0.00020</td><td>ESO2.2/GROND</td><td>0.33</td><td>K</td><td>18.5</td></tr><tr><td>2561</td><td>10.676272</td><td>-21.774082</td><td>148.9</td><td>13.8</td><td>10.14</td><td>0.274</td><td>0.00002</td><td>ESO2.2/GROND</td><td>0.33</td><td>r</td><td>20.0</td></tr><tr><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td>VLT/HAWKI</td><td>1.34</td><td>K</td><td>22.4</td></tr><tr><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td>VLT/HAWKI</td><td>6.28</td><td>K</td><td>22.5</td></tr><tr><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td>WHT/ACAM</td><td>0.27</td><td>r</td><td>21.2</td></tr><tr><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td>LT/IO:O</td><td>0.19</td><td>r</td><td>20.3</td></tr></tbody></table>
<table><caption>Table B.2. Wide-field observations.</caption><thead><tr><th>Facility</th><th>Tile ID</th><th>∫<sub>tile</sub> P<sub>2D</sub> dΩ<br/>[%]</th><th>Σ<sub>tile</sub> P<sub>gal</sub><br/>[%]</th><th>RA<br/>[deg]</th><th>Dec<br/>[deg]</th><th>MJD</th><th>Filter</th><th>Limiting mag<br/>[AB]</th></tr></thead><tbody><tr><td rowspan="10">VST</td><td rowspan="10">0</td><td rowspan="10">1.692</td><td rowspan="10">1.578</td><td rowspan="10">10.573625</td><td rowspan="10">-23.734472</td><td>58711.36</td><td>r</td><td>21.90</td></tr><tr><td>58715.36</td><td>r</td><td>21.60</td></tr><tr><td>58717.36</td><td>r</td><td>22.00</td></tr><tr><td>58724.36</td><td>r</td><td>22.10</td></tr><tr><td>58711.36</td><td>r</td><td>20.90</td></tr><tr><td>58715.36</td><td>r</td><td>21.50</td></tr><tr><td>58717.36</td><td>r</td><td>22.10</td></tr><tr><td>58724.36</td><td>r</td><td>22.20</td></tr><tr><td>58711.36</td><td>r</td><td>22.00</td></tr><tr><td>58715.36</td><td>r</td><td>22.10</td></tr><tr><td>58717.36</td><td>r</td><td>21.80</td></tr><tr><td>58724.36</td><td>r</td><td>21.90</td></tr><tr><td>58710.36</td><td>r</td><td>20.10</td></tr><tr><td>58710.36</td><td>r</td><td>20.60</td></tr><tr><td>58710.36</td><td>r</td><td>21.90</td></tr><tr><td>58711.36</td><td>r</td><td>21.80</td></tr><tr><td>58715.36</td><td>r</td><td>22.10</td></tr><tr><td>58717.36</td><td>r</td><td>22.10</td></tr><tr><td>58724.36</td><td>r</td><td>22.30</td></tr><tr><td>58710.36</td><td>r</td><td>20.70</td></tr><tr><td>58711.36</td><td>r</td><td>21.00</td></tr><tr><td>58714.36</td><td>r</td><td>21.70</td></tr><tr><td>58717.36</td><td>r</td><td>21.90</td></tr><tr><td>58725.36</td><td>r</td><td>22.00</td></tr><tr><td>58710.36</td><td>r</td><td>20.60</td></tr><tr><td>58711.36</td><td>r</td><td>20.90</td></tr><tr><td>58714.36</td><td>r</td><td>21.60</td></tr><tr><td>58717.36</td><td>r</td><td>22.10</td></tr><tr><td>58725.36</td><td>r</td><td>22.20</td></tr><tr><td>58711.36</td><td>r</td><td>22.20</td></tr><tr><td>58714.36</td><td>r</td><td>21.50</td></tr><tr><td>58717.36</td><td>r</td><td>22.10</td></tr><tr><td>58725.36</td><td>r</td><td>22.40</td></tr><tr><td>58710.36</td><td>r</td><td>21.50</td></tr><tr><td>58711.36</td><td>r</td><td>21.80</td></tr><tr><td>58715.36</td><td>r</td><td>22.10</td></tr><tr><td>58717.36</td><td>r</td><td>21.20</td></tr><tr><td>58725.36</td><td>r</td><td>22.00</td></tr><tr><td>58710.36</td><td>r</td><td>20.20</td></tr><tr><td>58710.36</td><td>r</td><td>21.00</td></tr><tr><td>58711.36</td><td>r</td><td>22.20</td></tr><tr><td>58715.36</td><td>r</td><td>22.20</td></tr><tr><td>58717.36</td><td>r</td><td>22.30</td></tr><tr><td>58724.36</td><td>r</td><td>22.20</td></tr><tr><td>58715.36</td><td>r</td><td>21.30</td></tr><tr><td>58717.36</td><td>r</td><td>21.60</td></tr><tr><td>58725.36</td><td>r</td><td>22.20</td></tr><tr><td>58710.36</td><td>r</td><td>21.40</td></tr><tr><td>58711.36</td><td>r</td><td>21.90</td></tr><tr><td>58714.36</td><td>r</td><td>21.70</td></tr><tr><td>58717.36</td><td>r</td><td>21.80</td></tr><tr><td>58725.36</td><td>r</td><td>22.00</td></tr><tr><td>58710.36</td><td>r</td><td>21.00</td></tr><tr><td>58711.36</td><td>r</td><td>21.80</td></tr><tr><td>58714.36</td><td>r</td><td>21.70</td></tr><tr><td>58717.36</td><td>r</td><td>21.90</td></tr><tr><td>58725.36</td><td>r</td><td>21.80</td></tr><tr><td>58711.36</td><td>r</td><td>21.80</td></tr><tr><td>58714.36</td><td>r</td><td>21.40</td></tr><tr><td>58717.36</td><td>r</td><td>21.70</td></tr><tr><td>58725.36</td><td>r</td><td>21.90</td></tr><tr><td>58711.36</td><td>r</td><td>22.00</td></tr><tr><td>58714.36</td><td>r</td><td>21.90</td></tr><tr><td>58717.36</td><td>r</td><td>21.40</td></tr><tr><td>58725.36</td><td>r</td><td>22.10</td></tr></tbody></table>

Notes. For each facility we report a list of unique tiles, identified by a Tile ID and by the RA and Dec of the tile centre, along with the LALInference sky location probability density ($P_{2D}$) integrated over the solid angle subtended by the tile, and the sum of the galaxy probabilities $P_{gal}$ over galaxies that fall within the tile (see Fig. B.1 for a comparison plot). For each tile, we then report the MJD, Filter and Limiting magnitude of the observations carried out of these tiles.
<table><caption>**Table B.2.** continued.</caption><thead><tr><th>Facility</th><th>Tile ID</th><th>∫<sub>tile</sub> P<sub>2D</sub> dΩ [%]</th><th>Σ<sub>tile</sub> P<sub>gal</sub> [%]</th><th>RA [deg]</th><th>Dec [deg]</th><th>MJD</th><th>Filter</th><th>Limiting mag [AB]</th></tr></thead><tbody><tr><td rowspan="12">VST</td><td rowspan="12">19</td><td rowspan="12">0.936</td><td rowspan="12">2.050</td><td rowspan="12">13.441833</td><td rowspan="12">-22.736889</td><td>58711.36</td><td>r</td><td>21.90</td></tr><tr><td>58715.36</td><td>r</td><td>21.90</td></tr><tr><td>58717.36</td><td>r</td><td>22.00</td></tr><tr><td>58725.36</td><td>r</td><td>22.00</td></tr><tr><td>58710.36</td><td>r</td><td>21.50</td></tr><tr><td>58715.36</td><td>r</td><td>21.70</td></tr><tr><td>58717.36</td><td>r</td><td>21.70</td></tr><tr><td>58725.36</td><td>r</td><td>21.80</td></tr><tr><td>58710.36</td><td>r</td><td>21.10</td></tr><tr><td>58715.36</td><td>r</td><td>21.90</td></tr><tr><td>58717.36</td><td>r</td><td>21.60</td></tr><tr><td>58725.36</td><td>r</td><td>22.10</td></tr><tr><td>58711.36</td><td>r</td><td>21.80</td></tr><tr><td>58714.36</td><td>r</td><td>21.40</td></tr><tr><td>58717.36</td><td>r</td><td>21.70</td></tr><tr><td>58725.36</td><td>r</td><td>22.10</td></tr><tr><td>58710.36</td><td>r</td><td>19.90</td></tr><tr><td>58711.36</td><td>r</td><td>20.90</td></tr><tr><td>58714.36</td><td>r</td><td>21.90</td></tr><tr><td>58717.36</td><td>r</td><td>22.00</td></tr><tr><td>58725.36</td><td>r</td><td>22.10</td></tr><tr><td>58711.36</td><td>r</td><td>21.90</td></tr><tr><td>58714.36</td><td>r</td><td>21.40</td></tr><tr><td>58717.36</td><td>r</td><td>21.80</td></tr><tr><td>58725.36</td><td>r</td><td>21.90</td></tr><tr><td>58711.36</td><td>r</td><td>21.80</td></tr><tr><td>58714.36</td><td>r</td><td>21.80</td></tr><tr><td>58717.36</td><td>r</td><td>21.70</td></tr><tr><td>58725.36</td><td>r</td><td>22.10</td></tr><tr><td>58714.36</td><td>r</td><td>21.80</td></tr><tr><td>58714.36</td><td>r</td><td>21.80</td></tr><tr><td>58717.36</td><td>r</td><td>21.90</td></tr><tr><td>58725.36</td><td>r</td><td>22.00</td></tr><tr><td>58711.36</td><td>r</td><td>22.00</td></tr><tr><td>58714.36</td><td>r</td><td>21.30</td></tr><tr><td>58717.36</td><td>r</td><td>21.70</td></tr><tr><td>58725.36</td><td>r</td><td>22.00</td></tr><tr><td>58710.36</td><td>r</td><td>20.80</td></tr><tr><td>58710.36</td><td>r</td><td>21.90</td></tr><tr><td>58712.32</td><td>K</td><td>21.00</td></tr><tr><td>58719.37</td><td>K</td><td>21.00</td></tr><tr><td>58750.36</td><td>K</td><td>21.00</td></tr><tr><td>58712.28</td><td>K</td><td>21.00</td></tr><tr><td>58719.33</td><td>K</td><td>21.00</td></tr><tr><td>58751.22</td><td>K</td><td>21.00</td></tr><tr><td>58712.36</td><td>K</td><td>21.00</td></tr><tr><td>58719.41</td><td>K</td><td>21.00</td></tr><tr><td>58751.34</td><td>K</td><td>21.00</td></tr><tr><td>58712.31</td><td>K</td><td>21.00</td></tr><tr><td>58719.36</td><td>K</td><td>21.00</td></tr><tr><td>58751.23</td><td>K</td><td>21.00</td></tr><tr><td>58712.27</td><td>K</td><td>21.00</td></tr><tr><td>58719.32</td><td>K</td><td>21.00</td></tr><tr><td>58751.21</td><td>K</td><td>21.00</td></tr><tr><td>58751.31</td><td>K</td><td>21.00</td></tr><tr><td>58712.34</td><td>K</td><td>21.00</td></tr><tr><td>58712.39</td><td>K</td><td>21.00</td></tr><tr><td>58719.40</td><td>K</td><td>21.00</td></tr><tr><td>58751.33</td><td>K</td><td>21.00</td></tr><tr><td>58712.33</td><td>K</td><td>21.00</td></tr><tr><td>58719.38</td><td>K</td><td>21.00</td></tr><tr><td>58750.37</td><td>K</td><td>21.00</td></tr><tr><td>58712.29</td><td>K</td><td>21.00</td></tr><tr><td>58719.35</td><td>K</td><td>21.00</td></tr><tr><td>58751.32</td><td>K</td><td>21.00</td></tr><tr><td>58714.22</td><td>K</td><td>21.00</td></tr><tr><td>58720.40</td><td>K</td><td>21.00</td></tr><tr><td>58752.28</td><td>K</td><td>21.00</td></tr><tr><td>58711.20</td><td>K</td><td>21.00</td></tr><tr><td>58719.23</td><td>K</td><td>21.00</td></tr><tr><td>58750.32</td><td>K</td><td>21.00</td></tr></tbody></table>
<table><caption>Table B.2. continued.</caption><thead><tr><th>Facility</th><th>Tile ID</th><th>∫<sub>tile</sub> P<sub>2D</sub> dΩ<br/>[%]</th><th>Σ<sub>tile</sub> P<sub>gal</sub><br/>[%]</th><th>RA<br/>[deg]</th><th>Dec<br/>[deg]</th><th>MJD</th><th>Filter</th><th>Limiting mag<br/>[AB]</th></tr></thead><tbody><tr><td rowspan="3">VISTA</td><td rowspan="3">10</td><td rowspan="3">5.010</td><td rowspan="3">1.799</td><td rowspan="3">14.489817</td><td rowspan="3">-25.399150</td><td>58712.17</td><td>K</td><td>21.00</td></tr><tr><td>58720.23</td><td>K</td><td>21.00</td></tr><tr><td>58750.12</td><td>K</td><td>21.00</td></tr><tr><td rowspan="3">VISTA</td><td rowspan="3">11</td><td rowspan="3">1.144</td><td rowspan="3">0.817</td><td rowspan="3">14.972704</td><td rowspan="3">-24.306950</td><td>58712.22</td><td>K</td><td>21.00</td></tr><tr><td>58720.27</td><td>K</td><td>21.00</td></tr><tr><td>58751.26</td><td>K</td><td>21.00</td></tr><tr><td rowspan="3">VISTA</td><td rowspan="3">12</td><td rowspan="3">2.520</td><td rowspan="3">2.361</td><td rowspan="3">12.988658</td><td rowspan="3">-27.491250</td><td>58714.21</td><td>K</td><td>21.00</td></tr><tr><td>58720.28</td><td>K</td><td>21.00</td></tr><tr><td>58752.20</td><td>K</td><td>21.00</td></tr><tr><td rowspan="3">VISTA</td><td rowspan="3">13</td><td rowspan="3">10.706</td><td rowspan="3">17.574</td><td rowspan="3">13.227671</td><td rowspan="3">-26.491250</td><td>58711.19</td><td>K</td><td>21.00</td></tr><tr><td>58719.22</td><td>K</td><td>21.00</td></tr><tr><td>58750.31</td><td>K</td><td>21.00</td></tr><tr><td rowspan="3">VISTA</td><td rowspan="3">14</td><td rowspan="3">5.098</td><td rowspan="3">3.471</td><td rowspan="3">14.252400</td><td rowspan="3">-24.441050</td><td>58712.40</td><td>K</td><td>21.00</td></tr><tr><td>58720.38</td><td>K</td><td>21.00</td></tr><tr><td>58751.37</td><td>K</td><td>21.00</td></tr><tr><td rowspan="3">VISTA</td><td rowspan="3">15</td><td rowspan="3">20.441</td><td rowspan="3">19.345</td><td rowspan="3">13.004717</td><td rowspan="3">-25.399150</td><td>58711.40</td><td>K</td><td>21.00</td></tr><tr><td>58720.22</td><td>K</td><td>21.00</td></tr><tr><td>58750.11</td><td>K</td><td>21.00</td></tr><tr><td rowspan="3">VISTA</td><td rowspan="3">16</td><td rowspan="3">2.230</td><td rowspan="3">2.282</td><td rowspan="3">14.070646</td><td rowspan="3">-23.348850</td><td>58714.18</td><td>K</td><td>21.00</td></tr><tr><td>58720.37</td><td>K</td><td>21.00</td></tr><tr><td>58752.18</td><td>K</td><td>21.00</td></tr><tr><td rowspan="3">VISTA</td><td rowspan="3">17</td><td rowspan="3">8.540</td><td rowspan="3">17.738</td><td rowspan="3">12.256171</td><td rowspan="3">-26.491250</td><td>58711.18</td><td>K</td><td>21.00</td></tr><tr><td>58719.21</td><td>K</td><td>21.00</td></tr><tr><td>58750.29</td><td>K</td><td>21.00</td></tr><tr><td rowspan="3">VISTA</td><td rowspan="3">18</td><td rowspan="3">12.904</td><td rowspan="3">7.413</td><td rowspan="3">13.439304</td><td rowspan="3">-24.306950</td><td>58712.21</td><td>K</td><td>21.00</td></tr><tr><td>58720.26</td><td>K</td><td>21.00</td></tr><tr><td>58751.25</td><td>K</td><td>21.00</td></tr><tr><td rowspan="3">VISTA</td><td rowspan="3">19</td><td rowspan="3">16.476</td><td rowspan="3">7.644</td><td rowspan="3">13.266400</td><td rowspan="3">-24.441050</td><td>58712.37</td><td>K</td><td>21.00</td></tr><tr><td>58720.35</td><td>K</td><td>21.00</td></tr><tr><td>58751.36</td><td>K</td><td>21.00</td></tr><tr><td rowspan="3">VISTA</td><td rowspan="3">20</td><td rowspan="3">0.990</td><td rowspan="3">0.634</td><td rowspan="3">11.821058</td><td rowspan="3">-27.491250</td><td>58714.20</td><td>K</td><td>21.00</td></tr><tr><td>58719.18</td><td>K</td><td>21.00</td></tr><tr><td>58752.05</td><td>K</td><td>21.00</td></tr><tr><td rowspan="3">VISTA</td><td rowspan="3">21</td><td rowspan="3">21.098</td><td rowspan="3">17.478</td><td rowspan="3">12.748617</td><td rowspan="3">-25.399150</td><td>58752.10</td><td>K</td><td>21.00</td></tr><tr><td>58752.12</td><td>K</td><td>21.00</td></tr><tr><td>58711.24</td><td>K</td><td>21.00</td></tr><tr><td rowspan="18">VISTA</td><td rowspan="18">22</td><td rowspan="18">5.813</td><td rowspan="18">8.108</td><td rowspan="18">13.273946</td><td rowspan="18">-23.348850</td><td>58711.26</td><td>K</td><td>21.00</td></tr><tr><td>58711.27</td><td>K</td><td>21.00</td></tr><tr><td>58711.28</td><td>K</td><td>21.00</td></tr><tr><td>58711.29</td><td>K</td><td>21.00</td></tr><tr><td>58711.30</td><td>K</td><td>21.00</td></tr><tr><td>58711.35</td><td>K</td><td>21.00</td></tr><tr><td>58711.36</td><td>K</td><td>21.00</td></tr><tr><td>58711.38</td><td>K</td><td>21.00</td></tr><tr><td>58711.39</td><td>K</td><td>21.00</td></tr><tr><td>58719.26</td><td>K</td><td>21.00</td></tr><tr><td>58719.27</td><td>K</td><td>21.00</td></tr><tr><td>58719.29</td><td>K</td><td>21.00</td></tr><tr><td>58719.42</td><td>K</td><td>21.00</td></tr><tr><td>58720.15</td><td>K</td><td>21.00</td></tr><tr><td>58720.16</td><td>K</td><td>21.00</td></tr><tr><td>58720.17</td><td>K</td><td>21.00</td></tr><tr><td>58720.18</td><td>K</td><td>21.00</td></tr><tr><td>58720.19</td><td>K</td><td>21.00</td></tr><tr><td>58720.21</td><td>K</td><td>21.00</td></tr><tr><td>58750.35</td><td>K</td><td>21.00</td></tr><tr><td>58751.07</td><td>K</td><td>21.00</td></tr><tr><td>58751.08</td><td>K</td><td>21.00</td></tr><tr><td>58751.11</td><td>K</td><td>21.00</td></tr><tr><td>58751.12</td><td>K</td><td>21.00</td></tr><tr><td>58751.13</td><td>K</td><td>21.00</td></tr><tr><td>58751.14</td><td>K</td><td>21.00</td></tr><tr><td>58751.16</td><td>K</td><td>21.00</td></tr><tr><td>58751.17</td><td>K</td><td>21.00</td></tr><tr><td>58751.18</td><td>K</td><td>21.00</td></tr><tr><td>58751.19</td><td>K</td><td>21.00</td></tr><tr><td rowspan="3">VISTA</td><td rowspan="3">23</td><td rowspan="3">18.996</td><td rowspan="3">16.287</td><td rowspan="3">12.246717</td><td rowspan="3">-25.399150</td><td>58714.17</td><td>K</td><td>21.00</td></tr><tr><td>58720.36</td><td>K</td><td>21.00</td></tr><tr><td>58751.38</td><td>K</td><td>21.00</td></tr><tr><td rowspan="3">VISTA</td><td rowspan="3">23</td><td rowspan="3">18.996</td><td rowspan="3">16.287</td><td rowspan="3">12.246717</td><td rowspan="3">-25.399150</td><td>58711.23</td><td>K</td><td>21.00</td></tr><tr><td>58719.25</td><td>K</td><td>21.00</td></tr><tr><td>58751.10</td><td>K</td><td>21.00</td></tr></tbody></table>
<table><caption>Table B.2. continued.</caption><thead><tr><th>Facility</th><th>Tile ID</th><th>∫<sub>tile</sub> P<sub>2D</sub> dΩ [%]</th><th>Σ<sub>tile</sub> P<sub>gal</sub> [%]</th><th>RA [deg]</th><th>Dec [deg]</th><th>MJD</th><th>Filter</th><th>Limiting mag [AB]</th></tr></thead><tbody><tr><td rowspan="12">VISTA</td><td rowspan="4">24</td><td rowspan="4">2.131</td><td rowspan="4">1.639</td><td rowspan="4">11.191771</td><td rowspan="4">-26.491250</td><td>58711.17</td><td>K</td><td>21.00</td></tr><tr><td>58719.14</td><td>K</td><td>21.00</td></tr><tr><td>58750.27</td><td>K</td><td>21.00</td></tr><tr><td>58750.28</td><td>K</td><td>21.00</td></tr><tr><td rowspan="4">25</td><td rowspan="4">16.608</td><td rowspan="4">9.775</td><td rowspan="4">11.966904</td><td rowspan="4">-24.306950</td><td>58712.20</td><td>K</td><td>21.00</td></tr><tr><td>58720.24</td><td>K</td><td>21.00</td></tr><tr><td>58750.15</td><td>K</td><td>21.00</td></tr><tr><td>58711.22</td><td>K</td><td>21.00</td></tr><tr><td rowspan="4">26</td><td rowspan="4">4.070</td><td rowspan="4">3.787</td><td rowspan="4">10.946017</td><td rowspan="4">-25.399150</td><td>58719.16</td><td>K</td><td>21.00</td></tr><tr><td>58750.33</td><td>K</td><td>21.00</td></tr><tr><td>58712.25</td><td>K</td><td>21.00</td></tr><tr><td>58720.33</td><td>K</td><td>21.00</td></tr><tr><td rowspan="4">27</td><td rowspan="4">6.280</td><td rowspan="4">9.353</td><td rowspan="4">11.613746</td><td rowspan="4">-23.214750</td><td>58751.28</td><td>K</td><td>21.00</td></tr><tr><td>58712.18</td><td>K</td><td>21.00</td></tr><tr><td>58719.17</td><td>K</td><td>21.00</td></tr><tr><td>58750.14</td><td>K</td><td>21.00</td></tr><tr><td rowspan="4">28</td><td rowspan="4">6.106</td><td rowspan="4">5.785</td><td rowspan="4">10.890504</td><td rowspan="4">-24.306950</td><td>58750.16</td><td>K</td><td>21.00</td></tr><tr><td>58712.24</td><td>K</td><td>21.00</td></tr><tr><td>58719.20</td><td>K</td><td>21.00</td></tr><tr><td>58751.27</td><td>K</td><td>21.00</td></tr><tr><td rowspan="2">29</td><td rowspan="2">2.701</td><td rowspan="2">2.745</td><td rowspan="2">10.778046</td><td rowspan="2">-23.214750</td><td>58710.58</td><td>z</td><td>19.65</td></tr><tr><td>58710.59</td><td>i</td><td>20.30</td></tr><tr><td rowspan="2">0</td><td rowspan="2">0.020</td><td rowspan="2">0.000</td><td rowspan="2">15.543640</td><td rowspan="2">-26.588360</td><td>58710.57</td><td>z</td><td>19.21</td></tr><tr><td>58710.57</td><td>i</td><td>20.26</td></tr><tr><td rowspan="2">1</td><td rowspan="2">0.067</td><td rowspan="2">0.000</td><td rowspan="2">15.195180</td><td rowspan="2">-26.989770</td><td>58710.59</td><td>i</td><td>20.26</td></tr><tr><td>58710.59</td><td>z</td><td>19.90</td></tr><tr><td rowspan="2">2</td><td rowspan="2">0.032</td><td rowspan="2">0.000</td><td rowspan="2">15.550540</td><td rowspan="2">-26.186890</td><td>58710.59</td><td>z</td><td>19.41</td></tr><tr><td>58710.59</td><td>i</td><td>20.07</td></tr><tr><td rowspan="2">3</td><td rowspan="2">0.001</td><td rowspan="2">0.000</td><td rowspan="2">16.007700</td><td rowspan="2">-25.389260</td><td>58710.57</td><td>z</td><td>20.23</td></tr><tr><td>58710.57</td><td>i</td><td>19.77</td></tr><tr><td rowspan="2">4</td><td rowspan="2">0.050</td><td rowspan="2">0.000</td><td rowspan="2">14.751030</td><td rowspan="2">-27.396890</td><td>58710.58</td><td>z</td><td>20.04</td></tr><tr><td>58710.59</td><td>i</td><td>20.48</td></tr><tr><td rowspan="2">5</td><td rowspan="2">0.058</td><td rowspan="2">0.000</td><td rowspan="2">15.188230</td><td rowspan="2">-26.588360</td><td>58710.59</td><td>i</td><td>20.27</td></tr><tr><td>58710.59</td><td>z</td><td>19.95</td></tr><tr><td rowspan="2">6</td><td rowspan="2">0.035</td><td rowspan="2">0.000</td><td rowspan="2">15.557400</td><td rowspan="2">-25.785410</td><td>58710.58</td><td>i</td><td>20.44</td></tr><tr><td>58710.59</td><td>z</td><td>20.08</td></tr><tr><td rowspan="2">7</td><td rowspan="2">0.174</td><td rowspan="2">0.085</td><td rowspan="2">14.745580</td><td rowspan="2">-26.995480</td><td>58710.58</td><td>i</td><td>20.44</td></tr><tr><td>58710.59</td><td>z</td><td>20.08</td></tr><tr><td rowspan="2">8</td><td rowspan="2">0.155</td><td rowspan="2">0.074</td><td rowspan="2">14.303600</td><td rowspan="2">-27.802590</td><td>58710.59</td><td>z</td><td>19.87</td></tr><tr><td>58710.59</td><td>z</td><td>20.06</td></tr><tr><td rowspan="2">9</td><td rowspan="2">0.058</td><td rowspan="2">0.064</td><td rowspan="2">15.181320</td><td rowspan="2">-26.186890</td><td>58710.59</td><td>i</td><td>20.38</td></tr><tr><td>58710.59</td><td>z</td><td>20.20</td></tr><tr><td rowspan="2">10</td><td rowspan="2">0.013</td><td rowspan="2">0.000</td><td rowspan="2">15.564210</td><td rowspan="2">-25.383940</td><td>58710.58</td><td>i</td><td>20.20</td></tr><tr><td>58710.58</td><td>z</td><td>20.03</td></tr><tr><td rowspan="2">11</td><td rowspan="2">0.218</td><td rowspan="2">0.301</td><td rowspan="2">14.740170</td><td rowspan="2">-26.593980</td><td>58710.58</td><td>i</td><td>20.51</td></tr><tr><td>58710.59</td><td>z</td><td>20.23</td></tr><tr><td rowspan="2">12</td><td rowspan="2">0.057</td><td rowspan="2">0.000</td><td rowspan="2">13.898950</td><td rowspan="2">-28.191420</td><td>58710.58</td><td>z</td><td>19.93</td></tr><tr><td>58710.58</td><td>z</td><td>19.90</td></tr><tr><td rowspan="2">13</td><td rowspan="2">0.165</td><td rowspan="2">0.014</td><td rowspan="2">14.299680</td><td rowspan="2">-27.401250</td><td>58710.58</td><td>z</td><td>19.74</td></tr><tr><td>58710.58</td><td>i</td><td>20.55</td></tr><tr><td rowspan="2">14</td><td rowspan="2">0.106</td><td rowspan="2">0.003</td><td rowspan="2">15.174460</td><td rowspan="2">-25.785400</td><td>58710.59</td><td>z</td><td>20.26</td></tr><tr><td>58710.59</td><td>i</td><td>20.26</td></tr><tr><td rowspan="2">15</td><td rowspan="2">0.013</td><td rowspan="2">0.000</td><td rowspan="2">15.570980</td><td rowspan="2">-24.982530</td><td>58710.58</td><td>i</td><td>20.26</td></tr><tr><td>58710.59</td><td>z</td><td>20.10</td></tr><tr><td rowspan="2">16</td><td rowspan="2">0.252</td><td rowspan="2">0.045</td><td rowspan="2">14.295780</td><td rowspan="2">-26.999770</td><td>58710.58</td><td>i</td><td>20.34</td></tr><tr><td>58710.59</td><td>z</td><td>20.20</td></tr><tr><td rowspan="2">17</td><td rowspan="2">0.158</td><td rowspan="2">0.101</td><td rowspan="2">14.734800</td><td rowspan="2">-26.192410</td><td>58710.58</td><td>i</td><td>20.63</td></tr><tr><td>58710.58</td><td>z</td><td>20.46</td></tr><tr><td rowspan="2">18</td><td rowspan="2">0.168</td><td rowspan="2">0.018</td><td rowspan="2">13.850540</td><td rowspan="2">-27.805550</td><td>58710.58</td><td>z</td><td>20.03</td></tr><tr><td>58710.58</td><td>i</td><td>20.54</td></tr><tr><td rowspan="2">19</td><td rowspan="2">0.078</td><td rowspan="2">0.247</td><td rowspan="2">15.167650</td><td rowspan="2">-25.383940</td><td>58710.59</td><td>z</td><td>20.43</td></tr><tr><td>58710.59</td><td>z</td><td>19.83</td></tr><tr><td rowspan="2">20</td><td rowspan="2">0.004</td><td rowspan="2">0.000</td><td rowspan="2">13.441750</td><td rowspan="2">-28.593110</td><td>58710.57</td><td>z</td><td>19.74</td></tr><tr><td>58710.57</td><td>z</td><td>19.79</td></tr><tr><td rowspan="2">21</td><td rowspan="2">0.016</td><td rowspan="2">0.000</td><td rowspan="2">15.577690</td><td rowspan="2">-24.581220</td><td>58710.58</td><td>i</td><td>20.67</td></tr><tr><td>58710.59</td><td>z</td><td>20.48</td></tr><tr><td rowspan="2">22</td><td rowspan="2">0.370</td><td rowspan="2">0.142</td><td rowspan="2">14.291920</td><td rowspan="2">-26.598190</td><td>58710.58</td><td>i</td><td>20.48</td></tr><tr><td>58710.59</td><td>z</td><td>20.73</td></tr><tr><td rowspan="2">23</td><td rowspan="2">0.189</td><td rowspan="2">0.015</td><td rowspan="2">14.729460</td><td rowspan="2">-25.790820</td><td>58710.58</td><td>i</td><td>20.48</td></tr><tr><td>58710.59</td><td>z</td><td>20.48</td></tr><tr><td rowspan="2">24</td><td rowspan="2">0.343</td><td rowspan="2">0.318</td><td rowspan="2">13.848190</td><td rowspan="2">-27.404160</td><td>58710.58</td><td>z</td><td>20.08</td></tr><tr><td>58710.58</td><td>z</td><td>20.14</td></tr><tr><td rowspan="2">25</td><td rowspan="2">0.050</td><td rowspan="2">0.000</td><td rowspan="2">13.442660</td><td rowspan="2">-28.191420</td><td>58710.58</td><td>i</td><td>20.64</td></tr><tr><td>58710.58</td><td>z</td><td>20.54</td></tr><tr><td rowspan="2">26</td><td rowspan="2">0.052</td><td rowspan="2">0.000</td><td rowspan="2">15.160890</td><td rowspan="2">-24.982530</td><td>58710.58</td><td>z</td><td>20.54</td></tr><tr><td>58710.58</td><td>z</td><td>19.75</td></tr><tr><td rowspan="2">27</td><td rowspan="2">0.012</td><td rowspan="2">0.000</td><td rowspan="2">15.584370</td><td rowspan="2">-24.180040</td><td>58710.58</td><td>z</td><td>20.19</td></tr><tr><td>58710.59</td><td>i</td><td>20.76</td></tr><tr><td rowspan="2">28</td><td rowspan="2">0.373</td><td rowspan="2">0.363</td><td rowspan="2">14.288080</td><td rowspan="2">-26.196550</td><td>58710.58</td><td>i</td><td>20.50</td></tr><tr><td>58710.58</td><td>z</td><td>20.50</td></tr></tbody></table>