SAGUARO: Time-domain Infrastructure for the Fourth Gravitational-wave Observing Run and Beyond

We present upgraded infrastructure for Searches after Gravitational Waves Using ARizona Observatories (SAGUARO) during LIGO, Virgo, and KAGRA's fourth gravitational-wave (GW) observing run (O4). These upgrades implement many of the lessons we learned after a comprehensive analysis of potential electromagnetic counterparts to the GWs discovered during the previous observing run. We have developed a new web-based target and observation manager (TOM) that allows us to coordinate sky surveys, vet potential counterparts, and trigger follow-up observations from one centralized portal. The TOM includes software that aggregates all publicly available information on the light curves and possible host galaxies of targets, allowing us to rule out potential contaminants like active galactic nuclei, variable stars, solar-system objects, and preexisting supernovae, as well as to assess the viability of any plausible counterparts. We have also upgraded our image-subtraction pipeline by assembling deeper reference images and training a new neural network-based real-bogus classifier. These infrastructure upgrades will aid coordination by enabling the prompt reporting of observations, discoveries, and analysis to the GW follow-up community, and put SAGUARO in an advantageous position to discover kilonovae in the remainder of O4 and beyond. Many elements of our open-source software stack have broad utility beyond multimessenger astronomy, and will be particularly relevant in the "big data" era of transient discoveries by the Vera C. Rubin Observatory.Comment: submitted to AAS Journal

Short GRB Host Galaxies I: Photometric and Spectroscopic Catalogs, Host Associations, and Galactocentric Offsets

We present a comprehensive optical and near-infrared census of the fields of 90 short gamma-ray bursts (GRBs) discovered in 2005-2021, constituting all short GRBs for which host galaxy associations are feasible ($\approx$ 60% of the total Swift short GRB population). We contribute 245 new multi-band imaging observations across 49 distinct GRBs and 25 spectra of their host galaxies. Supplemented by literature and archival survey data, the catalog contains 335 photometric and 40 spectroscopic data sets. The photometric catalog reaches $3\sigma$ depths of $\gtrsim 24-27$ mag and $\gtrsim 23-26$ mag for the optical and near-infrared bands, respectively. We identify host galaxies for 84 bursts, in which the most robust associations make up 54% (49/90) of events, while only a small fraction, 6.7%, have inconclusive host associations. Based on new spectroscopy, we determine 17 host spectroscopic redshifts with a range of $z\approx 0.15-1.6$ and find that $\approx$ 25-44% of Swift short GRBs originate from $z>1$. We also present the galactocentric offset catalog for 83 short GRBs. Taking into account the large range of individual measurement uncertainties, we find a median of projected offset of $\approx 7.9$ kpc, for which the bursts with the most robust associations have a smaller median of $\approx 4.9$ kpc. Our catalog captures more high-redshift and low-luminosity hosts, and more highly-offset bursts than previously found, thereby diversifying the population of known short GRB hosts and properties. In terms of locations and host luminosities, the populations of short GRBs with and without detectable extended emission are statistically indistinguishable. This suggests that they arise from the same progenitors, or from multiple progenitors which form and evolve in similar environments. All of the data products are available on the BRIGHT website.Comment: 53 pages, 9 figures, 6 tables, submitte

The Demographics, Stellar Populations, and Star Formation Histories of Fast Radio Burst Host Galaxies: Implications for the Progenitors

We present a comprehensive catalog of observations and stellar population properties for 23 highly secure host galaxies of fast radio bursts (FRBs). Our sample comprises six repeating FRBs and 17 apparent non-repeaters. We present 82 new photometric and eight new spectroscopic observations of these hosts. Using stellar population synthesis modeling and employing non-parametric star formation histories (SFHs), we find that FRB hosts have a median stellar mass of $\approx 10^{9.9}\,M_{\odot}$, mass-weighted age $\approx 5.1$ Gyr, and ongoing star formation rate $\approx 1.3\,M_{\odot}$ yr$^{-1}$ but span wide ranges in all properties. Classifying the hosts by degree of star formation, we find that 87% (20/23 hosts) are star-forming, two are transitioning, and one is quiescent. The majority trace the star-forming main sequence of galaxies, but at least three FRBs in our sample originate in less active environments (two non-repeaters and one repeater). Across all modeled properties, we find no statistically significant distinction between the hosts of repeaters and non-repeaters. However, the hosts of repeating FRBs generally extend to lower stellar masses, and the hosts of non-repeaters arise in more optically luminous galaxies. While four of the galaxies with the most clear and prolonged rises in their SFHs all host repeating FRBs, demonstrating heightened star formation activity in the last $\lesssim 100$ Myr, one non-repeating host shows this SFH as well. Our results support progenitor models with short delay channels (i.e., magnetars formed via core-collapse supernova) for most FRBs, but the presence of some FRBs in less active environments suggests a fraction form through more delayed channels.Comment: 52 pages, 32 figures, 6 tables, submitte

A Radio Flare in the Long-Lived Afterglow of the Distant Short GRB 210726A: Energy Injection or a Reverse Shock from Shell Collisions?

We present the discovery of the radio afterglow of the short $\gamma$-ray burst (GRB) 210726A, localized to a galaxy at a photometric redshift of $z\sim 2.4$. While radio observations commenced $\lesssim 1~$day after the burst, no radio emission was detected until $\sim11$~days. The radio afterglow subsequently brightened by a factor of $\sim 3$ in the span of a week, followed by a rapid decay (a ``radio flare''). We find that a forward shock afterglow model cannot self-consistently describe the multi-wavelength X-ray and radio data, and underpredicts the flux of the radio flare by a factor of $\approx 5$. We find that the addition of substantial energy injection, which increases the isotropic kinetic energy of the burst by a factor of $\approx 4$, or a reverse shock from a shell collision are viable solutions to match the broad-band behavior. At $z\sim 2.4$, GRB\,210726A is among the highest redshift short GRBs discovered to date as well as the most luminous in radio and X-rays. Combining and comparing all previous radio afterglow observations of short GRBs, we find that the majority of published radio searches conclude by $\lesssim 10~$days after the burst, potentially missing these late rising, luminous radio afterglows.Comment: 28 pages, 10 figures, submitted to Ap

Heavy element production in a compact object merger observed by JWST

The mergers of binary compact objects such as neutron stars and black holes are of central interest to several areas of astrophysics, including as the progenitors of gamma-ray bursts (GRBs) 1, sources of high-frequency gravitational waves (GWs) 2 and likely production sites for heavy-element nucleosynthesis by means of rapid neutron capture (the r-process) 3. Here we present observations of the exceptionally bright GRB 230307A. We show that GRB 230307A belongs to the class of long-duration GRBs associated with compact object mergers 4–6 and contains a kilonova similar to AT2017gfo, associated with the GW merger GW170817 (refs. 7–12). We obtained James Webb Space Telescope (JWST) mid-infrared imaging and spectroscopy 29 and 61 days after the burst. The spectroscopy shows an emission line at 2.15 microns, which we interpret as tellurium (atomic mass A = 130) and a very red source, emitting most of its light in the mid-infrared owing to the production of lanthanides. These observations demonstrate that nucleosynthesis in GRBs can create r-process elements across a broad atomic mass range and play a central role in heavy-element nucleosynthesis across the Universe

A kilonova following a long-duration gamma-ray burst at 350 Mpc

Full list of authors: Rastinejad, Jillian C.; Gompertz, Benjamin P.; Levan, Andrew J.; Fong, Wen-Fai; Nicholl, Matt; Lamb, Gavin P.; Malesani, Daniele B.; Nugent, Anya E.; Oates, Samantha R.; Tanvir, Nial R.; Postigo, Antonio de Ugarte; Kilpatrick, Charles D.; Moore, Christopher J.; Metzger, Brian D.; Ravasio, Maria Edvige; Rossi, Andrea; Schroeder, Genevieve; Jencson, Jacob; Sand, David J.; Smith, Nathan; Fernandez, Jose Feliciano Agui; Berger, Edo; Blanchard, Peter K.; Chornock, Ryan; Cobb, Bethany E.; De Pasquale, Massimiliano; Fynbo, Johan P. U.; Izzo, Luca; Kann, D. Alexander; Laskar, Tanmoy; Marini, Ester; Paterson, Kerry; Escorial, Alicia Rouco; Sears, Huei M.; Thone, Christina C.Gamma-ray bursts (GRBs) are divided into two populations1,2; long GRBs that derive from the core collapse of massive stars (for example, ref. 3) and short GRBs that form in the merger of two compact objects4,5. Although it is common to divide the two populations at a gamma-ray duration of 2 s, classification based on duration does not always map to the progenitor. Notably, GRBs with short (≲2 s) spikes of prompt gamma-ray emission followed by prolonged, spectrally softer extended emission (EE-SGRBs) have been suggested to arise from compact object mergers6,7,8. Compact object mergers are of great astrophysical importance as the only confirmed site of rapid neutron capture (r-process) nucleosynthesis, observed in the form of so-called kilonovae9,10,11,12,13,14. Here we report the discovery of a possible kilonova associated with the nearby (350 Mpc), minute-duration GRB 211211A. The kilonova implies that the progenitor is a compact object merger, suggesting that GRBs with long, complex light curves can be spawned from merger events. The kilonova of GRB 211211A has a similar luminosity, duration and colour to that which accompanied the gravitational wave (GW)-detected binary neutron star (BNS) merger GW170817 (ref. 4). Further searches for GW signals coincident with long GRBs are a promising route for future multi-messenger astronomy. © 2022, The Author(s), under exclusive licence to Springer Nature Limited.The Fong group at Northwestern acknowledges support by the National Science Foundation under grant nos. AST-1814782 and AST-1909358 and CAREER grant no. AST-2047919. W.F. gratefully acknowledges support by the David and Lucile Packard Foundation. A.J.L. and D.B.M. are supported by the European Research Council (ERC) under the European Union’s Horizon 2020 research and innovation programme (grant agreement no. 725246). M.N. and B.P.G. are supported by the ERC under the European Union’s Horizon 2020 research and innovation programme (grant agreement no. 948381). M.N. acknowledges a Turing Fellowship. G.P.L. is supported by the UK Science and Technology Facilities Council grant ST/S000453/1. A.R. and E.M. acknowledge support from the INAF research project ‘LBT - Supporto Arizona Italia’. J.F.A.F. acknowledges support from the Spanish Ministerio de Ciencia, Innovación y Universidades through the grant PRE2018-086507. D.A.K. and J.F.A.F. acknowledge support from Spanish National Research Project RTI2018-098104-J-I00 (GRBPhot). W. M. Keck Observatory and MMT Observatory access was supported by Northwestern University and the Center for Interdisciplinary Exploration and Research in Astrophysics (CIERA). Some of the data presented herein were obtained at the W. M. Keck Observatory, which is operated as a scientific partnership among the California Institute of Technology, the University of California and the National Aeronautics and Space Administration (NASA). The Observatory was made possible by the generous financial support of the W. M. Keck Foundation. We wish to recognize and acknowledge the very important cultural role and reverence that the summit of Maunakea has always had within the indigenous Hawaiian community. We are most fortunate to have the opportunity to conduct observations from this mountain. Observations reported here were obtained at the MMT Observatory, a joint facility of the University of Arizona and the Smithsonian Institution. On the basis of observations obtained at the international Gemini Observatory (programme ID GN2021B-Q-109), a programme of NOIRLab, which is managed by the Association of Universities for Research in Astronomy (AURA) under a cooperative agreement with the National Science Foundation on behalf of the Gemini Observatory partnership: the National Science Foundation (United States), National Research Council (Canada), Agencia Nacional de Investigación y Desarrollo (Chile), Ministerio de Ciencia, Tecnología e Innovación (Argentina), Ministério da Ciência, Tecnologia, Inovações e Comunicações (Brazil) and Korea Astronomy and Space Science Institute (Republic of Korea). Processed using the Gemini IRAF package and DRAGONS (Data Reduction for Astronomy from Gemini Observatory North and South). This work made use of data supplied by the UK Swift Science Data Centre at the University of Leicester. The National Radio Astronomy Observatory is a facility of the National Science Foundation operated under cooperative agreement by Associated Universities, Inc. This research is based on observations made with the NASA/ESA Hubble Space Telescope obtained from the Space Telescope Science Institute, which is operated by the AURA, Inc., under NASA contract NAS 5-26555. These observations are associated with programme no. 16923. This work is partly based on observations made with the Gran Telescopio Canarias, installed at the Spanish Observatorio del Roque de los Muchachos of the Instituto de Astrofísica de Canarias, on the island of La Palma. Partly based on observations collected at the Calar Alto Astronomical Observatory, operated jointly by Instituto de Astrofísica de Andalucía (CSIC) and Junta de Andalucía. Partly based on observations made with the Nordic Optical Telescope, under programme 64-502, owned in collaboration by the University of Turku and Aarhus University, and operated jointly by Aarhus University, the University of Turku and the University of Oslo, representing Denmark, Finland and Norway, respectively, the University of Iceland and Stockholm University at the Observatorio del Roque de los Muchachos, La Palma, Spain, of the Instituto de Astrofísica de Canarias. The LBT is an international collaboration among institutions in the United States, Italy and Germany. LBT Corporation partners are: The University of Arizona on behalf of the Arizona Board of Regents; Istituto Nazionale di Astrofisica, Italy; LBT Beteiligungsgesellschaft, Germany, representing the Max Planck Society, The Leibniz Institute for Astrophysics Potsdam and Heidelberg University; The Ohio State University, representing OSU, University of Notre Dame, University of Minnesota and University of Virginia.Peer reviewe

SAGUARO: Time-domain Infrastructure for the Fourth Gravitational-wave Observing Run and Beyond

We present upgraded infrastructure for Searches After Gravitational waves Using ARizona Observatories (SAGUARO) during LIGO, Virgo, and KAGRA’s fourth gravitational-wave (GW) observing run (O4). These upgrades implement many of the lessons we learned after a comprehensive analysis of potential electromagnetic counterparts to the GWs discovered during the previous observing run. We have developed a new web-based target and observation manager (TOM) that allows us to coordinate sky surveys, vet potential counterparts, and trigger follow-up observations from one centralized portal. The TOM includes software that aggregates all publicly available information on the light curves and possible host galaxies of targets, allowing us to rule out potential contaminants like active galactic nuclei, variable stars, solar system objects, and preexisting supernovae, as well as to assess the viability of any plausible counterparts. We have also upgraded our image-subtraction pipeline by assembling deeper reference images and training a new neural-network-based real–bogus classifier. These infrastructure upgrades will aid coordination by enabling the prompt reporting of observations, discoveries, and analysis to the GW follow-up community, and put SAGUARO in an advantageous position to discover kilonovae in the remainder of O4 and beyond. Many elements of our open-source software stack have broad utility beyond multimessenger astronomy, and will be particularly relevant in the “big data” era of transient discoveries by the Vera C. Rubin Observatory

A long-duration gamma-ray burst of dynamical origin from the nucleus of an ancient galaxy

The majority of long-duration (&gt;2 s) gamma-ray bursts (GRBs) arise from the collapse of massive stars, with a small proportion created from the merger of compact objects. Most of these systems form via standard stellar evolution pathways. However, a fraction of GRBs may result from dynamical interactions in dense environments. These channels could also contribute substantially to the samples of compact object mergers detected as gravitational wave sources. Here we report the case of GRB 191019A, a long GRB (a duration of T 90 = 64.4 ± 4.5 s), which we pinpoint close (⪅100 pc projected) to the nucleus of an ancient (&gt;1 Gyr old) host galaxy at z = 0.248. The lack of evidence for star formation and deep limits on any supernova emission disfavour a massive star origin. The most likely route for progenitor formation is via dynamical interactions in the dense nucleus of the host. The progenitor, in this case, could be a compact object merger. These may form in dense nuclear clusters or originate in a gaseous disc around the supermassive black hole. Identifying, to the best of our knowledge, a first example of a dynamically produced GRB demonstrates the role that such bursts may have in probing dense environments and constraining dynamical fractions in gravitational wave populations.</p

A kilonova following a long-duration gamma-ray burst at 350 Mpc

Here, we report the discovery of a kilonova associated with the nearby (350 Mpc) minute-duration GRB 211211A. In tandem with deep optical limits that rule out the presence of an accompanying supernova to $M_I > -13$ mag at 17.7 days post-burst, the identification of a kilonova confirms that this burst's progenitor was a compact object merger. While the spectrally softer tail in GRB 211211A's gamma-ray light curve is reminiscent of previous extended emission short GRBs (EE-SGRBs), its prompt, bright spikes last $\gtrsim 12$ s, separating it from past EE-SGRBs. GRB 211211A's kilonova has a similar luminosity, duration and color to AT2017gfo, the kilonova found in association with the gravitational wave (GW)-detected binary neutron star (BNS) merger GW170817. We find that the merger ejected $\approx 0.04 M_{\odot}$ of r-process-rich material, and is consistent with the merger of two neutron stars (NSs) with masses close to the canonical $1.4 M_{\odot}$. This discovery implies that GRBs with long, complex light curves can be spawned from compact object merger events and that a population of kilonovae following GRBs with durations $\gg 2$ s should be accounted for in calculations of the NS merger r-process contribution and rate. At 350 Mpc, the current network of GW interferometers at design sensitivity would have detected the merger precipitating GRB 211211A, had it been operating at the time of the event. Further searches for GW signals coincident with long GRBs are therefore a promising route for future multi-messenger astronomy.Comment: Submitted. 69 pages, 11 figures, 3 table