Showing 1–40 of 40 results for author: Ahumada, T

Searching for electromagnetic emission in an AGN from the gravitational wave binary black hole merger candidate S230922g

Authors: Tomás Cabrera, Antonella Palmese, Lei Hu, Brendan O'Connor, K. E. Saavik Ford, Barry McKernan, Igor Andreoni, Tomás Ahumada, Ariel Amsellem, Malte Busmann, Peter Clark, Michael W. Coughlin, Ekaterine Dadiani, Veronica Diaz, Matthew J. Graham, Daniel Gruen, Keerthi Kunnumkai, Jake Postiglione, Julian S. Sommer, Francisco Valdes

Abstract: We carried out long-term monitoring of the LIGO/Virgo/KAGRA binary black hole (BBH) merger candidate S230922g in search of electromagnetic emission from the interaction of the merger remnant with an embedding active galactic nuclei (AGN) accretion disk. Using a dataset primarily composed of wide-field imaging from the Dark Energy Camera (DECam) and supplemented by additional photometric and spectr… ▽ More We carried out long-term monitoring of the LIGO/Virgo/KAGRA binary black hole (BBH) merger candidate S230922g in search of electromagnetic emission from the interaction of the merger remnant with an embedding active galactic nuclei (AGN) accretion disk. Using a dataset primarily composed of wide-field imaging from the Dark Energy Camera (DECam) and supplemented by additional photometric and spectroscopic resources, we searched ~ 70% of the sky area probability for transient phenomena, and discovered 6 counterpart candidates. One especially promising candidate - AT 2023aagj - exhibited temporally varying asymmetric components in spectral broad line regions, a feature potentially indicative of an off-center event such as a BBH merger. This represents the first live search and multiwavelength, photometric, and spectroscopic monitoring of a GW BBH optical counterpart candidate in the disk of an AGN. △ Less

Submitted 15 July, 2024; originally announced July 2024.

Comments: 25 pages, 9 figures. Comments welcome

Searching for gravitational wave optical counterparts with the Zwicky Transient Facility: summary of O4a

Authors: Tomás Ahumada, Shreya Anand, Michael W. Coughlin, Vaidehi Gupta, Mansi M. Kasliwal, Viraj R. Karambelkar, Robert D. Stein, Gaurav Waratkar, Vishwajeet Swain, Theophile Jegou du Laz, Akash Anumarlapudi, Igor Andreoni, Mattia Bulla, Gokul P. Srinivasaragavan, Andrew Toivonen, Avery Wold, Eric C. Bellm, S. Bradley Cenko, David L. Kaplan, Jesper Sollerman, Varun Bhalerao, Daniel Perley, Anirudh Salgundi, Aswin Suresh, K-Ryan Hinds , et al. (27 additional authors not shown)

Abstract: During the first half of the fourth observing run (O4a) of the International Gravitational Wave Network (IGWN), the Zwicky Transient Facility (ZTF) conducted a systematic search for kilonova (KN) counterparts to binary neutron star (BNS) and neutron star-black hole (NSBH) merger candidates. Here, we present a comprehensive study of the five high-significance (FAR < 1 per year) BNS and NSBH candida… ▽ More During the first half of the fourth observing run (O4a) of the International Gravitational Wave Network (IGWN), the Zwicky Transient Facility (ZTF) conducted a systematic search for kilonova (KN) counterparts to binary neutron star (BNS) and neutron star-black hole (NSBH) merger candidates. Here, we present a comprehensive study of the five high-significance (FAR < 1 per year) BNS and NSBH candidates in O4a. Our follow-up campaigns relied on both target-of-opportunity observations (ToO) and re-weighting of the nominal survey schedule to maximize coverage. We describe the toolkit we have been developing, Fritz, an instance of SkyPortal, instrumental in coordinating and managing our telescope scheduling, candidate vetting, and follow-up observations through a user-friendly interface. ZTF covered a total of 2841 deg$^2$ within the skymaps of the high-significance GW events, reaching a median depth of g~20.2 mag. We circulated 15 candidates, but found no viable KN counterpart to any of the GW events. Based on the ZTF non-detections of the high-significance events in O4a, we used a Bayesian approach, nimbus, to quantify the posterior probability of KN model parameters that are consistent with our non-detections. Our analysis favors KNe with initial absolute magnitude fainter than -16 mag. The joint posterior probability of a GW170817-like KN associated with all our O4a follow-ups was 64%. Additionally, we use a survey simulation software, simsurvey, to determine that our combined filtered efficiency to detect a GW170817-like KN is 36%, when considering the 5 confirmed astrophysical events in O3 (1 BNS and 4 NSBH), along with our O4a follow-ups. Following Kasliwal et al. (2020), we derived joint constraints on the underlying KN luminosity function based on our O3 and O4a follow-ups, determining that no more than 76% of KNe fading at 1 mag/day can peak at a magnitude brighter than -17.5 mag. △ Less

Submitted 20 May, 2024; originally announced May 2024.

Comments: submitted

Kilonova parameters estimation with LSST at Vera C. Rubin Observatory

Authors: Fabio Ragosta, Tomas Ahumada, Silvia Piranomonte, Igor Andreoni, Andrea Melandri, Alberto Colombo, Michael W. Coughlin

Abstract: The up-coming Vera Rubin Observatory's Legacy Survey of Space and Time (LSST) opens a new opportunity to rapidly survey the southern Sky at optical wavelenghts (\ie ugrizy bands). In this study, we aim to test the possibility of using LSST observations to constrain the mass and velocity of different KN ejecta components from the observation of a combined set of light curves from GRB afterglows and… ▽ More The up-coming Vera Rubin Observatory's Legacy Survey of Space and Time (LSST) opens a new opportunity to rapidly survey the southern Sky at optical wavelenghts (\ie ugrizy bands). In this study, we aim to test the possibility of using LSST observations to constrain the mass and velocity of different KN ejecta components from the observation of a combined set of light curves from GRB afterglows and KNæ. We used a sample of simulated light curves from the aforementioned events as they would have been seen during the LSST survey to study how observing strategies' choice impacts the parameters' estimation. We found that the observing strategy design to be the best compromise between light curves coverage, observed filters and fit reliability involves a high number of visits with long gap pairs of about 4 hours every 2 nights in the same or different filters. The features of the observing strategy will allow us to recognize the different stages of the light curve evolution and gather observations in at least 3 filters. △ Less

Submitted 20 March, 2024; originally announced March 2024.

Dramatic rebrightening of the type-changing stripped-envelope supernova SN 2023aew

Authors: Yashvi Sharma, Jesper Sollerman, Shrinivas R. Kulkarni, Takashi J. Moriya, Steve Schulze, Stan Barmentloo, Michael Fausnaugh, Avishay Gal-Yam, Anders Jerkstrand, Tomás Ahumada, Eric C. Bellm, Kaustav K. Das, Andrew Drake, Christoffer Fremling, Saarah Hall, K. R. Hinds, Theophile Jegou du Laz, Viraj Karambelkar, Mansi M. Kasliwal, Frank J. Masci, Adam A. Miller, Guy Nir, Daniel A. Perley, Josiah N. Purdum, Yu-Jing Qin , et al. (10 additional authors not shown)

Abstract: Multi-peaked supernovae with precursors, dramatic light-curve rebrightenings, and spectral transformation are rare, but are being discovered in increasing numbers by modern night-sky transient surveys like the Zwicky Transient Facility (ZTF). Here, we present the observations and analysis of SN 2023aew, which showed a dramatic increase in brightness following an initial luminous (-17.4 mag) and lo… ▽ More Multi-peaked supernovae with precursors, dramatic light-curve rebrightenings, and spectral transformation are rare, but are being discovered in increasing numbers by modern night-sky transient surveys like the Zwicky Transient Facility (ZTF). Here, we present the observations and analysis of SN 2023aew, which showed a dramatic increase in brightness following an initial luminous (-17.4 mag) and long (~100 days) unusual first peak (possibly precursor). SN 2023aew was classified as a Type IIb supernova during the first peak but changed its type to resemble a stripped-envelope supernova (SESN) after the marked rebrightening. We present comparisons of SN 2023aew's spectral evolution with SESN subtypes and argue that it is similar to SNe Ibc during its main peak. P-Cygni Balmer lines are present during the first peak, but vanish during the second peak's photospheric phase, before H$α$ resurfaces again during the nebular phase. The nebular lines ([O I], [Ca II], Mg I], H$α$) exhibit a double-peaked structure which hints towards a clumpy or non-spherical ejecta. We analyze the second peak in the light curve of SN 2023aew and find it to be broader than normal SESNe as well as requiring a very high $^{56}$Ni mass to power the peak luminosity. We discuss the possible origins of SN 2023aew including an eruption scenario where a part of the envelope is ejected during the first peak which also powers the second peak of the light curve through SN-CSM interaction. △ Less

Submitted 5 February, 2024; originally announced February 2024.

Comments: 22 pages, 11 figures, 5 tables

AT2019pim: A Luminous Orphan Afterglow from a Moderately Relativistic Outflow

Authors: Daniel A. Perley, Anna Y. Q. Ho, Michael Fausnaugh, Gavin P. Lamb, Mansi M. Kasliwal, Tomas Ahumada, Shreya Anand, Igor Andreoni, Eric Bellm, Varun Bhalerao, Bryce Bolin, Thomas G. Brink, Eric Burns, S. Bradley Cenko, Alessandra Corsi, Alexei V. Filippenko, Dmitry Frederiks, Adam Goldstein, Rachel Hamburg, Rahul Jayaraman, Peter G. Jonker, Erik C. Kool, Shrinivas Kulkarni, Harsh Kumar, Russ Laher , et al. (12 additional authors not shown)

Abstract: Classical gamma-ray bursts (GRBs) have two distinct emission episodes: prompt emission from ultra-relativistic ejecta and afterglow from shocked circumstellar material. While both components are extremely luminous in known GRBs, a variety of scenarios predict the existence of luminous afterglow emission with little or no associated high-energy prompt emission. We present AT 2019pim, the first secu… ▽ More Classical gamma-ray bursts (GRBs) have two distinct emission episodes: prompt emission from ultra-relativistic ejecta and afterglow from shocked circumstellar material. While both components are extremely luminous in known GRBs, a variety of scenarios predict the existence of luminous afterglow emission with little or no associated high-energy prompt emission. We present AT 2019pim, the first secure example of this phenomenon to be identified. Serendipitously discovered during follow-up observations of a gravitational-wave trigger and located in a contemporaneous TESS sector, it is hallmarked by a fast-rising (t ~ 2 hr), luminous (M_UV,peak ~ -24.4 mag) optical transient with accompanying luminous X-ray and radio emission. No gamma-ray emission consistent with the time and location of the transient was detected by Fermi-GBM or by Konus, placing strong limits on an accompanying GRB. We investigate several independent observational aspects of the afterglow in the context of constraints on relativistic motion and find all of them are consistent with an initial Lorentz factor of Gamma_0 ~ 30-50, significantly lower than in any well-observed GRB and consistent with the theoretically-predicted "dirty fireball" scenario in which the high-energy prompt emission is stifled by pair production. However, we cannot rule out a structured jet model in which only the line-of-sight material was ejected at low-Gamma, off-axis from a classical high-Gamma jet core. This event represents a milestone in orphan afterglow searches, demonstrating that luminous afterglows with weak or no detectable gamma-ray radiation exist in nature and can be discovered by high-cadence optical surveys. △ Less

Submitted 29 January, 2024; originally announced January 2024.

Comments: Submitted to MNRAS

The 2022-2023 accretion outburst of the young star V1741 Sgr

Authors: Michael A. Kuhn, Lynne A. Hillenbrand, Michael S. Connelley, R. Michael Rich, Bart Staels, Adolfo S. Carvalho, Philip W. Lucas, Christoffer Fremling, Viraj R. Karambelkar, Ellen Lee, Tomás Ahumada, Emille E. O. Ishida, Kishalay De, Rafael S. de Souza, Mansi Kasliwal

Abstract: V1741 Sgr (= SPICY 71482/Gaia22dtk) is a Classical T Tauri star on the outskirts of the Lagoon Nebula. After at least a decade of stability, in mid-2022, the optical source brightened by ~3 mag over two months, remained bright until early 2023, then dimmed erratically over the next four months. This event was monitored with optical and infrared spectroscopy and photometry. Spectra from the peak (O… ▽ More V1741 Sgr (= SPICY 71482/Gaia22dtk) is a Classical T Tauri star on the outskirts of the Lagoon Nebula. After at least a decade of stability, in mid-2022, the optical source brightened by ~3 mag over two months, remained bright until early 2023, then dimmed erratically over the next four months. This event was monitored with optical and infrared spectroscopy and photometry. Spectra from the peak (October 2022) indicate an EX Lup-type (EXor) accretion outburst, with strong emission from H I, He I, and Ca II lines and CO bands. At this stage, spectroscopic absorption features indicated a temperature of T ~ 4750 K with low-gravity lines (e.g., Ba II and Sr II). By April 2023, with the outburst beginning to dim, strong TiO absorption appeared, indicating a cooler T ~ 3600 K temperature. However, once the source had returned to its pre-outburst flux in August 2023, the TiO absorption and the CO emission disappeared. When the star went into outburst, the source's spectral energy distribution became flatter, leading to bluer colours at wavelengths shorter than ~1.6 microns and redder colours at longer wavelengths. The brightening requires a continuum emitting area larger than the stellar surface, likely from optically thick circumstellar gas with cooler surface layers producing the absorption features. Additional contributions to the outburst spectrum may include blue excess from hotspots on the stellar surface, emission lines from diffuse gas, and reprocessed emission from the dust disc. Cooling of the circumstellar gas would explain the appearance of TiO, which subsequently disappeared once this gas had faded and the stellar spectrum reemerged. △ Less

Submitted 17 January, 2024; originally announced January 2024.

Comments: Accepted for publication in MNRAS; 17 pages, 16 figures, and 2 tables

Windows on the Universe: Establishing the Infrastructure for a Collaborative Multi-messenger Ecosystem

Authors: The 2023 Windows on the Universe Workshop White Paper Working Group, T. Ahumada, J. E. Andrews, S. Antier, E. Blaufuss, P. R. Brady, A. M. Brazier, E. Burns, S. B. Cenko, P. Chandra, D. Chatterjee, A. Corsi, M. W. Coughlin, D. A. Coulter, S. Fu, A. Goldstein, L. P. Guy, E. J. Hooper, S. B. Howell, T. B. Humensky, J. A. Kennea, S. M. Jarrett, R. M. Lau, T. R. Lewis, L. Lu , et al. (21 additional authors not shown)

Abstract: In this White Paper, we present recommendations for the scientific community and funding agencies to foster the infrastructure for a collaborative multi-messenger and time-domain astronomy (MMA/TDA) ecosystem. MMA/TDA is poised for breakthrough discoveries in the coming decade. In much the same way that expanding beyond the optical bandpass revealed entirely new and unexpected discoveries, cosmic… ▽ More In this White Paper, we present recommendations for the scientific community and funding agencies to foster the infrastructure for a collaborative multi-messenger and time-domain astronomy (MMA/TDA) ecosystem. MMA/TDA is poised for breakthrough discoveries in the coming decade. In much the same way that expanding beyond the optical bandpass revealed entirely new and unexpected discoveries, cosmic messengers beyond light (i.e., gravitational waves, neutrinos, and cosmic rays) open entirely new windows to answer some of the most fundamental questions in (astro)physics: heavy element synthesis, equation of state of dense matter, particle acceleration, etc. This field was prioritized as a frontier scientific pursuit in the 2020 Decadal Survey on Astronomy and Astrophysics via its "New Windows on the Dynamic Universe" theme. MMA/TDA science presents technical challenges distinct from those experienced in other disciplines. Successful observations require coordination across myriad boundaries -- different cosmic messengers, ground vs. space, international borders, etc. -- all for sources that may not be well localized, and whose brightness may be changing rapidly with time. Add that all of this work is undertaken by real human beings, with distinct backgrounds, experiences, cultures, and expectations, that often conflict. To address these challenges and help MMA/TDA realize its full scientific potential in the coming decade (and beyond), the second in a series of community workshops sponsored by the U.S. National Science Foundation (NSF) and NASA titled "Windows on the Universe: Establishing the Infrastructure for a Collaborative Multi-Messenger Ecosystem" was held on October 16-18, 2023 in Tucson, AZ. Here we present the primary recommendations from this workshop focused on three key topics -- hardware, software, and people and policy. [abridged] △ Less

Submitted 3 April, 2024; v1 submitted 3 January, 2024; originally announced January 2024.

Comments: Workshop white paper

Characterizing the Ordinary Broad-lined Type Ic SN 2023pel from the Energetic GRB 230812B

Authors: Gokul P. Srinivasaragavan, Vishwajeet Swain, Brendan M. O'Connor, Shreya Anand, Tomás Ahumada, Daniel A. Perley, Robert Stein, Jesper Sollerman, Christoffer Fremling, S. Bradley Cenko, Sarah Antier, Nidhal Guessoum, Thomas Hussenot-Desenonges, Patrice Hello, Stephen Lesage, Erica Hammerstein, M. Coleman Miller, Igor Andreoni, Varun Bhalerao, Joshua S. Bloom, Anirban Dutta, Avishay Gal-Yam, K-Ryan Hinds, Amruta D. Jaodand, Mansi M. Kasliwal , et al. (17 additional authors not shown)

Abstract: We report observations of the optical counterpart of the long gamma-ray burst (LGRB) GRB 230812B, and its associated supernova (SN) SN 2023pel. The proximity ($z = 0.36$) and high energy ($E_{γ, \rm{iso}} \sim 10^{53}$ erg) make it an important event to study as a probe of the connection between massive star core-collapse and relativistic jet formation. With a phenomenological power-law model for… ▽ More We report observations of the optical counterpart of the long gamma-ray burst (LGRB) GRB 230812B, and its associated supernova (SN) SN 2023pel. The proximity ($z = 0.36$) and high energy ($E_{γ, \rm{iso}} \sim 10^{53}$ erg) make it an important event to study as a probe of the connection between massive star core-collapse and relativistic jet formation. With a phenomenological power-law model for the optical afterglow, we find a late-time flattening consistent with the presence of an associated SN. SN 2023pel has an absolute peak $r$-band magnitude of $M_r = -19.46 \pm 0.18$ mag (about as bright as SN 1998bw) and evolves on quicker timescales. Using a radioactive heating model, we derive a nickel mass powering the SN of $M_{\rm{Ni}} = 0.38 \pm 0.01$ $\rm{M_\odot}$, and a peak bolometric luminosity of $L_{\rm{bol}} \sim 1.3 \times 10^{43}$ $\rm{erg}$ $\rm{s^{-1}}$. We confirm SN 2023pel's classification as a broad-lined Type Ic SN with a spectrum taken 15.5 days after its peak in $r$ band, and derive a photospheric expansion velocity of $v_{\rm{ph}} = 11,300 \pm 1,600$ $\rm{km}$ $\rm{s^{-1}}$ at that phase. Extrapolating this velocity to the time of maximum light, we derive the ejecta mass $M_{\rm{ej}} = 1.0 \pm 0.6$ $\rm{M_\odot}$ and kinetic energy $E_{\rm{KE}} = 1.3^{+3.3}_{-1.2} \times10^{51}$ $\rm{erg}$. We find that GRB 230812B/SN 2023pel has SN properties that are mostly consistent with the overall GRB-SN population. The lack of correlations found in the GRB-SN population between SN brightness and $E_{γ, \rm{iso}}$ for their associated GRBs, across a broad range of 7 orders of magnitude, provides further evidence that the central engine powering the relativistic ejecta is not coupled to the SN powering mechanism in GRB-SN systems. △ Less

Submitted 9 December, 2023; v1 submitted 22 October, 2023; originally announced October 2023.

Comments: 21 pages, 7 Figures, 1 Table, Accepted to ApJ Letters

Multi-band analyses of the bright GRB 230812B and the associated SN2023pel

Authors: T. Hussenot-Desenonges, T. Wouters, N. Guessoum, I. Abdi, A. Abulwfa, C. Adami, J. F. Agüí Fernández, T. Ahumada, V. Aivazyan, D. Akl, S. Anand, C. M. Andrade, S. Antier, S. A. Ata, P. D'Avanzo, Y. A. Azzam, A. Baransky, S. Basa, M. Blazek, P. Bendjoya, S. Beradze, P. Boumis, M. Bremer, R. Brivio, V. Buat , et al. (87 additional authors not shown)

Abstract: GRB~230812B is a bright and relatively nearby ($z =0.36$) long gamma-ray burst (GRB) that has generated significant interest in the community and has thus been observed over the entire electromagnetic spectrum. We report over 80 observations in X-ray, ultraviolet, optical, infrared, and sub-millimeter bands from the GRANDMA (Global Rapid Advanced Network for Multi-messenger Addicts) network of obs… ▽ More GRB~230812B is a bright and relatively nearby ($z =0.36$) long gamma-ray burst (GRB) that has generated significant interest in the community and has thus been observed over the entire electromagnetic spectrum. We report over 80 observations in X-ray, ultraviolet, optical, infrared, and sub-millimeter bands from the GRANDMA (Global Rapid Advanced Network for Multi-messenger Addicts) network of observatories and from observational partners. Adding complementary data from the literature, we then derive essential physical parameters associated with the ejecta and external properties (i.e. the geometry and environment) of the GRB and compare with other analyses of this event. We spectroscopically confirm the presence of an associated supernova, SN2023pel, and we derive a photospheric expansion velocity of v $\sim$ 17$\times10^3$ km s$^{-1}$. We analyze the photometric data first using empirical fits of the flux and then with full Bayesian Inference. We again strongly establish the presence of a supernova in the data, with a maximum (pseudo-)bolometric luminosity of $5.75 \times 10^{42}$ erg/s, at $15.76^{+0.81}_{-1.21}$ days (in the observer frame) after the trigger, with a half-max time width of 22.0 days. We compare these values with those of SN1998bw, SN2006aj, and SN2013dx. Our best-fit model favours a very low density environment ($\log_{10}({n_{\rm ISM}/{\rm cm}^{-3}}) = -2.38^{+1.45}_{-1.60}$) and small values for the jet's core angle $θ_{\rm core} = 1.54^{+1.02}_{-0.81} \ \rm{deg}$ and viewing angle $θ_{\rm obs} = 0.76^{+1.29}_{-0.76} \ \rm{deg}$. GRB 230812B is thus one of the best observed afterglows with a distinctive supernova bump. △ Less

Submitted 17 February, 2024; v1 submitted 22 October, 2023; originally announced October 2023.

Machine-directed gravitational-wave counterpart discovery

Authors: Niharika Sravan, Matthew J. Graham, Michael W. Coughlin, Tomas Ahumada, Shreya Anand

Abstract: Joint observations in electromagnetic and gravitational waves shed light on the physics of objects and surrounding environments with extreme gravity that are otherwise unreachable via siloed observations in each messenger. However, such detections remain challenging due to the rapid and faint nature of counterparts. Protocols for discovery and inference still rely on human experts manually inspect… ▽ More Joint observations in electromagnetic and gravitational waves shed light on the physics of objects and surrounding environments with extreme gravity that are otherwise unreachable via siloed observations in each messenger. However, such detections remain challenging due to the rapid and faint nature of counterparts. Protocols for discovery and inference still rely on human experts manually inspecting survey alert streams and intuiting optimal usage of limited follow-up resources. Strategizing an optimal follow-up program requires adaptive sequential decision-making given evolving light curve data that (i) maximizes a global objective despite incomplete information and (ii) is robust to stochasticity introduced by detectors/observing conditions. Reinforcement learning (RL) approaches allow agents to implicitly learn the physics/detector dynamics and the behavior policy that maximize a designated objective through experience. To demonstrate the utility of such an approach for the kilonova follow-up problem, we train a toy RL agent for the goal of maximizing follow-up photometry for the true kilonova among several contaminant transient light curves. In a simulated environment where the agent learns online, it achieves 3x higher accuracy compared to a random strategy. However, it is surpassed by human agents by up to a factor of 2. This is likely because our hypothesis function (Q that is linear in state-action features) is an insufficient representation of the optimal behavior policy. More complex agents could perform at par or surpass human experts. Agents like these could pave the way for machine-directed software infrastructure to efficiently respond to next generation detectors, for conducting science inference and optimally planning expensive follow-up observations, scalably and with demonstrable performance guarantees. △ Less

Submitted 18 July, 2023; originally announced July 2023.

Comments: Submitted to the Astrophysical Journal; Comments welcome!

Updated observing scenarios and multi-messenger implications for the International Gravitational-wave Network's O4 and O5

Authors: R. Weizmann Kiendrebeogo, Amanda M. Farah, Emily M. Foley, Abigail Gray, Nina Kunert, Anna Puecher, Andrew Toivonen, R. Oliver VandenBerg, Shreya Anand, Tomás Ahumada, Viraj Karambelkar, Michael W. Coughlin, Tim Dietrich, S. Zacharie Kam, Peter T. H. Pang, Leo P. Singer, Niharika Sravan

Abstract: An advanced LIGO and Virgo's third observing run brought another binary neutron star merger (BNS) and the first neutron-star black hole mergers. While no confirmed kilonovae were identified in conjunction with any of these events, continued improvements of analyses surrounding GW170817 allow us to project constraints on the Hubble Constant ($H_0$), the Galactic enrichment from $r$-process nucleosy… ▽ More An advanced LIGO and Virgo's third observing run brought another binary neutron star merger (BNS) and the first neutron-star black hole mergers. While no confirmed kilonovae were identified in conjunction with any of these events, continued improvements of analyses surrounding GW170817 allow us to project constraints on the Hubble Constant ($H_0$), the Galactic enrichment from $r$-process nucleosynthesis, and ultra-dense matter possible from forthcoming events. Here, we describe the expected constraints based on the latest expected event rates from the international gravitational-wave network (IGWN) and analyses of GW170817. We show the expected detection rate of gravitational waves and their counterparts, as well as how sensitive potential constraints are to the observed numbers of counterparts. We intend this analysis as support for the community when creating scientifically driven electromagnetic follow-up proposals. During the next observing run O4, we predict an annual detection rate of electromagnetic counterparts from BNS of $0.43^{+0.58}_{-0.26}$ ($1.97^{+2.68}_{-1.2}$) for the Zwicky Transient Facility (Rubin Observatory). △ Less

Submitted 12 December, 2023; v1 submitted 15 June, 2023; originally announced June 2023.

Journal ref: The Astrophysical Journal , 2023, 958, 158

A data science platform to enable time-domain astronomy

Authors: Michael W. Coughlin, Joshua S. Bloom, Guy Nir, Sarah Antier, Theophile Jegou du Laz, Stéfan van der Walt, Arien Crellin-Quick, Thomas Culino, Dmitry A. Duev, Daniel A. Goldstein, Brian F. Healy, Viraj Karambelkar, Jada Lilleboe, Kyung Min Shin, Leo P. Singer, Tomas Ahumada, Shreya Anand, Eric C. Bellm, Richard Dekany, Matthew J. Graham, Mansi M. Kasliwal, Ivona Kostadinova, R. Weizmann Kiendrebeogo, Shrinivas R. Kulkarni, Sydney Jenkins , et al. (28 additional authors not shown)

Abstract: SkyPortal is an open-source software package designed to efficiently discover interesting transients, manage follow-up, perform characterization, and visualize the results. By enabling fast access to archival and catalog data, cross-matching heterogeneous data streams, and the triggering and monitoring of on-demand observations for further characterization, a SkyPortal-based platform has been oper… ▽ More SkyPortal is an open-source software package designed to efficiently discover interesting transients, manage follow-up, perform characterization, and visualize the results. By enabling fast access to archival and catalog data, cross-matching heterogeneous data streams, and the triggering and monitoring of on-demand observations for further characterization, a SkyPortal-based platform has been operating at scale for 2 yr for the Zwicky Transient Facility Phase II community, with hundreds of users, containing tens of millions of time-domain sources, interacting with dozens of telescopes, and enabling community reporting. While SkyPortal emphasizes rich user experiences (UX) across common frontend workflows, recognizing that scientific inquiry is increasingly performed programmatically, SkyPortal also surfaces an extensive and well-documented API system. From backend and frontend software to data science analysis tools and visualization frameworks, the SkyPortal design emphasizes the re-use and leveraging of best-in-class approaches, with a strong extensibility ethos. For instance, SkyPortal now leverages ChatGPT large-language models (LLMs) to automatically generate and surface source-level human-readable summaries. With the imminent re-start of the next-generation of gravitational wave detectors, SkyPortal now also includes dedicated multi-messenger features addressing the requirements of rapid multi-messenger follow-up: multi-telescope management, team/group organizing interfaces, and cross-matching of multi-messenger data streams with time-domain optical surveys, with interfaces sufficiently intuitive for the newcomers to the field. (abridged) △ Less

Submitted 14 June, 2023; v1 submitted 28 April, 2023; originally announced May 2023.

Comments: Accepted to ApJS

Collapsars as Sites of r-process Nucleosynthesis: Systematic Near-Infrared Follow-up of Type Ic-BL Supernovae

Authors: Shreya Anand, Jennifer Barnes, Sheng Yang, Mansi M. Kasliwal, Michael W. Coughlin, Jesper Sollerman, Kishalay De, Christoffer Fremling, Alessandra Corsi, Anna Y. Q. Ho, Arvind Balasubramanian, Conor Omand, Gokul P. Srinivasaragavan, S. Bradley Cenko, Tomas Ahumada, Igor Andreoni, Aishwarya Dahiwale, Kaustav Kashyap Das, Jacob Jencson, Viraj Karambelkar, Harsh Kumar, Brian D. Metzger, Daniel Perley, Nikhil Sarin, Tassilo Schweyer , et al. (19 additional authors not shown)

Abstract: One of the open questions following the discovery of GW170817 is whether neutron star mergers are the only astrophysical sites capable of producing $r$-process elements. Simulations have shown that 0.01-0.1M$_\odot$ of $r$-process material could be generated in the outflows originating from the accretion disk surrounding the rapidly rotating black hole that forms as a remnant to both neutron star… ▽ More One of the open questions following the discovery of GW170817 is whether neutron star mergers are the only astrophysical sites capable of producing $r$-process elements. Simulations have shown that 0.01-0.1M$_\odot$ of $r$-process material could be generated in the outflows originating from the accretion disk surrounding the rapidly rotating black hole that forms as a remnant to both neutron star mergers and collapsing massive stars associated with long-duration gamma-ray bursts (collapsars). The hallmark signature of $r$-process nucleosynthesis in the binary neutron star merger GW170817 was its long-lasting near-infrared emission, thus motivating a systematic photometric study of the light curves of broadlined stripped-envelope (Ic-BL) supernovae (SNe) associated with collapsars. We present the first systematic study of 25 SNe Ic-BL -- including 18 observed with the Zwicky Transient Facility and 7 from the literature -- in the optical/near-infrared bands to determine what quantity of $r$-process material, if any, is synthesized in these explosions. Using semi-analytic models designed to account for $r$-process production in SNe Ic-BL, we perform light curve fitting to derive constraints on the $r$-process mass for these SNe. We also perform independent light curve fits to models without $r$-process. We find that the $r$-process-free models are a better fit to the light curves of the objects in our sample. Thus we find no compelling evidence of $r$-process enrichment in any of our objects. Further high-cadence infrared photometric studies and nebular spectroscopic analysis would be sensitive to smaller quantities of $r$-process ejecta mass or indicate whether all collapsars are completely devoid of $r$-process nucleosynthesis. △ Less

Submitted 12 February, 2024; v1 submitted 17 February, 2023; originally announced February 2023.

Comments: 37 pages, 8 figures, 4 tables

Journal ref: ApJ, 962, 68 (2024)

A very luminous jet from the disruption of a star by a massive black hole

Authors: Igor Andreoni, Michael W. Coughlin, Daniel A. Perley, Yuhan Yao, Wenbin Lu, S. Bradley Cenko, Harsh Kumar, Shreya Anand, Anna Y. Q. Ho, Mansi M. Kasliwal, Antonio de Ugarte Postigo, Ana Sagues-Carracedo, Steve Schulze, D. Alexander Kann, S. R. Kulkarni, Jesper Sollerman, Nial Tanvir, Armin Rest, Luca Izzo, Jean J. Somalwar, David L. Kaplan, Tomas Ahumada, G. C. Anupama, Katie Auchettl, Sudhanshu Barway , et al. (56 additional authors not shown)

Abstract: Tidal disruption events (TDEs) are bursts of electromagnetic energy released when supermassive black holes (SMBHs) at the centers of galaxies violently disrupt a star that passes too close. TDEs provide a new window to study accretion onto SMBHs; in some rare cases, this accretion leads to launching of a relativistic jet, but the necessary conditions are not fully understood. The best studied jett… ▽ More Tidal disruption events (TDEs) are bursts of electromagnetic energy released when supermassive black holes (SMBHs) at the centers of galaxies violently disrupt a star that passes too close. TDEs provide a new window to study accretion onto SMBHs; in some rare cases, this accretion leads to launching of a relativistic jet, but the necessary conditions are not fully understood. The best studied jetted TDE to date is Swift J1644+57, which was discovered in gamma-rays, but was too obscured by dust to be seen at optical wavelengths. Here we report the optical discovery of AT2022cmc, a rapidly fading source at cosmological distance (redshift z=1.19325) whose unique lightcurve transitioned into a luminous plateau within days. Observations of a bright counterpart at other wavelengths, including X-rays, sub-millimeter, and radio, supports the interpretation of AT2022cmc as a jetted TDE containing a synchrotron "afterglow", likely launched by a SMBH with spin $a \gtrsim 0.3$. Using 4 years of Zwicky Transient Facility (ZTF) survey data, we calculate a rate of $0.02 ^{+ 0.04 }_{- 0.01 }$ Gpc$^{-3}$ yr$^{-1}$ for on-axis jetted TDEs based on the luminous, fast-fading red component, thus providing a measurement complementary to the rates derived from X-ray and radio observations. Correcting for the beaming angle effects, this rate confirms that about 1% of TDEs have relativistic jets. Optical surveys can use AT2022cmc as a prototype to unveil a population of jetted TDEs. △ Less

Submitted 29 November, 2022; originally announced November 2022.

Comments: Accepted for publication in Nature

Deep Synoptic Array science I: discovery of the host galaxy of FRB 20220912A

Authors: Vikram Ravi, Morgan Catha, Ge Chen, Liam Connor, Jakob T. Faber, James W. Lamb, Gregg Hallinan, Charlie Harnach, Greg Hellbourg, Rick Hobbs, David Hodge, Mark Hodges, Casey Law, Paul Rasmussen, Kritti Sharma, Myles B. Sherman, Jun Shi, Dana Simard, Reynier Squillace, Sander Weinreb, David P. Woody, Nitika Yadlapalli, Tomas Ahumada, Dillon Dong, Christoffer Fremling , et al. (3 additional authors not shown)

Abstract: We report the detection and interferometric localization of the repeating fast radio burst (FRB) source FRB 20220912A during commissioning observations with the Deep Synoptic Array (DSA-110). Two bursts were detected from FRB 20220912A, one each on 2022 October 18 and 2022 October 25. The best-fit position is (R.A. J2000, decl. J2000) = (23:09:04.9, +48:42:25.4), with a 90% confidence error ellips… ▽ More We report the detection and interferometric localization of the repeating fast radio burst (FRB) source FRB 20220912A during commissioning observations with the Deep Synoptic Array (DSA-110). Two bursts were detected from FRB 20220912A, one each on 2022 October 18 and 2022 October 25. The best-fit position is (R.A. J2000, decl. J2000) = (23:09:04.9, +48:42:25.4), with a 90% confidence error ellipse of $\pm2$ arcsec and $\pm1$ arcsec in right ascension and declination respectively. The two bursts have disparate polarization properties and temporal profiles. We find a Faraday rotation measure that is consistent with the low value of $+0.6$ rad m$^{-2}$ reported by CHIME/FRB. The DSA-110 localization overlaps with the galaxy PSO J347.2702+48.7066 at a redshift $z=0.0771$, which we identify as the likely host. PSO J347.2702$+$48.7066 has a stellar mass of approximately $10^{10}M_{\odot}$, modest internal dust extinction, and a star-formation rate likely in excess of $0.1\,M_{\odot}$ yr$^{-1}$. The host-galaxy contribution to the dispersion measure is likely $\lesssim50$ pc cm$^{-3}$. The FRB 20220912A source is therefore likely viewed along a tenuous plasma column through the host galaxy. △ Less

Submitted 16 November, 2022; originally announced November 2022.

Comments: 10 pages, 7 figures, 2 tables, submitted to AAS Journals

The discovery and characterization of a kilometre sized asteroid inside the orbit of Venus

Authors: Bryce T. Bolin, T. Ahumada, P. van Dokkum, C. Fremling, M. Granvik, K. K. Hardegree-Ullman, Y. Harikane, J. N. Purdum, E. Serabyn, J. Southworth, C. Zhai

Abstract: Near-Earth asteroid population models predict the existence of bodies located inside the orbit of Venus. Despite searches up to the end of 2019, none had been found. We report discovery and follow-up observations of (594913) 'Ayló'chaxnim, an asteroid with an orbit entirely interior to Venus. (594913) 'Ayló'chaxnim has an aphelion distance of ~0.65 au, is ~2 km in diameter and is red in colour. Th… ▽ More Near-Earth asteroid population models predict the existence of bodies located inside the orbit of Venus. Despite searches up to the end of 2019, none had been found. We report discovery and follow-up observations of (594913) 'Ayló'chaxnim, an asteroid with an orbit entirely interior to Venus. (594913) 'Ayló'chaxnim has an aphelion distance of ~0.65 au, is ~2 km in diameter and is red in colour. The detection of such a large asteroid inside the orbit of Venus is surprising given their rarity according to near-Earth asteroid population models. As the first officially numbered and named asteroid located entirely within the orbit of Venus, we propose that the class of interior to Venus asteroids be referred to as 'Ayló'chaxnim asteroids. △ Less

Submitted 15 August, 2022; originally announced August 2022.

Comments: 7 pages, 8 figures, 1 table. Published in Monthly Notices of the Royal Astronomical Society: Letters (2022). arXiv admin note: text overlap with arXiv:2009.04125

From Data to Software to Science with the Rubin Observatory LSST

Authors: Katelyn Breivik, Andrew J. Connolly, K. E. Saavik Ford, Mario Jurić, Rachel Mandelbaum, Adam A. Miller, Dara Norman, Knut Olsen, William O'Mullane, Adrian Price-Whelan, Timothy Sacco, J. L. Sokoloski, Ashley Villar, Viviana Acquaviva, Tomas Ahumada, Yusra AlSayyad, Catarina S. Alves, Igor Andreoni, Timo Anguita, Henry J. Best, Federica B. Bianco, Rosaria Bonito, Andrew Bradshaw, Colin J. Burke, Andresa Rodrigues de Campos , et al. (75 additional authors not shown)

Abstract: The Vera C. Rubin Observatory Legacy Survey of Space and Time (LSST) dataset will dramatically alter our understanding of the Universe, from the origins of the Solar System to the nature of dark matter and dark energy. Much of this research will depend on the existence of robust, tested, and scalable algorithms, software, and services. Identifying and developing such tools ahead of time has the po… ▽ More The Vera C. Rubin Observatory Legacy Survey of Space and Time (LSST) dataset will dramatically alter our understanding of the Universe, from the origins of the Solar System to the nature of dark matter and dark energy. Much of this research will depend on the existence of robust, tested, and scalable algorithms, software, and services. Identifying and developing such tools ahead of time has the potential to significantly accelerate the delivery of early science from LSST. Developing these collaboratively, and making them broadly available, can enable more inclusive and equitable collaboration on LSST science. To facilitate such opportunities, a community workshop entitled "From Data to Software to Science with the Rubin Observatory LSST" was organized by the LSST Interdisciplinary Network for Collaboration and Computing (LINCC) and partners, and held at the Flatiron Institute in New York, March 28-30th 2022. The workshop included over 50 in-person attendees invited from over 300 applications. It identified seven key software areas of need: (i) scalable cross-matching and distributed joining of catalogs, (ii) robust photometric redshift determination, (iii) software for determination of selection functions, (iv) frameworks for scalable time-series analyses, (v) services for image access and reprocessing at scale, (vi) object image access (cutouts) and analysis at scale, and (vii) scalable job execution systems. This white paper summarizes the discussions of this workshop. It considers the motivating science use cases, identified cross-cutting algorithms, software, and services, their high-level technical specifications, and the principles of inclusive collaborations needed to develop them. We provide it as a useful roadmap of needs, as well as to spur action and collaboration between groups and individuals looking to develop reusable software for early LSST science. △ Less

Submitted 4 August, 2022; originally announced August 2022.

Comments: White paper from "From Data to Software to Science with the Rubin Observatory LSST" workshop

The long-active afterglow of GRB 210204A: Detection of the most delayed flares in a Gamma-Ray Burst

Authors: Harsh Kumar, Rahul Gupta, Divita Saraogi, Tomás Ahumada, Igor Andreoni, G. C. Anupama, Amar Aryan, Sudhanshu Barway, Varun Bhalerao, Poonam Chandra, Michael W. Coughlin, Dimple, Anirban Dutta, Ankur Ghosh, Anna Y. Q. Ho, E. C. Kool, Amit Kumar, Michael S. Medford, Kuntal Misra, Shashi B. Pandey, Daniel A. Perley, Reed Riddle, Amit Kumar Ror, Jason M. Setiadi, Yuhan Yao

Abstract: We present results from extensive broadband follow-up of GRB 210204A over the period of thirty days. We detect optical flares in the afterglow at 7.6 x 10^5 s and 1.1 x 10^6 s after the burst: the most delayed flaring ever detected in a GRB afterglow. At the source redshift of 0.876, the rest-frame delay is 5.8 x 10^5 s (6.71 d). We investigate possible causes for this flaring and conclude that th… ▽ More We present results from extensive broadband follow-up of GRB 210204A over the period of thirty days. We detect optical flares in the afterglow at 7.6 x 10^5 s and 1.1 x 10^6 s after the burst: the most delayed flaring ever detected in a GRB afterglow. At the source redshift of 0.876, the rest-frame delay is 5.8 x 10^5 s (6.71 d). We investigate possible causes for this flaring and conclude that the most likely cause is a refreshed shock in the jet. The prompt emission of the GRB is within the range of typical long bursts: it shows three disjoint emission episodes, which all follow the typical GRB correlations. This suggests that GRB 210204A might not have any special properties that caused late-time flaring, and the lack of such detections for other afterglows might be resulting from the paucity of late-time observations. Systematic late-time follow-up of a larger sample of GRBs can shed more light on such afterglow behaviour. Further analysis of the GRB 210204A shows that the late time bump in the light curve is highly unlikely due to underlying SNe at redshift (z) = 0.876 and is more likely due to the late time flaring activity. The cause of this variability is not clearly quantifiable due to the lack of multi-band data at late time constraints by the bad weather conditions. The flare of GRB 210204A is the latest flare detected to date. △ Less

Submitted 15 April, 2022; originally announced April 2022.

Comments: 18 pages, 12 figures

Neutrino follow-up with the Zwicky Transient Facility: Results from the first 24 campaigns

Authors: Robert Stein, Simeon Reusch, Anna Franckowiak, Marek Kowalski, Jannis Necker, Sven Weimann, Mansi M. Kasliwal, Jesper Sollerman, Tomas Ahumada, Pau Amaro-Seoane, Shreya Anand, Igor Andreoni, Eric C. Bellm, Joshua S. Bloom, Michael Coughlin, Kishalay De, Christoffer Fremling, Suvi Gezari, Matthew Graham, Steven L. Groom, George Helou, David L. Kaplan, Viraj Karambelkar, Albert K. H. Kong, Erik C. Kool , et al. (11 additional authors not shown)

Abstract: The Zwicky Transient Facility (ZTF) performs a systematic neutrino follow-up program, searching for optical counterparts to high-energy neutrinos with dedicated Target-of-Opportunity (ToO) observations. Since first light in March 2018, ZTF has taken prompt observations for 24 high-quality neutrino alerts from the IceCube Neutrino Observatory, with a median latency of 12.2 hours from initial neutri… ▽ More The Zwicky Transient Facility (ZTF) performs a systematic neutrino follow-up program, searching for optical counterparts to high-energy neutrinos with dedicated Target-of-Opportunity (ToO) observations. Since first light in March 2018, ZTF has taken prompt observations for 24 high-quality neutrino alerts from the IceCube Neutrino Observatory, with a median latency of 12.2 hours from initial neutrino detection. From two of these campaigns, we have already reported tidal disruption event (TDE) AT 2019dsg and likely TDE AT 2019fdr as probable counterparts, suggesting that TDEs contribute >7.8% of the astrophysical neutrino flux. We here present the full results of our program through to December 2021. No additional candidate neutrino sources were identified by our program, allowing us to place the first constraints on the underlying optical luminosity function of astrophysical neutrino sources. Transients with optical absolutes magnitudes brighter that $-21$ can contribute no more than 87% of the total, while transients brighter than $-22$ can contribute no more than 58% of the total, neglecting the effect of extinction and assuming they follow the star formation rate. These are the first observational constraints on the neutrino emission of bright populations such as superluminous supernovae. None of the neutrinos were coincident with bright optical AGN flares comparable to that observed for TXS 0506+056/IC170922A, with such optical blazar flares producing no more than 26% of the total neutrino flux. We highlight the outlook for electromagnetic neutrino follow-up programs, including the expected potential for the Rubin Observatory. △ Less

Submitted 4 April, 2024; v1 submitted 31 March, 2022; originally announced March 2022.

Journal ref: MNRAS, 521, 4, 2023, 5046-5063

In search of short gamma-ray burst optical counterpart with the Zwicky Transient Facility

Authors: Tomás Ahumada, Shreya Anand, Michael W. Coughlin, Igor Andreoni, Erik C. Kool, Harsh Kumar, Simeon Reusch, Ana Sagués-Carracedo, Robert Stein, S. Bradley Cenko, Mansi M. Kasliwal, Leo P. Singer, Rachel Dunwoody, Joseph Mangan, Varun Bhalerao, Mattia Bulla, Eric Burns, Matthew J. Graham, David L. Kaplan, Daniel Perley, Mouza Almualla, Joshua S. Bloom, Virginia Cunningham, Kishalay De, Pradip Gatkine , et al. (24 additional authors not shown)

Abstract: The Fermi Gamma-ray Burst Monitor (GBM) triggers on-board in response to $\sim$ 40 short gamma-ray bursts (SGRBs) per year; however, their large localization regions have made the search for optical counterparts a challenging endeavour. We have developed and executed an extensive program with the wide field of view of the Zwicky Transient Facility (ZTF) camera, mounted on the Palomar 48 inch Oschi… ▽ More The Fermi Gamma-ray Burst Monitor (GBM) triggers on-board in response to $\sim$ 40 short gamma-ray bursts (SGRBs) per year; however, their large localization regions have made the search for optical counterparts a challenging endeavour. We have developed and executed an extensive program with the wide field of view of the Zwicky Transient Facility (ZTF) camera, mounted on the Palomar 48 inch Oschin telescope (P48), to perform target-of-opportunity (ToO) observations on 10 Fermi-GBM SGRBs during 2018 and 2020-2021. Bridging the large sky areas with small field of view optical telescopes in order to track the evolution of potential candidates, we look for the elusive SGRB afterglows and kilonovae (KNe) associated with these high-energy events. No counterpart has yet been found, even though more than 10 ground based telescopes, part of the Global Relay of Observatories Watching Transients Happen (GROWTH) network, have taken part in these efforts. The candidate selection procedure and the follow-up strategy have shown that ZTF is an efficient instrument for searching for poorly localized SGRBs, retrieving a reasonable number of candidates to follow-up and showing promising capabilities as the community approaches the multi-messenger era. Based on the median limiting magnitude of ZTF, our searches would have been able to retrieve a GW170817-like event up to $\sim$ 200 Mpc and SGRB afterglows to z = 0.16 or 0.4, depending on the assumed underlying energy model. Future ToOs will expand the horizon to z = 0.2 and 0.7 respectively. △ Less

Submitted 22 March, 2022; originally announced March 2022.

Candidate Tidal Disruption Event AT2019fdr Coincident with a High-Energy Neutrino

Authors: Simeon Reusch, Robert Stein, Marek Kowalski, Sjoert van Velzen, Anna Franckowiak, Cecilia Lunardini, Kohta Murase, Walter Winter, James C. A. Miller-Jones, Mansi M. Kasliwal, Marat Gilfanov, Simone Garrappa, Vaidehi S. Paliya, Tomas Ahumada, Shreya Anand, Cristina Barbarino, Eric C. Bellm, Valery Brinnel, Sara Buson, S. Bradley Cenko, Michael W. Coughlin, Kishalay De, Richard Dekany, Sara Frederick, Avishay Gal-Yam , et al. (24 additional authors not shown)

Abstract: The origins of the high-energy cosmic neutrino flux remain largely unknown. Recently, one high-energy neutrino was associated with a tidal disruption event (TDE). Here we present AT2019fdr, an exceptionally luminous TDE candidate, coincident with another high-energy neutrino. Our observations, including a bright dust echo and soft late-time X-ray emission, further support a TDE origin of this flar… ▽ More The origins of the high-energy cosmic neutrino flux remain largely unknown. Recently, one high-energy neutrino was associated with a tidal disruption event (TDE). Here we present AT2019fdr, an exceptionally luminous TDE candidate, coincident with another high-energy neutrino. Our observations, including a bright dust echo and soft late-time X-ray emission, further support a TDE origin of this flare. The probability of finding two such bright events by chance is just 0.034%. We evaluate several models for neutrino production and show that AT2019fdr is capable of producing the observed high-energy neutrino, reinforcing the case for TDEs as neutrino sources. △ Less

Submitted 10 June, 2022; v1 submitted 17 November, 2021; originally announced November 2021.

Comments: 20 pages, 6 figures, 6 tables

Journal ref: Phys. Rev. Lett. 128, 221101 (2022)

Inferring kilonova population properties with a hierarchical Bayesian framework I : Non-detection methodology and single-event analyses

Authors: Siddharth R. Mohite, Priyadarshini Rajkumar, Shreya Anand, David L. Kaplan, Michael W. Coughlin, Ana Sagués-Carracedo, Muhammed Saleem, Jolien Creighton, Patrick R. Brady, Tomás Ahumada, Mouza Almualla, Igor Andreoni, Mattia Bulla, Matthew J. Graham, Mansi M. Kasliwal, Stephen Kaye, Russ R. Laher, Kyung Min Shin, David L. Shupe, Leo P. Singer

Abstract: We present ${\tt nimbus}$ : a hierarchical Bayesian framework to infer the intrinsic luminosity parameters of kilonovae (KNe) associated with gravitational-wave (GW) events, based purely on non-detections. This framework makes use of GW 3-D distance information and electromagnetic upper limits from multiple surveys for multiple events, and self-consistently accounts for finite sky-coverage and pro… ▽ More We present ${\tt nimbus}$ : a hierarchical Bayesian framework to infer the intrinsic luminosity parameters of kilonovae (KNe) associated with gravitational-wave (GW) events, based purely on non-detections. This framework makes use of GW 3-D distance information and electromagnetic upper limits from multiple surveys for multiple events, and self-consistently accounts for finite sky-coverage and probability of astrophysical origin. The framework is agnostic to the brightness evolution assumed and can account for multiple electromagnetic passbands simultaneously. Our analyses highlight the importance of accounting for model selection effects, especially in the context of non-detections. We show our methodology using a simple, two-parameter linear brightness model, taking the follow-up of GW190425 with the Zwicky Transient Facility (ZTF) as a single-event test case for two different prior choices of model parameters -- (i) uniform/uninformative priors and (ii) astrophysical priors based on surrogate models of Monte Carlo radiative transfer simulations of KNe. We present results under the assumption that the KN is within the searched region to demonstrate functionality and the importance of prior choice. Our results show consistency with ${\tt simsurvey}$ -- an astronomical survey simulation tool used previously in the literature to constrain the population of KNe. While our results based on uniform priors strongly constrain the parameter space, those based on astrophysical priors are largely uninformative, highlighting the need for deeper constraints. Future studies with multiple events having electromagnetic follow-up from multiple surveys should make it possible to constrain the KN population further. △ Less

Submitted 15 November, 2021; v1 submitted 15 July, 2021; originally announced July 2021.

Comments: 18 pages, 5 figures, accepted for publication in ApJ

Discovery and confirmation of the shortest gamma ray burst from a collapsar

Authors: Tomas Ahumada, Leo P. Singer, Shreya Anand, Michael W. Coughlin, Mansi M. Kasliwal, Geoffrey Ryan, Igor Andreoni, S. Bradley Cenko, Christoffer Fremling, Harsh Kumar, Peter T. H. Pang, Eric Burns, Virginia Cunningham, Simone Dichiara, Tim Dietrich, Dmitry S. Svinkin, Mouza Almualla, Alberto J. Castro-Tirado, Kishalay De, Rachel Dunwoody, Pradip Gatkine, Erica Hammerstein, Shabnam Iyyani, Joseph Mangan, Dan Perley , et al. (32 additional authors not shown)

Abstract: Gamma-ray bursts (GRBs) are among the brightest and most energetic events in the universe. The duration and hardness distribution of GRBs has two clusters, now understood to reflect (at least) two different progenitors. Short-hard GRBs (SGRBs; T90 <2 s) arise from compact binary mergers, while long-soft GRBs (LGRBs; T90 >2 s) have been attributed to the collapse of peculiar massive stars (collapsa… ▽ More Gamma-ray bursts (GRBs) are among the brightest and most energetic events in the universe. The duration and hardness distribution of GRBs has two clusters, now understood to reflect (at least) two different progenitors. Short-hard GRBs (SGRBs; T90 <2 s) arise from compact binary mergers, while long-soft GRBs (LGRBs; T90 >2 s) have been attributed to the collapse of peculiar massive stars (collapsars). The discovery of SN 1998bw/GRB 980425 marked the first association of a LGRB with a collapsar and AT 2017gfo/GRB 170817A/GW170817 marked the first association of a SGRB with a binary neutron star merger, producing also gravitational wave (GW). Here, we present the discovery of ZTF20abwysqy (AT2020scz), a fast-fading optical transient in the Fermi Satellite and the InterPlanetary Network (IPN) localization regions of GRB 200826A; X-ray and radio emission further confirm that this is the afterglow. Follow-up imaging (at rest-frame 16.5 days) reveals excess emission above the afterglow that cannot be explained as an underlying kilonova (KN), but is consistent with being the supernova (SN). Despite the GRB duration being short (rest-frame T90 of 0.65 s), our panchromatic follow-up data confirms a collapsar origin. GRB 200826A is the shortest LGRB found with an associated collapsar; it appears to sit on the brink between a successful and a failed collapsar. Our discovery is consistent with the hypothesis that most collapsars fail to produce ultra-relativistic jets. △ Less

Submitted 13 May, 2021; v1 submitted 11 May, 2021; originally announced May 2021.

Comments: submitted to Nature Astronomy

Fast-transient Searches in Real Time with ZTFReST: Identification of Three Optically-discovered Gamma-ray Burst Afterglows and New Constraints on the Kilonova Rate

Authors: Igor Andreoni, Michael W. Coughlin, Erik C. Kool, Mansi M. Kasliwal, Harsh Kumar, Varun Bhalerao, Ana Sagués Carracedo, Anna Y. Q. Ho, Peter T. H. Pang, Divita Saraogi, Kritti Sharma, Vedant Shenoy, Eric Burns, Tomás Ahumada, Shreya Anand, Leo P. Singer, Daniel A. Perley, Kishalay De, U. C. Fremling, Eric C. Bellm, Mattia Bulla, Arien Crellin-Quick, Tim Dietrich, Andrew Drake, Dmitry A. Duev , et al. (10 additional authors not shown)

Abstract: While optical surveys regularly discover slow transients like supernovae on their own, the most common way to discover extragalactic fast transients, fading away in a few nights, is via follow-up observations of gamma-ray burst and gravitational-wave triggers. However, wide-field surveys have the potential to also identify rapidly fading transients independently of such external triggers. The volu… ▽ More While optical surveys regularly discover slow transients like supernovae on their own, the most common way to discover extragalactic fast transients, fading away in a few nights, is via follow-up observations of gamma-ray burst and gravitational-wave triggers. However, wide-field surveys have the potential to also identify rapidly fading transients independently of such external triggers. The volumetric survey speed of the Zwicky Transient Facility (ZTF) makes it sensitive to faint and fast-fading objects as kilonovae, the optical counterparts to binary neutron stars and neutron star-black hole mergers, out to almost 200Mpc. We introduce an open-source software infrastructure, the ZTF REaltime Search and Triggering, ZTFReST, designed to identify kilonovae and fast optical transients in ZTF data. Using the ZTF alert stream combined with forced photometry, we have implemented automated candidate ranking based on their photometric evolution and fitting to kilonova models. Automated triggering of follow-up systems, such as Las Cumbres Observatory, has also been implemented. In 13 months of science validation, we found several extragalactic fast transients independent of any external trigger (though some counterparts were identified later), including at least one supernova with post-shock cooling emission, two known afterglows with an associated gamma-ray burst, two known afterglows without any known gamma-ray counterpart, and three new fast-declining sources (ZTF20abtxwfx, ZTF20acozryr, and ZTF21aagwbjr) that are likely associated with GRB200817A, GRB201103B, and GRB210204A. However, we have not found any objects which appear to be kilonovae; therefore, we constrain the rate of GW170817-like kilonovae to $R < 900$Gpc$^{-3}$yr$^{-1}$. A framework such as ZTFReST could become a prime tool for kilonova and fast transient discovery with the Vera C. Rubin Observatory. △ Less

Submitted 13 April, 2021; originally announced April 2021.

Comments: Submitted for publication

Time-series and Phasecurve Photometry of Episodically-Active Asteroid (6478) Gault in a Quiescent State Using APO, GROWTH, P200 and ZTF

Authors: Josiah N. Purdum, Zhong-Yi Lin, Bryce T. Bolin, Kritti Sharma, Philip I. Choi, Varun Bhalerao, Harsh Kumar, Robert Quimby, Joannes C. Van Roestel, Chengxing Zhai, Yanga R. Fernandez, Josef Hanuš, Carey M. Lisse, Dennis Bodewits, Christoffer Fremling, Nathan Ryan Golovich, Chen-Yen Hsu, Wing-Huen Ip, Chow-Choong Ngeow, Navtej S. Saini, Michael Shao, Yuhan Yao, Tomás Ahumada, Shreya Anand, Igor Andreoni , et al. (27 additional authors not shown)

Abstract: We observed Episodically Active Asteroid (6478) Gault in 2020 with multiple telescopes in Asia and North America and have found that it is no longer active after its recent outbursts at the end of 2018 and start of 2019. The inactivity during this apparation allowed us to measure the absolute magnitude of Gault of H_r = 14.63 +/- 0.02, G_r = 0.21 +/- 0.02 from our secular phasecurve observations.… ▽ More We observed Episodically Active Asteroid (6478) Gault in 2020 with multiple telescopes in Asia and North America and have found that it is no longer active after its recent outbursts at the end of 2018 and start of 2019. The inactivity during this apparation allowed us to measure the absolute magnitude of Gault of H_r = 14.63 +/- 0.02, G_r = 0.21 +/- 0.02 from our secular phasecurve observations. In addition, we were able to constrain Gault's rotation period using time-series photometric lightcurves taken over 17 hours on multiple days in 2020 August, September and October. The photometric lightcurves have a repeating $\lesssim$0.05 magnitude feature suggesting that (6478) Gault has a rotation period of ~2.5 hours and may have a semi-spherical or top-like shape, much like Near-Earth Asteroids Ryugu and Bennu. The rotation period of ~2.5 hours is near to the expected critical rotation period for an asteroid with the physical properties of (6478) Gault suggesting that its activity observed over multiple epochs is due to surface mass shedding from its fast rotation spun up by the Yarkovsky-O'Keefe-Radzievskii-Paddack effect. △ Less

Submitted 27 March, 2021; v1 submitted 25 February, 2021; originally announced February 2021.

Comments: 23 pages, 18 figures; Accepted by ApJ Letters

Optimizing Serendipitous Detections of Kilonovae: Cadence and Filter Selection

Authors: Mouza Almualla, Shreya Anand, Michael W. Coughlin, Tim Dietrich, Nidhal Guessoum, Ana Sagués Carracedo, Tomás Ahumada, Igor Andreoni, Sarah Antier, Eric C. Bellm, Mattia Bulla, Leo P. Singer

Abstract: The rise of multi-messenger astronomy has brought with it the need to exploit all available data streams and learn more about the astrophysical objects that fall within its breadth. One possible avenue is the search for serendipitous optical/near-infrared counterparts of gamma-ray bursts (GRBs) and gravitational-wave (GW) signals, known as kilonovae. With surveys such as the Zwicky Transient Facil… ▽ More The rise of multi-messenger astronomy has brought with it the need to exploit all available data streams and learn more about the astrophysical objects that fall within its breadth. One possible avenue is the search for serendipitous optical/near-infrared counterparts of gamma-ray bursts (GRBs) and gravitational-wave (GW) signals, known as kilonovae. With surveys such as the Zwicky Transient Facility (ZTF), which observes the sky with a cadence of ~ three days, the existing counterpart locations are likely to be observed; however, due to the significant amount of sky to explore, it is difficult to search for these fast-evolving candidates. Thus, it is beneficial to optimize the survey cadence for realtime kilonova identification and enable further photometric and spectroscopic observations. We explore how the cadence of wide field-of-view surveys like ZTF can be improved to facilitate such identifications. We show that with improved observational choices, e.g., the adoption of three epochs per night on a ~ nightly basis, and the prioritization of redder photometric bands, detection efficiencies improve by about a factor of two relative to the nominal cadence. We also provide realistic hypothetical constraints on the kilonova rate as a form of comparison between strategies, assuming that no kilonovae are detected throughout the long-term execution of the respective observing plan. These results demonstrate how an optimal use of ZTF increases the likelihood of kilonova discovery independent of GWs or GRBs, thereby allowing for a sensitive search with less interruption of its nominal cadence through Target of Opportunity programs. △ Less

Submitted 16 April, 2021; v1 submitted 20 November, 2020; originally announced November 2020.

Comments: 10 pages, 8 figures

Optical follow-up of the neutron star-black hole mergers S200105ae and S200115j

Authors: Shreya Anand, Michael W. Coughlin, Mansi M. Kasliwal, Mattia Bulla, Tomás Ahumada, Ana Sagués Carracedo, Mouza Almualla, Igor Andreoni, Robert Stein, Francois Foucart, Leo P. Singer, Jesper Sollerman, Eric C. Bellm, Bryce Bolin, M. D. Caballero-García, Alberto J. Castro-Tirado, S. Bradley Cenko, Kishalay De, Richard G. Dekany, Dmitry A. Duev, Michael Feeney, Christoffer Fremling, Daniel A. Goldstein, V. Zach Golkhou, Matthew J. Graham , et al. (24 additional authors not shown)

Abstract: LIGO and Virgo's third observing run (O3) revealed the first neutron star-black hole (NSBH) merger candidates in gravitational waves. These events are predicted to synthesize r-process elements creating optical/near-IR "kilonova" (KN) emission. The joint gravitational-wave (GW) and electromagnetic detection of an NSBH merger could be used to constrain the equation of state of dense nuclear matter,… ▽ More LIGO and Virgo's third observing run (O3) revealed the first neutron star-black hole (NSBH) merger candidates in gravitational waves. These events are predicted to synthesize r-process elements creating optical/near-IR "kilonova" (KN) emission. The joint gravitational-wave (GW) and electromagnetic detection of an NSBH merger could be used to constrain the equation of state of dense nuclear matter, and independently measure the local expansion rate of the universe. Here, we present the optical follow-up and analysis of two of the only three high-significance NSBH merger candidates detected to date, S200105ae and S200115j, with the Zwicky Transient Facility (ZTF). ZTF observed $\sim$\,48\% of S200105ae and $\sim$\,22\% of S200115j's localization probabilities, with observations sensitive to KNe brighter than $-$17.5\,mag fading at 0.5\,mag/day in g- and r-bands; extensive searches and systematic follow-up of candidates did not yield a viable counterpart. We present state-of-the-art KN models tailored to NSBH systems that place constraints on the ejecta properties of these NSBH mergers. We show that with depths of $\rm m_{\rm AB}\approx 22$ mag, attainable in meter-class, wide field-of-view survey instruments, strong constraints on ejecta mass are possible, with the potential to rule out low mass ratios, high BH spins, and large neutron star radii. △ Less

Submitted 14 September, 2020; originally announced September 2020.

Preliminary estimates of the Zwicky Transient Facility 'Ayló'chaxnim asteroid population completeness

Authors: B. T. Bolin, T. Ahumada, P. van Dokkum, C. Fremling, K. K. Hardegree-Ullman, J. N. Purdum, E. Serabyn, J. Southworth

Abstract: Near-Earth asteroids (NEAs) are organized into five main classes: Amor, Apollo, Aten, Atira and 'Ayló'chaxnim. Asteroids belonging to the 'Ayló'chaxnim class are located entirely within the orbit of Venus making them difficult to detect by ground-based observatories. The first-known asteroid of this class, (594913) 'Ayló'chaxnim, was discovered by the Zwicky Transient Facility (ZTF) in 2020 Januar… ▽ More Near-Earth asteroids (NEAs) are organized into five main classes: Amor, Apollo, Aten, Atira and 'Ayló'chaxnim. Asteroids belonging to the 'Ayló'chaxnim class are located entirely within the orbit of Venus making them difficult to detect by ground-based observatories. The first-known asteroid of this class, (594913) 'Ayló'chaxnim, was discovered by the Zwicky Transient Facility (ZTF) in 2020 January during a twilight search for asteroids at small solar elongations that ran between September 2019 and January 2020. Due to its large diameter of $\sim$2 km, the discovery of (594913) 'Ayló'chaxnim is surprising because contemporary NEA population models predict a scarcity of asteroids of this size located inside the orbit of Venus. To compare the discovery of (594913) 'Ayló'chaxnim by ZTF with the predictions of NEA population models, we estimated the ZTF survey completeness at detecting 'Ayló'chaxnim asteroids and the number of 'Ayló'chaxnim asteroids expected to have been discovered by simulating observations of synthetic 'Ayló'chaxnim asteroids. We find that the 'Ayló'chaxnim population completeness of the survey is $\sim$18$\%$ and there is only a 5$\%$ probability that a single 'Ayló'chaxnim asteroid would have been discovered. Given the small chance for (594913) 'Ayló'chaxnim to have been discovered, its presence is either a statistical fluke or it implies that asteroid population models may need to be revised. △ Less

Submitted 19 January, 2023; v1 submitted 9 September, 2020; originally announced September 2020.

Comments: 9 pages, 2 figures. Accepted for publication in Icarus

Characterization of Temporarily-Captured Minimoon 2020 CD$_3$ by Keck Time-resolved Spectrophotometry

Authors: Bryce T. Bolin, Christoffer Fremling, Timothy R. Holt, Matthew J. Hankins, Tomás Ahumada, Shreya Anand, Varun Bhalerao, Kevin B. Burdge, Chris M. Copperwheat, Michael Coughlin, Kunal P. Deshmukh, Kishalay De, Mansi M. Kasliwal, Alessandro Morbidelli, Josiah N. Purdum, Robert Quimby, Dennis Bodewits, Chan-Kao Chang, Wing-Huen Ip, Chen-Yen Hsu, Russ R. Laher, Zhong-Yi Lin, Carey M. Lisse, Frank J. Masci, Chow-Choong Ngeow , et al. (20 additional authors not shown)

Abstract: We present time-resolved visible spectrophotometry of minimoon 2020 CD$_3$, the second asteroid known to become temporarily captured by the Earth-Moon system's gravitational field. The spectrophotometry was taken with Keck I/LRIS between wavelengths 434 nm and 912 nm in $B$, $g$, $V$, $R$, $I$ and RG850 filters as it was leaving the Earth-Moon system on 2020 March 23 UTC. The spectrophotometry of… ▽ More We present time-resolved visible spectrophotometry of minimoon 2020 CD$_3$, the second asteroid known to become temporarily captured by the Earth-Moon system's gravitational field. The spectrophotometry was taken with Keck I/LRIS between wavelengths 434 nm and 912 nm in $B$, $g$, $V$, $R$, $I$ and RG850 filters as it was leaving the Earth-Moon system on 2020 March 23 UTC. The spectrophotometry of 2020 CD$_3$ most closely resembles the spectra of V-type asteroids and some Lunar rock samples with a reddish slope of ~18$\%$/100 nm between 434 nm and 761 nm corresponding to colors of $g$-$r$ = 0.62$\pm$0.08, $r$-$i$ = 0.21 $\pm$ 0.06 and an absorption band at ~900 nm corresponding to $i$-$z$ = -0.54$\pm$0.10. Combining our measured 31.9$\pm$0.1 absolute magnitude with an albedo of 0.35 typical for V-type asteroids, we determine 2020 CD$_3$'s diameter to be ~0.9$\pm$0.1 m making it the first minimoon and one of the smallest asteroids to be spectrally studied. We use our time-series photometry to detect periodic lightcurve variations with a $<$10$^{-4}$ false alarm probability corresponding to a lightcurve period of ~573 s and a lightcurve amplitude of ~1 mag implying 2020 CD$_3$ possesses a $b/a$ axial ratio of ~2.5. In addition, we extend the observational arc of 2020 CD$_3$ to 37 days between 2020 February 15 UTC and 2020 March 23 UTC. From the improved orbital solution for 2020 CD$_3$, we estimate its likely duration of its capture to be ~2 y, and we measure the non-gravitation perturbation on its orbit due to radiation pressure with an area-to-mass ratio of 6.9$\pm$2.4$\times$10$^{-4}$ m$^2$/kg implying a density of 2.3$\pm$0.8 g/cm$^3$, broadly compatible with the densities of other meter-scale asteroids and Lunar rock. We searched for pre-discovery detections of 2020 CD$_3$ in the ZTF archive as far back as 2018 October, but were unable to locate any positive detections. △ Less

Submitted 12 August, 2020; originally announced August 2020.

Comments: 26 pages, 8 figures, 3 Tables, accepted for publication in ApJL

Constraining the Kilonova Rate with Zwicky Transient Facility Searches Independent of Gravitational Wave and Short Gamma-ray Burst Triggers

Authors: Igor Andreoni, Erik C. Kool, Ana Sagues Carracedo, Mansi M. Kasliwal, Mattia Bulla, Tomas Ahumada, Michael W. Coughlin, Shreya Anand, Jesper Sollerman, Ariel Goobar, David L. Kaplan, Tegan T. Loveridge, Viraj Karambelkar, Jeff Cooke, Ashot Bagdasaryan, Eric C. Bellm, S. Bradley Cenko, David O. Cook, Kishalay De, Richard Dekany, Alexandre Delacroix, Andrew Drake, Dmitry A. Duev, Christoffer Fremling, V. Zach Golkhou , et al. (12 additional authors not shown)

Abstract: The first binary neutron star merger, GW170817, was accompanied by a radioactivity-powered optical/infrared transient called a kilonova. To date, no compelling kilonova has been found during optical surveys of the sky, independent of gravitational-wave triggers. In this work, we searched the first 23 months of the Zwicky Transient Facility (ZTF) data stream for candidate kilonovae in the form of r… ▽ More The first binary neutron star merger, GW170817, was accompanied by a radioactivity-powered optical/infrared transient called a kilonova. To date, no compelling kilonova has been found during optical surveys of the sky, independent of gravitational-wave triggers. In this work, we searched the first 23 months of the Zwicky Transient Facility (ZTF) data stream for candidate kilonovae in the form of rapidly evolving transients. We combined ZTF alert queries with forced point-spread-function photometry and nightly flux stacking to increase our sensitivity to faint and fast transients. Automatic queries yielded $>11,200$ candidates, 24 of which passed quality checks and strict selection criteria based on a grid of kilonova models tailored for both binary neutron star and neutron star-black hole mergers. None of the candidates in our sample was deemed a possible kilonova after thorough vetting, catalog cross-matching, and study of their color evolution. The sources that passed our selection criteria are dominated by Galactic cataclysmic variables. In addition, we identified two fast transients at high Galactic latitude, one of which is the confirmed afterglow of long-duration GRB190106A, and the other is a possible cosmological afterglow. Using a survey simulation code, we constrained the kilonova rate for a range of models including top-hat and linearly decaying light curves and synthetic light curves obtained with radiative transfer simulations. For prototypical GW170817-like kilonovae, we constrain the rate to be $R < 1775$ Gpc$^{-3}$ yr$^{-1}$ at 95% confidence level by requiring at least 2 high-significance detections. By assuming a population of kilonovae with the same geometry and composition of GW170817 observed under a uniform viewing angle distribution, we obtained a constraint on the rate of $R < 4029$ Gpc$^{-3}$ yr$^{-1}$. △ Less

Submitted 31 July, 2020; originally announced August 2020.

Comments: Submitted for publication in ApJ

Kilonova Luminosity Function Constraints based on Zwicky Transient Facility Searches for 13 Neutron Star Mergers

Authors: Mansi M. Kasliwal, Shreya Anand, Tomas Ahumada, Robert Stein, Ana Sagues Carracedo, Igor Andreoni, Michael W. Coughlin, Leo P. Singer, Erik C. Kool, Kishalay De, Harsh Kumar, Mouza AlMualla, Yuhan Yao, Mattia Bulla, Dougal Dobie, Simeon Reusch, Daniel A. Perley, S. Bradley Cenko, Varun Bhalerao, David L. Kaplan, Jesper Sollerman, Ariel Goobar, Christopher M. Copperwheat, Eric C. Bellm, G. C. Anupama , et al. (78 additional authors not shown)

Abstract: We present a systematic search for optical counterparts to 13 gravitational wave (GW) triggers involving at least one neutron star during LIGO/Virgo's third observing run. We searched binary neutron star (BNS) and neutron star black hole (NSBH) merger localizations with the Zwicky Transient Facility (ZTF) and undertook follow-up with the Global Relay of Observatories Watching Transients Happen (GR… ▽ More We present a systematic search for optical counterparts to 13 gravitational wave (GW) triggers involving at least one neutron star during LIGO/Virgo's third observing run. We searched binary neutron star (BNS) and neutron star black hole (NSBH) merger localizations with the Zwicky Transient Facility (ZTF) and undertook follow-up with the Global Relay of Observatories Watching Transients Happen (GROWTH) collaboration. The GW triggers had a median localization of 4480 deg^2, median distance of 267 Mpc and false alarm rates ranging from 1.5 to 1e-25 per yr. The ZTF coverage had a median enclosed probability of 39%, median depth of 20.8mag, and median response time of 1.5 hr. The O3 follow-up by the GROWTH team comprised 340 UVOIR photometric points, 64 OIR spectra, and 3 radio. We find no promising kilonova (radioactivity-powered counterpart) and we convert the upper limits to constrain the underlying kilonova luminosity function. Assuming that all kilonovae are at least as luminous as GW170817 at discovery (-16.1mag), we calculate our joint probability of detecting zero kilonovae is only 4.2%. If we assume that all kilonovae are brighter than -16.6mag (extrapolated peak magnitude of GW170817) and fade at 1 mag/day (similar to GW170817), the joint probability of zero detections is 7%. If we separate the NSBH and BNS populations, the joint probability of zero detections, assuming all kilonovae are brighter than -16.6mag, is 9.7% for NSBH and 7.9% for BNS mergers. Moreover, <57% (<89%) of putative kilonovae could be brighter than -16.6mag assuming flat (fading) evolution, at 90% confidence. If we further account for the online terrestrial probability for each GW trigger, we find that <68% of putative kilonovae could be brighter than -16.6mag. Comparing to model grids, we find that some kilonovae must have Mej < 0.03 Msun or Xlan>1e-4 or phi>30deg to be consistent with our limits. (Abridged) △ Less

Submitted 19 June, 2020; originally announced June 2020.

Comments: Submitted to ApJ

A tidal disruption event coincident with a high-energy neutrino

Authors: Robert Stein, Sjoert van Velzen, Marek Kowalski, Anna Franckowiak, Suvi Gezari, James C. A. Miller-Jones, Sara Frederick, Itai Sfaradi, Michael F. Bietenholz, Assaf Horesh, Rob Fender, Simone Garrappa, Tomás Ahumada, Igor Andreoni, Justin Belicki, Eric C. Bellm, Markus Böttcher, Valery Brinnel, Rick Burruss, S. Bradley Cenko, Michael W. Coughlin, Virginia Cunningham, Andrew Drake, Glennys R. Farrar, Michael Feeney , et al. (33 additional authors not shown)

Abstract: Cosmic neutrinos provide a unique window into the otherwise-hidden mechanism of particle acceleration in astrophysical objects. A flux of high-energy neutrinos was discovered in 2013, and the IceCube Collaboration recently reported the likely association of one high-energy neutrino with a flare from the relativistic jet of an active galaxy pointed towards the Earth. However a combined analysis of… ▽ More Cosmic neutrinos provide a unique window into the otherwise-hidden mechanism of particle acceleration in astrophysical objects. A flux of high-energy neutrinos was discovered in 2013, and the IceCube Collaboration recently reported the likely association of one high-energy neutrino with a flare from the relativistic jet of an active galaxy pointed towards the Earth. However a combined analysis of many similar active galaxies revealed no excess from the broader population, leaving the vast majority of the cosmic neutrino flux unexplained. Here we present the likely association of a radio-emitting tidal disruption event, AT2019dsg, with a second high-energy neutrino. AT2019dsg was identified as part of our systematic search for optical counterparts to high-energy neutrinos with the Zwicky Transient Facility. The probability of finding any coincident radio-emitting tidal disruption event by chance is 0.5%, while the probability of finding one as bright in bolometric energy flux as AT2019dsg is 0.2%. Our electromagnetic observations can be explained through a multi-zone model, with radio analysis revealing a central engine, embedded in a UV photosphere, that powers an extended synchrotron-emitting outflow. This provides an ideal site for PeV neutrino production. Assuming that the association is genuine, our observations suggest that tidal disruption events with mildly-relativistic outflows contribute to the cosmic neutrino flux. △ Less

Submitted 22 September, 2021; v1 submitted 11 May, 2020; originally announced May 2020.

Comments: Title and text modified during journal review. Version accepted for publication in Nature Astronomy

Journal ref: Nat Astron 5, 510-518 (2021)

DECam-GROWTH Search for the Faint and Distant Binary Neutron Star and Neutron Star-Black Hole Mergers in O3a

Authors: Shreya Anand, Igor Andreoni, Daniel A. Goldstein, Mansi M. Kasliwal, Tomás Ahumada, Jennifer Barnes, Joshua S. Bloom, Mattia Bulla, S. Bradley Cenko, Jeff Cooke, Michael W. Coughlin, Peter E. Nugent, Leo P. Singer

Abstract: Synoptic searches for the optical counterpart to a binary neutron star (BNS) or neutron star-black hole (NSBH) merger can pose significant challenges towards the discovery of kilonovae and performing multi-messenger science. In this work, we describe the advantage of a global multi-telescope network towards this end, with a particular focus on the key and complementary role the Dark Energy Camera… ▽ More Synoptic searches for the optical counterpart to a binary neutron star (BNS) or neutron star-black hole (NSBH) merger can pose significant challenges towards the discovery of kilonovae and performing multi-messenger science. In this work, we describe the advantage of a global multi-telescope network towards this end, with a particular focus on the key and complementary role the Dark Energy Camera (DECam) plays in multi-facility follow-up. We describe the Global Relay of Observatories Watching Transients Happen (GROWTH) Target-of-Opportunity (ToO) Marshal, a common web application we built to ingest events, plan observations, search for transient candidates, and retrieve performance summary statistics for all of the telescopes in our network. Our infrastructure enabled us to conduct observations of two events during O3a, S190426c and S190510g. Furthermore, our analysis of deep DECam observations of S190814bv conducted by the DESGW team, and access to a variety of global follow-up facilities allowed us to place meaningful constraints on the parameters of the kilonova and the merging binary. We emphasize the importance of a global telescope network in conjunction with a power telescope like DECam in performing searches for the counterparts to gravitational-wave sources. △ Less

Submitted 11 March, 2020; originally announced March 2020.

Comments: 9 pages, 4 figures, submitted to Rev. Mex. A. A. (AstroRob 2019 conference proceedings)

GROWTH on S190814bv: Deep Synoptic Limits on the Optical/Near-Infrared Counterpart to a Neutron Star-Black Hole Merger

Authors: Igor Andreoni, Daniel A. Goldstein, Mansi M. Kasliwal, Peter E. Nugent, Rongpu Zhou, Jeffrey A. Newman, Mattia Bulla, Francois Foucart, Kenta Hotokezaka, Ehud Nakar, Samaya Nissanke, Geert Raaijmakers, Joshua S. Bloom, Kishalay De, Jacob E. Jencson, Charlotte Ward, Tomás Ahumada, Shreya Anand, David A. H. Buckley, Maria D. Caballero-García, Alberto J. Castro-Tirado, Christopher M. Copperwheat, Michael W. Coughlin, S. Bradley Cenko, Mariusz Gromadzki , et al. (27 additional authors not shown)

Abstract: On 2019 August 14, the Advanced LIGO and Virgo interferometers detected the high-significance gravitational wave (GW) signal S190814bv. The GW data indicated that the event resulted from a neutron star--black hole (NSBH) merger, or potentially a low-mass binary black hole merger. Due to the low false alarm rate and the precise localization (23 deg$^2$ at 90\%), S190814bv presented the community wi… ▽ More On 2019 August 14, the Advanced LIGO and Virgo interferometers detected the high-significance gravitational wave (GW) signal S190814bv. The GW data indicated that the event resulted from a neutron star--black hole (NSBH) merger, or potentially a low-mass binary black hole merger. Due to the low false alarm rate and the precise localization (23 deg$^2$ at 90\%), S190814bv presented the community with the best opportunity yet to directly observe an optical/near-infrared counterpart to a NSBH merger. To search for potential counterparts, the GROWTH collaboration performed real-time image subtraction on 6 nights of public Dark Energy Camera (DECam) images acquired in the three weeks following the merger, covering $>$98\% of the localization probability. Using a worldwide network of follow-up facilities, we systematically undertook spectroscopy and imaging of optical counterpart candidates. Combining these data with a photometric redshift catalog, we ruled out each candidate as the counterpart to S190814bv and we placed deep, uniform limits on the optical emission associated with S190814bv. For the nearest consistent GW distance, radiative transfer simulations of NSBH mergers constrain the ejecta mass of S190814bv to be $M_\mathrm{ej} < 0.04$~$M_{\odot}$ at polar viewing angles, or $M_\mathrm{ej} < 0.03$~$M_{\odot}$ if the opacity is $κ< 2$~cm$^2$g$^{-1}$. Assuming a tidal deformability for the neutron star at the high end of the range compatible with GW170817 results, our limits would constrain the BH spin component aligned with the orbital momentum to be $ χ< 0.7$ for mass ratios $Q < 6$, with weaker constraints for more compact neutron stars. We publicly release the photometry from this campaign at http://www.astro.caltech.edu/~danny/static/s190814bv. △ Less

Submitted 31 December, 2019; v1 submitted 29 October, 2019; originally announced October 2019.

Comments: Accepted to ApJ with no major changes

GROWTH on S190425z: Searching thousands of square degrees to identify an optical or infrared counterpart to a binary neutron star merger with the Zwicky Transient Facility and Palomar Gattini IR

Authors: Michael W. Coughlin, Tomás Ahumada, Shreya Anand, Kishalay De, Matthew J. Hankins, Mansi M. Kasliwal, Leo P. Singer, Eric C. Bellm, Igor Andreoni, S. Bradley Cenko, Jeff Cooke, Christopher M. Copperwheat, Alison M. Dugas, Jacob E. Jencson, Daniel A. Perley, Po-Chieh Yu, Varun Bhalerao, Harsh Kumar, Joshua S. Bloom, G. C. Anupama, Michael C. B. Ashley, Ashot Bagdasaryan, Rahul Biswas, David A. H. Buckley, Kevin B. Burdge , et al. (54 additional authors not shown)

Abstract: The third observing run by LVC has brought the discovery of many compact binary coalescences. Following the detection of the first binary neutron star merger in this run (LIGO/Virgo S190425z), we performed a dedicated follow-up campaign with the Zwicky Transient Facility (ZTF) and Palomar Gattini-IR telescopes. The initial skymap of this single-detector gravitational wave (GW) trigger spanned most… ▽ More The third observing run by LVC has brought the discovery of many compact binary coalescences. Following the detection of the first binary neutron star merger in this run (LIGO/Virgo S190425z), we performed a dedicated follow-up campaign with the Zwicky Transient Facility (ZTF) and Palomar Gattini-IR telescopes. The initial skymap of this single-detector gravitational wave (GW) trigger spanned most of the sky observable from Palomar Observatory. Covering 8000 deg$^2$ of the initial skymap over the next two nights, corresponding to 46\% integrated probability, ZTF system achieved a depth of $\approx$\,21 $m_\textrm{AB}$ in $g$- and $r$-bands. Palomar Gattini-IR covered 2200 square degrees in $J$-band to a depth of 15.5\,mag, including 32\% integrated probability based on the initial sky map. The revised skymap issued the following day reduced these numbers to 21\% for the Zwicky Transient Facility and 19\% for Palomar Gattini-IR. We narrowed 338,646 ZTF transient "alerts" over the first two nights of observations to 15 candidate counterparts. Two candidates, ZTF19aarykkb and ZTF19aarzaod, were particularly compelling given that their location, distance, and age were consistent with the GW event, and their early optical lightcurves were photometrically consistent with that of kilonovae. These two candidates were spectroscopically classified as young core-collapse supernovae. The remaining candidates were ruled-out as supernovae. Palomar Gattini-IR did not identify any viable candidates with multiple detections only after merger time. We demonstrate that even with single-detector GW events localized to thousands of square degrees, systematic kilonova discovery is feasible. △ Less

Submitted 4 October, 2019; v1 submitted 29 July, 2019; originally announced July 2019.

GROWTH on S190510g: DECam Observation Planning and Follow-Up of a Distant Binary Neutron Star Merger Candidate

Authors: Igor Andreoni, Daniel A. Goldstein, Shreya Anand, Michael W. Coughlin, Leo P. Singer, Tomás Ahumada, Michael Medford, Erik C. Kool, Sara Webb, Mattia Bulla, Joshua S. Bloom, Mansi M. Kasliwal, Peter E. Nugent, Ashot Bagdasaryan, Jennifer Barnes, David O. Cook, Jeff Cooke, Dmitry A. Duev, U. Christoffer Fremling, Pradip Gatkine, V. Zach Golkhou, Albert K. H. Kong, Ashish Mahabal, Jorge Martínez-Palomera, Duo Tao , et al. (1 additional authors not shown)

Abstract: The first two months of the third Advanced LIGO and Virgo observing run (2019 April-May) showed that distant gravitational wave (GW) events can now be readily detected. Three candidate mergers containing neutron stars (NS) were reported in a span of 15 days, all likely located more than 100 Mpc away. However, distant events such as the three new NS mergers are likely to be coarsely localized, whic… ▽ More The first two months of the third Advanced LIGO and Virgo observing run (2019 April-May) showed that distant gravitational wave (GW) events can now be readily detected. Three candidate mergers containing neutron stars (NS) were reported in a span of 15 days, all likely located more than 100 Mpc away. However, distant events such as the three new NS mergers are likely to be coarsely localized, which highlights the importance of facilities and scheduling systems that enable deep observations over hundreds to thousands of square degrees to detect the electromagnetic counterparts. On 2019-05-10 02:59:39.292 UT the GW candidate S190510g was discovered and initially classified as a BNS merger with 98% probability. The GW event was localized within an area of 3462 deg2, later refined to 1166 deg2 (90%) at a distance of 227 +- 92 Mpc. We triggered Target of Opportunity observations with the Dark Energy Camera (DECam), a wide-field optical imager mounted at the prime focus of the 4m Blanco Telescope at CTIO in Chile. This Letter describes our DECam observations and our real-time analysis results, focusing in particular on the design and implementation of the observing strategy. Within 24 hours of the merger time, we observed 65% of the total enclosed probability of the final skymap with an observing efficiency of 94%. We identified and publicly announced 13 candidate counterparts. S190510g was re-classified 1.7 days after the merger, after our observations were completed, with a "binary neutron star merger" probability reduced from 98% to 42% in favor of a "terrestrial" classification. △ Less

Submitted 22 July, 2019; v1 submitted 31 May, 2019; originally announced June 2019.

Comments: 4 figures, 3 tables, accepted for publication in ApJL

GROWTH on S190426c. II. Real-Time Search for a Counterpart to the Probable Neutron Star-Black Hole Merger using an Automated Difference Imaging Pipeline for DECam

Authors: Daniel A. Goldstein, Igor Andreoni, Peter E. Nugent, Mansi M. Kasliwal, Michael W. Coughlin, Shreya Anand, Joshua S. Bloom, Jorge Martínez-Palomera, Keming Zhang, Tomás Ahumada, Ashot Bagdasaryan, Jeff Cooke, Kishalay De, Dmitry A. Duev, U. Christoffer Fremling, Pradip Gatkine, Matthew Graham, Eran O. Ofek, Leo P. Singer, Lin Yan

Abstract: The discovery of a transient kilonova following the gravitational-wave event GW170817 highlighted the critical need for coordinated rapid and wide-field observations, inference, and follow-up across the electromagnetic spectrum. In the Southern hemisphere, the Dark Energy Camera (DECam) on the Blanco 4-m telescope is well-suited to this task, as it is able to cover wide-fields quickly while still… ▽ More The discovery of a transient kilonova following the gravitational-wave event GW170817 highlighted the critical need for coordinated rapid and wide-field observations, inference, and follow-up across the electromagnetic spectrum. In the Southern hemisphere, the Dark Energy Camera (DECam) on the Blanco 4-m telescope is well-suited to this task, as it is able to cover wide-fields quickly while still achieving the depths required to find kilonovae like the one accompanying GW170817 to $\sim$500 Mpc, the binary neutron star horizon distance for current generation of LIGO/Virgo collaboration (LVC) interferometers. Here, as part of the multi-facility followup by the Global Relay of Observatories Watching Transients Happen (GROWTH) collaboration, we describe the observations and automated data movement, data reduction, candidate discovery, and vetting pipeline of our target-of-opportunity DECam observations of S190426c, the first possible neutron star--black hole merger detected via gravitational waves. Starting 7.5hr after S190426c, over 11.28\,hr of observations, we imaged an area of 525\,deg$^2$ ($r$-band) and 437\,deg$^2$ ($z$-band); this was 16.3\% of the total original localization probability and nearly all of the probability density visible from the Southern hemisphere. The machine-learning based pipeline was optimized for fast turnaround, delivering transient candidates for human vetting within 17 minutes, on average, of shutter closure. We reported nine promising counterpart candidates 2.5 hours before the end of our observations. Our observations yielded no detection of a bona fide counterpart to $m_z = 22.5$ and $m_r = 22.9$ at the 5$σ$ level of significance, consistent with the refined LVC positioning. We view these observations and rapid inferencing as an important real-world test for this novel end-to-end wide-field pipeline. △ Less

Submitted 16 May, 2019; originally announced May 2019.

Comments: submitted to AAS journals

2900 square degree search for the optical counterpart of short gamma-ray burst GRB 180523B with the Zwicky Transient Facility

Authors: Michael W. Coughlin, Tomas Ahumada, S. Bradley Cenko, Virginia Cunningham, Shaon Ghosh, Leo P. Singer, Eric C. Bellm, Eric Burns, Kishalay De, Adam Goldstein, V. Zach Golkhou, David L. Kaplan, Mansi M. Kasliwal, Daniel A. Perley, Jesper Sollerman, Ashot Bagdasaryan, Richard G. Dekany, Dmitry A. Duev, Michael Feeney, Matthew J. Graham, David Hale, Shri R. Kulkarni, Thomas Kupfer, Russ R. Laher, Ashish Mahabal , et al. (9 additional authors not shown)

Abstract: There is significant interest in the models for production of short gamma-ray bursts. Until now, the number of known short gamma-ray bursts with multi-wavelength afterglows has been small. While the {\it Fermi} Gamma-Ray Burst Monitor detects many gamma-ray bursts relative to the Neil Gehrels {\it Swift} Observatory, the large localization regions makes the search for counterparts difficult. With… ▽ More There is significant interest in the models for production of short gamma-ray bursts. Until now, the number of known short gamma-ray bursts with multi-wavelength afterglows has been small. While the {\it Fermi} Gamma-Ray Burst Monitor detects many gamma-ray bursts relative to the Neil Gehrels {\it Swift} Observatory, the large localization regions makes the search for counterparts difficult. With the Zwicky Transient Facility recently achieving first light, it is now fruitful to use its combination of depth ($m_\textrm{AB} \sim 20.6$), field of view ($\approx$ 47 square degrees), and survey cadence (every $\sim 3$ days) to perform Target of Opportunity observations. We demonstrate this capability on GRB 180523B, which was recently announced by the {\it Fermi} Gamma-Ray Burst Monitor as a short gamma-ray burst. ZTF imaged $\approx$ 2900\,square degrees of the localization region, resulting in the coverage of 61.6\,\% of the enclosed probability over 2 nights to a depth of $m_\textrm{AB} \sim 20.5$. We characterized 14 previously unidentified transients, and none were found to be consistent with a short gamma-ray burst counterpart. This search with the Zwicky Transient Facility shows it is an efficient camera for searching for coarsely-localized short gamma-ray burst and gravitational-wave counterparts, allowing for a sensitive search with minimal interruption to its nominal cadence. △ Less

Submitted 31 January, 2019; originally announced January 2019.

The Kitt Peak Electron Multiplying CCD demonstrator

Authors: Michael W. Coughlin, Richard G. Dekany, Dmitry A. Duev, Michael Feeney, S. R. Kulkarni, Reed Riddle, Tomás Ahumada, Kevin Burdge, Alison M. Dugas, Christoffer U. Fremling, Gregg Hallinan, Thomas A. Prince, Jan van Roestel

Abstract: The Kitt Peak Electron Multiplying CCD (EMCCD) demonstrator is a new instrument that has been developed for use at the Kitt Peak National Observatory's 84-inch telescope. The EMCCD enables single-band optical imaging in the Sloan \textit{g} and \textit{r} bands and Johnson \textit{UVRI} filters. The EMCCD is chosen for its sub-electron effective read noise using large multiplicative gains. As a re… ▽ More The Kitt Peak Electron Multiplying CCD (EMCCD) demonstrator is a new instrument that has been developed for use at the Kitt Peak National Observatory's 84-inch telescope. The EMCCD enables single-band optical imaging in the Sloan \textit{g} and \textit{r} bands and Johnson \textit{UVRI} filters. The EMCCD is chosen for its sub-electron effective read noise using large multiplicative gains. As a result, frame rates of greater than 1\,Hz are possible. The field-of-view is 4.4$^\prime$ $\times$ 4.4$^\prime$ and the pixel size is 0.259$^{\prime\prime}$. This camera, coupled with a fully roboticized telescope, is ideal for follow-up of short period, white dwarf binary candidates, as well as short duration transient and periodic sources identified by large field-of-view all-sky surveys such as the Zwicky Transient Facility. △ Less

Submitted 14 January, 2019; originally announced January 2019.

The Fast, Luminous Ultraviolet Transient AT2018cow: Extreme Supernova, or Disruption of a Star by an Intermediate-Mass Black Hole?

Authors: Daniel A. Perley, Paolo A. Mazzali, Lin Yan, S. Bradley Cenko, Suvi Gezari, Kirsty Taggart, Nadia Blagorodnova, Christoffer Fremling, Brenna Mockler, Avinash Singh, Nozomu Tominaga, Masaomi Tanaka, Alan M. Watson, Tomás Ahumada, G. C. Anupama, Chris Ashall, Rosa L. Becerra, David Bersier, Varun Bhalerao, Joshua S. Bloom, Nathaniel R. Butler, Chris Copperwheat, Michael W. Coughlin, Kishalay De, Andrew J. Drake , et al. (38 additional authors not shown)

Abstract: Wide-field optical surveys have begun to uncover large samples of fast (t_rise < 5d), luminous (M_peak < -18), blue transients. While commonly attributed to the breakout of a supernova shock into a dense wind, the great distances to the transients of this class found so far have hampered detailed investigation of their properties. We present photometry and spectroscopy from a comprehensive worldwi… ▽ More Wide-field optical surveys have begun to uncover large samples of fast (t_rise < 5d), luminous (M_peak < -18), blue transients. While commonly attributed to the breakout of a supernova shock into a dense wind, the great distances to the transients of this class found so far have hampered detailed investigation of their properties. We present photometry and spectroscopy from a comprehensive worldwide campaign to observe AT2018cow (ATLAS18qqn), the first fast-luminous optical transient to be found in real time at low redshift. Our first spectra (<2 days after discovery) are entirely featureless. A very broad absorption feature suggestive of near-relativistic velocities develops between 3-8 days, then disappears. Broad emission features of H and He develop after >10 days. The spectrum remains extremely hot throughout its evolution, and the photospheric radius contracts with time (receding below R<10^14 cm after 1 month). This behaviour does not match that of any known supernova, although a relativistic jet within a fallback supernova could explain some of the observed features. Alternatively, the transient could originate from the disruption of a star by an intermediate-mass black hole, although this would require long-lasting emission of highly super-Eddington thermal radiation. In either case, AT2018cow suggests that the population of fast luminous transients represents a new class of astrophysical event. Intensive follow-up of this event in its late phases, and of any future events found at comparable distance, will be essential to better constrain their origins. △ Less

Submitted 23 November, 2018; v1 submitted 2 August, 2018; originally announced August 2018.

Comments: Corrected Figure 8 / Table 4 to use final fits. Includes machine-readable photometry table (hopefully for real this time)