Hostname: page-component-8448b6f56d-wq2xx Total loading time: 0 Render date: 2024-04-19T05:30:39.453Z Has data issue: false hasContentIssue false

DESAlert: Enabling Real-Time Transient Follow-Up with Dark Energy Survey Data

Published online by Cambridge University Press:  30 September 2016

A. Poci
Affiliation:
Department of Physics and Astronomy, Macquarie University, NSW 2109, Australia
K. Kuehn*
Affiliation:
Australian Astronomical Observatory, North Ryde, NSW 2113, Australia
T. Abbott
Affiliation:
Cerro Tololo Inter-American Observatory, National Optical Astronomy Observatory, Casilla 603, La Serena, Chile
F. B. Abdalla
Affiliation:
Department of Physics & Astronomy, University College London, Gower Street, London, WC1E 6BT, UK
S. Allam
Affiliation:
Fermi National Accelerator Laboratory, P. O. Box 500, Batavia, IL 60510, USA
A.H. Bauer
Affiliation:
Institut de Ciències de l’Espai, IEEC-CSIC, Campus UAB, Facultat de Ciències, Torre C5 par-2, 08193 Bellaterra, Barcelona, Spain
A. Benoit-Lévy
Affiliation:
Department of Physics & Astronomy, University College London, Gower Street, London, WC1E 6BT, UK
E. Bertin
Affiliation:
Institut d’Astrophysique de Paris, Univ. Pierre et Marie Curie & CNRS UMR7095, F-75014 Paris, France
D. Brooks
Affiliation:
Department of Physics & Astronomy, University College London, Gower Street, London, WC1E 6BT, UK
P. J. Brown
Affiliation:
George P. and Cynthia Woods Mitchell Institute for Fundamental Physics and Astronomy, and Department of Physics and Astronomy, Texas A&M University, College Station, TX 77843, USA
E. Buckley-Geer
Affiliation:
Fermi National Accelerator Laboratory, P. O. Box 500, Batavia, IL 60510, USA
D. L. Burke
Affiliation:
Kavli Institute for Particle Astrophysics & Cosmology, P. O. Box 2450, Stanford University, Stanford, CA 94305, USA SLAC National Accelerator Laboratory, Menlo Park, CA 94025, USA
A. Carnero Rosell
Affiliation:
Laboratório Interinstitucional de e-Astronomia - LIneA, Rua Gal. José Cristino 77, Rio de Janeiro, RJ - 20921-400, Brazil Observatório Nacional, Rua Gal. José Cristino 77, Rio de Janeiro, RJ - 20921-400, Brazil
M. Carrasco Kind
Affiliation:
Department of Astronomy, University of Illinois, 1002 W. Green Street, Urbana, IL 61801, USA National Center for Supercomputing Applications, 1205 West Clark St., Urbana, IL 61801, USA
R. Covarrubias
Affiliation:
National Center for Supercomputing Applications, 1205 West Clark St., Urbana, IL 61801, USA
L. N. da Costa
Affiliation:
Laboratório Interinstitucional de e-Astronomia - LIneA, Rua Gal. José Cristino 77, Rio de Janeiro, RJ - 20921-400, Brazil Observatório Nacional, Rua Gal. José Cristino 77, Rio de Janeiro, RJ - 20921-400, Brazil
C. B. D’Andrea
Affiliation:
Institute of Cosmology & Gravitation, University of Portsmouth, Portsmouth, PO1 3FX, UK
D. L. DePoy
Affiliation:
George P. and Cynthia Woods Mitchell Institute for Fundamental Physics and Astronomy, and Department of Physics and Astronomy, Texas A&M University, College Station, TX 77843, USA
S. Desai
Affiliation:
Department of Physics, Ludwig-Maximilians-Universität, Scheinerstr. 1, 81679 München, Germany
J. P. Dietrich
Affiliation:
Department of Physics, Ludwig-Maximilians-Universität, Scheinerstr. 1, 81679 München, Germany Excellence Cluster Universe, Boltzmannstr. 2, 85748 Garching, Germany
C. E Cunha
Affiliation:
Kavli Institute for Particle Astrophysics & Cosmology, P. O. Box 2450, Stanford University, Stanford, CA 94305, USA
T. F. Eifler
Affiliation:
Department of Physics and Astronomy, University of Pennsylvania, Philadelphia, PA 19104, USA Jet Propulsion Laboratory, California Institute of Technology, 4800 Oak Grove Dr., Pasadena, CA 91109, USA
J. Estrada
Affiliation:
Fermi National Accelerator Laboratory, P. O. Box 500, Batavia, IL 60510, USA
A. E. Evrard
Affiliation:
Department of Physics, University of Michigan, Ann Arbor, MI 48109, USA
A. Fausti Neto
Affiliation:
Laboratório Interinstitucional de e-Astronomia - LIneA, Rua Gal. José Cristino 77, Rio de Janeiro, RJ - 20921-400, Brazil
D. A. Finley
Affiliation:
Fermi National Accelerator Laboratory, P. O. Box 500, Batavia, IL 60510, USA
B. Flaugher
Affiliation:
Fermi National Accelerator Laboratory, P. O. Box 500, Batavia, IL 60510, USA
P. Fosalba
Affiliation:
Institut de Ciències de l’Espai, IEEC-CSIC, Campus UAB, Facultat de Ciències, Torre C5 par-2, 08193 Bellaterra, Barcelona, Spain
J. Frieman
Affiliation:
Fermi National Accelerator Laboratory, P. O. Box 500, Batavia, IL 60510, USA Kavli Institute for Cosmological Physics, University of Chicago, Chicago, IL 60637, USA
D. Gerdes
Affiliation:
Department of Physics, University of Michigan, Ann Arbor, MI 48109, USA
D. Gruen
Affiliation:
Max Planck Institute for Extraterrestrial Physics, Giessenbachstrasse, 85748 Garching, Germany University Observatory Munich, Scheinerstrasse 1, 81679 Munich, Germany
R. A. Gruendl
Affiliation:
Department of Astronomy, University of Illinois, 1002 W. Green Street, Urbana, IL 61801, USA National Center for Supercomputing Applications, 1205 West Clark St., Urbana, IL 61801, USA
K. Honscheid
Affiliation:
Center for Cosmology and Astro-Particle Physics, The Ohio State University, Columbus, OH 43210, USA Department of Physics, The Ohio State University, Columbus, OH 43210, USA
D. James
Affiliation:
Cerro Tololo Inter-American Observatory, National Optical Astronomy Observatory, Casilla 603, La Serena, Chile
N. Kuropatkin
Affiliation:
Fermi National Accelerator Laboratory, P. O. Box 500, Batavia, IL 60510, USA
O. Lahav
Affiliation:
Department of Physics & Astronomy, University College London, Gower Street, London, WC1E 6BT, UK
T. S. Li
Affiliation:
George P. and Cynthia Woods Mitchell Institute for Fundamental Physics and Astronomy, and Department of Physics and Astronomy, Texas A&M University, College Station, TX 77843, USA
M. March
Affiliation:
Department of Physics and Astronomy, University of Pennsylvania, Philadelphia, PA 19104, USA
J. Marshall
Affiliation:
George P. and Cynthia Woods Mitchell Institute for Fundamental Physics and Astronomy, and Department of Physics and Astronomy, Texas A&M University, College Station, TX 77843, USA
K. W. Merritt
Affiliation:
Fermi National Accelerator Laboratory, P. O. Box 500, Batavia, IL 60510, USA
C.J. Miller
Affiliation:
Department of Physics, University of Michigan, Ann Arbor, MI 48109, USA Department of Astronomy, University of Michigan, Ann Arbor, MI 48109, USA
R. C. Nichol
Affiliation:
Institute of Cosmology & Gravitation, University of Portsmouth, Portsmouth, PO1 3FX, UK
B. Nord
Affiliation:
Fermi National Accelerator Laboratory, P. O. Box 500, Batavia, IL 60510, USA
R. Ogando
Affiliation:
Laboratório Interinstitucional de e-Astronomia - LIneA, Rua Gal. José Cristino 77, Rio de Janeiro, RJ - 20921-400, Brazil Observatório Nacional, Rua Gal. José Cristino 77, Rio de Janeiro, RJ - 20921-400, Brazil
A. A. Plazas
Affiliation:
Jet Propulsion Laboratory, California Institute of Technology, 4800 Oak Grove Dr., Pasadena, CA 91109, USA Brookhaven National Laboratory, Bldg 510, Upton, NY 11973, USA
A. K. Romer
Affiliation:
Department of Physics and Astronomy, Pevensey Building, University of Sussex, Brighton, BN1 QH, UK
A. Roodman
Affiliation:
Kavli Institute for Particle Astrophysics & Cosmology, P. O. Box 2450, Stanford University, Stanford, CA 94305, USA SLAC National Accelerator Laboratory, Menlo Park, CA 94025, USA
E. S. Rykoff
Affiliation:
Kavli Institute for Particle Astrophysics & Cosmology, P. O. Box 2450, Stanford University, Stanford, CA 94305, USA SLAC National Accelerator Laboratory, Menlo Park, CA 94025, USA
M. Sako
Affiliation:
Department of Physics and Astronomy, University of Pennsylvania, Philadelphia, PA 19104, USA
E. Sanchez
Affiliation:
Centro de Investigaciones Energéticas, Medioambientales y Tecnológicas (CIEMAT), Madrid, Spain
V. Scarpine
Affiliation:
Fermi National Accelerator Laboratory, P. O. Box 500, Batavia, IL 60510, USA
M. Schubnell
Affiliation:
Department of Physics, University of Michigan, Ann Arbor, MI 48109, USA
I. Sevilla
Affiliation:
Department of Astronomy, University of Illinois, 1002 W. Green Street, Urbana, IL 61801, USA Centro de Investigaciones Energéticas, Medioambientales y Tecnológicas (CIEMAT), Madrid, Spain
C. Smith
Affiliation:
Cerro Tololo Inter-American Observatory, National Optical Astronomy Observatory, Casilla 603, La Serena, Chile
M. Soares-Santos
Affiliation:
Fermi National Accelerator Laboratory, P. O. Box 500, Batavia, IL 60510, USA
F. Sobreira
Affiliation:
Fermi National Accelerator Laboratory, P. O. Box 500, Batavia, IL 60510, USA Laboratório Interinstitucional de e-Astronomia - LIneA, Rua Gal. José Cristino 77, Rio de Janeiro, RJ - 20921-400, Brazil
E. Suchyta
Affiliation:
Center for Cosmology and Astro-Particle Physics, The Ohio State University, Columbus, OH 43210, USA Department of Physics, The Ohio State University, Columbus, OH 43210, USA
M. E. C. Swanson
Affiliation:
National Center for Supercomputing Applications, 1205 West Clark St., Urbana, IL 61801, USA
G. Tarle
Affiliation:
Department of Physics, University of Michigan, Ann Arbor, MI 48109, USA
J. Thaler
Affiliation:
Department of Physics, University of Illinois, 1110 W. Green St., Urbana, IL 61801, USA
R. C. Thomas
Affiliation:
Lawrence Berkeley National Laboratory, 1 Cyclotron Road, Berkeley, CA 94720, USA
D. Tucker
Affiliation:
Fermi National Accelerator Laboratory, P. O. Box 500, Batavia, IL 60510, USA
A. R. Walker
Affiliation:
Cerro Tololo Inter-American Observatory, National Optical Astronomy Observatory, Casilla 603, La Serena, Chile
W. Wester
Affiliation:
Fermi National Accelerator Laboratory, P. O. Box 500, Batavia, IL 60510, USA
(The DES Collaboration)
Affiliation:
Fermi National Accelerator Laboratory, P. O. Box 500, Batavia, IL 60510, USA
Rights & Permissions [Opens in a new window]

Abstract

The Dark Energy Survey is undertaking an observational programme imaging 1/4 of the southern hemisphere sky with unprecedented photometric accuracy. In the process of observing millions of faint stars and galaxies to constrain the parameters of the dark energy equation of state, the Dark Energy Survey will obtain pre-discovery images of the regions surrounding an estimated 100 gamma-ray bursts over 5 yr. Once gamma-ray bursts are detected by, e.g., the Swift satellite, the DES data will be extremely useful for follow-up observations by the transient astronomy community. We describe a recently-commissioned suite of software that listens continuously for automated notices of gamma-ray burst activity, collates information from archival DES data, and disseminates relevant data products back to the community in near-real-time. Of particular importance are the opportunities that non-public DES data provide for relative photometry of the optical counterparts of gamma-ray bursts, as well as for identifying key characteristics (e.g., photometric redshifts) of potential gamma-ray burst host galaxies. We provide the functional details of the DESAlert software, and its data products, and we show sample results from the application of DESAlert to numerous previously detected gamma-ray bursts, including the possible identification of several heretofore unknown gamma-ray burst hosts.

Type
Research Article
Copyright
Copyright © Astronomical Society of Australia 2016 

1 INTRODUCTION

The Dark Energy Survey (DES) is an observational programme covering 5 000 deg2 of the southern sky, utilising the DECam instrument (Flaugher et al., Reference Flaugher2015) on the Blanco 4-m telescope at Cerro Tololo Interamerican Observatory, from 2013 August to 2018 February (Abbott et al., Reference Abbott2005). Over the 525 nights of the survey, the DES will observe in five filters broadly similar to the SDSS griz filter set (Gunn et al., Reference Gunn1998), but with some important differences—particularly, higher quantum efficiency at near-infrared wavelengths and the additional Y filter (see Figure 1). DES will reach a photometric accuracy of 1–2% for its two interleaved surveys—the wide-field survey covering the full footprint and the supernova survey covering smaller regions with increased cadence (Diehl et al., Reference Diehl2014)—and it will have significant overlap with other wide-area surveys, such as the Sloan Digital Sky Survey’s Stripe 82 (Abazajian et al., Reference Abazajian2009), the Vista Hemisphere Survey (McMahon et al., Reference McMahon2013), and the South Pole Telescope Survey (Ruhl et al., Reference Ruhl2004). The four Key Science programmes of the DES comprise observations of SNe Ia, large-scale galaxy clustering, galaxy clusters, and weak gravitational lensing; together these four probes will be used to measure the dark energy equation of state with unprecedented precision. The DES data have many uses beyond these primary science goals, however. In this work, we describe a service to provide data products to the transient observational community, particularly related to gamma-ray bursts located within the DES footprint.

Figure 1. Throughput as a function of wavelength of the DECam optical train, including the various filters the DES uses. These throughputs are calculated relative to the use of no filter at 9 000 ${\text{\AA}}$ .

Once a region of the sky has been observed by the DES, those observations will be useful whenever a transient (e.g., a gamma-ray burst) is detected in the same region. The DESAlert system is modelled upon the SkyAlert system (Williams et al., Reference Williams2009), and it bears similarities to the SDSS transient notification system (Cool et al., Reference Cool2006), though it focusses exclusively on data produced by the DES, and (at least initially) it relies on a single source for triggers. When a VOEvent notice (Seaman et al., Reference Seaman2011) is disseminated based on data from the Burst Alert Telescope (BAT) or X-Ray Telescope (XRT) aboard the Swift satellite (Barthelmy et al., Reference Barthelmy1999), the DESAlert system parses the notice for temporal and positional information, and then searches the DES data archives to find all observations of that region. DESAlert then provides finder images of the region derived from DES data, as well as a subset of data derived from DES observations of all nearby stars and galaxies. The finder images show other objects near the GRB, whilst the catalogue of stars provides nearby standards for the purpose of immediate relative photometry. Meanwhile, the galaxy catalogue provides critical information on potential host galaxies—especially magnitude and photometric redshift—for the given GRB. All of these data products will be of use to those who need to make decisions regarding the allocation of scarce resources for follow-up observations of these GRBs. For example, if a potential host has an extremely low (or extremely high) redshift, observers may decide that the GRB warrants further study; thus, they would begin follow-up observations as soon as possible. Alternately, if a host galaxy has unusual colours (as described by the multi-band DES observations), that may be indicative of unusual metallicity or dust content, again prompting observers to allocate observational resources to follow-up the detection of these transients.

By the conclusion of Year 2 in 2015 February, DES had observed nearly all of its survey footprint, with multi-epoch imagery covering the vast majority of that area—see Figure 2. Year 3 saw the completion of coverage of the entire footprint, with subsequent data-taking increasing the number of observed epochs (and thus the effective co-add survey depth) at each position within the footprint. Given the area covered, and assuming randomly distributed GRBs detected by Swift at a rate of ~ 100 yr−1, we expect 10–20 GRBs annually to have DES pre-discovery images that are amenable to analysis and publication via the DESAlert algorithm. We encourage all interested observers to take advantage of these data products provided to the astronomical community.

Figure 2. The DES footprint, with current coverage in the i band shown. Grey is the full survey area, whilst areas covered by 1, 2, 3, 4, or 5+ observations are coloured red, orange, light blue, dark blue, and green, respectively.

2 THE DESALERT ALGORITHM

The DESAlert algorithm is shown schematically in Figure 3. Though it is similar in effect to other VOEvent-focussed software such as dakota or comet (Swinbank, Reference Swinbank2014), in order to maintain minimal dependency on external code, it does not rely on these software packages—it does, however, makes use of the voeventlib libraries (Williams et al., Reference Williams, Peck, Seaman and Comeron2012). The functional code for DESAlert is written in python and sql, and will be made available online via the Astrophysics Source Code Library (Allen et al., Reference Allen2012). It encapsulates all functions necessary to complete the following steps:

  • Listen for VOEvent notices.

  • Select events confirmed to be GRBs by Swift.

  • Parse temporal and position data from notices.

  • Query the DES database for objects and images in a pre-determined region surrounding the GRB.

  • Extract relevant archival DES data for stars and galaxies.

  • Derive data products from extracted data.

  • Publish DESAlert data products as VOEvent notices, as well as to the DESAlert webpage.

Figure 3. The DESAlert algorithm flowchart.

First, the software opens a socket connection and receives (XML-formatted) messages from the VOEvent server as a series of packets, built up byte-by-byte until the full message is received. Error checking ensures that malformed (or improperly received) messages are deleted before the algorithm returns to listen mode. The XML is then parsed to determine the type of message (e.g., imalive, test, or observation). The first of these is necessary for the socket connection to be maintained at the client end, and it results in an identical reply from the client (imalive) so that the socket connection is likewise maintained at the server end. All other types except observation are discarded. Observations, however, have relevant data for each GRB extracted, including GRB position (RA, Dec), burst time, and detecting instrument. In order to ensure near-real-time response from DESAlert to initial GRB detections (almost always by the BAT instrument), we process the initial notice for each GRB. Typically, subsequent VOEvent notices from Swift’s XRT and/or UVOT instruments refine the position (or position uncertainty) of the GRB; in that case, we process the first such subsequent notice in the same fashion as the initial notice, with new and improved information superceding the old. In general, the XRT or UVOT position is sufficiently accurate (with uncertainties of only a few arcseconds) that processing of further notices is not required, so they are discarded.

Based on the extracted GRB data, the python code then calls a custom-made jython-based command-line database interface to query the DES archival catalogue of sources in order to find any stars or galaxies within a box ±60 arcsec in RA and Dec of the BAT-determined GRB position, or within a ±30 arcsec box for XRT-determined position. These objects have previously been extracted from the processed (and co-added) DES Multi-Extension FITS images, and have had a great many parameters determined by the Data Management pipeline (Desai et al., Reference Desai2012) running SourceExtractor (Bertin & Arnouts, Reference Bertin and Arnouts1996), including magnitudes (and uncertainties) for every observed filter, as well as object classification (either star or galaxy). Once the appropriate data are acquired from the DES database they are then stored, via sql Insert commands, in a local database connected to the machine running the DESAlert python code, and are made available to the DESAlert in several ways, as detailed in Section 3. A second query of the DES database finds all g, r, and i images containing the GRB position, from which the nearby objects have been extracted.

Whilst access to this type of data via publicly available online catalogues may be routine for GRB follow-up efforts, DESAlert has the advantage of accessing recent (i.e., proprietary) DES data, which includes objects with better than 2% photometric accuracy down to magnitudes fainter than g = 25. This is significantly fainter than the limits of other catalogues, e.g., the Digitized Sky Survey or 2MASS (Skrutskie at al., Reference Skrutskie2006), often used for this purpose. The DESAlert magnitude limits also compare favorably to those of GROND (Greiner et al., Reference Greiner2008), whilst the expected photometric accuracy of DES is better, and DES also makes use of Y-band imagery.

3 DESALERT DATA PRODUCTS

Once the software implementing the DESAlert algorithm has selected the relevant image segments and composed its catalogue of nearby objects, it disseminates the information of potential relevance to those seeking to follow-up the GRB detection with further observations. The primary data products produced by DESAlert are based on pre-discovery images of the regions around Swift-detected GRBs. The initial products are (XML-formatted) VOEvent notices, listing:

  • basic GRB parameters, including position, time, and discovering instrument;

  • closest stars to the GRB, including positions and magnitudes;

  • closest galaxies to the GRB, including positions, magnitudes, photometric redshifts, and chance alignment probabilities;

  • links to the finder images (in .fz and .jpg formats);

  • links to the DESAlert webpage containing a database with additional information for all objects near the GRB.

Finder images provided by DESAlert are processed and co-added DES images, each of which covers an area on the sky of approximately 0.75° on a side—although depending on the position of the GRB and its associated uncertainty, several images may be provided to completely cover the nominal search area. These images are provided in g, r, and i filters, when available. Links to the original images for each filter are published in the VOEvent notice, and also provided on the DESAlert webpageFootnote 1 .

In order to faciliate relative photometric measurements of the optical counterparts of GRBs, the VOEvent notices provide magnitudes (and uncertainties) in each filter in a range of magnitudes (typically 16–26, depending on filter, with 1–2% photometric accuracy) for several of the nearest stars within the search box centred on the GRB position. Faint stars are far more likely to be positioned closest to any given GRB (and thus included in the XML notice), but information for all stars in the full region is stored in a similar format in the ancillary data available on the DESAlert website.

Likely of greatest interest to those desiring to follow up GRB observations are the DES pre-discovery catalogues of potential GRB host galaxies. The VOEvent notices provide positions, magnitudes (and uncertainties), photometric redshift information (when available), and chance alignment probabilities for the nearest galaxies within the search box centred on the GRB position. Information for all galaxies in this same region is stored in a similar format in the ancillary data available on the DESAlert website. For the subset of galaxies with no photometric redshift determination, we estimate the redshift ‘on the fly’ with an empirical method based upon their relative gri magnitudes (Lopes, Reference Lopes2007). We only report physically realistic (i.e., positive) calculated values in the DESAlert data products (failures of the estimator are assigned a value of ‘-9999’). Because of the limited applicability of this empirical method, we stress that these values are approximations that will be supplanted by more precise determinations from the DES Data Management pipeline as they become available.

Host galaxies of long GRBs are expected to be precisely co-located with the GRB, whilst short GRBs may be several arcseconds distant from their original host. If no galaxy is found within 10 arcsec (for XRT-determined positions), the DESAlert VOEvent notice instead reports the 10σ galaxy detection limit derived (Rykoff, Rozo, & Keisler, Reference Rykoff, Rozo and Keisler2016) from the Year 1 co-added observations – g, r, i = 23.4±0.2, 23.1±0.2, 22.5±0.2, respectively, over the vast majority of the footprint (though not at the edges where there is less overlap amongst observations). A galaxy undetected even down to these limits (especially as the Survey progresses and the limits are extended even fainter) may indicate a high-redshift GRB (and host)—a particularly intriguing target for follow-up observations. Meanwhile, the DESAlert database still includes all objects found within the full search area.

4 TESTING THE DESALERT ALGORITHM

In addition to being submitted to a ‘code review’ by an experienced software engineer, each step of the DESAlert code outlined above has been individually tested, whilst end-to-end tests of the software show that all variations in expected VOEvent inputs are handled properly and output correct data products, even in rare and subtlely challenging cases (e.g., when a GRB is very close to RA=0.0, objects are returned with both RA > 0 and RA < 360). Worth noting also for desalert users are ‘edge’ cases where a GRB is near the border of the DES footprint (cf. Figure 2). Obviously, in this case, DESAlert can only return objects observed by DES, so the output object list may not be a true reflection of all objects that would have been contained within the search box. Therefore, a warning of potential incompleteness in the object catalogue is included with the VOEvent notice.

By the conclusion of Year 1 in 2014 February, DES had observed of order 2 500 deg2 across the southern sky. A variety of tests were performed with the DESAlert algorithm using both the Year 1 and Science Verification (Year 0) data (Diehl et al., Reference Diehl2014). First, we simulated VOEvent notices with systematically varying positions throughout the DES footprint—though of course the real-time nature of DESAlert was not exercised in this way, we nevertheless could determine in a statistical way the expected impact of DESAlert. Of the 45 simulated notices input to DESAlert, 23 were found to have archival DES data matching the locations, each of which had of order 500 objects (stars or galaxies) within 0.1° of the nominal GRB position. This yield is consistent with statistical expectations from a survey region that is 50% covered; with the data from Year 1 to Year 3 (and the start of Year 4) now populating the DES database, the future yield for DESAlert should be significantly higher. Next, to test the real-time processing capabilities of DESAlert, a Swift VOEvent notice of a newly discovered GRB was received and processed, with the software taking approximately 1 min to execute the full algorithm—parsing the notice, querying the relevant DES databases, building the VOEvent notice, and ‘submitting’ it (albeit to an internal recipient rather than the public VOEvent network).

5 RESULTS FROM THE APPLICATION OF DESALERT TO PAST BURSTS

To test the effectiveness of DESAlert specifically in finding GRB host galaxies, we compare positions of DESAlert galaxies with a selection of GRBs previously detected by Swift, along with a number of GRBs detected by BeppoSAX (Boella et al., Reference Boella1997) and the satellites of the InterPlanetary Network (Hurley, Reference Hurley1998) such as HETE-II (Ricker et al., Reference Ricker2002). Many (though not all) of these bursts have spectroscopically confirmed host galaxies drawn from The Optically Unbiased GRB Host (TOUGH) Survey (Hjorth et al., Reference Hjorth2012; Jakobsson et al., Reference Jakobsson2012), as well as GHostS, the GRB Host Studies (Savaglio, Glazebrook, & Le Borgne, Reference Savaglio, Glazebrook and Le Borgne2006). Of the 45 GRBs with XRT or other precise positions for which desalert searches returned nearby objects with photometric redshifts, 16 were deemed sufficiently close (within 5 arcsec) to a galaxy to warrant more detailed examination. These are GRBs 000210, 050219B, 060614A, 061007A, 071227A, 080514B, 080916A, 081109A, 090827A, 091018A, 110206A, 120701A, 140413A, 140928A, 151111A, and 160422A.

Of these 16, seven have spectroscopically confirmed redshifts, whilst nine do not. The galaxy closest to GRB000210 (only 0.21 arcsec away—smaller than the pixel scale of DECam) has a photometric redshift derived by DES of $z_{\text{phot}}$ = 0.844 (a difference of 0.3% from the spectroscopic redshift of $z_{\text{spec}}=0.846$ ). Furthermore, though the filter response curves are different, DES photometry is consistent with previously reported observations of the host galaxy (Christensen, Hjorth, & Gorosabel Reference Christensen, Hjorth and Gorosabel2004). Next, the galaxy closest to GRB081109A (2.93 arcsec away—well within the 4.8 arcsec uncertainty reported by XRT) has a photometric redshift derived by DES of $z_{\text{phot}}$ = 0.906 (a difference of 8.0% from the spectroscopic redshift of $z_{\text{spec}}=0.979$ ). DES photometry is consistent with the upper limits on host magnitude from Swift UVOT (Kuin, Reference Kuin and Immler2008), the REM Telescope (D’Avanzo, Reference D’Avanzo2008) and the Faulkes Telescope-South (Guidorzi, Reference Guidorzi2008), though all have significantly brighter (more than one magnitude) limits than the DES observations. Third, the galaxy closest to GRB091018A (2.16 arcsec away—close, but outside the reported 0.6 arcsec position uncertainty from XRT) has a photometric redshift derived by DES of $z_{\text{phot}}=1.01$ (a difference of 3.9% from the spectroscopic redshift of $z_{\text{spec}}=0.971$ ). Once again, DES photometry matches well with observations of the host galaxy by GROND, Gemini South/GMOS, and Faulkes Telescope-South (Wiersema et al., Reference Wiersema2012). The galaxies closest to the other four GRBs have photometric redshifts that do not reflect the spectroscopic redshifts as accurately (differences of 20% to nearly 50%). See Table 1 for details of all of these candidate host galaxies, as well as the galaxies closest to the GRBs without spectroscopically confirmed redshifts.

Table 1. Parameters of GRB-possible host galaxy matches, including GRB ID, galaxy position, galaxy magnitudes, angular separation from GRB, chance alignment probability, photometric redshift, and spectroscopic redshift.

To determine the general chance alignment of these galaxies with the GRB positions as determined by BAT, we start by finding the total number of galaxies in 80 arcsec × 80 arcsec searchboxes around a large subsample of all GRBs that were tested with DESAlert. Of those searches that returned nearby objects (whether with photo-z’s or not), there were 454 galaxies in the regions around 65 GRBs, totaling 416 000 arcsec2—approximately one galaxy per 910 arcsec2—see Figure 4. The probability of a chance alignment (P $_{\text{chance}}$ ) within a single 3.97 arcsec × 3.97 arcsec box (corresponding to the largest angular separation for any of the potential matches) is simply the ratio of the two areas; that is, P $_{\text{chance}}$ =1.72%. Of course, more than a single search was performed; to determine the chance probability of no matches in 45 searches, we calculate P 45 = [1−P $_{\text{chance}}$ ]45 = 0.458. The chance probability of one or more matches for a single GRB, then, is 1 − P 45 = 0.542, whilst the chance probability of 16 independent matches is P 16 = (0.542)16 = 5.47 × 10−5. That is, the probability that at least one of the 16 matches is the true host of the corresponding GRB is 99.9945%.

Figure 4. (Left) Number of galaxies as a function of magnitude returned by DESAlert from 80 arcsec × 80 arcsec search boxes around 65 GRBs with BAT positions. (Right) The cumulative number of galaxies equal to (or brighter than) each given magnitude. The total number of objects classified as galaxies within these 416 000 arcsec2 is 454, resulting in (on average) approximately one galaxy per 910 arcsec2.

We also performed a more refined chance alignment analysis with a subsample of 38 GRBs with XRT positions, taking into account the individual galaxies’ magnitudes and separations from the nominal GRB positions. The total area searched was 34 200 arcsec2, within which 89 galaxies were found with magnitudes 20 < g < 26 and angular separations < 30 arcsec. We calculate that an angular separation of 0.216 arcsec for a $g<25{\text{th}}$ mag galaxy (cf. GRB00210) has a chance probability of only 0.011%. Similarly, angular separations of 2.93 and 2.16 arcsec for g $<24{\text{th}}$ mag galaxies (cf. GRBs 081109A and 091018A) have a chance probability of 1.78 and 0.968%, respectively, though none of these probabilities factor in the entire ensemble of searches and matches. Figure 5 shows the cumulative chance alignment probability as a function of angular separation and galaxy magnitude. For those wishing to make decisions regarding allocation of valuable follow-up observational resources, suitable cut-offs in probability space can be straightforwardly derived from this information—specifically, we note that any bright ( $g<22{\text{nd}}$ mag) galaxy has a very low probability of chance alignment at any separation up to at least 20 arcsec, whilst any galaxy has a very low probability of being within 5 arcsec, regardless of magnitude. We also note that these probabilities are, if anything, underestimates of the significance of alignments, since the positions used in the analysis are, by construction, correlated with GRBs (and thus with the actual host galaxies), so any arbitrary point on the sky is less likely to have a galaxy nearby than these positions.

Figure 5. Probability of chance alignment of galaxies within the specified angular separation from the nominal GRB position and brighter than the specified magnitude, derived from $30\,\text{arcsec} \times 30\,\text{arcsec}$ search boxes around 38 GRBs. The total number of objects classified as galaxies within these 34 200 arcsec2 is 89, resulting in (on average) approximately one galaxy per 384 arcsec2.

Given the low probabilities associated with the observed galaxy alignments, the correspondence between the spectroscopic and photometric redshifts, and the similarity of the DES and non-DES photometry, we claim that DESAlert has successfully found the host for GRB000210, and has with high probability found the hosts of GRB081109A and GRB091018A. Though in the absence of spectroscopic redshift information, we cannot say with certainty that the additional galaxy matches are indeed the hosts of GRBs 050219B, 080514B, 090827A, 110206A, 120701A, 140413A, 140928A, 151111A, or 160422A, the low probabilities of chance alignment alone provide intriguing evidence pointing in that direction. Most importantly, we show that DESAlert can provide ideal candidates for real-time follow-up observations that could confirm or refute the status of proposed GRB host galaxies.

After completing the search of archival bursts, we tested the DESAlert algorithm with bursts detected by Swift during the Year 1 observations of DES. In Appendices A and B, we present an example of the data products related to a single burst (GRB131105A) that are provided to the community by DESAlert—specifically, the VOEvent XML notice and Finder Image, respectively. Table 2 shows a subset of the relevant data for selected stars and galaxies near this GRB as well. Table entries include object position, grizY magnitudes, the ${spread}\_{model}$ value, object separation from the nominal GRB positon (in arcseconds), object classification, chance alignment probability, photometric redshift, and object number as defined internally by DESAlert (where the same number indicates additional observational epochs of the same object). ${Spread}\_{model}$ is a neural-network based star/galaxy classifier within SourceExtractor; values greater than 0.003 correspond to galaxies (Desai et al., Reference Desai2012). We calculate the mean of the ${spread}\_{model}$ values for all filters, with each weighted by the inverse square of the uncertainty in the observations for that filter, yielding a value that relies most heavily upon the most accurate observation, but still incorporates all available data. As with other GRBs we tested, the photometric redshifts are calculated following the method described above in Section 3—for those observations with physically realistic values (i.e., $z_{\text{phot}}$ > 0), the values derived from repeated observations of the same object are broadly consistent (generally within ±10–20% of one another). All stars and galaxies in the region around the GRB are stored in a similar fashion in the ancillary data available on the DESAlert website.

Table 2. Parameters of selected galaxies and stars near GRB131105A, including position, magnitudes, star/galaxy classification ( ${Spread}\_{Model}$ ) value, angular separation from GRB, star/galaxy classification, chance alignment probability, photometric redshift, and DESAlert object ID (useful particularly for identifying duplicate observations of the same object).

6 FUTURE DEVELOPMENT OF DESALERT

Though DESAlert is fully functional and has made several confirmed (and proposed) matches to GRB hosts, we will continue to improve the functionality of DESAlert through future code releases. Dubbed DESAlert++, the next iteration of the DESAlert algorithm will incorporate several new aspects. In particular, we are investigating alternate search methods to including transients with larger position uncertainties, such as Fermi-detected GRBs with position uncertainties close to 1° (Atwood et al., Reference Atwood2009). One particularly promising category of transients for inclusion in DESAlert++ is Fast Radio Bursts (FRBs, Champion et al., Reference Champion2015). Though their positions are not well known (for example, the beam size of the Parkes radio telescope that has detected numerous FRBs is 14 arcmin in diameter), their redshift can be significantly constrained—for example, Petroff et al. (Reference Petroff2015) have determined that z < 0.5 for FRB 140514. Despite the probability of many galaxies being observed by DES in a search box more than 10 arcmin across, photometric redshift matches could positively identify optical counterparts (or hosts) of FRBs, or at the very least facilitate statistical approaches to identifying likely sources (similar to the method used above to identify multiple potential new GRB host galaxies).

Next, multi-messenger transient detections (e.g., gravitational wave observations) stemming from the Astrophysical Multimessenger Observatory Network (Smith et al., Reference Smith2013) provide yet another potentially interesting avenue of further development of DESAlert. Finally, we are exploring the possibility of including resolved stellar transients such as disappearing red supergiants (Kochanek et al., Reference Kochanek2008) or galactic stellar flares—in such cases, nearby stars (rather than galaxies) become the interesting targets for which the transient follow-up community can utilise the data products released by DESAlert.

7 CONCLUSIONS

DESAlert is an algorithm implemented in python and sql to receive automated notices of GRB parameters from VOEvent Notification triggers, and to provide the astronomical community with finder images as well as catalogues of nearby stars and galaxies (with relevant quantities such as positions, magnitudes, and photometric redshifts). The details of the algorithm have been laid out in the previous Sections, and a sample image and catalogue shown to familiarise readers with the DESAlert data products. Based on tests using historical GRBs, we confirm the effectiveness of DESAlert by finding multiple previously discovered GRB hosts (matching in both position and redshift, and consistent with independent photometric observations), and further by detecting several potential GRB hosts (matching in position, but with no spectroscopic redshift to confirm the DES-derived photometric redshift).

During the 5-yr lifetime of the DES, we expect to provide data products for of order 100 Swift-detected GRBs—though the DESAlert system is expected to function well beyond the formal lifetime of the DES itself. Extensions to DESAlert (dubbed DESAlert++) are currently being explored based upon other GRB-detecting instruments such as Fermi (or based on the incorporation of other transients such as FRBs into the DESAlert algorithm). To estimate the annual rate of GRBs from all sources likely to be within the DES footprint, we searched the SkyAlert database for unique GRB positions from all sources from one year; of the 230 bursts in the database, 34 (roughly 3 per month) fall within the DES footprint and could potentially trigger DESAlert++.

We encourage all observers interested in follow-up observations of transient astrophysical phenomena such as GRBs to avail themselves of the VOEvent notices and web-based data archive provided by DESAlert.

ACKNOWLEDGEMENTS

The authors thank Simon O’Toole for assistance in setup and maintenance of the DESAlert Data Archive at the Australian Astronomical Observatory, as well as Nuria P.F. Lorente for a comprehensive review of the DESAlert code, and Dr. Tayyaba Zafar for providing archival GRB data for testing DESAlert. We also thank the anonymous referee for many insightful comments that prompted numerous improvements in the manuscript and the functionality of the DESAlert code. KK is grateful to the Astronomical Society of Australia for convening an Early Career Researcher Writing Workshop, during which a significant portion of this text was written. The authors also recognise the Macquarie University PACE programme and Professor Mike Ireland for facilitating the undergraduate research opportunities, as well as Raymond Lam for early efforts towards DESAlert. This research has made use of the ghosts database (www.grbhosts.org), which is partly funded by Spitzer/NASA grant RSA Agreement No. 1287913, and has also made use of the simbad database, operated at CDS, Strasbourg, France. We are grateful for the extraordinary contributions of our CTIO colleagues and the DECam Construction, Commissioning and Science Verification teams in achieving the excellent instrument and telescope conditions that have made this work possible. The success of this project also relies critically on the expertise and dedication of the DES Data Management group.

This paper has gone through internal review by the DES collaboration. Funding for the DES Projects has been provided by the U.S. Department of Energy, the U.S. National Science Foundation, the Ministry of Science and Education of Spain, the Science and Technology Facilities Council of the United Kingdom, the Higher Education Funding Council for England, the National Center for Supercomputing Applications at the University of Illinois at Urbana–Champaign, the Kavli Institute of Cosmological Physics at the University of Chicago, the Center for Cosmology and Astro-Particle Physics at the Ohio State University, the Mitchell Institute for Fundamental Physics and Astronomy at Texas A&M University, Financiadora de Estudos e Projetos, Fundação Carlos Chagas Filho de Amparo à Pesquisa do Estado do Rio de Janeiro, Conselho Nacional de Desenvolvimento Científico e Tecnológico and the Ministério da Ciência e Tecnologia, the Deutsche Forschungsgemeinschaft and the Collaborating Institutions in the Dark Energy Survey.

The DES data management system is supported by the National Science Foundation under Grant Number AST-1138766. The DES participants from Spanish institutions are partially supported by MINECO under grants AYA2012-39559, ESP2013-48274, FPA2013-47986, and Centro de Excelencia Severo Ochoa SEV-2012-0234, some of which include ERDF funds from the European Union.

The Collaborating Institutions are Argonne National Laboratory, the University of California at Santa Cruz, the University of Cambridge, Centro de Investigaciones Energeticas, Medioambientales y Tecnologicas-Madrid, the University of Chicago, University College London, the DES-Brazil Consortium, the Eidgenössische Technische Hochschule (ETH) Zürich, Fermi National Accelerator Laboratory, the University of Edinburgh, the University of Illinois at Urbana-Champaign, the Institut de Ciencies de l’Espai (IEEC/CSIC), the Institut de Fisica d’Altes Energies, Lawrence Berkeley National Laboratory, the Ludwig-Maximilians Universität and the associated Excellence Cluster Universe, the University of Michigan, the National Optical Astronomy Observatory, the University of Nottingham, The Ohio State University, the University of Pennsylvania, the University of Portsmouth, SLAC National Accelerator Laboratory, Stanford University, the University of Sussex, and Texas A&M University.

APPENDIX A: SAMPLE DESALERT VOEVENT NOTICE

APPENDIX B: EXAMPLE FINDER IMAGE

Figure B1. This is an example of (a zoomed-in segment of) a GRB Finder Image. The initial Swift-BAT 90% error circle with 1.5 arcmin radius is shown (large red circle), as is the refined Swift-XRT 90% error circle with 1.4 arcsec radius (small white circle).

Footnotes

* Data for this paper is deposited here: https://www.aao.gov.au/DESalert.

References

Abazajian, K., et al. 2009, ApJS, 182, 543 CrossRefGoogle Scholar
Abbott, T., et al. 2005, White Paper Submitted to the Dark Energy Task Force (arXiv:astro-ph/0510346)Google Scholar
Allen, A., et al. 2012, arXiv:1212.1916Google Scholar
Atwood, W. B., et al. 2009, ApJ, 697, 1071 Google Scholar
Barthelmy, S. D., et al. 2005, SSRv, 120, 143 Google Scholar
Bertin, E., & Arnouts, S. A&AS, 117, 393 Google Scholar
Boella, G., et al. 1997, A&AS, 122, 299 Google Scholar
Champion, D. J., et al. 2016, MNRAS, 460, L30 Google Scholar
Christensen, L., Hjorth, J., & Gorosabel, J. 2004, A&A, 425, 913 Google Scholar
Cool, R., et al. 2006, PASP, 118, 733 Google Scholar
D’Avanzo, P., et al. 2008, GCN Circular #8501Google Scholar
Desai, S., et al. 2012, ApJ, 757, 83 CrossRefGoogle Scholar
Diehl, H. T., et al. 2014, in Proc. SPIE 9149, Observatory Operations: Strategies, Processes and Systems V, eds. A. B. Peck, C. R. Benn, and R. L. Seaman (Amsterdam: SPIE), 91490V Google Scholar
Flaugher, B., et al. 2015, AJ, 150, 150F CrossRefGoogle Scholar
Greiner, J., et al. 2008, PASP, 120, 405 CrossRefGoogle Scholar
Guidorzi, C., et al. 2008, GCN Circular #8508Google Scholar
Gunn, J. E., et al. 1998, AJ, 116, 3040 CrossRefGoogle Scholar
Hjorth, , et al. 2012, ApJ, 756, 187 Google Scholar
Hurley, K. 1998, ATel, #19Google Scholar
Jakobsson, , et al. 2012, ApJ, 752, 62 CrossRefGoogle Scholar
Kochanek, C. S. et al. 2008, ApJ, 684, 1336 Google Scholar
Kuin, N. P. M., & Immler, S. 2008, GCN Circular #8504Google Scholar
Lopes, P. A. A. 2007, MNRAS, 380, 1608 Google Scholar
McMahon, R. G., et al. 2013, ESO Messenger, 154, 35 Google Scholar
Petroff, E., et al. 2015, MNRAS, 447, 246P Google Scholar
Ricker, G., et al. 2002, ApJ, 571, L127 CrossRefGoogle Scholar
Ruhl, J. E., et al. 2004, Proc. SPIE, 5498, 11 Google Scholar
Rykoff, E. S., Rozo, E., & Keisler, R. 2016, AJ, submittedGoogle Scholar
Savaglio, S., Glazebrook, K., & Le Borgne, D. 2006, AIP Conf. Proc., 836, 540 Google Scholar
Seaman, R., et al. 2011, IVOA Recommendation (arXiv:1110.0523)Google Scholar
Skrutskie, M. F., et al. 2006, AJ, 131, 1163 Google Scholar
Smith, M. W. E., et al. 2013, APh, V45, 56 Google Scholar
Swinbank, J. 2014, A&C, 7–8, 12 (arXiv:astro-ph/1409.4805; see also http://comet.transientskp.org/)Google Scholar
Wiersema, K., et al. 2012, MNRAS, 426, 2Google Scholar
Williams, R. D., et al. 2009, ASPC, 411, 115W Google Scholar
Williams, R. D., et al. 2012, in Proc SPIE 8448, Responding to the Event Deluge, eds. Peck, A. B., Seaman, R. L., Comeron, F. (Amsterdam: SPIE), 84480RGoogle Scholar
Figure 0

Figure 1. Throughput as a function of wavelength of the DECam optical train, including the various filters the DES uses. These throughputs are calculated relative to the use of no filter at 9 000 ${\text{\AA}}$.

Figure 1

Figure 2. The DES footprint, with current coverage in the i band shown. Grey is the full survey area, whilst areas covered by 1, 2, 3, 4, or 5+ observations are coloured red, orange, light blue, dark blue, and green, respectively.

Figure 2

Figure 3. The DESAlert algorithm flowchart.

Figure 3

Table 1. Parameters of GRB-possible host galaxy matches, including GRB ID, galaxy position, galaxy magnitudes, angular separation from GRB, chance alignment probability, photometric redshift, and spectroscopic redshift.

Figure 4

Figure 4. (Left) Number of galaxies as a function of magnitude returned by DESAlert from 80 arcsec × 80 arcsec search boxes around 65 GRBs with BAT positions. (Right) The cumulative number of galaxies equal to (or brighter than) each given magnitude. The total number of objects classified as galaxies within these 416 000 arcsec2 is 454, resulting in (on average) approximately one galaxy per 910 arcsec2.

Figure 5

Figure 5. Probability of chance alignment of galaxies within the specified angular separation from the nominal GRB position and brighter than the specified magnitude, derived from $30\,\text{arcsec} \times 30\,\text{arcsec}$ search boxes around 38 GRBs. The total number of objects classified as galaxies within these 34 200 arcsec2 is 89, resulting in (on average) approximately one galaxy per 384 arcsec2.

Figure 6

Table 2. Parameters of selected galaxies and stars near GRB131105A, including position, magnitudes, star/galaxy classification (${Spread}\_{Model}$) value, angular separation from GRB, star/galaxy classification, chance alignment probability, photometric redshift, and DESAlert object ID (useful particularly for identifying duplicate observations of the same object).

Figure 7

Figure B1. This is an example of (a zoomed-in segment of) a GRB Finder Image. The initial Swift-BAT 90% error circle with 1.5 arcmin radius is shown (large red circle), as is the refined Swift-XRT 90% error circle with 1.4 arcsec radius (small white circle).