Search

We apply the color-magnitude intercept calibration method (CMAGIC) to the Nearby Supernova Factory SNe Ia spectrophotometric data set. The currently existing CMAGIC parameters are the slope and intercept of a straight line fit to the linear region in the color-magnitude diagram, which occurs over a span of approximately 30 days after maximum brightness. We define a new parameter, ω 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 linear region, β XY, in the CMAGIC diagram and ω 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, ω 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 1D simulations, we show that the correlation between the size of the “bump” feature and β XY can be understood as a result of chemical mixing due to large-scale Rayleigh-Taylor instabilities.

We calibrate spectrophotometric optical spectra of 32 stars commonly used as standard stars, referenced to 14 stars already on the Hubble Space Telescope-based CALSPEC flux system. Observations of CALSPEC and non-CALSPEC stars were obtained with the SuperNova Integral Field Spectrograph over the wavelength range 3300-9400 Å 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 presubmission 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 mag 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 cross-checks, and 14 reference stars, our results are the best calibration yet achieved with an integral-field spectrograph, and among the best calibrated surveys.

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 autoencoder that is interpreted probabilistically after training using a normalizing flow. 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 multistage 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 the 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 ± 0.010 mag, which corresponds to 0.074 ± 0.010 mag if peculiar velocity contributions are removed. Trained models and codes are released at https://github.com/georgestein/suPAErnova.

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 SNe from an untargeted, rolling survey, reaching 20.8, 20.6, and 20.3 mag in g r, and i band, respectively. With an FoV of 47 deg2, ZTF discovered > 3000 SNe Ia in a little over 2.5 yr. 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 2 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 d before maximum light, about 10 d earlier than the median in the literature. Furthermore, six SNe from our sample are at DL < 80 Mpc, for which host galaxy distances can be obtained in the JAMES WEBB SPACE TELESCOPE 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 50 SNe for which we can obtain calibrator distances, key for per cent level distance scale measurements.

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α-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 identify 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 ± 0.010 mag.

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 SALT2.4 light-curve stretch, which is 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σ. 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, however, and thereby the distances that are derived based on SN that are affected by observational selection effects. We highlight that this bias will increase with surveys covering increasingly larger redshift ranges, which is particularly important for the Large Synoptic Survey Telescope.

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 0.007 mag, which corresponds to 0.084 0.009 mag if peculiar velocity contributions are removed and to 0.073 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 0.045 mag. Standardization using the Twins Embedding also significantly decreases host-galaxy correlations. We recover a host mass step of 0.040 0.020 mag compared to 0.092 0.026 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.

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 Å. 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 can be applied to the standardization of SNe Ia for cosmology.

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 0.007 mag, which corresponds to 0.084 0.009 mag if peculiar velocity contributions are removed and to 0.073 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 0.045 mag. Standardization using the Twins Embedding also significantly decreases host-galaxy correlations. We recover a host mass step of 0.040 0.020 mag compared to 0.092 0.026 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.

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 Å. 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 can be applied to the standardization of SNe Ia for cosmology.

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 SALT2.4 light-curve stretch, which is 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σ. 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, however, and thereby the distances that are derived based on SN that are affected by observational selection effects. We highlight that this bias will increase with surveys covering increasingly larger redshift ranges, which is particularly important for the Large Synoptic Survey Telescope.

As part of an on-going effort to identify, understand and correct for astrophysics biases in the standardization of Type Ia supernovae (SN Ia) for cosmology, we have statistically classified a large sample of nearby SNe Ia into those that are located in predominantly younger or older environments. This classification is based on the specific star formation rate measured within a projected distance of 1 kpc from each SN location (LsSFR). This is an important refinement compared to using the local star formation rate directly, as it provides a normalization for relative numbers of available SN progenitors and is more robust against extinction by dust. We find that the SNe Ia in predominantly younger environments are ΔY = 0.163 ± 0.029 mag (5.7σ) fainter than those in predominantly older environments after conventional light-curve standardization. This is the strongest standardized SN Ia brightness systematic connected to the host-galaxy environment measured to date. The well-established step in standardized brightnesses between SNe Ia in hosts with lower or higher total stellar masses is smaller, at ΔM = 0.119 ± 0.032 mag (4.5σ), for the same set of SNe Ia. When fit simultaneously, the environment-age offset remains very significant, with ΔY = 0.129 ± 0.032 mag (4.0σ), while the global stellar mass step is reduced to ΔM = 0.064 ± 0.029 mag (2.2σ). Thus, approximately 70% of the variance from the stellar mass step is due to an underlying dependence on environment-based progenitor age. Also, we verify that using the local star formation rate alone is not as powerful as LsSFR at sorting SNe Ia into brighter and fainter subsets. Standardization that only uses the SNe Ia in younger environments reduces the total dispersion from 0.142 ± 0.008 mag to 0.120 ± 0.010 mag. We show that as environment-ages evolve with redshift, a strong bias, especially on the measurement of the derivative of the dark energy equation of state, can develop. Fortunately, data that measure and correct for this effect using our local specific star formation rate indicator, are likely to be available for many next-generation SN Ia cosmology experiments.

As part of an on-going effort to identify, understand and correct for astrophysics biases in the standardization of Type Ia supernovae (SN Ia) for cosmology, we have statistically classified a large sample of nearby SNe Ia into those that are located in predominantly younger or older environments. This classification is based on the specific star formation rate measured within a projected distance of 1 kpc from each SN location (LsSFR). This is an important refinement compared to using the local star formation rate directly, as it provides a normalization for relative numbers of available SN progenitors and is more robust against extinction by dust. We find that the SNe Ia in predominantly younger environments are ΔY = 0.163 ± 0.029 mag (5.7σ) fainter than those in predominantly older environments after conventional light-curve standardization. This is the strongest standardized SN Ia brightness systematic connected to the host-galaxy environment measured to date. The well-established step in standardized brightnesses between SNe Ia in hosts with lower or higher total stellar masses is smaller, at ΔM = 0.119 ± 0.032 mag (4.5σ), for the same set of SNe Ia. When fit simultaneously, the environment-age offset remains very significant, with ΔY = 0.129 ± 0.032 mag (4.0σ), while the global stellar mass step is reduced to ΔM = 0.064 ± 0.029 mag (2.2σ). Thus, approximately 70% of the variance from the stellar mass step is due to an underlying dependence on environment-based progenitor age. Also, we verify that using the local star formation rate alone is not as powerful as LsSFR at sorting SNe Ia into brighter and fainter subsets. Standardization that only uses the SNe Ia in younger environments reduces the total dispersion from 0.142 ± 0.008 mag to 0.120 ± 0.010 mag. We show that as environment-ages evolve with redshift, a strong bias, especially on the measurement of the derivative of the dark energy equation of state, can develop. Fortunately, data that measure and correct for this effect using our local specific star formation rate indicator, are likely to be available for many next-generation SN Ia cosmology experiments.

High-cadence transient surveys are able to capture supernovae closer to their first light than ever before. Applying analytical models to such early emission, we can constrain the progenitor stars’ properties. In this paper, we present observations of SN 2018fif (ZTF 18abokyfk). The supernova was discovered close to first light and monitored by the Zwicky Transient Facility (ZTF) and the Neil Gehrels Swift Observatory. Early spectroscopic observations suggest that the progenitor of SN 2018fif was surrounded by relatively small amounts of circumstellar material compared to all previous cases. This particularity, coupled with the high-cadence multiple-band coverage, makes it a good candidate to investigate using shock-cooling models. We employ the SOPRANOS code, an implementation of the model by Sapir & Waxman and its extension to early times by Morag et al. Compared with previous implementations, SOPRANOS has the advantage of including a careful account of the limited temporal validity domain of the shock-cooling model as well as allowing usage of the entirety of the early UV data. We find that the progenitor of SN 2018fif was a large red supergiant with a radius of R = 744.0-+128.0183.0 R☉ and an ejected mass of Mej = 9.3-+5.80.4 M☉. Our model also gives information on the explosion epoch, the progenitor’s inner structure, the shock velocity, and the extinction. The distribution of radii is double-peaked, with smaller radii corresponding to lower values of the extinction, earlier recombination times, and a better match to the early UV data. If these correlations persist in future objects, denser spectroscopic monitoring constraining the time of recombination, as well as accurate UV observations (e.g., with ULTRASAT), will help break the extinction/radius degeneracy and independently determine both.

© 2020. The American Astronomical Society. All rights reserved.. We present a survey of the early evolution of 12 Type IIn supernovae (SNe IIn) at ultraviolet and visible light wavelengths. We use this survey to constrain the geometry of the circumstellar material (CSM) surrounding SN IIn explosions, which may shed light on their progenitor diversity. In order to distinguish between aspherical and spherical CSM, we estimate the blackbody radius temporal evolution of the SNe IIn of our sample, following the method introduced by Soumagnac et al. We find that higher-luminosity objects tend to show evidence for aspherical CSM. Depending on whether this correlation is due to physical reasons or to some selection bias, we derive a lower limit between 35% and 66% for the fraction of SNe IIn showing evidence for aspherical CSM. This result suggests that asphericity of the CSM surrounding SNe IIn is common - consistent with data from resolved images of stars undergoing considerable mass loss. It should be taken into account for more realistic modeling of these events.

© 2020. The American Astronomical Society. All rights reserved.. We present a survey of the early evolution of 12 Type IIn supernovae (SNe IIn) at ultraviolet and visible light wavelengths. We use this survey to constrain the geometry of the circumstellar material (CSM) surrounding SN IIn explosions, which may shed light on their progenitor diversity. In order to distinguish between aspherical and spherical CSM, we estimate the blackbody radius temporal evolution of the SNe IIn of our sample, following the method introduced by Soumagnac et al. We find that higher-luminosity objects tend to show evidence for aspherical CSM. Depending on whether this correlation is due to physical reasons or to some selection bias, we derive a lower limit between 35% and 66% for the fraction of SNe IIn showing evidence for aspherical CSM. This result suggests that asphericity of the CSM surrounding SNe IIn is common - consistent with data from resolved images of stars undergoing considerable mass loss. It should be taken into account for more realistic modeling of these events.

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 (Mg ≈ 18.5 mag at peak) yet featured very high absorption velocities (v ≈ 15,000 km s-1 for Si ii λ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 56Ni is relatively low (we measure M56Ni=0.31±0.05,M⊙) 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 56Ni 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.

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 (Mg ≈ 18.5 mag at peak) yet featured very high absorption velocities (v ≈ 15,000 km s-1 for Si ii λ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 56Ni is relatively low (we measure M56Ni=0.31±0.05,M⊙) 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 56Ni 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.

© P.-F. Léget et al. 2020. Context. 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. Aims. 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. Methods. This model was constructed from SNe Ia spectral properties and spectrophotometric data from the Nearby Supernova Factory collaboration. In a first step, a principal component analysis-like method was used on spectral features measured at maximum light, which allowed us to extract the intrinsic properties of SNe Ia. Next, the intrinsic properties were used to extract the average extinction curve. Third, an interpolation using Gaussian processes facilitated using data taken at different epochs during the lifetime of an SN Ia and then projecting the data on a fixed time grid. Finally, the three steps were combined to build the SED model as a function of time and wavelength. This is the SUGAR model. Results. 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 and the second correlates with their calcium lines. The addition of these parameters, as well as the high quality of 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. Conclusions. The performance of this model makes it an excellent SED model for experiments like the Zwicky Transient Facility, the Large Synoptic Survey Telescope, or the Wide Field Infrared Survey Telescope.

© P.-F. Léget et al. 2020. Context. 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. Aims. 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. Methods. This model was constructed from SNe Ia spectral properties and spectrophotometric data from the Nearby Supernova Factory collaboration. In a first step, a principal component analysis-like method was used on spectral features measured at maximum light, which allowed us to extract the intrinsic properties of SNe Ia. Next, the intrinsic properties were used to extract the average extinction curve. Third, an interpolation using Gaussian processes facilitated using data taken at different epochs during the lifetime of an SN Ia and then projecting the data on a fixed time grid. Finally, the three steps were combined to build the SED model as a function of time and wavelength. This is the SUGAR model. Results. 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 and the second correlates with their calcium lines. The addition of these parameters, as well as the high quality of 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. Conclusions. The performance of this model makes it an excellent SED model for experiments like the Zwicky Transient Facility, the Large Synoptic Survey Telescope, or the Wide Field Infrared Survey Telescope.

Early-time observations of Type Ia supernovae (SNe Ia) are essential to constrain the properties of their progenitors. In this paper, we present high-quality light curves of 127 SNe Ia discovered by the Zwicky Transient Facility (ZTF) in 2018. We describe our method to perform forced point-spread function photometry, which can be applied to other types of extragalactic transients. With a planned cadence of six observations per night (three g + three r), all of the 127 SNe Ia are detected in both g and r bands more than 10 days (in the rest frame) prior to the epoch of g-band maximum light. The redshifts of these objects range from z = 0.0181 to 0.165; the median redshift is 0.074. Among the 127 SNe, 50 are detected at least 14 days prior to maximum light (in the rest frame), with a subset of nine objects being detected more than 17 days before g-band peak. This is the largest sample of young SNe Ia collected to date; it can be used to study the shape and color evolution of the rising light curves in unprecedented detail. We discuss six peculiar events in this sample: one 02cx-like event ZTF18abclfee (SN 2018crl), one Ia-CSM SN ZTF18aaykjei (SN 2018cxk), and four objects with possible super-Chandrasekhar mass progenitors: ZTF18abhpgje (SN 2018eul), ZTF18abdpvnd (SN 2018dvf), ZTF18aawpcel (SN 2018cir), and ZTF18abddmrf (SN 2018dsx).

The third observing run by LVC has brought the discovery of many compact binary coalescences. Following the detection of the first binary neutron star merger in this run (LIGO/Virgo S190425z), we performed a dedicated follow-up campaign with the Zwicky Transient Facility (ZTF) and Palomar Gattini-IR telescopes. The initial skymap of this single-detector gravitational wave (GW) trigger spanned most of the sky observable from Palomar Observatory. Covering 8000 deg2 of the initial skymap over the next two nights, corresponding to 46% integrated probability, ZTF system achieved a depth of ≈21 m AB in g-and r-bands. Palomar Gattini-IR covered 2200 square degrees in J-band to a depth of 15.5 mag, including 32% integrated probability based on the initial skymap. The revised skymap issued the following day reduced these numbers to 21% for the ZTF and 19% for Palomar Gattini-IR. We narrowed 338,646 ZTF transient "alerts" over the first two nights of observations to 15 candidate counterparts. Two candidates, ZTF19aarykkb and ZTF19aarzaod, were particularly compelling given that their location, distance, and age were consistent with the GW event, and their early optical light curves were photometrically consistent with that of kilonovae. These two candidates were spectroscopically classified as young core-collapse supernovae. The remaining candidates were ruled out as supernovae. Palomar Gattini-IR did not identify any viable candidates with multiple detections only after merger time. We demonstrate that even with single-detector GW events localized to thousands of square degrees, systematic kilonova discovery is feasible.

In the upcoming decade cadenced wide-field imaging surveys will increase thenumber of identified $z

In the upcoming decade cadenced wide-field imaging surveys will increase thenumber of identified $z

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 Very Large Telescope 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 ionization state of the emitting iron plasma is low, pointing to low ejecta temperatures and high densities. The optical light curves are characterized 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.

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 Very Large Telescope 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 ionization state of the emitting iron plasma is low, pointing to low ejecta temperatures and high densities. The optical light curves are characterized 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.

The Zwicky Transient Facility (ZTF) is a new optical time-domain survey that uses the Palomar 48 inch Schmidt telescope. A custom-built wide-field camera provides a 47 deg 2 field of view and 8 s readout time, yielding more than an order of magnitude improvement in survey speed relative to its predecessor survey, the Palomar Transient Factory. We describe the design and implementation of the camera and observing system. The ZTF data system at the Infrared Processing and Analysis Center provides near-real-time reduction to identify moving and varying objects. We outline the analysis pipelines, data products, and associated archive. Finally, we present on-sky performance analysis and first scientific results from commissioning and the early survey. ZTF’s public alert stream will serve as a useful precursor for that of the Large Synoptic Survey Telescope.

SN Ia 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 SN Ia spectral time series that capture a greater amount of supernova diversity than is possible with the models that are currently customary. These are entitled SuperNova Empirical MOdels (SNEMO; 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: SNEMO2, a two-component model for comparison with current Type Ia models; 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 an unexplained dispersion of 0.097 mag is found for the total set of training and validation supernovae); and SNEMO15, a comprehensive 15-component model that maximizes the amount of spectral time-series behavior captured.

SN Ia 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 SN Ia spectral time series that capture a greater amount of supernova diversity than is possible with the models that are currently customary. These are entitled SuperNova Empirical MOdels (SNEMO; 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: SNEMO2, a two-component model for comparison with current Type Ia models; 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 an unexplained dispersion of 0.097 mag is found for the total set of training and validation supernovae); and SNEMO15, a comprehensive 15-component model that maximizes the amount of spectral time-series behavior captured.

We present the discovery and measurements of a gravitationally lensed supernova (SN) behind the galaxy cluster MOO J1014+0038. Based on multi-band Hubble Space Telescope and Very Large Telescope (VLT) photometry of the supernova, and VLT spectroscopy of the host galaxy, we find a 97.5% probability that this SN is a SN Ia, and a 2.5% chance of a CC SN. Our typing algorithm combines the shape and color of the light curve with the expected rates of each SN type in the host galaxy. With a redshift of 2.2216, this is the highest redshift SN Ia discovered with a spectroscopic host-galaxy redshift. A further distinguishing feature is that the lensing cluster, at redshift 1.23, is the most distant to date to have an amplified SN. The SN lies in the middle of the color and light-curve shape distributions found at lower redshift, disfavoring strong evolution to z = 2.22. We estimate an amplification due to gravitational lensing of (1.10 0.23 mag) - compatible with the value estimated from the weak-lensing-derived mass and the mass-concentration relation from ΛCDM simulations - making it the most amplified SN Ia discovered behind a galaxy cluster.

We present the discovery and measurements of a gravitationally lensed supernova (SN) behind the galaxy cluster MOO J1014+0038. Based on multi-band Hubble Space Telescope and Very Large Telescope (VLT) photometry of the supernova, and VLT spectroscopy of the host galaxy, we find a 97.5% probability that this SN is a SN Ia, and a 2.5% chance of a CC SN. Our typing algorithm combines the shape and color of the light curve with the expected rates of each SN type in the host galaxy. With a redshift of 2.2216, this is the highest redshift SN Ia discovered with a spectroscopic host-galaxy redshift. A further distinguishing feature is that the lensing cluster, at redshift 1.23, is the most distant to date to have an amplified SN. The SN lies in the middle of the color and light-curve shape distributions found at lower redshift, disfavoring strong evolution to z = 2.22. We estimate an amplification due to gravitational lensing of (1.10 0.23 mag) - compatible with the value estimated from the weak-lensing-derived mass and the mass-concentration relation from ΛCDM simulations - making it the most amplified SN Ia discovered behind a galaxy cluster.

Context. Type Ia supernovae (SNe Ia) are widely used to measure the expansion of the Universe. To perform such measurements the luminosity and cosmological redshift (z) of the SNe Ia have to be determined. The uncertainty on z includes an unknown peculiar velocity, which can be very large for SNe Ia in the virialized cores of massive clusters. Aims. We determine which SNe Ia exploded in galaxy clusters using 145 SNe Ia from the Nearby Supernova Factory. We then study how the correction for peculiar velocities of host galaxies inside the clusters improves the Hubble residuals. Methods. We found 11 candidates for membership in clusters. We applied the biweight technique to estimate the redshift of a cluster. Then, we used the galaxy cluster redshift instead of the host galaxy redshift to construct the Hubble diagram. Results. For SNe Ia inside galaxy clusters, the dispersion around the Hubble diagram when peculiar velocities are taken into account is smaller compared with a case without peculiar velocity correction, which has a wRMS = 0.130 ± 0.038 mag instead of wRMS = 0.137 ± 0.036 mag. The significance of this improvement is 3.58σ. If we remove the very nearby Virgo cluster member SN2006X (z < 0.01) from the analysis, the significance decreases to 1.34σ. The peculiar velocity correction is found to be highest for the SNe Ia hosted by blue spiral galaxies. Those SNe Ia have high local specific star formation rates and smaller stellar masses, which is seemingly counter to what might be expected given the heavy concentration of old, massive elliptical galaxies in clusters. Conclusions. As expected, the Hubble residuals of SNe Ia associated with massive galaxy clusters improve when the cluster redshift is taken as the cosmological redshift of the supernova. This fact has to be taken into account in future cosmological analyses in order to achieve higher accuracy for cosmological redshift measurements. We provide an approach to do so.

Context. Observations of type Ia supernovae (SNe Ia) can be used to derive accurate cosmological distances through empirical standardization techniques. Despite this success neither the progenitors of SNe Ia nor the explosion process are fully understood. The U-band region has been less well observed for nearby SNe, due to technical challenges, but is the most readily accessible band for high-redshift SNe. Aims. Using spectrophotometry from the Nearby Supernova Factory, we study the origin and extent of U-band spectroscopic variations in SNe Ia and explore consequences for their standardization and the potential for providing new insights into the explosion process. Methods. We divide the U-band spectrum into four wavelength regions λ(uNi), λ(uTi), λ(uSi) and λ(uCa). Two of these span the Ca H&K λλ 3934, 3969 complex. We employ spectral synthesis using SYNAPPS to associate the two bluer regions with Ni/Co and Ti. Results. The flux of the uTi feature is an extremely sensitive temperature/luminosity indicator, standardizing the SN peak luminosity to 0.116 ± 0.011 mag root mean square (RMS). A traditional SALT2.4 fit on the same sample yields a 0.135 mag RMS. Standardization using uTi also reduces the difference in corrected magnitude between SNe originating from different host galaxy environments. Early U-band spectra can be used to probe the Ni+Co distribution in the ejecta, thus offering a rare window into the source of light curve power. The uCa flux further improves standardization, yielding a 0.086 ± 0.010 mag RMS without the need to include an additional intrinsic dispersion to reach χ2/dof ∼ 1. This reduction in RMS is partially driven by an improved standardization of Shallow Silicon and 91T-like SNe.

Context. Observations of type Ia supernovae (SNe Ia) can be used to derive accurate cosmological distances through empirical standardization techniques. Despite this success neither the progenitors of SNe Ia nor the explosion process are fully understood. The U-band region has been less well observed for nearby SNe, due to technical challenges, but is the most readily accessible band for high-redshift SNe. Aims. Using spectrophotometry from the Nearby Supernova Factory, we study the origin and extent of U-band spectroscopic variations in SNe Ia and explore consequences for their standardization and the potential for providing new insights into the explosion process. Methods. We divide the U-band spectrum into four wavelength regions λ(uNi), λ(uTi), λ(uSi) and λ(uCa). Two of these span the Ca H&K λλ 3934, 3969 complex. We employ spectral synthesis using SYNAPPS to associate the two bluer regions with Ni/Co and Ti. Results. The flux of the uTi feature is an extremely sensitive temperature/luminosity indicator, standardizing the SN peak luminosity to 0.116 ± 0.011 mag root mean square (RMS). A traditional SALT2.4 fit on the same sample yields a 0.135 mag RMS. Standardization using uTi also reduces the difference in corrected magnitude between SNe originating from different host galaxy environments. Early U-band spectra can be used to probe the Ni+Co distribution in the ejecta, thus offering a rare window into the source of light curve power. The uCa flux further improves standardization, yielding a 0.086 ± 0.010 mag RMS without the need to include an additional intrinsic dispersion to reach χ2/dof ∼ 1. This reduction in RMS is partially driven by an improved standardization of Shallow Silicon and 91T-like SNe.

Context. Type Ia supernovae (SNe Ia) are widely used to measure the expansion of the Universe. To perform such measurements the luminosity and cosmological redshift (z) of the SNe Ia have to be determined. The uncertainty on z includes an unknown peculiar velocity, which can be very large for SNe Ia in the virialized cores of massive clusters. Aims. We determine which SNe Ia exploded in galaxy clusters using 145 SNe Ia from the Nearby Supernova Factory. We then study how the correction for peculiar velocities of host galaxies inside the clusters improves the Hubble residuals. Methods. We found 11 candidates for membership in clusters. We applied the biweight technique to estimate the redshift of a cluster. Then, we used the galaxy cluster redshift instead of the host galaxy redshift to construct the Hubble diagram. Results. For SNe Ia inside galaxy clusters, the dispersion around the Hubble diagram when peculiar velocities are taken into account is smaller compared with a case without peculiar velocity correction, which has a wRMS = 0.130 ± 0.038 mag instead of wRMS = 0.137 ± 0.036 mag. The significance of this improvement is 3.58σ. If we remove the very nearby Virgo cluster member SN2006X (z < 0.01) from the analysis, the significance decreases to 1.34σ. The peculiar velocity correction is found to be highest for the SNe Ia hosted by blue spiral galaxies. Those SNe Ia have high local specific star formation rates and smaller stellar masses, which is seemingly counter to what might be expected given the heavy concentration of old, massive elliptical galaxies in clusters. Conclusions. As expected, the Hubble residuals of SNe Ia associated with massive galaxy clusters improve when the cluster redshift is taken as the cosmological redshift of the supernova. This fact has to be taken into account in future cosmological analyses in order to achieve higher accuracy for cosmological redshift measurements. We provide an approach to do so.

Aims. The scientific yield of current and future optical surveys is increasingly limited by systematic uncertainties in the flux calibration. This is the case for type Ia supernova (SN Ia) cosmology programs, where an improved calibration directly translates into improved cosmological constraints. Current methodology rests on models of stars. Here we aim to obtain flux calibration that is traceable to state-of-the-art detector-based calibration. Methods. We present the SNIFS Calibration Apparatus (SCALA), a color (relative) flux calibration system developed for the SuperNova integral field spectrograph (SNIFS), operating at the University of Hawaii 2.2 m (UH 88) telescope. Results. By comparing the color trend of the illumination generated by SCALA during two commissioning runs, and to previous laboratory measurements, we show that we can determine the light emitted by SCALA with a long-term repeatability better than 1%. We describe the calibration procedure necessary to control for system aging. We present measurements of the SNIFS throughput as estimated by SCALA observations. Conclusions. The SCALA calibration unit is now fully deployed at the UH 88 telescope, and with it color-calibration between 4000 Å and 9000 Å is stable at the percent level over a one-year baseline.

Aims. The scientific yield of current and future optical surveys is increasingly limited by systematic uncertainties in the flux calibration. This is the case for type Ia supernova (SN Ia) cosmology programs, where an improved calibration directly translates into improved cosmological constraints. Current methodology rests on models of stars. Here we aim to obtain flux calibration that is traceable to state-of-the-art detector-based calibration. Methods. We present the SNIFS Calibration Apparatus (SCALA), a color (relative) flux calibration system developed for the SuperNova integral field spectrograph (SNIFS), operating at the University of Hawaii 2.2 m (UH 88) telescope. Results. By comparing the color trend of the illumination generated by SCALA during two commissioning runs, and to previous laboratory measurements, we show that we can determine the light emitted by SCALA with a long-term repeatability better than 1%. We describe the calibration procedure necessary to control for system aging. We present measurements of the SNIFS throughput as estimated by SCALA observations. Conclusions. The SCALA calibration unit is now fully deployed at the UH 88 telescope, and with it color-calibration between 4000 Å and 9000 Å is stable at the percent level over a one-year baseline.

Correcting Type Ia Supernova brightnesses for extinction by dust has proven to be a vexing problem. Here we study the dust foreground to the highly reddened SN 2012cu, which is projected onto a dust lane in the galaxy NGC 4772. The analysis is based on multi-epoch, spectrophotometric observations spanning from 3300-9200 Å, obtained by the Nearby Supernova Factory. Phase-matched comparison of the spectroscopically twinned SN 2012cu and SN 2011fe across 10 epochs results in the best-fit color excess of (E (B - V), RMS) = (1.00, 0.03) and total-to-selective extinction ratio of (RV, RMS) = (2.95, 0.08) toward SN 2012cu within its host galaxy. We further identify several diffuse interstellar bands and compare the 5780 Å band with the dust-to-band ratio for the Milky Way (MW). Overall, we find the foreground dust-extinction properties for SN 2012cu to be consistent with those of the MW. Furthermore, we find no evidence for significant time variation in any of these extinction tracers. We also compare the dust extinction curve models of Cardelli et al., ODonnell, and Fitzpatrick, and find the predictions of Fitzpatrick fit SN 2012cu the best. Finally, the distance to NGC4772, the host of SN 2012cu, at a redshift of z = 0.0035, often assigned to the Virgo Southern Extension, is determined to be 16.6 ± 1.1 Mpc. We compare this result with distance measurements in the literature.

Correcting Type Ia Supernova brightnesses for extinction by dust has proven to be a vexing problem. Here we study the dust foreground to the highly reddened SN 2012cu, which is projected onto a dust lane in the galaxy NGC 4772. The analysis is based on multi-epoch, spectrophotometric observations spanning from 3300-9200 Å, obtained by the Nearby Supernova Factory. Phase-matched comparison of the spectroscopically twinned SN 2012cu and SN 2011fe across 10 epochs results in the best-fit color excess of (E (B - V), RMS) = (1.00, 0.03) and total-to-selective extinction ratio of (RV, RMS) = (2.95, 0.08) toward SN 2012cu within its host galaxy. We further identify several diffuse interstellar bands and compare the 5780 Å band with the dust-to-band ratio for the Milky Way (MW). Overall, we find the foreground dust-extinction properties for SN 2012cu to be consistent with those of the MW. Furthermore, we find no evidence for significant time variation in any of these extinction tracers. We also compare the dust extinction curve models of Cardelli et al., ODonnell, and Fitzpatrick, and find the predictions of Fitzpatrick fit SN 2012cu the best. Finally, the distance to NGC4772, the host of SN 2012cu, at a redshift of z = 0.0035, often assigned to the Virgo Southern Extension, is determined to be 16.6 ± 1.1 Mpc. We compare this result with distance measurements in the literature.

We introduce a method for identifying "twin" Type Ia supernovae (SNe Ia) and using them to improve distance measurements. This novel approach to SN Ia standardization is madeossible by spectrophotometric time series observations from the Nearby Supernova Factory (SNfactory). We begin with a well-measured set of SNe, findairs whose spectra match well across the entire optical window, and then test whether this leads to a smaller dispersion in their absolute brightnesses. This analysis is completed in a blinded fashion, ensuring that decisions made in implementing the method do not inadvertently bias the result. We find thatairs of SNe with more closely matched spectra indeed have reduced brightness dispersion. We are able to standardize this initial set of SNfactory SNe to 0.083 0.012 mag, implying a dispersion of 0.072 0.010 mag in the absence ofeculiar velocities. We estimate that with larger numbers of comparison SNe, e.g., using the final SNfactory spectrophotometric data set as a reference, this method will be capable of standardizing high-redshift SNe to within 0.06-0.07 mag. These results imply that at least 3/4 of the variance in Hubble residuals in current SN cosmology analyses is due toreviously unaccounted-for astrophysical differences among the SNe.

We introduce a method for identifying "twin" Type Ia supernovae (SNe Ia) and using them to improve distance measurements. This novel approach to SN Ia standardization is madeossible by spectrophotometric time series observations from the Nearby Supernova Factory (SNfactory). We begin with a well-measured set of SNe, findairs whose spectra match well across the entire optical window, and then test whether this leads to a smaller dispersion in their absolute brightnesses. This analysis is completed in a blinded fashion, ensuring that decisions made in implementing the method do not inadvertently bias the result. We find thatairs of SNe with more closely matched spectra indeed have reduced brightness dispersion. We are able to standardize this initial set of SNfactory SNe to 0.083 0.012 mag, implying a dispersion of 0.072 0.010 mag in the absence ofeculiar velocities. We estimate that with larger numbers of comparison SNe, e.g., using the final SNfactory spectrophotometric data set as a reference, this method will be capable of standardizing high-redshift SNe to within 0.06-0.07 mag. These results imply that at least 3/4 of the variance in Hubble residuals in current SN cosmology analyses is due toreviously unaccounted-for astrophysical differences among the SNe.

We estimate systematic errors due to K-corrections in standard photometric analyses of high-redshift Type Iasupernovae. Errors due to K-correction occur when the spectral template model underlying the light curve fitterpoorly represents the actual supernova spectral energy distribution, meaning that the distance modulus cannot berecovered accurately. In order to quantify this effect, synthetic photometry is performed on artificially redshiftedspectrophotometric data from 119 low-redshift supernovae from the Nearby Supernova Factory, and the resultinglight curves are fit with a conventional light curve fitter. We measure the variation in the standardized magnitudethat would be fit for a given supernova if located at a range of redshifts and observed with various filter setscorresponding to current and future supernova surveys. We find significant variation in the measurements of thesame supernovae placed at different redshifts regardless of filters used, which causes dispersion greater than?0.05 mag for measurements of photometry using the Sloan-like filters and a bias that corresponds to a 0.03 shift inw when applied to an outside data set. To test the result of a shift in supernova population or environment at higherredshifts, we repeat our calculations with the addition of a reweighting of the supernovae as a function of redshiftand find that this strongly affects the results and would have repercussions for cosmology. We discuss possiblemethods to reduce the contribution of the K-correction bias and uncertainty.

We estimate systematic errors due to K-corrections in standard photometric analyses of high-redshift Type Iasupernovae. Errors due to K-correction occur when the spectral template model underlying the light curve fitterpoorly represents the actual supernova spectral energy distribution, meaning that the distance modulus cannot berecovered accurately. In order to quantify this effect, synthetic photometry is performed on artificially redshiftedspectrophotometric data from 119 low-redshift supernovae from the Nearby Supernova Factory, and the resultinglight curves are fit with a conventional light curve fitter. We measure the variation in the standardized magnitudethat would be fit for a given supernova if located at a range of redshifts and observed with various filter setscorresponding to current and future supernova surveys. We find significant variation in the measurements of thesame supernovae placed at different redshifts regardless of filters used, which causes dispersion greater than?0.05 mag for measurements of photometry using the Sloan-like filters and a bias that corresponds to a 0.03 shift inw when applied to an outside data set. To test the result of a shift in supernova population or environment at higherredshifts, we repeat our calculations with the addition of a reweighting of the supernovae as a function of redshiftand find that this strongly affects the results and would have repercussions for cosmology. We discuss possiblemethods to reduce the contribution of the K-correction bias and uncertainty.

We present a sample of normal Type Ia supernovae (SNe Ia) from the Nearby Supernova Factory data set with spectrophotometry at sufficiently late phases to estimate the ejected mass using the bolometric light curve.Wemeasure 56 Ni masses from the peak bolometric luminosity, then compare the luminosity in the 56Co-decay tail to the expected rate of radioactive energy release from ejecta of a given mass. We infer the ejected mass in a Bayesian context using a semi-analytic model of the ejecta, incorporating constraints from contemporary numerical models as priors on the density structure and distribution of 56Ni throughout the ejecta. We find a strong correlation between ejected mass and light-curve decline rate, and consequently 56Ni mass, with ejected masses in our data ranging from 0.9 to 1.4 M⊙. Most fast-declining (SALT2 x1 < -1) normal SNe Ia have significantly sub-Chandrasekhar ejected masses in our fiducial analysis. © 2014 The Authors Published by Oxford University Press on behalf of the Royal Astronomical Society.