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

A Decade of Gamma-Ray Bursts Observed by Fermi-LAT: The Second GRB Catalog

M. Ajello<fnref n="1" />, M. Arimoto<fnref n="2" />, M. Axelsson<fnref n="3" />,<fnref n="4" /> ID, L. Baldini<fnref n="5" />, G. Barbiellini<fnref n="6" />,<fnref n="7" />, D. Bastieri<fnref n="8" />,<fnref n="9" />, R. Bellazzini<fnref n="10" />, P. N. Bhat<fnref n="11" />, E. Bissaldi<fnref n="12" />,<fnref n="13" /> ID, R. D. Blandford<fnref n="14" />, R. Bonino<fnref n="15" />,<fnref n="16" />, J. Bonnell<fnref n="17" />,<fnref n="18" />, E. Bottacini<fnref n="14" />,<fnref n="19" />, J. Bregeon<fnref n="20" />, P. Bruel<fnref n="21" />, R. Buehler<fnref n="22" />, R. A. Cameron<fnref n="14" />, R. Caputo<fnref n="23" />, P. A. Caraveo<fnref n="24" />, E. Cavazzuti<fnref n="25" />, S. Chen<fnref n="8" />,<fnref n="19" />, C. C. Cheung<fnref n="26" />, G. Chiaro<fnref n="24" />, S. Ciprini<fnref n="27" />,<fnref n="28" />, D. Costantin<fnref n="29" />, M. Crnogorcevic<fnref n="18" />, S. Cutini<fnref n="30" />, M. Dainotti<fnref n="14" />, F. D'Ammando<fnref n="31" />,<fnref n="32" />, P. de la Torre Luque<fnref n="12" />, F. de Palma<fnref n="15" />, A. Desai<fnref n="1" />, R. Desiante<fnref n="15" />, N. Di Lalla<fnref n="5" />, L. Di Venere<fnref n="12" />,<fnref n="13" />, F. Fana Dirirsa<fnref n="33" />, S. J. Fegan<fnref n="21" />, A. Franckowiak<fnref n="22" />, Y. Fukazawa<fnref n="34" />, S. Funk<fnref n="35" />, P. Fusco<fnref n="12" />,<fnref n="13" />, F. Gargano<fnref n="13" />, D. Gasparrini<fnref n="28" />,<fnref n="30" />, N. Giglietto<fnref n="12" />,<fnref n="13" />, F. Giordano<fnref n="12" />,<fnref n="13" />, M. Giroletti<fnref n="31" />, D. Green<fnref n="36" />, I. A. Grenier<fnref n="37" />, J. E. Grove<fnref n="26" />, S. Guiriec<fnref n="17" />,<fnref n="38" />, E. Hays<fnref n="17" />, J. W. Hewitt<fnref n="39" /> ID, D. Horan<fnref n="21" />, G. Jóhannesson<fnref n="40" />,<fnref n="41" />, D. Kocevski<fnref n="17" />, M. Kuss<fnref n="10" />, L. Latronico<fnref n="15" />, J. Li<fnref n="22" />, F. Longo<fnref n="6" />,<fnref n="7" />, F. Loparco<fnref n="12" />,<fnref n="13" />, M. N. Lovellette<fnref n="26" />, P. Lubrano<fnref n="30" />, S. Maldera<fnref n="15" /> ID, A. Manfreda<fnref n="5" />, G. Martí-Devesa<fnref n="42" />, M. N. Mazziotta<fnref n="13" />, I. Mereu<fnref n="43" />, M. Meyer<fnref n="14" />, P. F. Michelson<fnref n="14" />, N. Mirabal<fnref n="17" />,<fnref n="44" />, W. Mitthumsiri<fnref n="45" />, T. Mizuno<fnref n="46" />, M. E. Monzani<fnref n="14" />, E. Moretti<fnref n="47" />, A. Morselli<fnref n="27" />, I. V. Moskalenko<fnref n="14" />, M. Negro<fnref n="15" />,<fnref n="16" />, E. Nuss<fnref n="20" />, M. Ohno<fnref n="34" />, N. Omodei<fnref n="14" /> ID, M. Orienti<fnref n="13" />, E. Orlando<fnref n="14" />, M. Palatiello<fnref n="6" />,<fnref n="7" />, V. S. Paliya<fnref n="22" />, D. Paneque<fnref n="36" />, M. Persic<fnref n="6" />,<fnref n="48" />, M. Pesce-Rollins<fnref n="10" />, V. Petrosian<fnref n="14" />, F. Piron<fnref n="20" />, S. Poolakkil<fnref n="11" />, H. Poon<fnref n="34" />, T. A. Porter<fnref n="14" />, G. Principe<fnref n="35" />, J. L. Racusin<fnref n="17" />, S. Rainò<fnref n="12" />,<fnref n="13" />, R. Rando<fnref n="8" />,<fnref n="9" />, M. Razzano<fnref n="10" />,⁶¹, S. Razzaque<fnref n="33" />, A. Reimer<fnref n="14" />,<fnref n="42" />, O. Reimer<fnref n="14" />,<fnref n="42" />, T. Reposeur<fnref n="49" />, F. Ryde<fnref n="4" />,<fnref n="50" />, D. Serini<fnref n="12" />, C. Sgrò<fnref n="10" />, E. J. Siskind⁵¹, E. Sonbas⁵², G. Spandre<fnref n="10" />, P. Spinelli<fnref n="12" />,<fnref n="13" />, D. J. Suson⁵³, H. Tajima<fnref n="14" />,⁵⁴, M. Takahashi<fnref n="36" />, D. Tak<fnref n="17" />,<fnref n="18" />, J. B. Thayer<fnref n="14" />, D. F. Torres⁵⁵,⁵⁶, E. Troja<fnref n="17" />,<fnref n="18" />, J. Valverde<fnref n="21" />, P. Veres<fnref n="11" />, G. Vianello<fnref n="14" /> ID, A. von Kienlin⁵⁷, K. Wood⁵⁸, M. Yassine<fnref n="6" />,<fnref n="7" />, S. Zhu⁵⁹, and S. Zimmer<fnref n="42" />,⁶⁰

<footnote n="1">¹Department of Physics and Astronomy, Clemson University, Kinard Lab of Physics, Clemson, SC 29634-0978, USA</footnote>

<footnote n="2">²Faculty of Mathematics and Physics, Institute of Science and Engineering, Kanazawa University, Kakuma, Kanazawa, Ishikawa 920-1192, Japan</footnote>

<footnote n="3">³Department of Physics and Oskar Klein Center for Cosmoparticle Physics, Stockholm University, SE-106 91 Stockholm, Sweden; magnusa@astro.su.se</footnote>

<footnote n="4">⁴Department of Physics, KTH Royal Institute of Technology, AlbaNova, SE-106 91 Stockholm, Sweden</footnote>

<footnote n="5">⁵Università di Pisa and Istituto Nazionale di Fisica Nucleare, Sezione di Pisa I-56127 Pisa, Italy</footnote>

<footnote n="6">⁶Istituto Nazionale di Fisica Nucleare, Sezione di Trieste, I-34127 Trieste, Italy</footnote>

<footnote n="7">⁷Dipartimento di Fisica, Università di Trieste, I-34127 Trieste, Italy</footnote>

<footnote n="8">⁸Istituto Nazionale di Fisica Nucleare, Sezione di Padova, I-35131 Padova, Italy</footnote>

<footnote n="9">⁹Dipartimento di Fisica e Astronomia “G. Galilei,” Università di Padova, I-35131 Padova, Italy</footnote>

<footnote n="10">¹⁰Istituto Nazionale di Fisica Nucleare, Sezione di Pisa, I-56127 Pisa, Italy</footnote>

<footnote n="11">¹¹Center for Space Plasma and Aeronomic Research (CSPAR), University of Alabama in Huntsville, Huntsville, AL 35899, USA</footnote>

<footnote n="12">¹²Dipartimento di Fisica “M. Merlin” dell’Università e del Politecnico di Bari, I-70126 Bari, Italy; elisabetta.bissaldi@ba.infn.it</footnote>

<footnote n="13">¹³Istituto Nazionale di Fisica Nucleare, Sezione di Bari, I-70126 Bari, Italy</footnote>

<footnote n="14">¹⁴W.W. Hansen Experimental Physics Laboratory, Kavli Institute for Particle Astrophysics and Cosmology, Department of Physics and SLAC National Accelerator Laboratory, Stanford University, Stanford, CA 94305, USA; nicola.omodei@stanford.edu, giacomov@slac.stanford.edu</footnote>

<footnote n="15">¹⁵Istituto Nazionale di Fisica Nucleare, Sezione di Torino, I-10125 Torino, Italy</footnote>

<footnote n="16">¹⁶Dipartimento di Fisica, Università degli Studi di Torino, I-10125 Torino, Italy</footnote>

<footnote n="17">¹⁷NASA Goddard Space Flight Center, Greenbelt, MD 20771, USA</footnote>

<footnote n="18">¹⁸Department of Astronomy, University of Maryland, College Park, MD 20742, USA</footnote>

<footnote n="19">¹⁹Department of Physics and Astronomy, University of Padova, Vicolo Osservatorio 3, I-35122 Padova, Italy</footnote>

<footnote n="20">²⁰Laboratoire Univers et Particules de Montpellier, Université Montpellier, CNRS/IN2P3, F-34095 Montpellier, France</footnote>

<footnote n="21">²¹Laboratoire Leprince-Ringuet, École polytechnique, CNRS/IN2P3, F-91128 Palaiseau, France</footnote>

<footnote n="22">²²Deutsches Elektronen Synchrotron DESY, D-15738 Zeuthen, Germany</footnote>

<footnote n="23">²³Center for Research and Exploration in Space Science and Technology (CRESST) and NASA Goddard Space Flight Center, Greenbelt, MD 20771, USA</footnote>

<footnote n="24">²⁴INAF-Istituto di Astrofisica Spaziale e Fisica Cosmica Milano, via E. Bassini 15, I-20133 Milano, Italy</footnote>

<footnote n="25">²⁵Italian Space Agency, Via del Politecnico snc, I-00133 Roma, Italy</footnote>

<footnote n="26">²⁶Space Science Division, Naval Research Laboratory, Washington, DC 20375-5352, USA</footnote>

<footnote n="27">²⁷Istituto Nazionale di Fisica Nucleare, Sezione di Roma “Tor Vergata,” I-00133 Roma, Italy</footnote>

<footnote n="28">²⁸Space Science Data Center—Agenzia Spaziale Italiana, Via del Politecnico, snc, I-00133, Roma, Italy</footnote>

<footnote n="29">²⁹University of Padua, Department of Statistical Science, Via 8 Febbraio, 2, I-35122 Padova, Italy</footnote>

<footnote n="30">³⁰Istituto Nazionale di Fisica Nucleare, Sezione di Perugia, I-06123 Perugia, Italy</footnote>

<footnote n="31">³¹INAF Istituto di Radioastronomia, I-40129 Bologna, Italy</footnote>

<footnote n="32">³²Dipartimento di Astronomia, Università di Bologna, I-40127 Bologna, Italy</footnote>

<footnote n="33">³³Department of Physics, University of Johannesburg, P.O. Box 524, Auckland Park 2006, South Africa</footnote>

<footnote n="34">³⁴Department of Physical Sciences, Hiroshima University, Higashi-Hiroshima 739-8526, Japan</footnote>

<footnote n="35">³⁵Friedrich-Alexander-Universität Erlangen-Nürnberg, Erlangen Centre for Astroparticle Physics, Erwin-Rommel-Str. 1, D-91058 Erlangen, Germany</footnote>

<footnote n="36">³⁶Max-Planck-Institut für Physik, D-80805 München, Germany</footnote>

<footnote n="37">³⁷AIM, CEA, CNRS, Université Paris-Saclay, Université Paris Diderot, Sorbonne Paris Cité, F-91191 Gif-sur-Yvette, France</footnote>

<footnote n="38">³⁸The George Washington University, Department of Physics, 725 21st St., NW, Washington, DC 20052, USA</footnote>

<footnote n="39">³⁹University of North Florida, Department of Physics, 1 UNF Drive, Jacksonville, FL 32224, USA</footnote>

<footnote n="40">⁴⁰Science Institute, University of Iceland, IS-107 Reykjavik, Iceland</footnote>

<footnote n="41">⁴¹Nordita, Royal Institute of Technology and Stockholm University, Roslagstullsbacken 23, SE-106 91 Stockholm, Sweden</footnote>

<footnote n="42">⁴²Institut für Astro- und Teilchenphysik und Institut für Theoretische Physik, Leopold-Franzens-Universität Innsbruck, A-6020 Innsbruck, Austria</footnote>

<footnote n="43">⁴³Dipartimento di Fisica, Università degli Studi di Perugia, I-06123 Perugia, Italy</footnote>

<footnote n="44">⁴⁴Department of Physics and Center for Space Sciences and Technology, University of Maryland Baltimore County, Baltimore, MD 21250, USA</footnote>

<footnote n="45">⁴⁵Department of Physics, Faculty of Science, Mahidol University, Bangkok 10400, Thailand</footnote>

<footnote n="46">⁴⁶Hiroshima Astrophysical Science Center, Hiroshima University, Higashi-Hiroshima 739-8526, Japan</footnote>

<footnote n="47">⁴⁷Institut de Física d'Altes Energies (IFAE), Edifici Cn, Universitat Autònoma de Barcelona (UAB), E-08193 Bellaterra (Barcelona), Spain</footnote>

<footnote n="48">⁴⁸Osservatorio Astronomico di Trieste, Istituto Nazionale di Astrofisica, I-34143 Trieste, Italy</footnote>

<footnote n="49">⁴⁹Centre d'Études Nucléaires de Bordeaux Gradignan, IN2P3/CNRS, Université Bordeaux 1, BP120, F-33175 Gradignan Cedex, France</footnote>

<footnote n="50">⁵⁰The Oskar Klein Centre for Cosmoparticle Physics, AlbaNova, SE-106 91 Stockholm, Sweden</footnote>
Abstract

The Large Area Telescope (LAT) aboard the Fermi spacecraft routinely observes high-energy emission from gamma-ray bursts (GRBs). Here we present the second catalog of LAT-detected GRBs, covering the first 10 yr of operations, from 2008 to 2018 August 4. A total of 186 GRBs are found; of these, 91 show emission in the range 30–100 MeV (17 of which are seen only in this band) and 169 are detected above 100 MeV. Most of these sources were discovered by other instruments (Fermi/GBM, Swift/BAT, AGILE, INTEGRAL) or reported by the Interplanetary Network (IPN); the LAT has independently triggered on four GRBs. This catalog presents the results for all 186 GRBs. We study onset, duration, and temporal properties of each GRB, as well as spectral characteristics in the 100 MeV–100 GeV energy range. Particular attention is given to the photons with the highest energy. Compared with the first LAT GRB catalog, our rate of detection is significantly improved. The results generally confirm the main findings of the first catalog: the LAT primarily detects the brightest GBM bursts, and the high-energy emission shows delayed onset as well as longer duration. However, in this work we find delays exceeding 1 ks and several GRBs with durations over 10 ks. Furthermore, the larger number of LAT detections shows that these GRBs not only cover the high-fluence range of GBM-detected GRBs but also sample lower fluences. In addition, the greater number of detected GRBs with redshift estimates allows us to study their properties in both the observer and rest frames. Comparison of the observational results with theoretical predictions reveals that no model is currently able to explain all results, highlighting the role of LAT observations in driving theoretical models.

Key words: catalogs – gamma-ray burst – general – methods: data analysis

Supporting material: FITS file

1. Introduction

Observations by the Fermi Gamma-Ray Space Telescope in the 10 yr since it was placed into orbit on 2008 June 11 have allowed for the opportunity to study the broadband properties of gamma-ray bursts (GRBs) over an unprecedented energy range. Its two scientific instruments, the Large Area Telescope (LAT; Atwood et al. 2009) and the Gamma-Ray Burst Monitor (GBM; Meegan et al. 2009), provide the combined capability of probing emission from GRBs over seven decades in energy. These groundbreaking observations have helped to characterize the highest-energy emission from these events (Abdo et al. 2009a, 2009e, 2010a; Ackermann et al. 2010a, 2010b, 2011), furthered our understanding of the emission processes associated with GRBs (Ryde et al. 2010; Axelsson et al. 2012; Preece et al. 2014; Ahlgren et al. 2015; Arimoto et al. 2016; Burgess et al. 2016; Moretti & Axelsson 2016), helped place constraints on the Lorentz invariance of the speed of light (Abdo et al. 2009b, 2009c; Shao et al. 2010; Nemiroff et al. 2012; Vasileiou et al. 2013) and the gamma-ray opacity of the universe (Abdo et al. 2010b; Desai et al. 2017), and motivated revisions in our basic understanding of collisionless relativistic shock physics (Ackermann et al. 2014).

The first Fermi-LAT GRB catalog was published in 2013 (Ackermann et al. 2013b, hereafter 1FLGC). It is a compilation of the 35 GRBs, including 30 long-duration (>2 s) and 5 short-duration (<2 s) GRBs, detected in the period 2008 August

through 2011 July (the first GRB included is GRB 080825C, and the last is GRB 110731A). Of these GRBs, 28 were detected with standard analysis techniques at energies above 100 MeV, while 7 GRBs were detected at only <100 MeV using the LAT Low-Energy (LLE) technique. It established three main observational characteristics of the high-energy emission:

1. Additional spectral components: Many of the bright GRBs observed by Fermi-LAT were not well fit with a single Band function (Band et al. 1993). In particular, an additional power-law (PL) component was required to account for the high-energy data of four bursts.

2. Delayed onset: The emission above 100 MeV was systematically delayed compared with the emission seen at lower energies, in the keV–MeV energy range. Delays of up to 40 s have been detected, with a few seconds being a typical value.

3. Extended duration: The emission above 100 MeV also systematically lasted longer than the keV–MeV prompt emission. The flux generally followed a PL decay with time, $F \propto t^{-\alpha}$, with $\alpha$ close to -1.

The 1FLGC also left some open questions regarding GRB properties at high energy, which we plan to address in the 2FLGC:

1. Hyperfluent GRBs: Four GRBs of the 1FLGC hinted at the possibility of a different class of high-energy fluence, significantly greater than the average fluence of the other GRBs. Because of the small sample, the hyperfluent class of GRBs was not significant and its confirmation was left for subsequent observations.

<footnote n="51">51 NYCB Real-Time Computing Inc., Lattingtown, NY 11560-1025, USA</footnote>

<footnote n="52">52 Adiyaman University, 02040 Adiyaman, Turkey</footnote>

<footnote n="53">53 Purdue University Northwest, Hammond, IN 46323, USA</footnote>

<footnote n="54">54 Solar-Terrestrial Environment Laboratory, Nagoya University, Nagoya 464-8601, Japan</footnote>

<footnote n="55">55 Institute of Space Sciences (CSICIEEC), Campus UAB, Carrer de Magrans s/n, E-08193 Barcelona, Spain</footnote>

<footnote n="56">56 Institució Catalana de Recerca i Estudis Avançats (ICREA), E-08010 Barcelona, Spain</footnote>

<footnote n="57">57 Max-Planck Institut für extraterrestrische Physik, D-85748 Garching, Germany</footnote>

<footnote n="58">58 Praxis Inc., Alexandria, VA 22303, resident at Naval Research Laboratory, Washington, DC 20375, USA</footnote>

<footnote n="59">59 Albert-Einstein-Institut, Max-Planck-Institut für Gravitationsphysik, D-30167 Hannover, Germany</footnote>

<footnote n="60">60 University of Geneva, Département de physique nucléaire et corpusculaire (DPNC), 24 quai Ernest-Ansermet, CH-1211 Genève 4, Switzerland
Received 2019 March 28; accepted 2019 April 25; published 2019 June 13</footnote>

<footnote n="61">61 Funded by contract FIRB-2012-RBFR12PM1F from the Italian Ministry of Education, University and Research (MIUR).</footnote>
2. Late-time temporal decay index: The distribution of the late-time temporal decay index in the 1FLGC was clustered around -1, supporting the hypothesis of an adiabatically expanding fireball as a common origin of the extended emission. Although this scenario could explain all observations, only nine GRBs had enough data to allow the decay index to be determined.

3. Breaks in the late-time light curve: Three GRBs in the 1FLGC showed breaks in the temporal decay at late times, similar to the breaks observed in the X-ray light curves. Exploration of this feature was limited by the small sample and the relatively low significance of the breaks.

Since the publication of the 1FLGC, a number of improvements have been made with regard to LAT data processing. The two major changes concern the development of a new event analysis. Since launch, the LAT event classes have undergone a number of versions (or “passes”), and the latest “Pass 8” analysis constitutes a major improvement on the previous versions. In addition, a new detection algorithm has been developed running in parallel over a range of different timescales. Taken together, this has increased the detection efficiency by over 60% (Vianello et al. 2015), and in particular allows the detection of fainter high-energy GRB counterparts.

This paper presents the second catalog of GRBs detected by the *Fermi*–LAT (2FLGC), covering a 10 yr period, from 2008 August to 2018 August. During this time, the GBM triggered on 2357 GRBs, approximately half of which were in the field of view (FOV) of the LAT at the time of trigger. LAT counterparts are searched for in ground processing, following external triggers provided by the GBM, as well as by other instruments. The LAT instrument is also capable of detecting GRBs through an onboard trigger search algorithm. This is a very rare occurrence and has so far only happened 4 times (see Section 4). In addition, continual blind searches are performed as part of the standard ground processing, to look for untriggered events. These efforts are further described in Ackermann et al. (2016) and Ajello et al. (2018).

In the 2FLGC, we have performed an entirely new, standardized analysis to look for LAT counterparts to all GRB triggers reported during the first decade of *Fermi* operations. In Section 2 we describe the data used in this study, giving first a short description of the *Fermi* instruments, followed by the description of data cuts and sample selection. This is followed in Section 3 by a detailed description of the analysis methods used for detection and localization. We also present the methodology that we followed to characterize the temporal and spectral properties of the detected GRBs. The 2FLGC is focused on the LAT-only properties of the bursts in our sample, and we do not perform any joint spectral analysis of LAT and GBM data. In Sections 4 and 5 we present and discuss our results. Finally, we examine the theoretical implications of our results in Section 6.

## 2. Data Preparation

The 2FLGC presents analysis done with data from the LAT, i.e., covering energies above 30 MeV for GRBs detected through LLE, and from 100 MeV to 300 GeV for GRBs detected with the standard analysis. However, the GBM provides the vast majority of the GRB triggers and is an integral part of the GRB observations made by *Fermi*. We

therefore begin with a brief description of both instruments. This is followed by a more thorough description of the LAT data used for the analysis. Finally, we present the selection of GRB triggers that formed the seed of the catalog.

### 2.1. Instrument Overview

The LAT is a pair production telescope sensitive to γ-rays in the energy range from ~30 MeV to more than 300 GeV. The instrument and its on-orbit calibrations are described in detail in Atwood et al. (2009) and Abdo et al. (2009d). The telescope consists of a 4 × 4 array of identical towers, each including a tracker of 18 x-y silicon strip detector planes interleaved with tungsten foils, followed by an 8.6 radiation length imaging calorimeter made with CsI(Tl) scintillation crystals with a hodoscopic layout. This array is surrounded by a segmented anticoincidence detector made of 89 plastic scintillator tiles that identifies and rejects charged particle background events with an efficiency above 99.97% (Ackermann et al. 2012a).

Whether or not an event is observable by the LAT is primarily defined by two angles: the angle $\zeta$ with respect to the spacecraft zenith, and the viewing angle $\theta$ from the LAT boresight. The LAT performance—including the dependence of the effective area on energy and $\theta$—is presented on the official *Fermi*–LAT performance web page.<fnref n="62" /> In the analysis performed in this catalog, we do not make any explicit cuts on the angle $\theta$; however, the exposure will drop very quickly for $\theta$ greater than ~75°. The wide FOV (~2.4 sr at 1 GeV) of the LAT, its high observing efficiency (obtained by keeping the FOV on the sky with scanning observations), its broad energy range, its large effective area, its low dead time per event (~27 µs), its efficient background rejection, and its good angular resolution (the 68% containment angle of the point-spread function is 0.8 at 1 GeV) are vastly improved in comparison with those of previous instruments such as EGRET (Esposito et al. 1999). As a result, the LAT provides more GRB detections, higher statistics per detection, and more accurate localizations.

The GBM is composed of 12 sodium iodide (NaI) and two bismuth germanate (BGO) detectors sensitive in the 8 keV–1 MeV and 250 keV–40 MeV energy ranges, respectively. The NaI detectors are arranged in groups of three at each of the four edges of the spacecraft, and the two BGO detectors are placed symmetrically on opposite sides of the spacecraft, resulting in an FOV of 9.5 sr. Triggering and localization are determined from the NaI detectors, while spectroscopy is performed using both the NaI and BGO detectors. The GBM flight software continually monitors the detector rates and triggers when a statistically significant rate increase occurs in two or more NaI detectors. A combination of 28 timescales and energy ranges are currently tested, with the first combination that exceeds the predefined threshold (generally $4.5\sigma$) being considered the triggering timescale. Localization is performed using the relative event rates of detectors with different orientations with respect to the source and is typically accurate to a few degrees (statistical uncertainty). An additional systematic uncertainty has been characterized as a core-plus-tail model, with 90% of GRBs having a 3.7° uncertainty and a small tail suffering a larger than 10° systematic uncertainty (Connaughton et al. 2015). The GBM covers roughly four decades in energy and provides a bridge from the low energies (below ~1 MeV),

<footnote n="62">⁶² http://www.slac.stanford.edu/exp/glast/groups/canda/lat_Performance.htm</footnote>
where most of the GRB emission takes place, to the less studied energy range that is accessible to the LAT.

For GRBs exceeding a preset threshold for peak flux or fluence in the GBM, an autonomous repoint request (ARR) is sent to the spacecraft flight software. If the GBM request is accepted, a special LAT observation mode is initiated. This will keep the GBM flight software location in the LAT FOV for an extended period of time, typically ~2.5 hr, subject to observational constraints.

On 2018 March 16, one of the solar array drive assemblies on *Fermi* suffered a malfunction. This led to the LAT and GBM being switched off, and normal science operations were only resumed on April 13. However, as of writing one of the solar panels remains in a fixed position, and a modified rocking strategy has been adopted. As a result, the LAT rocks between the northern and southern sky every week, as opposed to every orbit as before. A further impact is that ARRs have been disabled.

## 2.2. Sample Selection

The 2FLGC presents the results of a search for high-energy counterparts of GRBs that triggered space instruments and have an available public localization. We considered in particular bursts detected by the GBM, the *Swift* Burst Alert Telescope (BAT; Barthelmy et al. 2005) on board the *Neil Gehrels Swift Observatory* (Gehrels et al. 2004), and the *International Gamma-Ray Astrophysics Laboratory* Soft Gamma-Ray Imager (INTEGRAL/ISGRI; Lebrun et al. 2003) on board the INTEGRAL satellite (Winkler et al. 2003), or reported by the Interplanetary Network (IPN<sup>63</sup>). In the 10 yr period covered by the 2FLGC, there were 2357 GBM GRB triggers, 876 *Swift*/BAT GRB triggers, 65 INTEGRAL/ISGRI GRB triggers, and 83 events reported by the IPN (a few through private communication by the PI K. Hurley). We also considered seven bursts contained in the first catalog of GRBs detected by the *Astrorivelatore Gamma a Immagini Leggero* mini-calorimeter (AGILE/MCAL; Galli et al. 2013) and not contained in any other list. After accounting for bursts that triggered more than one instrument, we have a total of 3044 independent GRBs. We use the localizations provided by the GBM, unless a better localization is provided by one of the instruments on board *Swift*, namely, the BAT, the X-ray Telescope (XRT; Burrows et al. 2005), or the UV-Optical Telescope (UVOT; Roming et al. 2005), or by the IPN. All these latter positions are distributed via the Gamma-Ray Burst Coordinates Network (GCN),<sup>64</sup> while the GBM-only localizations are reported in the *Fermi*-GBM online GRB Catalog<sup>65</sup> (hereafter FGGC; see Bhat et al. 2016, for more details).

## 2.3. LAT Data Cuts and Temporal Selection

For the standard analysis, we use Pass 8 LAT data with energies between 100 MeV and 100 GeV, selecting the time interval around each trigger from 600 s before to 100 ks after the GRB trigger time, and defining a standard region of interest (ROI) around the trigger location of 12°. In order to reduce the contamination of the Earth limb, in some dedicated cases (13 GRBs) we define a smaller ROI with a radius of 8°. It is worth

noting that as a final check we look for high-energy events over a larger energy range (up to 300 GeV), as discussed in Section 4.8.

We then use *gtmktime* to select only those time intervals when the center of the ROI has a zenith angle $\zeta < 97°$ (i.e., every point of the ROI has $\zeta < 105°$). For bursts with an initial value of $\zeta > 90°$, we increase our selection to include all time intervals when the center of the ROI has a zenith angle $\zeta < 102°$. This allows us to study the prompt emission of those GRBs that started close to the limb of the Earth. The choice of the event class depends on the timescale on which we detect the signal from the GRB and is described in Section 3.2.

## 2.4. LLE Data

The LLE technique is an analysis method designed to study bright transient phenomena, such as GRBs and solar flares, in the 30 MeV–1 GeV energy range. The LAT Collaboration developed this analysis using a different approach than the one used in the standard photon analysis, which is based on sophisticated classification procedures (a detailed description of the standard analysis can be found in Atwood et al. 2009; Ackermann et al. 2012a). The idea behind LLE is to maximize the effective area below ~1 GeV by relaxing the standard analysis requirement on background rejection.

The basic LLE selection is based on a few simple requirements on the event topology in the three LAT subdetectors. First of all, an event passing the LLE selection must have at least one reconstructed track in the tracker and therefore an estimate of the direction of the incoming photon. Second, we require that the reconstructed energy of the event be nonzero.

We use the information provided by the flight software in LLE to efficiently select events that are gamma-ray like. With the release of Pass 8 data, we have also improved the LLE selection. For events with an incident angle $\theta < 40°$, we require that no anticoincidence tiles are in “veto” condition (to suppress charged particle contamination), while for angles greater than 40°, we allow a maximum of two tiles in “veto” condition, but no tracker hits can be found in proximity of the anticoincidence hits. This condition helps prevent suppression of large incident angle gamma-rays due to secondary electrons or positrons interacting with the anticoincidence detector downstream. In order to reduce the number of photons originating from the Earth limb in our LLE sample, we only keep reconstructed events with a zenith angle $\zeta < 90°$ or $\zeta < 100°$ (depending on the location of the GRB). Finally, we explicitly include in the selection a cut on the ROI, i.e., the position in the sky of the transient source we are observing. In other words, the localization of the source is embedded in the event selection, and therefore for a given analysis the LLE data are tailored to a particular location in the sky. The response of the detector for the LLE class is encoded in a response matrix, which is generated using a dedicated Monte Carlo simulation for each GRB, and is saved in the standard HEASARC RMF file format.<sup>66</sup> LLE data and the relative response are made available for any transient signal (GRB or solar flare) detected with a significance above $4\sigma$ through the HEASARC FERMILLE website.<sup>67</sup>

<footnote n="63"><sup>63</sup> http://www.ssl.berkeley.edu/ipn3/</footnote>

<footnote n="64"><sup>64</sup> https://gcn.gsfc.nasa.gov/gcn3_archive.html</footnote>

<footnote n="65"><sup>65</sup> https://heasarc.gsfc.nasa.gov/W3Browse/fermi/fermigrbrst.html</footnote>

<footnote n="66"><sup>66</sup> Described at http://heasarc.gsfc.nasa.gov/docs/heasarc/caldb/docs/memos/cal_gen_92_002/cal_gen_92_002.html#Sec:RMF-format.</footnote>

<footnote n="67"><sup>67</sup> http://heasarc.gsfc.nasa.gov/W3Browse/fermi/fermille.html</footnote>
## 2.5. Low-energy Data Used for Comparisons

While we do not perform any joint spectral fitting with GBM data in this work, comparisons with the sample of GBM-detected GRBs are both highly interesting and inevitable in order to characterize our sample. For this purpose we use the official data from the FGGC.

In order to perform the comparisons, we use the standard GRB properties that characterize the GRB emission, in particular the onset time and duration (both calculated by GBM in the 50–300 keV energy range), values of peak fluxes ($F_P$, calculated in 1024 ms and 64 ms intervals for long and short GRBs, respectively), fluences ($F$, calculated over the time interval used in the GBM spectral analysis, which might not always be coincident with the burst duration), and the spectral parameters of the best-fit model derived by the GBM in the 10–1000 keV energy band. The best fit is determined from the four standard spectral models tested in the GBM time-integrated spectral catalog analysis (see Gruber et al. 2014, for more details), namely, the simple PL, the smoothly broken PL (SBPL), the phenomenological Band function (BAND; Band et al. 1993), and the Comptonized model (COMP), which is a subset of the Band function. The latter three models are characterized by a low-energy spectral index $\alpha$, by a high-energy spectral index $\beta$, which in the case of the COMP model goes to minus infinity, and by the energy $E_{\text{peak}}$, which describes the peak of the $\nu F_{\nu}$ distribution.

The classification of GRBs into long and short classes is primarily derived from the low-energy duration as measured by GBM and follows the standard rule that long and short bursts are longer and shorter than 2 s, respectively (Kouveliotou et al. 1993). For bursts that did not trigger GBM and are not included in the FGGC, we use the durations calculated by the Konus-Wind instrument in the 20 keV–5 MeV energy range, which have been published in GCN Circulars.

## 3. Analysis Methods and Procedures

As presented in Section 1, significant improvements have been made to the analysis techniques since the 1FLGC. The two developments with highest impact are the “Pass 8” event analysis and a redesigned detection algorithm. “Pass 8” has been rebuilt from the ground up with respect to previous versions (“Pass 6” was used for the first catalog), resulting in increased sensitivity. It is thoroughly described in Atwood et al. (2013a). In this section, we focus on the improvements to the triggered search and detection algorithms for GRBs used to produce the 2FLGC.

### 3.1. Analysis Sequence

The procedure followed to produce the catalog is summarized in Figure 1. Each individual step of the analysis is described in detail in the following subsections.

Each input trigger was first run through a detection algorithm as described in Section 3.2. All triggers that passed the criteria for detection in any time window were placed in the list of potential candidates for analysis. This list was manually inspected following the procedure in Section 3.5. Only detections retained after this inspection were passed to the analysis pipeline. Here, analyses were performed to determine a number of key properties for each GRB, such as onset, duration, and spectral parameters. The steps of this analysis are described in Section 3.6.

### 3.2. The LAT Transient Factory

The LAT Transient Factory (LTF) algorithm was introduced after the publication of the 1FLGC and is presented in Vianello et al. (2015). LTF has been running continuously ever since, looking in real time for GRB counterparts in LAT data. When compared to the old algorithm on the same data set as the first catalog, it returned 50% more GRBs. The adoption of “Pass 8” data has led to a further 10% improvement.

LTF is based on the application of an unbinned maximum likelihood technique. This analysis starts by selecting all the photons detected by the LAT above 100 MeV in a circular ROI with radius $r$ and center ($\alpha_{ROI}, \delta_{ROI}$) in a time window starting at the trigger time $T_0$. Then, the presence of a new point source at a position ($\alpha_t, \delta_t$) is tested by using the likelihood-ratio test (LRT). The null hypothesis for the test is represented by a baseline likelihood model including all point sources from the LAT source catalog (Acero et al. 2015) with all parameters fixed, as well as the Galactic and isotropic diffuse emission templates$^{68}$ provided by the Fermi-LAT Collaboration (Acero et al. 2016), with the normalizations left free to vary. The alternative hypothesis is represented by the baseline model plus the new point source (“test source”), modeled with a PL spectrum with free index and normalization. The LRT uses as the test statistic (TS) twice the logarithm of the maximum of the likelihood function for the alternative hypothesis ($L_1$) divided by the maximum of the likelihood function for the null hypothesis ($L_0$):

$$TS = 2(\log L_1 - \log L_0). \qquad (1)$$

Detailed instructions on how to perform an unbinned likelihood analysis using the Fermi Science Tools can be found on the Fermi website.$^{69}$

For each trigger, a search is performed on five time windows starting at the trigger time $T_0$ and ending 10, 100, 500, 4000, and 10,000 s after $T_0$. This selection slightly differs from the standard LTF real-time analysis, which consists of 10 searches running in parallel over time intervals logarithmically spaced from the trigger time to 10 ks after that, as stated in Vianello et al. (2015). For the 10 and 100 s time windows, we use the Pass 8 P8R2_TRANSIENT020E_V6 event class and the corresponding response functions; for the longer time windows, the event class P8R2_TRANSIENT010E_V6 is used. In each time window the LTF starts from the input coordinates and trigger time as measured by the triggering instrument (see Section 2.2) and performs the following steps:

1. Finding map: we consider an ROI with radius $r$ centered on the input position and a square grid of side $\delta x$ inscribed in the ROI with a spacing of 0.8°. The size of the grid is fixed according to the triggering instrument and its typical localization accuracy (statistical + systematic), as well as the typical size of the LAT PSF. Specifically, $\delta x = 32°$ for triggers localized by the GBM that are dominated by systematic uncertainties (Connaughton et al. 2015), $\delta x = 1.6°$ for Swift and INTEGRAL triggers, and $\delta x = 5°$ for IPN triggers. The radius of the ROI is chosen as $r = \delta x/2 + 12°$ in order to have enough data around each point in the grid for performing an LRT

<footnote>$^{68}$ gll_iem_v06.fits and iso_P8R3_SOURCE_V2.txt; https://fermi.gsfc.nasa.gov/ssc/data/access/lat/BackgroundModels.html</footnote>

<footnote>$^{69}$ https://fermi.gsfc.nasa.gov/ssc/data/analysis/scitools/likelihood_tutorial.html</footnote>
<figure><img src="image_2.png" /><figcaption>**Figure 1.** Flow diagram representation of the analysis pipeline. The LTF detection algorithm is highlighted in yellow, and the catalog analysis pipeline, which is executed after manually checking the candidate seeds, is highlighted in blue. The LLE analysis pipeline, which is executed on all input triggers, is highlighted in green. All input candidates are analyzed by the LTF/LLE pipelines; only those passing the detection criteria are retained. Inspection in the manual checks can also lead to rejection of further cases (see text for details).</figcaption></figure>

test (see below). In order to reduce the contamination from the Earth limb—a bright source of γ-rays—all events with $\zeta < 100°$ are filtered out. The effect of this selection is taken into account when computing the exposure by the tool *gtltcube*. We then use the LRT as described above to test for the presence of a source at each position of the grid having at least three photons within $10°$. The latter requirement filters out points without any photon cluster around them, in order to reduce the computational cost. The point $(\alpha_{max}, \delta_{max})$ in the grid providing the maximum of the TS is considered the best guess for the position of the new transient and marked for further analysis.

2. *TS and position refinement:* we consider an ROI centered on $(\alpha_{max}, \delta_{max})$ with a radius of $8°$, and we perform an LRT as described above considering only the time intervals within the time window when the border of the ROI is at a zenith angle smaller than $105°$ ("good time intervals"). This is a different way of reducing the contamination from the Earth limb that is more effective than the one used in the previous step, but it can only be applied on small ROIs. We then use the tool *gtfindsrc* to search for the maximum of the likelihood under the alternative hypothesis (i.e., when the test source is added to the model), varying the position of the test source and profiling the other free parameters. The position $(\alpha_1, \delta_1)$ yielding the maximum of the likelihood is considered the new putative position for the candidate counterpart.

3. *Final TS and position:* the previous step is repeated using an ROI centered on $(\alpha_1, \delta_1)$ that yields the final TS (TS<sub>GRB</sub>). The tool *gtfindsrc* is run again returning the final estimate of the localization uncertainty.

4. *Photon-by-photon assignment of probability:* we run the tool *gtsrcprob* using the final optimized likelihood model under the alternative hypothesis. This tool assigns to each detected photon the probability of belonging to the test source, i.e., to the candidate counterpart. We then measure the number of photons $N_p > 0.9$ having a probability larger than 90% of belonging to the candidate counterpart.

The final products of LTF are five sets of results, one for each time window. In order to consider a counterpart detected, we consider in particular TS<sub>GRB</sub> and $N_p>0.9$, as explained in the next section.

### 3.3. Detection Threshold and False Discovery Rate

A classic result from Chernoff (1954) states that under the null hypothesis the TS of a single LRT as applied in LTF is a random variable that is zero half of the time and is distributed as $\chi^2$ with 1 degree of freedom the other half. This result was confirmed by Monte Carlo simulation in Mattox et al. (1996). Under these circumstances, the significance of the detection (z-score) is $\sqrt{\text{TS}}$; thus, a threshold of TS = 25 corresponds to a detection >5σ for one LRT.

As described in the previous section, LTF consists of multiple LRT procedures and the trial factor needs to be taken into account. The effective number of trials for one time window is, however, difficult to determine because the trials are not independent. Furthermore, we also need to account for the number of time windows and for the number of triggers searched.

To account for the number of triggers searched, we use the procedure proposed by Benjamini & Hochberg (1995). It assumes independent trials and is simple: all the p-values $p_i$ for
all the $n_{tr}$ searches, with $i=1...n_{tr}$, are sorted in increasing order. We then find $k$ so that $p_k$ is the largest $p$-value where $p_k \le (k/n_{tr})\alpha$, where $\alpha$ is the error probability for one test. All the triggers with $i=1...k$ are considered detected. In practice, we first compute through Monte Carlo simulations the effective number of trials for one time window $n_{trw}$. The value of $n_{trw}$ is different depending on the instrument that generated the trigger, and it reflects the size of the finding map (see previous section), so it is larger for larger finding maps. We find $n_{trw} = 110$ for GBM triggers, $n_{trw} = 12.28$ for IPN triggers, and $n_{trw} = 1.25$ for Swift and INTEGRAL triggers. We then compute the post-trial $p$-value for one time window by using the binomial distribution as $p_w = 1 - (1 - p_{pr})^{n_{trw}}$, where $p_{pr}$ is the $p$-value coming from the LRT applied to the current time window. There is also another independent trial factor $n_{et}$, which we consider to be equal to the number of time windows where we effectively searched for a counterpart. Note that this is a slightly conservative approach, as the time windows are not independent and thus $n_{et}$ is in reality a little smaller. The number $n_{et}$ can vary from 0 to 5 depending on how many time windows had an exposure larger than zero after our data cuts. For example, if the trigger was never in the FOV or it was always at a zenith angle larger than our cut during a time window, this will not constitute a search, and it will not contribute to $n_{et}$. We can now compute the $p$-value for a GRB corrected for both the spatial and the time trials as $p = 1 - (1 - p_{w,min})^{n_{et}}$, where $p_{w,min}$ is the minimum $p_w$ corresponding to the maximum final TS$_{GRB}$ found by LTF among the timescales searched. We then apply the Benjamini & Hochberg (1995) procedure using these $p$-values and correct for the number of triggers searched as explained above.

We further apply the quality cut $n_{p>0.9} \ge 3$ to the list of detections, i.e., we require at least three photons with a probability larger than 90% of belonging to the GRB. This neutralizes the effect of isolated high-energy photons ($\ge 10$ GeV) within the search region that tend to return high TS values during the unbinned analysis but also very hard spectra, with photon indices close to 0. Moreover, three photons are required in order to have both the normalization and the photon index free during the likelihood maximization under the alternative hypothesis and still have at least 1 degree of freedom.

### 3.4. The Bayesian Blocks Burst Detection Algorithm for LLE Data

In order to detect GRB counterparts in LLE data, we use a counting analysis based on the well-known Bayesian blocks (BB) algorithm of Scargle et al. (2013). The BB algorithm is capable of dividing a time series into intervals of constant rate, opening a new block only when the rate of events changes in a statistically significant way. In particular, we use the unbinned version of the algorithm, which presents as the only parameter the probability $p_0$ of opening a new block when the rate is constant (false positive). However, before we can apply the BB algorithm, we need to introduce a preprocessing step to account for the time-varying background in LLE data. Otherwise, the BB algorithm will find many blocks following the variation in the event rate owing to the variations in the background.

We start by fitting a polynomial function to the data in two off-pulse time windows, one before and one after the trigger window, as shown in the left panel of Figure 2 for GRB 131014A. The trigger window is defined on the basis of the duration measured by the GBM in the 50–300 keV energy

range. Then, we exploit the fact that a nonuniform Poisson process with expected value $\lambda(t)$ can be converted to a uniform one by transforming the time reference system to

$$t' = \int_{0}^{t} \lambda(t) dt. \qquad (2)$$

In our case $\lambda(t)$ is the expected number of events of the time-varying background as modeled by the polynomial. We then transpose the time series of LLE events to the $t'$ time reference system and apply the standard BB algorithm. If more than one block is found in the trigger window, it means that there is a rate change on top of what is predicted by the background model. We then deem the transient detected, and we transform back to the original time reference system to yield the time interval of the detection. For all searches we use $p_0 = 10^{-3}$. The analysis is illustrated in the right panel of Figure 2, which shows the final BB representation of the light curve for GRB 131014A.

### 3.5. Final Manual Checks

The normalization of the Galactic and diffuse templates is allowed to vary in the analysis, but the parameters for the point sources are kept fixed to the values found in the updated Fermi-LAT point-source list.<fnref n="70" /> This means that the algorithm can potentially pick up non-GRB sources, such as AGN flares. Moreover, GRBs detected far from zenith might also suffer from Earth limb contamination. Therefore, a final manual check of all potential detections was also performed. The list of detections derived by the LTF pipeline previously described was divided into random subsets of bursts that were assigned for analysis to the members of the Fermi-LAT GRB team. Each putative event was independently cross-checked by two people, whose task was to either confirm or reject the detection. In case of agreement, the classification was seen as final; otherwise, the case was reviewed by a third person.

The manual checks included a series of tasks to be carried out. First, the LTF results were evaluated in each of the five temporal intervals, taking into account (1) the number of photons with probability $>90\%$, (2) the distance to the nearest known source, (3) the localization error, (4) the spectral index, and (5) the final TS value. We identified many “simple” cases, in which both the number of detected photons was high and the final TS was above 80 in several time intervals, and no other known sources were present in the ROI. These candidates were marked as confirmed with no further inspection.

Intermediate cases that needed deeper investigations included (1) cases where the final TS in all time intervals was close to the threshold, (2) cases where only three to four high-energy photons were detected, (3) cases where a bright source (an AGN, solar flares, etc.) was at an angle $<1°$ from the GRB candidate in the ROI, and (4) cases where a high TS value was obtained only by integrating over the longest timescale (from 0 to 10,000 s; see Section 3.2). In order to check for other active sources in the ROI around the time of the GRB trigger, we looked for flaring blazars using the publicly available FAVA tool<fnref n="71" /> (Ciprinì et al. 2013; Abdollahi et al. 2017), and for solar activity we checked the Solar Monitor public pages.<fnref n="72" />

<footnote n="70">⁷⁰ https://fermi.gsfc.nasa.gov/ssc/data/access/lat/fl8y/</footnote>

<footnote n="71">⁷¹ https://fermi.gsfc.nasa.gov/ssc/data/access/lat/FAVA/</footnote>

<footnote n="72">⁷² http://www.solarmonitor.org</footnote>
<figure><img src="image_3.png" /></figure>

<figure><img src="image_4.png" /><figcaption>**Figure 2.** Left panel: LLE light curve of GRB 131014 in the 30 MeV–1 GeV energy range. The polynomial fit to the background is marked by a red line. Right panel: BB representation of the light curve of GRB 131014. All light curves are centered around the GBM trigger time $T_0 = 403420143.2$ (in MET).</figcaption></figure>

<table><caption>**Table 1**

Definition of the Four Time Intervals (with Start and Stop Times) Used in the Time-integrated Spectral Analysis</caption><thead><tr><td>Name</td><td>Interval</td><td>Description</td></tr></thead><tbody><tr><td>GBM</td><td>[T<sub>GBM,05</sub>, T<sub>GBM,95</sub>]</td><td>GRB duration measured by GBM in the 50–300 keV energy range</td></tr><tr><td>LTF</td><td>[T<sub>0</sub>, T<sub>max</sub><sup>LTF</sup>]</td><td>Time interval showing the highest TS value as calculated by the LAT transient factory, starting from the GRB trigger time</td></tr><tr><td>LAT</td><td>[T<sub>LAT,0</sub>, T<sub>LAT,1</sub>]</td><td>GRB duration measured by LAT by performing a time-resolved likelihood analysis in the 100 MeV–100 GeV energy range</td></tr><tr><td>EXT</td><td>[T<sub>GBM,95</sub>, T<sub>LAT,1</sub>]</td><td>Interval from end of GBM to end of LAT duration</td></tr></tbody></table>

In case of particularly uncertain candidates, we performed an ad hoc likelihood analysis, similar to the one performed by the LTF pipeline, but running on dedicated time intervals that might differ from the catalog ones.

Through these manual checks, ~15%–20% of the examined cases were rejected as not connected to a GRB. The remaining events were processed in the dedicated analysis pipeline.

### 3.6. Catalog Analysis Description

In this section we describe the analysis steps we performed on each GRB of the final sample. The idea is to perform an automated analysis, which is implemented in a series of `python` scripts that are used to control the various steps. The analysis is based on `ScienceTools v11r05p03`, available for download at the Fermi Science Support Center.<sup>73</sup>

#### 3.6.1. Time-integrated Likelihood Analysis

We perform an unbinned likelihood analysis in four different time intervals. The “GBM” time interval represents the GRB duration as given by $T_{GBM,90}$ reported in the FGGC. $T_{GBM,90}$ is the interval during which the instrument measures from 5% to 95% of the total GRB flux in the 50–300 keV energy range (i.e., from $T_{GBM,05}$ to $T_{GBM,95}$). The “LTF” interval corresponds to the time interval in the LAT Transient Factory analysis where the highest TS was found. The “LAT” interval encompasses the signal detected by the LAT, as defined in **3.6.2. The “EXT” interval is defined as the time interval including LAT emission (if any) after the $T_{GBM,95}$. Table 1 summarizes the definition of these intervals.**

#### 3.6.2. Time-resolved Likelihood Analysis

In order to perform time-resolved likelihood analysis, we have developed an algorithm for adaptively binning the LAT events. Starting from the result of the analysis in the “LTF” time window, we apply `gtsrcprob` to calculate the probability of each LAT event to be associated with the GRB source. Starting with preselected logarithmically spaced time bins (48 bins from 0.01 s to 50 ks after the GBM trigger), we merge them until at least three events with probability >0.9 are present in each final bin. In practice, we have 3 degrees of freedom ($N_{dof}$): 2 associated with the PL describing the GRB, and 1 with the normalization of the isotropic diffuse component. The normalization of the Galactic model has been fixed to its nominal value (1). We require at least $N_{dof}$ events with probability >0.9 in every bin. In this way, we optimize the duration of the time intervals in order to always have enough photons to perform the fit. Once we have identified the time bins, we perform unbinned likelihood analysis in each bin, calculating the value of flux or, in the case of a TS value <10, we calculate the flux upper limit (95%) by profiling the likelihood function.

In the 1FLGC, the duration in the LAT was calculated based on the concept of $T_{GBM,90}$, i.e., the time during which 90% of the flux is collected. As the LAT observes each photon individually, this requires the simulation of light curves. In this analysis, we instead use a technique based on the individual photons intrinsic to the LAT. The total duration of the signal in the LAT, defined as $T_{LAT,100}$, is estimated starting from the results of the time-resolved analysis. The LAT onset time $T_{LAT,0}$ corresponds to the time when the first photon with probability $p > 0.9$ to be associated with the GRB is detected,

<footnote n="73">73 https://fermi.gsfc.nasa.gov/ssc/data/analysis/</footnote>
<figure><img src="image_5.png" /><figcaption>**Figure 3.** Temporal extended emission for two bright LAT GRBs, the long GRB 080916C (left panel) and the short GRB 090510 (right panel). Blue points show the flux in each time bin, while black arrows mark upper limits. The green vertical dashed lines indicate the first and the last LAT-detected event, while the vertical dashed black line marks the end of the GBM emission ($T_{GBM,95}$). Shaded gray areas mark intervals when the GRB is outside the FOV. The red vertical markers at the bottom of each panel indicate the arrival times of the events with probability >0.9 to be associated with the GRB, with the green cross being the event with maximum energy. The horizontal red line indicates the estimated duration of the GRB ($T_{LAT,100}$).</figcaption></figure>

while $T_{LAT,1}$ corresponds to the last event with $p > 0.9$. $T_{LAT,100}$ of the signal is simply $T_{LAT,1} - T_{LAT,0}$. These are also the quantities that define the “LAT” time interval, as previously discussed (see Table 1).

In order to correctly estimate the uncertainty on $T_{LAT,1}$ ($\delta T_{LAT,1}$) for an event with $n$ detected photons with probability $p > 0.9$, we define $\Delta t_{n-1,n}$ as the time interval between the second-to-last and the last event. Assuming Poisson statistics, the probability to measure an event between $t$ and $t+dt$ is $P(t, t+dt) = \lambda dt$, where $\lambda$ is the rate: in our case $\lambda = 2/\Delta t_{n-1,n}$. Therefore, we conservatively compute the uncertainty as $\delta T_{LAT,1} = 1/\lambda = \Delta t_{n-1,n}/2$. Similarly, considering the first two events with probability $p > 0.9$, we define the uncertainty on $T_{LAT,0}$ as $\delta T_{LAT,0} = \Delta t_{1,2}/2$. The error on $T_{LAT,100}$ follows using standard error propagation.

In order to better illustrate this analysis, Figure 3 shows two light curves of two bright GRBs, the long GRB 080916C (left panel) and the short GRB 090510 (right panel). For the long burst, the arrival time of the last event is substantially smaller than the end of the bin of the last detection. This could indicate that in the last bin the GRB is only marginally detected. For most of the other bursts, the arrival time of the last event is very close to the end of the last bin with positive detection, as in the case of GRB 090510.

### 3.6.3. Calculation of Energetics

In addition to reporting the flux and fluence of each GRB, for the subset of GRBs with measured redshift $z$ we also calculate their total radiated energy ($E_{iso}$). Starting from the measured spectrum of each burst, this is done by using the best-fit model over a specific energy range and by assuming that the energy emitted by a GRB at the source in the cosmological source frame is isotropically radiated. The isotropic radiated energy is

defined by following expression:

$$E_{iso} = \frac{4\pi d_L^2}{1+z} S(E_1, E_2, z), \quad (3)$$

where $d_L$ is the luminosity distance and $S(E_1, E_2, z)$ is the fluence integrated between the minimum energy $E_1$ and the maximum energy $E_2$. It can be expressed as

$$S(E_1, E_2, z) = T_{100} \int_{E_1/(1+z)}^{E_2/(1+z)} E N(E) dE. \quad (4)$$

Here $N(E)$ describes the best-fit spectral model, and $T_{100}$ represents the total duration of the burst as defined in the previous section. LAT data are always fit with a simple PL model in the energy range 100 MeV–10 GeV, i.e.,

$$N(E) = AE^{\Gamma}. \quad (5)$$

Finally, assuming a spatially flat universe $\Lambda$CDM model with $\Omega_{\Lambda} = 0.714$, $\Omega_M = 0.286$, and $H_0 = 69.6$ km s$^{-1}$ Mpc$^{-1}$ (Bennett et al. 2014; Planck Collaboration et al. 2016), the luminosity distance is given by (Weinberg 1972)

$$d_L(z, \Omega_{\Lambda}, \Omega_M) = (1+z) \frac{c}{H_0} \int_0^z \frac{dz'}{\sqrt{\Omega_M(1+z')^3 + \Omega_{\Lambda}}} \quad (6)$$

### 3.6.4. Localization

The LTF algorithm described in Section 3.2 returns the position of the GRB, as well as its detection probability. The steps of the procedure include refining the source location, and the position given in the final step is taken as the definitive one; no further optimization is performed.
<figure><img src="image_6.png" /><figcaption>**Figure 4.** Left panel: distribution of the TS for 14 short (red) and 155 long (blue) GRBs detected by the LTF algorithm. Right panel: TS values for long and short bursts as a function of the angle $\theta$ at the trigger time. Bursts that triggered an ARR are marked with a triangle.</figcaption></figure>

### 3.6.5. LLE Light Curve and Duration

The BB burst detection algorithm described in Section 3.4 provides a way of binning the data taking into account background fluctuations: blocks are defined only when an intrinsic rate variation above the background is detected, as opposed to an absolute variation. In our analysis we therefore define the onset of the LLE signal ($T_{\text{LLE,05}}$) as the starting time of the first block above background. Similarly, we define the $T_{\text{LLE,95}}$ as the ending time of the last block above background. The LLE duration ($T_{\text{LLE,90}}$) is simply defined as $T_{\text{LLE,95}}-T_{\text{LLE,05}}$.

## 4. Results

In the following subsections, we examine the main results of our analysis. The focus will be on the properties of the overall population, rather than a presentation of individual GRBs.

### 4.1. LAT Detections

This 10 yr catalog comprises 186 detections, 17 short GRBs (sGRBs), and 169 long GRBs (lGRBs). Adopting the analysis methods described in Section 3, we detect 169 GRBs with our likelihood analysis above 100 MeV. Of these, 155 are lGRBs and 14 are sGRBs. The distribution of the TS obtained by the LTF algorithm is shown in the left panel of Figure 4. The distribution peaks at relatively low values of TS (~30) and then smoothly falls with increasing TS value. Only a handful of GRBs (~5%) form a tail at very high TS (above 1000).

Using the LLE technique, 91 GRBs are found below 100 MeV. Out of those, 85 are lGRBs and 6 are sGRBs. Moreover, 17 of these GRBs (of which 2 are sGRBs) are found only with the LLE technique and are not detected at higher energies with the LAT standard analysis chain.

Of all 3044 triggers in our initial list, we thus detect ~6% at high energies with the LAT. About 18% of the LAT-detected bursts were outside the nominal FOV of the LAT at the time of the trigger. We note that the position in the sky of events that were outside the FOV at trigger time may have entered the FOV at a later time. Moreover, in 10 yr 220 triggers initiated an ARR of the satellite, a small fraction (<10%) of which are

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

caused by other sources, such as solar flares or particle events. Of these ARRs, 83 successfully resulted in LAT detection. The distribution of the LTF TS values as a function of $\theta$ at the trigger time is shown in the right panel of Figure 4. The highest TS values are seen for GRBs with $\theta < 50^{\circ}$.

Furthermore, this catalog includes four GRBs that triggered the LAT directly: one short burst, GRB 090510, and three long ones, GRB 131108A, GRB 160509A, and GRB 160821A. This underscores that onboard LAT GRB detections are relatively rare, implying exceptional brightness in high-energy gamma-rays. It is worth noting that the very bright GRB 130427A did not result in LAT onboard trigger, since the GBM had triggered and issued an ARR on the first emission peak (Preece et al. 2014), which was very bright at low energies but not particularly strong above 100 MeV.

Figure 5 shows the position in equatorial coordinates of 2357 GBM GRB triggers (gray symbols) detected over the 10 yr period of the catalog. A total of 160 lGRBs and 16 sGRBs also detected by LAT are marked by blue and red asterisks, respectively. In fact, out of the 169 likelihood-detected GRBs, 10 did not trigger the GBM instrument. Of these, two GRBs triggered Swift-BAT, namely, GRB 081203A and GRB 130907A, while six GRBs were reported by the IPN: GRB 090427A, GRB 110518A, GRB 120911B, GRB 140825A, the short GRB 160702A, and GRB 180526A. While GRB 120911B did not trigger GBM, it was the only burst to be later found in on-ground analysis of GBM data and announced by Gruber et al. (2012). Furthermore, we report for the first time the detection of GRB 100213C and GRB 111210B. These triggers were reported via private IPN communication, as stated in Section 2.2.

A list of all LAT detections is given in Table 2. For each event, we state the trigger date and time (both in UT and in MET), the final LTF localization with error, the off-axis and zenith angles at trigger time, whether an ARR was issued, the likelihood TS value and LLE significance, the redshift, and the references to the corresponding GCN circulars published by the LAT Collaboration. For those events detected only with the
<figure><img src="image_8.png" /><figcaption>**Figure 5.** Sky distribution of 2357 GBM-triggered GRBs (from 2008 July 14 to 2018 July 31) in equatorial coordinates (gray asterisks). Blue (red) asterisks indicate 160 (16) long (short) LAT-detected GRBs included in the 2FLGC over the same time period.</figcaption></figure>

LLE technique, we report the best possible localization of the burst as determined by, e.g., GBM or *Swift*.

In our sample, 34 GRBs have a measured redshift (19%), as compared to 10 (29%) in 1FLGC. For comparison, the fraction of *Swift*-detected bursts with redshift is ~29%.<fnref n="74" /> The smaller fraction of LAT bursts with a measured redshift in the 2FLGC with respect to the 1FLGC is not surprising, as ~50 new GRBs were discovered by our analysis, which have not been previously reported to the community. In addition, the improvements to the analysis techniques enable us to detect fainter GBRs, which are more difficult targets for follow-up observations.

On average, the (90% containment, statistical only) uncertainty in LAT detections is 0.36 with a range from 0.04 to 2.0. In order to assess the LAT location accuracy, we also checked for joint detections by Fermi-LAT and *Swift* and found that 75 bursts (~40%) have a BAT position, while 67 bursts (~36%) have an XRT position. By comparing LAT and *Swift/XRT* localizations of the codetected GRBs, we find that ~70% of the *Swift* localizations are inside the LAT 90% confidence region. The majority of the remaining XRT positions are only marginally outside the LAT region, indicating that the LAT localization error is slightly underestimated (~0.1).

### 4.1.1. Comparison with the First LAT GRB Catalog

The changes and improvements in the 2FLGC mean that the results reported here will differ from those in the 1FLGC. In the time interval of the 1FLGC, 2008 August to 2011 July (3 yr), we now recover more events: instead of 28 standard likelihood detections, we now have 50 detections. Three of these new detections are short GRBs, namely, GRB 081102B, GRB 090228A, and GRB 110728A. Four of the new detections come from non-GBM triggers.

The 1FLGC included 21 GRBs also detected with the LLE technique below 100 MeV. During the same period, we now find 25 LLE detections. Four of those—GRB 090531B, GRB 100225A, GRB 101123A, and GRB 110529A—are LLE-only bursts as reported also in the 1FLGC, with the first and the last one being short GRBs. The total number of LLE-only

detections is lower with respect to the 1FLGC, where we retrieved seven LLE-only bursts. Indeed, this is not surprising, since we now detect more events with the likelihood analysis thanks to Pass 8 and to the improved LTF pipeline.

As a result of the new analysis, we do not include in the current catalog two events that were included in the 1FLGC: GRB 091208B and GRB 110709A. Both GRBs were long, with estimated LAT durations of ~40 s; however, only three photons were detected for each GRB, and their detection was marked as marginal. The highest-energy photon in GRB 091208B was 1.2 GeV, while GRB 110709A had no detected emission above 500 MeV. By selecting Pass 8 data and applying the new detection algorithm, the significance of these two detections further decreased, thus resulting in their exclusion from the 2FLGC.

### 4.1.2. LAT Detections after 2011 July

We have also cross-checked the current catalog with the LAT detections that were publicly announced through GCN Circulars in the time period from 2011 July until 2018 August. Using the standardized catalog analysis described in Section 3.2, we now detect 31 previously unreported GRBs, for which no GCN has been issued. As expected, this is a much smaller relative increase than during the period of the 1FLGC, since Pass 8 data and the improved detection algorithm have been used since 2015.

Moreover, we do not retrieve eight GRBs that have previously been publicly announced by the Fermi-LAT Collaboration, namely, GRB 120916A (GCN 13777), GRB 130206A (GCN 14190), GRB 131018B (GCN 15357), GRB 140329A (GCN 16047), GRB 150127A (GCN 17356), GRB 150724B (GCN 18065), GRB 161202A (GCN 20229), and GRB 170810A (GCN 21452). In general, these are all GRBs that at the time of detection were reported with low significance, or with few photons. All these cases were analyzed at the time of the GCN writing either on ad hoc time intervals chosen by the burst advocates or on the 10 real-time LTF temporal windows. These differ from the five fixed time intervals chosen for the catalog analysis presented in Section 3.2, thus leading to different results. GRB 130206A and GRB 150127A were previously reported through GCNs as marginal LLE detections, both with a significance <3σ, again not matching the current catalog requirements.

<footnote n="74">⁷⁴ https://swift.gsfc.nasa.gov/archive/grb_table/</footnote>
<figure><img src="image_9.png" /><figcaption>**Figure 6.** Temporal extended emission for two GRBs detected by LAT at very late times, namely, GRB 080818B (left panel) and GRB 160503A (right panel). Markers and colors are the same as in Figure 3.</figcaption></figure>

<figure><img src="image_10.png" /><figcaption>**Figure 7.** (a) $T_{\text{LAT,0}}$ and (b) $T_{\text{LAT,1}}$ calculated in the 100 MeV–100 GeV energy range vs. the same quantities calculated in the 50–300 keV energy range. Panel (c) shows the onset times ($T_{\text{LAT,0}}$) in the 100 MeV–100 GeV vs. the durations ($T_{\text{GBM,90}}$) in the 50–300 keV energy range. The solid line denotes where values are equal. Blue and red circles represent long and short GRBs, respectively. In panel (c), we additionally mark the GRBs that were outside the FOV at trigger time with a thick orange contour.</figcaption></figure>

## 4.2. LAT Onset Times and Duration

In the following paragraphs we discuss the temporal properties of the bursts in our sample. As presented in Section 2.5, the classification of GRBs into long and short classes is derived from the low-energy duration as measured by GBM in the 50–300 keV energy band. The LAT durations are calculated in the 100 MeV–10 GeV energy range. Table 3 summarizes the various temporal characteristics of the GRBs in our catalog. This includes the values of $T_{\text{GBM,05}}$, $T_{\text{GBM,95}}$, and $T_{\text{GBM,90}}$ for GBM; $T_{\text{LLE,05}}$, $T_{\text{LLE,95}}$, and $T_{\text{LLE,90}}$ for LLE; and $T_{\text{LAT,0}}$, $T_{\text{LAT,1}}$, and $T_{\text{LAT,100}}$ for the LAT. Two GRBs, GRB 100213C and GRB 111210B, were reported only by the IPN through private communication: we do not provide any duration information for those. We mark all non-GBM durations in Table 3 with an asterisk in the $T_{\text{GBM,05}}$ column.

For some GRBs, the LAT detection of the first photon occurs at very late times. This could be due to high-energy photons not

being emitted during the initial phase, but it could also be due to observational constraints, where the GRB location is outside the FOV for long intervals. This is illustrated for two GRBs in Figure 6. In both panels, blue points are photon flux measurements, while upper bounds are displayed as black arrows. In the left panel, the first shaded gray area marking when GRB 080818B was outside the FOV spans almost 10 ks ($T_{\text{LAT,0}} = 9.0 \pm 0.6$ ks). The estimated duration of the burst, $T_{\text{LAT,100}} = 500 \pm 200$ s, is almost not visible owing to the late time of the detection. Similarly, in the right panel, the first detection of GRB 160503A occurs at 5.3 ks, again after a period of several kiloseconds where the burst was first not detected and then outside the FOV. In this case, the duration was $T_{\text{LAT,100}} \sim 18$ ks.

In Figure 7(a) we compare the onset times estimated in the LAT energy band (100 MeV–100 GeV) with the ones estimated in the GBM energy band (50–300 keV). A negative $T_{\text{GBM,05}}$ value in the low-energy band means that the burst
<figure><img src="image_11.png" /><figcaption>**Figure 8.** Composite light curve of GRB 160625B: summed GBM/NaI detectors (first two panels), GBM/BGO (third panel), LLE (fourth panel), and LAT rates above 100 MeV (bottom panel). The gray shaded area indicates the $T_{GBM,90}$ calculated in the 50–300 keV energy range. The dashed red lines mark the GBM trigger time. A zoom-in around trigger time is shown in two insets for each energy range of the NaI detectors, highlighting the peak detected at low energies that caused the trigger in GBM. This peak does not show any counterpart at higher energies. The association of the low-energy ($<250$ keV) late-time emission (from 500 to 800 s post-trigger) with the main GRB emission episode has been cross-checked by the GBM Team (FGGC).</figcaption></figure>

onset occurred before the trigger time. In general, we notice that the high-energy emission starts significantly later with respect to the low-energy one, for both long and short bursts. Burst durations are compared in panel (b) of the same figure. Here, the end of the signal at high energies ($T_{LAT,1}$) appears to be significantly later than the one measured in the GBM energy band. Both these characteristics were already reported in the 1FLGC. Our results confirm and strongly support the claim that when high-energy emission is observed in GRBs, this emission is delayed and lasts longer compared to that in the low-energy band.

In Figure 7(c) we show the onset time ($T_{LAT,0}$) of the high-energy emission versus the burst duration ($T_{GBM,90}$) in the 50–300 keV range. It is worth noting that the $T_{LAT,0}$ of the majority of GRBs (both long and short ones) occurs before the prompt emission measured by the GBM is over. Events that were outside the nominal LAT FOV (~75°) at the time of the GBM trigger are marked with thick orange contours. They compose the majority of GRBs where the onset of the high-energy emission came after the low-energy emission had faded, indicating that most such events are due to observational bias. This effect will be further investigated below.

As shown in panels (a) and (b) of Figure 7, there are just a few outliers that have high-energy emission that is not delayed and/or has shorter duration compared to the

low-energy band. However, since the procedure to calculate onset times and durations differs between the two energy ranges, we caution that further analysis is needed before strong conclusions can be drawn about individual GRBs. The difference is in most cases less than a few seconds. The most prominent outlier to the right of the line is GRB 160625B, where the GRB $T_{GBM,05}$ is ~190 s, whereas $T_{LAT,0}$ is ~25 s. However, this burst showed three emission episodes spread over a period of more than ~10 minutes, as shown in Figure 8. The first one triggered the GBM, a second one 3 minutes later resulted in LAT onboard trigger, and then the GBM triggered again 10 minutes after the first trigger. It is thus not surprising that the $T_{GBM,05}$ is much greater than the arrival time of the first LAT photon.

Short GRBs in general have more similar onset times in LAT and GBM. They also exhibit shorter durations in the high-energy range, although they last still significantly (generally more than an order of magnitude) longer than at lower energies. The short GRB 170127C is the short burst with the longest-lasting high-energy duration, more than 2 ks.

In our sample, 16 GRBs have high-energy emission lasting over 5 ks, and four have durations over 10 ks, namely, GRB 160623A (~35 ks), GRB 130427A (~34 ks), GRB 140810A (~18 ks), and GRB 160503A (~18 ks). Figure 9 shows the temporal extended emission for the two longest bursts, GRB 130427A in the left
<figure><img src="image_12.png" /><figcaption>**Figure 9.** Temporal extended emission for the two longest LAT GRBs, GRB 130427A (left panel) and GRB 160623A (right panel). For the bright GRB 130427A we display the results of the fit for a simple PL model with temporal index $\alpha$ (red line) and with a BPL with indices $\alpha_1$ and $\alpha_2$ (green line). For GRB 160623A, only the result of the fit with a simple PL is shown.</figcaption></figure>

<figure><img src="image_13.png" /><figcaption>**Figure 10.** Left panel: difference of the $T_{\text{LAT,0}}$ and the $T_{\text{GBM,05}}$ with respect to the incident angle of the GRB at the time of the trigger. The yellow region highlights the GRBs that were outside the LAT FOV ($\theta \sim 75^\circ$) at the time of the GBM trigger. Right panel: zenith angle ($\zeta$) vs. incident angle ($\theta$). The symbols are colored as a function of the difference of the $T_{\text{LAT,0}}$ and the $T_{\text{GBM,05}}$ (left panel, y-axis). Circles mark long GRBs, while stars mark short GRBs. Bursts that triggered an ARR are marked with red contours. The vertical dashed line marks the LAT FOV, while the horizontal line marks the zenith angle cut used in the analysis ($\zeta_{\text{MAX}} = 105^\circ$). In general, LAT-detected events at $\zeta > 105^\circ$ and $\theta > 75^\circ$ are seen at very late times or thanks to ARRs.</figcaption></figure>

panel and GRB 160623A in the right panel. In each panel, we also indicate the fit results to the temporal decay, giving the corresponding model parameters in the top right corner. This will be further discussed in Section 4.7.

As already mentioned, a possible bias in the estimation of the onset time in the LAT is related to the initial position of the GRB at the time of the GBM trigger. For a GRB outside the nominal LAT FOV at trigger time, the first significant detection would happen only when the GRB reenters the FOV. We further illustrate this effect in Figure 10. In the left panel the delay of the LAT onset time with respect to the GBM one is plotted as a

function of the incident angle of the GRB, while in the right panel we plot the zenith angle as a function of the incident angle. It is evident that all GRBs that were outside the LAT FOV at trigger time have a large delay ($\gtrsim 100$ s) with respect to the GBM trigger, which corresponds to the time needed for the GRB to reenter the LAT FOV. On the other hand, we also measure significant delays for GRBs that were in the FOV at the time of the GBM trigger, supporting the intrinsic nature of the delay of the high-energy component.

In the right panel of Figure 10, we highlight GRBs that resulted in an ARR, marking each symbol with a red contour.
<figure><img src="image_14.png" /><figcaption>**Figure 11.** Left panel: off-axis angle for GRB 170127C as a function of time since the GBM trigger. The horizontal red line corresponds to $\theta = 75^\circ$. Right panel: temporal extended emission for GRB 170127C. The blue circle is the photon flux measurement with a significance of $\text{TS} > 10$, while, for lower values of $\text{TS}$, upper bounds are displayed as black arrows. The horizontal red line indicates the estimated duration of the burst, and the thick parts of it correspond to the uncertainty on the $T_{\text{LAT},0}$ and $T_{\text{LAT},100}$ parameters. The red vertical markers indicate the arrival time of each photon with probability >90% to be associated with the GRB, with the green cross being the event with maximum energy.</figcaption></figure>

<figure><img src="image_15.png" /><figcaption>**Figure 12.** Onset times ($T_{\text{LLE},05}$; left panel) and durations ($T_{\text{LLE},90}$; right panel) calculated using LLE data in the 30 MeV–1 GeV energy range vs. the corresponding quantities ($T_{\text{GBM},05}$ and $T_{\text{GBM},90}$) calculated using GBM data in the 50–300 keV energy range. The solid green line denotes where values are equal. Blue and red circles represent long and short GRBs, respectively.</figcaption></figure>

Most of the bursts for whom an ARR was issued were in the LAT FOV at the time of the trigger, whereas in seven cases the GRBs were outside the LAT FOV and the detection happened only at later times. To better illustrate this effect, we display the case of the short GRB 170127C in Figure 11. This burst was at a zenith angle of $\zeta \sim 94^\circ$ and at an off-axis angle $\theta \sim 142^\circ$ when it triggered the GBM (it is the outlier sGRB seen to the far right in both panels of Figure 10). The trigger resulted in an ARR, and the spacecraft slewed to move the location of the burst close to the center of the FOV (at $\sim 17^\circ$). This can be seen in the left panel of Figure 11, where the blue circles show the evolution of $\theta$ as a function of time after the trigger. In the right panel, we show the photon flux light curve resulting from the time-dependent analysis. The first detection is at $\sim 400$ s, well beyond the end of the GBM signal ($T_{\text{GBM},95} = 0.13$ s).

### 4.3. LLE Onset and Duration

If we restrict our considerations to the LLE analysis, where the bulk of the emission is in the energy range from 30 MeV to 1 GeV, we see that the left panel of Figure 12 shows how the onset times are relatively similar to the onset times as measured by the GBM. Here, two GRBs are not shown: GRB 120624B and GRB 150513A. Both GRBs triggered *Swift* before they triggered GBM, 257 s (Barthelmy et al. 2012) and 157 s (Kocevski et al. 2015) before the GBM trigger time, respectively. As a result, since all our calculations are referred to GBM trigger times, $T_{\text{LLE},05}$ is negative and omitted from the figure (see Table 3).

In contrast to the emission above 100 MeV, the right panel of Figure 12 indicates that the duration of the signal in LLE is systematically shorter than the duration of the signal in the GBM, as was seen also in the 1FLGC. If we assume that the
<figure><img src="image_16.png" /><figcaption>**Figure 13.** Left panel: distribution of energy fluence calculated in the 10–1000 keV energy range for 178 bursts detected by the LAT compared to the entire sample of 2357 GRBs detected by GBM over the same time period. Middle panel: distribution of peak photon flux in the 10–1000 keV energy range for the same sample of LAT and GBM-detected populations. Right panel: GBM peak photon flux (10–1000 keV), as measured on a 1024 ms timescale, vs. the energy fluence (10–1000 keV) derived using a spectral model fit to a single spectrum over the entire duration of the burst.</figcaption></figure>

LLE emission is dominated by the same emission episodes as in the GBM, we can infer that the pulses that make up the time profile of the prompt emission are systematically shorter in the LLE range than at lower energies. This behavior has previously been reported by Norris et al. (1996) and Norris (2002) using BATSE data, as well as for several LAT-observed GRBs (e.g., Axelsson et al. 2012; Bissaldi et al. 2017; Vianello et al. 2018).

### 4.4. Comparison to the GBM Population

Since the majority of our triggers come from the GBM, and the GBM has observed nearly all GRBs in our sample, we examine how the LAT-detected bursts are drawn from the general GBM population covering the same 10 yr time period. For this comparison, we extracted the peak photon flux, as measured on a 1024 ms timescale, and energy fluence measured by the GBM in the 10–1000 keV energy range from the FGGC. Here the GBM fluence is derived from the parameters of the best-fit spectral model applied to GBM data over a time interval where the signal-to-noise ratio (S/N) exceeds a predefined value (S/N > 3.5; see Gruber et al. 2014, for more details). This requirement ensures that there are enough counts to perform a spectral fit, but as a result, the time interval does not always coincide with $T_{GBM,90}$. Note that eight GRBs, two triggered by Swift and six by the IPN, were not detected by the GBM and are omitted from this comparison and from the following figures.

Figure 13 shows the distribution of the energy fluence (left panel) and of the peak photon flux (middle panel) for 178 bursts detected by the LAT compared to the entire sample of 2357 GRBs detected by GBM over the same time period. Here we have also made a distinction between short and long bursts for both the LAT (16 sGRBs and 162 IGRBs) and GBM (400 sGRBs and 1957 IGRBs) populations, showing a bifurcation in the range of flux and fluence values covered by these two classes of bursts. The right panel shows the peak photon flux plotted against the energy fluence for the LAT bursts compared again to the entire GBM burst catalog. Again, we separate short and long bursts for both the LAT and GBM populations.

These comparisons show that although the majority of the LAT-detected GRBs come from the GBM-detected bursts with the highest peak flux and fluence, they cover a large range. LAT-detected short (long) bursts are present with a fluence

$>10^{-7}$ erg cm$^{-2}$ ($>8 \times 10^{-7}$ erg cm$^{-2}$) and with a peak flux $>0.8$ photons cm$^{-2}$ s$^{-1}$ ($>1.5$ photons cm$^{-2}$ s$^{-1}$). The LAT-detected long GRBs cover more than two orders of magnitude in both distributions, and the prominence of bright GRBs is even less pronounced in the short GRB sample. The spread is also evident from the right panel of Figure 13, where the cluster of LAT events is only slightly shifted with respect to the GBM one. The burst with the highest fluence (and flux) is GRB 130427A. It is worth noting that Figure 13 does not include any selection on the $\theta$ angle.

### 4.5. Flux, Fluences, and Photon Indices from the Time-integrated Analysis

The results of the likelihood analysis are summarized in Table 4. For each time window, we report the number of detected and predicted LAT events in the ROI, the resulting TS, the spectral index obtained using a PL fit, and the LAT flux and fluence calculated in the 100 MeV–100 GeV energy range. For 34 GRBs with known redshift we also report the total radiated energy ($E_{iso}$).

Figure 14 shows the distributions of fluxes (left panels) and fluences (right panels) as a function of the measured duration of the signal in the “GBM” (top row) and “LAT” (bottom row) time windows. LAT fluxes decrease with increasing burst duration in both time windows, as expected. In the “GBM” time window, the LAT fluence seems to be clustered around a value of $10^{-6}$ erg cm$^{-2}$ for the majority of IGRBs (regardless of duration), while sGRBs show slightly lower values. Both groups have bursts that are very much brighter than the average. At late times, there is instead a tendency for the fluence to increase with duration. The same conclusion can be drawn from the fluence values in the “LAT” time window, where most of the values are distributed around $\sim 5 \times 10^{-6}$ erg cm$^{-2}$ and there is a less evident spread toward higher values.

Comparing our results to Figure 11 in the 1FLGC, we find that the four “hyperfluent” GRBs are no longer outliers. Instead, they are part of a continuous distribution. The range in both flux and fluence has also increased dramatically as compared to the sample in the 1FLGC.

In Figure 15 we then compare the LAT fluence calculated in the 100 MeV–100 GeV energy range during the “GBM” time window with the GBM fluence calculated between 10 keV and
<figure><img src="image_17.png" /><figcaption>**Figure 14.** Fluxes (left panels) and fluences (right panels) calculated in the 100 MeV–100 GeV energy range vs. GRB durations. Top panels show fluxes and fluences evaluated in the “GBM” time window vs. durations calculated in the 50–300 keV energy range ($T_{\text{GBM,90}}$). Bottom panels show fluxes and fluences evaluated in the “LAT” time window vs. durations calculated in the 100 MeV–100 GeV energy range ($T_{\text{LAT,100}}$). Blue and red circles represent long and short GRBs, respectively.</figcaption></figure>

1 MeV (left panel) and with the LAT fluence calculated in the same energy range during the “EXT” time window (right panel). In the left panel, it can be noted that the “GBM” time window is dominated by the low-energy emission, with the 100 MeV–100 GeV energy range contributing only a small fraction of the emission for the majority of long GRBs. Indeed, most events are clustered to the left of the solid and dashed lines, which indicate equality and a factor of 10 less, respectively. For short GRBs this difference seems less pronounced, and several lie close to the solid line of equality. Comparing the “GBM” and “EXT” time windows in the right

panel, the points are instead much closer to the line of equality, suggesting that the high-energy emission in the two time windows is comparable. As in Figure 14, the four “hyperfluent” GRBs of the 1FLGC (GRB 080916C, GRB 090510, GRB 090902B, and GRB 090926A) are no longer outliers.

In Figure 16 we also compare the photon index measured by the LAT during the “GBM” time window ($\Gamma_{\text{GBM}}$) and the “EXT” time window ($\Gamma_{\text{EXT}}$). Both indices are plotted as a function of the GRB duration as calculated in the 50–300 keV energy range (top panels) and of the LAT flux calculated in the “GBM” (bottom left panel) and “EXT” (bottom right panel)
<figure><img src="image_18.png" /><figcaption>Figure 15. Left panel: fluences calculated in the 10–1000 keV energy range vs. fluences calculated in the 100 MeV–100 GeV energy range. All values are estimated in the “GBM” time window. Right panel: fluences calculated in the 100 MeV–100 GeV energy range evaluated in the “EXT” time window vs. the same quantities evaluated in the “GBM” time window. The solid green lines denote where values are equal. The dashed and dot-dashed green lines are shifted by factors of 10 and 100, respectively.</figcaption></figure>

time windows. The photon index shows no sign of being correlated with either the GBM duration or the flux in either time window, and it is similar for long and short GRBs. The value is indeed similar between the two time windows but is slightly harder in the “EXT” window. In the “GBM” time window, the values of the photon index are more scattered, with a mean value of $\Gamma_{GBM} = -2.49$ and a 10th (90th) percentile of $-3.22$ ($-1.86$). In the “EXT” time window, the values are more uniform, with a mean of $\Gamma_{EXT} = -2.03$ and a 10th (90th) percentile of $-2.45$ ($-1.6$). For comparison we recall the same values reported in the 1FLGC: $-2.08 \pm 0.04$ in the “GBM” time window and $-2.00 \pm 0.04$ in the “EXT” time window. While the latter is in agreement with the current value, the photon index during the “GBM” time window was much harder than the one we derive in the 2FLGC. Interestingly, it showed a weak inverse correlation with the duration of the burst (see Figure 26 of the 1FLGC). This correlation is now less evident in the larger sample of bursts, but the $\Gamma_{GBM}$ distribution still underlines the agreement with previous findings that the spectra of short-duration GRBs tend to be harder. No clear trend can be seen in the comparison with flux (bottom panels), except a slight tendency for low-flux GRBs to show harder spectra when looking in the “EXT” window.

## 4.6. Energetics

In order to more closely study the energetics of the bursts in this catalog, and to put the detections in a wider context, we focus on the GRBs with known redshift. We decided to compare our sample to other bursts with measured redshift detected by Swift and GBM. As of the end of 2018 July, Swift-BAT has detected 1246 GRBs, of which ~35% have a measured redshift. In the case of the GBM-detected GRBs, only ~5% have measured redshift. The redshift distributions of 405 bursts detected by Swift-BAT (gray histogram), 116 bursts detected by GBM (cyan histogram), and 34 bursts detected by

LAT (blue histogram) are shown together in Figure 17. We see no obvious difference between the three distributions.

We next compare the isotropic radiated energy ($E_{iso}$) and the bolometric gamma-ray peak luminosity ($L_{iso}$) of LAT-detected GRBs to the same quantities in the Swift and GBM samples. The values for $E_{iso}$ are computed according to Equation (3) in the 1 keV–10 MeV energy range. In the case of GBM-detected GRBs, we adopt the fluence listed in the FGGC as computed from the best-fit spectral model, which is usually calculated on a slightly different time interval with respect to the burst $T_{GBM,90}$, according to the burst brightness.

In order to compute $E_{iso}$ of Swift-detected events (with no GBM observation), we used the parameters of the best-fit spectral models obtained in the 15–350 keV energy range reported in the Swift-BAT online catalog<fnref n="75" /> (see Lien et al. 2016, for more details). These are calculated over a time interval corresponding to a duration that contains 100% of the burst emission. For both the GBM and BAT $E_{iso}$ calculation, we only consider bursts for which the spectral parameters are globally well constrained (see Gruber et al. 2014). Thus, we find 116 (405) GBM (BAT) GRBs that satisfy these criteria, out of which 108 (376) are IGRBs and 8 (29) are sGRBs. The LAT sample comprises 32 IGRBs (2 of the 34 were not detected by the GBM, as previously discussed) and only one sGRB (090510).

We also calculate the isotropic luminosity $L_{iso}$, which takes into account the GRB prompt emission spectrum and is defined in a 1 s time interval centered around the time of the peak flux. It can be expressed as

$$L_{iso} = 4\pi d_L^2 P(E_1, E_2, z), \quad (7)$$

<footnote n="75">⁷⁵ https://swift.gsfc.nasa.gov/results/batgrbcat/</footnote>
<figure><img src="image_19.png" /><figcaption>Figure 16. Photon indices Γ<sub>GBM</sub> (left panels) and Γ<sub>EXT</sub> (right panels) calculated in the 100 MeV–100 GeV energy range vs. duration calculated in the 50–300 keV energy range (top panels) and LAT 100 MeV–100 GeV fluxes (bottom panels). These are calculated in the “GBM” (bottom left) and “EXT” (bottom right) time windows, respectively. The green dashed lines denote the mean values and the dotted lines the 10th and 90th percentiles of each distribution. Blue and red circles represent long and short GRBs, respectively.</figcaption></figure>

where P(E<sub>1</sub>, E<sub>2</sub>, z) represents the bolometric peak flux, defined as

$$P(E_1, E_2, z) = \int_{E_1/(1+z)}^{E_2/(1+z)} E N(E) dE. \qquad (8)$$

As with E<sub>iso</sub>, L<sub>iso</sub> is computed in the 1 keV–10 MeV energy range, for GBM-detected GRBs, using the 1 s peak flux of the best-fit model as reported in the FGGC. We again consider only GRBs whose time-integrated spectra are well defined, as reported in the GBM and Swift-BAT GRB catalogs. This leaves us with 394 BAT GRBs and with the same number (116) of

GBM GRBs. The slightly lower number of BAT GRBs is expected, as the time interval (and thereby the number of photon counts) is smaller.

Figure 18 shows the distribution of E<sub>iso</sub> (left panel) and L<sub>iso</sub> (right panel) as a function of redshift. Swift-BAT and GBM bursts are indicated by gray and black points, respectively, with long (short) bursts marked with filled (open) symbols. LAT long and short bursts are marked with the standard blue and red circles, respectively, used in this paper. LAT-detected GRBs populate the top portion of both distributions, as was previously seen in the 1FLGC. At that time, this figure only contained nine LAT-detected GRBs with redshift. It is
<figure><img src="image_20.png" /><figcaption>**Figure 17.** Redshift distribution of 34 GRBs detected by LAT (blue histogram), 405 GRBs detected by Swift-BAT (gray histogram), and 116 GRBs detected by GBM (cyan histogram).</figcaption></figure>

worth noting that quite a few bursts have a moderate 1 keV–10 MeV $E_{iso}$ ($\lesssim 10^{-53}$ erg) yet have nevertheless been detected by the LAT.

### 4.7. Time-resolved Light Curves

We now turn to the temporal decay of the high-energy extended emission. Using the analysis described in Section 3.6.2, we were able to determine the evolution of the flux as a function of time for 115 long and 11 short GRBs in our sample. This is shown in the two panels of Figure 19, displaying the temporal decay of long (left panel) and short (right panel) bursts separately. Each event is marked with a different color. The light curves of both sGRBs and lGRBs show a fairly large spread in the observer frame.

In order to determine the corresponding temporal decay index, we perform a fit of all the light curves maximizing the $\chi^2$, with two different spectral models, namely, (1) a simple PL,

$$F(t) = F_0 \left(\frac{t}{T_0}\right)^{-\alpha}, \qquad (9)$$

where $\alpha$ is the temporal decay index, $T_0$ is the GRB trigger time, and $F_0$ is the normalization flux; and (2) a BPL,

$$F(t) \propto t^{-p} \begin{cases} p = \alpha_1 & \text{for } t < T_b \\ p = \alpha_2 & \text{for } t \ge T_b, \end{cases} \qquad (10)$$

with index $\alpha_1$ for times before the break time $T_b$ and index $\alpha_2$ afterward. If there are at least three flux points (with $TS > 10$) in the light curve after the $T_{GBM,95}$, we fit a PL, and if there are at least four flux points, we also try a BPL.

The results of the fits are presented in Table 5. By fitting the flux temporal decay with a BPL, we find a significant improvement in 12 cases. We show three examples in Figure 20 (GRB 090510, GRB 150627A, and GRB 180720B; a fourth, GRB 130427A, has already been shown in the left panel of Figure 9). The BPL fit is indicated with a solid green line, and the corresponding fit values are given in the upper right corner of each panel. The PL fit is shown with a red solid line for comparison. In all but two cases the light curves manifest a steep-to-shallow decay, while for

GRB 171120A and GRB 180720B (right panel in Figure 20) the decay steepens after the break.

The distribution of late-time temporal decay indices ($\alpha$ or $\alpha_2$) is displayed in Figure 21, together with a Gaussian fit (black dashed line) to the distribution. The distribution comprises 88 GRBs, including 86 long and 2 short bursts. Among the long (short) GRBs, 77 (1) are best fit with the PL model, while 11 (1) prefer a BPL model. This is a large increase compared to the 1FLGC, where only nine GRBs had enough data to allow the decay index to be determined, ranging from 0.8 (for GRB 090902B) to 1.8 (for GRB 080916C) and with a mean value of $\sim 1.1$. We now find a mean value of $0.99 \pm 0.04$ with a standard deviation of $0.80 \pm 0.07$, still in agreement with the results presented in the 1FLGC.

These results and their interpretation will be discussed in more detail in Sections 5.6 and 6.3, but it is worth noting that there are several cases in which a BPL would likely be required if the data in the “GBM” time window had been included. Furthermore, several GRBs show features in the light curve that deviate from both a PL and BPL. Four light curves exemplifying both these cases are displayed in Figure 22.

In the figure, GRB 131108A (top left) exemplifies how some bursts display strong variability in the LAT $> 100$ MeV light curve during the GBM emission. GRB 140206B (top right) instead shows the possible presence of a late-time high-energy pulse. GRB 160625B (bottom left) shows a strong pulse with a very sharp decay. GRB 160816A (bottom right panel) would clearly require a break to accommodate the data in the “GBM” time window. However, in the “EXT” time window used for the fits, a PL model is statistically preferred. These peculiar features are all observed for the first time in the 2FLGC. Some light curves, like GRB 140206B (top right panel in Figure 22), show a sharp step between the end of the “GBM” and the beginning of the “EXT” time window (marked by the vertical black dashed line). In these cases, a fit of the complete data set would again favor a BPL model rather than the PL one currently used in the “EXT” time window. In a few even more extreme cases, the light curve could not be fit at all in the “EXT” time window, since there is only a single point in the light curve after the “GBM” time window.

In Figure 23, we show the 100 MeV–100 GeV luminosity evolution for the 34 GRBs in our sample with measured redshift. Among those, there is only one short burst, namely, GRB 090510. The three panels of the figure show first the light curves in the observer frame (left), the evolution as a function of time in the source frame (middle), and finally the luminosity divided by isotropic energy ($E_{iso}$; right), calculated in the 1 keV–10 MeV energy range. Each correction brings the light curves closer together, and in the right panel there is a remarkable alignment of all the GRBs. This analysis was done following the one presented by Nava et al. (2014), where a similar result was found. It is worth noting that in the right panel of Figure 23 one of the GRBs, GRB 160623A, does not align with the others, indicating a possible outlier. This burst was occulted by Earth for a large part of its duration, and the GBM trigger occurred $\sim 50$ s after the start of the GRB based on the Konus-Wind light curve. This likely leads us to underestimate the total energy release $E_{iso}$ of the burst and thus to overestimate the normalization of the light curve. In the right panel we also include a linear fit to all 34 GRBs, indicated by the solid line. The decay index is $1.25 \pm 0.03$. For
<figure><img src="image_21.png" /><figcaption>**Figure 18.** Left panel: isotropic radiated energy vs. redshift. Right panel: isotropic luminosity vs. redshift. Blue/red circles indicate the LAT long/short GRBs, and black and gray circles indicate GBM/Swift-BAT GRBs, with short bursts always marked by open symbols.</figcaption></figure>

<figure><img src="image_22.png" /><figcaption>**Figure 19.** Flux calculated in the 100 MeV–100 GeV energy range vs. elapsed time after trigger for 115 long (left panel) and 11 short (right panel) bursts. Each color represents a separate GRB.</figcaption></figure>

comparison, we also show a dashed line with decay index 10/7 (see further Section 5.6).

4.8. *High-energy Events*

The highest-energy GRB photon ever recorded by *Fermi* thus far is a 94.1 GeV event connected with GRB 130427A (Ackermann et al. 2014). While displaying photon energies of a few hundred MeV is a common feature among the LAT-detected GRBs, higher energies are relatively rare. Table 6 summarizes the highest-energy photon characteristics for each burst in our sample. It lists the total number of photons detected with probability >90% of belonging to the burst, as well as the energy, arrival time, and probability of the highest-energy

photon detected in the “GBM” time window. We also list the same quantities calculated in the time-resolved analysis.

Figure 24 shows the fraction of GRBs detected above selected energy thresholds (250 MeV, 500 MeV, 1 GeV, 5 GeV, 10 GeV, 50 GeV). A sharp drop from ~70% to ~30% is seen at 5 GeV. There are three GRBs with emission above 50 GeV (2%), namely, GRB 130427A (95 GeV), GRB 140928A (52 GeV), and GRB 160509A (52 GeV).

Our sample of 34 GRBs with measured redshift also allows us to study the source-frame-corrected energies. This is shown as the dashed line in Figure 24. This distribution shows a more gradual decrease with energy: almost 80% of the included GRBs have a maximum source-frame photon energy above 5 GeV, and ~12% (four GRBs) above 100 GeV. The highest
<figure><img src="image_23.png" /><figcaption>**Figure 20.** Temporally extended emission for three GRBs that are best fit by a BPL model (a fourth one, GRB 130427A, is shown in the left panel of Figure 9). GRB names are shown in the title in each panel. The red vertical markers indicate the arrival times of the events with probability >0.9 to be associated with the GRB, with the green cross being the event with maximum energy. The green vertical dashed lines indicate the first and the last event. The vertical dashed black line marks the $T_{GBM,95}$. The PL and BPL fits are indicated with solid red and green lines, respectively. Fit parameters are given in the top right corner of each panel. The asterisk highlights the value selected as the temporal decay index in the analysis (which, in these cases, always corresponds to $\alpha_2$ of the BPL model).</figcaption></figure>

<figure><img src="image_24.png" /><figcaption>**Figure 21.** Distribution of the temporal decay indices for 86 long (blue histogram) and 2 short (red histogram) bursts, calculated in the “ETX” time window. A Gaussian fit of the IGRB indices (mean $0.99 \pm 0.04$; standard deviation $0.80 \pm 0.07$) is superimposed on the distributions (dashed black line).</figcaption></figure>

source-frame energy is a 147 GeV photon from GRB 080916C, at $z = 4.35$ (Atwood et al. 2013b). The figure also includes a linear fit to the bin centers of the source-frame distribution. The fit is remarkably good, showing that the fraction of GRBs decreases as $A \times \log(E/1 \text{ MeV}) + B$, where $A = -49 \pm 4$ and $B = 266 \pm 21$.

Figure 25 shows the energy of the highest-energy photon in each GRB as a function of arrival time (left panel). In the right panel of the figure, the arrival time is shown as a fraction of $T_{GBM,90}$, calculated in the 50–300 keV range. No clear pattern can be distinguished, with long and short bursts overlapping in the right panel. The one clear outlier in the right panel is the short GRB 170127C, where the highest-energy photon (500 MeV) was detected almost 3 ks after the trigger. This GRB was outside the LAT FOV at $T_0$, and the data therefore

only cover the time from around 300 s after the trigger (see also Section 4.2 and Figures 10 and 11).

In order to check for more high-energy photons to be detected beyond the standard energy range (>100 GeV), we performed an additional analysis up to 300 GeV. No LAT detection was made above 100 GeV. Table 7 presents the 29 GRBs in the sample detected above >10 GeV. All high-energy photons with a probability >90% and with an observed energy >10 GeV are reported for each burst. As in Table 6, we specify the photon energy (listed in decreasing order), its arrival time, the GRB redshift, and the source-frame-corrected energy ($E_{sf}$). Values of $E_{sf} > 100$ GeV are marked in bold. The only short burst listed in this table is GRB 090510: in this case, a ~30 GeV photon was detected 830 ms after the GBM trigger time. GRB 130427A holds the record with 17 photons detected above 10 GeV, with the highest event ever detected from a GRB (the 94.1 GeV photon) observed 243 s post-trigger. The second burst with the most HE photons is GRB 090902B, with seven photons detected above 10 GeV. There are two bursts where a high-energy photon is detected at very late times (>10 ks): GRB 130427A (~34 ks) and GRB 160623A (~12 ks).

## 5. Discussion

We will now discuss our results and compare them to previous results, in particular what was seen in the 1FLGC. The discussion generally follows the outline of Section 4, starting with the LAT sample as a whole, then continuing with a comparison to the GBM sample. Finally, we will consider the energetics, temporal decay, and possibilities for detections at very high energy (VHE). The aim of this section is to put our results in a wider context; broader implications in the framework of theoretical models will be discussed in Section 6.

### 5.1. Detectability of LAT Bursts and LAT Detection Rate

Before launch it was estimated that the LAT would detect 10–12 GRBs per year above 100 MeV (Band et al. 2009). The results from the first GRB catalog showed 28 GRBs detected in the first 3 yr of the mission, slightly below expectations. The current work instead shows that the LAT has exceeded
<figure><img src="image_25.png" /><figcaption>**Figure 22.** Temporal extended emission for four GRBs showing a peculiar behavior (the name of each GRB is shown in the title of each panel). Markers are the same as in Figure 20.</figcaption></figure>

expectations, with 169 GRBs detected above 100 MeV in 10 yr. This is in large part due to the continuous improvements in event analysis and detection algorithms.

It is interesting also to look specifically at the highest energies of the LAT-detected GRBs. As can be seen from Figure 24, only a small fraction of GRBs are seen in the upper energy range of the LAT. For example, 20% of the GRBs have detected emission above 10 GeV, corresponding to ~3 GRBs per year. We will further discuss the occurrence of these highest-energy events in Section 5.5.2.

Another interesting aspect to discuss is that while the highest-energy photon most often arrives after the $T_{\text{GBM,90}}$ (right panel of Figure 25), there is also a large fraction where it arrives during this interval. Interestingly, this fraction seems to vary slightly with energy: below 1 GeV more than half of the events arrive during $T_{\text{GBM,90}}$, whereas 70% of events arrive later at energies above 10 GeV. We also note that there is little difference between the short and long GRB populations

regarding normalized arrival time, although the maximum energy of sGRBs appears on average a bit lower. This could link the highest-energy photons to a process common in the two scenarios, and independent from the progenitor, such as external shocks in the circumburst medium. This will be further discussed in Section 6.

As of writing, the observing pattern adopted after the solar array malfunction (see Section 2.1) could potentially lead to a loss in GRB detections, as events occurring in the unobserved hemisphere of the sky may not be monitored until long after the event. However, studying the detections reported in this catalog, no such event was detected even during the previous observing strategy. The exception to this are the ARRs, as they could bring bright events into the FOV. It is thus likely that the detection rate will be slightly impacted while ARRs remain disabled. Even if the LAT detects late-time emission, without ARRs there is the risk of losing high-energy prompt emission.
<figure><img src="image_26.png" /><figcaption>**Figure 23.** Isotropic luminosity $L_{iso}$ calculated in the 100 MeV–100 GeV rest-frame energy range for the 34 GRBs in our sample with measured redshift. The left panel shows $L_{iso}$ vs. the time elapsed since the trigger. GRB 090510 is the only short GRB with known redshift, and it is marked with black stars. In the middle panel $L_{iso}$ is plotted against the time elapsed since the trigger in the rest frame. In the right panel, we show $L_{iso}$ divided by the isotropic energy (calculated in the 1 keV–10 MeV rest-frame energy range) vs. the time elapsed since the trigger in the rest frame. We mark GRB 160623A, which represents an outlier of the distribution (see text for explanation), with orange crosses. The solid line shows a linear fit, giving a decay index of $1.25 \pm 0.03$. For comparison, we also show a dashed line with decay index 10/7.</figcaption></figure>

<figure><img src="image_27.png" /><figcaption>**Figure 24.** Fraction of GRBs with the highest-energy photon detected above selected threshold energies (250 MeV, 500 MeV, 1 GeV, 5 GeV, 10 GeV, 50 GeV; green solid line). The distribution of the source-frame-corrected energies for the redshift sample is indicated with the dashed green line. The dashed black line denotes a linear fit to the values corresponding to the center of each bin.</figcaption></figure>

## 5.2. Onset and Duration of the High-energy Emission

As in the 1FLGC, we can firmly establish the general trend that high-energy emission from GRBs tends to have delayed onset and longer duration as compared to emission at lower energies. However, Figure 7(c) also shows that when high-energy emission is detected, it starts during the prompt phase in >60% of the cases. The majority of the other GRBs were outside the LAT FOV at trigger time, meaning that the fraction is likely even higher. This is an interesting result that provides valuable input to models of the emission mechanisms and will be discussed further in Section 6.

Figure 7 shows how varied the difference between $T_{GBM,0.5}$ (calculated in the 50–300 keV range) and $T_{LAT,0}$ (100 MeV–10 GeV) can be. In some cases, the LAT emission is completely contemporaneous with the GBM. In other cases, the LAT emission starts hundreds or even thousands of seconds later. Considering only bursts that were in the LAT FOV at the time of trigger, there are a number of cases where the high-energy emission came much later than the one at lower energies. For example, GRB 160503A shown in the right panel of Figure 6 was at $\theta = 25.1^{\circ}$ at $T_{LAT,0}$, and it remained in the FOV for over 2 ks without any high-energy emission being seen. The first detection instead came much later, at 5 ks. These extreme delays are much longer than the ones seen in the 1FLGC and represent a new result in the 2FLGC.

In addition to greater delays, we now report much longer durations. In the 1FLGC, the longest duration reported was >800 s for GRB 090902B. In the 2FLGC, many GRBs have durations of order $10^3$ s, with the longest duration being 35 ks (GRB 160623A). In general, the durations have increased also for most bursts contained in the 1FLGC, likely due to better sensitivity as a result of Pass 8 (we note, however, that the duration estimates were made using a different technique as described in Section 3.6.2, so the numbers are not directly comparable).

## 5.3. Comparison with GBM Flux and Fluence Distribution

As already presented in Section 4, the LAT-detected GRBs tend to sample the upper range of the GBM flux and fluence distribution (Figure 13). At the high end of the GBM fluence distribution, the LAT detects a high fraction of the GRBs; above $\sim 10^{-4}$ erg cm$^{-2}$ the two distributions practically overlap (left panel). The few additional bursts seen by the GBM could be explained by some GRBs being outside the LAT FOV or at high zenith angles. The effect of the $\theta$ angle is further investigated in Ajello et al. (2018), who show that this is the main factor determining LAT detectability. The bias toward high GBM flux and fluence for LAT-detected bursts is therefore not related to a difference in sensitivity between the instruments.
<figure><img src="image_28.png" /><figcaption>**Figure 25.** Energy vs. arrival time for the highest-energy photon of each GRB. In the right panel, the arrival time is normalized to the duration ($T_{90}$) calculated in the 50–300 keV energy range (indicated by the dashed vertical line). Blue and red circles represent long and short GRBs, respectively.</figcaption></figure>

Looking at the high end of the flux distribution (middle panel of Figure 13), the fraction of GRBs not seen by the LAT is higher. This shows that the flux might not be such a good indicator for high-energy emission. While the fluence is a measure of the total energy output, the flux simply shows the “strength” of the peak. The light curves also show clear differences, and peaks at low energies are not necessarily mirrored at high energies. An example of this is GRB 180728A, where the GBM peak flux was about 230 photons $cm^{-2}s^{-1}$ (Veres et al. 2018), but no detection was made at high energy even though the burst was inside the LAT FOV ($\theta \sim 35^{\circ}$). This is likely due to the low value of $E_{peak}$, 80 keV; indeed, the flux was dominated by energies below 50 keV.

However, the LAT has also detected GRBs that have relatively low fluence in the GBM. These outliers are predominantly short GRBs, where the low fluence is naturally explained by the short duration. The fluence distribution of sGRBs in the GBM is also shifted to lower values overall. Also for sGRBs there is a tendency for the LAT to sample the higher fluence end of the GBM population, but this is much less marked than for the IGRBs.

The result that the LAT detections are biased toward the brighter GRBs was clear already in the 1FLGC. However, with the larger sample, the picture presented here becomes more nuanced. The sGRBs show that high-energy emission can be produced even at lower fluence, raising the question of why not more IGRBs are detected. There may be differences in the emission mechanisms or environments between the classes that explain why low-fluence sGRBs are more likely to be detected in high-energy emission than IGRBs of similar fluence. For instance, IGRBs are expected to have a denser circumburst medium, as they are coupled to massive stars that have strong stellar winds.

The possibility of GeV emission from sGRBs is particularly interesting in the light of GW170817 and the associated sGRB 170817A (Abbott et al. 2017a, 2017b; Goldstein et al. 2017). This event was outside the LAT FOV and had a fairly low

fluence of $2.8 \times 10^{-7} \text{ erg cm}^{-2}$. However, Figure 13 shows that the LAT has detected GRBs with a similar fluence. This is very promising in view of the upcoming observation period scheduled to start in spring 2019, showing a strong potential for LAT detections of similar events (see also Ajello et al. 2018).

### 5.4. Origin of Emission below 100 MeV

In our sample almost twice as many GRBs are detected above 100 MeV (169) as in the 30–100 MeV LLE range (91). Several studies have found evidence for a separate spectral component behind the emission above 100 MeV, and high-energy cutoffs between the LLE and LAT energy ranges have also been seen (e.g., Vianello et al. 2018). The behavior of the emission below 100 MeV further shows more similar temporal behavior to the GBM range than do the data above 100 MeV (see Figures 7 and 12). The LLE-only GRBs in our sample would then be the result of the low-energy emission being strong enough (and/or $E_{peak}$ at high enough energy) to extend into the LLE range. As the burst evolves, $E_{peak}$ moves to lower energies, and the emission in the LLE range will therefore appear to fade before that at lower energies. This explains the fact that the duration at 30–100 MeV is almost always shorter than the one measured by the GBM. This picture also explains the fact that the LLE emission appears to start earlier than the GBM emission in a few cases (left panel in Figure 12). As the duration is shorter in the LLE range, this will naturally make the value of $T_{LLE,05}$ shorter as well. Again, it is just a sign that the emission in the GBM range lasts much longer than the one in the LLE range.

For GRBs detected in both LAT and LLE energy ranges, the duration is generally shorter in the 30–100 MeV band. While direct comparison of the $T_{LLE,90}$ and $T_{LAT,100}$ should be done with caution, the large differences in duration here clearly point to an intrinsic origin rather than observational bias. For example, in the LLE range the effective area is up to a factor of 2–3 lower than that above 1 GeV; however, the $T_{LAT,100}$ can be more than an order of magnitude longer than the $T_{LLE,90}$.
Although there is much to suggest that the LAT emission is often the result of a separate component, we caution that there are also cases where the GRB spectrum is seen in the full range from keV to GeV and well fit by a single component extending also to energies above 100 MeV (e.g., Abdo et al. 2009c; Axelsson et al. 2012). So while the dominance of LAT-only detections supports a separate emission process >100 MeV in general, proper spectral analysis must be made to draw conclusions about an individual burst.

### 5.5. Energetics in the Prompt and Afterglow Phases

We first investigate the energy output in different time windows. Figure 14 allows us to qualitatively compare the energy output during the total duration of each GRB as measured in the 50–300 keV (GBM) and 100 MeV–10 GeV (LAT) ranges. In the “GBM” duration window, the flux (top left panel) is inversely proportional to the duration. A possible interpretation is that a limited amount of total energy is available for any given GRB. A longer duration then means that the average flux will decrease. Comparing the fluences (top right panel), long and short GRBs are clearly separated, with all but one of the sGRBs having fluence below $10^{-6}$ erg cm$^{-2}$. Within each group, fluence is independent of duration, suggesting that the separation arises as a result of different energy budgets, as might be expected if the groups arise from two different progenitor scenarios.

In contrast, long and short GRBs show no clear separation in the “LAT” duration window (bottom panels). While the populations occupy slightly different regions of the plots, there is a smooth transition between them. This may indicate that the emission in this window is dominated by the afterglow and that this process is similar regardless of progenitor. In the right panel, the fluence shows a tendency to increase with duration. This could have many explanations: a variable energy budget, meaning that more energetic afterglows last longer; an effect of differences in the circumburst medium, leading to varying radiative efficiency; or variations in the viewing angle with respect to the jet. We caution that there is likely also an observational bias: GRBs with favorable conditions can be studied longer; these will then have both a higher number of detected photons (i.e., greater fluence) and a longer measured duration.

Looking closer at the energetics, we can compare the LAT 100 MeV–10 GeV fluence to that of the 10 keV–1 MeV range in the GBM. We make this comparison in the “GBM” time window (left panel of Figure 15), which can be seen as a proxy for the prompt phase of the burst. We again stress that no joint spectral fits have been made and that the fluences are calculated based on different spectral fits in each range: the GBM fits are taken from the FGGC, while a PL fit is used for the LAT range (see Section 3.6). In the 1FLGC, it was found that the high-energy fluences from the joint GBM-LAT spectral fit agreed with those from LAT-only analysis, with a small discrepancy seen only for the brightest GRBs. Therefore, we do not expect any significant bias in our results, but a more thorough study needs to be performed before detailed conclusions can be drawn. Such analyses will be presented in future publications.

As found in the 1FLGC, the bulk of the population emits most of its radiation in the GBM energy range, with the high-energy emission reaching ≲20% of that in the lower band. However, given the larger sample in the 2FLGC, it is clear that this mainly applies to long GRBs. The short GRBs are more

clustered around the line of equality, and several have higher fluence in the 100 MeV–10 GeV range than in the 10 keV–1 MeV one. We therefore conclude that short GRBs in the prompt phase tend to emit a relatively larger fraction of their energy at higher energies, as compared to long GRBs.

In the next step we compare the relative energy output in the 100 MeV–10 GeV range during the “GBM” (prompt) and “EXT” (afterglow) time windows of each burst (right panel of Figure 15). Similar to the results found in the 1FLGC, most GRBs are clustered around the line of equality, meaning that in the LAT energy range comparable amounts of energy are released in the prompt and afterglow phases. For long GRBs there is a large fraction where the fluence during the afterglow phase is greater than that during the prompt phase. The short GRBs instead tend to have equal fluence in the two time windows.

In the 1FLGC, there appeared to be a correlation between the photon index and GBM duration, with shorter bursts tending to have harder spectra. This was coupled to the general property that short GRBs on average have harder spectra. Figure 16 shows that this overall trend is not held up in the current larger sample. Furthermore, although the short GRBs on average have slightly harder spectra in the “GBM” time window (top left panel), this difference is not seen in the “EXT” window (top right panel). This could merely be an indication that the GRB spectra are sometimes “contaminated” by emission from lower energies, i.e., the $\beta$ value of the Band function. This value is typically just below -2, and the distribution has a long tail to lower values (Gruber et al. 2014). The larger number of IGRBs means that we are more likely to see the effect in this population. Comparing the spectral index between the “GBM” and “EXT” windows, the small difference ($\Delta\Gamma \sim 0.5$) might indicate that we are seeing the early afterglow already in the “GBM” window. This possibility will be discussed further in Section 6.

The bottom left panel of Figure 16 shows that the spectral index is not dependent on flux in the “GBM” window. Instead, the bottom right panel would appear to suggest a correlation between flux and photon index in the “EXT” time window, such that stronger emission is coupled to softer spectra. This is especially obvious in bursts with lower flux. However, we caution that this is likely an observational bias. Weaker bursts will only be detected if they have harder spectra.

#### 5.5.1. Intrinsic Energetics

In the 1FLGC, there were 10 LAT-detected GRBs with measured redshift. In the 2FLGC, the number has increased to 34, allowing us to better study their intrinsic properties. Already Figure 13 suggests that the LAT GRBs are among the brightest. Comparing with the Swift and GBM samples in Figure 18, it is clear that this is an intrinsic property and that the LAT preferentially detects brighter GRBs, regardless of redshift. The two panels in Figure 18 show that in order for the LAT to detect a GRB at high redshift, such a burst must be intrinsically very bright. For example, all detected GRBs at $z > 1$ have $E_{\text{iso}} > 10^{52}$ erg, while detecting a GRB at $z > 4$ requires $E_{\text{iso}} > 10^{54}$ erg. After 10 yr, GRB080916C is still the most distant GRB observed by the LAT, at $z = 4.35$.
<figure><img src="image_29.png" /><figcaption>**Figure 26.** Energy of the highest-energy photon calculated in the rest frame vs. $E_{\text{iso}}$ in the 1 keV–10 MeV (left panel) and in the 100 MeV–10 GeV (right panel) energy ranges for each GRB with redshift estimation. Blue and red circles represent long and short GRBs, respectively.</figcaption></figure>

5.5.2. *The Highest-energy Photons*

We now turn to the highest-energy photons detected by the LAT, as presented in Figure 24. Although less than 80% of GRBs do not reach above 5 GeV in the observer frame, the figure indicates that higher energies are more common in the source frame, with ~15% of GRBs reaching source-frame-corrected energies above 100 GeV. We note, however, that as there is only one sGRB with measured redshift, the rest-frame distribution is dominated by IGRBs. Interestingly, the fraction of GRBs detected drops smoothly as the threshold energy in the rest frame increases. It is tempting to connect this behavior to the underlying spectral distribution, such that we are seeing the limit determined by the intrinsic spectral shape, which seems to be similar for all GRBs (see Figure 16). In bright bursts more high-energy emission is produced, allowing GeV emission to be observed. Faint bursts will produce too little GeV emission, and the LAT will only detect MeV photons.

Although it is clear that the high-energy emission in general
can last a long time, the highest-energy photon in a given GRB
in some cases arrives relatively soon after the trigger. For
example, this photon arrives within <2 s from the trigger time
for ~50% of the sGRBs in the left panel of Figure 25. The
right panel of Figure 25 also shows that among the GRBs with
a maximum photon energy >10 GeV, the highest-energy
photon is in some cases detected before the emission in the
50–300 keV energy range is over. Particle acceleration in
GRBs must be efficient in order to produce such high-energy
gamma-rays within such a short time. Considering internal
opacity of the jet outflow also leads us to conclude that a high
bulk Lorentz factor is necessary in order for us to be able to
detect these photons. Indeed, these considerations have been
used to estimate a Lorentz factor well above 500 in several
LAT-detected bursts (see, e.g., Abdo et al. 2009a, 2009c).

In the 1FLGC, there was a trend for the photons with highest
source-frame energy to appear in the GRBs with highest 1 keV–
10 MeV $E_{\text{iso}}$ (Figure 21 in the 1FLGC), and the one short GRB
with redshift (GRB 090510) did not follow the pattern seen in long
GRBs. However, these results were based on only nine GRBs.
With the greater statistics in the 2FLGC, the trend not only persists

but also is extended over a greater range, as shown in the left panel
of Figure 26. This indicates that the maximum rest-frame energy
produced is simply a function of the total energy output. The right
panel instead shows the relation between the highest energy and
$E_{\text{iso}}$ calculated during the “GBM” time window in the LAT energy
range (100 MeV–100 GeV). In this range, the correlation is less
obvious and the points are more scattered. As the energy output in
the LAT range is much smaller than at lower energies (see the left
panel in Figure 15), it is perhaps not surprising that the low-energy
$E_{\text{iso}}$ is a more reliable estimator of the energy budget. However,
there are a small number of GRBs that deviate from the general
trend, showing low rest-frame energies despite a large $E_{\text{iso}}$. More
detailed studies are required in order to understand this behavior.
In both panels the region of GRB 090510, which is still the only
short GRB with a redshift measurement, is populated also by long
GRBs. This again points to similar conditions for long and short
GRBs. Finally, the high-energy gamma-rays from GRBs offer a
valuable way to probe more physics than the burst itself. For
example, the interaction of >10 GeV gamma-rays from sources at
cosmological distance with optical and UV photons of the
extragalactic background light (EBL) causes absorption of gamma-
rays, modifying the high-energy part of the spectrum. High-
redshift GRBs can thereby serve as a probe of EBL opacity (Desai
et al. 2017). For example, EBL attenuation could contribute to the
differences between the observer and rest-frame distributions in
Figure 24.

5.6. *High-energy Emission Temporal Profile*

In the 1FLGC the study of the temporal decay of the
extended emission led to the hypothesis that this is part of an
early afterglow component very similar to the afterglow
component seen in the X-ray band. The stable decay index
with a value around $\alpha = 1$ for the small sample of bursts in the
1FLGC was taken as support for adiabatic expansion within the
fireball scenario. The index reported here remains centered
around $\alpha = 0.99$. If the extended component is fast-cooling
afterglow emission from either a radiative or an adiabatic
fireball in a constant-density environment, the decay indices
(given a photon index of $-2$) would be $10/7$ or $1$, respectively
(Katz & Piran 1997; Sari 1997; Ghisellini et al. 2010). The
<figure><img src="image_30.png" /><figcaption>**Figure 27.** Light curve of GRB 160509A, together with the fit using the Willingale et al. (2007) model (solid line). The black circle indicates the end time of the plateau emission and its corresponding flux.</figcaption></figure>

peak in the decay index distribution would thus primarily suggest the adiabatic scenario. Nevertheless, the range of measured indices allows also for a radiative fireball in at least a few individual bursts.

The luminosity distribution seen in Figure 23 shows that the GRBs with known redshift all have similar decay indices. Yet although the slope of the luminosity decay is similar among the GRBs, its distribution versus the elapsed time since the trigger time is scattered (left panel). In the middle panel, the time axis for each light curve has been corrected to the time in the rest frame, which reduces the spread. If the luminosity is divided by the isotropic energy of each burst, measured in the 1 keV–10 MeV range, the resulting light curves cluster even more (right panel in Figure 23), as reported by Ghisellini et al. (2010) and Nava et al. (2014). With higher statistics with respect to these studies (which were limited to a few LAT GRBs observed up until 2013), the clustering points to a strong correlation between the isotropic luminosity of the high-energy emission and the prompt isotropic energy calculated in the range 1 keV–10 MeV.

In the analysis performed in Sections 3.6.2 and 4.7, some light curves are found to be best fit with a spectral break. Such features may be related to the X-ray afterglow plateaus, which are well established, but their nature is still debated. To investigate this, we analyze all LAT GRBs with $TS \ge 64$, with a sufficient number of photons. This choice of TS is dictated by the need of observing the γ-ray afterglow for a sufficiently long time to allow us to characterize the nature of the deviation from a PL. We then determine whether the light curves can be fitted within the phenomenological Willingale et al. (2007) model, used to fit the X-ray plateaus (the details of the analysis will be presented in M. Dainotti et al. 2019, in preparation).

We find that four GRBs (GRB 160509A, GRB 090902B, GRB 150627A, GRB 160816A) show a late-time flattening reminiscent of the X-ray plateaus (their light curves can be seen in Figures 20 and 22). The light curve of GRB 160509A, together with the Willingale et al. (2007) model, is presented in Figure 27. From this resemblance we infer that there is indication of the existence of a plateau in a few Fermi-LAT light curves. Our recovered values for $T_γ$, which is the time at the end of the plateau emission, are smaller than those typically found in the X-ray regime (from hundreds to a few thousand seconds; Willingale et al. 2007). This suggests that more

energetic bursts have smaller plateau durations. Although a one-to-one comparison is difficult owing to the lack of corresponding X-ray observations for those bursts, it is nevertheless interesting to mention this trend for future reference and possible new observations.

## 5.7. Prospects for GRB Detections at VHE

The detection of VHE emission from GRB 190114C reported by the MAGIC Collaboration (Mirzoyan 2019) opened up a new era in the observations of GRBs. The bright source was seen across the electromagnetic spectrum, from radio to VHE. While studies of this burst are certain to yield interesting results, more VHE detections will be necessary to build an understanding of the mechanisms behind such emission.

The Cherenkov Telescope Array (CTA; Actis et al. 2011; The CTA Consortium et al. 2019) will mark an important milestone in the VHE range, covering the whole sky with unprecedented sensitivity. The CTA can steer to a GRB location in less than 30 s after receiving a trigger notice from a space-based instrument (see, e.g., Gilmore et al. 2013; Inoue et al. 2013). In an optimistic case, one can therefore assume that the CTA is on target within ~50 s after the onset of the burst.

The GRB sample in the 2FLGC allows us to make a simple and phenomenological estimate of the detectability by VHE instruments. For this exercise, we use the observed LAT spectra as input for CTA simulations, considering only the 33 IGRBs with redshift measurements. We performed an ad hoc spectral analysis in a 2 hr interval starting at $T_{LAT,0} +50$ s in the energy range from 100 MeV to 100 GeV. The LAT data are fit using a PL spectrum. The output spectral parameters (flux, $F_0$, and photon index, $\Gamma$) were then used as input for the CTA tools software package `ctools`,<fnref n="76" /> version 1.5.2, with the instrument response functions (IRFs) Prod2⁷⁷ referred to the CTA configuration in the Southern Hemisphere site. The CTA simulation requires as input the GRB location and spectral model and gives as output a TS map showing whether the source is detected or not. The extrapolation into the CTA range assumes either a PL or a cutoff PL (CUTPL) spectral model. The background is already taken into account in the used IRFs. We performed simulations over an energy range of 50 GeV–50 TeV, assuming that the CTA received the LAT detection notice with several different time delays: 0, 0.5, 1, 2, 4, 6, and 10 hr. We considered a 2 hr CTA observational interval following the notice time, assuming a flux temporal decay of $F = F_0 t^{-\alpha}$, where $\alpha$ is the temporal decay index taken from Table 4. For the 0 hr delay we used all 33 GRBs, while for all the other analyses the number of cases was limited to the 27 GRBs where $\alpha$ could be constrained.

An important factor to consider in the VHE range (>50 GeV) is that the gamma-ray spectrum can be attenuated by γγ absorption due to the annihilation of high-energy photons with photons from the EBL. In order to include the EBL attenuation in our simulation, the cutoff energy of the CUTPL model was chosen as the critical energy defined by Ackermann et al. (2011). For the GRBs at redshift >1 (20 GRBs), we fixed the critical energy at 30 GeV (following Figure 2 of Ackermann et al. 2012b), while for the others at redshift <1 (13 GRBs) the

<footnote n="76">76 http://cta.irap.omp.eu/ctools/users/index.html</footnote>

<footnote n="77">77 https://www.cta-observatory.org/science/cta-performance/</footnote>
<figure><img src="image_31.png" /><figcaption>**Figure 28.** TS evolution with respect to the time delay (in hours) of the CTA alert reception, obtained from the analysis with a PL model, with respect to the beginning of the LAT-detected extended emission. The figure shows the TS values for 27 GRBs. For clarity, TS values less than 1 are fixed to 1.</figcaption></figure>

critical energy was chosen following Ackermann et al. (2012b): $E_{crit}(z) = 170(1+z)^{-2.38}$ GeV. Spectral analyses of the CTA simulated data were then performed using both PL and CUTPL spectral models.

The analyses performed without any time delay (0 hr) reveal that the CTA is able to detect all 33 GRBs with high TS, as shown in Figure 28 for the PL case. This shows that CTA telescopes would readily detect GRBs given ideal observational conditions during night time, i.e., assuming that the CTA was already on target or started observing immediately after a GRB trigger. However, the LAT data generally require more than 6 hr to be downlinked and analyzed on the ground. Figure 28 also shows the TS evolution with respect to the time delay between the LAT detection and the CTA alert reception time in the case of a PL model. The TS values decrease rapidly with time, and given a 10 hr delay from the LAT detection, the CTA can detect only two GRBs among 27 (with $TS > 25$). In the case of the CUTPL model, all GRBs are detected only in the most optimistic observing condition, and the TS values fall below 25 for all 27 GRBs already with a 0.5 hr delay. These results show the strong need for rapid alerts and precise localizations to enable follow-up observations at VHEs.

We stress that the analysis presented here is the result of a very simplified scenario using a single spectral component. Extending this study to give more reliable estimates of the detectability with VHE instruments will require both more realistic models for the high-energy emission mechanisms in GRBs and using proper theoretical models of EBL.

## 5.8. LAT sGRBs and Implications for GW Observations

The recent detection of GW170817 and the coincident GRB 1701817A clearly showed a coupling between merging neutron stars, gravitational waves, and short GRBs. This has made the rapid search for counterparts also in the MeV–GeV range a high priority for the LAT. The LAT Collaboration has previously reported on follow-up studies of several GW triggers: GW150914 (Ackermann et al. 2016), LVT 151012

and GW151226 (Racusin et al. 2017), and GW170817 (Ajello et al. 2018). The first three triggers were associated with mergers of two black holes, so no electromagnetic counterpart was expected. In the last case, the location could only be observed by the LAT starting 1100 s after the GBM trigger time, and again no detection was made. Nevertheless, these events have shown the importance of a dedicated pipeline to search for high-energy counterparts for a larger class of transient events than just GRBs. The LTF pipeline described in Section 3.2 has formed the baseline for these searches. In addition, dedicated searches have been developed to cover various time intervals (also very long ones) around the GW triggers, which are thoroughly described in Vianello et al. (2017). A particular challenge is the large area given by the LIGO probability map (several hundred square degrees), much larger than for a typical GRB.

As already said, due to observational constraints, the LAT could observe the location of GW170817/GRB 170817A only after more than 1 ks. In the sample of 14 sGRBs presented here, most events are only detected up to a few seconds after the GBM trigger. However, GRB 160702A was not seen until ~1900 s after the trigger and lasted for about 300 s. Also GRB 170127C has an unusually long duration: from ~650 to ~2900 s after the trigger. Neither of these events was especially bright (flux < $10^{-5}$ erg cm$^{-2}$ s$^{-1}$), showing that the LAT has the potential to find counterparts even to GW events that occur outside the FOV. Overall, the LAT has found high-energy counterparts to 5% of the sGRBs seen by the GBM. A critical discriminating factor that has been carefully addressed by Ajello et al. (2018) is the fact that the LAT can detect fainter sGRBs the earlier the observation starts, ideally within 100–200 s post-trigger. They calculated that the LAT can observe (either a detection or upper bound) ~23% of the full-sky sGRB population within 100 s from their GBM trigger. Moreover, assuming that the LAT will have the same efficiency for GRB/GW triggers and a rate of joint GBM/GW events of 1 (2) per year, they obtained a ~5% (~10%) probability of detecting one or more GRB/GW events with the LAT in 1 yr, respectively. This was derived assuming that GRB/GW events will be representative of the entire GBM-detected sGRB population when observed in gamma-rays. Indeed, two other important factors affecting the detectability of such events by the LAT might be the off-axis viewing angle and the GRB distance.

Recently, von Kienlin et al. (2019) presented a search for GRBs in the entire 10 yr Fermi-GBM sample that show similar characteristics to GRB 170817A, i.e., a nonthermal pulse followed by a thermal component. They identified 13 candidates and concluded that the observed similarities likely arose as a result of the same processes that shaped the gamma-ray signal of GRB 170817A. We examined their sGRB sample and found that GRB 120915A has also been detected by the LAT and is included in the 2FLGC. The short duration in the LAT ($T_{LAT,100} = 0.6 \pm 0.3$ s) clearly ties the high-energy emission to the prompt phase. This again shows that the LAT is likely to detect such events as long as they are in the FOV at the time of the GW trigger, as in the case of GRB 120915A (with $\theta = 11^\circ$).

# 6. Theoretical Implications

Observations of GRBs with Fermi-LAT have significantly broadened our understanding of the nature of high-energy
emission from these powerful transients and have led to renewed theoretical activities to model this emission (see, e.g., Gehrels & Razzaque 2013, for a review). Already the 1FLGC hinted at several new prominent characteristics of the high-energy emission: delayed onset compared to the GBM trigger; a hard spectrum, often requiring an additional component when modeling together with GBM data in the prompt phase; variable emission during the early prompt phase; and temporally extended emission beyond the prompt phase without measurable short-timescale variation. Although no joint GBM-LAT spectral fits have been performed for the 2FLGC, with its greatly enhanced sample the results presented here have firmly established at least the temporal properties as being general to the population as a whole. Modeling these characteristics gives us clues to the particle acceleration sites, environment, and radiation processes, while measurements of fluence, redshift, and energetics help us to classify the LAT-detected GRBs in the context of the broader GRB population.

Simultaneous multiwavelength observations of a handful of LAT-detected GRBs, e.g., GRB 090510 (Ackermann et al. 2010a), GRB 110731A (Ackermann et al. 2013a), and GRB 130427A (Ackermann et al. 2014; Maselli et al. 2014), have further helped constrain the origin of a few high-energy emission features, such as the temporally extended emission. Regardless of the origin of high-energy emission, merely the fact that we detect it implies that the bulk Lorentz factor of the GRB fireball (Lithwick & Sari 2001; Granot et al. 2008; Hascoët et al. 2012) is large. In particular, the fireball can be completely transparent to high-energy emission for sufficiently large bulk Lorentz factors (Razzaque et al. 2004). Furthermore, detection of $\gtrsim 10$ GeV photons from GRBs has constrained models of the extragalactic background light (Abdo et al. 2010b) and the models leading to violation of Lorentz invariance (Abdo et al. 2009b). Finally, the detection of the short GRB 170817A by Fermi-GBM (Goldstein et al. 2017) in coincidence with the gravitational wave event GW170817 (Abbott et al. 2017a, 2017b, 2017c) has proven that at least some short GRBs originate from mergers of binary neutron stars.

The search for the origin of GeV emission from GRBs has intensified since Fermi was launched and the LAT has detected an increasing number of GRBs. The results presented here make these efforts even more promising, as the large sample of the 2FLGC shows that high-energy emission can be seen also from some fainter GRBs (see Figure 13). In the following, we will examine how GeV emission from GRBs can be used to test theoretical predictions. These examples are not exhaustive—rather, they are meant to illustrate various ways in which LAT observations can inform theoretical models.

## 6.1. Delayed Onset and Long Duration: Challenges for Emission Models

Synchrotron self-Compton (SSC) models for GeV emission are a natural extension of the synchrotron model for keV–MeV emission from GRBs, in the context of the internal shock scenario (see, e.g., Bošnjak et al. 2009; Wang et al. 2009; Ackermann et al. 2010a). However, in a simple SSC model the Compton emission can be delayed compared to synchrotron emission by at most the flux variability timescale ($\lesssim 1$ s for long GRBs). Such models therefore cannot account for the $\gg 1$ s delay observed in a large fraction of LAT GRBs (see Figure 7). SSC emission from late internal shocks, with much longer

variability time at GeV than in MeV, can explain delayed LAT onset in some cases (see, e.g., Toma et al. 2009; Daigne et al. 2011).

A possible alternative for GeV emission during the prompt phase is hadronic emission, either proton-synchrotron radiation and radiation from $e^{\pm}$ pair cascades (Razzaque et al. 2010) or photohadronic interactions and associated cascade radiation (Asano et al. 2009). In contrast to many SSC models, hadronic models can account for delayed onset of GeV emission via a hard PL component with its flux increasing over time owing to proton acceleration and cooling of cascade particles. Hadronic models, however, require $\gtrsim 100$ times more energy in the hadrons than in leptons, which is challenging to produce in GRBs. Furthermore, neither SSC nor hadronic emission mechanisms active during the prompt phase are able to explain the extended GeV emission beyond the keV–MeV $T_{GBM,90}$.

The delayed onset, smooth light curve, and temporally extended emission $\gtrsim 100$ MeV have led to the development of an early afterglow interpretation of high-energy radiation from the external forward shock (Kumar & Barniol Duran 2009; Ackermann et al. 2010a; Corsi et al. 2010; Ghisellini et al. 2010; Razzaque 2010). The PL spectrum and temporally extended emission are then interpreted as synchrotron radiation from the external shock. Modeling of simultaneous multi-wavelength (optical to GeV) data has shown this approach to generally fit observations for several GRBs, such as GRB 090510, GRB 110731A, and GRB 130427A. However, detection of $\gtrsim 10$ GeV photons at late times from a number of GRBs (see Figure 25) is not easy to explain with the synchrotron interpretation (nor are the VHE photons from GRB 190114C seen by MAGIC). The maximum photon energy from electron-synchrotron radiation is $h\nu_{max} \approx 0.24 \phi^{-1} (1+z)^{-1}\Gamma(t)$ GeV, where $\phi^{-1} \lesssim 1$ is the acceleration efficiency for electrons and $\Gamma(t)$ is the bulk Lorentz factor (see, e.g., Razzaque 2010).

In light of the difficulties in simultaneously matching both the delayed onset and long duration of GeV emission, alternative models have also been proposed. For instance, several authors have pointed out that the prompt MeV emission should cause the creation of $e^{\pm}$ pairs in the medium ahead of the blast wave (e.g., Thompson & Madau 2000; Mészáros et al. 2001; Beloborodov 2002). This pair loading led Beloborodov (2005) to suggest that GeV emission arises owing to inverse Compton scattering of the prompt MeV radiation streaming through the external medium.

In Beloborodov et al. (2014), this scenario was shown to predict a delayed onset of $\sim 1–10$ s, in agreement with the majority of LAT GRBs. However, it is unclear whether also the much longer ($>100$ s) delays seen in a few GRBs can be explained. This would require either a very low Lorentz factor of the blast wave (not easy to reconcile with GeV emission) or a suppression of the blast-wave luminosity until it has reached a very large radius. Likewise, while the model can explain a GeV duration several times greater than that at lower energies, it may face difficulties for the LAT GRBs where the GeV duration is more than an order of magnitude longer (see panel (c) of Figure 7). Observational tests so far have also predicted rather large densities of the circumburst medium (comparable to stellar winds of $\sim 10^{-5} M_{\odot} \text{yr}^{-1}$), which may prove challenging as more sGRBs are detected, as binary neutron stars are not expected to have a dense surrounding medium.
## 6.2. Fluence and Energetics of LAT GRBs: Probing the Central Engine

LAT GRBs constitute a subsample of the GBM GRBs with a flux distribution peaking at a higher value than the complete GBM sample (see Figure 13). This is also true for the fluence as evident from the same figure. On the other hand, the redshift distribution of the LAT GRBs is rather typical of Swift-detected GRBs (see Figure 17). As a result, the isotropic energy release and peak luminosity (Figure 18) are also typically higher for the LAT GRBs across the redshifts. These results have implications for the GRB central engines.

The nature of the GRB central engine is a fundamental question in astrophysics. A newly born magnetar is one of the two leading hypotheses (the other is a newly born black hole). The vacuum/force-free electromagnetic spin-down power of a neutron star with surface (equatorial) dipole magnetic field $B_{dip}$ and spin period $P$ is given by (Spitkovsky 2006)

$$ \dot{E}_{\text{FF}} \approx \frac{4\pi^4 B_{\text{dip}}^2 R^6}{P^4 c^3} \sim 10^{49} \left( \frac{B_{\text{dip}}}{10^{15} \text{ G}} \right)^2 \times \left( \frac{P}{1 \text{ ms}} \right)^{-4} \left( \frac{R}{10 \text{ km}} \right)^6 \text{ erg s}^{-1}, \quad (11) $$

where $R$ is the neutron star radius. The maximum isotropic γ-ray power is then ~$10^{49}-10^{51}$ erg s⁻¹, after taking into account a jet-beaming factor of up to ~200 for a ~6° jet opening angle. In comparison, the measured peak luminosity for a number of LAT-detected GRBs exceeds these values, e.g., GRB 080916C (~$6 \times 10^{53}$ erg s⁻¹) and GRB 090510 and GRB 160625B (~$4 \times 10^{53}$ erg s⁻¹). The observations of LAT GRBs therefore significantly constrain the parameter space of the models with a magnetar as the GRB central engine (see, e.g., Duncan & Thompson 1992; Usov 1992; Metzger et al. 2011).

## 6.3. Closure Relations: Test of the Standard External Forward Shock Model

As described in Section 6.1, a possible explanation for the temporally extended emission is synchrotron emission from an external forward shock developed by collisions between outgoing ejecta and the circumburst medium. According to the model, a PL distribution of electrons radiates synchrotron emission, and the spectrum and corresponding light curve are described as a series of BPLs (e.g., Sari et al. 1998; Chevalier & Li 2000; Zhang et al. 2006; Gao et al. 2013). The indices of each segment of both spectrum and corresponding light curve are either constant or a function of the electron spectral index, and the relation between the spectral and temporal indices is called the “closure relation.” The closure relations have different forms depending on the physical conditions, such as the electron cooling regime, surrounding environment profile, electron spectral index, and jet geometry.

A systematic analysis of 59 GRBs testing a set of closure relations against temporal and spectral indices computed from LAT extended emission was performed in Tak et al. (2019, submitted). Among the analyzed GRBs, 81% of them can be classified by a set of closure relations. A third of these allow both interstellar medium (ISM) and wind environments, while the remainder favor an ISM environment. The observed spectra and light curves of the extended emission in LAT GRBs are thus generally consistent with the standard external forward

shock model in most cases, but there is still a significant fraction that cannot be explained by this model. The GRBs that do not match any closure relation tend to show a slower temporal decay (~$\alpha < 1$). This implies the need for continuous energy injection or other physical sources to sustain the fluxes.

As already mentioned in Section 5.6, some light curves show a break or even hints of a plateau at late times. Such features can cause problems for the models used above, and instead a more phenomenological approach needs to be adopted. Using the parameters derived from fitting the light curves with the Willingale et al. (2007) model, we tested also these GRBs against the closure relations. Unfortunately, due to the large uncertainties, no conclusive results could be reached at least at this stage. Thus, further analysis and future observations are required to cast light on this topic. Additional tests in other theoretical frameworks (e.g., the models proposed by Beloborodov et al. 2014; Ito et al. 2014) will be the object of a forthcoming paper.

## 7. Summary and Conclusions

We have here presented the second Fermi-LAT GRB catalog. Compared to the 1FLGC, there are several changes and improvements to the analysis, beyond the three times longer timescale:

1. Pass 8 data have replaced Pass 6 data.

2. Over the past 2 yr, a new detection algorithm has been developed at energies above 100 MeV.

3. A BB algorithm has been employed to look for detections at 30–100 MeV.

4. All available triggers (both GBM and other space-based missions) were used in the search for LAT counterparts.

All potential detections have been processed by a dedicated and standardized analysis. This differs from both the analysis performed for the 1FLGC and the standard LAT real-time analysis. In total we find 186 detections, a much larger sample than the 35 GRBs of the 1FLGC.

The results presented in the 2FLGC have strengthened some of the conclusions from the 1FLGC, most notably the temporal characteristics. A major difference from the 1FLGC is a new definition of the duration, where we consider the first and last detected photon. Nevertheless, the properties of later onset and longer duration for the high-energy emission remain unchanged. Other trends suggested by the 1FLGC, most notably regarding the energetics and differences between long and short GRBs, have been shown not to hold with the larger sample in the 2FLGC.

While this catalog describes the temporal properties of the detections at both keV and MeV energy ranges, the spectral analysis is limited to the LAT range (100 MeV–100 GeV). This has shown to be problematic in some cases, such as when testing for breaks in the spectrum. These results are not surprising, as the division between the energy ranges is purely artificial. Nevertheless, a joint analysis is beyond the scope of this paper and will instead be the focus of a dedicated study presented in a separate publication.

The larger sample in this catalog is also reflected in the 34 GRBs with measured redshift. This has allowed us to better investigate the energetics also in the rest frame, probing both the highest-energy photons and the closure relations. The remarkable alignment seen in the decay of the luminosity further points to a strong correlation between isotropic
luminosity in the rest-frame MeV–GeV range and isotropic energy as measured at lower energies.

The 2FLGC presents the results of the first 10 yr of Fermi-LAT; it represents an initial step in understanding the processes behind high-energy emission from gamma-ray bursts. It is our hope that the Fermi data presented here will lead to a large number of follow-up analyses and continue to drive the evolution of theoretical models. Within the Fermi-LAT Collaboration, several such studies are already underway: focused studies on the closure relations, energetics, and the connection to VHE emission. The large sample of the 2FLGC will also allow more complete analysis of the detectability of MeV emission from GRBs, as well as detailed comparison between the energy output in the LAT range during the prompt versus afterglow phases. The collaboration is also in the process of setting up an online repository, similar to the analysis of the 2FLGC, which will automatically be updated with each new GRB detection. The past few years have truly shown the strength of Fermi observations, and its importance will continue for many years to come. This is especially true as the era of multimessenger studies of GRBs has just begun.

The Fermi-LAT Collaboration acknowledges generous ongoing support from a number of agencies and institutes that have supported both the development and the operation of the LAT, as well as scientific data analysis. These include the National Aeronautics and Space Administration and the Department of Energy in the United States; the Commissariat à l'Energie Atomique and the Centre National de la Recherche Scientifique/Institut National de Physique Nucléaire et de Physique des Particules in France; the Agenzia Spaziale Italiana and the Istituto Nazionale di Fisica Nucleare in Italy; the Ministry of Education, Culture, Sports, Science and Technology (MEXT), High Energy Accelerator Research Organization (KEK), and Japan Aerospace Exploration Agency (JAXA) in Japan; and the K. A. Wallenberg Foundation, the Swedish Research Council, and the Swedish National Space Agency in Sweden.

Additional support for science analysis during the operations phase is gratefully acknowledged from the Istituto Nazionale di Astrofisica in Italy and the Centre National d'Études Spatiales in France. This work performed in part under DOE Contract DE-AC02-76SF00515. M. Axelsson gratefully acknowledges funding from the European Union's Horizon 2020 research and innovation program under the Marie Skłodowska-Curie grant agreement No. 734303 (NEWS).

# Appendix A
## Tables

### A1. Sample of LAT Bursts

The complete list of the Fermi-LAT 10 yr catalog GRBs is presented in Table 2. The first column gives the GRB name as it appeared in the corresponding GCN. The following two columns display the GRB trigger date (in UTC) and time (in Mission Elapsed Time, or MET, in seconds from 2001 January 1 UT 00:00:00). Columns (4) through (6) report the best reconstructed position using only Fermi-LAT data. Coordinates are expressed in the equatorial frame of epoch J2000. The localization error (Column (6)) is the error obtained from the LTF analysis. The boresight angle $\theta$ and the zenith angle $\zeta$ at trigger time are given in Columns (7) and (8). Column (9) indicates the presence or

absence of an ARR. The detection significance of the LTF and LLE analysis are listed in Columns (10) and (11). The LTF TS value corresponds to the highest TS obtained during one of the four time windows listed in Table 1. LLE-only events, for which no position calculation was possible, report the best location available for that event. The reference to the burst location is then given in Column (12). When available, we report the redshift of the host galaxy in Column (13), together with its corresponding reference (Column (14)).

### A2. Comparison of Duration Estimation

Table 3 presents the different burst duration estimations. The GRB class (listed in Column (2); L = long, S = short) is derived from the GBM duration in the 50–300 keV energy range (long bursts have $T_{GBM,90} > 2$ s, and short bursts have $T_{GBM,90} < 2$ s). All GBM values ($T_{GBM,05}$, $T_{GBM,95}$, and $T_{GBM,90}$; Columns (3)–(5)) are taken from the FGGC. The LLE onset times and durations ($T_{LLE,05}$, $T_{LLE,95}$, and $T_{LLE,90}$), listed in Columns (6)–(8), are calculated in the 30 MeV–1 GeV energy range. The LAT onset time $T_{LAT,0}$, end time $T_{LAT,1}$, and duration $T_{LAT,100}$, calculated in the 100 MeV–100 GeV energy range, are given in Columns (9)–(11).

### A3. Likelihood Analysis in Different Time Window

Table 4 reports the likelihood analysis results in the different time windows. The definition of the time windows is explained in Table 1. The time range of each interval in which we perform the analysis is given in Column (2). Columns (3) and (4) display the number of selected events in the ROI and the number of events predicted by the best-fit model obtained after likelihood maximization, respectively. The TS is given in Column (5). Note that only intervals with $TS > 20$ are displayed in the table. Columns (6) through (9) contain the value of the spectral index, flux, and fluence in each interval and, when possible, the isotropic energy calculated in the 100 MeV–100 GeV rest-frame energy band.

### A4. Temporally Extended High-energy Emission

Table 5 displays the temporally extended high-energy emission analysis results. Columns (2) and (3) give the peak flux and peak flux time. These are followed by the fit parameters of the simple PL ($\alpha$; Column (4)) and BPL ($\alpha_1$, $\alpha_2$, and $T_b$; Columns (5)–(7)). Column (8) reports the best-fit decay index, indicating whether the SPL or the BPL is the preferred model. The last column indicates the total duration of the extended emission.

### A5. Highest-energy Events

For each GRB, Table 6 shows the highest-energy photon detected in the GBM time window (Columns (2)–(5)), as well as the highest-energy photon found overall in the time-resolved analysis (Columns (6)–(9)). For each time span the values given are the total number of photons found with probability >90% to be associated with the GRB, the energy of the most energetic photon, its arrival time, and the probability to be associated with the GRB.

Table 7 lists the 29 GRBs from which at least one photon with energy greater than 10 GeV has been detected. For each such photon, the energy and arrival time are given. For GRBs with a measured redshift, this is listed along with the corresponding source-frame energy ($E_{sf}$) of each photon.
<table><caption>Table 2
Sample of Fermi-LAT GRBs, from 2008 to 2018 August</caption><thead><tr><th>GRB Name</th><th>Date<br>(UT)</th><th>Trigger Time<br>(MET)</th><th>R.A.<br>(deg, J2000)</th><th>Decl.<br>(deg, J2000)</th><th>Loc. Err.<br>(deg)</th><th>θ<br>(deg)</th><th>ζ<br>(deg)</th><th>ARR</th><th>LTF TS</th><th>LLE Sig.<br>(σ)</th><th>Loc. GCN</th><th>Redshift</th><th>GCN</th></tr></thead><tbody><tr><td>080818B</td><td>2008 Aug 18 22:40:49</td><td>240792050.080</td><td>109.46</td><td>15.0</td><td>0.47</td><td>56.4</td><td>75.0</td><td>0</td><td>26</td><td>...</td><td></td><td>...</td><td></td></tr><tr><td>080825C</td><td>2008 Aug 25 14:13:48</td><td>241366429.110</td><td>234.83</td><td>-3.84</td><td>0.65</td><td>60.2</td><td>26.3</td><td>0</td><td>56</td><td>8</td><td></td><td>...</td><td></td></tr><tr><td>080916C</td><td>2008 Sep 16 00:12:45</td><td>243216766.620</td><td>119.89</td><td>-56.7</td><td>0.07</td><td>49.0</td><td>83.8</td><td>0</td><td>2200</td><td>46</td><td></td><td>4.35</td><td>Greiner et al. (2009)</td></tr><tr><td>081006</td><td>2008 Oct 6 14:29:34</td><td>244996175.170</td><td>135.0</td><td>-61.95</td><td>0.38</td><td>10.8</td><td>40.3</td><td>0</td><td>97</td><td>...</td><td></td><td>...</td><td></td></tr><tr><td>081009</td><td>2008 Oct 9 03:20:58</td><td>245215259.060</td><td>251.57</td><td>18.92</td><td>0.28</td><td>69.6</td><td>55.4</td><td>0</td><td>50</td><td>...</td><td></td><td>...</td><td></td></tr><tr><td>081024B</td><td>2008 Oct 24 21:22:40</td><td>246576161.860</td><td>323.01</td><td>20.84</td><td>0.31</td><td>19.2</td><td>30.6</td><td>0</td><td>83</td><td>8</td><td></td><td>...</td><td></td></tr><tr><td>081102B</td><td>2008 Nov 2 08:45:00</td><td>247308301.510</td><td>212.95</td><td>30.33</td><td>0.58</td><td>48.2</td><td>47.6</td><td>0</td><td>26</td><td>...</td><td></td><td>...</td><td></td></tr><tr><td>081122A</td><td>2008 Nov 22 12:28:12</td><td>249049693.210</td><td>338.82</td><td>37.57</td><td>0.27</td><td>21.5</td><td>25.9</td><td>1</td><td>38</td><td>...</td><td></td><td>...</td><td></td></tr><tr><td>081203A</td><td>2008 Dec 3 13:57:11</td><td>250005432.570</td><td>233.49</td><td>63.96</td><td>0.39</td><td>78.0</td><td>106.1</td><td>0</td><td>30</td><td>...</td><td></td><td>2.10</td><td>Landsman et al. (2008)</td></tr><tr><td>081224</td><td>2008 Dec 24 21:17:55</td><td>251846276.410</td><td>201.1</td><td>72.67</td><td>0.33</td><td>15.4</td><td>49.2</td><td>1</td><td>29</td><td>...</td><td></td><td>...</td><td></td></tr><tr><td>090102</td><td>2009 Jan 2 02:55:30</td><td>252557732.850</td><td>127.56</td><td>33.46</td><td>0.65</td><td>85.5</td><td>78.7</td><td>0</td><td>20</td><td>...</td><td></td><td>1.55</td><td>de Ugarte Postigo et al. (2009)</td></tr><tr><td>090217</td><td>2009 Feb 17 04:56:42</td><td>256539404.560</td><td>205.01</td><td>-8.39</td><td>0.46</td><td>33.7</td><td>46.2</td><td>0</td><td>100</td><td>16</td><td></td><td>...</td><td></td></tr><tr><td>090227A</td><td>2009 Feb 27 07:25:57</td><td>257412359.000</td><td>4.56</td><td>-44.63</td><td>0.34</td><td>21.6</td><td>45.4</td><td>0</td><td>33</td><td>6</td><td></td><td>...</td><td></td></tr><tr><td>090227B</td><td>2009 Feb 27 18:31:01</td><td>257452263.407</td><td>11.8</td><td>32.2</td><td>1.80</td><td>72.7</td><td>86.6</td><td>1</td><td>...</td><td>15</td><td>FGGC</td><td>...</td><td></td></tr><tr><td>090228A</td><td>2009 Feb 28 04:53:20</td><td>257489602.920</td><td>98.6</td><td>-28.79</td><td>0.21</td><td>23.8</td><td>56.6</td><td>1</td><td>41</td><td>...</td><td></td><td>...</td><td></td></tr><tr><td>090323</td><td>2009 Mar 23 00:02:42</td><td>259459364.630</td><td>190.75</td><td>17.07</td><td>0.09</td><td>56.5</td><td>66.2</td><td>1</td><td>230</td><td>18</td><td></td><td>3.57</td><td>Chornock et al. (2009)</td></tr><tr><td>090328</td><td>2009 Mar 28 09:36:46</td><td>259925808.510</td><td>90.67</td><td>-41.95</td><td>0.12</td><td>65.7</td><td>72.4</td><td>1</td><td>290</td><td>20</td><td></td><td>0.74</td><td>Cenko et al. (2009)</td></tr><tr><td>090427A</td><td>2009 Apr 27 23:26:32</td><td>262567594.000</td><td>235.86</td><td>-14.05</td><td>0.11</td><td>134.0</td><td>135.4</td><td>0</td><td>40</td><td>...</td><td></td><td>...</td><td></td></tr><tr><td>090510</td><td>2009 May 10 00:22:59</td><td>263607781.970</td><td>333.57</td><td>-26.62</td><td>0.05</td><td>13.5</td><td>52.3</td><td>1</td><td>2600</td><td>24</td><td></td><td>0.90</td><td>Rau et al. (2009)</td></tr><tr><td>090531B</td><td>2009 May 31 18:35:56</td><td>265487758.492</td><td>252.06</td><td>-36.05</td><td>...</td><td>21.9</td><td>49.2</td><td>0</td><td>...</td><td>10</td><td>FGGC</td><td>...</td><td></td></tr><tr><td>090626</td><td>2009 Jun 26 04:32:08</td><td>267683530.880</td><td>170.19</td><td>-33.56</td><td>0.29</td><td>22.6</td><td>10.4</td><td>0</td><td>100</td><td>9</td><td></td><td>...</td><td></td></tr><tr><td>090720B</td><td>2009 Jul 20 17:02:56</td><td>269802178.900</td><td>203.638</td><td>-51.198</td><td>0.99</td><td>56.3</td><td>94.2</td><td>1</td><td>35</td><td>...</td><td></td><td>...</td><td></td></tr><tr><td>090902B</td><td>2009 Sep 2 11:05:08</td><td>273582310.320</td><td>264.95</td><td>27.33</td><td>0.04</td><td>50.3</td><td>37.6</td><td>1</td><td>3600</td><td>8</td><td></td><td>1.82</td><td>Cucchiara et al. (2009c)</td></tr><tr><td>090926A</td><td>2009 Sep 26 04:20:26</td><td>275631628.990</td><td>353.56</td><td>-66.37</td><td>0.06</td><td>48.0</td><td>80.4</td><td>1</td><td>3200</td><td>29</td><td></td><td>2.11</td><td>Malesani et al. (2009)</td></tr><tr><td>091003</td><td>2009 Oct 3 04:35:45</td><td>276237347.580</td><td>251.54</td><td>36.62</td><td>0.10</td><td>12.2</td><td>62.1</td><td>1</td><td>170</td><td>...</td><td></td><td>0.90</td><td>Cucchiara et al. (2009b)</td></tr><tr><td>091031</td><td>2009 Oct 31 12:00:28</td><td>278683230.850</td><td>71.81</td><td>-57.73</td><td>0.23</td><td>24.6</td><td>39.4</td><td>0</td><td>110</td><td>17</td><td></td><td>...</td><td></td></tr><tr><td>091120</td><td>2009 Nov 20 04:34:40</td><td>280384482.230</td><td>226.7</td><td>-21.36</td><td>0.18</td><td>48.1</td><td>25.7</td><td>1</td><td>38</td><td>...</td><td></td><td>...</td><td></td></tr><tr><td>091127</td><td>2009 Nov 27 23:25:45</td><td>281057147.490</td><td>36.65</td><td>-19.14</td><td>0.31</td><td>25.8</td><td>54.5</td><td>1</td><td>34</td><td>...</td><td></td><td>0.49</td><td>Cucchiara et al. (2009a)</td></tr><tr><td>100116A</td><td>2010 Jan 16 21:31:00</td><td>285370262.240</td><td>305.0</td><td>14.51</td><td>0.06</td><td>25.9</td><td>36.3</td><td>0</td><td>130</td><td>31</td><td></td><td>...</td><td></td></tr><tr><td>100213C</td><td>2010 Feb 13 12:44:41</td><td>287757883.000</td><td>282.45</td><td>15.15</td><td>0.06</td><td>146.4</td><td>147.2</td><td>0</td><td>45</td><td>...</td><td></td><td>...</td><td></td></tr><tr><td>100225A</td><td>2010 Feb 25 02:45:31</td><td>288758733.150</td><td>311.93</td><td>-60.51</td><td>0.34</td><td>57.0</td><td>83.7</td><td>0</td><td>47</td><td>10</td><td></td><td>...</td><td></td></tr><tr><td>100325A</td><td>2010 Mar 25 06:36:08</td><td>291191770.030</td><td>330.26</td><td>-26.99</td><td>0.87</td><td>7.8</td><td>52.9</td><td>0</td><td>30</td><td>...</td><td></td><td>...</td><td></td></tr><tr><td>100414A</td><td>2010 Apr 14 02:20:21</td><td>292904423.990</td><td>192.15</td><td>8.56</td><td>0.09</td><td>69.6</td><td>29.0</td><td>1</td><td>190</td><td>...</td><td></td><td>1.37</td><td>Cucchiara & Fox (2010)</td></tr><tr><td>100423B</td><td>2010 Apr 23 05:51:25</td><td>293694687.750</td><td>120.18</td><td>13.36</td><td>0.19</td><td>47.9</td><td>32.8</td><td>0</td><td>31</td><td>...</td><td></td><td>...</td><td></td></tr><tr><td>100511A</td><td>2010 May 11 00:49:56</td><td>295231798.230</td><td>111.38</td><td>-0.99</td><td>0.07</td><td>48.0</td><td>26.2</td><td>1</td><td>93</td><td>...</td><td></td><td>...</td><td></td></tr><tr><td>100620A</td><td>2010 Jun 20 02:51:29</td><td>298695091.120</td><td>87.68</td><td>-50.66</td><td>1.24</td><td>25.2</td><td>73.3</td><td>0</td><td>29</td><td>...</td><td></td><td>...</td><td></td></tr><tr><td>100724B</td><td>2010 Jul 24 00:42:06</td><td>301624928.090</td><td>123.59</td><td>75.86</td><td>0.10</td><td>50.0</td><td>96.2</td><td>1</td><td>110</td><td>40</td><td></td><td>...</td><td></td></tr><tr><td>100728A</td><td>2010 Jul 28 02:17:30</td><td>301976252.610</td><td>88.79</td><td>-15.41</td><td>0.18</td><td>60.5</td><td>36.9</td><td>1</td><td>63</td><td>...</td><td></td><td>1.57</td><td>Kruehler et al. (2010)</td></tr><tr><td>100826A</td><td>2010 Aug 26 22:58:22</td><td>304556304.900</td><td>277.71</td><td>-25.27</td><td>1.13</td><td>73.0</td><td>85.1</td><td>1</td><td>29</td><td>21</td><td></td><td>...</td><td></td></tr><tr><td>101014A</td><td>2010 Oct 14 04:11:52</td><td>308722314.620</td><td>26.5</td><td>-48.54</td><td>0.09</td><td>55.0</td><td>102.6</td><td>1</td><td>67</td><td>12</td><td></td><td>...</td><td></td></tr><tr><td>101107A</td><td>2010 Nov 7 00:16:25</td><td>310781787.120</td><td>170.27</td><td>23.88</td><td>0.20</td><td>36.6</td><td>28.4</td><td>0</td><td>47</td><td>...</td><td></td><td>...</td><td></td></tr><tr><td>101123A</td><td>2010 Nov 23 22:51:34</td><td>312245496.973</td><td>135.16</td><td>1.91</td><td>1.00</td><td>77.4</td><td>97.8</td><td>0</td><td>...</td><td>8</td><td>FGGC</td><td>...</td><td></td></tr><tr><td>101227B</td><td>2010 Dec 27 09:45:06</td><td>315135908.570</td><td>240.21</td><td>-15.93</td><td>0.32</td><td>13.7</td><td>37.3</td><td>0</td><td>27</td><td>...</td><td></td><td>...</td><td></td></tr><tr><td>110120A</td><td>2011 Jan 20 15:59:39</td><td>317231981.230</td><td>61.62</td><td>-12.02</td><td>0.20</td><td>14.6</td><td>31.5</td><td>1</td><td>67</td><td>6</td><td></td><td>...</td><td></td></tr><tr><td>110123A</td><td>2011 Jan 23 19:17:45</td><td>317503067.050</td><td>240.61</td><td>28.37</td><td>0.22</td><td>55.4</td><td>91.9</td><td>1</td><td>62</td><td>...</td><td></td><td>...</td><td></td></tr></tbody></table>
<table><caption>Table 2
(Continued)</caption><thead><tr><th>GRB Name</th><th>Date<br>(UT)</th><th>Trigger Time<br>(MET)</th><th>R.A.<br>(deg, J2000)</th><th>Decl.<br>(deg, J2000)</th><th>Loc. Err.<br>(deg)</th><th>θ<br>(deg)</th><th>ζ<br>(deg)</th><th>ARR</th><th>LTF TS</th><th>LLE Sig.<br>(σ)</th><th>Loc. GCN</th><th>Redshift</th><th>GCN</th></tr></thead><tbody><tr><td>110213A</td><td>2011 Feb 13 05:17:11</td><td>319267033.280</td><td>42.98</td><td>49.11</td><td>0.22</td><td>98.6</td><td>53.2</td><td>0</td><td>28</td><td>...</td><td></td><td>1.46</td><td>Milne & Cenko (2011)</td></tr><tr><td>110328B</td><td>2011 Mar 28 12:29:19</td><td>323008161.190</td><td>121.66</td><td>41.52</td><td>0.16</td><td>36.0</td><td>43.3</td><td>0</td><td>50</td><td>29</td><td></td><td>...</td><td></td></tr><tr><td>110428A</td><td>2011 Apr 28 09:18:30</td><td>325675112.400</td><td>5.42</td><td>65.01</td><td>0.18</td><td>34.6</td><td>84.3</td><td>0</td><td>52</td><td>...</td><td></td><td>...</td><td></td></tr><tr><td>110518A</td><td>2011 May 18 20:38:10</td><td>327443892.000</td><td>67.19</td><td>-34.27</td><td>0.19</td><td>165.9</td><td>124.1</td><td>0</td><td>40</td><td>...</td><td></td><td>...</td><td></td></tr><tr><td>110529A</td><td>2011 May 29 00:48:42</td><td>328322924.870</td><td>118.33</td><td>67.91</td><td>1.50</td><td>30.0</td><td>79.9</td><td>1</td><td>...</td><td>15</td><td>FGGC</td><td>...</td><td></td></tr><tr><td>110625A</td><td>2011 Jun 25 21:08:18</td><td>330728900.230</td><td>286.38</td><td>6.4</td><td>0.13</td><td>87.9</td><td>51.4</td><td>1</td><td>55</td><td>...</td><td></td><td>...</td><td></td></tr><tr><td>110721A</td><td>2011 Jul 21 04:47:43</td><td>332916465.760</td><td>333.53</td><td>-38.58</td><td>0.12</td><td>40.2</td><td>82.9</td><td>1</td><td>220</td><td>61</td><td></td><td>...</td><td></td></tr><tr><td>110728A</td><td>2011 Jul 28 01:20:22</td><td>333508824.810</td><td>173.57</td><td>4.34</td><td>0.62</td><td>41.7</td><td>16.8</td><td>0</td><td>43</td><td>...</td><td></td><td>...</td><td></td></tr><tr><td>110731A</td><td>2011 Jul 31 11:09:29</td><td>333803371.950</td><td>280.42</td><td>-28.65</td><td>0.24</td><td>2.9</td><td>50.7</td><td>1</td><td>390</td><td>20</td><td></td><td>2.83</td><td>Tanvir et al. (2011)</td></tr><tr><td>110903A</td><td>2011 Sep 3 02:39:34</td><td>336710376.420</td><td>197.26</td><td>58.85</td><td>0.07</td><td>44.3</td><td>93.9</td><td>0</td><td>49</td><td>...</td><td></td><td>...</td><td></td></tr><tr><td>110921B</td><td>2011 Sep 21 21:52:45</td><td>338334767.090</td><td>5.77</td><td>-26.59</td><td>0.39</td><td>33.9</td><td>21.2</td><td>1</td><td>50</td><td>7</td><td></td><td>...</td><td></td></tr><tr><td>111210B</td><td>2011 Dec 10 08:37:39</td><td>345199061.000</td><td>122.51</td><td>44.09</td><td>0.82</td><td>28.3</td><td>24.1</td><td>0</td><td>34</td><td>...</td><td></td><td>...</td><td></td></tr><tr><td>120107A</td><td>2012 Jan 7 09:12:15</td><td>347620337.410</td><td>245.01</td><td>-70.59</td><td>0.34</td><td>46.8</td><td>96.7</td><td>0</td><td>30</td><td>...</td><td></td><td>...</td><td></td></tr><tr><td>120226A</td><td>2012 Feb 26 20:54:17</td><td>351982459.030</td><td>302.97</td><td>48.98</td><td>0.68</td><td>37.1</td><td>42.2</td><td>1</td><td>37</td><td>12</td><td></td><td>...</td><td></td></tr><tr><td>120316A</td><td>2012 Mar 16 00:11:02</td><td>353549464.560</td><td>57.35</td><td>-56.41</td><td>0.31</td><td>9.9</td><td>48.7</td><td>0</td><td>57</td><td>5</td><td></td><td>...</td><td></td></tr><tr><td>120328B</td><td>2012 Mar 28 06:26:20</td><td>354608782.953</td><td>228.14</td><td>22.8</td><td>1.00</td><td>68.5</td><td>106.6</td><td>0</td><td>...</td><td>8</td><td>FGGC</td><td>...</td><td></td></tr><tr><td>120420B</td><td>2012 Apr 20 20:35:13</td><td>356646915.070</td><td>108.77</td><td>15.17</td><td>0.29</td><td>109.7</td><td>76.8</td><td>0</td><td>29</td><td>...</td><td></td><td>...</td><td></td></tr><tr><td>120526A</td><td>2012 May 26 07:16:40</td><td>359709402.760</td><td>53.82</td><td>-35.23</td><td>0.11</td><td>109.7</td><td>62.9</td><td>0</td><td>110</td><td>...</td><td></td><td>...</td><td></td></tr><tr><td>120624B</td><td>2012 Jun 24 22:23:54</td><td>362269436.930</td><td>170.91</td><td>8.94</td><td>0.21</td><td>68.2</td><td>26.6</td><td>1</td><td>450</td><td>9</td><td></td><td>2.20</td><td>de Ugarte Postigo et al. (2013a)</td></tr><tr><td>120709A</td><td>2012 Jul 9 21:11:40</td><td>363561103.370</td><td>318.8</td><td>-50.42</td><td>0.40</td><td>27.8</td><td>40.5</td><td>0</td><td>84</td><td>9</td><td></td><td>...</td><td></td></tr><tr><td>120711A</td><td>2012 Jul 11 02:44:53</td><td>363667496.290</td><td>94.64</td><td>-70.88</td><td>0.13</td><td>134.5</td><td>85.0</td><td>1</td><td>140</td><td>...</td><td></td><td>1.41</td><td>Tanvir et al. (2012)</td></tr><tr><td>120729A</td><td>2012 Jul 29 10:56:12</td><td>365252175.660</td><td>13.26</td><td>49.91</td><td>0.38</td><td>82.6</td><td>32.7</td><td>0</td><td>35</td><td>...</td><td></td><td>0.80</td><td>Tanvir & Ball (2012)</td></tr><tr><td>120830A</td><td>2012 Aug 30 07:07:03</td><td>368003226.540</td><td>88.59</td><td>-28.79</td><td>0.62</td><td>38.7</td><td>75.3</td><td>0</td><td>44</td><td>...</td><td></td><td>...</td><td></td></tr><tr><td>120911B</td><td>2012 Sep 11 06:26:12</td><td>369037575.000</td><td>172.66</td><td>-37.5</td><td>0.60</td><td>66.1</td><td>22.2</td><td>0</td><td>340</td><td>18</td><td></td><td>...</td><td></td></tr><tr><td>120915A</td><td>2012 Sep 15 00:00:41</td><td>369360044.640</td><td>240.95</td><td>57.04</td><td>1.76</td><td>10.8</td><td>59.2</td><td>0</td><td>32</td><td>...</td><td></td><td>...</td><td></td></tr><tr><td>120919B</td><td>2012 Sep 19 01:14:23</td><td>369710066.070</td><td>302.73</td><td>-37.73</td><td>0.14</td><td>22.9</td><td>65.9</td><td>0</td><td>30</td><td>...</td><td></td><td>...</td><td></td></tr><tr><td>121011A</td><td>2012 Oct 11 11:15:25</td><td>371646928.697</td><td>260.2</td><td>41.14</td><td>0.01</td><td>55.9</td><td>84.8</td><td>0</td><td>...</td><td>8</td><td>GCN13845</td><td>...</td><td></td></tr><tr><td>121029A</td><td>2012 Oct 29 08:24:19</td><td>373191862.930</td><td>229.07</td><td>-25.56</td><td>0.21</td><td>61.2</td><td>42.0</td><td>1</td><td>40</td><td>...</td><td></td><td>...</td><td></td></tr><tr><td>121123B</td><td>2012 Nov 23 10:35:55</td><td>375359758.710</td><td>28.62</td><td>-24.06</td><td>0.50</td><td>78.8</td><td>76.4</td><td>0</td><td>27</td><td>...</td><td></td><td>...</td><td></td></tr><tr><td>121225B</td><td>2012 Dec 25 10:00:53</td><td>378122456.583</td><td>310.45</td><td>-34.83</td><td>1.50</td><td>64.3</td><td>90.7</td><td>0</td><td>...</td><td>8</td><td>FGGC</td><td>...</td><td></td></tr><tr><td>130305A</td><td>2013 Mar 5 11:39:11</td><td>384176354.369</td><td>116.765</td><td>52.025</td><td>0.06</td><td>41.5</td><td>99.6</td><td>0</td><td>...</td><td>19</td><td>GCN14257</td><td>...</td><td></td></tr><tr><td>130310A</td><td>2013 Mar 10 20:09:41</td><td>384638984.500</td><td>142.38</td><td>-17.1</td><td>0.47</td><td>75.6</td><td>49.8</td><td>1</td><td>47</td><td>11</td><td></td><td>...</td><td></td></tr><tr><td>130325A</td><td>2013 Mar 25 04:51:54</td><td>385879917.300</td><td>122.97</td><td>-19.15</td><td>0.18</td><td>47.3</td><td>76.1</td><td>1</td><td>60</td><td>...</td><td></td><td>...</td><td></td></tr><tr><td>130327B</td><td>2013 Mar 27 08:24:04</td><td>386065447.050</td><td>218.26</td><td>-69.76</td><td>0.21</td><td>46.3</td><td>86.2</td><td>1</td><td>250</td><td>...</td><td></td><td>...</td><td></td></tr><tr><td>130427A</td><td>2013 Apr 27 07:47:06</td><td>388741629.420</td><td>173.14</td><td>27.73</td><td>0.04</td><td>46.4</td><td>62.8</td><td>1</td><td>4200</td><td>15</td><td></td><td>0.34</td><td>Xu et al. (2013)</td></tr><tr><td>130502B</td><td>2013 May 2 07:51:11</td><td>389173874.760</td><td>66.65</td><td>70.83</td><td>0.32</td><td>46.4</td><td>96.4</td><td>0</td><td>330</td><td>8</td><td></td><td>...</td><td></td></tr><tr><td>130504C</td><td>2013 May 4 23:28:57</td><td>389402940.520</td><td>91.55</td><td>3.78</td><td>0.12</td><td>47.6</td><td>52.3</td><td>1</td><td>91</td><td>29</td><td></td><td>...</td><td></td></tr><tr><td>130518A</td><td>2013 May 18 13:54:37</td><td>390578080.530</td><td>355.59</td><td>47.53</td><td>0.19</td><td>40.3</td><td>82.4</td><td>1</td><td>140</td><td>19</td><td></td><td>2.49</td><td>Sanchez-Ramirez et al. (2013)</td></tr><tr><td>130606B</td><td>2013 Jun 6 11:55:33</td><td>392212536.630</td><td>218.57</td><td>-22.15</td><td>0.11</td><td>90.5</td><td>10.9</td><td>1</td><td>87</td><td>...</td><td></td><td>...</td><td></td></tr><tr><td>130702A</td><td>2013 Jul 2 00:05:23</td><td>394416326.080</td><td>216.39</td><td>16.15</td><td>0.62</td><td>76.7</td><td>78.9</td><td>0</td><td>26</td><td>...</td><td></td><td>0.14</td><td>Leloudas et al. (2013)</td></tr><tr><td>130821A</td><td>2013 Aug 21 16:10:28</td><td>398794231.010</td><td>314.23</td><td>-11.52</td><td>0.18</td><td>34.9</td><td>49.2</td><td>1</td><td>170</td><td>16</td><td></td><td>...</td><td></td></tr><tr><td>130828A</td><td>2013 Aug 28 07:20:00</td><td>399367203.150</td><td>259.72</td><td>28.21</td><td>0.29</td><td>39.9</td><td>45.5</td><td>1</td><td>150</td><td>9</td><td></td><td>...</td><td></td></tr><tr><td>130907A</td><td>2013 Sep 7 21:41:13</td><td>400282876.090</td><td>215.92</td><td>45.7</td><td>0.26</td><td>105.4</td><td>57.2</td><td>0</td><td>30</td><td>...</td><td></td><td>1.24</td><td>de Ugarte Postigo et al. (2013c)</td></tr></tbody></table>
<table><caption>Table 2
(Continued)</caption><thead><tr><th>GRB Name</th><th>Date<br>(UT)</th><th>Trigger Time<br>(MET)</th><th>R.A.<br>(deg, J2000)</th><th>Decl.<br>(deg, J2000)</th><th>Loc. Err.<br>(deg)</th><th>θ<br>(deg)</th><th>ζ<br>(deg)</th><th>ARR</th><th>LTF TS</th><th>LLE Sig.<br>(σ)</th><th>Loc. GCN</th><th>Redshift</th><th>GCN</th></tr></thead><tbody><tr><td>131014A</td><td>2013 Oct 14 05:09:00</td><td>403420143.200</td><td>100.78</td><td>-20.76</td><td>0.81</td><td>72.5</td><td>76.0</td><td>1</td><td>98</td><td>17</td><td></td><td>...</td><td></td></tr><tr><td>131029A</td><td>2013 Oct 29 23:20:48</td><td>404781651.580</td><td>200.85</td><td>48.03</td><td>0.18</td><td>59.8</td><td>94.1</td><td>0</td><td>110</td><td>...</td><td></td><td>...</td><td></td></tr><tr><td>131108A</td><td>2013 Nov 8 20:41:55</td><td>405636118.760</td><td>156.81</td><td>9.57</td><td>0.22</td><td>20.4</td><td>21.9</td><td>1</td><td>1200</td><td>37</td><td></td><td>2.40</td><td>de Ugarte Postigo et al.<br>(2013b)</td></tr><tr><td>131209A</td><td>2013 Dec 9 13:07:56</td><td>408287279.960</td><td>137.06</td><td>-33.9</td><td>1.16</td><td>19.8</td><td>53.6</td><td>0</td><td>29</td><td>...</td><td></td><td>...</td><td></td></tr><tr><td>131216A</td><td>2013 Dec 16 01:56:32</td><td>408851795.063</td><td>91.59</td><td>-35.5</td><td>2.24</td><td>44.6</td><td>77.2</td><td>1</td><td>...</td><td>7</td><td>FGGC</td><td>...</td><td></td></tr><tr><td>131231A</td><td>2013 Dec 31 04:45:16</td><td>410157919.080</td><td>10.55</td><td>-1.72</td><td>0.05</td><td>37.2</td><td>39.9</td><td>1</td><td>240</td><td>18</td><td></td><td>0.64</td><td>Xu et al. (2014b)</td></tr><tr><td>140102A</td><td>2014 Jan 2 21:17:37</td><td>410390260.810</td><td>212.02</td><td>1.47</td><td>0.24</td><td>47.2</td><td>44.7</td><td>1</td><td>55</td><td>11</td><td></td><td>...</td><td></td></tr><tr><td>140104B</td><td>2014 Jan 4 17:32:03</td><td>410549526.150</td><td>218.76</td><td>-8.98</td><td>0.18</td><td>17.6</td><td>36.5</td><td>0</td><td>62</td><td>...</td><td></td><td>...</td><td></td></tr><tr><td>140110A</td><td>2014 Jan 10 06:18:37</td><td>411027520.940</td><td>29.12</td><td>-36.5</td><td>0.25</td><td>36.1</td><td>64.6</td><td>0</td><td>200</td><td>13</td><td></td><td>...</td><td></td></tr><tr><td>140124A</td><td>2014 Jan 24 12:38:31</td><td>412259914.190</td><td>59.27</td><td>36.86</td><td>0.33</td><td>19.4</td><td>65.6</td><td>0</td><td>26</td><td>...</td><td></td><td>...</td><td></td></tr><tr><td>140206B</td><td>2014 Feb 6 06:36:12</td><td>413361375.840</td><td>315.24</td><td>-8.79</td><td>0.08</td><td>46.2</td><td>86.0</td><td>1</td><td>300</td><td>39</td><td></td><td>...</td><td></td></tr><tr><td>140219A</td><td>2014 Feb 19 19:46:32</td><td>414531995.240</td><td>158.2</td><td>7.26</td><td>0.85</td><td>96.2</td><td>113.4</td><td>0</td><td>26</td><td>...</td><td></td><td>...</td><td></td></tr><tr><td>140323A</td><td>2014 Mar 23 10:22:53</td><td>417262976.120</td><td>356.52</td><td>-79.92</td><td>0.18</td><td>30.3</td><td>68.4</td><td>0</td><td>40</td><td>...</td><td></td><td>...</td><td></td></tr><tr><td>140402A</td><td>2014 Apr 2 00:10:06</td><td>418090209.990</td><td>207.66</td><td>5.97</td><td>0.16</td><td>13.7</td><td>46.0</td><td>0</td><td>93</td><td>8</td><td></td><td>...</td><td></td></tr><tr><td>140416A</td><td>2014 Apr 16 01:26:40</td><td>419304403.930</td><td>35.74</td><td>42.1</td><td>0.10</td><td>138.8</td><td>94.8</td><td>0</td><td>34</td><td>...</td><td></td><td>...</td><td></td></tr><tr><td>140523A</td><td>2014 May 23 03:05:57</td><td>422507160.630</td><td>133.09</td><td>24.83</td><td>0.38</td><td>55.4</td><td>26.4</td><td>0</td><td>180</td><td>...</td><td></td><td>...</td><td></td></tr><tr><td>140528A</td><td>2014 May 28 20:05:22</td><td>423000325.750</td><td>280.59</td><td>-59.46</td><td>0.36</td><td>46.2</td><td>66.3</td><td>1</td><td>31</td><td>...</td><td></td><td>...</td><td></td></tr><tr><td>140619B</td><td>2014 Jun 19 11:24:40</td><td>424869883.520</td><td>132.66</td><td>-9.65</td><td>0.20</td><td>38.9</td><td>25.8</td><td>0</td><td>210</td><td>10</td><td></td><td>...</td><td></td></tr><tr><td>140723A</td><td>2014 Jul 23 01:36:30</td><td>427772193.730</td><td>211.05</td><td>-3.57</td><td>0.51</td><td>54.4</td><td>38.6</td><td>0</td><td>110</td><td>8</td><td></td><td>...</td><td></td></tr><tr><td>140729A</td><td>2014 Jul 29 00:36:53</td><td>428287016.710</td><td>193.9</td><td>15.42</td><td>0.26</td><td>26.5</td><td>40.9</td><td>0</td><td>69</td><td>...</td><td></td><td>...</td><td></td></tr><tr><td>140810A</td><td>2014 Aug 10 18:46:10</td><td>429389173.090</td><td>118.93</td><td>27.4</td><td>0.07</td><td>123.5</td><td>89.3</td><td>1</td><td>79</td><td>...</td><td></td><td>...</td><td></td></tr><tr><td>140825A</td><td>2014 Aug 25 06:55:25</td><td>430642528.000</td><td>89.01</td><td>-12.16</td><td>0.57</td><td>113.4</td><td>112.0</td><td>0</td><td>19</td><td>...</td><td></td><td>...</td><td></td></tr><tr><td>140928A</td><td>2014 Sep 28 10:29:53</td><td>433592996.550</td><td>43.76</td><td>-55.92</td><td>0.10</td><td>110.0</td><td>111.3</td><td>0</td><td>57</td><td>...</td><td></td><td>...</td><td></td></tr><tr><td>141012A</td><td>2014 Oct 12 18:33:17</td><td>434831600.670</td><td>289.92</td><td>-53.0</td><td>0.27</td><td>47.2</td><td>94.7</td><td>0</td><td>110</td><td>...</td><td></td><td>...</td><td></td></tr><tr><td>141028A</td><td>2014 Oct 28 10:54:46</td><td>436186489.780</td><td>322.78</td><td>-0.64</td><td>0.27</td><td>36.4</td><td>48.1</td><td>1</td><td>130</td><td>31</td><td></td><td>2.33</td><td>Xu et al. (2014a)</td></tr><tr><td>141102A</td><td>2014 Nov 2 12:51:39</td><td>436625502.260</td><td>209.07</td><td>-46.73</td><td>0.68</td><td>52.3</td><td>53.5</td><td>1</td><td>34</td><td>...</td><td></td><td>...</td><td></td></tr><tr><td>141113A</td><td>2014 Nov 13 08:17:43</td><td>437559466.500</td><td>182.32</td><td>77.38</td><td>0.71</td><td>53.7</td><td>86.4</td><td>0</td><td>37</td><td>...</td><td></td><td>...</td><td></td></tr><tr><td>141207A</td><td>2014 Dec 7 19:11:21</td><td>439672284.100</td><td>159.8</td><td>3.95</td><td>0.20</td><td>61.1</td><td>36.8</td><td>0</td><td>260</td><td>15</td><td></td><td>...</td><td></td></tr><tr><td>141221B</td><td>2014 Dec 21 21:31:48</td><td>440890311.780</td><td>99.12</td><td>-66.9</td><td>0.33</td><td>38.7</td><td>84.3</td><td>0</td><td>31</td><td>...</td><td></td><td>...</td><td></td></tr><tr><td>141222A</td><td>2014 Dec 22 07:08:57</td><td>440924940.390</td><td>178.0</td><td>-57.31</td><td>0.08</td><td>45.9</td><td>93.5</td><td>0</td><td>150</td><td>8</td><td></td><td>...</td><td></td></tr><tr><td>150118B</td><td>2015 Jan 18 09:48:17</td><td>443267300.770</td><td>240.23</td><td>-35.93</td><td>0.56</td><td>65.3</td><td>91.7</td><td>1</td><td>28</td><td>6</td><td></td><td>...</td><td></td></tr><tr><td>150202B</td><td>2015 Feb 2 23:59:08</td><td>444614351.270</td><td>86.27</td><td>59.05</td><td>0.33</td><td>17.2</td><td>66.9</td><td>1</td><td>26</td><td>20</td><td></td><td>...</td><td></td></tr><tr><td>150210A</td><td>2015 Feb 10 22:26:24</td><td>445299987.290</td><td>112.28</td><td>13.57</td><td>0.43</td><td>53.7</td><td>10.0</td><td>0</td><td>92</td><td>14</td><td></td><td>...</td><td></td></tr><tr><td>150314A</td><td>2015 Mar 14 04:54:50</td><td>448001693.860</td><td>126.12</td><td>63.79</td><td>0.72</td><td>45.1</td><td>91.9</td><td>1</td><td>30</td><td>6</td><td></td><td>1.76</td><td>de Ugarte Postigo et al.<br>(2015a)</td></tr><tr><td>150403A</td><td>2015 Apr 3 21:54:10</td><td>449790853.960</td><td>311.75</td><td>-62.46</td><td>0.23</td><td>55.3</td><td>105.1</td><td>1</td><td>46</td><td>17</td><td></td><td>2.06</td><td>Pugliese et al. (2015)</td></tr><tr><td>150416A</td><td>2015 Apr 16 18:33:25</td><td>450902008.965</td><td>58.75</td><td>52.96</td><td>1.93</td><td>68.6</td><td>90.7</td><td>0</td><td>...</td><td>9</td><td>FGGC</td><td>...</td><td></td></tr><tr><td>150510A</td><td>2015 May 10 03:19:53</td><td>452920796.740</td><td>16.42</td><td>4.91</td><td>0.33</td><td>56.2</td><td>11.4</td><td>1</td><td>53</td><td>12</td><td></td><td>...</td><td></td></tr><tr><td>150513A</td><td>2015 May 13 20:33:15</td><td>453241998.250</td><td>49.06</td><td>-22.91</td><td>0.04</td><td>51.8</td><td>16.5</td><td>0</td><td>...</td><td>8</td><td>GCN17810</td><td>...</td><td></td></tr><tr><td>150514A</td><td>2015 May 14 18:35:05</td><td>453321308.350</td><td>74.92</td><td>-60.86</td><td>0.15</td><td>39.2</td><td>59.1</td><td>0</td><td>31</td><td>...</td><td></td><td>0.81</td><td>de Ugarte Postigo et al.<br>(2015b)</td></tr><tr><td>150523A</td><td>2015 May 23 09:29:48</td><td>454066191.080</td><td>115.34</td><td>-45.43</td><td>0.08</td><td>24.8</td><td>54.7</td><td>1</td><td>220</td><td>8</td><td></td><td>...</td><td></td></tr><tr><td>150627A</td><td>2015 Jun 27 04:23:23</td><td>457071806.680</td><td>117.5</td><td>-51.55</td><td>0.06</td><td>74.7</td><td>97.9</td><td>1</td><td>240</td><td>9</td><td></td><td>...</td><td></td></tr><tr><td>150702A</td><td>2015 Jul 2 23:56:38</td><td>457574202.550</td><td>52.59</td><td>-56.88</td><td>0.54</td><td>117.5</td><td>88.1</td><td>0</td><td>33</td><td>...</td><td></td><td>...</td><td></td></tr><tr><td>150820A</td><td>2015 Aug 20 21:07:15</td><td>461797639.642</td><td>258.65</td><td>-48.33</td><td>15.30</td><td>23.0</td><td>66.6</td><td>0</td><td>...</td><td>6</td><td>FGGC</td><td>...</td><td></td></tr></tbody></table>