Showing 1–50 of 301 results for author: Coughlin, M W

Candidate strongly-lensed Type Ia supernovae in the Zwicky Transient Facility archive

Authors: A. Townsend, J. Nordin, A. Sagués Carracedo, M. Kowalski, N. Arendse, S. Dhawan, A. Goobar, J. Johansson, E. Mörtsell, S. Schulze, I. Andreoni, E. Fernández, A. G. Kim, P. E. Nugent, F. Prada, M. Rigault, N. Sarin, D. Sharma, E. C. Bellm, M. W. Coughlin, R. Dekany, S. L. Groom, L. Lacroix, R. R. Laher, R. Riddle , et al. (39 additional authors not shown)

Abstract: Gravitationally lensed Type Ia supernovae (glSNe Ia) are unique astronomical tools for studying cosmological parameters, distributions of dark matter, the astrophysics of the supernovae and the intervening lensing galaxies themselves. Only a few highly magnified glSNe Ia have been discovered by ground-based telescopes, such as the Zwicky Transient Facility (ZTF), but simulations predict the existe… ▽ More Gravitationally lensed Type Ia supernovae (glSNe Ia) are unique astronomical tools for studying cosmological parameters, distributions of dark matter, the astrophysics of the supernovae and the intervening lensing galaxies themselves. Only a few highly magnified glSNe Ia have been discovered by ground-based telescopes, such as the Zwicky Transient Facility (ZTF), but simulations predict the existence of a fainter, undetected population. We present a systematic search in the ZTF archive of alerts from 1 June 2019 to 1 September 2022. Using the AMPEL platform, we developed a pipeline that distinguishes candidate glSNe Ia from other variable sources. Initial cuts were applied to the ZTF alert photometry before forced photometry was obtained for the remaining candidates. Additional cuts were applied to refine the candidates based on their light curve colours, lens galaxy colours, and the resulting parameters from fits to the SALT2 SN Ia template. Candidates were also cross-matched with the DESI spectroscopic catalogue. Seven transients passed all the cuts and had an associated galaxy DESI redshift, which we present as glSN Ia candidates. While superluminous supernovae (SLSNe) cannot be fully rejected, two events, ZTF19abpjicm and ZTF22aahmovu, are significantly different from typical SLSNe and their light curves can be modelled as two-image glSN Ia systems. From this two-image modelling, we estimate time delays of 22 $\pm$ 3 and 34 $\pm$ 1 days for the two events, respectively, which suggests that we have uncovered a population with longer time delays. The pipeline is efficient and sensitive enough to parse full alert streams. It is currently being applied to the live ZTF alert stream to identify and follow-up future candidates while active. This pipeline could be the foundation for glSNe Ia searches in future surveys, like the Vera C. Rubin Observatory's Legacy Survey of Space and Time. △ Less

Submitted 28 May, 2024; originally announced May 2024.

Comments: 21 pages, 15 figures

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

An Online Framework for Fitting Fast Transient Lightcurves

Authors: Tyler Barna, Brandon Reed, Igor Andreoni, Michael W. Coughlin, Tim Dietrich, Steven L. Groom, Theophile Jegou du Laz, Peter T. H Pang, Josiah N. Purdum, Ben Rusholme

Abstract: The identification of extragalactic fast optical transients (eFOTs) as potential multi-messenger sources is one of the main challenges in time-domain astronomy. However, recent developments have allowed for probes of rapidly-evolving transients. With the increasing number of alert streams from optical time-domain surveys, the next paradigm is building technologies to rapidly identify the most inte… ▽ More The identification of extragalactic fast optical transients (eFOTs) as potential multi-messenger sources is one of the main challenges in time-domain astronomy. However, recent developments have allowed for probes of rapidly-evolving transients. With the increasing number of alert streams from optical time-domain surveys, the next paradigm is building technologies to rapidly identify the most interesting transients for follow-up. One effort to make this possible is the fitting of objects to a variety of eFOT lightcurve models such as kilonovae and $γ$-ray burst (GRB) afterglows. In this work, we describe a new framework designed to efficiently fit transients to light curve models and flag them for further follow-up. We describe the pipeline's workflow and a handful of performance metrics, including the nominal sampling time for each model. We highlight as examples ZTF20abwysqy, the shortest long gamma ray burst discovered to date, and ZTF21abotose, a core-collapse supernova initially identified as a potential kilonova candidate. △ Less

Submitted 26 April, 2024; originally announced April 2024.

Comments: Accepted to MNRAS, 11 pages, 9 figures

The Lunar Gravitational-wave Antenna: Mission Studies and Science Case

Authors: Parameswaran Ajith, Pau Amaro Seoane, Manuel Arca Sedda, Riccardo Arcodia, Francesca Badaracco, Enis Belgacem, Stefano Benetti, Alexey Bobrick, Alessandro Bonforte, Elisa Bortolas, Valentina Braito, Marica Branchesi, Adam Burrows, Enrico Cappellaro, Roberto Della Ceca, Chandrachur Chakraborty, Shreevathsa Chalathadka Subrahmanya, Michael W. Coughlin, Stefano Covino, Andrea Derdzinski, Aayushi Doshi, Maurizio Falanga, Stefano Foffa, Alessia Franchini, Alessandro Frigeri , et al. (58 additional authors not shown)

Abstract: The Lunar Gravitational-wave Antenna (LGWA) is a proposed array of next-generation inertial sensors to monitor the response of the Moon to gravitational waves (GWs). Given the size of the Moon and the expected noise produced by the lunar seismic background, the LGWA would be able to observe GWs from about 1 mHz to 1 Hz. This would make the LGWA the missing link between space-borne detectors like L… ▽ More The Lunar Gravitational-wave Antenna (LGWA) is a proposed array of next-generation inertial sensors to monitor the response of the Moon to gravitational waves (GWs). Given the size of the Moon and the expected noise produced by the lunar seismic background, the LGWA would be able to observe GWs from about 1 mHz to 1 Hz. This would make the LGWA the missing link between space-borne detectors like LISA with peak sensitivities around a few millihertz and proposed future terrestrial detectors like Einstein Telescope or Cosmic Explorer. In this article, we provide a first comprehensive analysis of the LGWA science case including its multi-messenger aspects and lunar science with LGWA data. We also describe the scientific analyses of the Moon required to plan the LGWA mission. △ Less

Submitted 14 April, 2024; originally announced April 2024.

A machine-learning pipeline for real-time detection of gravitational waves from compact binary coalescences

Authors: Ethan Marx, William Benoit, Alec Gunny, Rafia Omer, Deep Chatterjee, Ricco C. Venterea, Lauren Wills, Muhammed Saleem, Eric Moreno, Ryan Raikman, Ekaterina Govorkova, Dylan Rankin, Michael W. Coughlin, Philip Harris, Erik Katsavounidis

Abstract: The promise of multi-messenger astronomy relies on the rapid detection of gravitational waves at very low latencies ($\mathcal{O}$(1\,s)) in order to maximize the amount of time available for follow-up observations. In recent years, neural-networks have demonstrated robust non-linear modeling capabilities and millisecond-scale inference at a comparatively small computational footprint, making them… ▽ More The promise of multi-messenger astronomy relies on the rapid detection of gravitational waves at very low latencies ($\mathcal{O}$(1\,s)) in order to maximize the amount of time available for follow-up observations. In recent years, neural-networks have demonstrated robust non-linear modeling capabilities and millisecond-scale inference at a comparatively small computational footprint, making them an attractive family of algorithms in this context. However, integration of these algorithms into the gravitational-wave astrophysics research ecosystem has proven non-trivial. Here, we present the first fully machine learning-based pipeline for the detection of gravitational waves from compact binary coalescences (CBCs) running in low-latency. We demonstrate this pipeline to have a fraction of the latency of traditional matched filtering search pipelines while achieving state-of-the-art sensitivity to higher-mass stellar binary black holes. △ Less

Submitted 27 March, 2024; originally announced March 2024.

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.

SN 2023zaw: an ultra-stripped, nickel-poor supernova from a low-mass progenitor

Authors: Kaustav K. Das, Christoffer Fremling, Mansi M. Kasliwal, Steve Schulze, Jesper Sollerman, Viraj Karambelkar, Sam Rose, Shreya Anand, Igor Andreoni, Marie Aubert, Sean J. Brennan, S. Bradley Cenko, Michael W. Coughlin, B. O'Connor, Kishalay De, Jim Fuller, Matthew Graham, Erica Hammerstein, Annastasia Haynie, K-Ryan Hinds, Io Kleiser, S. R. Kulkarni, Zeren Lin, Chang Liu, Ashish A. Mahabal , et al. (12 additional authors not shown)

Abstract: We present SN 2023zaw $-$ a sub-luminous ($\mathrm{M_r} = -16.7$ mag) and rapidly-evolving supernova ($\mathrm{t_{1/2,r}} = 4.9$ days), with the lowest nickel mass ($\approx0.002$ $\mathrm{M_\odot}$) measured among all stripped-envelope supernovae discovered to date. The photospheric spectra are dominated by broad He I and Ca NIR emission lines with velocities of $\sim10\ 000 - 12\ 000$… ▽ More We present SN 2023zaw $-$ a sub-luminous ($\mathrm{M_r} = -16.7$ mag) and rapidly-evolving supernova ($\mathrm{t_{1/2,r}} = 4.9$ days), with the lowest nickel mass ($\approx0.002$ $\mathrm{M_\odot}$) measured among all stripped-envelope supernovae discovered to date. The photospheric spectra are dominated by broad He I and Ca NIR emission lines with velocities of $\sim10\ 000 - 12\ 000$ $\mathrm{km\ s^{-1}}$. The late-time spectra show prominent narrow He I emission lines at $\sim$1000$\ \mathrm{km\ s^{-1}}$, indicative of interaction with He-rich circumstellar material. SN 2023zaw is located in the spiral arm of a star-forming galaxy. We perform radiation-hydrodynamical and analytical modeling of the lightcurve by fitting with a combination of shock-cooling emission and nickel decay. The progenitor has a best-fit envelope mass of $\approx0.2$ $\mathrm{M_\odot}$ and an envelope radius of $\approx50$ $\mathrm{R_\odot}$. The extremely low nickel mass and low ejecta mass ($\approx0.5$ $\mathrm{M_\odot}$) suggest an ultra-stripped SN, which originates from a mass-losing low mass He-star (ZAMS mass $<$ 10 $\mathrm{M_\odot}$) in a close binary system. This is a channel to form double neutron star systems, whose merger is detectable with LIGO. SN 2023zaw underscores the existence of a previously undiscovered population of extremely low nickel mass ($< 0.005$ $\mathrm{M_\odot}$) stripped-envelope supernovae, which can be explored with deep and high-cadence transient surveys. △ Less

Submitted 14 March, 2024; v1 submitted 12 March, 2024; originally announced March 2024.

Comments: Submitted to ApJL. Host galaxy analysis added. Comments are welcome

Ultralight vector dark matter search using data from the KAGRA O3GK run

Authors: The LIGO Scientific Collaboration, the Virgo Collaboration, the KAGRA Collaboration, A. G. Abac, R. Abbott, H. Abe, I. Abouelfettouh, F. Acernese, K. Ackley, C. Adamcewicz, S. Adhicary, N. Adhikari, R. X. Adhikari, V. K. Adkins, V. B. Adya, C. Affeldt, D. Agarwal, M. Agathos, O. D. Aguiar, I. Aguilar, L. Aiello, A. Ain, P. Ajith, T. Akutsu, S. Albanesi , et al. (1778 additional authors not shown)

Abstract: Among the various candidates for dark matter (DM), ultralight vector DM can be probed by laser interferometric gravitational wave detectors through the measurement of oscillating length changes in the arm cavities. In this context, KAGRA has a unique feature due to differing compositions of its mirrors, enhancing the signal of vector DM in the length change in the auxiliary channels. Here we prese… ▽ More Among the various candidates for dark matter (DM), ultralight vector DM can be probed by laser interferometric gravitational wave detectors through the measurement of oscillating length changes in the arm cavities. In this context, KAGRA has a unique feature due to differing compositions of its mirrors, enhancing the signal of vector DM in the length change in the auxiliary channels. Here we present the result of a search for $U(1)_{B-L}$ gauge boson DM using the KAGRA data from auxiliary length channels during the first joint observation run together with GEO600. By applying our search pipeline, which takes into account the stochastic nature of ultralight DM, upper bounds on the coupling strength between the $U(1)_{B-L}$ gauge boson and ordinary matter are obtained for a range of DM masses. While our constraints are less stringent than those derived from previous experiments, this study demonstrates the applicability of our method to the lower-mass vector DM search, which is made difficult in this measurement by the short observation time compared to the auto-correlation time scale of DM. △ Less

Submitted 5 March, 2024; originally announced March 2024.

Comments: 20 pages, 5 figures

Report number: LIGO-P2300250

Correlated 0.01Hz-40Hz seismic and Newtonian noise and its impact on future gravitational-wave detectors

Authors: Kamiel Janssens, Guillaume Boileau, Nelson Christensen, Nick van Remortel, Francesca Badaracco, Benjamin Canuel, Alessandro Cardini, Andrea Contu, Michael W. Coughlin, Jean-Baptiste Decitre, Rosario De Rosa, Matteo Di Giovanni, Domenico D'Urso, Stéphane Gaffet, Carlo Giunchi, Jan Harms, Soumen Koley, Valentina Mangano, Luca Naticchioni, Marco Olivieri, Federico Paoletti, Davide Rozza, Dylan O. Sabulsky, Shahar Shani-Kadmiel, Lucia Trozzo

Abstract: We report correlations in underground seismic measurements with horizontal separations of several hundreds of meters to a few kilometers in the frequency range 0.01Hz to 40Hz. These seismic correlations could threaten science goals of planned interferometric gravitational-wave detectors such as the Einstein Telescope as well as atom interferometers such as MIGA and ELGAR. We use seismic measuremen… ▽ More We report correlations in underground seismic measurements with horizontal separations of several hundreds of meters to a few kilometers in the frequency range 0.01Hz to 40Hz. These seismic correlations could threaten science goals of planned interferometric gravitational-wave detectors such as the Einstein Telescope as well as atom interferometers such as MIGA and ELGAR. We use seismic measurements from four different sites, i.e. the former Homestake mine (USA) as well as two candidate sites for the Einstein Telescope, Sos Enattos (IT) and Euregio Maas-Rhein (NL-BE-DE) and the site housing the MIGA detector, LSBB (FR). At all sites, we observe significant coherence for at least 50% of the time in the majority of the frequency region of interest. Based on the observed correlations in the seismic fields, we predict levels of correlated Newtonian noise from body waves. We project the effect of correlated Newtonian noise from body waves on the capabilities of the triangular design of the Einstein Telescope's to observe an isotropic gravitational-wave background (GWB) and find that, even in case of the most quiet site, its sensitivity will be affected up to $\sim$20Hz. The resolvable amplitude of a GWB signal with a negatively sloped power-law behaviour would be reduced by several orders of magnitude. However, the resolvability of a power-law signal with a slope of e.g. $α=0$ ($α=2/3$) would be more moderately affected by a factor $\sim$ 6-9 ($\sim$3-4) in case of a low noise environment. Furthermore, we bolster confidence in our results by showing that transient noise features have a limited impact on the presented results. △ Less

Submitted 27 February, 2024; originally announced February 2024.

An overview of existing and new nuclear and astrophysical constraints on the equation of state of neutron-rich dense matter

Authors: Hauke Koehn, Henrik Rose, Peter T. H. Pang, Rahul Somasundaram, Brendan T. Reed, Ingo Tews, Adrian Abac, Oleg Komoltsev, Nina Kunert, Aleksi Kurkela, Michael W. Coughlin, Brian F. Healy, Tim Dietrich

Abstract: Through continuous progress in nuclear theory and experiment and an increasing number of neutron-star observations, a multitude of information about the equation of state (EOS) for matter at extreme densities is available. Here, we apply these different pieces of data individually to a broad set of physics-agnostic candidate EOSs and analyze the resulting constraints. Specifically, we make use of… ▽ More Through continuous progress in nuclear theory and experiment and an increasing number of neutron-star observations, a multitude of information about the equation of state (EOS) for matter at extreme densities is available. Here, we apply these different pieces of data individually to a broad set of physics-agnostic candidate EOSs and analyze the resulting constraints. Specifically, we make use of information from chiral effective field theory, perturbative quantum chromodynamics, as well as data from heavy-ion collisions and the PREX-II and CREX experiments. We also investigate the impact of current mass and radius measurements of neutron stars, such as radio timing measurements of heavy pulsars, NICER data, and other X-ray observations. We augment these by reanalyses of the gravitational-wave (GW) signal GW170817, its associated kilonova AT2017gfo and gamma-ray burst afterglow, the GW signal GW190425, and the GRB211211A afterglow, where we use improved models for the tidal waveform and kilonova light curves. Additionally, we consider the postmerger fate of GW170817 and its consequences for the EOS. This large and diverse set of constraints is eventually combined in numerous ways to explore limits on quantities such as the typical neutron-star radius, the maximum neutron-star mass, the nuclear symmetry-energy parameters, and the speed of sound. Based on the priors from our EOS candidate set, we find the radius of the canonical 1.4 M$_\odot$ neutron star to be $R_{1.4}= 12.27_{-0.94}^{+0.83}$ km and the TOV mass $M_{\rm TOV}= 2.26_{-0.22}^{+0.45}$ M$_\odot$ at 95% credibility, when including those constraints where systematic uncertainties are deemed small. A less conservative approach, combining all the presented constraints, similarly yields $R_{1.4}= 12.20_{-0.50}^{+0.53}$ km and $M_{\rm TOV}= 2.31_{-0.20}^{+0.08}$ M$_\odot$. △ Less

Submitted 6 February, 2024; originally announced February 2024.

Comments: 49 pages, 32 figures, webinterface for custom constraint combinations in https://enlil.gw.physik.uni-potsdam.de/eos_constraints

Report number: LA-UR-24-20420

The Zwicky Transient Facility Bright Transient Survey. III. $\texttt{BTSbot}$: Automated Identification and Follow-up of Bright Transients with Deep Learning

Authors: Nabeel Rehemtulla, Adam A. Miller, Theophile Jegou Du Laz, Michael W. Coughlin, Christoffer Fremling, Daniel A. Perley, Yu-Jing Qin, Jesper Sollerman, Ashish A. Mahabal, Russ R. Laher, Reed Riddle, Ben Rusholme, Shrinivas R. Kulkarni

Abstract: The Bright Transient Survey (BTS) aims to obtain a classification spectrum for all bright ($m_\mathrm{peak}\,\leq\,18.5\,$mag) extragalactic transients found in the Zwicky Transient Facility (ZTF) public survey. BTS critically relies on visual inspection ("scanning") to select targets for spectroscopic follow-up, which, while effective, has required a significant time investment over the past… ▽ More The Bright Transient Survey (BTS) aims to obtain a classification spectrum for all bright ($m_\mathrm{peak}\,\leq\,18.5\,$mag) extragalactic transients found in the Zwicky Transient Facility (ZTF) public survey. BTS critically relies on visual inspection ("scanning") to select targets for spectroscopic follow-up, which, while effective, has required a significant time investment over the past $\sim5$ yr of ZTF operations. We present $\texttt{BTSbot}$, a multi-modal convolutional neural network, which provides a bright transient score to individual ZTF detections using their image data and 25 extracted features. $\texttt{BTSbot}$ is able to eliminate the need for daily human scanning by automatically identifying and requesting spectroscopic follow-up observations of new bright transient candidates. $\texttt{BTSbot}$ recovers all bright transients in our test split and performs on par with scanners in terms of identification speed (on average, $\sim$1 hour quicker than scanners). We also find that $\texttt{BTSbot}$ is not significantly impacted by any data shift by comparing performance across a concealed test split and a sample of very recent BTS candidates. $\texttt{BTSbot}$ has been integrated into Fritz and $\texttt{Kowalski}$, ZTF's first-party marshal and alert broker, and now sends automatic spectroscopic follow-up requests for the new transients it identifies. During the month of October 2023, $\texttt{BTSbot}$ selected 296 sources in real-time, 93% of which were real extragalactic transients. With $\texttt{BTSbot}$ and other automation tools, the BTS workflow has produced the first fully automatic end-to-end discovery and classification of a transient, representing a significant reduction in the human-time needed to scan. Future development has tremendous potential for creating similar models to identify and request follow-up observations for specific types of transients. △ Less

Submitted 26 January, 2024; originally announced January 2024.

Comments: 26 pages, 12 figures; to be submitted to ApJ; comments welcome

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

Optimizing Likelihood-free Inference using Self-supervised Neural Symmetry Embeddings

Authors: Deep Chatterjee, Philip C. Harris, Maanas Goel, Malina Desai, Michael W. Coughlin, Erik Katsavounidis

Abstract: Likelihood-free inference is quickly emerging as a powerful tool to perform fast/effective parameter estimation. We demonstrate a technique of optimizing likelihood-free inference to make it even faster by marginalizing symmetries in a physical problem. In this approach, physical symmetries, for example, time-translation are learned using joint-embedding via self-supervised learning with symmetry… ▽ More Likelihood-free inference is quickly emerging as a powerful tool to perform fast/effective parameter estimation. We demonstrate a technique of optimizing likelihood-free inference to make it even faster by marginalizing symmetries in a physical problem. In this approach, physical symmetries, for example, time-translation are learned using joint-embedding via self-supervised learning with symmetry data augmentations. Subsequently, parameter inference is performed using a normalizing flow where the embedding network is used to summarize the data before conditioning the parameters. We present this approach on two simple physical problems and we show faster convergence in a smaller number of parameters compared to a normalizing flow that does not use a pre-trained symmetry-informed representation. △ Less

Submitted 11 December, 2023; originally announced December 2023.

Comments: Accepted for Machine Learning and the Physical Sciences Workshop (submission 69) at NeurIPS 2023; for codes, see https://github.com/ML4GW/summer-projects-2023/blob/neurips-2023/symmetry-informed-flows/README.md

The ZTF Source Classification Project: III. A Catalog of Variable Sources

Authors: Brian F. Healy, Michael W. Coughlin, Ashish A. Mahabal, Theophile Jegou du Laz, Andrew Drake, Matthew J. Graham, Lynne A. Hillenbrand, Jan van Roestel, Paula Szkody, LeighAnna Zielske, Mohammed Guiga, Muhammad Yusuf Hassan, Jill L. Hughes, Guy Nir, Saagar Parikh, Sungmin Park, Palak Purohit, Umaa Rebbapragada, Draco Reed, Daniel Warshofsky, Avery Wold, Joshua S. Bloom, Frank J. Masci, Reed Riddle, Roger Smith

Abstract: The classification of variable objects provides insight into a wide variety of astrophysics ranging from stellar interiors to galactic nuclei. The Zwicky Transient Facility (ZTF) provides time series observations that record the variability of more than a billion sources. The scale of these data necessitates automated approaches to make a thorough analysis. Building on previous work, this paper re… ▽ More The classification of variable objects provides insight into a wide variety of astrophysics ranging from stellar interiors to galactic nuclei. The Zwicky Transient Facility (ZTF) provides time series observations that record the variability of more than a billion sources. The scale of these data necessitates automated approaches to make a thorough analysis. Building on previous work, this paper reports the results of the ZTF Source Classification Project (SCoPe), which trains neural network and XGBoost machine learning (ML) algorithms to perform dichotomous classification of variable ZTF sources using a manually constructed training set containing 170,632 light curves. We find that several classifiers achieve high precision and recall scores, suggesting the reliability of their predictions for 209,991,147 light curves across 77 ZTF fields. We also identify the most important features for XGB classification and compare the performance of the two ML algorithms, finding a pattern of higher precision among XGB classifiers. The resulting classification catalog is available to the public, and the software developed for SCoPe is open-source and adaptable to future time-domain surveys. △ Less

Submitted 11 April, 2024; v1 submitted 30 November, 2023; originally announced December 2023.

Comments: 23 pages, 12 figures, accepted to ApJS. Catalog available at https://zenodo.org/doi/10.5281/zenodo.8410825

Minutes-duration Optical Flares with Supernova Luminosities

Authors: Anna Y. Q. Ho, Daniel A. Perley, Ping Chen, Steve Schulze, Vik Dhillon, Harsh Kumar, Aswin Suresh, Vishwajeet Swain, Michael Bremer, Stephen J. Smartt, Joseph P. Anderson, G. C. Anupama, Supachai Awiphan, Sudhanshu Barway, Eric C. Bellm, Sagi Ben-Ami, Varun Bhalerao, Thomas de Boer, Thomas G. Brink, Rick Burruss, Poonam Chandra, Ting-Wan Chen, Wen-Ping Chen, Jeff Cooke, Michael W. Coughlin , et al. (52 additional authors not shown)

Abstract: In recent years, certain luminous extragalactic optical transients have been observed to last only a few days. Their short observed duration implies a different powering mechanism from the most common luminous extragalactic transients (supernovae) whose timescale is weeks. Some short-duration transients, most notably AT2018cow, display blue optical colours and bright radio and X-ray emission. Seve… ▽ More In recent years, certain luminous extragalactic optical transients have been observed to last only a few days. Their short observed duration implies a different powering mechanism from the most common luminous extragalactic transients (supernovae) whose timescale is weeks. Some short-duration transients, most notably AT2018cow, display blue optical colours and bright radio and X-ray emission. Several AT2018cow-like transients have shown hints of a long-lived embedded energy source, such as X-ray variability, prolonged ultraviolet emission, a tentative X-ray quasiperiodic oscillation, and large energies coupled to fast (but subrelativistic) radio-emitting ejecta. Here we report observations of minutes-duration optical flares in the aftermath of an AT2018cow-like transient, AT2022tsd (the "Tasmanian Devil"). The flares occur over a period of months, are highly energetic, and are likely nonthermal, implying that they arise from a near-relativistic outflow or jet. Our observations confirm that in some AT2018cow-like transients the embedded energy source is a compact object, either a magnetar or an accreting black hole. △ Less

Submitted 16 November, 2023; originally announced November 2023.

Comments: 79 pages, 3 figures (main text) + 7 figures (extended data) + 2 figures (supplementary information). Published online in Nature on 15 November 2023

Bayesian real-time classification of multi-messenger electromagnetic and gravitational-wave observations

Authors: Marina Berbel, Miquel Miravet-Tenés, Sushant Sharma Chaudhary, Simone Albanesi, Marco Cavaglià, Lorena Magaña Zertuche, Dimitra Tseneklidou, Yanyan Zheng, Michael W. Coughlin, Andrew Toivonen

Abstract: Because of the electromagnetic radiation produced during the merger, compact binary coalescences with neutron stars may result in multi-messenger observations. In order to follow up on the gravitational-wave signal with electromagnetic telescopes, it is critical to promptly identify the properties of these sources. This identification must rely on the properties of the progenitor source, such as t… ▽ More Because of the electromagnetic radiation produced during the merger, compact binary coalescences with neutron stars may result in multi-messenger observations. In order to follow up on the gravitational-wave signal with electromagnetic telescopes, it is critical to promptly identify the properties of these sources. This identification must rely on the properties of the progenitor source, such as the component masses and spins, as determined by low-latency detection pipelines in real time. The output of these pipelines, however, might be biased, which could decrease the accuracy of parameter recovery. Machine learning algorithms are used to correct this bias. In this work, we revisit this problem and discuss two new implementations of supervised machine learning algorithms, K-Nearest Neighbors and Random Forest, which are able to predict the presence of a neutron star and post-merger matter remnant in low-latency compact binary coalescence searches across different search pipelines and data sets. Additionally, we present a novel approach for calculating the Bayesian probabilities for these two metrics. Instead of metric scores derived from binary machine learning classifiers, our scheme is designed to provide the astronomy community well-defined probabilities. This would deliver a more direct and easily interpretable product to assist electromagnetic telescopes in deciding whether to follow up on gravitational-wave events in real time. △ Less

Submitted 31 October, 2023; originally announced November 2023.

Comments: 13 pages, 8 figures

Ready for O4 II: GRANDMA Observations of Swift GRBs during eight-weeks of Spring 2022

Authors: I. Tosta e Melo, J. -G. Ducoin, Z. Vidadi, C. Andrade, V. Rupchandani, S. Agayeva, J. Abdelhadi, L. Abe, O. Aguerre-Chariol, V. Aivazyan, S. Alishov, S. Antier, J. -M. Bai, A. Baransky, S. Bednarz, Ph. Bendjoya, Z. Benkhaldoun, S. Beradze, M. A. Bizouard, U. Bhardwaj, M. Blazek, M. Boër, E. Broens, O. Burkhonov, N. Christensen , et al. (84 additional authors not shown)

Abstract: We present a campaign designed to train the GRANDMA network and its infrastructure to follow up on transient alerts and detect their early afterglows. In preparation for O4 II campaign, we focused on GRB alerts as they are expected to be an electromagnetic counterpart of gravitational-wave events. Our goal was to improve our response to the alerts and start prompt observations as soon as possible… ▽ More We present a campaign designed to train the GRANDMA network and its infrastructure to follow up on transient alerts and detect their early afterglows. In preparation for O4 II campaign, we focused on GRB alerts as they are expected to be an electromagnetic counterpart of gravitational-wave events. Our goal was to improve our response to the alerts and start prompt observations as soon as possible to better prepare the GRANDMA network for the fourth observational run of LIGO-Virgo-Kagra (which started at the end of May 2023), and future missions such as SM. To receive, manage and send out observational plans to our partner telescopes we set up dedicated infrastructure and a rota of follow-up adcates were organized to guarantee round-the-clock assistance to our telescope teams. To ensure a great number of observations, we focused on Swift GRBs whose localization errors were generally smaller than the GRANDMA telescopes' field of view. This allowed us to bypass the transient identification process and focus on the reaction time and efficiency of the network. During 'Ready for O4 II', 11 Swift/INTEGRAL GRB triggers were selected, nine fields had been observed, and three afterglows were detected (GRB 220403B, GRB 220427A, GRB 220514A), with 17 GRANDMA telescopes and 17 amateur astronomers from the citizen science project Kilonova-Catcher. Here we highlight the GRB 220427A analysis where our long-term follow-up of the host galaxy allowed us to obtain a photometric redshift of $z=0.82\pm0.09$, its lightcurve elution, fit the decay slope of the afterglows, and study the properties of the host galaxy. △ Less

Submitted 26 October, 2023; originally announced October 2023.

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.

The first systematically identified repeating partial tidal disruption event

Authors: Jean J. Somalwar, Vikram Ravi, Yuhan Yao, Muryel Guolo, Matthew Graham, Erica Hammerstein, Wenbin Lu, Matt Nicholl, Yashvi Sharma, Robert Stein, Sjoert van Velzen, Eric C. Bellm, Michael W. Coughlin, Steven L. Groom, Frank J. Masci, Reed Riddle

Abstract: Tidal disruption events (TDEs) occur when a star enters the tidal radius of a supermassive black hole (SMBH). If the star only grazes the tidal radius, a fraction of the stellar mass will be accreted in a partial TDE (pTDE). The remainder can continue orbiting and may re-disrupted at pericenter, causing a repeating pTDE. pTDEs may be as or more common than full TDEs (fTDEs), yet few are known. In… ▽ More Tidal disruption events (TDEs) occur when a star enters the tidal radius of a supermassive black hole (SMBH). If the star only grazes the tidal radius, a fraction of the stellar mass will be accreted in a partial TDE (pTDE). The remainder can continue orbiting and may re-disrupted at pericenter, causing a repeating pTDE. pTDEs may be as or more common than full TDEs (fTDEs), yet few are known. In this work, we present the discovery of the first repeating pTDE from a systematically-selected sample, AT\,2020vdq. AT\,2020vdq was originally identified as an optically- and radio-flaring TDE. Around $3$ years after its discovery, it rebrightened dramatically and rapidly in the optical. The optical flare was remarkably fast and luminous compared to previous TDEs. It was accompanied by extremely broad (${\sim}0.1c$) optical/UV spectral features and faint X-ray emission ($L_X \sim 3\times10^{41}$\,erg\,s$^{-1}$), but no new radio-emitting component. Based on the transient optical/UV spectral features and the broadband light curve, we show that AT\,2020vdq is a repeating pTDE. We then use it to constrain TDE models; in particular, we favor a star originally in a very tight binary system that is tidally broken apart by the Hills mechanism. We also constrain the repeating pTDE rate to be $10^{-6}$ to $10^{-5}$ yr$^{-1}$ galaxy$^{-1}$, with uncertainties dominated by the unknown distribution of pTDE repeat timescales. In the Hills framework, this means the binary fraction in the galactic nucleus is of the order few percent. △ Less

Submitted 5 October, 2023; originally announced October 2023.

Comments: 24 pages, 13 figures, submitted to ApJ

GWAK: Gravitational-Wave Anomalous Knowledge with Recurrent Autoencoders

Authors: Ryan Raikman, Eric A. Moreno, Ekaterina Govorkova, Ethan J Marx, Alec Gunny, William Benoit, Deep Chatterjee, Rafia Omer, Muhammed Saleem, Dylan S Rankin, Michael W Coughlin, Philip C Harris, Erik Katsavounidis

Abstract: Matched-filtering detection techniques for gravitational-wave (GW) signals in ground-based interferometers rely on having well-modeled templates of the GW emission. Such techniques have been traditionally used in searches for compact binary coalescences (CBCs), and have been employed in all known GW detections so far. However, interesting science cases aside from compact mergers do not yet have ac… ▽ More Matched-filtering detection techniques for gravitational-wave (GW) signals in ground-based interferometers rely on having well-modeled templates of the GW emission. Such techniques have been traditionally used in searches for compact binary coalescences (CBCs), and have been employed in all known GW detections so far. However, interesting science cases aside from compact mergers do not yet have accurate enough modeling to make matched filtering possible, including core-collapse supernovae and sources where stochasticity may be involved. Therefore the development of techniques to identify sources of these types is of significant interest. In this paper, we present a method of anomaly detection based on deep recurrent autoencoders to enhance the search region to unmodeled transients. We use a semi-supervised strategy that we name Gravitational Wave Anomalous Knowledge (GWAK). While the semi-supervised nature of the problem comes with a cost in terms of accuracy as compared to supervised techniques, there is a qualitative advantage in generalizing experimental sensitivity beyond pre-computed signal templates. We construct a low-dimensional embedded space using the GWAK method, capturing the physical signatures of distinct signals on each axis of the space. By introducing signal priors that capture some of the salient features of GW signals, we allow for the recovery of sensitivity even when an unmodeled anomaly is encountered. We show that regions of the GWAK space can identify CBCs, detector glitches and also a variety of unmodeled astrophysical sources. △ Less

Submitted 20 September, 2023; originally announced September 2023.

An Optically-Discovered Outburst from XTE J1859+226

Authors: Eric C. Bellm, Yuankun Wang, Jan van Roestel, Rebecca A. Phillipson, Michael W. Coughlin, John A. Tomsick, Steven L. Groom, Brian Healy, Josiah Purdum, Ben Rusholme, Jesper Sollerman, Peter Bealo, Stefano Lora, Eddy Muyllaert, Ivo Peretto, Erik J. Schwendeman

Abstract: Using the Zwicky Transient Facility, in 2021 February we identified the first known outburst of the Black Hole X-ray Transient XTE J1859+226 since its discovery in 1999. The outburst was visible at X-ray, UV, and optical wavelengths for less than 20 days, substantially shorter than its 320-day full outburst in 1999, and the observed peak luminosity was two orders of magnitude lower. Its peak bolom… ▽ More Using the Zwicky Transient Facility, in 2021 February we identified the first known outburst of the Black Hole X-ray Transient XTE J1859+226 since its discovery in 1999. The outburst was visible at X-ray, UV, and optical wavelengths for less than 20 days, substantially shorter than its 320-day full outburst in 1999, and the observed peak luminosity was two orders of magnitude lower. Its peak bolometric luminosity was only $2\times 10^{35}$ erg s$^{-1}$, implying an Eddington fraction of about $3\times10^{-4}$. The source remained in the hard spectral state throughout the outburst. From optical spectroscopy measurements we estimate an outer disk radius of 10$^{11}$ cm. The low observed X-ray luminosity is not sufficient to irradiate the entire disk, but we observe a surprising exponential decline in the X-ray lightcurve. These observations highlight the potential of optical and infrared (O/IR) synoptic surveys to discover low-luminosity activity from X-ray transients. △ Less

Submitted 19 September, 2023; originally announced September 2023.

Comments: 12 pages, 6 figures, accepted for publication in ApJ

The Progenitor Star of SN 2023ixf: A Massive Red Supergiant with Enhanced, Episodic Pre-Supernova Mass Loss

Authors: Yu-Jing Qin, Keming Zhang, Joshua Bloom, Jesper Sollerman, Erez A. Zimmerman, Ido Irani, Steve Schulze, Avishay Gal-Yam, Mansi Kasliwal, Michael W. Coughlin, Daniel A. Perley, Christoffer Fremling, Shrinivas Kulkarni

Abstract: We identify the progenitor star of SN 2023ixf in the nearby galaxy Messier 101 using Keck/NIRC2 adaptive optics imaging and pre-explosion HST/ACS images. The supernova position, localized with diffraction-spike pattern and high precision relative astrometry, unambiguously coincides with a single progenitor candidate of m_F814W=24.96(-0.04)(+0.05). Forced photometry further recovers 2-sigma detecti… ▽ More We identify the progenitor star of SN 2023ixf in the nearby galaxy Messier 101 using Keck/NIRC2 adaptive optics imaging and pre-explosion HST/ACS images. The supernova position, localized with diffraction-spike pattern and high precision relative astrometry, unambiguously coincides with a single progenitor candidate of m_F814W=24.96(-0.04)(+0.05). Forced photometry further recovers 2-sigma detections in the F673N and F675W bands and imposes robust flux limits on the bluer bands. Given the reported infrared excess and semi-regular variability of the progenitor, we fit a time-dependent spectral energy distribution (SED) model of a dusty red supergiant (RSG) to a combined dataset of HST photometry, as well as ground-based near-infrared and Spitzer/IRAC [3.6], [4.5] photometry from the literature. The progenitor closely resembles a RSG of T_eff=3343+/-27 K and logL=5.10+/-0.02, with a 0.11+/-0.01 dex (25.2+/-1.7 per cent) variation over the mean luminosity at a period of P=1128.3+/-6.5 days, heavily obscured by a dust envelope with an optical depth of tau=2.83+/-0.03 at 1 micron (or A_V=10.28+/-0.11 mag). Such observed signatures match a post-main sequence star of 18.1(-1.2)(+0.7) Msun, close to the most massive SN II progenitor, with a pulsation-enhanced mass-loss rate of M_dot=(3.58+/-0.15) x 10^(-4) Msun/yr. The dense and confined circumstellar material is likely ejected during the last episode of radial pulsation before the explosion. Notably, we find strong evidence for periodic variation of tau (or both T_eff and tau) along with luminosity, a necessary assumption to reproduce the wavelength dependence of the variability, which implies dust sublimation and condensation during radial pulsations. Given the observed SED, partial dust obscuration remains a possible scenario, but any unobstructed binary companion over 7.1 Msun can be ruled out. △ Less

Submitted 18 September, 2023; originally announced September 2023.

Comments: submitted

A Joint Fermi-GBM and Swift-BAT Analysis of Gravitational-Wave Candidates from the Third Gravitational-wave Observing Run

Authors: C. Fletcher, J. Wood, R. Hamburg, P. Veres, C. M. Hui, E. Bissaldi, M. S. Briggs, E. Burns, W. H. Cleveland, M. M. Giles, A. Goldstein, B. A. Hristov, D. Kocevski, S. Lesage, B. Mailyan, C. Malacaria, S. Poolakkil, A. von Kienlin, C. A. Wilson-Hodge, The Fermi Gamma-ray Burst Monitor Team, M. Crnogorčević, J. DeLaunay, A. Tohuvavohu, R. Caputo, S. B. Cenko , et al. (1674 additional authors not shown)

Abstract: We present Fermi Gamma-ray Burst Monitor (Fermi-GBM) and Swift Burst Alert Telescope (Swift-BAT) searches for gamma-ray/X-ray counterparts to gravitational wave (GW) candidate events identified during the third observing run of the Advanced LIGO and Advanced Virgo detectors. Using Fermi-GBM on-board triggers and sub-threshold gamma-ray burst (GRB) candidates found in the Fermi-GBM ground analyses,… ▽ More We present Fermi Gamma-ray Burst Monitor (Fermi-GBM) and Swift Burst Alert Telescope (Swift-BAT) searches for gamma-ray/X-ray counterparts to gravitational wave (GW) candidate events identified during the third observing run of the Advanced LIGO and Advanced Virgo detectors. Using Fermi-GBM on-board triggers and sub-threshold gamma-ray burst (GRB) candidates found in the Fermi-GBM ground analyses, the Targeted Search and the Untargeted Search, we investigate whether there are any coincident GRBs associated with the GWs. We also search the Swift-BAT rate data around the GW times to determine whether a GRB counterpart is present. No counterparts are found. Using both the Fermi-GBM Targeted Search and the Swift-BAT search, we calculate flux upper limits and present joint upper limits on the gamma-ray luminosity of each GW. Given these limits, we constrain theoretical models for the emission of gamma-rays from binary black hole mergers. △ Less

Submitted 25 August, 2023; originally announced August 2023.

Low-latency gravitational wave alert products and their performance at the time of the fourth LIGO-Virgo-KAGRA observing run

Authors: Sushant Sharma Chaudhary, Andrew Toivonen, Gaurav Waratkar, Geoffrey Mo, Deep Chatterjee, Sarah Antier, Patrick Brockill, Michael W. Coughlin, Reed Essick, Shaon Ghosh, Soichiro Morisaki, Pratyusava Baral, Amanda Baylor, Naresh Adhikari, Patrick Brady, Gareth Cabourn Davies, Tito Dal Canton, Marco Cavaglià, Jolien Creighton, Sunil Choudhary, Yu-Kuang Chu, Patrick Clearwater, Luke Davis, Thomas Dent, Marco Drago , et al. (28 additional authors not shown)

Abstract: Multi-messenger searches for BNS and NSBH mergers are currently one of the most exciting areas of astronomy. The search for joint electromagnetic and neutrino counterparts to GWs has resumed with O4. To support this effort, public semi-automated data products are sent in near real-time and include localization and source properties to guide complementary observations. In preparation for O4, we hav… ▽ More Multi-messenger searches for BNS and NSBH mergers are currently one of the most exciting areas of astronomy. The search for joint electromagnetic and neutrino counterparts to GWs has resumed with O4. To support this effort, public semi-automated data products are sent in near real-time and include localization and source properties to guide complementary observations. In preparation for O4, we have conducted a study using a simulated population of compact binaries and a MDC in the form of a real-time replay to optimize and profile the software infrastructure and scientific deliverables. End-to-end performance was tested, including data ingestion, running online search pipelines, performing annotations, and issuing alerts to the astrophysics community. We present an overview of the low-latency infrastructure and the performance of the data products that are now being released during O4 based on the MDC. We report the expected median latency for the preliminary alert of full bandwidth searches (29.5s) and show consistency and accuracy of released data products using the MDC. For the first time, we report the expected median latency for triggers from early warning searches (-3.1s), which are new in O4 and target neutron star mergers during inspiral phase. This paper provides a performance overview for LVK low-latency alert infrastructure and data products using the MDC and serves as a useful reference for the interpretation of O4 detections. △ Less

Submitted 27 May, 2024; v1 submitted 8 August, 2023; originally announced August 2023.

Journal ref: PNAS 121 (18) e2316474121 (2024)

Search for Eccentric Black Hole Coalescences during the Third Observing Run of LIGO and Virgo

Authors: The LIGO Scientific Collaboration, the Virgo Collaboration, the KAGRA Collaboration, A. G. Abac, R. Abbott, H. Abe, F. Acernese, K. Ackley, C. Adamcewicz, S. Adhicary, N. Adhikari, R. X. Adhikari, V. K. Adkins, V. B. Adya, C. Affeldt, D. Agarwal, M. Agathos, O. D. Aguiar, I. Aguilar, L. Aiello, A. Ain, P. Ajith, T. Akutsu, S. Albanesi, R. A. Alfaidi , et al. (1750 additional authors not shown)

Abstract: Despite the growing number of confident binary black hole coalescences observed through gravitational waves so far, the astrophysical origin of these binaries remains uncertain. Orbital eccentricity is one of the clearest tracers of binary formation channels. Identifying binary eccentricity, however, remains challenging due to the limited availability of gravitational waveforms that include effect… ▽ More Despite the growing number of confident binary black hole coalescences observed through gravitational waves so far, the astrophysical origin of these binaries remains uncertain. Orbital eccentricity is one of the clearest tracers of binary formation channels. Identifying binary eccentricity, however, remains challenging due to the limited availability of gravitational waveforms that include effects of eccentricity. Here, we present observational results for a waveform-independent search sensitive to eccentric black hole coalescences, covering the third observing run (O3) of the LIGO and Virgo detectors. We identified no new high-significance candidates beyond those that were already identified with searches focusing on quasi-circular binaries. We determine the sensitivity of our search to high-mass (total mass $M>70$ $M_\odot$) binaries covering eccentricities up to 0.3 at 15 Hz orbital frequency, and use this to compare model predictions to search results. Assuming all detections are indeed quasi-circular, for our fiducial population model, we place an upper limit for the merger rate density of high-mass binaries with eccentricities $0 < e \leq 0.3$ at $0.33$ Gpc$^{-3}$ yr$^{-1}$ at 90\% confidence level. △ Less

Submitted 7 August, 2023; originally announced August 2023.

Comments: 24 pages, 5 figures

Report number: LIGO-P2300080

Chemical Distribution of the Dynamical Ejecta in the Neutron Star Merger GW170817

Authors: Shreya Anand, Peter T. H. Pang, Mattia Bulla, Michael W. Coughlin, Tim Dietrich, Brian Healy, Thomas Hussenot-Desenonges, Theophile Jegou du Laz, Mansi M. Kasliwal, Nina Kunert, Ivan Markin, Kunal Mooley, Vsevolod Nedora, Anna Neuweiler

Abstract: GW170817 and its associated electromagnetic counterpart AT2017gfo continue to be a treasure trove as observations and modeling continue. Recent precision astrometry of AT2017gfo with the Hubble Space Telescope combined with previous constraints from Very Long Baseline Interferometry (VLBI) constraints narrowed down the inclination angle to 19-25 deg (90\% confidence). This paper explores how the i… ▽ More GW170817 and its associated electromagnetic counterpart AT2017gfo continue to be a treasure trove as observations and modeling continue. Recent precision astrometry of AT2017gfo with the Hubble Space Telescope combined with previous constraints from Very Long Baseline Interferometry (VLBI) constraints narrowed down the inclination angle to 19-25 deg (90\% confidence). This paper explores how the inclusion of precise inclination information can reveal new insights about the ejecta properties, in particular, about the composition of the dynamical ejecta of AT2017gfo. Our analysis relies on updated kilonova modeling, which includes state-of-the-art heating rates, thermalization efficiencies, and opacities and is parameterized by $\bar{Y}_{\rm e,dyn}$, the average electron fraction of the dynamical ejecta component. Using this model, we incorporate the latest inclination angle constraint of AT2017gfo into a light curve fitting framework to derive updated parameter estimates. Our results suggest that the viewing angle of the observer is pointed towards the lanthanide-poor ($Y_{\rm e,dyn}\gtrsim0.25$), squeezed polar dynamical ejecta component, which can explain the early blue emission observed in the light curve of AT2017gfo. In contrast to a recent claim of spherical ejecta powering AT2017gfo, our study indicates that the composition of the dynamical ejecta has a strong angular dependence, with a lanthanide-rich ($Y_{\rm e,dyn}\lesssim0.25$), tidal component distributed around the merger plane with a half-opening angle of $35^\circ$. The inclination angle constraint reduces $\bar{Y}_{\rm e,dyn}$ from $0.24$ to $0.22$, with values $0.17\lesssim Y_{\rm e, dyn} \lesssim0.41$ enabling the robust production of $r$-process elements up to the $3^{\rm rd}$ peak in the tidal dynamical ejecta. △ Less

Submitted 20 July, 2023; originally announced July 2023.

Comments: 8 pages, 6 figures

Enabling Kilonova Science with Nancy Grace Roman Space Telescope

Authors: Igor Andreoni, Michael W. Coughlin, Alexander W. Criswell, Mattia Bulla, Andrew Toivonen, Leo P. Singer, Antonella Palmese, E. Burns, Suvi Gezari, Mansi M. Kasliwal, R. Weizmann Kiendrebeogo, Ashish Mahabal, Takashi J. Moriya, Armin Rest, Dan Scolnic, Robert A. Simcoe, Jamie Soon, Robert Stein, Tony Travouillon

Abstract: Binary neutron star mergers and neutron star-black hole mergers are multi-messenger sources that can be detected in gravitational waves and in electromagnetic radiation. The low electron fraction of neutron star merger ejecta favors the production of heavy elements such as lanthanides and actinides via rapid neutron capture (r-process). The decay of these unstable nuclei powers an infrared-bright… ▽ More Binary neutron star mergers and neutron star-black hole mergers are multi-messenger sources that can be detected in gravitational waves and in electromagnetic radiation. The low electron fraction of neutron star merger ejecta favors the production of heavy elements such as lanthanides and actinides via rapid neutron capture (r-process). The decay of these unstable nuclei powers an infrared-bright transient called a "kilonova". The discovery of a population of kilonovae will allow us to determine if neutron star mergers are the dominant sites for r-process element nucleosynthesis, constrain the equation of state of nuclear matter, and make independent measurements of the Hubble constant. The Nancy Grace Roman Space Telescope (Roman) will have a unique combination of depth, near-infrared sensitivity, and wide field of view. These characteristics will enable Roman's discovery of GW counterparts that will be missed by optical telescopes, such as kilonova that are associated with large distances, high lanthanide fractions, high binary mass-ratios, large dust extinction in the line of sight, or that are observed from equatorial viewing angles. Our analysis suggests to (i) make available a rapid (about 1 week) Target of Opportunity mode for GW follow-up; (ii) include observations of the High Latitude Time-Domain survey footprint in at least two filters (preferably the F158 and F213 filters) with a cadence of < 8 days; (iii) operate in synergy with Rubin Observatory. Following these recommendations, we expect that 1-6 kilonovae can be identified by Roman via ToO observations of well localized (A < 10 sq. deg., 90% C.I.) neutron star mergers during 1.5 years of the LIGO-Virgo-KAGRA fifth (or about 4-21 in during the sixth) observing run. A sample of 5-40 serendipitously discovered kilonovae can be collected in a 5-year high latitude survey. △ Less

Submitted 4 October, 2023; v1 submitted 18 July, 2023; originally announced July 2023.

Comments: accepted for publication in Astroparticle Physics

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!

$\texttt{BTSbot}$: A Multi-input Convolutional Neural Network to Automate and Expedite Bright Transient Identification for the Zwicky Transient Facility

Authors: Nabeel Rehemtulla, Adam A. Miller, Michael W. Coughlin, Theophile Jegou du Laz

Abstract: The Bright Transient Survey (BTS) relies on visual inspection ("scanning") to select sources for accomplishing its mission of spectroscopically classifying all bright extragalactic transients found by the Zwicky Transient Facility (ZTF). We present $\texttt{BTSbot}$, a multi-input convolutional neural network, which provides a bright transient score to individual ZTF detections using their image d… ▽ More The Bright Transient Survey (BTS) relies on visual inspection ("scanning") to select sources for accomplishing its mission of spectroscopically classifying all bright extragalactic transients found by the Zwicky Transient Facility (ZTF). We present $\texttt{BTSbot}$, a multi-input convolutional neural network, which provides a bright transient score to individual ZTF detections using their image data and 14 extracted features. $\texttt{BTSbot}$ eliminates the need for scanning by automatically identifying and requesting follow-up observations of new bright ($m\,<18.5\,\mathrm{mag}$) transient candidates. $\texttt{BTSbot}$ outperforms BTS scanners in terms of completeness (99% vs. 95%) and identification speed (on average, 7.4 hours quicker). △ Less

Submitted 14 July, 2023; originally announced July 2023.

Comments: Accepted at the ICML 2023 Workshop on Machine Learning for Astrophysics

Early Spectroscopy and Dense Circumstellar Medium Interaction in SN 2023ixf

Authors: K. Azalee Bostroem, Jeniveve Pearson, Manisha Shrestha, David J. Sand, Stefano Valenti, Saurabh W. Jha, Jennifer E. Andrews, Nathan Smith, Giacomo Terreran, Elizabeth Green, Yize Dong, Michael Lundquist, Joshua Haislip, Emily T. Hoang, Griffin Hosseinzadeh, Daryl Janzen, Jacob E. Jencson, Vladimir Kouprianov, Emmy Paraskeva, Nicolas E. Meza Retamal, Daniel E. Reichart, Iair Arcavi, Alceste Z. Bonanos, Michael W. Coughlin, Ross Dobson , et al. (31 additional authors not shown)

Abstract: We present the optical spectroscopic evolution of SN~2023ixf seen in sub-night cadence spectra from 1.18 to 14 days after explosion. We identify high-ionization emission features, signatures of interaction with material surrounding the progenitor star, that fade over the first 7 days, with rapid evolution between spectra observed within the same night. We compare the emission lines present and the… ▽ More We present the optical spectroscopic evolution of SN~2023ixf seen in sub-night cadence spectra from 1.18 to 14 days after explosion. We identify high-ionization emission features, signatures of interaction with material surrounding the progenitor star, that fade over the first 7 days, with rapid evolution between spectra observed within the same night. We compare the emission lines present and their relative strength to those of other supernovae with early interaction, finding a close match to SN~2020pni and SN~2017ahn in the first spectrum and SN~2014G at later epochs. To physically interpret our observations we compare them to CMFGEN models with confined, dense circumstellar material around a red supergiant progenitor from the literature. We find that very few models reproduce the blended \NC{} emission lines observed in the first few spectra and their rapid disappearance thereafter, making this a unique diagnostic. From the best models, we find a mass-loss rate of $10^{-3}-10^{-2}$ \mlunit{}, which far exceeds the mass-loss rate for any steady wind, especially for a red supergiant in the initial mass range of the detected progenitor. These mass-loss rates are, however, similar to rates inferred for other supernovae with early circumstellar interaction. Using the phase when the narrow emission features disappear, we calculate an outer dense radius of circumstellar material $R_\mathrm{CSM, out}\sim5\times10^{14}~\mathrm{cm}$ and a mean circumstellar material density of $ρ=5.6\times10^{-14}~\mathrm{g\,cm^{-3}}$. This is consistent with the lower limit on the outer radius of the circumstellar material we calculate from the peak \Halpha{} emission flux, $R_\text{CSM, out}\gtrsim9\times10^{13}~\mathrm{cm}$. △ Less

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

Comments: Published in ApJL

Journal ref: The Astrophysical Journal Letters, Volume 956, Issue 1, id.L5, 17 pp., Oct 2023

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

QoQ: a Q-transform based test for Gravitational Wave transient events

Authors: Siddharth Soni, Ethan Marx, Erik Katsavounidis, Reed Essick, G. S. Cabourn Davies, Patrick Brockill, Michael W. Coughlin, Shaon Ghosh, Patrick Godwin

Abstract: The observation of transient gravitational waves is hindered by the presence of transient noise, colloquially referred to as glitches. These glitches can often be misidentified as gravitational waves by searches for unmodeled transients using the excess-power type of methods and sometimes even excite template waveforms for compact binary coalescences while using matched filter techniques. They thu… ▽ More The observation of transient gravitational waves is hindered by the presence of transient noise, colloquially referred to as glitches. These glitches can often be misidentified as gravitational waves by searches for unmodeled transients using the excess-power type of methods and sometimes even excite template waveforms for compact binary coalescences while using matched filter techniques. They thus create a significant background in the searches. This background is more critical in getting identified promptly and efficiently within the context of real-time searches for gravitational-wave transients. Such searches are the ones that have enabled multi-messenger astrophysics with the start of the Advanced LIGO and Advanced Virgo data taking in 2015 and they will continue to enable the field for further discoveries. With this work we propose and demonstrate the use of a signal-based test that quantifies the fidelity of the time-frequency decomposition of the putative signal based on first principles on how astrophysical transients are expected to be registered in the detectors and empirically measuring the instrumental noise. It is based on the Q-transform and a measure of the occupancy of the corresponding time-frequency pixels over select time-frequency volumes; we call it ``QoQ''. Our method shows a 40% reduction in the number of retraction of public alerts that were issued by the LIGO-Virgo-KAGRA collaborations during the third observing run with negligible loss in sensitivity. Receiver Operator Characteristic measurements suggest the method can be used in online and offline searches for transients, reducing their background significantly. △ Less

Submitted 14 May, 2023; originally announced May 2023.

Comments: 39 Figures, 5 Tables

Performance of the low-latency GstLAL inspiral search towards LIGO, Virgo, and KAGRA's fourth observing run

Authors: Becca Ewing, Rachael Huxford, Divya Singh, Leo Tsukada, Chad Hanna, Yun-Jing Huang, Prathamesh Joshi, Alvin K. Y. Li, Ryan Magee, Cody Messick, Alex Pace, Anarya Ray, Surabhi Sachdev, Shio Sakon, Ron Tapia, Shomik Adhicary, Pratyusava Baral, Amanda Baylor, Kipp Cannon, Sarah Caudill, Sushant Sharma Chaudhary, Michael W. Coughlin, Bryce Cousins, Jolien D. E. Creighton, Reed Essick , et al. (18 additional authors not shown)

Abstract: GstLAL is a stream-based matched-filtering search pipeline aiming at the prompt discovery of gravitational waves from compact binary coalescences such as the mergers of black holes and neutron stars. Over the past three observation runs by the LIGO, Virgo, and KAGRA (LVK) collaboration, the GstLAL search pipeline has participated in several tens of gravitational wave discoveries. The fourth observ… ▽ More GstLAL is a stream-based matched-filtering search pipeline aiming at the prompt discovery of gravitational waves from compact binary coalescences such as the mergers of black holes and neutron stars. Over the past three observation runs by the LIGO, Virgo, and KAGRA (LVK) collaboration, the GstLAL search pipeline has participated in several tens of gravitational wave discoveries. The fourth observing run (O4) is set to begin in May 2023 and is expected to see the discovery of many new and interesting gravitational wave signals which will inform our understanding of astrophysics and cosmology. We describe the current configuration of the GstLAL low-latency search and show its readiness for the upcoming observation run by presenting its performance on a mock data challenge. The mock data challenge includes 40 days of LIGO Hanford, LIGO Livingston, and Virgo strain data along with an injection campaign in order to fully characterize the performance of the search. We find an improved performance in terms of detection rate and significance estimation as compared to that observed in the O3 online analysis. The improvements are attributed to several incremental advances in the likelihood ratio ranking statistic computation and the method of background estimation. △ Less

Submitted 13 July, 2023; v1 submitted 9 May, 2023; originally announced May 2023.

Comments: 19 pages, 21 figures

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

Early-time spectroscopic modelling of the transitional Type Ia Supernova 2021rhu with TARDIS

Authors: Luke Harvey, Kate Maguire, Mark R. Magee, Mattia Bulla, Suhail Dhawan, Steve Schulze, Jesper Sollerman, Maxime Deckers, Georgios Dimitriadis, Simeon Reusch, Mathew Smith, Jacco Terwel, Michael W. Coughlin, Frank Masci, Josiah Purdum, Alexander Reedy, Estelle Robert, Avery Wold

Abstract: An open question in SN Ia research is where the boundary lies between 'normal' Type Ia supernovae (SNe Ia) that are used in cosmological measurements and those that sit off the Phillips relation. We present the spectroscopic modelling of one such '86G-like' transitional SN Ia, SN 2021rhu, that has recently been employed as a local Hubble Constant calibrator using a tip of the red-giant branch meas… ▽ More An open question in SN Ia research is where the boundary lies between 'normal' Type Ia supernovae (SNe Ia) that are used in cosmological measurements and those that sit off the Phillips relation. We present the spectroscopic modelling of one such '86G-like' transitional SN Ia, SN 2021rhu, that has recently been employed as a local Hubble Constant calibrator using a tip of the red-giant branch measurement. We detail its modelling from -12 d until maximum brightness using the radiative-transfer spectral-synthesis code tardis. We base our modelling on literature delayed-detonation and deflagration models of Chandrasekhar mass white dwarfs, as well as the double-detonation models of sub-Chandrasekhar mass white dwarfs. We present a new method for 'projecting' abundance profiles to different density profiles for ease of computation. Due to the small velocity extent and low outer densities of the W7 profile, we find it inadequate to reproduce the evolution of SN 2021rhu as it fails to match the high-velocity calcium components. The host extinction of SN 2021rhu is uncertain but we use modelling with and without an extinction correction to set lower and upper limits on the abundances of individual species. Comparing these limits to literature models we conclude that the spectral evolution of SN 2021rhu is also incompatible with double-detonation scenarios, lying more in line with those resulting from the delayed detonation mechanism (although there are some discrepancies, in particular a larger titanium abundance in SN 2021rhu compared to the literature). This suggests that SN 2021rhu is likely a lower luminosity, and hence lower temperature, version of a normal SN Ia. △ Less

Submitted 21 April, 2023; v1 submitted 20 April, 2023; originally announced April 2023.

Comments: 25 pages, 22 figures, accepted for publication in MNRAS

Search for gravitational-lensing signatures in the full third observing run of the LIGO-Virgo network

Authors: The LIGO Scientific Collaboration, the Virgo Collaboration, the KAGRA Collaboration, R. Abbott, H. Abe, F. Acernese, K. Ackley, S. Adhicary, N. Adhikari, R. X. Adhikari, V. K. Adkins, V. B. Adya, C. Affeldt, D. Agarwal, M. Agathos, O. D. Aguiar, L. Aiello, A. Ain, P. Ajith, T. Akutsu, S. Albanesi, R. A. Alfaidi, C. Alléné, A. Allocca, P. A. Altin , et al. (1670 additional authors not shown)

Abstract: Gravitational lensing by massive objects along the line of sight to the source causes distortions of gravitational wave-signals; such distortions may reveal information about fundamental physics, cosmology and astrophysics. In this work, we have extended the search for lensing signatures to all binary black hole events from the third observing run of the LIGO--Virgo network. We search for repeated… ▽ More Gravitational lensing by massive objects along the line of sight to the source causes distortions of gravitational wave-signals; such distortions may reveal information about fundamental physics, cosmology and astrophysics. In this work, we have extended the search for lensing signatures to all binary black hole events from the third observing run of the LIGO--Virgo network. We search for repeated signals from strong lensing by 1) performing targeted searches for subthreshold signals, 2) calculating the degree of overlap amongst the intrinsic parameters and sky location of pairs of signals, 3) comparing the similarities of the spectrograms amongst pairs of signals, and 4) performing dual-signal Bayesian analysis that takes into account selection effects and astrophysical knowledge. We also search for distortions to the gravitational waveform caused by 1) frequency-independent phase shifts in strongly lensed images, and 2) frequency-dependent modulation of the amplitude and phase due to point masses. None of these searches yields significant evidence for lensing. Finally, we use the non-detection of gravitational-wave lensing to constrain the lensing rate based on the latest merger-rate estimates and the fraction of dark matter composed of compact objects. △ Less

Submitted 17 April, 2023; originally announced April 2023.

Comments: 28 pages, 11 figures

Report number: LIGO-P2200031

SN 2019odp: A Massive Oxygen-Rich Type Ib Supernova

Authors: T. Schweyer, J. Sollerman, A. Jerkstrand, M. Ergon, T. -W. Chen, C. M. B. Omand, S. Schulze, M. W. Coughlin, I. Andreoni, C. Fremling, A. Rau, Y. Sharma, N. L. Strotjohann, L. Yan, M. J. Graham, M. M. Kasliwal, R. R. Laher, J. Purdum, P. Rosnet, B. Rusholme, R. Smith

Abstract: We present and analyze observations of the Type Ib supernova (SN) 2019odp (a.k.a ZTF19abqwtfu) covering epochs within days of the explosion to the nebular phase at 360 d post-explosion. We discuss them in the context of recombination cooling emission for the early excess emission and consider progenitor models based on the nebular phase spectra. Our observations include photometric observations ma… ▽ More We present and analyze observations of the Type Ib supernova (SN) 2019odp (a.k.a ZTF19abqwtfu) covering epochs within days of the explosion to the nebular phase at 360 d post-explosion. We discuss them in the context of recombination cooling emission for the early excess emission and consider progenitor models based on the nebular phase spectra. Our observations include photometric observations mainly in the optical and low to medium-resolution spectroscopic observations covering the complete observable time-range. We expand on existing methods to derive oxygen mass estimates from nebular phase spectroscopy. Our spectroscopic observations confirm the presence of He in the SN ejecta and we thus (re)classify it as a Type Ib supernova. From the pseudo-bolometric lightcurve we estimate a high ejecta mass $M_\text{ej} \sim 4 - 7~M_\odot$. The high ejecta mass, large nebular [O I]/[Ca II] line flux ratio ($1.2 - 1.9$) and an oxygen mass above $\gtrapprox 0.5\, M_\odot$ point towards a progenitor with pre-explosion mass higher than $18\,M_\odot$. The compact nature of the progenitor ($\lesssim 10\,R_\odot$) suggests a Wolf-Rayet (WR) star as progenitor. △ Less

Submitted 24 March, 2023; originally announced March 2023.

Comments: Submitted to A&A

A Sensitive Search for Supernova Emission Associated with the Extremely Energetic and Nearby GRB 221009A

Authors: Gokul P. Srinivasaragavan, Brendan O' Connor, S. Bradley Cenko, Alexander J. Dittmann, Sheng Yang, Jesper Sollerman, G. C. Anupama, Sudhanshu Barway, Varun Bhalerao, Harsh Kumar, Vishwajeet Swain, Erica Hammerstein, Isiah Holt, Shreya Anand, Igor Andreoni, Michael W. Coughlin, Simone Dichiara, Avishay Gal-Yam, M. Coleman Miller, Jaime Soon, Roberto Soria, Joseph Durbak, James H. Gillanders, Sibasish Laha, Anna M. Moore , et al. (2 additional authors not shown)

Abstract: We report observations of the optical counterpart of the long gamma-ray burst (LGRB) GRB 221009A. Due to the extreme rarity of being both nearby ($z = 0.151$) and highly energetic ($E_{γ,\mathrm{iso}} \geq 10^{54}$ erg), GRB 221009A offers a unique opportunity to probe the connection between massive star core collapse and relativistic jet formation across a very broad range of $γ$-ray properties.… ▽ More We report observations of the optical counterpart of the long gamma-ray burst (LGRB) GRB 221009A. Due to the extreme rarity of being both nearby ($z = 0.151$) and highly energetic ($E_{γ,\mathrm{iso}} \geq 10^{54}$ erg), GRB 221009A offers a unique opportunity to probe the connection between massive star core collapse and relativistic jet formation across a very broad range of $γ$-ray properties. Adopting a phenomenological power-law model for the afterglow and host galaxy estimates from high-resolution Hubble Space Telescope imaging, we use Bayesian model comparison techniques to determine the likelihood of an associated SN contributing excess flux to the optical light curve. Though not conclusive, we find moderate evidence ($K_{\rm{Bayes}}=10^{1.2}$) for the presence of an additional component arising from an associated supernova, SN 2022xiw, and find that it must be substantially fainter ($<$ 67% as bright at the 99% confidence interval) than SN 1998bw. Given the large and uncertain line-of-sight extinction, we attempt to constrain the supernova parameters ($M_{\mathrm{Ni}}$, $M_{\mathrm{ej}}$, and $E_{\mathrm{KE}}$) under several different assumptions with respect to the host galaxy's extinction. We find properties that are broadly consistent with previous GRB-associated SNe: $M_{\rm{Ni}}=0.05$ - $0.25 \, \rm{M_\odot}$, $M_{\rm{ej}}=3.5$ - $11.1 \, \rm{M_\odot}$, and $E_{\rm{KE}} = (1.6$ - $5.2) \times 10^{52} \, \rm{erg}$. We note that these properties are weakly constrained due to the faintness of the supernova with respect to the afterglow and host emission, but we do find a robust upper limit on the $M_{\rm{Ni}}$ of $M_{\rm{Ni}}<0.36\, \rm{M_\odot}$. Given the tremendous range in isotropic gamma-ray energy release exhibited by GRBs (7 orders of magnitude), the SN emission appears to be decoupled from the central engine in these systems. △ Less

Submitted 28 May, 2023; v1 submitted 22 March, 2023; originally announced March 2023.

Comments: 18 pages, accepted to ApJL, 4 tables, 5 figures. Updated abstract in Preview

Insight-HXMT and GECAM-C observations of the brightest-of-all-time GRB 221009A

Authors: Zheng-Hua An, S. Antier, Xing-Zi Bi, Qing-Cui Bu, Ce Cai, Xue-Lei Cao, Anna-Elisa Camisasca, Zhi Chang, Gang Chen, Li Chen, Tian-Xiang Chen, Wen Chen, Yi-Bao Chen, Yong Chen, Yu-Peng Chen, Michael W. Coughlin, Wei-Wei Cui, Zi-Gao Dai, T. Hussenot-Desenonges, Yan-Qi Du, Yuan-Yuan Du, Yun-Fei Du, Cheng-Cheng Fan, Filippo Frontera, He Gao , et al. (153 additional authors not shown)

Abstract: GRB 221009A is the brightest gamma-ray burst ever detected since the discovery of this kind of energetic explosions. However, an accurate measurement of the prompt emission properties of this burst is very challenging due to its exceptional brightness. With joint observations of \textit{Insight}-HXMT and GECAM-C, we made an unprecedentedly accurate measurement of the emission during the first… ▽ More GRB 221009A is the brightest gamma-ray burst ever detected since the discovery of this kind of energetic explosions. However, an accurate measurement of the prompt emission properties of this burst is very challenging due to its exceptional brightness. With joint observations of \textit{Insight}-HXMT and GECAM-C, we made an unprecedentedly accurate measurement of the emission during the first $\sim$1800 s of GRB 221009A, including its precursor, main emission (ME, which dominates the burst in flux), flaring emission and early afterglow, in the hard X-ray to soft gamma-ray band from $\sim$ 10 keV to $\sim$ 6 MeV. Based on the GECAM-C unsaturated data of the ME, we measure a record-breaking isotropic equivalent energy ($E_{\rm iso}$) of $\bf \sim 1.5 \times 10^{55}$ erg, which is about eight times the total rest-mass energy of the Sun. The early afterglow data require a significant jet break between 650 s and 1100 s, most likely at $\sim950$ s from the afterglow starting time $T_{AG}$, which corresponds to a jet opening angle of $\sim {0.7^\circ} \ (η_γn)^{1/8}$, where $n$ is the ambient medium density in units of $\rm cm^{-3}$ and $η_γ$ is the ratio between $γ$-ray energy and afterglow kinetic energy. The beaming-corrected total $γ$-ray energy $E_γ$ is $\sim 1.15 \times10^{51} \ (η_γn)^{1/4}$ erg, which is typical for long GRBs. These results suggest that this GRB may have a special central engine, which could launch and collimate a very narrowly beamed jet with an ordinary energy budget, leading to exceptionally luminous gamma-ray radiation per unit solid angle. Alternatively, more GRBs might have such a narrow and bright beam, which are missed by an unfavorable viewing angle or have been detected without distance measurement. △ Less

Submitted 3 March, 2023; v1 submitted 2 March, 2023; originally announced March 2023.

Comments: Submitted to National Science Review. This paper is under press embargo, contact the corresponding author for details

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)

GRANDMA and HXMT Observations of GRB 221009A -- the Standard-Luminosity Afterglow of a Hyper-Luminous Gamma-Ray Burst

Authors: D. A. Kann, S. Agayeva, V. Aivazyan, S. Alishov, C. M. Andrade, S. Antier, A. Baransky, P. Bendjoya, Z. Benkhaldoun, S. Beradze, D. Berezin, M. Boër, E. Broens, S. Brunier, M. Bulla, O. Burkhonov, E. Burns, Y. Chen, Y. P. Chen, M. Conti, M. W. Coughlin, W. W. Cui, F. Daigne, B. Delaveau, H. A. R. Devillepoix , et al. (91 additional authors not shown)

Abstract: GRB 221009A is the brightest Gamma-Ray Burst (GRB) detected in more than 50 years of study. In this paper, we present observations in the X-ray and optical domains after the GRB obtained by the GRANDMA Collaboration (which includes observations from more than 30 professional and amateur telescopes) and the Insight-HXMT Collaboration. We study the optical afterglow with empirical fitting from GRAND… ▽ More GRB 221009A is the brightest Gamma-Ray Burst (GRB) detected in more than 50 years of study. In this paper, we present observations in the X-ray and optical domains after the GRB obtained by the GRANDMA Collaboration (which includes observations from more than 30 professional and amateur telescopes) and the Insight-HXMT Collaboration. We study the optical afterglow with empirical fitting from GRANDMA+HXMT data, augmented with data from the literature up to 60 days. We then model numerically, using a Bayesian approach, the GRANDMA and HXMT-LE afterglow observations, that we augment with Swift-XRT and additional optical/NIR observations reported in the literature. We find that the GRB afterglow, extinguished by a large dust column, is most likely behind a combination of a large Milky-Way dust column combined with moderate low-metallicity dust in the host galaxy. Using the GRANDMA+HXMT-LE+XRT dataset, we find that the simplest model, where the observed afterglow is produced by synchrotron radiation at the forward external shock during the deceleration of a top-hat relativistic jet by a uniform medium, fits the multi-wavelength observations only moderately well, with a tension between the observed temporal and spectral evolution. This tension is confirmed when using the extended dataset. We find that the consideration of a jet structure (Gaussian or power-law), the inclusion of synchrotron self-Compton emission, or the presence of an underlying supernova do not improve the predictions, showing that the modelling of GRB22109A will require going beyond the most standard GRB afterglow model. Placed in the global context of GRB optical afterglows, we find the afterglow of GRB 221009A is luminous but not extraordinarily so, highlighting that some aspects of this GRB do not deviate from the global known sample despite its extreme energetics and the peculiar afterglow evolution. △ Less

Submitted 27 March, 2023; v1 submitted 13 February, 2023; originally announced February 2023.

Comments: Accepted to ApJL for the special issue, 37 pages, 23 pages main text, 6 tables, 13 figures

Open data from the third observing run of LIGO, Virgo, KAGRA and GEO

Authors: The LIGO Scientific Collaboration, the Virgo Collaboration, the KAGRA Collaboration, R. Abbott, H. Abe, F. Acernese, K. Ackley, S. Adhicary, N. Adhikari, R. X. Adhikari, V. K. Adkins, V. B. Adya, C. Affeldt, D. Agarwal, M. Agathos, O. D. Aguiar, L. Aiello, A. Ain, P. Ajith, T. Akutsu, S. Albanesi, R. A. Alfaidi, A. Al-Jodah, C. Alléné, A. Allocca , et al. (1719 additional authors not shown)

Abstract: The global network of gravitational-wave observatories now includes five detectors, namely LIGO Hanford, LIGO Livingston, Virgo, KAGRA, and GEO 600. These detectors collected data during their third observing run, O3, composed of three phases: O3a starting in April of 2019 and lasting six months, O3b starting in November of 2019 and lasting five months, and O3GK starting in April of 2020 and lasti… ▽ More The global network of gravitational-wave observatories now includes five detectors, namely LIGO Hanford, LIGO Livingston, Virgo, KAGRA, and GEO 600. These detectors collected data during their third observing run, O3, composed of three phases: O3a starting in April of 2019 and lasting six months, O3b starting in November of 2019 and lasting five months, and O3GK starting in April of 2020 and lasting 2 weeks. In this paper we describe these data and various other science products that can be freely accessed through the Gravitational Wave Open Science Center at https://gwosc.org. The main dataset, consisting of the gravitational-wave strain time series that contains the astrophysical signals, is released together with supporting data useful for their analysis and documentation, tutorials, as well as analysis software packages. △ Less

Submitted 7 February, 2023; originally announced February 2023.

Comments: 27 pages, 3 figures

Report number: LIGO-P2200316

Search for subsolar-mass black hole binaries in the second part of Advanced LIGO's and Advanced Virgo's third observing run

Authors: The LIGO Scientific Collaboration, the Virgo Collaboration, the KAGRA Collaboration, R. Abbott, H. Abe, F. Acernese, K. Ackley, S. Adhicary, N. Adhikari, R. X. Adhikari, V. K. Adkins, V. B. Adya, C. Affeldt, D. Agarwal, M. Agathos, O. D. Aguiar, L. Aiello, A. Ain, P. Ajith, T. Akutsu, S. Albanesi, R. A. Alfaidi, C. Alléné, A. Allocca, P. A. Altin , et al. (1680 additional authors not shown)

Abstract: We describe a search for gravitational waves from compact binaries with at least one component with mass 0.2 $M_\odot$ -- $1.0 M_\odot$ and mass ratio $q \geq 0.1$ in Advanced LIGO and Advanced Virgo data collected between 1 November 2019, 15:00 UTC and 27 March 2020, 17:00 UTC. No signals were detected. The most significant candidate has a false alarm rate of 0.2 $\mathrm{yr}^{-1}$. We estimate t… ▽ More We describe a search for gravitational waves from compact binaries with at least one component with mass 0.2 $M_\odot$ -- $1.0 M_\odot$ and mass ratio $q \geq 0.1$ in Advanced LIGO and Advanced Virgo data collected between 1 November 2019, 15:00 UTC and 27 March 2020, 17:00 UTC. No signals were detected. The most significant candidate has a false alarm rate of 0.2 $\mathrm{yr}^{-1}$. We estimate the sensitivity of our search over the entirety of Advanced LIGO's and Advanced Virgo's third observing run, and present the most stringent limits to date on the merger rate of binary black holes with at least one subsolar-mass component. We use the upper limits to constrain two fiducial scenarios that could produce subsolar-mass black holes: primordial black holes (PBH) and a model of dissipative dark matter. The PBH model uses recent prescriptions for the merger rate of PBH binaries that include a rate suppression factor to effectively account for PBH early binary disruptions. If the PBHs are monochromatically distributed, we can exclude a dark matter fraction in PBHs $f_\mathrm{PBH} \gtrsim 0.6$ (at 90% confidence) in the probed subsolar-mass range. However, if we allow for broad PBH mass distributions we are unable to rule out $f_\mathrm{PBH} = 1$. For the dissipative model, where the dark matter has chemistry that allows a small fraction to cool and collapse into black holes, we find an upper bound $f_{\mathrm{DBH}} < 10^{-5}$ on the fraction of atomic dark matter collapsed into black holes. △ Less

Submitted 26 January, 2024; v1 submitted 2 December, 2022; originally announced December 2022.

Comments: https://dcc.ligo.org/P2200139

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

Needle in a Bayes Stack: a Hierarchical Bayesian Method for Constraining the Neutron Star Equation of State with an Ensemble of Binary Neutron Star Post-merger Remnants

Authors: Alexander W. Criswell, Jesse Miller, Noah Woldemariam, Theodoros Soultanis, Andreas Bauswein, Katerina Chatziioannou, Michael W. Coughlin, Galin Jones, Vuk Mandic

Abstract: Binary neutron star (BNS) post-merger gravitational-wave emission can occur in the aftermath of a BNS merger -- provided the system avoids prompt collapse to a black hole -- as a quasistable hypermassive remnant experiences quadrupolar oscillations and non-axisymmetric deformations. The post-merger gravitational-wave spectrum possesses a characteristic peak frequency that has been shown to be depe… ▽ More Binary neutron star (BNS) post-merger gravitational-wave emission can occur in the aftermath of a BNS merger -- provided the system avoids prompt collapse to a black hole -- as a quasistable hypermassive remnant experiences quadrupolar oscillations and non-axisymmetric deformations. The post-merger gravitational-wave spectrum possesses a characteristic peak frequency that has been shown to be dependent on the binary chirp mass and the neutron star equation of state (EoS), rendering post-merger gravitational waves a powerful tool for constraining neutron star composition. Unfortunately, the BNS post-merger signal is emitted at high ($\gtrsim 1.5$ kHz) frequencies, where ground-based gravitational wave detectors suffer from reduced sensitivity. It is therefore unlikely that post-merger signals will be detected with sufficient signal-to-noise ratio (SNR) until the advent of next-generation detectors. However, by employing empirical relations derived from numerical relativity simulations, we can combine information across an ensemble of BNS mergers, allowing us to obtain EoS constraints with many low-SNR signals. We present a hierarchical Bayesian method for deriving constraints on $R_{1.6}$, the radius of a 1.6$\mathrm{M_{\odot}}$ neutron star, through an ensemble analysis of binary neutron star mergers. We apply this method to simulations of the next two LIGO-Virgo-KAGRA observing runs, O4 and O5, as well as an extended 4-year run at A+ sensitivity, demonstrating the potential of our approach to yield EoS information from the post-merger signal with current-generation detectors. The A+ 4-year scenario is predicted to improve the constraint on $R_{1.6}$ from the currently available multimessenger-based 95\% credible interval (C.I.) uncertainty of $R_{1.6}=12.07^{+0.98}_{-0.77}$ km to $R_{1.6}=11.91^{+0.80}_{-0.56}$ km, a 22% reduction of the 95% C.I. width. △ Less

Submitted 8 March, 2023; v1 submitted 9 November, 2022; originally announced November 2022.

Comments: 22 pages, 15 figures, 4 tables

Journal ref: Phys. Rev. D 107, 043021 (2023)

Search for gravitational-wave transients associated with magnetar bursts in Advanced LIGO and Advanced Virgo data from the third observing run

Authors: The LIGO Scientific Collaboration, the Virgo Collaboration, the KAGRA Collaboration, R. Abbott, H. Abe, F. Acernese, K. Ackley, N. Adhikari, R. X. Adhikari, V. K. Adkins, V. B. Adya, C. Affeldt, D. Agarwal, M. Agathos, K. Agatsuma, N. Aggarwal, O. D. Aguiar, L. Aiello, A. Ain, P. Ajith, T. Akutsu, S. Albanesi, R. A. Alfaidi, A. Allocca, P. A. Altin , et al. (1645 additional authors not shown)

Abstract: Gravitational waves are expected to be produced from neutron star oscillations associated with magnetar giant flares and short bursts. We present the results of a search for short-duration (milliseconds to seconds) and long-duration ($\sim$ 100 s) transient gravitational waves from 13 magnetar short bursts observed during Advanced LIGO, Advanced Virgo and KAGRA's third observation run. These 13 bu… ▽ More Gravitational waves are expected to be produced from neutron star oscillations associated with magnetar giant flares and short bursts. We present the results of a search for short-duration (milliseconds to seconds) and long-duration ($\sim$ 100 s) transient gravitational waves from 13 magnetar short bursts observed during Advanced LIGO, Advanced Virgo and KAGRA's third observation run. These 13 bursts come from two magnetars, SGR 1935$+$2154 and Swift J1818.0$-$1607. We also include three other electromagnetic burst events detected by Fermi GBM which were identified as likely coming from one or more magnetars, but they have no association with a known magnetar. No magnetar giant flares were detected during the analysis period. We find no evidence of gravitational waves associated with any of these 16 bursts. We place upper bounds on the root-sum-square of the integrated gravitational-wave strain that reach $2.2 \times 10^{-23}$ $/\sqrt{\text{Hz}}$ at 100 Hz for the short-duration search and $8.7 \times 10^{-23}$ $/\sqrt{\text{Hz}}$ at $450$ Hz for the long-duration search, given a detection efficiency of 50%. For a ringdown signal at 1590 Hz targeted by the short-duration search the limit is set to $1.8 \times 10^{-22}$ $/\sqrt{\text{Hz}}$. Using the estimated distance to each magnetar, we derive upper bounds on the emitted gravitational-wave energy of $3.2 \times 10^{43}$ erg ($7.3 \times 10^{43}$ erg) for SGR 1935$+$2154 and $8.2 \times 10^{42}$ erg ($2.8 \times 10^{43}$ erg) for Swift J1818.0$-$1607, for the short-duration (long-duration) search. Assuming isotropic emission of electromagnetic radiation of the burst fluences, we constrain the ratio of gravitational-wave energy to electromagnetic energy for bursts from SGR 1935$+$2154 with available fluence information. The lowest of these ratios is $3 \times 10^3$. △ Less

Submitted 19 October, 2022; originally announced October 2022.

Comments: 30 pages with appendices, 5 figures, 10 tables

Report number: LIGO-P2100387

Model-based cross-correlation search for gravitational waves from the low-mass X-ray binary Scorpius X-1 in LIGO O3 data

Authors: The LIGO Scientific Collaboration, the Virgo Collaboration, the KAGRA Collaboration, R. Abbott, H. Abe, F. Acernese, K. Ackley, S. Adhicary, N. Adhikari, R. X. Adhikari, V. K. Adkins, V. B. Adya, C. Affeldt, D. Agarwal, M. Agathos, O. D. Aguiar, L. Aiello, A. Ain, P. Ajith, T. Akutsu, S. Albanesi, R. A. Alfaidi, C. Alléné, A. Allocca, P. A. Altin , et al. (1670 additional authors not shown)

Abstract: We present the results of a model-based search for continuous gravitational waves from the low-mass X-ray binary Scorpius X-1 using LIGO detector data from the third observing run of Advanced LIGO, Advanced Virgo and KAGRA. This is a semicoherent search which uses details of the signal model to coherently combine data separated by less than a specified coherence time, which can be adjusted to bala… ▽ More We present the results of a model-based search for continuous gravitational waves from the low-mass X-ray binary Scorpius X-1 using LIGO detector data from the third observing run of Advanced LIGO, Advanced Virgo and KAGRA. This is a semicoherent search which uses details of the signal model to coherently combine data separated by less than a specified coherence time, which can be adjusted to balance sensitivity with computing cost. The search covered a range of gravitational-wave frequencies from 25Hz to 1600Hz, as well as ranges in orbital speed, frequency and phase determined from observational constraints. No significant detection candidates were found, and upper limits were set as a function of frequency. The most stringent limits, between 100Hz and 200Hz, correspond to an amplitude h0 of about 1e-25 when marginalized isotropically over the unknown inclination angle of the neutron star's rotation axis, or less than 4e-26 assuming the optimal orientation. The sensitivity of this search is now probing amplitudes predicted by models of torque balance equilibrium. For the usual conservative model assuming accretion at the surface of the neutron star, our isotropically-marginalized upper limits are close to the predicted amplitude from about 70Hz to 100Hz; the limits assuming the neutron star spin is aligned with the most likely orbital angular momentum are below the conservative torque balance predictions from 40Hz to 200Hz. Assuming a broader range of accretion models, our direct limits on gravitational-wave amplitude delve into the relevant parameter space over a wide range of frequencies, to 500Hz or more. △ Less

Submitted 2 January, 2023; v1 submitted 6 September, 2022; originally announced September 2022.

Comments: 19 pages, Open Access Journal PDF

Report number: LIGO-P2100110-v13

Journal ref: The Astrophysical Journal Letters, 941, L30 (2022)

GROWTH on S190426c II: GROWTH-India Telescope search for an optical counterpart with a custom image reduction and candidate vetting pipeline

Authors: Harsh Kumar, Varun Bhalerao, G. C. Anupama, Sudhanshu Barway, Michael W. Coughlin, Kishalay De, Kunal Deshmukh, Anirban Dutta, Daniel A Goldstein, Adeem Jassani, Simran Joharle, Viraj Karambelker, Maitreya Khandagale, Brajesh Kumar, Divita Saraogi, Yashvi Sharma, Vedant Shenoy, Leo singer, Avinash Singh, Gaurav Waratkar

Abstract: S190426c / GW190426_152155 was the first probable neutron star - black hole merger candidate detected by the LIGO-Virgo Collaboration. We undertook a tiled search for optical counterparts of this event using the 0.7m GROWTH-India Telescope. Over a period of two weeks, we obtained multiple observations over a 22.1 deg^2 area, with a 17.5% probability of containing the source location. Initial effor… ▽ More S190426c / GW190426_152155 was the first probable neutron star - black hole merger candidate detected by the LIGO-Virgo Collaboration. We undertook a tiled search for optical counterparts of this event using the 0.7m GROWTH-India Telescope. Over a period of two weeks, we obtained multiple observations over a 22.1 deg^2 area, with a 17.5% probability of containing the source location. Initial efforts included obtaining photometry of sources reported by various groups, and a visual search for sources in all galaxies contained in the region. Subsequently, we have developed an image subtraction and candidate vetting pipeline with ~ 94% efficiency for transient detection. Processing the data with this pipeline, we find several transients, but none that are compatible with kilonova models. We present the details of our observations, working of our pipeline, results from the search, and our interpretations of the non-detections that will work as a pathfinder during the O4 run of LVK. △ Less

Submitted 5 September, 2022; originally announced September 2022.

Comments: 12 pages, 6 figures, Accepted for publication in Monthly Notices of the Royal Astronomical Society

Correlated 1-1000 Hz magnetic field fluctuations from lightning over earth-scale distances and their impact on gravitational wave searches

Authors: Kamiel Janssens, Matthew Ball, Robert M. S. Schofield, Nelson Christensen, Raymond Frey, Nick van Remortel, Sharan Banagiri, Michael W. Coughlin, Anamaria Effler, Mark Gołkowski, Jerzy Kubisz, Michał Ostrowski

Abstract: We report Earth-scale distance magnetic correlations from lightning strokes in the frequency range 1-1000 Hz at several distances ranging from 1100 to 9000 km. Noise sources which are correlated on Earth-scale distances can affect future searches for gravitational-wave signals with ground-based gravitational-wave interferometric detectors. We consider the impact of correlations from magnetic field… ▽ More We report Earth-scale distance magnetic correlations from lightning strokes in the frequency range 1-1000 Hz at several distances ranging from 1100 to 9000 km. Noise sources which are correlated on Earth-scale distances can affect future searches for gravitational-wave signals with ground-based gravitational-wave interferometric detectors. We consider the impact of correlations from magnetic field fluctuations on gravitational-wave searches due to Schumann resonances ($<$50 Hz) as well as higher frequencies ($>$100 Hz). We demonstrate that individual lightning strokes are a likely source for the observed correlations in the magnetic field fluctuations at gravitational-wave observatories and discuss some of their characteristics. Furthermore, we predict their impact on searches for an isotropic gravitational-wave background, as well as for searches looking for short-duration transient gravitational waves, both unmodeled signals (bursts) as well as modeled signals (compact binary coalescence). Whereas the recent third observing run by LIGO and Virgo was free of an impact from correlated magnetic field fluctuations, future runs could be affected. For example, at current magnetic coupling levels, neutron star inspirals in third generation detectors are likely to be contaminated by multiple correlated lightning glitches. We suggest that future detector design should consider reducing lightning coupling by, for example, reducing the lightning-induced beam tube currents that pass through sensitive magnetic coupling regions in current detectors. We also suggest that the diurnal and seasonal variation in lightning activity may be useful in discriminating between detector correlations that are produced by gravitational waves and those produced by lightning. △ Less

Submitted 1 September, 2022; originally announced September 2022.