Search

A supermassive black hole (SMBH) ejected from the potential well of its host galaxy via gravitational wave recoil carries important information about the mass ratio and spin alignment of the pre-merger SMBH binary. Such a recoiling SMBH may be detectable as an active galactic nucleus (AGN) broad-line region offset by up to 10 kpc from a disturbed host galaxy. We describe a novel methodology using forward modeling with The Tractor to search for such offset AGNs in a sample of 5493 optically variable AGNs detected with the Zwicky Transient Facility (ZTF). We present the discovery of nine AGNs that may be spatially offset from their host galaxies and are candidates for recoiling SMBHs. Five of these offset AGNs exhibit double-peaked broad Balmer lines, which may have arisen from unobscured accretion disk emission, and four show radio emission indicative of a relativistic jet. The fraction of double-peaked emitters in our spatially offset AGN sample is significantly larger than the 16% double-peaked emitter fraction observed for ZTF AGNs overall. In our sample of variable AGNs we also identified 52 merging galaxies, including a new spectroscopically confirmed dual AGN. Finally, we detected the dramatic rebrightening of SDSS 1133, a previously discovered variable object and recoiling SMBH candidate, in ZTF. The flare was accompanied by the reemergence of strong P Cygni line features, indicating that SDSS 1133 may be an outbursting luminous blue variable star.

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

Colors of Type Ia supernovae (SNe Ia) in the first few days after explosion provide a potential discriminant between different models. In this paper, we present g-r colors of 65 SNe Ia discovered within 5 days from first light by the Zwicky Transient Facility in 2018, a sample that is about three times larger than that in the literature. We find that g-r colors are intrinsically rather homogeneous at early phases, with about half of the dispersion attributable to photometric uncertainties (σnoise ∼ σ int ∼ 0.18 mag). Colors are nearly constant starting from 6 days after first light (g-r ∼-0.15 mag), while the time evolution at earlier epochs is characterized by a continuous range of slopes, from events rapidly transitioning from redder to bluer colors (slope of ∼-0.25 mag day-1) to events with a flatter evolution. The continuum in the slope distribution is in good agreement both with models requiring some amount of 56Ni mixed in the outermost regions of the ejecta and with "double-detonation"models having thin helium layers MHe=0.01 M⊙) and varying carbon-oxygen core masses. At the same time, six events show evidence for a distinctive "red bump"signature predicted by double-detonation models with larger helium masses. We finally identify a significant correlation between the early-time g-r slopes and supernova brightness, with brighter events associated to flatter color evolution (p-value = 0.006). The distribution of slopes, however, is consistent with being drawn from a single population, with no evidence for two components as claimed in the literature based on B-V colors.

While it is clear that Type Ia supernovae (SNe) are the result of thermonuclear explosions in C/O white dwarfs (WDs), a great deal remains uncertain about the binary companion that facilitates the explosive disruption of the WD. Here, we present a comprehensive analysis of a large, unique data set of 127 SNe Ia with exquisite coverage by the Zwicky Transient Facility (ZTF). High-cadence (six observations per night) ZTF observations allow us to measure the SN rise time and examine its initial evolution. We develop a Bayesian framework to model the early rise as a power law in time, which enables the inclusion of priors in our model. For a volume-limited subset of normal SNe Ia, we find that the mean power-law index is consistent with 2 in the r ZTF-band (α r=2.01± 0.02), as expected in the expanding fireball model. There are, however, individual SNe that are clearly inconsistent with α r=2. We estimate a mean rise time of 18.9 days (with a range extending from ∼15 to 22 days), though this is subject to the adopted prior. We identify an important, previously unknown, bias whereby the rise times for higher-redshift SNe within a flux-limited survey are systematically underestimated. This effect can be partially alleviated if the power-law index is fixed to α = 2, in which case we estimate a mean rise time of 21.7 days (with a range from ∼18 to 23 days). The sample includes a handful of rare and peculiar SNe Ia. Finally, we conclude with a discussion of lessons learned from the ZTF sample that can eventually be applied to observations from the Vera C. Rubin Observatory.

High-cadence transient surveys are able to capture supernovae closer to their first light than ever before. Applying analytical models to such early emission, we can constrain the progenitor stars’ properties. In this paper, we present observations of SN 2018fif (ZTF 18abokyfk). The supernova was discovered close to first light and monitored by the Zwicky Transient Facility (ZTF) and the Neil Gehrels Swift Observatory. Early spectroscopic observations suggest that the progenitor of SN 2018fif was surrounded by relatively small amounts of circumstellar material compared to all previous cases. This particularity, coupled with the high-cadence multiple-band coverage, makes it a good candidate to investigate using shock-cooling models. We employ the SOPRANOS code, an implementation of the model by Sapir & Waxman and its extension to early times by Morag et al. Compared with previous implementations, SOPRANOS has the advantage of including a careful account of the limited temporal validity domain of the shock-cooling model as well as allowing usage of the entirety of the early UV data. We find that the progenitor of SN 2018fif was a large red supergiant with a radius of R = 744.0-+128.0183.0 R☉ and an ejected mass of Mej = 9.3-+5.80.4 M☉. Our model also gives information on the explosion epoch, the progenitor’s inner structure, the shock velocity, and the extinction. The distribution of radii is double-peaked, with smaller radii corresponding to lower values of the extinction, earlier recombination times, and a better match to the early UV data. If these correlations persist in future objects, denser spectroscopic monitoring constraining the time of recombination, as well as accurate UV observations (e.g., with ULTRASAT), will help break the extinction/radius degeneracy and independently determine both.

© 2020. The American Astronomical Society. All rights reserved.. We present a survey of the early evolution of 12 Type IIn supernovae (SNe IIn) at ultraviolet and visible light wavelengths. We use this survey to constrain the geometry of the circumstellar material (CSM) surrounding SN IIn explosions, which may shed light on their progenitor diversity. In order to distinguish between aspherical and spherical CSM, we estimate the blackbody radius temporal evolution of the SNe IIn of our sample, following the method introduced by Soumagnac et al. We find that higher-luminosity objects tend to show evidence for aspherical CSM. Depending on whether this correlation is due to physical reasons or to some selection bias, we derive a lower limit between 35% and 66% for the fraction of SNe IIn showing evidence for aspherical CSM. This result suggests that asphericity of the CSM surrounding SNe IIn is common - consistent with data from resolved images of stars undergoing considerable mass loss. It should be taken into account for more realistic modeling of these events.

Early observations of Type Ia supernovae (SNe Ia) provide essential clues for understanding the progenitor system that gave rise to the terminal thermonuclear explosion. We present exquisite observations of SN 2019yvq, the second observed SN Ia, after iPTF 14atg, to display an early flash of emission in the ultraviolet (UV) and optical. Our analysis finds that SN 2019yvq was unusual, even when ignoring the initial flash, in that it was moderately underluminous for an SN Ia (Mg ≈ 18.5 mag at peak) yet featured very high absorption velocities (v ≈ 15,000 km s-1 for Si ii λ6355 at peak). We find that many of the observational features of SN 2019yvq, aside from the flash, can be explained if the explosive yield of radioactive 56Ni is relatively low (we measure M56Ni=0.31±0.05,M⊙) and it and other iron-group elements are concentrated in the innermost layers of the ejecta. To explain both the UV/optical flash and peak properties of SN 2019yvq we consider four different models: interaction between the SN ejecta and a nondegenerate companion, extended clumps of 56Ni in the outer ejecta, a double-detonation explosion, and the violent merger of two white dwarfs. Each of these models has shortcomings when compared to the observations; it is clear additional tuning is required to better match SN 2019yvq. In closing, we predict that the nebular spectra of SN 2019yvq will feature either H or He emission, if the ejecta collided with a companion, strong [Ca ii] emission, if it was a double detonation, or narrow [O i] emission, if it was due to a violent merger.

We present detailed observations of ZTF18abukavn (SN2018gep), discovered in high-cadence data from the Zwicky Transient Facility as a rapidly rising (1.4 ± 0.1 mag hr-1) and luminous (Mg,peak = -20 mag) transient. It is spectroscopically classified as a broad-lined stripped-envelope supernova (Ic-BL SN). The high peak luminosity (Lbol ≳ 3 × 1044 erg s-1), the short rise time (trise = 3 days in g band), and the blue colors at peak (g-r ∼ -0.4) all resemble the high-redshift Ic-BL iPTF16asu, as well as several other unclassified fast transients. The early discovery of SN2018gep (within an hour of shock breakout) enabled an intensive spectroscopic campaign, including the highest-temperature (Teff ≳ 40,000 K) spectra of a stripped-envelope SN. A retrospective search revealed luminous (Mg ∼ Mr ≈ mag) emission in the days to weeks before explosion, the first definitive detection of precursor emission for a Ic-BL. We find a limit on the isotropic gamma-ray energy release E γ,iso < 4.9 × 10 48 erg, a limit on X-ray emission LX < 1040 erg s-1, and a limit on radio emission ν Lν ≲ 1037 erg s-1. Taken together, we find that the early (< 10 days) data are best explained by shock breakout in a massive shell of dense circumstellar material (0.02 M⊙) at large radii (3 × 1014 cm) that was ejected in eruptive pre-explosion mass-loss episodes. The late-time (> 10 days) light curve requires an additional energy source, which could be the radioactive decay of Ni-56.

Early-time observations of Type Ia supernovae (SNe Ia) are essential to constrain the properties of their progenitors. In this paper, we present high-quality light curves of 127 SNe Ia discovered by the Zwicky Transient Facility (ZTF) in 2018. We describe our method to perform forced point-spread function photometry, which can be applied to other types of extragalactic transients. With a planned cadence of six observations per night (three g + three r), all of the 127 SNe Ia are detected in both g and r bands more than 10 days (in the rest frame) prior to the epoch of g-band maximum light. The redshifts of these objects range from z = 0.0181 to 0.165; the median redshift is 0.074. Among the 127 SNe, 50 are detected at least 14 days prior to maximum light (in the rest frame), with a subset of nine objects being detected more than 17 days before g-band peak. This is the largest sample of young SNe Ia collected to date; it can be used to study the shape and color evolution of the rising light curves in unprecedented detail. We discuss six peculiar events in this sample: one 02cx-like event ZTF18abclfee (SN 2018crl), one Ia-CSM SN ZTF18aaykjei (SN 2018cxk), and four objects with possible super-Chandrasekhar mass progenitors: ZTF18abhpgje (SN 2018eul), ZTF18abdpvnd (SN 2018dvf), ZTF18aawpcel (SN 2018cir), and ZTF18abddmrf (SN 2018dsx).