We present observations of SN 2023xgo, a transitional Type Ibn/Icn SN, from-5.6 to 63 d relative to r-band peak. Early spectra show C III lambda 5696 emission like Type Icn SNe, shifting to Type Ibn features. The He I velocities (1800-10 000 km s(-1)) and pseudo-equivalent widths are among the highest in the Ibn/Icn class. The light curve declines at 0.14 mag d(-1) until 30 d, matching SNe Ibn/Icn but slower than fast transients. SN 2023xgo is the faintest in our SN Ibn sample ( M-r = -17.65 +/- 0 . 04) but shows typical colour and host properties. Semi-analytical modelling of the light curve suggests a compact CSM shell (similar to 10(12)-10(13) cm), mass-loss rate between 10(-4) and 10(-3) M-circle dot yr(-1) with CSM and ejecta masses of similar to 0.22 and 0.12 M-circle dot, respectively. Post-maximum light-curve, spectral modelling favours a similar to 3 M-circle dot helium star progenitor with extended (similar to 10(15) cm), stratified CSM (density exponent of 2.9) and mass-loss rate of 0 . 1-2 . 7 M-circle dot yr(-1). These two mass-loss regimes imply a radially varying CSM, shaped by asymmetry or changes in the progenitor's mass-loss over time. This mass-loss behaviour fits both binary and single-star evolution. Early Icn-like features stem from hot carbon ionization, fading to Ibn-like with cooling. SN 2023xgo thus offers rare insight into the connection between SNe Icn, Ibn, and SNe Ibn with ejecta signatures.

The observational link between long gamma-ray bursts (GRBs) and broad-lined stripped-envelope core-collapse supernovae (SNe Ic-BL) is well established. Significant progress has been made in constraining what fraction of SNe Ic-BL may power high- or low-luminosity GRBs when viewed at small off-axis angles. However, the GRB-SN connection still lacks a complete understanding in the broader context of massive-star evolution and explosion physics. Models predict a continuum of outcomes for the fastest ejecta, from choked to ultrarelativistic jets, and observations from radio to X-rays are key to probing these scenarios across a range of viewing angles and velocities. Here, we present results from a coordinated radio-to-X-ray campaign targeting nearby (z less than or similar to 0.1) SNe Ic-BL designed to explore this diversity. With eight new radio-monitored events and updated data for one previously observed SN, we further tighten constraints on the fraction of SNe Ic-BL as relativistic as SN 1998bw (GRB 980425). We identify SN 2024rjw as a new radio-loud event likely powered by strong interaction with circumstellar material, and add evidence supporting a similar interpretation for SN 2020jqm. We also establish new limits on the properties of radio-emitting ejecta with velocities consistent with cocoons from choked jets, highlighting SN 2022xxf as a promising cocoon-dominated candidate. These results refine our understanding of the continuum linking ordinary SNe Ic-BL, engine-driven explosions, and GRBs, and contribute to building a sample that will inform future multimessenger searches for electromagnetic counterparts to high-energy neutrinos.

The Zwicky Transient Facility Census of the Local Universe survey yielded a sample of 330 Type IIP supernovae (SNe) with well-constrained peak luminosities. In paper I, we measured their luminosity function and volumetric rate. Here (paper II), we present the largest systematic study of lightcurve properties for Type IIP SNe from a volume-limited survey, analyzing a selected subset of 129 events, including 16 low-luminosity Type IIP (LLIIP) SNe with Mr,peak >= -16 mag. We find that plateau slope correlates with peak brightness, with many LLIIP SNe showing positive slopes-suggesting smaller progenitor radii and distinct density profiles compared to brighter Type IIP SNe. The plateau duration shows only a weak dependence on peak brightness, likely suggesting binary interaction. One SN exhibits a plateau-to-tail drop of >3.5 mag, consistent with an electron-capture or failed SN with very low or zero nickel mass. We derive explosion and progenitor parameters of the entire Type IIP SN sample using semi-analytical and radiation-hydrodynamical models. Based on radiation-hydrodynamical model fitting, LLIIP SNe are characterized by low nickel masses (0.001-0.025 M-circle dot), low explosion energies (0.1-0.28 & times; 10(51 )erg), low ejecta masses ( 8.1(-1.7)(+0.8) M-circle dot), and ZAMS masses below 11 M-circle dot. In comparison, the full Type IIP SN sample spans a wider range with nickel masses (0.001-0.222 M-circle dot), explosion energies (0.10-4.43 & times; 10(51 )erg), ejecta masses (5.4-24.8 M-circle dot), and ZAMS masses (9.3-16.7 M-circle dot). We find strong correlations between peak brightness, explosion energy, and nickel mass that extend to the low-luminosity end. We conclude that LLIIP SNe represent the faint, low-energy end of the Type IIP population and originate from the lowest-mass core-collapse progenitors.

We present optical/UV photometric and spectroscopic observations, as well as X-ray and radio follow-up, of the extraordinary event AT2019cmw. With a peak bolometric luminosity of similar to 10(45.6) erg s(-1), it is one of the most luminous thermal transients ever discovered. Extensive spectroscopic follow-up post-peak showed only a featureless continuum throughout its evolution. This, combined with its nuclear location, blue colour at peak and lack of prior evidence of an AGN in its host lead us to interpret this event as a 'featureless' tidal disruption event (TDE). It displays photometric evolution atypical of most TDEs, cooling from similar to 30 to similar to 10kK in the first similar to 300 d post-peak, with potential implications for future photometric selection of candidate TDEs. No X-ray or radio emission is detected, placing constraints on the presence of on-axis jetted emission or a visible inner-accretion disc. Modelling the optical light curve with existing theoretical prescriptions, we find that AT2019cmw may be the result of the disruption of a star in the tens of solar masses by a supermassive black hole (SMBH). Combined with a lack of detectable star formation in its host galaxy, it could imply the existence of a localized region of star formation around the SMBH. This could provide a new window to probe nuclear star formation and the shape of the initial mass function (IMF) in close proximity to SMBHs out to relatively high redshifts.

We present the searches conducted with the Zwicky Transient Facility (ZTF) in response to S250206dm, a bona fide event with an online false alarm rate of one in 25 yr, detected by the International Gravitational Wave Network. Although the event is significant, the nature of the compact objects involved remains unclear, with at least one likely neutron star. ZTF covered 68% of the last refined Bilby localization region, though we did not identify any likely optical counterpart. We describe the ZTF strategy, potential candidates, and the observations that helped rule out candidates, including sources circulated by other collaborations. Similar to Ahumada et al., we perform a frequentist analysis, using simsurvey, as well as Bayesian analysis, using nimbus, to quantify the efficiency of our searches. We find that, given the nominal up-to-date distance to this event of 373 +/- 104 Mpc, our efficiencies are above 10% for KNe brighter than -17.5 absolute magnitude. Assuming the optical counterpart known as kilonova (KN) lies within the ZTF footprint, our limits constrain the brightest end of the KN parameter space. Through dedicated radiative transfer simulations of KNe from binary neutron star (BNS) and black hole-neutron star mergers, we exclude parts of the BNS KN parameter space. Up to 35% of the models with high wind ejecta mass (Mwind approximate to 0.13 M circle dot) are ruled out when viewed face-on ( cos theta obs=1.0 ). Finally, we present a joint analysis using the combined coverage from ZTF and the Gravitational Wave Multimessenger Dark Energy Camera Survey. The joint observations cover 73% of the Bilby localization region, and the combined efficiency has a stronger impact on rising and slowly fading models, allowing us to rule out 55% of the high-mass KN models viewed face-on.

We present our observations and analysis of SN 2023gpw, a hydrogen-rich superluminous supernova (SLSN II) with broad emission lines in its post-peak spectra. Unlike previously observed SLSNe II, its light curve suggests an abrupt drop during a solar conjunction between ∼80 and ∼180 d after the light curve peak, which is possibly analogous to a normal hydrogen-rich supernova (SN). Spectra taken at and before the peak show hydrogen and helium “flash” emission lines attributed to early interaction with a dense confined circumstellar medium (CSM). A well-observed ultraviolet excess appears as these lines disappear, also as a result of CSM interaction. The blackbody photosphere expands roughly at the same velocity throughout the observations, indicating little or no bulk deceleration. This velocity is much higher than what is seen in spectral lines, suggesting asymmetry in the ejecta. The high total radiated energy (≳9 × 1050 erg) and aforementioned lack of bulk deceleration in SN 2023gpw are difficult to reconcile with a neutrino-driven SN simply combined with efficient conversion from kinetic energy to emission through interaction. This suggests an additional energy source, such as a central engine. While magnetar-powered models qualitatively similar to SN 2023gpw exist, more modeling work is required to determine if they can reproduce the observed properties in combination with early interaction. The required energy might alternatively be provided by accretion onto a black hole created in the collapse of a massive progenitor star.

AT2022cmc is the first on-axis jetted tidal disruption event (TDE) to be discovered at optical wavelengths. The optically bright nature of AT2022cmc presents an unprecedented opportunity to place this jetted TDE in the context of the larger optically selected thermal TDE population and explore potential connections to optical TDE subclasses, particularly the class of luminous TDEs that lack optical spectral features. In this work we present late-time optical observations of AT2022cmc, both imaging and spectroscopy, that extend the optical dataset to ∼ 160 days from the first detection in the observed frame. The light curve clearly evolves from red to blue, which we interpret as a transition from a non-thermally dominated spectral energy distribution (SED) to thermally dominated SED. By accounting for the non-thermal emission evident in the optical SED at early times, we extract the properties of the thermal emission and compare to a sample of optically selected thermal TDEs. We find that the properties of AT2022cmc are consistent with previous correlations found for the evolution and properties of thermal TDEs, with the thermal properties of AT2022cmc aligning with the class of featureless and luminous TDEs. The confirmation of this similarity motivates the importance of prompt and multi-wavelength follow-up of featureless and luminous TDEs in order to further explore the connection they have with jetted TDEs.

In recent years, a class of stripped-envelope supernovae (SESNe) has emerged that show two distinct peaks in their light curves, where the first peak cannot be attributed to shock cooling emission. Such peculiar supernovae are often studied individually, explained by invoking some combination of powering mechanisms. However, they have seldom been discussed in the broader context of double-peaked SESNe. In this paper, we attempt to form a picture of the landscape of double-peaked SESNe and their powering mechanisms by adding two more objects-SN 2021uvy and SN 2022hgk. SN 2021uvy is a broad and luminous SN Ib with an unusually long rise of the first peak and constant color evolution with rising photospheric temperature during the second peak. Although its first peak is similar to that of SN 2019stc, their second peaks differ in properties, making it unique among double-peaked objects. SN 2022hgk shows striking photometric similarity to SN 2019cad and spectroscopic similarity to SN 2005bf, both of which have been suggested to be powered by a double-nickel distribution in their ejecta. We analyze their light curves and colors, compare them with a sample of other double-peaked published supernovae for which we have additional data, and analyze the light curve parameters of the sample. We observe a correlation (p-value similar to 0.025) between the peak absolute magnitudes of the first and second peaks. The sample shows variety in the photometric and spectroscopic properties, and thus no single definitive powering mechanism applies to the whole sample. However, sub-groups of similarity exist that can be explained by mechanisms like the double-nickel distribution, magnetar central engine, interaction, and fallback accretion. We also map out the duration between the peaks (Delta t21) versus the difference between peak absolute magnitudes (Delta M21) as a phase-space that could potentially delineate the most promising powering mechanisms for the double-peaked SESNe.

For over 25 yr, the origin of long-duration gamma-ray bursts (lGRBs) has been linked to the collapse of rotating massive stars. However, we have yet to pinpoint the stellar progenitor powering these transients. Moreover, the dominant engine powering the explosions remains open to debate. Observations of both lGRBs, supernovae associated with these GRBs, such as broad-line (BL) stripped-envelope (type Ic) supernovae (hereafter, Ic-BL), supernovae (SNe), and perhaps superluminous SNe, fast blue optical transients, and fast x-ray transients, may provide clues to both engines and progenitors. In this paper, we conduct a detailed study of the tight-binary formation scenario for lGRBs, comparing this scenario to other leading progenitor models. Combining this progenitor scenario with different lGRB engines, we can compare to existing data and make predictions for future observational tests. We find that the combination of the tight-binary progenitor scenario with the black hole accretion disk engine can explain lGRBs, low-luminosity GRBs, ultra-long GRBs, and Ic-BL. We discuss the various progenitor properties required for these different subclasses and note such systems would be future gravitational-wave merger sources. We show that the current literature on other progenitor-engine scenarios cannot explain all of these transient classes with a single origin, motivating additional work. We find that the tight-binary progenitor with a magnetar engine is excluded by existing observations. The observations can be used to constrain the properties of stellar evolution, the nature of the GRB, and the associated SN engines in lGRBs and Ic-BL. We discuss the future observations needed to constrain our understanding of these rare, but powerful, explosions.

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 1 9 − 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 ) × 1 0 4 Gpc − 3 yr − 1 for Type IIP SNe and ( 7 . 3 − 0.6 + 0.6 ) × 1 0 3 Gpc − 3 yr − 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.

Although all Type II supernovae (SNe) originate from massive stars possessing a hydrogen-rich envelope, their light-curve morphology is diverse, reflecting poorly characterized heterogeneity in the physical properties of their progenitor systems. Here, we present a detailed light-curve analysis of a magnitude-limited sample of 639 Type II SNe from the Zwicky Transient Facility Bright Transient Survey. Using Gaussian processes, we systematically measure empirical light-curve features (e.g. rise times, peak colours, and luminosities) in a robust sampling-independent manner. We focus on rise times as they are highly sensitive to pre-explosion progenitor properties, especially the presence of a dense circumstellar medium (CSM) shed by the progenitor in the years immediately pre-explosion. By correlating our feature measurements with physical parameters from an extensive grid of stella hydrodynamical models with varying progenitor properties (CSM structure, M-center dot, R-CSM, and M-ZAMS), we quantify the proportion of events with sufficient pre-explosion mass loss to significantly alter the initial light curve (roughly MCSM >= 10(-2.5)M(circle dot)) in a highly complete sample of 377 spectroscopically classified Type II SNe. We find that 67 +/- 6 per cent of observed SNe in our magnitude-limited sample show evidence for substantial CSM (M-CSM >= 10(-2.5)M(circle dot)) close to the progenitor (R-CSM < 10(15) cm) at the time of explosion. After applying a volumetric-correction, we find 36(-7)(+5) per cent of all Type II SN progenitors possess substantial CSM within 1015 cm at the time of explosion. This high fraction of progenitors with dense CSM, supported by photometric and spectroscopic evidence of previous SNe, reveals mass-loss rates significantly exceeding those measured in local group red supergiants or predicted by current theoretical models.

Astrophysical high-energy (TeV-PeV) neutrinos were first discovered in 2013, but their origin remains largely unknown. Here we present SN 2023uqf, a supernova found in coincidence with high-energy neutrino IC231004A, as part of a systematic optical follow-up program with the Zwicky Transient Facility. SN 2023uqf had a luminous and rapidly-evolving lightcurve, and spectroscopic observations indicated that the source was a Type Ibn supernova. Spectroscopic signatures confirm ongoing interaction between the supernova ejecta and a dense circumstellar medium, as expected for high-energy neutrino production in a core-collapse supernova. Given the rare nature of Type Ibn supernovae, SN 2023uqf is unlikely to have been discovered by chance over the course of our program (p=0.3

Classical gamma-ray bursts (GRBs) have two distinct emission episodes: prompt emission from ultrarelativistic ejecta and afterglow from shocked circumstellar material. While both components are extremely luminous in known GRBs, a variety of scenarios predict the existence of luminous afterglow emission with little or no associated high-energy prompt emission. We present AT 2019pim, the first spectroscopically confirmed afterglow with no observed high-energy emission to be identified. Serendipitously discovered during follow-up observations of a gravitational-wave trigger and located in a contemporaneous TESS sector, it is hallmarked by a fast-rising (t approximate to 2 h), luminous (M-UV,M-peak approximate to-24.4 mag) optical transient with accompanying luminous X-ray and radio emission. No gamma-ray emission consistent with the time and location of the transient was detected by Fermi-GBM or by Konus, placing constraining limits on an accompanying GRB. We investigate several independent observational aspects of the afterglow in the context of constraints on relativistic motion and find all of them are consistent with an initial Lorentz factor of Gamma(0)approximate to 10-30 for the on-axis material, significantly lower than in any well-observed GRB and consistent with the theoretically predicted 'dirty fireball' scenario in which the high-energy prompt emission is stifled by pair production. However, we cannot rule out a structured jet model in which only the line-of-sight material was ejected at low-Gamma, off-axis from a classical high-Gamma jet core, and an on-axis GRB with below-average gamma-ray efficiency also remains a possibility. This event represents a milestone in orphan afterglow searches, demonstrating that luminous optical afterglows lacking detected GRB counterparts can be identified and spectroscopically confirmed in real time.

In the past few years, the improved sensitivity and cadence of wide-field optical surveys have enabled the discovery of several afterglows without associated detected gamma-ray bursts (GRBs). We present the identification, observations, and multiwavelength modeling of a recent such afterglow (AT 2023lcr), and model three literature events (AT 2020blt, AT 2021any, and AT 2021lfa) in a consistent fashion. For each event, we consider the following possibilities as to why a GRB was not observed: (1) the jet was off-axis; (2) the jet had a low initial Lorentz factor; and (3) the afterglow was the result of an on-axis classical GRB (on-axis jet with physical parameters typical of the GRB population), but the emission was undetected by gamma-ray satellites. We estimate all physical parameters using afterglowpy and Markov Chain Monte Carlo methods from emcee . We find that AT 2023lcr, AT 2020blt, and AT 2021any are consistent with on-axis classical GRBs, and AT 2021lfa is consistent with both on-axis low Lorentz factor (Γ _0 ≈ 5–13) and off-axis ( θ _obs = 2 θ _jet ) high Lorentz factor (Γ _0 ≈ 100) jets.

Near-sun sky twilight observations allow for the detection of asteroids interior to the orbit of Venus (Aylos) and the Earth (Atiras) and comets. We present the results of observations with the Palomar 48-inch telescope (P48)/Zwicky Transient Facility (ZTF) camera in 30 s r-band exposures taken during evening astronomical twilight from 2019 Sep 20 to 2022 March 7 and during morning astronomical twilight sky from 2019 Sep 21 to 2022 Sep 29. More than 21,940 exposures were taken in evening astronomical twilight within 31 degrees and 66 degrees from the Sun with an r-band limiting magnitude between 18.0 and 20.8 (5th to 95th percentile), and more than 24,370 exposures were taken in morning astronomical twilight within 31 degrees and 65 degrees from the Sun with an r-band limiting magnitude between 18.2 and 20.9 (5th to 95th percentile). The morning and evening twilight pointings show a slight seasonal dependence in limiting magnitude and ability to point closer towards the Sun, with limiting magnitude improving by 0.5 magnitudes during the summer months and Sun-centric angular distances as small as 31-32 degrees during the spring and fall months. In total, the one Aylo, (594913) 'Ayl & oacute;'chaxnim, and 4 Atiras, 2020 OV1, 2021 BS1, 2021 PB2, and 2021 VR3, were discovered in evening and morning twilight observations. Additional twilight survey discoveries also include 6 long period comets: C/2020 T2, C/2020 V2, C/2021 D2, C/2021 E3, C/2022 E3 and C/2022 P3, and two short period comets: P/2021 N1 and P/2022 P2 using deep learning comet detection pipelines. The P48/ZTF twilight survey also recovered 11 known Atiras, one Aylo, three short period comes, two long period comets, one interstellar object, 45,536 Main Belt asteroids, and 265 near-Earth objects. Additionally, observations from the GROWTH network of telescopes were used to recover the Aylo, Atira, and comet discoveries made during the ZTF twilight survey. Lastly, we discuss the future twilight surveys for the discovery of Aylos such as with the Vera Rubin Observatory which will have a twilight survey starting in its first year of operations and will cover the sky as within 45 degrees from the Sun. Twilight surveys such as those by ZTF and future surveys will provide opportunities for the discovery of asteroids inside the orbits of the terrestrial planets that would otherwise be unavailable in conventional sky survey observations.

We present optical, radio, and X-ray observations of EP250108a/SN 2025kg, a broad-line Type Ic supernova (SN Ic-BL) accompanying an Einstein Probe (EP) fast X-ray transient at z = 0.176. EP250108a/SN 2025kg possesses a double-peaked optical light curve, and its spectrum transitions from a blue underlying continuum to a typical SN Ic-BL spectrum over time. We fit a radioactive decay model to the second peak of the optical light curve and find SN parameters that are consistent with the SN Ic-BL population, while its X-ray and radio properties are consistent with those of low-luminosity GRB (LLGRB) 060218/SN 2006aj. We explore three scenarios to understand the system’s multiwavelength emission: (a) SN ejecta interacting with an extended circumstellar medium (CSM), (b) the shocked cocoon of a collapsar-driven jet choked in its stellar envelope, and (c) the shocked cocoon of a collapsar-driven jet choked in an extended CSM. Models (b) and (c) can explain the optical light curve and are also consistent with the radio and X-ray observations. We favor model (c) because it can self-consistently explain both the X-ray prompt emission and first optical peak, but we do not rule out model (b). From the properties of the first peak in model (c), we find evidence that EP250108a/SN 2025kg interacts with an extended CSM and infer an envelope mass M _e ∼ 0.1 M _⊙ and radius R _e ∼ 4 × 10 ^13 cm. EP250108a/SN 2025kg’s multiwavelength properties make it a close analog to LLGRB 060218/SN 2006aj and highlight the power of early follow-up observations in mapping the environments of massive stars prior to core collapse.

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 similar to 1.3 x 10(16) cm, moving with a maximum velocity of 4275 km s(-1), and the shell of SN 2022xgc is located at similar to 0.8 x 10(16) cm, reaching up to 4400 km s(-1). These shells were expelled similar to 11 and similar to 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 similar to 7 - 9 M-circle dot. The ejecta masses inferred from the magnetar modeling are not consistent with the PPI scenario pointing toward stars > 50 M-circle dot 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 10(7.8) M-circle dot, a star formation rate of 0.96(-0.26)(+0.47) M-circle dot yr(-1), and a metallicity of similar to 0.2 Z(circle dot).

SN 2023ixf is one of the brightest core collapse supernovae of the 21st century and offers a rare opportunity to investigate the late stage of a supernova through nebular phase spectroscopy. We present four nebular phase spectra from day +291 to +413 after explosion. This is supplemented with high-cadence early phase spectroscopic observations and photometry covering the first 500 d to investigate explosion parameters. The narrow and blueshifted nebular oxygen emission lines are used to infer an ejected oxygen mass of < 0.65 M-circle dot, consistent with models of a relatively low mass (M-ZAMS < 15 M-circle dot) progenitor. An energy of 0.3-1.4 x 10(51) erg and a light curve powered by an initial Ni-56 mass of 0.049 +/- 0.005 M-circle dot appear consistent with a relatively standard Type II explosion, while an incomplete gamma-ray trapping (with time-scale of 240 +/- 4 d) suggests a lower ejecta mass. Assuming a typical explosion, the broad hydrogen and calcium profiles suggest a common origin within a lower mass, partially stripped envelope. Hydrogen emission broadens with time, indicating contribution from an additional power source at an extended distance; while the emergence of high-velocity (similar to 6000 km s(-1)) hydrogen emission features (beginning around day +200) may be explained by shock interaction with a dense hydrogen-rich region located at similar to 1.5 x 10(16) cm. Such envelope mass-loss for a low-mass progenitor may be explained through theoretical models of binary interaction.

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 CSM interaction is the main responsible for their features. However, other mechanisms can not be discarded. Large sample analysis can provide clarification. 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. We analyze the gri light curves of SLSNe II obtained through ZTF. We study peak absolute magnitudes and characteristic timescales. When possible we compute g-r colors, pseudo-bolometric light curves, and estimate lower limits for their total radiated energy. We also study the luminosity distribution of our sample and estimate the percentage of them that would be observable by the LSST. Finally, we compare our sample to other H-rich SNe and to H-poor SLSNe I. 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 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. We present the largest sample of SLSNe II light curves to date, comprising of 107 events. Their diversity could be explained by considering 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 multi-wavelength light curve coverage becomes paramount. LSST could miss 30 percent of the ZTF events in the its footprint in gri bands. Redder bands become important to construct complete samples.

Stars are initially powered by the fusion of hydrogen to helium. These ashes serve as fuel in a series of stages1-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-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-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.