Search

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

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

Gravitationally lensed Type Ia supernovae (glSNe Ia) are unique astronomical tools for studying cosmological parameters, distributions of dark matter, the astrophysics of the supernovae and the intervening lensing galaxies themselves. Only a few highly magnified glSNe Ia have been discovered by ground-based telescopes, such as the Zwicky Transient Facility (ZTF), but simulations predict the existence of a fainter, undetected population. We present a systematic search in the ZTF archive of alerts from 1 June 2019 to 1 September 2022. Using the AMPEL platform, we developed a pipeline that distinguishes candidate glSNe Ia from other variable sources. Initial cuts were applied to the ZTF alert photometry before forced photometry was obtained for the remaining candidates. Additional cuts were applied to refine the candidates based on their light curve colours, lens galaxy colours, and the resulting parameters from fits to the SALT2 SN Ia template. Candidates were also cross-matched with the DESI spectroscopic catalogue. Seven transients passed all the cuts and had an associated galaxy DESI redshift, which we present as glSN Ia candidates. While superluminous supernovae (SLSNe) cannot be fully rejected, two events, ZTF19abpjicm and ZTF22aahmovu, are significantly different from typical SLSNe and their light curves can be modelled as two-image glSN Ia systems. From this two-image modelling, we estimate time delays of 22 $\pm$ 3 and 34 $\pm$ 1 days for the two events, respectively, which suggests that we have uncovered a population with longer time delays. The pipeline is efficient and sensitive enough to parse full alert streams. It is currently being applied to the live ZTF alert stream to identify and follow-up future candidates while active. This pipeline could be the foundation for glSNe Ia searches in future surveys, like the Vera C. Rubin Observatory's Legacy Survey of Space and Time., Comment: 21 pages, 15 figures

We apply the color-magnitude intercept calibration method (CMAGIC) to the Nearby Supernova Factory SNe Ia spectrophotometric dataset. The currently existing CMAGIC parameters are the slope and intercept of a straight line fit to the first linear region in the color-magnitude diagram, which occurs over a span of approximately 30 days after maximum brightness. We define a new parameter, $\omega_{XY}$, the size of the ``bump'' feature near maximum brightness for arbitrary filters $X$ and $Y$. We find a significant correlation between the slope of the first linear region, $\beta_{XY, 1}$, in the CMAGIC diagram and $\omega_{XY}$. These results may be used to our advantage, as they are less affected by extinction than parameters defined as a function of time. Additionally, $\omega_{XY}$ is computed independently of templates. We find that current empirical templates are successful at reproducing the features described in this work, particularly SALT3, which correctly exhibits the negative correlation between slope and bump size seen in our data. In 1-D simulations, we show that the correlation between the size of the bump feature and $\beta_{XY, 1}$ can be understood as a result of chemical mixing due to large-scale Rayleigh-Taylor instabilities., Comment: 19 pages, 9 figures

Recent developments in time domain astronomy, like the Zwicky Transient Facility, have made possible a daily scan of the entire visible sky, leading to the discovery of hundreds of new transients every night. Among them, 10 to 15 are supernovae (SNe), which have to be classified prior to cosmological use. The Spectral Energy Distribution machine (SEDm), a low resolution Integral Field Spectrograph, has been designed, built, and operated to spectroscopically classify targets detected by the ZTF main camera. The current Pysedm pipeline is limited by contamination when the transient is too close to its host galaxy core; this can lead to an incorrect typing and ultimately bias the cosmological analyses, and affect the SN sample homogeneity in terms of local environment properties. We present a new scene modeler to extract the transient spectrum from its structured background, aiming at improving the typing efficiency of the SEDm. HyperGal is a fully chromatic scene modeler, which uses pre-transient photometric images to generate a hyperspectral model of the host galaxy; it is based on the CIGALE SED fitter used as a physically-motivated spectral interpolator. The galaxy model, complemented by a point source and a diffuse background component, is projected onto the SEDm spectro-spatial observation space and adjusted to observations. The full procedure is validated on 5000 simulated cubes. We introduce the contrast as the transient-to-total flux ratio at SN location. From estimated contrast distribution of real SEDm observations, we show that HyperGal correctly classifies ~95% of SNe Ia. Compared to the standard extraction method, HyperGal correctly classifies 10% more SNe Ia. The false positive rate is less than 2%, half as much as the standard extraction method. Assuming a similar contrast distribution for core-collapse SNe, HyperGal classifies 14% (11%) additional SNe II (Ibc)., Comment: 14 pages, 17 figures, submitted to Astronomy & Astrophysics

We construct a physically-parameterized probabilistic autoencoder (PAE) to learn the intrinsic diversity of type Ia supernovae (SNe Ia) from a sparse set of spectral time series. The PAE is a two-stage generative model, composed of an Auto-Encoder (AE) which is interpreted probabilistically after training using a Normalizing Flow (NF). We demonstrate that the PAE learns a low-dimensional latent space that captures the nonlinear range of features that exists within the population, and can accurately model the spectral evolution of SNe Ia across the full range of wavelength and observation times directly from the data. By introducing a correlation penalty term and multi-stage training setup alongside our physically-parameterized network we show that intrinsic and extrinsic modes of variability can be separated during training, removing the need for the additional models to perform magnitude standardization. We then use our PAE in a number of downstream tasks on SNe Ia for increasingly precise cosmological analyses, including automatic detection of SN outliers, the generation of samples consistent with the data distribution, and solving the inverse problem in the presence of noisy and incomplete data to constrain cosmological distance measurements. We find that the optimal number of intrinsic model parameters appears to be three, in line with previous studies, and show that we can standardize our test sample of SNe Ia with an RMS of $0.091 \pm 0.010$ mag, which corresponds to $0.074 \pm 0.010$ mag if peculiar velocity contributions are removed. Trained models and codes are released at \href{https://github.com/georgestein/suPAErnova}{github.com/georgestein/suPAErnova}, Comment: 23 pages, 8 Figures, 1 Table. Accepted to ApJ

We calibrate spectrophotometric optical spectra of 32 stars commonly used as standard stars, referenced to 14 stars already on the HST-based CALSPEC flux system. Observations of CALSPEC and non-CALSPEC stars were obtained with the SuperNova Integral Field Spectrograph over the wavelength range 3300 A to 9400 A as calibration for the Nearby Supernova Factory cosmology experiment. In total, this analysis used 4289 standard-star spectra taken on photometric nights. As a modern cosmology analysis, all pre-submission methodological decisions were made with the flux scale and external comparison results blinded. The large number of spectra per star allows us to treat the wavelength-by-wavelength calibration for all nights simultaneously with a Bayesian hierarchical model, thereby enabling a consistent treatment of the Type Ia supernova cosmology analysis and the calibration on which it critically relies. We determine the typical per-observation repeatability (median 14 mmag for exposures >~ 5 s), the Maunakea atmospheric transmission distribution (median dispersion of 7 mmag with uncertainty 1 mmag), and the scatter internal to our CALSPEC reference stars (median of 8 mmag). We also check our standards against literature filter photometry, finding generally good agreement over the full 12-magnitude range. Overall, the mean of our system is calibrated to the mean of CALSPEC at the level of ~ 3 mmag. With our large number of observations, careful crosschecks, and 14 reference stars, our results are the best calibration yet achieved with an integral-field spectrograph, and among the best calibrated surveys., Comment: Accepted for publication in ApJS

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

In the new era of time-domain surveys Type Ia supernovae are being caught sooner after explosion, which has exposed significant variation in their early light curves. Two driving factors for early time evolution are the distribution of nickel in the ejecta and the presence of flux excesses of various causes. We perform an analysis of the largest young SN Ia sample to date. We compare 115 SN Ia light curves from the Zwicky Transient Facility to the turtls model grid containing light curves of Chandrasekhar-mass explosions with a range of nickel masses, nickel distributions and explosion energies. We find that the majority of our observed light curves are well reproduced by Chandrasekhar-mass explosion models with a preference for highly extended nickel distributions. We identify six SNe Ia with an early-time flux excess in our g- and r-band data (four `blue' and two `red' flux excesses). We find an intrinsic rate of 18+/-11 per cent of early flux excesses in SNe Ia at z < 0.07, based on three detected flux excesses out of 30 (10 per cent) observed SNe Ia with a simulated efficiency of 57 per cent. This is comparable to rates of flux excesses in the literature but also accounts for detection efficiencies. Two of these events are mostly consistent with CSM interaction, while the other four have longer lifetimes in agreement with companion interaction and nickel-clump models. We find a higher frequency of flux excesses in 91T/99aa-like events (44+/-13 per cent).

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

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

Type Ia Supernovae (SNe Ia) are standardizable candles that allow us to measure the recent expansion rate of the Universe. Due to uncertainties in progenitor physics, potential astrophysical dependencies may bias cosmological measurements if not properly accounted for. The dependency of the intrinsic luminosity of SNe Ia with their host-galaxy environment is often used to standardize SNe Ia luminosity and is commonly parameterized as a step function. This functional form implicitly assumes two-populations of SNe Ia. In the literature, multiple environmental indicators have been considered, finding different, sometimes incompatible, step function amplitudes. We compare these indicators in the context of a two-populations model, based on their ability to distinguish the two populations. We show that local H$\alpha$-based specific star formation rate (lsSFR) and global stellar mass are better tracers than, for instance, host galaxy morphology. We show that tracer accuracy can explain the discrepancy between the observed SNe Ia step amplitudes found in the literature. Using lsSFR or global mass to distinguish the two populations can explain all other observations, though lsSFR is favoured. As lsSFR is strongly connected to age, our results favour a prompt and delayed population model. In any case, there exists two populations that differ in standardized magnitude by at least $0.121\pm0.010\,\mathrm{mag}$., Comment: Has been accepted for publication in Astronomy and Astrophysics (A&A)

We present SN 2020jfo, a Type IIP supernova in the nearby galaxy M61. Optical light curves from the Zwicky Transient Facility, complemented with data from Swift and near-IR photometry are presented. The 350-day duration bolometric light curve exhibits a relatively short (~ 65 days) plateau. This implies a moderate ejecta mass (~ 5 Msun). A series of spectroscopy is presented, including spectropolarimetric observations. The nebular spectra are dominated by Halpha but also reveal emission lines from oxygen and calcium. Comparisons to synthetic nebular spectra indicate an initial progenitor mass of about 12 Msun. Stable nickel is present in the nebular spectrum, with a super-solar Ni/Fe ratio. Several years of pre-discovery data are examined, but no signs of pre-cursor activity is found. Pre-explosion Hubble Space Telescope imaging reveals a probable progenitor star, detected only in the reddest band and is fainter than expected for stars in the 10 - 15 Msun range, in tension with the analysis of the LC and the nebular spectral modeling. We present two additional core-collapse SNe monitored by the ZTF, which also have nebular Halpha-dominated spectra. This illustrates how the absence or presence of interaction with circumstellar material affect both the LCs and in particular the nebular spectra. Type II SN 2020amv has a LC powered by CSM interaction, in particular after about 40 days when the LC is bumpy and slowly evolving. The late-time spectra show strong Halpha emission with a structure suggesting emission from a thin, dense shell. The evolution of the complex three-horn line profile is reminiscent of that observed for SN 1998S. SN 2020jfv has a poorly constrained early-time LC, but shows a transition from a hydrogen-poor Type IIb to a Type IIn, where the nebular spectrum after the light-curve rebrightening is dominated by Halpha, although with an intermediate line width., Comment: Paper on SN 2020jfo in M61, and on SNe 2020amv and 2020jfv. This is the version resubmitted to A&A after responding to first referee comments. 27 pages, 12 figures. Somewhat shortened abstract

We present observations and analysis of SN 2020cxd, a Low luminous (LL), long-lived Type IIP SN. This object was a clear outlier in the magnitude-limited SN sample recently presented by the ZTF Bright Transient Survey. We demonstrate that SN 2020cxd is an additional member of the group of LL SNe, and discuss the rarity of LL SNe in the context of the ZTF survey, and how further studies of these faintest members of the CC SN family might help understand the underlying initial mass function for stars that explode., Comment: 22 pages, 16 figures. This is the version resubmitted to A&A with replies to referee report

We study the spectral diversity of Type Ia supernovae (SNe Ia) at maximum light using high signal-to-noise spectrophotometry of 173 SNe Ia from the Nearby Supernova Factory. We decompose the diversity of these spectra into different extrinsic and intrinsic components, and we construct a nonlinear parameterization of the intrinsic diversity of SNe Ia that preserves pairings of "twin" SNe Ia. We call this parameterization the "Twins Embedding". Our methodology naturally handles highly nonlinear variability in spectra, such as changes in the photosphere expansion velocity, and uses the full spectrum rather than being limited to specific spectral line strengths, ratios or velocities. We find that the time evolution of SNe Ia near maximum light is remarkably similar, with 84.6% of the variance in common to all SNe Ia. After correcting for brightness and color, the intrinsic variability of SNe Ia is mostly restricted to specific spectral lines, and we find intrinsic dispersions as low as ~0.02 mag between 6600 and 7200 A. With a nonlinear three-dimensional model plus one dimension for color, we can explain 89.2% of the intrinsic diversity in our sample of SNe Ia, which includes several different kinds of "peculiar" SNe Ia. A linear model requires seven dimensions to explain a comparable fraction of the intrinsic diversity. We show how a wide range of previously-established indicators of diversity in SNe Ia can be recovered from the Twins Embedding. In a companion article, we discuss how these results an be applied to standardization of SNe Ia for cosmology., Comment: Accepted to ApJ

We show how spectra of Type Ia supernovae (SNe Ia) at maximum light can be used to improve cosmological distance estimates. In a companion article, we used manifold learning to build a three-dimensional parameterization of the intrinsic diversity of SNe Ia at maximum light that we call the "Twins Embedding". In this article, we discuss how the Twins Embedding can be used to improve the standardization of SNe Ia. With a single spectrophotometrically-calibrated spectrum near maximum light, we can standardize our sample of SNe Ia with an RMS of $0.101 \pm 0.007$ mag, which corresponds to $0.084 \pm 0.009$ mag if peculiar velocity contributions are removed and $0.073 \pm 0.008$ mag if a larger reference sample were obtained. Our techniques can standardize the full range of SNe Ia, including those typically labeled as peculiar and often rejected from other analyses. We find that traditional light curve width + color standardization such as SALT2 is not sufficient. The Twins Embedding identifies a subset of SNe Ia including but not limited to 91T-like SNe Ia whose SALT2 distance estimates are biased by $0.229 \pm 0.045$ mag. Standardization using the Twins Embedding also significantly decreases host-galaxy correlations. We recover a host mass step of $0.040 \pm 0.020$ mag compared to $0.092 \pm 0.024$ mag for SALT2 standardization on the same sample of SNe Ia. These biases in traditional standardization methods could significantly impact future cosmology analyses if not properly taken into account., Comment: Accepted to ApJ

Helium is expected to be present in the ejecta of some hydrogen-poor superluminous supernovae (SLSN-I). However, so far only one event, PTF10hgi has been identified with He features in its photospheric spectra (Quimby et al. 2018). We present the discovery of a new He-rich SLSN-I, ZTF19aawfbtg (SN2019hge) at $z=0.0866$. This event has more than 10 optical spectra at phases from $-41$ to $+103$\,days relative to the peak, most of which match well with that of PTF10hgi. Confirmation comes from a near-IR spectrum taken at $+34$ days, revealing He I features with P-Cygni profiles at 1.083 and 2.058$\mu$m. Using the optical spectra of PTF10hgi and SN2019hge as templates, we examine 70 SLSN-I discovered by ZTF in the first two years of operation and found additional five SLSN-I with distinct He-features. The excitation of He\,I atoms in normal core collapse supernovae requires non-thermal radiation, as proposed by previous studies. These He-rich events can not be explained by traditional $^{56}$Ni mixing model because of their blue spectra, high peak luminosity and slow rising time scales. Magnetar models offer a possible solution since pulsar winds naturally generate high energy particles as sources of non-thermal excitation. An alternative model is ejecta interaction with H-poor CSM which may be supported by the observed light curve undulations. These six SLSN-Ib appear to have relatively low-peak luminosities (rest-frame $M_g = -20.06\pm0.16$ mag)., Comment: 12 pages, 6 figures, Accepted for publication in ApJL. Matched with the published version

The detailed nature of type Ia supernovae (SNe Ia) remains uncertain, and as survey statistics increase, the question of astrophysical systematic uncertainties arises, notably that of the evolution of SN Ia populations. We study the dependence on redshift of the SN Ia light-curve stretch, a purely intrinsic SN property, to probe its potential redshift drift. The SN stretch has been shown to be strongly correlated with the SN environment, notably with stellar age tracers. We modeled the underlying stretch distribution as a function of redshift, using the evolution of the fraction of young and old SNe Ia as predicted using the SNfactory dataset, and assuming a constant underlying stretch distribution for each age population consisting of Gaussian mixtures. We tested our prediction against published samples that were cut to have marginal magnitude selection effects so that any observed change is indeed astrophysical and not observational in origin. In this first study, there are indications that the underlying SN Ia stretch distribution evolves as a function of redshift, and that the age drifting model is a better description of the data than any time-constant model, including the sample-based asymmetric distributions that are often used to correct Malmquist bias at a significance higher than 5 $\sigma$. The favored underlying stretch model is a bimodal one, composed of a high-stretch mode shared by both young and old environments, and a low-stretch mode that is exclusive to old environments. The precise effect of the redshift evolution of the intrinsic properties of a SN Ia population on cosmology remains to be studied. The astrophysical drift of the SN stretch distribution does affect current Malmquist bias corrections and hence the distances that are derived using SNe that are affected by observational selection effects. This bias increases with surveys covering larger redshift ranges., Comment: 10 pages, 8 figures; submitted to A&A and accepted on 17th March 2021. Replacing the previous version, adding section 2.1 with 1 more page and 1 more figure

Early observations of Type Ia supernovae (SNe$\,$Ia) provide essential clues for understanding the progenitor system that gave rise to the terminal thermonuclear explosion. We present exquisite observations of SN$\,$2019yvq, the second observed SN$\,$Ia, after iPTF$\,$14atg, to display an early flash of emission in the ultraviolet (UV) and optical. Our analysis finds that SN$\,$2019yvq was unusual, even when ignoring the initial flash, in that it was moderately underluminous for an SN$\,$Ia ($M_g \approx -18.5\,$mag at peak) yet featured very high absorption velocities ($v \approx 15,000\,\mathrm{km\,s}^{-1}$ for Si II $\lambda$6355 at peak). We find that many of the observational features of SN$\,$2019yvq, aside from the flash, can be explained if the explosive yield of radioactive $^{56}\mathrm{Ni}$ is relatively low (we measure $M_{^{56}\mathrm{Ni}} = 0.31 \pm 0.05\,M_\odot$) and it and other iron-group elements are concentrated in the innermost layers of the ejecta. To explain both the UV/optical flash and peak properties of SN$\,$2019yvq we consider four different models: interaction between the SN ejecta and a nondegenerate companion, extended clumps of $^{56}\mathrm{Ni}$ in the outer ejecta, a double-detonation explosion, and the violent merger of two white dwarfs. Each of these models has shortcomings when compared to the observations; it is clear additional tuning is required to better match SN$\,$2019yvq. In closing, we predict that the nebular spectra of SN$\,$2019yvq will feature either H or He emission, if the ejecta collided with a companion, strong [Ca II] emission, if it was a double detonation, or narrow [O I] emission, if it was due to a violent merger., Comment: 27 pages, 14 figures, accepted in ApJ

The Nearby Supernova Factory has made spectrophotometric observations of Type Ia supernovae since $2004$. This work presents an interim version of the data produced, including $210$ supernovae observed between $2004$ and $2013$., Comment: 5 pages

Type Ia Supernovae have yet again the opportunity to revolutionize the field of cosmology as the new generation of surveys are acquiring thousands of nearby SNeIa opening a new era in cosmology: the direct measurement of the growth of structure parametrized by $fD$. This method is based on the SNeIa peculiar velocities derived from the residual to the Hubble law as direct tracers of the full gravitational potential caused by large scale structure. With this technique, we could probe not only the properties of dark energy, but also the laws of gravity. In this paper we present the analytical framework and forecasts. We show that ZTF and LSST will be able to reach 5\% precision on $fD$ by 2027. Our analysis is not significantly sensitive to photo-typing, but known selection functions and spectroscopic redshifts are mandatory. We finally introduce an idea of a dedicated spectrograph that would get all the required information in addition to boost the efficiency to each SNeIa so that we could reach the 5\% precision within the first two years of LSST operation and the few percent level by the end of the survey., Comment: White paper CNRS/IN2P3 Prospectives 2020

The Zwicky Transient Facility (ZTF) is performing a three-day cadence survey of the visible Northern sky (~3$\pi$). The transient candidates found in this survey are announced via public alerts. As a supplementary product ZTF is also conducting a large spectroscopic campaign: the ZTF Bright Transient Survey (BTS). The goal of the BTS is to spectroscopically classify all extragalactic transients brighter than 18.5 mag at peak brightness and immediately announce those classifications to the public. Extragalactic discoveries from ZTF are predominantly Supernovae (SNe). The BTS is the largest flux-limited SN survey to date. Here we present a catalog of the761 SNe that were classified during the first nine months of the survey (2018 Apr. 1 to 2018 Dec. 31). The BTS SN catalog contains redshifts based on SN template matching and spectroscopic host galaxy redshifts when available. Based on this data we perform an analysis of the redshift completeness of local galaxy catalogs, dubbed as the Redshift Completeness Fraction (RCF; the number of SN host galaxies with known spectroscopic redshift prior to SN discovery divided by the total number of SN hosts). In total, we identify the host galaxies of 512 Type Ia supernovae, 227 of which have known spectroscopic redshifts, yielding an RCF estimate of $44\% \pm1\%$. We find a steady decrease in the RCF with increasing distance in the local universe. For z<0.05, or ~200 Mpc, we find RCF=0.6, which has important ramifications when searching for multimessenger astronomical events. Prospects for dramatically increasing the RCF are limited to new multi-fiber spectroscopic instruments, or wide-field narrowband surveys. We find that existing galaxy redshift catalogs are only $50\%$ complete at $r\approx16.9$ mag. Pushing this limit several magnitudes deeper will pay huge dividends when searching for electromagnetic counterparts to gravitational wave events., Comment: 22 pages, 10 figures. Submitted to ApJ

Type Ia Supernovae (SNe Ia) are widely used to measure the expansion of the Universe. Improving distance measurements of SNe Ia is one technique to better constrain the acceleration of expansion and determine its physical nature. This document develops a new SNe Ia spectral energy distribution (SED) model, called the SUpernova Generator And Reconstructor (SUGAR), which improves the spectral description of SNe Ia, and consequently could improve the distance measurements. This model is constructed from SNe Ia spectral properties and spectrophotometric data from The Nearby Supernova Factory collaboration. In a first step, a PCA-like method is used on spectral features measured at maximum light, which allows us to extract the intrinsic properties of SNe Ia. Next, the intrinsic properties are used to extract the average extinction curve. Third, an interpolation using Gaussian Processes facilitates using data taken at different epochs during the lifetime of a SN Ia and then projecting the data on a fixed time grid. Finally, the three steps are combined to build the SED model as a function of time and wavelength. This is the SUGAR model. The main advancement in SUGAR is the addition of two additional parameters to characterize SNe Ia variability. The first is tied to the properties of SNe Ia ejecta velocity, the second is correlated with their calcium lines. The addition of these parameters, as well as the high quality the Nearby Supernova Factory data, makes SUGAR an accurate and efficient model for describing the spectra of normal SNe Ia as they brighten and fade. The performance of this model makes it an excellent SED model for experiments like ZTF, LSST or WFIRST., Comment: 25 pages, 27 figures

As a candidate 'super-Chandrasekhar' or 09dc-like Type Ia supernova (SN Ia), SN 2012dn shares many characteristics with other members of this remarkable class of objects but lacks their extraordinary luminosity. Here, we present and discuss the most comprehensive optical data set of this SN to date, comprised of a densely sampled series of early-time spectra obtained within the Nearby Supernova Factory project, plus photometry and spectroscopy obtained at the VLT about 1 yr after the explosion. The light curves, colour curves, spectral time series and ejecta velocities of SN 2012dn are compared with those of other 09dc-like and normal SNe Ia, the overall variety within the class of 09dc-like SNe Ia is discussed, and new criteria for 09dc-likeness are proposed. Particular attention is directed to additional insight that the late-phase data provide. The nebular spectra show forbidden lines of oxygen and calcium, elements that are usually not seen in late-time spectra of SNe Ia, while the ionisation state of the emitting iron plasma is low, pointing to low ejecta temperatures and high densities. The optical light curves are characterised by an enhanced fading starting ~60 d after maximum and very low luminosities in the nebular phase, which is most readily explained by unusually early formation of clumpy dust in the ejecta. Taken together, these effects suggest a strongly perturbed ejecta density profile, which might lend support to the idea that 09dc-like characteristics arise from a brief episode of interaction with a hydrogen-deficient envelope during the first hours or days after the explosion., Comment: 16 pages, 11 figures, MNRAS in press, missing line added in Table 1

Both multi-messenger astronomy and new high-throughput wide-field surveys require flexible tools for the selection and analysis of astrophysical transients. We here introduce the Alert Management, Photometry and Evaluation of Lightcurves (AMPEL) system, an analysis framework designed for high-throughput surveys and suited for streamed data. AMPEL combines the functionality of an alert broker with a generic framework capable of hosting user-contributed code, that encourages provenance and keeps track of the varying information states that a transient displays. The latter concept includes information gathered over time and data policies such as access or calibration levels. We describe a novel ongoing real-time multi-messenger analysis using AMPEL to combine IceCube neutrino data with the alert streams of the Zwicky Transient Facility (ZTF). We also reprocess the first four months of ZTF public alerts, and compare the yields of more than 200 different transient selection functions to quantify efficiencies for selecting Type Ia supernovae that were reported to the Transient Name Server (TNS). We highlight three channels suitable for (1) the collection of a complete sample of extragalactic transients, (2) immediate follow-up of nearby transients and (3) follow-up campaigns targeting young, extragalactic transients. We confirm ZTF completeness in that all TNS supernovae positioned on active CCD regions were detected. AMPEL can assist in filtering transients in real time, running alert reaction simulations, the reprocessing of full datasets as well as in the final scientific analysis of transient data. This text introduces how users can design their own channels for inclusion in the AMPEL live instance that parses the ZTF stream and the real-time submission of high quality extragalactic supernova candidates to the TNS., Comment: Updated to match version accepted for publication in A&A

In the upcoming decade cadenced wide-field imaging surveys will increase the number of identified $z<0.3$ Type~Ia supernovae (SNe~Ia) from the hundreds to the hundreds of thousands. The increase in the number density and solid-angle coverage of SNe~Ia, in parallel with improvements in the standardization of their absolute magnitudes, now make them competitive probes of the growth of structure and hence of gravity. The peculiar velocity power spectrum is sensitive to the growth index $\gamma$, which captures the effect of gravity on the linear growth of structure through the relation $f=\Omega_M^\gamma$. We present the first projections for the precision in $\gamma$ for a range of realistic SN peculiar-velocity survey scenarios. In the next decade the peculiar velocities of SNe~Ia in the local $z<0.3$ Universe will provide a measure of $\gamma$ to $\pm 0.01$ precision that can definitively distinguish between General Relativity and leading models of alternative gravity., Comment: Decadal Survey White Paper submission

Current time domain facilities are discovering hundreds of new galactic and extra-galactic transients every week. Classifying the ever-increasing number of transients is challenging, yet crucial to further our understanding of their nature, discover new classes, or ensuring sample purity, for instance, for Supernova Ia cosmology. The Zwicky Transient Facility is one example of such a survey. In addition, it has a dedicated very-low resolution spectrograph, the SEDMachine, operating on the Palomar 60-inch telescope. This spectrograph's primary aim is object classification. In practice most, if not all, transients of interest brighter than ~19 mag are typed. This corresponds to approximately 10 to 15 targets a night. In this paper, we present a fully automated pipeline for the SEDMachine. This pipeline has been designed to be fast, robust, stable and extremely flexible. pysedm enables the fully automated spectral extraction of a targeted point source object in less than 5 minutes after the end of the exposure. The spectral color calibration is accurate at the few percent level. In the 19 weeks since pysedm entered production in early August of 2018, we have classified, among other objects, about 400 Type Ia supernovae and 140 Type II supernovae. We conclude that low resolution, fully automated spectrographs such as the `SEDMachine with pysedm' installed on 2-m class telescopes within the southern hemisphere could allow us to automatically and simultaneously type and obtain a redshift for most (if not all) bright transients detected by LSST within z<0.2, notably potentially all Type Ia Supernovae. In comparison to the current SEDM design, this would require higher spectral resolution (R~1000) and slightly improved throughput. With this perspective in mind, pysedm has been designed to easily be adaptable to any IFU-like spectrograph (see https://github.com/MickaelRigault/pysedm)., Comment: Submitted to Astronomy and Astrophysics

Type Ia supernova cosmology depends on the ability to fit and standardize observations of supernova magnitudes with an empirical model. We present here a series of new models of Type Ia Supernova spectral time series that capture a greater amount of supernova diversity than possible with the models that are currently customary. These are entitled SuperNova Empirical MOdels (\textsc{SNEMO}\footnote{https://snfactory.lbl.gov/snemo}). The models are constructed using spectrophotometric time series from $172$ individual supernovae from the Nearby Supernova Factory, comprising more than $2000$ spectra. Using the available observations, Gaussian Processes are used to predict a full spectral time series for each supernova. A matrix is constructed from the spectral time series of all the supernovae, and Expectation Maximization Factor Analysis is used to calculate the principal components of the data. K-fold cross-validation then determines the selection of model parameters and accounts for color variation in the data. Based on this process, the final models are trained on supernovae that have been dereddened using the Fitzpatrick and Massa extinction relation. Three final models are presented here: \textsc{SNEMO2}, a two-component model for comparison with current Type~Ia models; \textsc{SNEMO7}, a seven component model chosen for standardizing supernova magnitudes which results in a total dispersion of $0.100$~mag for a validation set of supernovae, of which $0.087$~mag is unexplained (a total dispersion of $0.113$~mag with unexplained dispersion of $0.097$~mag is found for the total set of training and validation supernovae); and \textsc{SNEMO15}, a comprehensive $15$ component model that maximizes the amount of spectral time series behavior captured., Comment: 51 page, 19 figures, accepted in ApJ