Search for Gravitational Waves Associated with Gamma-Ray Bursts Detected by Fermi and Swift during the LIGO–Virgo Run O3a

We search for gravitational-wave transients associated with gamma-ray bursts (GRBs) detected by the Fermi and Swift satellites during the first part of the third observing run of Advanced LIGO and Advanced Virgo (2019 April 1 15:00 UTC–2019 October 1 15:00 UTC). A total of 105 GRBs were analyzed using a search for generic gravitational-wave transients; 32 GRBs were analyzed with a search that specifically targets neutron star binary mergers as short GRB progenitors. We find no significant evidence for gravitational-wave signals associated with the GRBs that we followed up, nor for a population of unidentified subthreshold signals. We consider several source types and signal morphologies, and report for these lower bounds on the distance to each GRB.

Gamma-ray bursts (GRBs) are transient flashes of gamma radiation of cosmological origin observed at a rate of 1 per day (Nakar 2007).The interaction of matter with a compact central object, e.g., an accreting black hole (BH; Woosley 1993;Popham et al. 1999) or a magnetar (Usov 1992;Zhang & Meszaros 2001), is believed to drive highly relativistic jets which power the prompt emission of these astrophysical events.GRBs are broadly grouped into two classes -long and short GRBs -depending on the duration and spectral hardness of their prompt emission (Mazets et al. 1981;Norris et al. 1984;Kouveliotou et al. 1993).
Long, soft GRBs have durations 2 s and are firmly associated by optical observations to the collapse of massive stars (Galama et al. 1998;Hjorth et al. 2003;Stanek et al. 2003;Hjorth & Bloom 2012).Gravitational waves (GWs) will be radiated by the core-collapse process, (e.g., Fryer & New 2011).Several models of this process do not yield radiation that is detectable by the current generation of GW interferometers beyond Galactic distances (Abbott et al. 2020c).However, rotational instabilities and instabilities induced by the additional presence of an accretion disk as part of the GRB engine may enhance the GW emission, making it detectable even for extragalactic sources (van Putten 2001; Davies et al. 2002;Fryer et al. 2002;Kobayashi & Meszaros 2003;Shibata et al. 2003;Piro & Pfahl 2007;Corsi & Meszaros 2009;Romero et al. 2010;Gossan et al. 2016;Abbott et al. 2020c).
The unambiguous association (Abbott et al. 2017a) of neutron star (NS) binary merger GW170817 (Abbott et al. 2017b(Abbott et al. , 2019d) and short GRB 170817A (Goldstein et al. 2017;Savchenko et al. 2017) has confirmed that compact binary mergers of this kind can produce short GRBs.This milestone in multimessenger astronomy corroborated the idea first proposed in the 1980s (Blinnikov et al. 1984;Paczynski 1986;Eichler et al. 1989;Paczynski 1991;Narayan et al. 1992) that the progenitors of short GRBs are compact binaries containing NSs (for a review of proposed progenitors, see Lee & Ramirez-Ruiz 2007;Nakar 2007).Indirect evidence that had previously reinforced this idea was due to the observation of a possible kilonova associated with GRB 130603B (Berger et al. 2013;Tanvir et al. 2013), and to numerous studies of the environments of short GRBs (for reviews see Berger 2011Berger , 2014)), starting with the afterglow observation and host-galaxy association of GRB 050509B (Castro-Tirado et al. 2005;Gehrels et al. 2005;Bloom et al. 2006).
In addition to confirming the origin of some short GRBs, combining data from observations of GW170817 and GRB 170817A allowed for the inference of basic properties of short GRB jets.These include the isotropic equivalent luminosity of the jet, determined through a redshift measurement made possible by the optical follow-up of the GW localization (Abbott et al. 2017a;Goldstein et al. 2017), and the geometry of the GRB jets (Williams et al. 2018;Farah et al. 2020;Mogushi et al. 2019).The precise mechanism by which the jet is launched is still unknown, although it is typically believed to be either neutrino-driven or magnetically driven (Nakar 2007, but see also Liu et al. 2015, and references therein).Indeed, the scientific debate about the emission profile of the jet and the subsequent gammaray production mechanism of GRB 170817A is still ongoing (Hallinan et al. 2017;Kasliwal et al. 2017;Lamb & Kobayashi 2017;Troja et al. 2017;Gottlieb et al. 2018b;Lazzati et al. 2018;Gill & Granot 2018;Mooley et al. 2018;Zhang et al. 2018;Ghirlanda et al. 2019).It is generally believed that there are symmetric polar outflows of highly relativistic material that travel parallel to the total angular momentum of the binary system (Aloy et al. 2004;Kumar & Zhang 2014;Murguia-Berthier et al. 2017).These jets are thought to be collimated and roughly axisymmetric, emitting preferentially in a narrow opening angle due to a combination of outflow geometry and relativistic beaming.The data from extensive multi-wavelength observation campaigns that ran for nearly 20 months following the merger (Fong et al. 2019;Makhathini et al. 2020;Troja et al. 2020) are in agreement with a structured jet model, in which the energy and bulk Lorentz factor gradually decrease with angular distance from the jet symmetry axis (e.g., Dai & Gou 2001;Lipunov et al. 2001;Rossi et al. 2002;Zhang & Mészáros 2002;Ghirlanda et al. 2019;Beniamini et al. 2020).Further, according to one of the models proposed, as the jet drills through the surrounding merger ejecta it inflates a mildly relativistic cocoon due to interactions between the material at the edge of the jet and the ejecta (Lazzati et al. 2017;Gottlieb et al. 2018a).In this case, it is possible that the cocoon alone could produce the gamma-rays observed from GRB 170817A (Gottlieb et al. 2018b).Additional joint detections of GRBs and GWs will significantly aid in our understanding of the underlying energetics (Lamb & Kobayashi 2017;Wu & MacFadyen 2018;Burns et al. 2019), jet geometry (Farah et al. 2020;Mogushi et al. 2019;Biscoveanu et al. 2020;Hayes et al. 2020), and jet ignition mechanisms (Veres et al. 2018;Ciolfi et al. 2019;Zhang 2019) of binary neutron star (BNS) coalescences.
A targeted search for GWs in sky and time coincidence with GRBs enhances our potential of achieving such joint detections.In this paper we present our results for the targeted GW follow-up of GRBs reported during the first part of the third observing run of Advanced LIGO and Advanced Virgo (O3a) by the Fermi (Meegan et al. 2009) and Swift (Gehrels et al. 2004;Barthelmy et al. 2005) satellites.As in the first (Abbott et al. 2017c) and second (Abbott et al. 2017a(Abbott et al. , 2019b) ) observing runs, two searches with different assumptions about signal morphology are applied to the GW data: we process all GRBs with a search for generic GW transients (X-Pipeline; Sutton et al. 2010;Was et al. 2012, see Sec. 3.2 for details) and we follow up short GRBs with an additional, modeled search for BNS and neutron star-black hole (NSBH) GW inspiral signals (PyGRB; Harry & Fairhurst 2011;Williamson et al. 2014, see Sec. 3.1 for details).These searches were able to process 105 and 32 GRBs in O3a, respectively.
The scope of these targeted searches is to enhance our ability to detect GW signals in coincidence with GRBs with respect to all-sky searches for transient GW signals carried out by the LIGO Scientific & Virgo Collaboration (Abbott et al. 2019c;Abbott et al. 2021).These may lead to joint GW-GRB detections in the case of loud GW events, as for GW170817 and GRB 170817A, but the targeted searches we report on here aim at uncovering sub-threshold GW signals by exploiting the time and localization information of the GRBs themselves.The Fermi Gamma-ray Burst Monitor (GBM) team conducts an analogous effort when searching through GBM data for gamma-ray transients coincident with confirmed events and low-significance candidates reported by LIGO-Virgo offline analyses (Hamburg et al. 2020).Similarly, the Swift/Burst Alert Telescope (BAT) team has developed their own autonomous pipeline to enable subthreshold GRB searches for externally triggered events (Tohuvavohu et al. 2020).
This first part of the third observing run took place between 2019 April 1 15:00 UTC and 2019 October 1 15:00 UTC.Setting the false-alarm-rate threshold to two per year, 39 compact binary coalescence events were identified in O3a (Abbott et al. 2021).The majority of these have been classified as signals emitted by binary BH mergers; however, three events have the possibility of coming from a binary with at least one NS, that is, a potential short GRB progenitor.
1. GW190425 (Abbott et al. 2020a) was a compact binary coalescence with primary mass 2.0 +0.6 −0.3 and secondary mass 1.4 +0.3  −0.3 (all measurements quoted at the 90% credible level) and is therefore consistent with being the result of a BNS merger (Abbott et al. 2020a;Abbott et al. 2021).
2. GW190426 was the GW candidate event with the highest false-alarm rate in Abbott et al. (2021); assuming it is a real signal, its inferred component masses of 5.7 +3.9 −2.3 and 1.5 +0.8 −0.5 indicate that it may have originated from an NSBH, or a binary BH merger.
3. GW190814 (Abbott et al. 2020b) could have originated from an NSBH, or a binary BH merger, as it has a primary mass measurement of 23.2 +1.1 −1.0 and posterior support for a secondary mass 2.59 +0.08 −0.09 .This makes the secondary compact object either the lightest BH or the heaviest NS known to be in a compact binary system.
While there is considerable uncertainty in source type for all three of these events, GW190425 is the one for which the prospects of observing an associated GRB were most promising, as it is consistent with a BNS merger, rather than a binary BH merger or an NSBH merger with high or moderately high mass ratio.However, no confirmed electromagnetic or neutrino counterparts were observed in association with this event (Hosseinzadeh et al. 2019, Lundquist et al. 2019, Abbott et al. 2020a, Coughlin et al. 2020; see also Pozanenko et al. 2020) despite extensive searches, which are logged in the Gamma-ray Coordinates Network (GCN) Circular archive. 204There are a number of reasons for which an electromagnetic counterpart associated with GW190425 may not have been detected.First, the large area covered by the localization region of GW190425 determined from GW data (> 8000 deg 2 ) posed a considerable challenge for electromagnetic follow-up.45.4% of this localization region was occulted by the Earth for the Fermi satellite, so, if gamma-rays were emitted from the source, it is possible they were not detectable.Other gamma-ray observatories with lower sensitivities to short GRBs, such as INTEGRAL and KONUS-Wind, were covering relevant fractions of the localization region, however (Martin-Carrillo et al. 2019;Svinkin et al. 2019).Second, GRB jets are expected to be aligned with the total angular momentum of the binary system, and thus more easily detectable at small viewing angles.The binary inclination angle of GW190425 was poorly constrained, so it is possible that a jet from this system was formed but was oriented away from our line of sight.Additionally, the luminosity distance inferred for GW190425 (∼ 160 Mpc) was significantly larger than that for GW170817 (∼ 40 Mpc).GRB 170817A, which followed GW170817, was such an exceptionally faint short GRB (Abbott et al. 2017a) that its prompt emission photon flux would have dipped below the detection threshold for Fermi-GBM, had the source been farther than ∼ 75 Mpc (Abbott et al. 2017a;Goldstein et al. 2017), and by ∼ 100 Mpc it would become undetectable by Swift/BAT (Tohuvavohu et al. 2020).Thus, if emission from the system that produced GW190425 was similarly faint, it would not have been detectable by Swift/BAT or Fermi-GBM.Therefore, we do not necessarily expect a GRB detection to be associated with GW190425 due to its almost unconstrained inclination angle, large localization region, and distance, even if gamma-rays were emitted from this system.Scenarios like this one further motivate the need for GW follow-up analyses of GRB events which, by definition, constrain the sky localization and inclination angle of the progenitor.
In Section 2 we discuss the set of GRBs analyzed in this paper.In Section 3, we summarize the two targeted search methods used to follow up GRBs.Section 4 presents the results obtained with these two methods.We also consider each of the two sets of results collectively and quantify its consistency with the no-signal hypothesis.Finally, in Section 5 we provide our concluding remarks.

GRB SAMPLE
The sample of GRBs analyzed in this paper includes events circulated by the GCN, 205 complemented with information from the Swift/BAT catalog (Lien et al. 2016), 206 the online Swift GRBs Archive, 207 and the Fermi GBM Catalog. 208(Gruber et al. 2014;von Kienlin et al. 2014;Narayana Bhat et al. 2016) Once an alert detailing an event has been received via the GCN, the dedicated Vetting Automation and Literature Informed Database (Coyne 2015) is applied to find the latest GRB results by comparing the time and localization parameters with those in tables relating to each satellite, the published catalogs, and an automatic literature search.The GCN notices provide a set of 141 GRBs during the O3a data-taking period (2019 April 1 15:00 UTC -2019 October 1 15:00 UTC).
As mentioned in the Introduction, we carry out two searches with distinct assumptions about signal morphology (see Sec. 3 for details on both methods): a search for generic GW transients and a modeled search for GW signals from NS binary inspirals, i.e., BNS and NSBH.We do this because GRBs of different durations are expected to have different origins and therefore different GW signal morphologies.In particular, if a compact binary merger were to produce a GRB it would be expected to have a short duration.In order to specifically target such phenomena with the modeled search, we classify each GRB as long, short, or ambiguous.This classification relies on the measurement of the time interval over which 90% of the total backgroundsubtracted photon counts are observed (T 90 , with error |δT 90 |).When T 90 + |δT 90 | < 2 s the GRBs are labeled as short, when T 90 − |δT 90 | > 4 s the GRBs are labeled as long, the rest are labeled as ambiguous.The unmodeled search for generic transients is applied to GRBs of all classifications.All of the short and ambiguous GRBs are additionally analyzed with the modeled search in order to maximize the chances of uncovering any potential binary coalescence candidate.
The classification process results in 20 short GRBs, 108 long GRBs, and 13 ambiguous GRBs.As in Abbott et al. (2019b), we require a minimum amount of coincident data from at least two GW detectors around the time of a GRB for the generic unmodeled GW transient search to assess the significance of a GW candidate with sub-percent level accuracy (see Sec. 3.2 for technical details).This requirement is applied to GRBs of all classifications and results in 105 GRBs being analyzed with this method, out of the 141 GRBs recorded by Fermi and Swift during O3a.This amounts to 74.5%, a percentage of events that is compatible with the fraction of observing time during which at least two interferometers in the network were operating in observing mode (81.9 %; Abbott et al. 2021).Similarly, requirements from the modeled search (see Sec. 3.1 for technical details) set the minimum amount of data needed from at least one detector around the time of the GRBs.It leads to 32 short and ambiguous GRBs being analyzed with this method, 209 that is, 97.0% of the 33 possible ones.This value matches the fraction of observing time during which at least one interferometer in the network was operating in observing mode during O3a (96.9 %; Abbott et al. 2021).
Of the 141 Fermi and Swift GRBs in our sample, the vast majority do not have redshift measurements.Those that do are the ambiguous GRB 190627A at z = 1.942 (Japelj et al. 2019), and the two long GRBs 190719C and 190829A at z = 2.469 and z = 0.0785, respectively (Rossi et al. 2019;Valeev et al. 2019).All three fall beyond the detection range of our interferometers, and are not expected to produce measurable GW results.Regardless of availability of redshift information, however, we followed up as many GRBs as we could and we were indeed able to analyze these three cases.

SEARCH METHODS
We now provide a description of the two targeted search methods used in this paper.These are the same methods applied to GW data coincident with GRBs that occurred during the first (Abbott et al. 2017c) and second (Abbott et al. 2017a(Abbott et al. , 2019b) ) Advanced LIGO and Virgo observing runs.In Sec.3.1 we summarize the modeled search method that aims at uncovering subthreshold GW signals emitted by BNS and NSBH binaries (PyGRB; Harry & Fairhurst 2011;Williamson et al. 2014).In Sec.3.2 we discuss the search for generic GW transients (X-Pipeline; Sutton et al. 2010;Was et al. 2012).Results from these two searches are presented in Sec. 4.

Modeled search for binary mergers
This analysis searches for a GW signal compatible with the inspiral of a BNS or NSBH binary -collectively NS binaries -within 6 s of data associated with an observed short GRB.This stretch of data is the onsource window and runs from −5 s to +1 s around the start of the GRB emission (i.e., the GRB trigger time).The surrounding ∼30-90 minutes of data are split into 6 s off-source trials which are also analyzed in order to build a background.Around 30 minutes allows the modeled search to accurately estimate the power spectral density of the available instruments and ensures that it can assess at sub-percent level accuracy the significance of any candidate events found in the on-source window.All the data are processed using PyGRB (Harry GRB 190626254, GRB 190706B, GRB 190714573, GRB 190716917, GRB 190719113, GRB 190723309, GRB 190731943, GRB 190804792, GRB 190806675, GRB 190808498, GRB 190814837, GRB 190821A, GRB 190821716, GRB 190828614 .& Fairhurst 2011; Williamson et al. 2014), a coherent matched filtering pipeline that is part of the general open-source software PyCBC (Nitz et al. 2020) and has core elements in the LALSuite software library (LIGO Scientific Collaboration 2018).We scan each trial of data and the on-source window in the 30-1000 Hz frequency band using a predefined bank of waveform templates (Owen & Sathyaprakash 1999) created with a hybrid geometric-stochastic method (Capano et al. 2016;Dal Canton & Harry 2017) and using a phenomenological inspiral-merger-ringdown waveform model for nonprecessing point-particle binaries (IMRPhenomD; Husa et al. 2016;Khan et al. 2016). 210The waveform template bank includes waveforms corresponding to a range of masses ([1.0, 2.8]M for NSs, [1.0, 25.0]M for BHs) and dimensionless spin magnitudes ([0, 0.05] for NSs, [0, 0.998] for BHs) for aligned-spin, zero-eccentricity BNS or NSBH systems that may produce an electromagnetic counterpart via the tidal disruption of the NS (Pannarale & Ohme 2014).Aside from the updated sensitivity of our detectors, the only difference with respect to the second LIGO-Virgo observing run (Abbott et al. 2019b) is that the generation of the bank has been updated to apply more accurate physics to determine whether an NSBH system could produce an accretion disk from this disruption (Foucart et al. 2018).We only search for circularly polarized GWs, which may be emitted by binaries with inclinations of 0 • or 180 • : such systems have GW amplitudes that are consistent (Williamson et al. 2014) with those of binary progenitors with inclination angles over the full range of viewing angles that we expect for typical brightness GRBs ( 30 • ; Fong et al. 2015), such as those in our sample.
The strength of any potential signal is ranked via a coherent matched filter signal-to-noise ratio (S/N; Harry & Fairhurst 2011;Williamson et al. 2014) which is reweighted according to a χ 2 goodness-of-fit between the template that identified it and the signal itself.The significance of the latter is quantified as the probability of background alone producing such an event.This is evaluated by comparing the re-weighted S/N of the loudest trigger within the 6 s on-source to the distribution of the re-weighted S/Ns of the loudest triggers in the 6 s off-source trials.When data from more than one detector are available, this background S/N distribution is extended by generating additional off-source trials via time slides, that is, by combining data from detectors after introducing time shifts longer than the light-travel time across the network.Specifically, our time shifts are 6 s long, in order to match the width of the on-source window and the off-source trials.
In order to derive the sensitivity of this search to potential GRB sources, simulated signals are injected in software into the off-source data.The 90% exclusion distances, D 90 , are defined as the distances within which 90% of the injected simulated signals are recovered with a greater ranking statistic than the loudest on-source event.Three different astrophysical populations are considered: BNS binaries with generically oriented -i.e., precessing -spins, aligned spin NSBH binaries, and NSBH binaries with generically oriented spins.These simulated signals cover a portion of parameter space that extends beyond that covered by the template bank, as they include NS dimensionless spin values up to 0.4 and, for two families of injected signals, admit precession.As stated previously, the templates used to filter the data are produced using IMRPhenomD.In order to factor into the sensitivity assessment any potential loss due to uncertainties in GW signal modeling, the injected signals are not produced with the same model used for the templates.Precessing BNS signals are simulated using the TaylorT2 time-domain, post-Newtonian inspiral approximant (SpinTaylorT2; Sathyaprakash & Dhurandhar 1991;Blanchet et al. 1996;Mikoczi et al. 2005;Arun et al. 2009;Bohé et al. 2013Bohé et al. , 2015;;Mishra et al. 2016), while NSBH-injected waveforms are generated assuming a point-particle effective-one-body model tuned to numerical simulations which can allow for precession effects from misaligned spins (SEOBNRv3; Pan et al. 2014;Taracchini et al. 2014;Babak et al. 2017).The three populations used to build the injected signals are defined as in the first two LIGO-Virgo observing runs, to allow for direct comparisons (Abbott et al. 2017c(Abbott et al. , 2019b)).NS masses for the injections are taken between 1 M and 3 M from a normal distribution centered at 1.4 M with a standard deviation of 0.2 M (Kiziltan et al. 2013) and 0.4 M for BNS and NSBH systems, respectively.BH masses are taken to be between 3 M and 15 M from a normal distribution centered at 10 M with a standard deviation of 6 M .Spins are drawn uniformly in magnitude and, when applicable, with random orientation; the maximum allowed NS spin magnitude is 0.4, from the fastest observed pulsar spin (Hessels et al. 2006), while the maximum BH spin magnitude is set to 0.98, motivated by X-ray binary observations (e.g., Özel et al. 2010;Kreidberg et al. 2012;Miller & Miller 2014).Injected signals have a range of total inclinations from 0 • -30 • and 150 • -180 • whilst re-moving any systems which could not feasibly produce a short GRB (Pannarale & Ohme 2014).

Unmodelled search for generic transients
X-Pipeline looks for excess power that is coherent across the network of GW detectors and consistent with the sky localization and time window for each GRB.As in the first two observing runs, we use a search time window that begins 600 s before the GRB trigger time and ends 60 s after it, or at the T 90 time itself (whichever is larger).This window is long enough to encapsulate the time delay between GW emission from a progenitor and the GRB prompt emission (Koshut et al. 1995;Aloy et al. 2000;MacFadyen et al. 2001;Zhang et al. 2003;Lazzati 2005;Wang & Meszaros 2007;Burlon et al. 2008Burlon et al. , 2009;;Lazzati et al. 2009;Vedrenne & Atteia 2009).Our frequency range is restricted to the most sensitive band of the GW detectors, namely 20-500 Hz.While gravitational radiation from core-collapse supernovae is expected to contain frequency content above this band (Radice et al. 2019), detection of bursts above a few hundred hertz is not energetically favorable (see, e.g., Fig. 4 in Abbott et al. 2019a) and increasing the frequency upper limit also increases the computational cost.
The generic transient search pipeline coherently combines data from all detectors and produces timefrequency maps of this GW data stream.The maps are scanned for clusters of pixels with excess energy, referred to as events.The events obtained this way are first ranked according to a detection statistic based on energy and then subject to coherent consistency tests.These are based on correlations between data in different detectors and reject events associated with noise transients.The surviving event with the largest ranking statistic is taken to be the best candidate for a GW detection.Its significance is evaluated in the same way as in the modeled analysis, but with 660 s long off-source trials.In order to ensure that the significance is assessed at a sub-percent level, we require at least ∼ 1.5 hr of coincident data from at least two detectors around the time of a GRB.Non-Gaussian noise transients, or glitches, are handled as described in Abbott et al. (2019b).
Similarly to the modeled search, we quantify the sensitivity of the generic transient search by injecting simulated signals into off-source data in software and recovering them.Calibration errors are accounted for by jittering the amplitude and arrival time of the injections according to a Gaussian distribution representative of the calibration uncertainties in O3a (Abbott et al. 2017c).We report results obtained for four dis- tinct sets of circular sine-Gaussian (CSG) waveforms, with fixed quality factor Q = 9 and with central frequencies of 70, 100, 150, and 300 Hz (see Equation 1and Section 3.2 of Abbott et al. 2017c).These models are intended to represent the GWs from stellar collapses.In all four cases, we set the total radiated energy to E GW = 10 −2 M c 2 , a choice that is about an order of magnitude higher than the results presented in Abbott et al. (2020c) for the detectability of core-collapse supernovae.As optimistic representatives (Ott & Santamaría 2013) of longer-duration GW signals detectable by the unmodeled search, we use accretion disk instability (ADI) waveforms (van Putten 2001;van Putten et al. 2014).In these ADI models, instabilities form in a magnetically suspended torus around a rapidly spinning BH, causing GWs to be emitted.The model specifics and parameters used to generate the five families of ADI signals that we consider are the same as in Table 1 and Section 3.2 of Abbott et al. (2017c).

RESULTS
During O3a we used the generic transient method to follow up a total of 105 GRBs, whereas the modeled search was applied to the 32 GRB triggers classified as short or ambiguous.For all of the most GW-signallike triggers associated with the examined GRBs, the searches returned no significant probability of incompatibility with background alone (p-value).This indicates that no GW signal was uncovered in association with any of these GRBs.This is consistent with the estimated GW-GRB joint detection rate with Fermi-GBM of 0.07-1.80 per year reported in Abbott et al. (2019b) for the 2019-2020 LIGO-Virgo observing run.The most significant events found by the generic transient method and by the modeled search had p-values of 5.5 × 10 −3 (GRB 190804058) and 2.7 × 10 −2 (GRB 190601325), respectively.
Figures 1 and 2 show the cumulative distributions of p-values returned by the modeled search and the generic transient search, respectively.For cases in which no associated on-source trigger survived the analysis cuts of the modeled search, the associated p-value ranges between 1 -i.e. an upper bound on a probability -and the fraction of background trials for the GRB that also yielded no associated GW trigger.In both figures, the expected background distribution under the no-signal hypothesis is shown by the dashed line, and its 2σ limits are indicated by the two dotted lines.Both cumulative distributions are within the 2σ lines and therefore compatible with the no-signal hypothesis.These figures .Cumulative histograms of the 90% confidence exclusion distances, D90, for the binary neutron star (blue, thin line) and generically spinning neutron star-black hole (orange, thick line) signal models, shown for the sample of 32 short and ambiguous gamma-ray bursts that were followed up by the NS binary modeled search during O3a, none of which had an identified gravitational wave counterpart.For a given GRB event and signal model, D90 is the distance within which 90% of simulated signals inserted into off-source data are recovered with greater significance than the most significant on-source trigger.These simulated signals have inclinations θJN -the angle between the total angular momentum and the line of sight -drawn uniformly in sin θJN with θJN restricted to [0 indicate that the lowest p-value found by each search is compatible with the no-signal hypothesis.
Having found no GW signal associated with the GRBs followed up by our searches, we consider the set of modeled search results and the set of generic transient search results, collectively.We apply a weighted binomial test described in Abadie et al. (2012) to evaluate how consistent each set of results is collectively with the nosignal hypothesis.This test is conducted using the most significant 5% of p-values in the sample weighted by a prior probability of detection estimated using the network detector sensitivity at the time and location of each GRB.This final probability of observing this distribution of p-values given background alone, i.e. under the no-signal hypothesis, was 0.43 (0.31) for the modeled (generic transient) search method.Therefore, both searches gave no significant evidence for a population of unidentified subthreshold GW signals.For the analyses carried out in the first observing run of Advanced LIGO and Advanced Virgo, the combined p-values were 0.57 and 0.75 for the modeled and generic transient search, respectively (Abbott et al. 2017c); in the second observing run of Advanced LIGO and Advanced Virgo, they were 0.30 and 0.75 (Abbott et al. 2019b).
In Fig. 3, we show the cumulative 90% exclusion distances for the 32 short and ambiguous GRBs followed up with the modeled search.The lowest exclusion distance values (∼ 20 Mpc) were obtained for ambiguous GRB 190409901.This is due to the fact that only Virgo data were available for this GRB and that the sky location of this event was in a direction in which Virgo had ∼ 30% sensitivity with respect to an optimal sky location.For each of the three simulated signal classes, we quote the median of the 32 D 90 results in the top part of Table 1.All three values are 40 %-60 % times higher than those reported in Abbott et al. (2019b)  also place a 90% confidence level lower limit on the distance for each of the 105 GRBs analyzed by the generic transient search, assuming the various emission models discussed in Sec.3.2 (see also Abbott et al. 2017c) (Abbott et al. 2017c).These limits depend on the sensitivity of the detector network, which, in turn, varies over time and with sky location, and have been marginalized over errors introduced by detector calibration.For the ADI and the CSG models mentioned above, as well as for the other seven models used in the generic transient method search (see Sec. 3.2), we provide population median exclusion limits, D 90 , in Table 1.These vary roughly over one order of magnitude, which reflects the wide range of models used in the analysis.We report the D 90 values found for each GRB in the case of ADI model A simulated signals and CSG simulated signals with central frequency of 150 Hz in Table 2, at the end of this paper.For each event in the O3a sample that was localized with an error radius smaller than 0 • .5, we searched GLADE (Dálya et al. 2018) for galaxies within 200 Mpc.We then compared the angular separation between each GRB and galaxy, and recorded all separations less than or equal to twice the error radius for each GRB.Of the 141 events in our sample, four had nearby galaxies according to the definition above: GRB 190530430, GRB 190531840, GRB 190610A, and GRB 190731943.Data for our GW follow-up analysis were available only in the case of the short GRB 190610A, first observed by Swift/BAT (Evans et al. 2019) and localized to within a 90% error radius of 1 .9(Lien et al. 2016;Palmer et al. 2019).On the edge of its localization region, there is a nearby galaxy at a luminosity distance of approximately 165 Mpc (z = 0.037), as reported in GLADE (see Fig. 5). 212The angular separation between the center of the localization region and the nearby galaxy is at the 2.21σ level relative to the formal fit error, which is slightly less conservative than the quoted 90% localization derived from the S/N, and is consistent with expectations of angular offsets from a host galaxy at that distance (Fong & Berger 2013).
We did not find any GW signal associated with GRB 190610A in the data available from the two LIGO detectors (Virgo data were not in observing mode at that particular time).Our modeled search described in Sec.3.1, which uses an on-source window from −5 s to +1 s around the GRB trigger time, placed 90% confidence exclusion distances of 63 Mpc, 82 Mpc, and 114 Mpc for BNS binaries with generically oriented spins, NSBH binaries with generically oriented spins, and aligned spin NSBH binaries (see Sec. 3.1 for more details on these three populations).In general, a distance of 165 Mpc can be within the reach of our modeled search, but GRB 190610A was in a sky location such that the sensitivity of both detectors was less than 30% of what it would have been in an optimal sky location.

CONCLUSIONS
We carried out targeted analyses for GWs associated with Fermi and Swift GRBs reported during the O3a LIGO-Virgo observing run.In the case of short and ambiguous GRBs events (see Sec. 2), we ran a modeled search for NS binary merger signals (Harry & Fairhurst 2011;Williamson et al. 2014), while an unmodeled search for GW transient signals was performed for all GRBs (Sutton et al. 2010;Was et al. 2012).As a result of our analyses, we found no GW signal in association with the GRBs that we followed up.This is consistent with the previously predicted rate of coincident detections of 0.1-1.4 per year for the third observing run of Advanced LIGO and Advanced Virgo (Abbott et al. 2017a).Additionally, by carrying out a weighted binomial test, we found no strong evidence for a population of unidentified subthreshold GW signals in our results.We set lower bounds on the distances to the progenitors of all GRBs we analyzed for a number of emission models.These D 90 values are reported in Table 2, along with other information about each GRB that we considered; this includes timing, sky location, observing instrument, and GW detectors with available data.The 90% confidence level exclusion distances achieved in this run include the largest values published so far for some individual GRBs (see Abbott et al. 2017cAbbott et al. , 2019b)).Among the GRBs we analyzed is GRB 190610A, the sky localization of which included a nearby galaxy at a luminosity distance of 165 Mpc.We placed 90% confidence level exclusion distances lower than this value for NS binary merger GW signals and are therefore unable to rule out the possibility that GRB 190610A occurred in such galaxy.
the support of the NSF, STFC, INFN and CNRS for provision of computational resources.
We would like to thank all of the essential workers who put their health at risk during the COVID-19 pandemic, without whom we would not have been able to complete this work.
We would also like to thank Christian Malacaria and Aaron Tohuvavohu for providing useful comments that helped improve this paper.
Table 2. Information and Limits on Associated GW Emission for Each of the Fermi and Swift GRBs Followed Up during the LIGO-Virgo Run O3a.

Figure 1 .
Figure1.Cumulative distribution of loudest on-source event p-values for the neutron star binary modeled search in O3a.If the search reports no trigger in the on-source, we plot an upper limit on the p-value of 1 (open circles), and a lower limit equal to the fraction of off-source trials that contained no trigger (full circles).The dashed line indicates the expected distribution of p-values under the no-signal hypothesis, with the corresponding 2σ envelope marked by dotted lines.

Figure 2 .
Figure 2. Cumulative distribution of p-values from the unmodeled search for transient gravitational waves associated with 105 gamma-ray bursts.The dashed line represents the expected distribution under the no-signal hypothesis, with dotted lines indicating a 2σ deviation from this distribution.
Figure3.Cumulative histograms of the 90% confidence exclusion distances, D90, for the binary neutron star (blue, thin line) and generically spinning neutron star-black hole (orange, thick line) signal models, shown for the sample of 32 short and ambiguous gamma-ray bursts that were followed up by the NS binary modeled search during O3a, none of which had an identified gravitational wave counterpart.For a given GRB event and signal model, D90 is the distance within which 90% of simulated signals inserted into off-source data are recovered with greater significance than the most significant on-source trigger.These simulated signals have inclinations θJN -the angle between the total angular momentum and the line of sight -drawn uniformly in sin θJN with θJN restricted to [0 • , 30 • ] ∪ [150 • , 180 • ].
Figure 4. Cumulative histograms of the 90% confidence exclusion distances, D90, for accretion disk instability (ADI; van Putten 2001; van Putten et al. 2014) signal model A (orange, thin line) and circular sine-Gaussian (CSG) 150 Hz (Abbott et al. 2017c) model (green, thick line).For a given GRB and signal model this is the distance within which 90% of simulated signals inserted into off-source data are successfully recovered with a significance greater than the loudest on-source trigger.The median values for ADI-A and CSG-150 waveforms are 23 Mpc and 73 Mpc respectively.
Figure 4 shows the distribution of D 90 values for the ADI model A (van Putten 2001; van Putten et al. 2014) and for a CSG with central frequency of 150 Hz

Figure 5 .
Figure 5. Overlay of the estimated 90% Swift/BAT error radius for GRB 190610A (orange circle) on the sky.A galaxy at around 165 Mpc (Dálya et al. 2018) compatible with this localization is indicated by the blue crosshair.
The Satellite column lists the instrument the sky localization of which was used for GW analysis purposes.The Network column lists the GW detector network used in the analysis of each GRB: H1 = LIGO Hanford, L1 = LIGO Livingston, V1 = Virgo.A † denotes cases in which T90 > 60 s, so the on-source window of the generic transient search was extended to cover the GRB duration.For cases in which the generic transient search (Sec.3.2) and the neutron star binary search (Sec.3.1) used a different network, we report the network used by the latter in parentheses.Columns 8-12 display the 90% confidence exclusion distances to the GRB (D90) for several emission scenarios: BNS, generic and aligned spin NSBH, ADI-A, and CSG GW burst at 150 Hz with total radiated energy EGW = 10 −2 M c 2 .The first three are determined with the neutron star binary search, while the last two are calculated with the generic transient search.