The detection of the optical transient AT2017gfo proved that binary neutron star mergers are progenitors of kilonovae (KNe). Using a combination of numerical-relativity and radiative-transfer simulations, the community has developed sophisticated models for these transients for a wide portion of the expected parameter space. Using these simulations and surrogate models made from them, it has been possible to perform Bayesian inference of the observed signals to infer properties of the ejected matter. It has been pointed out that combining inclination constraints derived from the KN with gravitational-wave measurements increases the accuracy with which binary parameters can be estimated, in particular breaking the distance-inclination degeneracy from gravitational wave inference. To avoid bias from the unknown ejecta geometry, constraints on the inclination angle for AT2017gfo should be insensitive to the employed models. In this work, we compare different assumptions about the ejecta and radiative reprocesses used by the community and we investigate their impact on the parameter inference. While most inferred parameters agree, we find disagreement between posteriors for the inclination angle for different geometries that have been used in the current literature. According to our study, the inclusion of reprocessing of the photons between different ejecta types improves the modeling fits to AT2017gfo and, in some cases, affects the inferred constraints. Our study motivates the inclusion of large similar to 1-mag uncertainties in the KN models employed for Bayesian analysis to capture yet unknown systematics, especially when inferring inclination angles, although smaller uncertainties seem appropriate to capture model systematics for other intrinsic parameters. We can use this method to impose soft constraints on the ejecta geometry of the KN AT2017gfo.

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

The most common way to discover extragalactic fast transients, which fade within a few nights in the optical, is via follow-up of gamma-ray burst and gravitational-wave triggers. However, wide-field surveys have the potential to identify rapidly fading transients independently of such external triggers. The volumetric survey speed of the Zwicky Transient Facility (ZTF) makes it sensitive to objects as faint and fast fading as kilonovae, the optical counterparts to binary neutron star mergers, out to almost 200 Mpc. We introduce an open-source software infrastructure, the ZTF REaltime Search and Triggering, ZTFReST, designed to identify kilonovae and fast transients in ZTF data. Using the ZTF alert stream combined with forced point-spread-function photometry, we have implemented automated candidate ranking based on their photometric evolution and fitting to kilonova models. Automated triggering, with a human in the loop for monitoring, of follow-up systems has also been implemented. In 13 months of science validation, we found several extragalactic fast transients independently of any external trigger, including two supernovae with post-shock cooling emission, two known afterglows with an associated gamma-ray burst (ZTF20abbiixp, ZTF20abwysqy), two known afterglows without any known gamma-ray counterpart (ZTF20aajnksq, ZTF21aaeyldq), and three new fast-declining sources (ZTF20abtxwfx, ZTF20acozryr, ZTF21aagwbjr) that are likely associated with GRB200817A, GRB201103B, and GRB210204A. However, we have not found any objects that appear to be kilonovae. We constrain the rate of GW170817-like kilonovae to R < 900 Gpc(-3) yr(-1) (95% confidence). A framework such as ZTFReST could become a prime tool for kilonova and fast-transient discovery with the Vera Rubin Observatory.

We present the goals, strategy, and first results of the high-cadence Galactic plane survey using the Zwicky Transient Facility (ZTF). The goal of the survey is to unveil the Galactic population of short-period variable stars, including short-period binaries, and stellar pulsators with periods less than a few hours. Between 2018 June and 2019 January, we observed 64 ZTF fields resulting in 2990 deg² of high stellar density in the ZTF-r band along the Galactic plane. Each field was observed continuously for 1.5 to 6 h with a cadence of 40 sec. Most fields have between 200 and 400 observations obtained over 2–3  continuous nights. As part of this survey, we extract a total of ≈230 million individual objects with at least 80 epochs obtained during the high-cadence Galactic plane survey reaching an average depth of ZTF–r ≈ 20.5 mag. For four selected fields with 2–10 million individual objects per field, we calculate different variability statistics and find that ≈1–2  per cent of the objects are astrophysically variable over the observed period. We present a progress report on recent discoveries, including a new class of compact pulsators, the first members of a new class of Roche lobe filling hot subdwarf binaries as well as new ultracompact double white dwarfs and flaring stars. Finally, we present a sample of 12 new single-mode hot subdwarf B-star pulsators with pulsation amplitudes between ZTF–r = 20–76 mmag and pulsation periods between P = 5.8–16 min with a strong cluster of systems with periods ≈6 min. All of the data have now been released in either ZTF Data Release 3 or Data Release 4.

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

Joint multimessenger observations with gravitational waves and electromagnetic (EM) data offer new insights into the astrophysical studies of compact objects. The third Advanced LIGO and Advanced Virgo observing run began on 2019 April 1; during the 11 months of observation, there have been 14 compact binary systems candidates for which at least one component is potentially a neutron star. Although intensive follow-up campaigns involving tens of ground and space-based observatories searched for counterparts, no EM counterpart has been detected. Following on a previous study of the first six months of the campaign, we present in this paper the next five months of the campaign from 2019 October to 2020 March. We highlight two neutron star-black hole candidates (S191205ah and S200105ae), two binary neutron star candidates (S191213g and S200213t), and a binary merger with a possible neutron star and a `MassGap' component, S200115j. Assuming that the gravitational-wave (GW) candidates are of astrophysical origin and their location was covered by optical telescopes, we derive possible constraints on the matter ejected during the events based on the non-detection of counterparts. We find that the follow-up observations during the second half of the third observing run did not meet the necessary sensitivity to constrain the source properties of the potential GW candidate. Consequently, we suggest that different strategies have to be used to allow a better usage of the available telescope time. We examine different choices for follow-up surveys to optimize sky localization coverage versus observational depth to understand the likelihood of counterpart detection.

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 deg(2), 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 <57% (<89%) of putative kilonovae could be brighter than -16.6 mag assuming flat evolution (fading by 1 mag day(-1)) at the 90% confidence level. If we further take into account the online terrestrial probability for each GW trigger, we find that no more than <68% of putative kilonovae could be brighter than -16.6 mag. Comparing to model grids, we find that some kilonovae must have M-ej M, X-lan > 10(-4), or > 30 degrees 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 <25%.

The detection of GW170817 is revolutionizing many areas of astrophysics with the joint observation of gravitational waves and electromagnetic emissions. These multimessenger events provide a new approach to determine the Hubble constant, thus, they are a promising candidate for mitigating the tension between measurements of type-Ia supernovae via the local distance ladder and the cosmic microwave background. In addition to the “standard siren” provided by the gravitational-wave measurement, the kilonova itself has characteristics that allow one to improve existing measurements or to perform yet another, independent measurement of the Hubble constant without gravitational-wave information. Here, we employ standardization techniques borrowed from the type-Ia community and apply them to kilonovae, not using any information from the gravitational-wave signal. We use two versions of this technique, one derived from direct observables measured from the light curve, and the other based on inferred ejecta parameters, e.g., mass, velocity, and composition, for two different models. These lead to Hubble constant measurements of H0=109+49−35 km s−1 Mpc−1 for the measured analysis, and H0=85+22−17 km s−1 Mpc−1 and H0=79+23−15 km s−1 Mpc−1 for the inferred analyses. This measurement has error bars within ∼2 to the gravitational-wave measurements (H0=74+16−8 km s−1 Mpc−1), showing its promise as an independent constraint on H0.

We describe a dynamic science portal called the GROWTH Marshal that allows time-domain astronomers to define science programs; program filters to save sources from different discovery streams; coordinate follow-up with various robotic or classical telescopes; analyze the panchromatic follow-up data; and generate summary tables for publication. The GROWTH marshal currently serves 137 scientists, 38 science programs, and 67 telescopes. Every night, in real time, several science programs apply various customized filters to the 105 nightly alerts from the Zwicky Transient Facility. Here, we describe the schematic and explain the functionality of the various components of this international collaborative platform.