Search

We present a systematic search for optical counterparts to 13 gravitational wave (GW) triggers involving at least one neutron star during LIGO/Virgo's third observing run (O3). We searched binary neutron star (BNS) and neutron star black hole (NSBH) merger localizations with the Zwicky Transient Facility (ZTF) and undertook follow-up with the Global Relay of Observatories Watching Transients Happen (GROWTH) collaboration. The GW triggers had a median localization area of 4480 deg2, a median distance of 267 Mpc, and false-alarm rates ranging from 1.5 to 10-25 yr-1. The ZTF coverage in the g and r bands had a median enclosed probability of 39%, median depth of 20.8 mag, and median time lag between merger and the start of observations of 1.5 hr. The O3 follow-up by the GROWTH team comprised 340 UltraViolet/Optical/InfraRed (UVOIR) photometric points, 64 OIR spectra, and three radio images using 17 different telescopes. We find no promising kilonovae (radioactivity-powered counterparts), and we show how to convert the upper limits to constrain the underlying kilonova luminosity function. Initially, we assume that all GW triggers are bona fide astrophysical events regardless of false-alarm rate and that kilonovae accompanying BNS and NSBH mergers are drawn from a common population; later, we relax these assumptions. Assuming that all kilonovae are at least as luminous as the discovery magnitude of GW170817 (-16.1 mag), we calculate that our joint probability of detecting zero kilonovae is only 4.2%. If we assume that all kilonovae are brighter than-16.6 mag (the extrapolated peak magnitude of GW170817) and fade at a rate of 1 mag day-1 (similar to GW170817), the joint probability of zero detections is 7%. If we separate the NSBH and BNS populations based on the online classifications, the joint probability of zero detections, assuming all kilonovae are brighter than-16.6 mag, is 9.7% for NSBH and 7.9% for BNS mergers. Moreover, no more than 10-4, or φ > 30° to be consistent with our limits. We look forward to searches in the fourth GW observing run; even 17 neutron star mergers with only 50% coverage to a depth of-16 mag would constrain the maximum fraction of bright kilonovae to

The Advanced X-ray Imaging Satellite (AXIS) promises revolutionary science in the X-ray and multi-messenger time domain. AXIS will leverage excellent spatial resolution (<1.5 arcsec), sensitivity ( 80 × that of Swift), and a large collecting area (5–10× that of Chandra) across a 24-arcmin diameter field of view at soft X-ray energies (0.3–10.0 keV) to discover and characterize a wide range of X-ray transients from supernova-shock breakouts to tidal disruption events to highly variable supermassive black holes. The observatory's ability to localize and monitor faint X-ray sources opens up new opportunities to hunt for counterparts to distant binary neutron star mergers, fast radio bursts, and exotic phenomena like fast X-ray transients. AXIS will offer a response time of <2 h to community alerts, enabling studies of gravitational wave sources, high-energy neutrino emitters, X-ray binaries, magnetars, and other targets of opportunity. This white paper highlights some of the discovery science that will be driven by AXIS in this burgeoning field of time domain and multi-messenger astrophysics. This White Paper is part of a series commissioned for the AXIS Probe Concept Mission; additional AXIS White Papers can be found at the AXIS website. [ABSTRACT FROM AUTHOR]

On 2017 August 17, Advanced LIGO and Virgo observed GW170817, the first gravitational-wave (GW) signal from a binary neutron star merger. It was followed by a short-duration gamma-ray burst, GRB 170817A, and by a non-thermal afterglow emission. In this work, a combined simultaneous fit of the electromagnetic (EM, specifically, afterglow) and GW domains is implemented, both using the posterior distribution of a GW standalone analysis as prior distribution to separately process the EM data, and fitting the EM and GW domains simultaneously. These approaches coincide mathematically, as long as the actual posterior of the GW analysis, and not an approximation, is used as prior for the EM analysis. We treat the viewing angle, θ v , as shared parameter across the two domains. In the afterglow modelling with a Gaussian structured jet this parameter and the jet core angle, θc, are correlated, leading to high uncertainties on their values. The joint EM + GW analysis relaxes this degeneracy, reducing the uncertainty compared to an EM-only fit. We also apply our methodology to hypothetical GW170817-like events occurring in the next GW observing run at ∼140 and 70 Mpc. At 70 Mpc the existing EM degeneracy is broken, thanks to the inclusion of the GW domain in the analysis. At 140 Mpc, the EM-only fit cannot constrain θ v nor θc because of the lack of detections in the afterglow rising phase. Folding the GW data into the analysis leads to tighter constraints on θ v , still leaving θc unconstrained, requiring instruments with higher sensitivities, such as Athena. [ABSTRACT FROM AUTHOR]

Swift has now completed 18 years of mission, during which it discovered thousands of gamma-ray bursts as well as new classes of high-energy transient phenomena. Its first breakthrough result was the localization of short duration GRBs, which enabled for redshift measurements and kilonova searches. Swift, in synergy with the Hubble Space Telescope and a wide array of ground-based telescopes, provided the first tantalizing evidence of a kilonova in the aftermath of a short GRB. In 2017, Swift observations of the gravitational wave event GW170817 captured the early UV photons from the kilonova AT2017gfo, opening a new window into the physics of kilonovae. Since then, Swift has continued to expand the sample of known kilonovae, leading to the surprising discovery of a kilonova in a long duration GRB. This article will discuss recent advances in the study of kilonovae driven by the fundamental contribution of Swift. [ABSTRACT FROM AUTHOR]

In this paper we explore the scientific synergies between Athena and some of the key multi-messenger facilities that should be operative concurrently with Athena. These facilities include LIGO A+, Advanced Virgo+ and future detectors for ground-based observation of gravitational waves (GW), LISA for space-based observations of GW, IceCube and KM3NeT for neutrino observations, and CTA for very high energy observations. These science themes encompass pressing issues in astrophysics, cosmology and fundamental physics such as: the central engine and jet physics in compact binary mergers, accretion processes and jet physics in Super-Massive Binary Black Holes (SMBBHs) and in compact stellar binaries, the equation of state of neutron stars, cosmic accelerators and the origin of Cosmic Rays (CRs), the origin of intermediate and high-Z elements in the Universe, the Cosmic distance scale and tests of General Relativity and the Standard Model. Observational strategies for implementing the identified science topics are also discussed. A significant part of the sources targeted by multi-messenger facilities is of transient nature. We have thus also discussed the synergy of Athena with wide-field high-energy facilities, taking THESEUS as a case study for transient discovery. This discussion covers all the Athena science goals that rely on follow-up observations of high-energy transients identified by external observatories, and includes also topics that are not based on multi-messenger observations, such as the search for missing baryons or the observation of early star populations and metal enrichment at the cosmic dawn with Gamma-Ray Bursts (GRBs). [ABSTRACT FROM AUTHOR]

Gamma-ray bursts (GRBs) associated with gravitational wave events are, and will likely continue to be, viewed at a larger inclination than GRBs without gravitational wave detections. As demonstrated by the afterglow of GW 170817A, this requires an extension of the common GRB afterglow models, which typically assume emission from an on-axis top-hat jet. We present a characterization of the afterglows arising from structured jets, providing a framework covering both successful and choked jets. We compute new closure relations for decelerating structured jets and compare them with the established relations for energy injection and refreshed shock models. The temporal slope before the jet break is found to be a simple function of the ratio between the viewing angle and effective opening angle of the jet. A numerical model to calculate synthetic light curves and spectra is publicly available as the open-source Python package afterglowpy. [ABSTRACT FROM AUTHOR]