Search

The second data release of Type Ia supernovae (SNe Ia) observed by the Zwicky Transient Facility has provided a homogeneous sample of 3628 SNe Ia with photometric and spectral information. This unprecedented sample size enables us to better explore our currently tentative understanding of the dependence of host environment on SN Ia properties. In this paper, we make use of two-dimensional image decomposition to model the host galaxies of SNe Ia. We model elliptical galaxies as well as disk/spiral galaxies with or without central bulges and bars. This allows for the categorisation of SN Ia based on their morphological host environment, as well as the extraction of intrinsic galaxy properties corrected for both cosmological and atmospheric effects. We find that although this image decomposition technique leads to a significant bias towards elliptical galaxies in our final sample of galaxies, the overall results are robust. By successfully modelling 728 host galaxies, we find that the photometric properties of SNe Ia found in disks and in elliptical galaxies, correlate fundamentally differently with their host environment. We identified strong linear relations between light-curve stretch and our model-derived galaxy colour for both the elliptical (16.8$\sigma$) and disk (5.1$\sigma$) subpopulations of SNe Ia. Lower stretch SNe Ia are found in redder environments, which we identify as an age/metallicity effect. Within the subpopulation of SNe Ia found in disk containing galaxies, we find a significant linear trend (6.1$\sigma$) between light-curve stretch and model-derived local $r$-band surface brightness, which we link to the age/metallicity gradients found in disk galaxies. SN Ia colour shows little correlation with host environment as seen in the literature. We identify a possible dust effect in our model-derived surface brightness (3.3$\sigma$), for SNe Ia in disk galaxies., Comment: ZTF SN Ia DR2 release paper. Submitted to A&A (ZTF DR2 Special Issue)

Massive stars experience strong mass-loss, producing a dense, H-rich circumstellar medium (CSM). After the explosion, the collision and continued interaction of the supernova (SN) ejecta with the CSM power the light curve through the conversion of kinetic energy into radiation. When the interaction is strong, the light curve shows a broad peak and high luminosity lasting for a relatively long time. Also the spectral evolution is slower, compared to non-interacting SNe. Energetic shocks between the ejecta and the CSM create the ideal conditions for particle acceleration and production of high-energy (HE) neutrinos above 1 TeV. In this paper, we study four strongly-interacting Type IIn SNe: 2021acya, 2021adxl, 2022qml, and 2022wed to highlight their peculiar characteristics, derive the kinetic energy of the explosion and the characteristics of the CSM, infer clues on the possible progenitors and their environment and relate them to the production of HE neutrinos. The SNe analysed in this sample exploded in dwarf, star-forming galaxies and they are consistent with energetic explosions and strong interaction with the surrounding CSM. For SNe 2021acya and 2022wed, we find high CSM masses and mass-loss rates, linking them to very massive progenitors. For SN 2021adxl, the spectral analysis and less extreme CSM mass suggest a stripped-envelope massive star as possible progenitor. SN 2022qml is marginally consistent with being a Type Ia thermonuclear explosion embedded in a dense CSM. The mass-loss rates for all SNe are consistent with the expulsion of several solar masses of material during eruptive episodes in the last few decades before the explosion. Finally, we find that the SNe in our sample are marginally consistent with HE neutrino production., Comment: Submitted to A&A

We present photometric and spectroscopic observations of SN2020xga and SN2022xgc, two hydrogen-poor superluminous supernovae (SLSNe-I) at $z = 0.4296$ and $z = 0.3103$ respectively, that show an additional set of broad Mg II absorption lines, blueshifted by a few thousand km s$^{-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 CSM expelled shortly before the explosion. The peak rest-frame $g$-band magnitude of SN2020xga is $-22.30 \pm 0.04$ mag and of SN2022xgc is $-21.97 \pm 0.05$ mag, placing them among the brightest SLSNe-I. We use 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 SN2020xga resides at $\sim 1.3 \times 10^{16} \rm cm$ moving with a maximum velocity of $4275~\rm km~s^{-1}$, and the shell of SN2022xgc is located at $\sim 0.8 \times 10^{16} \rm cm$ reaching up to $4400~\rm km~s^{-1}$. These shells were expelled $\sim 11$ and $\sim 5$ months before explosion for SN2020xga and SN2022xgc respectively, possibly as a result of Luminous Blue Variable-like eruptions or pulsational pair instability (PPI) mass loss. We also analyze optical photometric data and model 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 $\sim 7-9 \rm~M_\odot$. The ejecta masses inferred from the magnetar modeling are not consistent with the PPI scenario pointing towards stars $> 50~\rm M_\odot$ He-core, hence alternative scenarios such as fallback accretion are discussed., Comment: 25 pages text, 8 pages appendix, 21 figures. Submitted to A&A

SN 2023tsz is a Type Ibn supernova (SNe Ibn) discovered in an extremely low-mass host. SNe Ibn are an uncommon subtype of stripped-envelope core-collapse SNe. They are characterised by narrow helium emission lines in their spectra and are believed to originate from the collapse of massive Wolf-Rayet (WR) stars, though their progenitor systems still remain poorly understood. In terms of energetics and spectrophotometric evolution, SN 2023tsz is largely a typical example of the class, although line profile asymmetries in the nebular phase are seen, which may indicate the presence of dust formation or unshocked circumstellar material. Intriguingly, SN 2023tsz is located in an extraordinarily low-mass host galaxy that is in the 2nd percentile for SESN host masses and star formation rates (SFR). The host has a radius of 1.0 kpc, a $g$-band absolute magnitude of $-12.73$, and an estimated metallicity of $\log(Z_{*}/Z_{\odot}$) = $-1.56$. The SFR and metallicity of the host galaxy raise questions about the progenitor of SN 2023tsz. The low SFR suggests that a star with sufficient mass to evolve into a WR would be uncommon in this galaxy. Further, the very low-metallicity is a challenge for single stellar evolution to enable H and He stripping of the progenitor and produce a SN Ibn explosion. The host galaxy of SN 2023tsz adds another piece to the ongoing puzzle of SNe Ibn progenitors, and demonstrates that they can occur in hosts too faint to be observed in contemporary sky surveys at a more typical SN Ibn redshift., Comment: 14 pages, 10 figures, submitted to MNRAS

Type Ia supernovae (SNe Ia) constitute an historical probe to derive cosmological parameters through the fit of the Hubble-Lema\^itre diagram, i.e. SN Ia distance modulus versus their redshift. In the era of precision cosmology, realistic simulation of SNe Ia for any survey entering in an Hubble-Lema\^itre diagram is a key tool to address observational systematics, like Malmquist bias. As the distance modulus of SNe Ia is derived from the fit of their light-curves, a robust simulation framework is required. In this paper, we present the performances of the simulation framework skysurvey to reproduce the the Zwicky Transient Facility (ZTF) SN Ia DR2 covering the first phase of ZTF running from April 2018 up to December 2020. The ZTF SN Ia DR2 sample correspond to almost 3000 classified SNe Ia of cosmological quality. First, a targeted simulation of the ZTF SN Ia DR2 was carried on to check the validity of the framework after some fine tuning of the observing conditions and instrument performance. Then, a realistic simulation has been run using observing ZTF logs and ZTF SN Ia DR2 selection criteria on simulated light-curves to demonstrate the ability of the simulation framework to match the ZTF SN Ia DR2 sample. Furthermore a redshift dependency of SALT2 light-curve parameters (stretch and colour) was conducted to deduce a volume limited sample, i.e. an unbiased SNe Ia sample, characterized with $z_{lim} \leq 0.06$. This volume limited sample of about 1000 SNe Ia is unique to carry on new analysis on standardization procedure with a precision never reached (those analysis are presented in companion papers)., Comment: 10 pages, 10 figures. Submitted to A & A

The Zwicky Transient Facility SN Ia Data Release 2 (ZTF SN Ia DR2) contains more than 3,000 Type Ia supernovae (SNe Ia), providing the largest homogeneous low-redshift sample of SNe Ia. Having at least one spectrum per event, this data collection is ideal for large-scale statistical studies of the photometric, spectroscopic and host-galaxy properties of SNe Ia, particularly of the more rare "peculiar" subclasses. In this paper, we first present the method we developed to spectroscopically classify the SNe in the sample, and the techniques we used to model their multi-band light curves and explore their photometric properties. We then show a method to distinguish between the "peculiar" subtypes and the normal SNe Ia. We also explore the properties of their host galaxies and estimate their relative rates, focusing on the "peculiar" subtypes and their connection to the cosmologically useful SNe Ia. Finally, we discuss the implications of our study with respect to the progenitor systems of the "peculiar" SN Ia events., Comment: 19 pages, 13 figures, submitted to Astronomy and Astrophysics

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., Comment: 35 pages. 27 figures. 9 tables

We present optical and near-infrared observations of two Type Ibn supernovae (SNe), SN 2018jmt and SN 2019cj. Their light curves have rise times of about 10 days, reaching an absolute peak magnitude of $M_g$(SN 2018jmt) = $-$19.07 $\pm$ 0.37 and $M_V$(SN 2019cj) = $-$18.94 $\pm$ 0.19 mag, respectively. The early-time spectra of SN 2018jmt are dominated by a blue continuum, accompanied by narrow (600$-$1000 km~s$^{-1}$) He I lines with P-Cygni profile. At later epochs, the spectra become more similar to those of the prototypical SN Ibn 2006jc. At early phases, the spectra of SN 2019cj show flash ionisation emission lines of C III, N III and He II superposed on a blue continuum. These features disappear after a few days, and then the spectra of SN 2019cj evolve similarly to those of SN 2018jmt. The spectra indicate that the two SNe exploded within a He-rich circumstellar medium (CSM) lost by the progenitors a short time before the explosion. We model the light curves of the two SNe Ibn to constrain the progenitor and the explosion parameters. The ejecta masses are consistent with either that expected for a canonical SN Ib ($\sim$ 2 M$_{\odot}$) or those from a massive WR star ($>$ $\sim$ 4 M$_{\odot}$), with the kinetic energy on the order of $10^{51}$ erg. The lower limit on the ejecta mass ($>$ $\sim$ 2 M$_{\odot}$) argues against a scenario involving a relatively low-mass progenitor (e.g., $M_{ZAMS}$ $\sim$ 10 M$_{\odot}$). We set a conservative upper limit of $\sim$0.1 M$_{\odot}$ for the $^{56}$Ni masses in both SNe. From the light curve modelling, we determine a two-zone CSM distribution, with an inner, flat CSM component, and an outer CSM with a steeper density profile. The physical properties of SN 2018jmt and SN 2019cj are consistent with those expected from the core collapse of relatively massive, stripped-envelope (SE) stars., Comment: 28 pages, 19 figures, Accepted for publication in Astronomy & Astrophysics

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. 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, then we discuss more in detail the high resolution spectrum collected for NGC 300 2008OT-1 with the Very Large Telescope equipped with UVES. Finally we analyse late time spectra of NGC 300 2008OT-1 and AT 2019ahd through comparisons with both synthetic and observed spectra. 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$\alpha$, H$\beta$ 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 ($\sim$30 km s$^{-1}$) progenitor wind. At late epochs, optical spectra of NGC 300 2008OT-1 and AT 2019ahd show broad ($\sim$2500 km s$^{-1}$) emission features at $\sim$6170 A and $\sim$7000 A 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., Comment: 28 pages, 19 figures. Submitted to A&A

More than 3000 spectroscopically confirmed Type Ia supernovae (SNe Ia) are presented in the Zwicky Transient Facility SN Ia Data Release 2 (ZTF DR2). In this paper, we detail the spectral properties of 482 SNe Ia near maximum light, up to a redshift limit of $z$ $\leq$ 0.06. We measure the velocities and pseudo-equivalent widths (pEW) of key spectral features (Si II $\lambda$5972 and Si II $\lambda$6355) and investigate the relation between the properties of the spectral features and the photometric properties from the SALT2 light-curve parameters as a function of spectroscopic sub-class. We discuss the non-negligible impact of host galaxy contamination on SN Ia spectral classifications, as well as investigate the accuracy of spectral template matching of the ZTF DR2 sample. We define a new subclass of underluminous SNe Ia (`04gs-like') that lie spectroscopically between normal SNe Ia and transitional 86G-like SNe Ia (stronger Si II $\lambda$5972 than normal SNe Ia but significantly weaker Ti II features than `86G-like' SNe). We model these `04gs-like' SN Ia spectra using the radiative-transfer spectral synthesis code tardis and show that cooler temperatures alone are unable to explain their spectra; some changes in elemental abundances are also required. However, the broad continuity in spectral properties seen from bright (`91T-like') to faint normal SN Ia, including the transitional and 91bg-like SNe Ia, suggests that variations within a single explosion model may be able to explain their behaviour.

Type Ia supernova (SN Ia) light curves have a secondary maximum that exists in the $r$, $i$, and near-infrared filters. The secondary maximum is relatively weak in the $r$ band, but holds the advantage that it is accessible, even at high redshift. We used Gaussian Process fitting to parameterise the light curves of 893 SNe Ia from the Zwicky Transient Facility's (ZTF) second data release (DR2), and we were able to extract information about the timing and strength of the secondary maximum. We found $>5\sigma$ correlations between the light curve decline rate ($\Delta m_{15}(g)$) and the timing and strength of the secondary maximum in the $r$ band. Whilst the timing of the secondary maximum in the $i$ band also correlates with $\Delta m_{15}(g)$, the strength of the secondary maximum in the $i$ band shows significant scatter as a function of $\Delta m_{15}(g)$. We found that the transparency timescales of 97 per cent of our sample are consistent with double detonation models, and that SNe Ia with small transparency timescales ($<$ 32 d) reside predominantly in locally red environments. We measured the total ejected mass for the normal SNe Ia in our sample using two methods, and both were consistent with medians of $1.3\ \pm \ 0.3$ and $1.2\ \pm\ 0.2$ solar masses. We find that the strength of the secondary maximum is a better standardisation parameter than the SALT light curve stretch ($x_1$). Finally, we identified a spectral feature in the $r$ band as Fe II, which strengthens during the onset of the secondary maximum. The same feature begins to strengthen at $<$ 3 d post maximum light in 91bg-like SNe. Finally, the correlation between $x_1$ and the strength of the secondary maximum was best fit with a broken line, with a split at $x_1^0\ =\ -0.5\ \pm\ 0.2$, suggestive of the existence of two populations of SNe Ia., Comment: 18 pages, 14 figures

The unprecedented statistics of detected Type Ia supernovae (SNe Ia) brought by the Zwicky Transient Facility enables us to probe the impact of the Large-Scale Structure on the properties of these objects. The goal of this paper is to explore the possible impact of the under-dense part of the large-scale structure on the intrinsic SALT2 light curve properties of SNe Ia and uncover possible biases in SN Ia analyses. With a volume-limited selection of ZTF-Cosmo-DR2 Type Ia supernovae overlapping with the SDSS-DR7 survey footprint, we investigate the distribution of their properties with regard to voids detected in the SDSS-DR7 galaxy sample. We further use Voronoi volumes as proxy for local density environments within the large-scale structure. We find a moderate dependency of the stretch toward the localisation around the void centre and none when considering colour. The local Voronoi volumes mostly affect the fraction of low/high stretch supernovae. With the current statistics available, we consider that the impact of high or low local density environment can be considered as a proxy for the colour of the host galaxy. Under-dense environments should not cause any biases in supernova analyses., Comment: 10 pages, 8 figures. Submitted to A&A

Type Ia supernova (SNIa) are excellent probes of local distance, and the increasing sample sizes of SNIa have driven an increased need to study the associated systematic uncertainties and improve the standardisation methods in preparation for the next generation of cosmological surveys into the dark energy equation-of-state $w$. We aim to probe the potential change in the SNIa standardisation parameter $c$ with redshift and the host-galaxy of the supernova. Improving the standardisation of SNIa brightnesses will require accounting for the relationship between the host and the SNIa, and potential shifts in the SNIa standardisation parameters with redshift will cause biases in the recovered cosmology. Here, we assemble a volume-limited sample of ~3000 likely SNIa across a redshift range of $z = 0.015$ to $z = 0.36$. This sample is fitted with changing mass and redshift bins to determine the relationship between intrinsic properties of SNe Ia and their redshift and host galaxy parameters. We then investigate the colour-luminosity parameter $\beta$ as a further test of the SNIa standardisation process. We find that the changing colour distribution of SNe Ia with redshift is driven by dust at a confidence of $>4\sigma$. Additionally, we show a strong correlation between the host galaxy mass and the colour-luminosity coefficient $\beta$ ($> 4\sigma$), even when accounting for the quantity of dust in a host galaxy.

Type Ia supernova (SN Ia) cosmology relies on the estimation of lightcurve parameters to derive precision distances that leads to the estimation of cosmological parameters. The empirical SALT2 lightcurve modeling that relies on only two parameters, a stretch x1, and a color c, has been used by the community for almost two decades. In this paper we study the ability of the SALT2 model to fit the nearly 3000 cosmology-grade SN Ia lightcurves from the second release of the Zwicky Transient Facility (ZTF) cosmology science working group. While the ZTF data was not used to train SALT2, the algorithm is modeling the ZTF SN Ia optical lightcurves remarkably well, except for lightcurve points prior to -10 d from maximum, where the training critically lacks statistics. We find that the lightcurve fitting is robust against the considered choice of phase-range, but we show the [-10; +40] d range to be optimal in terms of statistics and accuracy. We do not detect any significant features in the lightcurve fit residuals that could be connected to the host environment. Potential systematic population differences related to the SN Ia host properties might thus not be accountable for by the addition of extra lightcurve parameters. However, a small but significant inconsistency between residuals of blue- and red-SN Ia strongly suggests the existence of a phase-dependent color term, with potential implications for the use of SNe Ia in precision cosmology. We thus encourage modellers to explore this avenue and we emphasize the importance that SN Ia cosmology must include a SALT2 retraining to accurately model the lightcurves and avoid biasing the derivation of cosmological parameters., Comment: 10 pages, 9 figures. Submitted to Astronomy and Astrophysics

As Type Ia supernova cosmology transitions from a statistics dominated to a systematics dominated era, it is crucial to understand leftover unexplained uncertainties affecting their luminosity, such as the ones stemming from astrophysical biases. Indeed, SNe Ia are standardisable candles, whose absolute magnitude reach a 0.15~mag scatter once empirical correlations with their lightcurve stretch and colour and with their environment are accounted for. In this paper, we investigate how the standardisation process of SNe Ia depends on environment, to ultimately reduce their scatter in magnitude, focusing on colour standardisation. We use the volume-limited ZTF SN Ia DR2 sample, which offers unprecedented statistics for the low redshift ($z<0.06$) range. We first study the colour distribution, focusing on the effects of dust, to then select a dustless subsample of objects from low stellar mass environments and from the outskirts of their host galaxies. We then look at the colour-residuals relation and its associated parameter $\beta$. Finally, we investigate the colour dependency of the environment-dependent magnitude offsets (steps), to try to disentangle intrinsic and extrinsic colour origin. Our sample probes well the red tail of the colour distribution, up to $c=0.8$. The dustless sample exhibits a significantly lower red tail ($4.6\sigma$) in comparison to the whole sample. This suggests that reddening above $c\geq0.2$ is dominated by host interstellar dust absorption. Looking at the colour-residuals relation, we find it to be linear with lightcurve colour. We show hints of a potential evolution of $\beta$ with host stellar mass at a $2.5\sigma$ level. Finally, unlike recent claims from the literature, we see no evolution of steps as a function of lightcurve colour, suggesting that dust may not be the dominating mechanism responsible for the environmental dependency of SNe Ia magnitude., Comment: 10 pages, 7 figures, submitted to Astronomy and Astrophysics

Understanding Type Ia supernovae (SNe~Ia) and the empirical standardisation relations that make them excellent distance indicators is vital to improving cosmological constraints. SN~Ia ``siblings", i.e. two or more SNe~Ia in the same host or parent galaxy offer a unique way to infer the standardisation relations and their diversity across the population. We analyse a sample of 25 SN~Ia pairs, observed homogeneously by the Zwicky Transient Factory (ZTF) to infer the SNe~Ia light curve width-luminosity and colour-luminosity parameters $\alpha$ and $\beta$. Using the pairwise constraints from siblings, allowing for a diversity in the standardisation relations, we find $\alpha = 0.218 \pm 0.055 $ and $\beta = 3.084 \pm 0.312$, respectively, with a dispersion in $\alpha$ and $\beta$ of $\leq 0.195$ and $\leq 0.923$, respectively, at 95$\%$ C.L. While the median dispersion is large, the values within $\sim 1 \sigma$ are consistent with no dispersion. Hence, fitting for a single global standardisation relation, we find $\alpha = 0.228 \pm 0.029 $ and $\beta = 3.160 \pm 0.191$. We find a very small intrinsic scatter of the siblings sample $\sigma_{\rm int} \leq 0.10$ at 95\% C.L. compared to $\sigma_{\rm int} = 0.22 \pm 0.04$ when computing the scatter using the Hubble residuals without comparing them as siblings. Splitting the sample based on host galaxy stellar mass, we find that SNe~Ia in both subsamples have consistent $\alpha$ and $\beta$. The $\beta$ value is consistent with the value for the cosmological sample. However, we find a higher $\alpha$ by $\sim 2.5 - 3.5 \sigma$. The high $\alpha$ is driven by low $x_1$ pairs, potentially suggesting that the slow and fast declining SN~Ia have different slopes of the width-luminosity relation. We can confirm or refute this with increased statistics from near future time-domain surveys. (abridged), Comment: 12 pages, 7 figures, submitted to A&A

Understanding the impact of the astrophysical environment on Type Ia supernova (SN Ia) properties is crucial to minimize systematic uncertainties in cosmological analyses based on this probe. We investigate the dependence of the SN Ia SALT2.4 light-curve stretch on the distance from their nearest galaxy cluster to study a potential effect of the intracluster medium (ICM) environment on SN Ia intrinsic properties. We use the largest SN Ia sample to date and cross-match it with existing X-ray, Sunyaev-Zel'dovich, and optical cluster catalogs in order to study the dependence between stretch and distance to the nearest detected cluster from each SN Ia. We model the underlying stretch distribution with a Gaussian mixture with relative amplitudes that depend on redshift and cluster-centric distance. We find a significant improvement of the fit quality of the stretch distribution if we include the distance-dependant term in the model with a variation of the Akaike information criterion $\rm{\Delta AIC} = -10.2$. Because of the known correlation between galaxy age and distance from cluster center, this supports previous evidence that the age of the stellar population is the underlying driver of the bimodial shape of the SN Ia stretch distribution. We further compute the evolution of the fraction of quenched galaxies as a function of distance with respect to cluster center from our best-fit model of the SNe Ia stretch distribution and compare it to previous results obtained from $H\alpha$ line measurements, optical broadband photometry, and simulations. We find our estimate to be compatible with these results. The results of this work indicate that SNe Ia searches at high redshift targeted towards clusters to maximize detection probability should be considered with caution as the stretch distribution of the detected sample would be strongly biased towards the old sub-population of SNe Ia., Comment: Accepted in Astronomy & Astrophysics

To get distances, Type Ia Supernovae magnitudes are corrected for their correlation with lightcurve width and colour. Here we investigate how this standardisation is affected by the SN environment, with the aim to reduce scatter and improve standardisation. We first study the SN Ia stretch distribution, as well as its dependence on environment, as characterised by local and global (g-z) colour and stellar mass. We then look at the standardisation parameter $\alpha$, which accounts for the correlation between residuals and stretch, along with its environment dependence and linearity. We finally compute magnitude offsets between SNe in different astrophysical environments after colour and stretch standardisation, aka steps. This analysis is made possible due to the unprecedented statistics of the ZTF SN Ia DR2 volume-limited sample. The stretch distribution exhibits a bimodal behaviour, as previously found in literature. However, we find the distribution means to decrease with host stellar mass at a 9.0$\sigma$ significance. We demonstrate, at the 14.3$\sigma$ level, that the stretch-magnitude relation is non-linear, challenging the usual linear stretch-residuals relation. Fitting for a broken-$\alpha$ model, we indeed find two different slopes between stretch regimes ($x_1<-0.49\pm0.06$): $\alpha_{low}=0.28\pm0.01$ and $\alpha_{high}=0.09\pm0.01$, a $\Delta_{\alpha}=-0.19\pm0.01$ difference. As the relative proportion of SNe Ia in the high-/low-stretch modes evolves with redshift and environment, this implies that a linear $\alpha$ also evolves with redshift and environment. Concerning the environmental magnitude offset $\gamma$, we find it to be greater than 0.14 mag regardless of the considered environmental tracer used (local or global colour and stellar mass), all measured at the $\geq 6\sigma$ level, increased to $\sim0.18\pm0.01$ mag when accounting for the stretch-non linearity., Comment: 15 pages, 9 figures, submitted to Astronomy and Astrophysics

SNe Ia are used to determine the distance-redshift relation and build the Hubble diagram. Neglecting their host-galaxy peculiar velocities (PVs) may bias the measurement of cosmological parameters. The smaller the redshift, the larger the effect is. We use realistic simulations of SNe Ia observed by the Zwicky Transient Facility (ZTF) to investigate the effect of different methods to take into account PVs. We study the impact of neglecting galaxy PVs and their correlations in an analysis of the SNe Ia Hubble diagram. We find that it is necessary to use the PV full covariance matrix computed from the velocity power spectrum to take into account the sample variance. Considering the results we have obtained using simulations, we determine the PV systematic effects in the context of the ZTF DR2 SNe Ia sample. We determine the PV impact on the intercept of the Hubble diagram, $a_B$, which is directly linked to the measurement of $H_0$. We show that not taking into account PVs and their correlations results in a shift of the $H_0$ value of about $1.0$km.s$^{-1}$.Mpc$^{-1}$ and a slight underestimation of the $H_0$ error bar., Comment: 12 pages, 4 figures

The Zwicky Transient Facility (ZTF) was expected to detect more than one strong gravitationally-lensed supernova (glSN) per year, but only one event was identified in the first four years of the survey. This work investigates selection biases in the search strategy that could explain the discrepancy and revise discovery predictions. We present simulations of realistic lightcurves for lensed thermonuclear (glSNIa) and core-collapse supernova (glCCSN) explosions over a span of 5.33 years of the survey, utilizing the actual observation logs of ZTF. We find that the magnitude limit in spectroscopic screening significantly biases the selection towards highly magnified glSNe, for which the detection rates are consistent with the identification of a single object by ZTF. To reach the higher predicted rate of detections requires an optimization of the identification criteria for fainter objects. We find that around 1.36 (3.08) Type Ia SNe (CCSNe) are identifiable with the magnification method per year in ZTF, but when applying the magnitude cut of m < 19 mag, the detection rates decrease to 0.17 (0.32) per year. We compare our simulations with the previously found lensed Type Ia SNe, iPTF16geu and SN Zwicky, and conclude that considering the bias towards highly magnified events, the findings are within expectations in terms of detection rates and lensing properties of the systems. In addition, we provide a set of selection cuts based on simple observables to distinguish glSNe from regular, unlensed, supernovae to select potential candidates for spectroscopic and high-spatial resolution follow-up campaigns. We find optimal cuts in observed colours $g-r$, $g-i$, and $r-i$ as well as in the colour SALT2 fit parameter. The developed pipeline and the simulated lightcurves employed in this analysis can be found in the $LENSIT$ github repository., Comment: 12 pages, 14 figures

The supernova remnant SNR 0540-69.3, twin of the Crab Nebula, offers an excellent opportunity to study the continuum emission from a young pulsar and pulsar-wind nebula (PWN). We present observations taken with the VLT instruments MUSE and X-shooter in the wavelength range 3000-25,000 \r{A}, which allow us to study spatial variations of the optical spectra, along with the first near-infrared (NIR) spectrum of the source. We model the optical spectra with a power law (PL) $F_\nu\propto\nu^{-\alpha}$ and find clear spatial variations (including a torus-jet structure) in the spectral index across the PWN. Generally, we find spectral hardening toward the outer parts, from $\alpha\sim1.1$ to $\sim0.1$, which may indicate particle reacceleration by the PWN shock at the inner edge of the ejecta or alternatively time variability of the pulsar wind. The optical-NIR spectrum of the PWN is best described by a broken PL, confirming that several breaks are needed to model the full spectral energy distribution of the PWN, suggesting the presence of more than one particle population. Finally, subtracting the PWN contribution from the pulsar spectrum we find that the spectrum is best described with a broken-PL model with a flat and a positive spectral index, in contrast to the Crab pulsar that has a negative spectral index and no break in the optical. This might imply that pulsar differences propagate to the PWN spectra., Comment: 29 pages, 19 Figures, 1 Table. Accepted for publication in ApJ

Obtaining spectroscopic observations of the progenitors of core-collapse supernovae is often unfeasible due to an inherent lack of knowledge as to which stars will go supernova and when they will explode. In this letter, we present photometric and spectroscopic observations of the progenitor activity of SN 2023fyq in the preceding 150 days before the He-rich progenitor exploded as a Type Ibn supernova. The progenitor of SN 2023fyq shows an exponential rise in flux prior to core-collapse. Complex He I emission line features are observed, with a P-Cygni like profile, as well as an evolving broad base with velocities on the order of 10,000 km/s, possibly due to electron scattering. The luminosity and evolution of SN 2023fyq are consistent with a faint Type Ibn, reaching a peak r-band magnitude of 18.1 mag, although there is some uncertainty in the distance to the host, NGC 4388, located in the Virgo cluster. We present additional evidence of asymmetric He-rich material being present prior to the explosion of SN 2023fyq, as well as after, suggesting this material has survived the ejecta-CSM interaction. Broad [O I] and the Ca II triplet lines are observed at late phases, confirming that SN 2023fyq was a genuine supernova rather than a non-terminal interacting transient. SN 2023fyq provides insight into the final moments of a massive star's life, highlighting that the progenitor is likely highly unstable before core-collapse., Comment: 7 Pages, 5 Figures, accepted to A&A Letters

We present an extensive analysis of the optical and UV properties of AT2023clx, the closest TDE to date, that occurred in the nucleus of the interacting LINER galaxy, NGC3799 (z=0.01107). After correcting for the host reddening (E(B-V) = 0.179 mag), we find its peak absolute g-band magnitude to be -18.03{+/-}0.07 mag, and its peak bolometric luminosity to be L=(1.57{+/-}0.19)x10^43 erg/s. AT2023clx displays several distinctive features: first, it rose to peak within 10.4{+/-}2.5 days, making it the fastest rising TDE to date. Our SMBH mass estimate of M_BH ~ 10^6 Msol rules out the possibility of an intermediate-mass BH as the reason for the fast rise. Dense spectral follow-up reveals a blue continuum that cools slowly and broad Balmer and HeII lines as well as weak HeI 5876,6678 emission features that are typically seen in TDEs. The early, broad (width ~ 15000 km/s) profile of Ha matches theoretical expectations from an optically thick outflow. A flat Balmer decrement (~ 1.58) suggests that the lines are collisionally excited rather than being produced via photoionisation, in contrast to typical active galactic nuclei. A second distinctive feature, seen for the first time in TDE spectra, is a sharp, narrow emission peak at a rest wavelength of ~ 6353 A. This feature is clearly visible up to 10d post-peak; we attribute it to clumpy material preceding the bulk outflow, which manifests as a high-velocity component of Ha (-9584 km/s). Its third distinctive feature is the rapid cooling during the first ~ 20 days after peak, reflected as a break in the temperature evolution. Combining these findings, we propose a scenario for AT2023clx involving the disruption of a very low-mass star (<=0.1 Msol) with an outflow launched in our line of sight and with disruption properties that led to efficient circularisation and prompt accretion disc formation, observed through a low-density photosphere., Comment: Accepted for publication in Astronomy and Astrophysics (in production; official acceptance date: 28/06/2024)

SN 2021adxl is a slowly evolving, luminous, Type IIn supernova with asymmetric emission line profiles, similar to the well-studied SN 2010jl. We present extensive optical, near-ultraviolet, and near-infrared photometry and spectroscopy covering ~1.5 years post discovery. SN 2021adxl occurred in an unusual environment, atop a vigorously star-forming region that is offset from its host galaxy core. The appearance of Ly-alpha, O II, as well as the compact core, would classify the host of SN 2021adxl as a Blueberry galaxy, analogous to the higher redshift Green Pea galaxies. Using several abundance indicators, we find a metallicity of the explosion environment of only 10% solar, the lowest reported metallicity for a Type IIn SN environment. SN 2021adxl reaches a peak magnitude of r ~ -20.2 mag and since discovery, SN 2021adxl has faded by only ~4 magnitudes in the r band with a cumulative radiated energy of ~1.5e50 erg over 18 months. SN 2021adxl shows strong signs of interaction with a complex circumstellar medium, seen by the detection of X-rays, revealed by the detection of coronal emission lines, and through multi-component hydrogen and helium profiles. In order to further understand this interaction, we model the H{\alpha} profile using a Monte-Carlo electron scattering code. The blueshifted high-velocity component is consistent with emission from a radially thin, spherical shell resulting in the broad emission components due to electron scattering. Using the velocity evolution of this emitting shell, we find that the SN ejecta collide with circumstellar material of at least 5 Msun, assuming a steady-state mass-loss rate of 4-6e-3 Msun per year for the first ~200 days of evolution. Continuing the observations of SN 2021adxl may reveal signatures of dust formation or an infrared excess, similar to that seen for SN 2010jl., Comment: 21 pages with 27 Figures. Accepted for publication to A&A. MC line scattering code can be found at https://github.com/Astro-Sean/escatter

Observing a supernova explosion shortly after it occurs can reveal important information about the physics of stellar explosions and the nature of the progenitor stars of supernovae (SNe). When a star with a well-defined edge explodes in vacuum, the first photons to escape from its surface appear as a brief shock-breakout flare. The duration of this flare can extend to at most a few hours even for nonspherical breakouts from supergiant stars, after which the explosion ejecta should expand and cool. Alternatively, for stars exploding within a distribution of sufficiently dense optically thick circumstellar material, the first photons escape from the material beyond the stellar edge, and the duration of the initial flare can extend to several days, during which the escaping emission indicates photospheric heating. The difficulty in detecting SN explosions promptly after the event has so far limited data regarding supergiant stellar explosions mostly to serendipitous observations that, owing to the lack of ultraviolet (UV) data, were unable to determine whether the early emission is heating or cooling, and hence the nature of the early explosion event. Here, we report observations of SN 2023ixf in the nearby galaxy M101, covering the early days of the event. Using UV spectroscopy from the Hubble Space Telescope (HST) as well as a comprehensive set of additional multiwavelength observations, we trace the photometric and spectroscopic evolution of the event and are able to temporally resolve the emergence and evolution of the SN emission.

SN\,2020zbf is a hydrogen-poor superluminous supernova (SLSN) at $z = 0.1947$ that shows conspicuous \ion{C}{II} features at early times, in contrast to the majority of H-poor SLSNe. Its peak magnitude is $M_{\rm g}$ = $-21.2$~mag and its rise time ($\lesssim 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 \ion{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 $^{56}$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_{\rm ej}$ = 1.5~$\rm M_\odot$. 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 \ion{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~\AA,\ 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~$\rm M_\odot$. Modeling the spectral energy distribution of the host galaxy reveals a host mass of 10$^{8.7}$~$\rm M_\odot$, a star formation rate of 0.24$^{+0.41}_{-0.12}$~$\rm M_\odot$~yr$^{-1}$, and a metallicity of $\sim$ 0.4~$\rm Z_\odot$., Comment: Accepted in A&A. 26 pages, 22 figures

Hydrogen-rich Type II supernovae (SNe II) are the most frequently observed class of core-collapse SNe (CCSNe). However, most studies that analyse large samples of SNe II lack events with absolute peak magnitudes brighter than -18.5 mag at rest-frame optical wavelengths. Thanks to modern surveys, the detected number of such luminous SNe II (LSNe II) is growing. There exist several mechanisms that could produce luminous SNe II. The most popular propose either the presence of a central engine (a magnetar gradually spinning down or a black hole accreting fallback material) or the interaction of supernova ejecta with circumstellar material (CSM) that turns kinetic energy into radiation energy. In this work, we study the light curves and spectral series of a small sample of six LSNe II that show peculiarities in their H$\alpha$ profile, to attempt to understand the underlying powering mechanism. We favour an interaction scenario with CSM that is not dense enough to be optically thick to electron scattering on large scales -- thus, no narrow emission lines are observed. This conclusion is based on the observed light curve (higher luminosity, fast decline, blue colours) and spectral features (lack of persistent narrow lines, broad H$\alpha$ emission, lack of H$\alpha$ absorption, weak or nonexistent metal lines) together with comparison to other luminous events available in the literature. We add to the growing evidence that transients powered by ejecta-CSM interaction do not necessarily display persistent narrow emission lines., Comment: 27 pages, 16 figures, accepted for publication in MNRAS

We report on our study of supernova (SN) 2022xxf based on observations obtained during the first four months of its evolution. The light curves (LCs) display two humps of similar maximum brightness separated by 75 days, unprecedented for a broad-lined (BL) Type Ic supernova (SN IcBL). SN 2022xxf is the most nearby SN IcBL to date (in NGC 3705, $z = 0.0037$, at a distance of about 20 Mpc). Optical and near-infrared photometry and spectroscopy are used to identify the energy source powering the LC. Nearly 50 epochs of high signal-to-noise-ratio spectroscopy were obtained within 130 days, comprising an unparalleled dataset for a SN IcBL, and one of the best-sampled SN datasets to date. The global spectral appearance and evolution of SN 2022xxf points to typical SN Ic/IcBL, with broad features (up to $\sim14000$ km s$^{-1}$) and a gradual transition from the photospheric to the nebular phase. However, narrow emission lines (corresponding to $\sim1000-2500$ km s$^{-1}$) are present in the spectra from the time of the second rise, suggesting slower-moving circumstellar material (CSM). These lines are subtle, in comparison to the typical strong narrow lines of CSM-interacting SNe, for example, Type IIn, Ibn, and Icn, but some are readily noticeable at late times such as in Mg I $\lambda$5170 and [O I] $\lambda$5577. Unusually, the near-infrared spectra show narrow line peaks in a number of features formed by ions of O and Mg. We infer the presence of CSM that is free of H and He. We propose that the radiative energy from the ejecta-CSM interaction is a plausible explanation for the second LC hump. This interaction scenario is supported by the color evolution, which progresses to the blue as the light curve evolves along the second hump, and the slow second rise and subsequent rapid LC drop. (Abstract abridged), Comment: Accepted version

We present and analyze observations of the Type Ib supernova (SN) 2019odp (a.k.a ZTF19abqwtfu) covering epochs within days of the explosion to the nebular phase at 360 d post-explosion. We discuss them in the context of recombination cooling emission for the early excess emission and consider progenitor models based on the nebular phase spectra. Our observations include photometric observations mainly in the optical and low to medium-resolution spectroscopic observations covering the complete observable time-range. We expand on existing methods to derive oxygen mass estimates from nebular phase spectroscopy. Our spectroscopic observations confirm the presence of He in the SN ejecta and we thus (re)classify it as a Type Ib supernova. From the pseudo-bolometric lightcurve we estimate a high ejecta mass $M_\text{ej} \sim 4 - 7~M_\odot$. The high ejecta mass, large nebular [O I]/[Ca II] line flux ratio ($1.2 - 1.9$) and an oxygen mass above $\gtrapprox 0.5\, M_\odot$ point towards a progenitor with pre-explosion mass higher than $18\,M_\odot$. The compact nature of the progenitor ($\lesssim 10\,R_\odot$) suggests a Wolf-Rayet (WR) star as progenitor., Comment: Submitted to A&A

Context: The extreme luminosity of gamma-ray bursts (GRBs) makes them powerful beacons for studies of the distant Universe. The most luminous bursts are typically detected at moderate/high redshift, where the volume for seeing such rare events is maximized and the star-formation activity is greater than at z = 0. For distant events, not all observations are feasible, such as at TeV energies. Aims: Here we present a spectroscopic redshift measurement for the exceptional GRB 221009A, the brightest GRB observed to date with emission extending well into the TeV regime. Methods: We used the X-shooter spectrograph at the ESO Very Large Telescope (VLT) to obtain simultaneous optical to near-IR spectroscopy of the burst afterglow 0.5 days after the explosion. Results: The spectra exhibit both absorption and emission lines from material in a host galaxy at z = 0.151. Thus GRB 221009A was a relatively nearby burst with a luminosity distance of 745 Mpc. Its host galaxy properties (star-formation rate and metallicity) are consistent with those of LGRB hosts at low redshift. This redshift measurement yields information on the energy of the burst. The inferred isotropic energy release, $E_{\rm iso} > 5 \times 10^{54}$ erg, lies at the high end of the distribution, making GRB 221009A one of the nearest and also most energetic GRBs observed to date. We estimate that such a combination (nearby as well as intrinsically bright) occurs between once every few decades to once per millennium., Comment: 9 pages, 4 figures, submitted to Astronomy & Astrophysics

We present linear polarimetry for seven hydrogen-poor superluminous supernovae (SLSNe-I). For SN 2017gci, for which we present two epochs of spectropolarimetry at +3 d and +29 d post-peak in rest frame, accompanied by four epochs of imaging polarimetry up to +108 d. The spectropolarimetry at +3 d shows increasing polarisation degree P towards the redder wavelengths and exhibits signs of axial symmetry, but at +29 d P=0 throughout the spectrum implying that the photosphere of SN 2017gci evolved from a slightly aspherical configuration to a more spherical one in the first month post-peak. However, an increase of P to 0.5% at +55 d accompanied by a different orientation of the axial symmetry compared to +3 d implies the presence of additional sources of polarisation at this phase. The increase in polarisation is possibly caused by interaction with circumstellar matter as already suggested by a knee in the light curve and a possible detection of broad Ha emission. We also analysed the sample of all 16 SLSNe-I with polarimetry to date. The data taken during the early spectroscopic phase show consistently low P indicating spherical photospheres. No clear relation between the polarimetry and spectral phase was seen when the spectra resemble Type Ic SNe during the photospheric and nebular phases. The light curve decline rate also shows no clear relation with the polarisation properties. While only slow-evolving SLSNe-I have shown non-zero P, the fast-evolving ones have not been observed at sufficiently late times to conclude that none of them exhibit changing P. However, the four SLSNe-I with increasing polarisation degree also have irregular light curve declines. For up to half of them, the photometric, spectroscopic and polarimetric properties are affected by CSM interaction. As such CSM interaction clearly plays an important role in understanding the polarimetric evolution of SLSNe-I., Comment: Accepted to A&A on 18/03/2023

Thanks to the advent of large-scale optical surveys, a diverse set of flares from the nuclear regions of galaxies has recently been discovered. These include the disruption of stars by supermassive black holes at the centers of galaxies - nuclear transients known as tidal disruption events (TDEs). Active galactic nuclei (AGN) can show extreme changes in the brightness and emission line intensities, often referred to as changing-look AGN (CLAGN). Given the physical and observational similarities, the interpretation and distinction of nuclear transients as CLAGN or TDEs remains difficult. One of the obstacles of making progress in the field is the lack of well-sampled data of long-lived nuclear outbursts in AGN. Here, we study PS16dtm, a nuclear transient in a Narrow Line Seyfert 1 (NLSy1) galaxy, which has been proposed to be a TDE candidate. Our aim is to study the spectroscopic and photometric properties of PS16dtm, in order to better understand the outbursts originating in NLSy1 galaxies. Our extensive multiwavelength follow-up that spans around 2000 days includes photometry and spectroscopy in the UV/optical, as well as mid-infrared (MIR) and X-ray observations. Furthermore, we improved an existing semiempirical model in order to reproduce the spectra and study the evolution of the spectral lines. The UV/optical light curve shows a double peak at $\sim50$ and $\sim100$ days after the first detection, and it declines and flattens afterward, reaching preoutburst levels after 2000 days of monitoring. The MIR light curve rises almost simultaneously with the optical, but unlike the UV/optical which is approaching the preoutburst levels in the last epochs of our observations, the MIR emission is still rising at the time of writing. The optical spectra show broad Balmer features and the strongest broad Fe II emission ever detected in a nuclear transient. [abridged], Comment: Accepted for publication in A&A, 19 pages and 18 figures

If high-mass stars (>20-25 Msun) are the progenitors of stripped-envelope (SE) supernovae (SNe), their massive ejecta should lead to broad, long-duration lightcurves (LCs). Instead, literature samples of SE~SNe have reported relatively narrow LCs with ejecta masses between 1-4 Msun that favor progenitors <20-25 Msun. Working with the untargeted sample of (i)PTF SNe to better constrain their rates, we search for SE~SNe with broad LCs. Using a simple LC stretch compared to a template to measure broadness, we identified eight significantly broader Type~Ibc SNe after applying quantitative sample selection criteria. The LCs, colors, and spectra of these SNe are found to evolve more slowly relative to typical Type~Ibc SNe, proportional with the stretch. Bolometric LC modeling and their nebular spectra indicate high ejecta and nickel masses, assuming radioactive decay powering. Additionally, these objects are preferentially located in low-metallicity host galaxies with high star-formation rates, which may account for their massive progenitors, as well as their relative absence from the literature. Our study thus supports the link between broad LCs (as measured by stretch) and high-mass progenitor stars in SE~SNe with independent evidence from bolometric LC modeling, nebular spectra, host environment properties, and photometric evolution. In the first systematic search of its kind using an untargeted sample, we use the stretch distribution to identify a higher than previously appreciated fraction of SE~SNe with broad LCs (~13%). Correcting for Malmquist and LC duration observational biases, we conservatively estimate that a minimum of ~6% of SE~SNe are consistent with high-mass progenitors. This result has implications for the progenitor channels of SE~SNe, late stages of massive stellar evolution, oxygen fraction in the universe, and formation channels for stellar-mass black holes.(Abridged), Comment: 38 pages, 24 figures, 12 tables. Submitted to A&A

We present the discovery and analysis of SN\,2022oqm, a Type Ic supernova (SN) detected $<1$\,day after explosion. The SN rises to a blue and short-lived (2\,days) initial peak. Early-time spectral observations of SN\,2022oqm show a hot (40,000\,K) continuum with high-ionization C and O absorption features at velocities of 4000\,km\,s$^{-1}$, while its photospheric radius expands at 20,000\,\kms, indicating a pre-existing distribution of expanding C/O material. After $\sim2.5$\,days, both the spectrum and light curves evolve into those of a typical SN Ic, with line velocities of $\sim10,000$\,km\,s$^{-1}$, in agreement with the photospheric radius evolution. The optical light curves reach a second peak at $t\approx15$\,days. By $t=60$\,days, the spectrum of \oqm\ becomes nearly nebular, displaying strong \ion{Ca}{2} and [\ion{Ca}{2}] emission with no detectable [\ion{O}{1}], marking this event as Ca-rich. The early behavior can be explained by $10^{-3}$\,\msun\ of optically thin circumstellar material (CSM) surrounding either (1) a massive compact progenitor such as a Wolf-Rayet star, (2) a massive stripped progenitor with an extended envelope, or (3) a binary system with a white dwarf. We propose that the early-time light curve is powered by both interaction of the ejecta with the optically thin CSM and shock cooling (in the massive-star scenario). The observations can be explained by CSM that is optically thick to X-ray photons, is optically thick in the lines as seen in the spectra, and is optically thin to visible-light continuum photons that come either from downscattered X-rays or from the shock-heated ejecta. Calculations show that this scenario is self-consistent., Comment: 36 pages, 23 figures. Comments are welcome ido.irani@weizmann.ac.il or idoirani@gmail.com. Accepted to ApJ

We present the results from multi-wavelength observations of a transient discovered during the follow-up of S191213g, a gravitational wave (GW) event reported by the LIGO-Virgo Collaboration as a possible binary neutron star merger in a low latency search. This search yielded SN2019wxt, a young transient in a galaxy whose sky position (in the 80\% GW contour) and distance ($\sim$150\,Mpc) were plausibly compatible with the localisation uncertainty of the GW event. Initially, the transient's tightly constrained age, its relatively faint peak magnitude ($M_i \sim -16.7$\,mag) and the $r-$band decline rate of $\sim 1$\,mag per 5\,days appeared suggestive of a compact binary merger. However, SN2019wxt spectroscopically resembled a type Ib supernova, and analysis of the optical-near-infrared evolution rapidly led to the conclusion that while it could not be associated with S191213g, it nevertheless represented an extreme outcome of stellar evolution. By modelling the light curve, we estimated an ejecta mass of $\sim 0.1\,M_\odot$, with $^{56}$Ni comprising $\sim 20\%$ of this. We were broadly able to reproduce its spectral evolution with a composition dominated by helium and oxygen, with trace amounts of calcium. We considered various progenitors that could give rise to the observed properties of SN2019wxt, and concluded that an ultra-stripped origin in a binary system is the most likely explanation. Disentangling electromagnetic counterparts to GW events from transients such as SN2019wxt is challenging: in a bid to characterise the level of contamination, we estimated the rate of events with properties comparable to those of SN2019wxt and found that $\sim 1$ such event per week can occur within the typical GW localisation area of O4 alerts out to a luminosity distance of 500\,Mpc, beyond which it would become fainter than the typical depth of current electromagnetic follow-up campaigns., Comment: By the ENGRAVE collaboration (engrave-eso.org). 35 pages, 20 figures, final version accepted by A&A

SN\,2020qlb (ZTF20abobpcb) is a hydrogen-poor superluminous supernova (SLSN-I) that is among the most luminous (maximum M$_{g} = -22.25$ mag) and that has one of the longest rise times (77 days from explosion to maximum). We estimate the total radiated energy to be $>2.1\times10^{51}$ erg. SN\,2020qlb has a well-sampled light curve that exhibits clear near and post peak undulations, a phenomenon seen in other SLSNe, whose physical origin is still unknown. We discuss the potential power source of this immense explosion as well as the mechanisms behind its observed light curve undulations. We analyze photospheric spectra and compare them to other SLSNe-I. We constructed the bolometric light curve using photometry from a large data set of observations from the Zwicky Transient Facility (ZTF), Liverpool Telescope (LT), and Neil Gehrels Swift Observatory and compare it with radioactive, circumstellar interaction and magnetar models. Model residuals and light curve polynomial fit residuals are analyzed to estimate the undulation timescale and amplitude. We also determine host galaxy properties based on imaging and spectroscopy data, including a detection of the [O III]$\lambda$4363, auroral line, allowing for a direct metallicity measurement. We rule out the Arnett $^{56}$Ni decay model for SN\,2020qlb's light curve due to unphysical parameter results. Our most favored power source is the magnetic dipole spin-down energy deposition of a magnetar. Two to three near peak oscillations, intriguingly similar to those of SN\,2015bn, were found in the magnetar model residuals with a timescale of $32\pm6$ days and an amplitude of 6$\%$ of peak luminosity. We rule out centrally located undulation sources due to timescale considerations; and we favor the result of ejecta interactions with circumstellar material (CSM) density fluctuations as the source of the undulations., Comment: 22 pages, 25 figures, submitted to A&A

The merger of two compact objects of which at least one is a neutron star is signalled by transient electromagnetic emission in a kilonova (KN). This event is accompanied by gravitational waves and possibly other radiation messengers such as neutrinos or cosmic rays. The electromagnetic emission arises from the radioactive decay of heavy $r-$process elements synthesized in the material ejected during and after the merger. In this paper we show that the analysis of KNe light curves can provide cosmological distance measurements and constrain the properties of the ejecta. In this respect, MAAT, the new Integral Field Unit in the OSIRIS spectrograph on the $10.4$ m Gran Telescopio CANARIAS (GTC), is well suited for the study of KNe by performing absolute spectro-photometry over the entire 3600-10000 Angstron spectral range. Here, we study the most representative cases regarding the scientific interest of KNe from binary neutron stars, and we evaluate the observational prospects and performance of MAAT on the GTC to do the following: a) study the impact of the equation of state on the KN light curve, and determine to what extent bounds on neutron star (NS) radii or compactness deriving from KN peak magnitudes can be identified and b) measure the Hubble constant, $H_0$, with precision improved by up to 40$\%$, when both gravitational wave data and photometric-light curves are used. In this context we discuss how the equation of state, the viewing angle, and the distance affect the precision and estimated value of $H_0$., Comment: 21 pages, 23 figures, matches published version in A&A 666, A67 (2022)

Currently time-domain astronomy can scan the entire sky on a daily basis, discovering thousands of interesting transients every night. Classifying the ever-increasing number of new transients is one of the main challenges for the astronomical community. One solution that addresses this issue is the robotically controlled Spectral Energy Distribution Machine (SEDM) which supports the Zwicky Transient Facility (ZTF). SEDM with its pipeline PYSEDM demonstrates that real-time robotic spectroscopic classification is feasible. In an effort to improve the quality of the current SEDM data, we present here two new modules, BYECR and CONTSEP. The first removes contamination from cosmic rays, and the second removes contamination from non-target light. These new modules are part of the automated PYSEDM pipeline and fully integrated with the whole process. Employing BYECR and CONTSEP modules together automatically extracts more spectra than the current PYSEDM pipeline. Using SNID classification results, the new modules show an improvement in the classification rate and accuracy of 2.8% and 1.7%, respectively, while the strength of the cross-correlation remains the same. Improvements to the SEDM astrometry would further boost the improvement of the CONTSEP module. This kind of robotic follow-up with a fully automated pipeline has the potential to provide the spectroscopic classifications for the transients discovered by ZTF and also by the Rubin Observatory's Legacy Survey of Space and Time., Comment: 9 pages, 6 figures, 2 tables, accepted for publication in PASP

We present analysis of the light curves (LCs) of 77 hydrogen-poor superluminous supernovae (SLSNe-I) discovered during the Zwicky Transient Facility Phase-I operation. We find that the majority (67\%) of the sample can be fit equally well by both magnetar and ejecta-circumstellar medium (CSM) interaction plus $^{56}$Ni decay models. This implies that LCs alone can not unambiguously constrain the physical power sources for a SLSN-I. However, 23\% of the sample show inverted V-shape, steep declining LCs or features of long rise and fast post-peak decay, which are better described by the CSM+Ni model. The remaining 10\% of the sample favor the magnetar model. Moreover, our analysis shows that the LC undulations are quite common, with a fraction of $18-44\% $ in our gold sample. Among those strongly undulating events, about 62\% of them are found to be CSM-favored, implying that the undulations tend to occur in the CSM-favored events. Undulations show a wide range in energy and duration, with median values (and 1$\sigma$ errors) being as $1.7\%^{+1.5\%}_{-0.7\%}\,\rm E_{\rm rad,total}$ and $28.8^{+14.4}_{-9.1}$\,days, respectively. Our analysis of the undulation time scales suggests that intrinsic temporal variations of the central engine can explain half of the undulating events, while CSM interaction can account for the majority of the sample. Finally, all of the well-observed He-rich SLSNe-Ib have either strongly undulating LCs or the LCs are much better fit by the CSM+Ni model. These observations imply that their progenitor stars have not had enough time to lose all of the He-envelopes before supernova explosions, and H-poor CSM are likely to present in these events., Comment: 30 pages, 20 figures, Accepted by APJ

During the Zwicky Transient Facility (ZTF) Phase-I operation, 78 hydrogen-poor superluminous supernovae (SLSNe-I) were discovered in less than three years, making up the largest sample from a single survey. This paper (Paper I) presents the data, including the optical/ultraviolet light curves and classification spectra, while Paper II in this series will focus on the detailed analysis of the light curves and modeling. Our photometry is primarily taken by the ZTF in the $g,r,i$ bands, and with additional data from other ground-based facilities and Swift. The events of our sample cover a redshift range of $z = 0.06 - 0.67$, with a median and $1\sigma$ error (16\% and 84\% percentiles) $z_{\rm med} = 0.265^{+0.143}_{-0.135}$. The peak luminosity covers $-22.8\,{\rm mag} \leq M_{g,\rm peak} \leq -19.8$\,mag, with a median value of $-21.48^{+1.13}_{-0.61}$\,mag. Their light curves evolve slowly with the mean rest-frame rise time of $t_{\rm rise} = 41.9\pm17.8$\,days. The luminosity and time scale distributions suggest that low luminosity SLSNe-I with peak luminosity $\sim -20$\,mag or extremely fast rising events ($<10$\,days) exist but are rare. We confirm previous findings that slowly rising SLSNe-I also tend to fade slowly. The rest-frame color and temperature evolution show large scatters, suggesting that the SLSN-I population may have diverse spectral energy distributions. The peak rest-frame color shows a moderate correlation with the peak absolute magnitude, i.e. brighter SLSNe-I tend to have bluer colors. With optical and ultraviolet photometry, we construct bolometric luminosity and derive a bolometric correction relation generally applicable for converting $g,r$-band photometry to bolometric luminosity for SLSNe-I., Comment: 38 pages, 25 figures, Accepted by APJ

Supernova (SN) siblings -- two or more SNe in the same parent galaxy -- are useful tools for exploring progenitor stellar populations as well as properties of the host galaxies such as distance, star formation rate, dust extinction, and metallicity. Since the average SN rate for a Milky Way-type galaxy is just one per century, a large imaging survey is required to discover an appreciable sample of SN siblings. From the wide-field Zwicky Transient Facility (ZTF) Bright Transient Survey (BTS; which aims for spectroscopic completeness for all transients which peak brighter than $r{<}$18.5 mag) we present 10 SN siblings in 5 parent galaxies. For each of these families we analyze the SN's location within the host and its underlying stellar population, finding agreement with expectations that SNe from more massive progenitors are found nearer to their host core and in regions of more active star formation. We also present an analysis of the relative rates of core collapse and thermonuclear SN siblings, finding a significantly lower ratio than past SN sibling samples due to the unbiased nature of the ZTF., Comment: 14 pages, 4 tables, 7 figures, accepted to MNRAS

The final explosive fate of massive stars, and the nature of the compact remnants they leave behind (black holes and neutron stars), are major open questions in astrophysics. Many massive stars are stripped of their outer hydrogen envelopes as they evolve. Such Wolf-Rayet (W-R) stars emit strong and rapidly expanding (v_wind>1000 km/s) winds indicating a high escape velocity from the stellar surface. A fraction of this population is also helium depleted, with spectra dominated by highly-ionized emission lines of carbon and oxygen (Types WC/WO). Evidence indicates that the most commonly-observed supernova (SN) explosions that lack hydrogen and helium (Types Ib/Ic) cannot result from massive WC/WO stars, leading some to suggest that most such stars collapse directly into black holes without a visible supernova explosions. Here, we present observations of supernova SN 2019hgp, discovered about a day after explosion. The short rise time and rapid decline place it among an emerging population of rapidly-evolving transients (RETs). Spectroscopy reveals a rich set of emission lines indicating that the explosion occurred within a nebula composed of carbon, oxygen, and neon. Narrow absorption features show that this material is expanding at relatively high velocities (>1500 km/s) requiring a compact progenitor. Our observations are consistent with an explosion of a massive WC/WO star, and suggest that massive W-R stars may be the progenitors of some rapidly evolving transients., Comment: Unedited author version, Nature in press

Type Ia supernovae (SNe~Ia) in the nearby Hubble flow are excellent distance indicators in cosmology. The Zwicky Transient Facility (ZTF) has observed a large sample of supernovae from an untargeted, rolling survey, reaching $20.8, 20.6, 20.3$ mag in $g$ $r$, and $i$-band, respectively. With a FoV of 47 sq.deg, ZTF discovered $>$ 3000 SNe~Ia in a little over 2.5 years. Here, we report on the sample of 761 spectroscopically classified SNe~Ia from the first year of operations (DR1). The sample has a median redshift $\bar z =$ 0.057, nearly a factor of two higher than the current low-$z$ sample. Our sample has a total of 934 spectra, of which 632 were obtained with the robotic SEDm on Palomar P60. We assess the potential for precision cosmology for a total of 305 SNe with redshifts from host galaxy spectra. The sample is already comparable in size to the entire combined literature low-$z$ anchor sample. The median first detection is 13.5 days before maximum light, about 10 days earlier than the median in the literature. Furthermore, six SNe from our sample are at $D_L < 80$ Mpc, for which host galaxy distances can be obtained in the JWST era, such that we have calibrator and Hubble flow SNe observed with the same instrument. In the entire duration of ZTF-I, we have observed nearly fifty SNe for which we can obtain calibrator distances, key for percent level distance scale measurements., Comment: 16 pages, 9 figures, accepted for publication in MNRAS. Data and associated information for the sample can be found at https://github.com/ZwickyTransientFacility/ztfcosmodr

We present observations of three Core-collapse supernovae (CCSNe) in elliptical hosts, detected by the Zwicky Transient Facility Bright Transient Survey (BTS). SN 2019ape is a SN Ic that exploded in the main body of a typical elliptical galaxy. Its properties are consistent with an explosion of a regular SN Ic progenitor. A secondary g-band light curve peak could indicate interaction of the ejecta with circumstellar material (CSM). An H$\alpha$-emitting source at the explosion site suggests a residual local star formation origin. SN 2018fsh and SN 2020uik are SNe II which exploded in the outskirts of elliptical galaxies. SN 2020uik shows typical spectra for SNe II, while SN 2018fsh shows a boxy nebular H$\alpha$ profile, a signature of CSM interaction. We combine these 3 SNe with 7 events from the literature and analyze their hosts as a sample. We present multi-wavelength photometry of the hosts, and compare this to archival photometry of all BTS hosts. Using the spectroscopically complete BTS we conclude that $0.3\%^{+0.3}_{-0.1}$ of all CCSNe occur in elliptical galaxies. We derive star-formation rates and stellar masses for the host-galaxies and compare them to the properties of other SN hosts. We show that CCSNe in ellipticals have larger physical separations from their hosts compared to SNe Ia in elliptical galaxies, and discuss implications for star-forming activity in elliptical galaxies., Comment: Comments are welcome. Submitted to ApJ

The division between stripped-envelope supernovae (SE-SNe) and superluminous supernovae (SLSNe) is not well defined in either photometric or spectroscopic space. While a sharp luminosity threshold has been suggested, there remains an increasing number of transitional objects that reach this threshold without the spectroscopic signatures common to SLSNe. In this work we present data and analysis on four SNe transitional between SE-SNe and SLSNe; the He-poor SNe 2019dwa and 2019cri, and the He-rich SNe 2019hge and 2019unb. Each object displays long-lived and variable photometric evolution with luminosities around the SLSN threshold of $M_r < -19.8$ mag. Spectroscopically however, these objects are similar to SE-SNe, with line velocities lower than either SE-SNe and SLSNe, and thus represent an interesting case of rare transitional events., Comment: Accepted for publication in MNRAS

Stripped-envelope supernovae (SE SNe) of Type Ib and Type Ic are thought to result from explosions of massive stars having lost their outer envelopes. The favoured explosion mechanism is by core-collapse, with the shock later revived by neutrino heating. However, there is an upper limit to the amount of radioactive Nickel-56 that such models can accomplish. Recent literature point to a tension between the maximum luminosity from such simulations and observations. We use a well characterized sample of SE SNe from the Zwicky Transient Facility (ZTF) Bright Transient Survey (BTS) to collect a sample of spectroscopically classified normal Type Ibc SNe for which we use the ZTF light curves to determine the maximum luminosity. We cull the sample further based on data quality, light-curve shape, distance and colors. The methodology of the sample construction from this BTS sample can be used for many other future investigations. In total we use 129 Type Ib or Type Ic BTS SNe with an initial rough luminosity distribution peaked at Mr = -17.61 +- 0.72, and where 36% are apparently brighter than the theoretically predicted maximum brightness of Mr = -17.8. When we further cull this sample to ensure that the SNe are normal Type Ibc with good LC data within the Hubble flow, the sample of 94 objects has Mr = -17.64 +- 0.54. A main uncertainty in absolute magnitude determinations for SNe is the host galaxy extinction correction, but the reddened objects only get more luminous after corrections. If we simply exclude objects with red, unusual or uncertain colors, we are left with 14 objects at M = -17.90 +- 0.73, whereof a handful are most certainly brighter than the suggested theoretical limit. The main result of this study is thus that normal SNe Ibc do indeed reach luminosities above $10^{42.6}$ erg/s, apparently in conflict with existing explosion models., Comment: Submitted to A&A, this version includes revisions based on referee comments

We investigate the thermal emission and extinction from dust associated with the nearby superluminous supernova (SLSN) 2018bsz. Our dataset has daily cadence and simultaneous optical and near-infrared coverage up to ~ 100 days, together with late time (+ 1.7 yr) MIR observations. At 230 days after light curve peak the SN is not detected in the optical, but shows a surprisingly strong near-infrared excess, with r - J > 3 mag and r - Ks > 5 mag. The time evolution of the infrared light curve enables us to investigate if the mid-infrared emission is from newly formed dust inside the SN ejecta, from a pre-existing circumstellar envelope, or interstellar material heated by the radiation from the SN. We find the latter two scenarios can be ruled out, and a scenario where new dust is forming in the SN ejecta at epochs > 200 days can self-consistently reproduce the evolution of the SN flux. We can fit the spectral energy distribution well at +230 d with 5 x 10^-4 solar mass of carbon dust, increasing over the following several hundred days to 10^-2 solar mass by +535 d. SN 2018bsz is the first SLSN showing evidence for dust formation within the SN ejecta, and appears to form ten times more dust than normal core-collapse SNe at similar epochs. Together with their preference for low mass, low metallicity host galaxies, we suggest that SLSNe may be a significant contributor to dust formation in the early Universe., Comment: 14 pages and 7 figures in main text, 12 pages and 6 figures in appendix. The observational data will be updated once the paper is accepted

The distribution of ejecta in young supernova remnants offers a powerful observational probe of their explosions and progenitors. Here we present a 3D reconstruction of the ejecta in SNR 0540-69.3, which is an O-rich remnant with a pulsar wind nebula located in the LMC. We use observations from VLT/MUSE to study H\beta, [O III] \lambda \lambda 4959, 5007, H\alpha, [S II] \lambda \lambda 6717, 6731, [Ar III] \lambda 7136 and [S III] \lambda 9069. This is complemented by 2D spectra from VLT/X-shooter, which also cover [O II] \lambda \lambda 3726, 3729 and [Fe II] \lambda 12567. We identify three main emission components: (i) Clumpy rings in the inner nebula (<1000 km/s) with similar morphologies in all lines; (ii) Faint extended [O III] emission dominated by an irregular ring-like structure with radius ~1600 km/s and inclination ~40 \dg, but with maximal velocities reaching ~3000 km/s; and (iii) A blob of H\alpha and H\beta located southeast of the pulsar at velocities ~1500-3500 km/s. We analyze the geometry using a clump-finding algorithm and use the clumps in the [O III] ring to estimate an age of 1146 \pm 116 years. The observations favor an interpretation of the [O III] ring as ejecta, while the origin of the H-blob is more uncertain. An alternative explanation is that it is the blown-off envelope of a binary companion. From the detection of Balmer lines in the innermost ejecta we confirm that SNR 0540 was a Type II supernova and that hydrogen was mixed down to low velocities in the explosion., Comment: 19 pages, 14 figures, accepted for publication ApJ

We present a detailed analysis of SN 2020qmp, a nearby type IIP core-collapse supernova (CCSN), discovered by the Palomar Gattini-IR (PGIR) survey in the galaxy UGC07125. We illustrate how the multiwavelength study of this event helps our general understanding of stellar progenitors and circumstellar medium (CSM) interactions in CCSNe. We also highlight the importance of near-infrared (NIR) surveys for early detections of SNe in dusty environments. SN 2020qmp displays characteristic hydrogen lines in its optical spectra, as well as a plateau in its optical LC, hallmarks of a type IIP SN. We do not detect linear polarization during the plateau phase, with a 3 sigma upper limit of 0.78%. Through hydrodynamical LC modeling and an analysis of its nebular spectra, we estimate a progenitor mass of around 11 solar masses, and an explosion energy of around 0.8e51 erg. We find that the spectral energy distribution cannot be explained by a simple CSM interaction model, assuming a constant shock velocity and steady mass-loss rate, and the excess X-ray luminosity compared with the synchrotron radio luminosity suggests deviations from equipartition. Finally, we demonstrate the advantages of NIR surveys over optical surveys for the detection of dust-obscured CCSNe in the local universe. Specifically, our simulations show that the Wide-Field Infrared Transient Explorer will detect about 14 more CCSNe out of 75 expected in its footprint within 40 Mpc, over five years than an optical survey equivalent to the Zwicky Transient Facility would detect. We have determined or constrained the main properties of SN 2020qmp and of its progenitor, highlighting the value of multiwavelength follow-up observations of nearby CCSNe. We have also shown that forthcoming NIR surveys will finally enable us to do a nearly complete census of CCSNe in the local universe., Comment: 15 pages, 14 figures, Accepted in Astronomy & Astrophysics; Revised version after referee comments and language edits from journal