skip to main content
US FlagAn official website of the United States government
dot gov icon
Official websites use .gov
A .gov website belongs to an official government organization in the United States.
https lock icon
Secure .gov websites use HTTPS
A lock ( lock ) or https:// means you've safely connected to the .gov website. Share sensitive information only on official, secure websites.


Title: A Probabilistic Autoencoder for Type Ia Supernova Spectral Time Series
Abstract 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.  more » « less
Award ID(s):
1839217
PAR ID:
10392077
Author(s) / Creator(s):
; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; more » ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; « less
Date Published:
Journal Name:
The Astrophysical Journal
Volume:
935
Issue:
1
ISSN:
0004-637X
Page Range / eLocation ID:
5
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. ABSTRACT After correcting for their light-curve shape and colour, Type Ia supernovae (SNe Ia) are precise cosmological distance indicators. However, there remains a non-zero intrinsic scatter in the differences between measured distance and that inferred from a cosmological model (i.e. Hubble residuals or HRs), indicating that SN Ia distances can potentially be further improved. We use the open-source relational data base kaepora to generate composite spectra with desired average properties of phase, light-curve shape, and HR. At many phases, the composite spectra from two subsamples with positive and negative average HRs are significantly different. In particular, in all spectra from 9 d before to 15 d after peak brightness, we find that SNe with negative HRs have, on average, higher ejecta velocities (as seen in nearly every optical spectral feature) than SNe with positive HRs. At +4 d relative to B-band maximum, using a sample of 62 SNe Ia, we measure a 0.091 ± 0.035 mag (2.7σ) HR step between SNe with Si ii λ6355 line velocities ($$v_{Si\, rm{\small II}}$$) higher/lower than −11 000 km s−1 (the median velocity). After light-curve shape and colour correction, SNe with higher velocities tend to have underestimated distance moduli relative to a cosmological model. The intrinsic scatter in our sample reduces from 0.094 to 0.082 mag after making this correction. Using the Si ii λ6355 velocity evolution of 115 SNe Ia, we estimate that a velocity difference >500 km s−1 exists at each epoch between the positive-HR and negative-HR samples with 99.4 per cent confidence. Finally at epochs later than +37 d, we observe that negative-HR composite spectra tend to have weaker spectral features in comparison to positive-HR composite spectra. 
    more » « less
  2. Abstract We investigate the potential of using a sample of very high-redshift (2 ≲z≲ 6) (VHZ) Type Ia supernovae (SNe Ia) attainable by JWST on constraining cosmological parameters. At such high redshifts, the age of the universe is young enough that the VHZ SN Ia sample comprises the very first SNe Ia of the universe, with progenitors among the very first generation of low-mass stars that the universe has made. We show that the VHZ SNe Ia can be used to disentangle systematic effects due to the luminosity distance evolution with redshifts intrinsic to SN Ia standardization. Assuming that the systematic evolution can be described by a linear or logarithmic formula, we found that the coefficients of this dependence can be determined accurately and decoupled from cosmological models. Systematic evolution as large as 0.15 mag and 0.45 mag out toz= 5 can be robustly separated from popular cosmological models for linear and logarithmic evolution, respectively. The VHZ SNe Ia will lay the foundation for quantifying the systematic redshift evolution of SN Ia luminosity distance scales. When combined with SN Ia surveys at comparatively lower redshifts, the VHZ SNe Ia allow for the precise measurement of the history of the expansion of the universe fromz∼ 0 to the epoch approaching reionization. 
    more » « less
  3. ABSTRACT The Type Ia supernovae (SNe Ia) 2011by, hosted in NGC 3972, and 2011fe, hosted in M101, are optical ‘twins,’ having almost identical optical light-curve shapes, colours, and near-maximum-brightness spectra. However, SN 2011fe had significantly more ultraviolet (UV; 1600 < λ < 2500 Å) flux than SN 2011by before and at peak luminosity. Several theoretical models predict that SNe Ia with higher progenitor metallicity should (1) have additional UV opacity and thus lower UV flux; (2) have an essentially unchanged optical spectral-energy distribution; (3) have a similar optical light-curve shape; and (4) because of the excess neutrons, produce more stable Fe-group elements at the expense of radioactive 56Ni and thus have a lower peak luminosity. Following these predictions, Foley and Kirshner suggested that the difference in UV flux between SNe 2011by and 2011fe was the result of their progenitors having significantly different metallicities. They also measured a large, but insignificant, difference between the peak absolute magnitudes of the SNe (ΔMV, peak = 0.60 ± 0.36 mag), with SN 2011fe being more luminous. We present a new Cepheid-based distance to NGC 3972, substantially improving the precision of the distance measurement for SN 2011by. With these new data, we determine that the SNe have significantly different peak luminosities (ΔMV, peak = 0.335 ± 0.069 mag). Consequently, SN 2011fe produced 38 per cent more 56Ni than SN 2011by, consistent with predictions for progenitor metallicity differences for these SNe, although alternative models may also explain this difference. We discuss how progenitor metallicity differences can contribute to the intrinsic scatter for light-curve-shape-corrected SN luminosities, the use of ‘twin’ SNe for measuring distances, and implications for using SNe Ia for constraining cosmological parameters. 
    more » « less
  4. Abstract We present the largest and most homogeneous collection of near-infrared (NIR) spectra of Type Ia supernovae (SNe Ia): 339 spectra of 98 individual SNe obtained as part of the Carnegie Supernova Project-II. These spectra, obtained with the FIRE spectrograph on the 6.5 m Magellan Baade telescope, have a spectral range of 0.8–2.5μm. Using this sample, we explore the NIR spectral diversity of SNe Ia and construct a template of spectral time series as a function of the light-curve-shape parameter, color stretchsBV. Principal component analysis is applied to characterize the diversity of the spectral features and reduce data dimensionality to a smaller subspace. Gaussian process regression is then used to model the subspace dependence on phase and light-curve shape and the associated uncertainty. Our template is able to predict spectral variations that are correlated withsBV, such as the hallmark NIR features: Mgiiat early times and theH-band break after peak. Using this template reduces the systematic uncertainties inK-corrections by ∼90% compared to those from the Hsiao template. These uncertainties, defined as the meanK-correction differences computed with the color-matched template and observed spectra, are on the level of 4 × 10−4mag on average. This template can serve as the baseline spectral energy distribution for light-curve fitters and can identify peculiar spectral features that might point to compelling physics. The results presented here will substantially improve future SN Ia cosmological experiments, for both nearby and distant samples. 
    more » « less
  5. ABSTRACT We present improved photometric measurements for the host galaxies of 206 spectroscopically confirmed type Ia supernovae discovered by the Dark Energy Survey Supernova Program (DES-SN) and used in the first DES-SN cosmological analysis. For the DES-SN sample, when considering a 5D (z, x1, c, α, β) bias correction, we find evidence of a Hubble residual ‘mass step’, where SNe Ia in high-mass galaxies (>1010M⊙) are intrinsically more luminous (after correction) than their low-mass counterparts by $$\gamma =0.040\pm 0.019$$ mag. This value is larger by 0.031 mag than the value found in the first DES-SN cosmological analysis. This difference is due to a combination of updated photometric measurements and improved star formation histories and is not from host-galaxy misidentification. When using a 1D (redshift-only) bias correction the inferred mass step is larger, with $$\gamma =0.066\pm 0.020$$ mag. The 1D−5D γ difference for DES-SN is $$0.026\pm 0.009$$ mag. We show that this difference is due to a strong correlation between host galaxy stellar mass and the x1 component of the 5D distance-bias correction. Including an intrinsic correlation between the observed properties of SNe Ia, stretch and colour, and stellar mass in simulated SN Ia samples, we show that a 5D fit recovers γ with −9 mmag bias compared to a +2 mmag bias for a 1D fit. This difference can explain part of the discrepancy seen in the data. Improvements in modelling correlations between galaxy properties and SN is necessary to ensure unbiased precision estimates of the dark energy equation of state as we enter the era of LSST. 
    more » « less