No Arabic abstract
We present a systematic search for optical counterparts to 13 gravitational wave (GW) triggers involving at least one neutron star during LIGO/Virgos third observing run. We searched binary neutron star (BNS) and neutron star black hole (NSBH) merger localizations with the Zwicky Transient Facility (ZTF) and undertook follow-up with the Global Relay of Observatories Watching Transients Happen (GROWTH) collaboration. The GW triggers had a median localization of 4480 deg^2, median distance of 267 Mpc and false alarm rates ranging from 1.5 to 1e-25 per yr. The ZTF coverage had a median enclosed probability of 39%, median depth of 20.8mag, and median response time of 1.5 hr. The O3 follow-up by the GROWTH team comprised 340 UVOIR photometric points, 64 OIR spectra, and 3 radio. We find no promising kilonova (radioactivity-powered counterpart) and we convert the upper limits to constrain the underlying kilonova luminosity function. Assuming that all kilonovae are at least as luminous as GW170817 at discovery (-16.1mag), we calculate our joint probability of detecting zero kilonovae is only 4.2%. If we assume that all kilonovae are brighter than -16.6mag (extrapolated peak magnitude of GW170817) and fade at 1 mag/day (similar to GW170817), the joint probability of zero detections is 7%. If we separate the NSBH and BNS populations, the joint probability of zero detections, assuming all kilonovae are brighter than -16.6mag, is 9.7% for NSBH and 7.9% for BNS mergers. Moreover, <57% (<89%) of putative kilonovae could be brighter than -16.6mag assuming flat (fading) evolution, at 90% confidence. If we further account for the online terrestrial probability for each GW trigger, we find that <68% of putative kilonovae could be brighter than -16.6mag. Comparing to model grids, we find that some kilonovae must have Mej < 0.03 Msun or Xlan>1e-4 or phi>30deg to be consistent with our limits. (Abridged)
The discovery rate of fast radio bursts (FRBs) is increasing dramatically thanks to new radio facilities. Meanwhile, wide-field instruments such as the 47 deg$^2$ Zwicky Transient Facility (ZTF) survey the optical sky to study transient and variable sources. We present serendipitous ZTF observations of the CHIME repeating source FRB 180916.J0158+65, that was localized to a spiral galaxy 149 Mpc away and is the first FRB suggesting periodic modulation in its activity. While 147 ZTF exposures corresponded to expected high-activity periods of this FRB, no single ZTF exposure was at the same time as a CHIME detection. No $>3sigma$ optical source was found at the FRB location in 683 ZTF exposures, totalling 5.69 hours of integration time. We combined ZTF upper limits and expected repetitions from FRB 180916.J0158+65 in a statistical framework using a Weibull distribution, agnostic of periodic modulation priors. The analysis yielded a constraint on the ratio between the optical and radio fluences of $eta lesssim 200$, corresponding to an optical energy $E_{rm opt} lesssim 3 times 10^{46}$ erg for a fiducial 10 Jy ms FRB (90% confidence). A deeper (but less statistically robust) constraint of $eta lesssim 3$ can be placed assuming a rate of $r(>5$ Jy ms)= hr$^{-1}$ and $1.2pm 1.1$ FRB occurring during exposures taken in high-activity windows. The constraint can be improved with shorter per-image exposures and longer integration time, or observing FRBs at higher Galactic latitudes. This work demonstrated how current surveys can statistically constrain multi-wavelength counterparts to FRBs even without deliberately scheduled simultaneous radio observation.
The first binary neutron star merger, GW170817, was accompanied by a radioactivity-powered optical/infrared transient called a kilonova. To date, no compelling kilonova has been found during optical surveys of the sky, independent of gravitational-wave triggers. In this work, we searched the first 23 months of the Zwicky Transient Facility (ZTF) data stream for candidate kilonovae in the form of rapidly evolving transients. We combined ZTF alert queries with forced point-spread-function photometry and nightly flux stacking to increase our sensitivity to faint and fast transients. Automatic queries yielded $>11,200$ candidates, 24 of which passed quality checks and strict selection criteria based on a grid of kilonova models tailored for both binary neutron star and neutron star-black hole mergers. None of the candidates in our sample was deemed a possible kilonova after thorough vetting, catalog cross-matching, and study of their color evolution. The sources that passed our selection criteria are dominated by Galactic cataclysmic variables. In addition, we identified two fast transients at high Galactic latitude, one of which is the confirmed afterglow of long-duration GRB190106A, and the other is a possible cosmological afterglow. Using a survey simulation code, we constrained the kilonova rate for a range of models including top-hat and linearly decaying light curves and synthetic light curves obtained with radiative transfer simulations. For prototypical GW170817-like kilonovae, we constrain the rate to be $R < 1775$ Gpc$^{-3}$ yr$^{-1}$ at 95% confidence level by requiring at least 2 high-significance detections. By assuming a population of kilonovae with the same geometry and composition of GW170817 observed under a uniform viewing angle distribution, we obtained a constraint on the rate of $R < 4029$ Gpc$^{-3}$ yr$^{-1}$.
We present detailed simulations of black hole-neutron star (BH-NS) mergers kilonova and gamma-ray burst (GRB) afterglow and kilonova luminosity function, and discuss the detectability of electromagnetic (EM) counterpart in connection with gravitational wave (GW) detections, GW-triggered target-of-opportunity observations, and time-domain blind searches. The predicted absolute magnitude of the BH-NS kilonovae at $0.5,{rm days}$ after the merger falls in $[-10,-15.5]$. The simulated luminosity function contains the potential viewing-angle distribution information of the anisotropic kilonova emission. We simulate the GW detection rates, detectable distances and signal duration, for the future networks of 2nd/2.5th/3rd-generation GW detectors. BH-NSs tend to produce brighter kilonovae and afterglows if the BH has a higher aligned-spin, and a less massive NS with a stiffer EoS. The detectability of kilonova is especially sensitive to the BH spin. If BHs typically have low spins, the BH-NS EM counterparts are hard to discover. For the 2nd generation GW detector networks, a limiting magnitude of $m_{rm limit}sim23-24,{rm mag}$ is required to detect the kilonovae even if BH high spin is assumed. Thus, a plausible explanation for the lack of BH-NS associated kilonova detection during LIGO/Virgo O3 is that either there is no EM counterpart (plunging events), or the current follow-ups are too shallow. These observations still have the chance to detect the on-axis jet afterglow associated with an sGRB or an orphan afterglow. Follow-up observations can detect possible associated sGRB afterglows, from which kilonova signatures may be studied. For time-domain observations, a high-cadence search in redder filters is recommended to detect more BH-NS associated kilonovae and afterglows.
We present radiative transfer simulations for blue kilonovae hours after neutron star (NS) mergers by performing detailed opacity calculations for the first time. We calculate atomic structures and opacities of highly ionized elements (up to the tenth ionization) with atomic number Z = 20 - 56. We find that the bound-bound transitions of heavy elements are the dominant source of the opacities in the early phase (t < 1 day after the merger), and that the ions with a half-closed electron shell provide the highest contributions. The Planck mean opacity for lanthanide-free ejecta (with electron fraction of Ye = 0.30 - 0.40) can only reach around kappa ~ 0.5 - 1 cm^2 g^-1 at t = 0.1 day, whereas that increases up to kappa ~ 5 - 10 cm^2 g^-1 at t = 1 day. The spherical ejecta model with an ejecta mass of Mej = 0.05Msun gives the bolometric luminosity of ~ 2 x 10^42 erg s^-1 at t ~ 0.1 day. We confirm that the existing bolometric and multi-color data of GW170817 can be naturally explained by the purely radioactive model. The expected early UV signals reach 20.5 mag at t ~ 4.3 hours for sources even at 200 Mpc, which is detectable by the facilities such as Swift and the Ultraviolet Transient Astronomy Satellite (ULTRASAT). The early-phase luminosity is sensitive to the structure of the outer ejecta, as also pointed out by Kasen et al. (2017). Therefore, the early UV observations give strong constraints on the structure of the outer ejecta as well as the presence of a heating source besides r-process nuclei.
The Zwicky Transient Facility (ZTF) is conducting a nightly public survey of all 13 TESS northern sectors in 2019-2020. ZTF will observe the portions of the current TESS sectors visible from Palomar Observatory each night. Each ZTF pointing will have one exposure each with $g$ and $r$ filters, totaling two images per night. ZTF is also making additional nightly $g$- and $r$-band observations of denser stellar regions (e.g. near the Galactic Plane) to better facilitate variability studies of Galactic objects. The limiting magnitude of the ZTF observations is $approx$20.6 and ZTF saturates at magnitude $approx$13. ZTF will release data from TESS fields in three forms: nightly alerts distributed by established ZTF brokers, nightly alerts converted to JSON format are distributed via ZTFs bucket on Google Cloud as a tarball, and monthly photometric light curves also distributed via Google Cloud.