Aims. Stripped envelope (SE) supernovae are explosions of stars that have somehow lost most of their outer envelopes. We present the discovery and analyse the observations of the Type Ib supernova 2019odp (a.k.a. ZTF19abqwtfu) covering epochs within days of the explosion to late nebular phases at 360 d post-explosion.Methods. Our observations include an extensive set of photometric observations and low- to medium-resolution spectroscopic observations, both covering the complete observable time range. We analysed the data using analytic models for the recombination cooling emission of the early excess emission and the diffusion of the peak light curve. We expanded on existing methods to derive oxygen mass estimates from nebular phase spectroscopy, and briefly discuss progenitor models based on this analysis.Results. Our spectroscopic observations confirm the presence of He in the supernova ejecta and we thus (re)classify SN 2019odp as a Type Ib supernova. From the pseudo-bolometric light curve, we estimate a high ejecta mass of Mej ∼ 4 − 7 M⊙. The high ejecta mass, large nebular [O I]/[Ca II] line flux ratio (1.2 − 1.9), and an oxygen mass above ⪆0.5 M⊙ point towards a progenitor with a pre-explosion mass higher than 18 M⊙. Whereas a majority of analysed SE supernovae in the literature seem to have low ejecta masses, indicating stripping in a binary star system, SN 2019odp instead has parameters that are consistent with an origin in a single massive star. The compact nature of the progenitor (≲10 R⊙) suggests that a Wolf-Rayet star is the progenitor.

Stars with zero-age main sequence masses between 140 and 260 M_⊙ are thought to explode as pair-instability supernovae (PISNe). During their thermonuclear runaway, PISNe can produce up to several tens of solar masses of radioactive nickel, resulting in luminous transients similar to some superluminous supernovae (SLSNe). Yet, no unambiguous PISN has been discovered so far. SN 2018ibb is a hydrogen-poor SLSN at z = 0.166 that evolves extremely slowly compared to the hundreds of known SLSNe. Between mid 2018 and early 2022, we monitored its photometric and spectroscopic evolution from the UV to the near-infrared (NIR) with 2–10 m class telescopes. SN 2018ibb radiated > 3 × 10⁵¹ erg during its evolution, and its bolometric light curve reached > 2 × 10⁴⁴ erg s⁻¹ at its peak. The long-lasting rise of > 93 rest-frame days implies a long diffusion time, which requires a very high total ejected mass. The PISN mechanism naturally provides both the energy source (⁵⁶Ni) and the long diffusion time. Theoretical models of PISNe make clear predictions as to their photometric and spectroscopic properties. SN 2018ibb complies with most tests on the light curves, nebular spectra and host galaxy, and potentially all tests with the interpretation we propose. Both the light curve and the spectra require 25–44 M_⊙ of freshly nucleosynthesised ⁵⁶Ni, pointing to the explosion of a metal-poor star with a helium core mass of 120–130 M_⊙ at the time of death. This interpretation is also supported by the tentative detection of [Co II] λ 1.025 μm, which has never been observed in any other PISN candidate or SLSN before. We observe a significant excess in the blue part of the optical spectrum during the nebular phase, which is in tension with predictions of existing PISN models. However, we have compelling observational evidence for an eruptive mass-loss episode of the progenitor of SN 2018ibb shortly before the explosion, and our dataset reveals that the interaction of the SN ejecta with this oxygen-rich circumstellar material contributed to the observed emission. That may explain this specific discrepancy with PISN models. Powering by a central engine, such as a magnetar or a black hole, can be excluded with high confidence. This makes SN 2018ibb by far the best candidate for being a PISN, to date.

Several types of energetic supernovae, such as superluminous supernovae (SLSNe) and broad-line Ic supernovae (Ic-BL SNe), could be powered by the spin-down of a rapidly rotating magnetar. Currently, most models used to infer the parameters for potential magnetar-driven supernovae make several unsuitable assumptions that likely bias the estimated parameters. In this work, we present a new model for magnetar-driven supernovae that relaxes several of these assumptions and an inference workflow that enables accurate estimation of parameters from light curves of magnetar-driven supernovae. In particular, in this model, we include the dynamical evolution of the ejecta, coupling it to the energy injected by the magnetar itself while also allowing for non-dipole spin down. We show that the model can reproduce SLSN and Ic-BL SN light curves consistent with the parameter space from computationally expensive numerical simulations. We also show the results of parameter inference on four well-known example supernovae, demonstrating the model’s effectiveness at capturing the considerable diversity in magnetar-driven supernova light curves. The model fits each light curve well and recovers parameters broadly consistent with previous works. This model will allow us to explore the full diversity of magnetar-driven supernovae under one theoretical framework, more accurately characterize these supernovae from only photometric data, and make more accurate predictions of future multiwavelength emission to test the magnetar-driven scenario better.

One of the open questions following the discovery of GW170817 is whether neutron star (NS) mergers are the only astrophysical sites capable of producing r-process elements. Simulations have shown that 0.01–0.1 M_⊙ of r-process material could be generated in the outflows originating from the accretion disk surrounding the rapidly rotating black hole that forms as a remnant to both NS mergers and collapsing massive stars associated with long-duration gamma-ray bursts (collapsars). The hallmark signature of r-process nucleosynthesis in the binary NS merger GW170817 was its long-lasting near-infrared (NIR) emission, thus motivating a systematic photometric study of the light curves of broad-lined stripped-envelope (Ic-BL) supernovae (SNe) associated with collapsars. We present the first systematic study of 25 SNe Ic-BL—including 18 observed with the Zwicky Transient Facility and 7 from the literature—in the optical/NIR bands to determine what quantity of r-process material, if any, is synthesized in these explosions. Using semi-analytic models designed to account for r-process production in SNe Ic-BL, we perform light curve fitting to derive constraints on the r-process mass for these SNe. We also perform independent light curve fits to models without the r-process. We find that the r-process-free models are a better fit to the light curves of the objects in our sample. Thus, we find no compelling evidence of r-process enrichment in any of our objects. Further high-cadence infrared photometric studies and nebular spectroscopic analysis would be sensitive to smaller quantities of r-process ejecta mass or indicate whether all collapsars are completely devoid of r-process nucleosynthesis.

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 the r band of mag. Using spectra obtained around peak light, we compute expansion velocities from the Fe ii 5169 Å line for each event with high enough signal-to-noise ratio spectra, and find an average value of 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 , , and 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.

Fulfilling the rich promise of rapid advances in time-domain astronomy is only possible through confronting our observations with physical models and extracting the parameters that best describe what we see. Here, we introduce REDBACK; a Bayesian inference software package for electromagnetic transients. REDBACK provides an object-orientated PYTHON interface to over 12 different samplers and over 100 different models for kilonovae, supernovae, gamma-ray burst afterglows, tidal disruption events, engine-driven transients among other explosive transients. The models range in complexity from simple analytical and semi-analytical models to surrogates built upon numerical simulations accelerated via machine learning. REDBACK also provides a simple interface for downloading and processing data from various catalogues such as Swift and FINK. The software can also serve as an engine to simulate transients for telescopes such as the Zwicky Transient Facility and Vera Rubin with realistic cadences, limiting magnitudes, and sky coverage or a hypothetical user-constructed survey or a generic transient for target-of-opportunity observations with different telescopes. As a demonstration of its capabilities, we show how REDBACK can be used to jointly fit the spectrum and photometry of a kilonova, enabling a more powerful, holistic probe into the properties of a transient. We also showcase general examples of how REDBACK can be used as a tool to simulate transients for realistic surveys, fit models to real, simulated, or private data, multimessenger inference with gravitational waves, and serve as an end-to-end software toolkit for parameter estimation and interpreting the nature of electromagnetic transients.

One of the most studied objects in astronomy, the Crab Nebula, is the remnant of the historical supernova SN 1054. Historical observations of the supernova imply a typical supernova luminosity, but contemporary observations of the remnant imply a low explosion energy and low ejecta kinetic energy. These observations are incompatible with a standard Ni-56-powered supernova, hinting at an an alternate power source such as circumstellar interaction or a central engine. We examine SN 1054 using a pulsar-driven supernova model, similar to those used for superluminous supernovae. The model can reproduce the luminosity and velocity of SN 1054 for an initial spin period of similar to 14 ms and an initial dipole magnetic field of 10(14-15) G. We discuss the implications of these results, including the evolution of the Crab pulsar, the evolution of the remnant structure, formation of filaments, and limits on freely expanding ejecta. We discuss how our model could be tested further through potential light echo photometry and spectroscopy, as well as the modern analogues of SN 1054.

SN 2020zbf is a hydrogen-poor superluminous supernova (SLSN) at z = 0.1947 that shows conspicuous C II features at early times, in contrast to the majority of H-poor SLSNe. Its peak magnitude is M_g = −21.2 mag and its rise time (≲26.4 days from first light) places SN 2020zbf among the fastest rising type I SLSNe. We used spectra taken from ultraviolet (UV) to near-infrared wavelengths to identify spectral features. We paid particular attention to the C II lines as they present distinctive characteristics when compared to other events. We also analyzed UV and optical photometric data and modeled the light curves considering three different powering mechanisms: radioactive decay of ⁵⁶Ni, magnetar spin-down, and circumstellar medium (CSM) interaction. The spectra of SN 2020zbf match the model spectra of a C-rich low-mass magnetar-powered supernova model well. This is consistent with our light curve modeling, which supports a magnetar-powered event with an ejecta mass M_ej = 1.5 M_⊙. However, we cannot discard the CSM-interaction model as it may also reproduce the observed features. The interaction with H-poor, carbon-oxygen CSM near peak light could explain the presence of C II emission lines. A short plateau in the light curve around 35–45 days after peak, in combination with the presence of an emission line at 6580 Å, can also be interpreted as being due to a late interaction with an extended H-rich CSM. Both the magnetar and CSM-interaction models of SN 2020zbf indicate that the progenitor mass at the time of explosion is between 2 and 5 M_⊙. Modeling the spectral energy distribution of the host galaxy reveals a host mass of 10^8.7 M_⊙, a star formation rate of 0.24_−0.12^+0.41 M_⊙ yr⁻¹, and a metallicity of ∼0.4 Z_⊙.

Type Ia supernovae (SNe Ia) are thermonuclear explosions of degenerate white dwarf stars destabilized by mass accretion from a companion star¹, but the nature of their progenitors remains poorly understood. A way to discriminate between progenitor systems is through radio observations; a non-degenerate companion star is expected to lose material through winds² or binary interaction³ before explosion, and the supernova ejecta crashing into this nearby circumstellar material should result in radio synchrotron emission. However, despite extensive efforts, no type Ia supernova (SN Ia) has ever been detected at radio wavelengths, which suggests a clean environment and a companion star that is itself a degenerate white dwarf star^4,5. Here we report on the study of SN 2020eyj, a SN Ia showing helium-rich circumstellar material, as demonstrated by its spectral features, infrared emission and, for the first time in a SN Ia to our knowledge, a radio counterpart. On the basis of our modelling, we conclude that the circumstellar material probably originates from a single-degenerate binary system in which a white dwarf accretes material from a helium donor star, an often proposed formation channel for SNe Ia (refs. ^6,7). We describe how comprehensive radio follow-up of SN 2020eyj-like SNe Ia can improve the constraints on their progenitor systems.

The dividing line between gamma-ray bursts (GRBs) and ordinary stripped-envelope core-collapse supernovae (SNe) is yet to be fully understood. Observationally mapping the variety of ejecta outcomes (ultrarelativistic, mildly relativistic, or nonrelativistic) in SNe of Type Ic with broad lines (Ic-BL) can provide a key test to stellar explosion models. However, this requires large samples of the rare SN Ic-BL events with follow-up observations in the radio, where fast ejecta can be probed largely free of geometry and viewing angle effects. Here, we present the results of a radio (and X-ray) follow-up campaign of 16 SNe Ic-BL detected by the Zwicky Transient Facility (ZTF). Our radio campaign resulted in four counterpart detections and 12 deep upper limits. None of the events in our sample is as relativistic as SN 1998bw and we constrain the fraction of SN 1998bw-like explosions to <19% (3σ Gaussian equivalent), a factor of ≈2 smaller than previously established. We exclude relativistic ejecta with radio luminosity densities in between ≈5 × 10²⁷ erg s⁻¹ Hz⁻¹ and ≈10²⁹ erg s⁻¹ Hz⁻¹ at t ≳ 20 days since explosion for ≈60% of the events in our sample. This shows that SNe Ic-BL similar to the GRB-associated SNe 1998bw, 2003lw, and 2010bh, or to the relativistic SNe 2009bb and iPTF17cw, are rare. Our results also exclude an association of the SNe Ic-BL in our sample with largely off-axis GRBs with energies E ≳ 10⁵⁰ erg. The parameter space of SN 2006aj-like events (faint and fast-peaking radio emission) is, on the other hand, left largely unconstrained, and systematically exploring it represents a promising line of future research.