Aims. We investigate the photometric characteristics of a sample of intermediate-luminosity red transients (ILRTs), a class of elusive objects with peak luminosity between that of classical novae and standard supernovae. Our goal is to provide a stepping stone in the path to reveal the physical origin of such events, thanks to the analysis of the datasets collected. Methods. We present the multi-wavelength photometric follow-up of four ILRTs, namely NGC 300 2008OT-1, AT 2019abn, AT 2019ahd, and AT 2019udc. Through the analysis and modelling of their spectral energy distribution and bolometric light curves, we inferred the physical parameters associated with these transients. Results. All four objects display a single-peaked light curve which ends in a linear decline in magnitudes at late phases. A flux excess with respect to a single blackbody emission is detected in the infrared domain for three objects in our sample, a few months after maximum. This feature, commonly found in ILRTs, is interpreted as a sign of dust formation. Mid-infrared monitoring of NGC 300 2008OT-1 761 days after maximum allowed us to infer the presence of ∼10-3-10-5 M⊙ of dust, depending on the chemical composition and the grain size adopted. The late-time decline of the bolometric light curves of the considered ILRTs is shallower than expected for 56Ni decay, hence requiring an additional powering mechanism. James Webb Space Telescope observations of AT 2019abn prove that the object has faded below its progenitor luminosity in the mid-infrared domain, five years after its peak. Together with the disappearance of NGC 300 2008OT-1 in Spitzer images seven years after its discovery, this supports the terminal explosion scenario for ILRTs. With a simple semi-analytical model we tried to reproduce the observed bolometric light curves in the context of a few solar masses ejected at few 103 km s-1 and enshrouded in an optically thick circumstellar medium.

Aims. We investigate the spectroscopic characteristics of intermediate-luminosity Red Transients (ILRTs), a class of elusive objects with peak luminosity between that of classical novae and standard supernovae. Our goal is to provide a stepping stone in the path to unveiling the physical origin of these events based on the analysis of the collected datasets. Methods. We present the extensive optical and near-infrared (NIR) spectroscopic monitoring of four ILRTs, namely NGC 300 2008OT-1, AT 2019abn, AT 2019ahd and AT 2019udc. First we focus on the evolution of the most prominent spectral features observed in the low-resolution spectra. We then present a more detailed description of the high-resolution spectrum collected for NGC 300 2008OT-1 with the Very Large Telescope equipped with UVES. Finally, we describe our analysis of late-time spectra of NGC 300 2008OT-1 and AT 2019ahd through comparisons with both synthetic and observed spectra. Results. Balmer and Ca lines dominate the optical spectra, revealing the presence of slowly moving circumstellar medium (CSM) around the objects. The line luminosity of Hα, Hβ, and Ca II NIR triplet presents a double peaked evolution with time, possibly indicative of interaction between fast ejecta and the slow CSM. The high-resolution spectrum of NGC 300 2008OT-1 reveals a complex circumstellar environment, with the transient being surrounded by a slow (∼30 km s-1) progenitor wind. At late epochs, optical spectra of NGC 300 2008OT-1 and AT 2019ahd show broad (∼2500 km s-1) emission features at ∼6170 Å and ∼7000 Å which are unprecedented for ILRTs. We find that these lines originate most likely from the blending of several narrow lines, possibly of iron-peak elements.

We present photometric and spectroscopic observations of SN 2020xga and SN 2022xgc, two hydrogen-poor superluminous supernovae (SLSNe-I) at z=-0.4296 and z = 0.3103, respectively, which show an additional set of broad Mg II absorption lines, blueshifted by a few thousands kilometer second-1 with respect to the host galaxy absorption system. Previous work interpreted this as due to resonance line scattering of the SLSN continuum by rapidly expanding circumstellar material (CSM) expelled shortly before the explosion. The peak rest-frame g-band magnitude of SN 2020xga is -22.30 ± 0.04 mag and of SN 2022xgc is -21.97 ± 0.05 mag, placing them among the brightest SLSNe-I. We used high-quality spectra from ultraviolet to near-infrared wavelengths to model the Mg II line profiles and infer the properties of the CSM shells. We find that the CSM shell of SN 2020xga resides at ∼1.3×1016 cm, moving with a maximum velocity of 4275 km s-1, and the shell of SN 2022xgc is located at ∼0.8×1016 cm, reaching up to 4400 km s-1. These shells were expelled ∼11 and ∼5 months before the explosions of SN 2020xga and SN 2022xgc, respectively, possibly as a result of luminous-blue-variable-like eruptions or pulsational pair instability (PPI) mass loss. We also analyzed optical photometric data and modeled the light curves, considering powering from the magnetar spin-down mechanism. The results support very energetic magnetars, approaching the mass-shedding limit, powering these SNe with ejecta masses of ∼7-9M⊙. The ejecta masses inferred from the magnetar modeling are not consistent with the PPI scenario pointing toward stars > 50M⊙ He-core; hence, alternative scenarios such as fallback accretion and CSM interaction are discussed. Modeling the spectral energy distribution of the host galaxy of SN 2020xga reveals a host mass of 107.8 M⊙, a star formation rate of 0.96-0.26+0.47 M⊙ yr-1, and a metallicity of ∼0.2 Z⊙

Stars are initially powered by the fusion of hydrogen to helium. These ashes serve as fuel in a series of stages1, 2–3, transforming massive stars into a structure of shells. These are composed of natal hydrogen on the outside and consecutively heavier compositions inside, predicted to be dominated by He, C/O, O/Ne/Mg and O/Si/S (refs. 4,5). Silicon and sulfur are fused into iron, leading to the collapse of the core and either a supernova explosion or the formation of a black hole6, 7, 8–9. Stripped stars, in which the outer hydrogen layer has been removed and the internal He-rich or even the C/O layer below it is exposed10, provide evidence for this shell structure and the cosmic element production mechanism it reflects. The supernova types that arise from stripped stars embedded in shells of circumstellar material (CSM) confirm this scenario11, 12, 13, 14–15. However, direct evidence for the most interior shells, which are responsible for producing elements heavier than oxygen, is lacking. Here we report the discovery of the supernova (SN) 2021yfj resulting from a star stripped to its O/Si/S-rich layer. We directly observe a thick, massive Si/S-rich shell, expelled by the progenitor shortly before the supernova explosion. Exposing such an inner stellar layer is theoretically challenging and probably requires a rarely observed mass-loss mechanism. This rare supernova event reveals advanced stages of stellar evolution, forming heavier elements, including silicon, sulfur and argon, than those detected on the surface of any known class of massive stars.

We present panchromatic optical + near-infrared (NIR) + mid-infrared (MIR) observations of the intermediate-luminosity Type Iax supernova (SN Iax) 2024pxl and the extremely low-luminosity SN Iax 2024vjm. JWST observations provide unprecedented MIR spectroscopy of SN Iax, spanning from +11 to +42 day past maximum light. We detect forbidden emission lines in the MIR at these early times while the optical and NIR are dominated by permitted lines with an absorption component. Panchromatic spectra at early times can thus simultaneously show nebular and photospheric lines, probing both inner and outer layers of the ejecta. We identify spectral lines not seen before in SN Iax, including [Mg ii] 4.76 μm, [Mg ii] 9.71 μm, [Ne ii] 12.81 μm, and isolated O i 2.76 μm that traces unburned material. Forbidden emission lines of all species are centrally peaked with similar kinematic distributions, indicating that the ejecta are well mixed in both SN 2024pxl and SN 2024vjm, a hallmark of pure deflagration explosion models. Radiative transfer modeling of SN 2024pxl shows good agreement with a weak deflagration of a near-Chandrasekhar-mass white dwarf, but additional IR flux is needed to match the observations, potentially attributable to a surviving remnant. Similarly, we find SN 2024vjm is also best explained by a weak deflagration model, despite the large difference in luminosity between the two supernovae. Future modeling should push to even weaker explosions and include the contribution of a bound remnant. Our observations demonstrate the diagnostic power of panchromatic spectroscopy for unveiling explosion physics in thermonuclear supernovae.

We present the luminosity function and volumetric rate of a sample of Type IIP supernovae (SNe) from the Zwicky Transient Facility Census of the Local Universe survey (CLU). This is the largest sample of Type IIP SNe from a systematic volume-limited survey to-date. The final sample includes 330 Type IIP SNe and 36 low-luminosity Type II (LLIIP) SNe with M_r,peak > −16 mag, which triples the literature sample of LLIIP SNe. The fraction of LLIIP SNe is 19−4+3% of the total CLU Type IIP SNe population (8−2+1% of all core-collapse SNe). This implies that while LLIIP SNe likely represent the fate of core-collapse SNe of 8–12 M_⊙ progenitors, they alone cannot account for the fate of all massive stars in this mass range. To derive an absolute rate, we estimate the ZTF pipeline efficiency as a function of the apparent magnitude and the local surface brightness. We derive a volumetric rate of (3.9−0.4+0.4)×104Gpc−3yr−1 for Type IIP SNe and (7.3−0.6+0.6)×103Gpc−3yr−1 for LLIIP SNe. Now that the rate of LLIIP SNe is robustly derived, the unresolved discrepancy between core-collapse SN rates and star formation rates cannot be explained by LLIIP SNe alone.

Context. Hydrogen-rich superluminous supernovae (SLSNe II) are rare. The exact mechanism producing their extreme light curve peaks is not understood. Analysis of single events and small samples suggest that circumstellar material (CSM) interaction is the main mechanism responsible for the observed features. However, other mechanisms cannot be discarded. Large sample analysis can provide clarification.

Aims. We aim to characterize the light curves of a sample of 107 SLSNe II to provide valuable information that can be used to validate theoretical models.

Methods. We analyzed the gri light curves of SLSNe II obtained through ZTF. We studied the peak absolute magnitudes and characteristic timescales. When possible, we computed the g − r colors and pseudo-bolometric light curves, and estimated lower limits for their total radiated energy. We also studied the luminosity distribution of our sample and estimated the fraction that would be observable by the LSST. Finally, we compared our sample to other H-rich SNe and to H-poor SLSNe I.

Results. SLSNe II are heterogeneous. Their median peak absolute magnitude is ∼ − 20.3 mag in optical bands. Their rise can take from ∼two weeks to over three months, and their decline times range from ∼twenty days to over a year. We found no significant correlations between peak magnitude and timescales. SLSNe II tend to show fainter peaks, longer declines, and redder colors than SLSNe I.

Conclusions. We present the largest sample of SLSN II light curves to date, comprising 107 events. Their diversity could be explained by different CSM morphologies, although theoretical analysis is needed to explore alternative scenarios. Other luminous transients, such as active galactic nuclei, tidal disruption events or SNe Ia-CSM, can easily become contaminants. Thus, good multiwavelength light curve coverage becomes paramount. LSST could miss ∼30% of the ZTF events in its gri band footprint.

We present a comprehensive photometric and spectroscopic study of the Type IIP supernova (SN) 2018is. The V band luminosity and the expansion velocity at 50 days post-explosion are −15.1 ± 0.2 mag (corrected for AV = 1.34 mag) and 1400 km s−1, classifying it as a low-luminosity SN II. The recombination phase in the V band is shorter, lasting around 110 days, and exhibits a steeper decline (1.0 mag per 100 days) compared to most other low-luminosity SNe II. Additionally, the optical and near-infrared spectra display hydrogen emission lines that are strikingly narrow, even for this class. The Fe ii and Sc ii line velocities are at the lower end of the typical range for low-luminosity SNe II. Semi-analytical modelling of the bolometric light curve suggests an ejecta mass of ∼8 M, corresponding to a pre-supernova mass of ∼9.5 M, and an explosion energy of ∼0.40 × 1051 erg. Hydrodynamical modelling further indicates that the progenitor had a zero-age main sequence mass of 9 M, coupled with a low explosion energy of 0.19 × 1051 erg. The nebular spectrum reveals weak [O i] λλ6300,6364 lines, consistent with a moderate-mass progenitor, while features typical of Fe core-collapse events, such as He i, [C i], and Fe i, are indiscernible. However, the redder colours and low ratio of Ni to Fe abundance do not support an electron-capture scenario either. As a low-luminosity SN II with an atypically steep decline during the photospheric phase and remarkably narrow emission lines, SN 2018is contributes to the diversity observed within this population.

We present the long-term photometric and spectroscopic analysis of a transitioning SN IIn/Ibn from –10.8 d to 150.7 d post V-band maximum. SN 2021foa shows prominent He I lines comparable in strength to the H α line around peak, placing SN 2021foa between the SN IIn and SN Ibn populations. The spectral comparison shows that it resembles the SN IIn population at pre-maximum, becomes intermediate between SNe IIn/Ibn, and at post-maximum matches with SN IIn 1996al. The photometric evolution shows a precursor at –50 d and a light curve shoulder around 17 d. The peak luminosity and colour evolution of SN 2021foa are consistent with most SNe IIn and Ibn in our comparison sample. SN 2021foa shows the unique case of an SN IIn where the narrow P-Cygni in H α becomes prominent at 7.2 d. The H α profile consists of a narrow (500–1200 km s-1) component, intermediate width (3000–8000 km s-1) and broad component in absorption. Temporal evolution of the H α profile favours a disc-like CSM geometry. Hydrodynamical modelling of the light curve well reproduces a two-component CSM structure with different densities (ρ ∝ r-2–ρ ∝ r-5), mass-loss rates (10-3–10-1 M☉ yr-1) assuming a wind velocity of 1000 km s-1 and having a CSM mass of 0.18 M☉. The overall evolution indicates that SN 2021foa most likely originated from an LBV star transitioning to a WR star with the mass-loss rate increasing in the period from 5 to 0.5 yr before the explosion or it could be due to a binary interaction.

Type IIn supernovae (SNe) resembling SN 2009ip (09ip-like SNe) originate from the interaction between circumstellar material (CSM) and the ejecta. This subclass not only shares similar observational properties around the maximum but is also commonly characterized by a long-duration precursor before its maximum. Investigating the observed properties of the precursor provides constraints on the mass-loss history of the progenitor. We present observational data of SN 2023vbg, an 09ip-like type IIn SN that displayed unique observational properties compared to other 09ip-like SNe. SN 2023vbg showed a long-duration precursor at M_g ∼ −14 mag lasting for ∼100 days, followed by a bright bump at M_g ∼ −17 mag at 12–25 days before the maximum. The luminosity of the precursor is similar to those of other 09ip-like SNe, but the bright bump has not been observed in other cases. After reaching the peak luminosity, the light curve exhibited a relative smooth decline. While the Hα profile displays two velocity components (∼500 and 3000 km s⁻¹), a broad component observed in other 09ip-like SNe was not seen, but it may emerge later. We suggest that these properties are explained by the difference in the CSM structure as compared to other 09ip-like SNe; SN 2023vbg had an inner denser CSM component, as well as generally smooth CSM density distribution in a more extended scale than in the others. Such diversity of CSM likely reflects the diversity of pre-SN outbursts, which in turn may mirror the range of evolutionary pathways in the final stages of the progenitors.