Search for: All records

We identify the progenitor star of SN 2023ixf in Messier 101 using Keck/NIRC2 adaptive optics imaging and pre-explosion Hubble Space Telescope (HST)/Advanced Camera for Surveys (ACS) images. The supernova, localized with diffraction spikes and high-precision astrometry, unambiguously coincides with a progenitor candidate of $m_\text{F814W}=24.87\pm 0.05$ (AB). Given its reported infrared excess and semiregular variability, we fit a time-dependent spectral energy distribution (SED) model of a dusty red supergiant (RSG) to a combined data set of HST optical, ground-based near-infrared, and Spitzer Infrared Array Camera (IRAC) [3.6], [4.5] photometry. The progenitor resembles an RSG of $T_\text{eff}=3488\pm 39$ K and $\log (L/\mathrm{L}_\odot)=5.15\pm 0.02$, with a $0.13\pm 0.01$ dex ($31.1\pm 1.7$ per cent) luminosity variation at a period of $P=1144.7\pm 4.8$ d, obscured by a dusty envelope of $\tau =2.92\pm 0.02$ at $1\, \mu \text{m}$ in optical depth (or $A_\text{V}=8.43\pm 0.11$ mag). The signatures match a post-main-sequence star of $18.2_{-0.6}^{+1.3}\, \mathrm{M}_\odot$ in zero-age main-sequence mass, among the most massive SN II progenitor, with a pulsation-enhanced mass-loss rate of $\dot{M}=(4.32\pm 0.26)\times 10^{-4} \, \mathrm{M}_\odot \, \text{yr}^{-1}$. The dense and confined circumstellar material is ejected during the last episode of radial pulsation before the explosion. Notably, we find strong evidence for variations of $\tau$ or $T_\text{eff}$ along with luminosity, a necessary assumption to reproduce the wavelength-dependent variability, which implies periodic dust sublimation and condensation. Given the observed SED, partial dust obscuration remains possible, but any unobstructed binary companion over $5.6\, \mathrm{ M}_\odot$ can be ruled out.

The Bright Transient Survey (BTS) aims to obtain a classification spectrum for all bright (m_peak≤ 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 ∼5 yr of ZTF operations. We presentBTSbot, a multimodal convolutional neural network, which provides a bright transient score to individual ZTF detections using their image data and 25 extracted features.BTSbotis able to eliminate the need for daily human scanning by automatically identifying and requesting spectroscopic follow-up observations of new bright transient candidates.BTSbotrecovers all bright transients in our test split and performs on par with scanners in terms of identification speed (on average, ∼1 hr quicker than scanners). We also find thatBTSbotis not significantly impacted by any data shift by comparing performance across a concealed test split and a sample of very recent BTS candidates.BTSbothas been integrated intoFritzandKowalski, ZTF’s first-party marshal and alert broker, and now sends automatic spectroscopic follow-up requests for the new transients it identifies. Between 2023 December and 2024 May,BTSbotselected 609 sources in real time, 96% of which were real extragalactic transients. WithBTSbotand 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.

We present the tidal disruption event (TDE) AT2022lri, hosted in a nearby (≈144 Mpc) quiescent galaxy with a low-mass massive black hole (10⁴M_⊙<M_BH< 10⁶M_⊙). AT2022lri belongs to the TDE-H+He subtype. More than 1 Ms of X-ray data were collected with NICER, Swift, and XMM-Newton from 187 to 672 days after peak. The X-ray luminosity gradually declined from 1.5 × 10⁴⁴erg s⁻¹to 1.5 × 10⁴³erg s⁻¹and remains much above the UV and optical luminosity, consistent with a super-Eddington accretion flow viewed face-on. Sporadic strong X-ray dips atop a long-term decline are observed, with a variability timescale of ≈0.5 hr–1 days and amplitude of ≈2–8. When fitted with simple continuum models, the X-ray spectrum is dominated by a thermal disk component with inner temperature going from ∼146 to ∼86 eV. However, there are residual features that peak around 1 keV, which, in some cases, cannot be reproduced by a single broad emission line. We analyzed a subset of time-resolved spectra with two physically motivated models describing a scenario either where ionized absorbers contribute extra absorption and emission lines or where disk reflection plays an important role. Both models provide good and statistically comparable fits, show that the X-ray dips are correlated with drops in the inner disk temperature, and require the existence of subrelativistic (0.1–0.3c) ionized outflows. We propose that the disk temperature fluctuation stems from episodic drops of the mass accretion rate triggered by magnetic instabilities or/and wobbling of the inner accretion disk along the black hole’s spin axis.

Optical surveys have become increasingly adept at identifying candidate tidal disruption events (TDEs) in large numbers, but classifying these generally requires extensive spectroscopic resources. Here we presenttdescore, a simple binary photometric classifier that is trained using a systematic census of ∼3000 nuclear transients from the Zwicky Transient Facility (ZTF). The sample is highly imbalanced, with TDEs representing ∼2% of the total.tdescoreis nonetheless able to reject non-TDEs with 99.6% accuracy, yielding a sample of probable TDEs with recall of 77.5% for a precision of 80.2%.tdescoreis thus substantially better than any available TDE photometric classifier scheme in the literature, with performance not far from spectroscopy as a method for classifying ZTF nuclear transients, despite relying solely on ZTF data and multiwavelength catalog cross matching. In a novel extension, we use “Shapley additive explanations” to provide a human-readable justification for each individualtdescoreclassification, enabling users to understand and form opinions about the underlying classifier reasoning.tdescorecan serve as a model for photometric identification of TDEs with time-domain surveys, such as the upcoming Rubin observatory.

We study a magnitude-limited sample of 36 broad-lined type Ic supernovae (SNe Ic-BL) from the Zwicky Transient Facility Bright Transient Survey (detected between 2018 March and 2021 August), which is the largest systematic study of SNe Ic-BL done in literature thus far. We present the light curves (LCs) for each of the SNe and analyze the shape of the LCs to derive empirical parameters, along with the explosion epochs for every event. The sample has an average absolute peak magnitude in therband of${\overline{\mathit{M}}}_{\mathrm{r},\max }=-18.51\pm 0.15$mag. Using spectra obtained around peak light, we compute expansion velocities from the Feii5169 Å line for each event with high enough signal-to-noise ratio spectra, and find an average value of$\overline{v_{\mathrm{ph}}}=16,100\pm 1100$km s⁻¹. We also compute bolometric LCs, study the blackbody temperature and radii evolution over time, and derive the explosion properties of the SNe. The explosion properties of the sample have average values of${\overline{\mathit{M}}}_{\mathrm{Ni}}={0.37}_{-0.06}^{+0.08}{\mathit{M}}_{\odot }$,${\overline{\mathit{M}}}_{\mathrm{ej}}={2.45}_{-0.41}^{+0.47}{\mathit{M}}_{\odot }$, and${\overline{\mathit{E}}}_{\mathrm{K}}=({4.02}_{-1.00}^{+1.37})\times {10}^{51}$erg. Thirteen events have radio observations from the Very Large Array, with eight detections and five non-detections. We find that the populations that have radio detections and radio non-detections are indistinct from one another with respect to their optically inferred explosion properties, and there are no statistically significant correlations present between the events’ radio luminosities and optically inferred explosion properties. This provides evidence that the explosion properties derived from optical data alone cannot give inferences about the radio properties of SNe Ic-BL and likely their relativistic jet formation mechanisms.

Eruptive mass loss of massive stars prior to supernova (SN) explosion is key to understanding their evolution and end fate. An observational signature of pre-SN mass loss is the detection of an early, short-lived peak prior to the radioactive-powered peak in the lightcurve of the SN. This is usually attributed to the SN shock passing through an extended envelope or circumstellar medium. Such an early peak is common for double-peaked Type IIb SNe with an extended hydrogen envelope but uncommon for normal Type Ibc SNe with very compact progenitors. In this paper, we systematically study a sample of 14 double-peaked Type Ibc SNe out of 475 Type Ibc SNe detected by the Zwicky Transient Facility. The rate of these events is ∼3%–9% of Type Ibc SNe. A strong correlation is seen between the peak brightness of the first and the second peak. We perform a holistic analysis of this sample’s photometric and spectroscopic properties. We find that six SNe have ejecta mass less than 1.5M_⊙. Based on the nebular spectra and lightcurve properties, we estimate that the progenitor masses for these are less than ∼12M_⊙. The rest have an ejecta mass >2.4M_⊙and a higher progenitor mass. This sample suggests that the SNe with low progenitor masses undergo late-time binary mass transfer. Meanwhile, the SNe with higher progenitor masses are consistent with wave-driven mass loss or pulsation-pair instability-driven mass-loss simulations.

The origin of cosmic high-energy neutrinos remains largely unexplained. For high-energy neutrino alerts from IceCube, a coincidence with time-variable emission has been seen for three different types of accreting black holes: (1) a gamma-ray flare from a blazar (TXS 0506+056), (2) an optical transient following a stellar tidal disruption event (TDE; AT2019dsg), and (3) an optical outburst from an active galactic nucleus (AGN; AT2019fdr). For the latter two sources, infrared follow-up observations revealed a powerful reverberation signal due to dust heated by the flare. This discovery motivates a systematic study of neutrino emission from all supermassive black hole with similar dust echoes. Because dust reprocessing is agnostic to the origin of the outburst, our work unifies TDEs and high-amplitude flares from AGN into a population that we dub accretion flares. Besides the two known events, we uncover a third flare that is coincident with a PeV-scale neutrino (AT2019aalc). Based solely on the optical and infrared properties, we estimate a significance of 3.6σ for this association of high-energy neutrinos with three accretion flares. Our results imply that at least ∼10 per cent of the IceCube high-energy neutrino alerts could be due to accretion flares. This is surprising because the sum of the fluence of these flares is at least three orders of magnitude lower compared to the total fluence of normal AGN. It thus appears that the efficiency of high-energy neutrino production in accretion flares is increased compared to non-flaring AGN. We speculate that this can be explained by the high Eddington ratio of the flares.

During the first half of the fourth observing run (O4a) of the International Gravitational Wave Network, 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 (False Alarm Rate less than 1 yr⁻¹) BNS and NSBH candidates in O4a. Our follow-up campaigns relied on both target-of-opportunity observations and re-weighting of the nominal survey schedule to maximize coverage. We describe the toolkit we have been developing,Fritz, an instance ofSkyPortal, 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²within the skymaps of the high-significance GW events, reaching a median depth ofg≈ 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 events), along with our O4a follow-ups. Following Kasliwal et al., 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.

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 full outburst of 320 days in 1999, and the observed peak luminosity was 2 orders of magnitude lower. Its peak bolometric luminosity was only 2 × 10³⁵erg s⁻¹, implying an Eddington fraction of about 3 × 10⁻⁴. The source remained in the hard spectral state throughout the outburst. From optical spectroscopy measurements we estimate an outer disk radius of 10¹¹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 light curve. These observations highlight the potential of optical and infrared synoptic surveys to discover low-luminosity activity from X-ray transients.

Multipeaked supernovae with precursors, dramatic light-curve rebrightenings, and spectral transformation are rare, but are being discovered in increasing numbers by modern night-sky transient surveys like the Zwicky Transient Facility. Here, we present the observations and analysis of SN 2023aew, which showed a dramatic increase in brightness following an initial luminous (−17.4 mag) and long (∼100 days) unusual first peak (possibly precursor). SN 2023aew was classified as a Type IIb supernova during the first peak but changed its type to resemble a stripped-envelope supernova (SESN) after the marked rebrightening. We present comparisons of SN 2023aew’s spectral evolution with SESN subtypes and argue that it is similar to SNe Ibc during its main peak. P-Cygni Balmer lines are present during the first peak, but vanish during the second peak’s photospheric phase, before Hαresurfaces again during the nebular phase. The nebular lines ([Oi], [Caii], Mgi], Hα) exhibit a double-peaked structure that hints toward a clumpy or nonspherical ejecta. We analyze the second peak in the light curve of SN 2023aew and find it to be broader than that of normal SESNe as well as requiring a very high⁵⁶Ni mass to power the peak luminosity. We discuss the possible origins of SN 2023aew including an eruption scenario where a part of the envelope is ejected during the first peak and also powers the second peak of the light curve through interaction of the SN with the circumstellar medium.