skip to main content


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
NSF-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

    The unknown cause of the correlation between Type Ia supernova (SN Ia) Hubble residuals and their host-galaxy masses (the “mass step”) may bias cosmological parameter measurements. To better understand the mass step, we develop a SALT3 light-curve model for SN cosmology that uses the host-galaxy masses of 296 low-redshift SNe Ia to derive a spectral energy distribution–host-galaxy mass relationship. The resulting model has larger CaiiH and K, Caiinear-infrared triplet, and Siiiequivalent widths for SNe in low-mass host galaxies at 2.2–2.7σsignificance; this indicates higher explosion energies per unit mass in low-mass-hosted SNe. The model has phase-dependent changes in SN Ia colors as a function of host mass, indicating intrinsic differences in mean broadband light curves. Although the model provides a better fit to the SN data overall, it does not substantially reduce data–model residuals for a typical light curve in our sample nor does it significantly reduce Hubble residual dispersion. This is because we find that previous SALT models parameterized most host-galaxy dependencies with their first principal component, although they failed to model some significant spectral variations. Our new model is luminosity and cosmology independent, and applying it to data reduces the mass step by 0.021 ± 0.002 mag (uncertainty accounts for correlated data sets); these results indicate that ∼35% of the mass step can be attributed to luminosity-independent effects. This SALT model version could be trained using alternative host-galaxy properties and at different redshifts, and therefore will be a tool for understanding redshift-dependent correlations between SNe Ia and their host properties as well as their impact on cosmological parameter measurements.

     
    more » « less
  2. 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
  3. ABSTRACT

    While conventional Type Ia supernova (SN Ia) cosmology analyses rely primarily on rest-frame optical light curves to determine distances, SNe Ia are excellent standard candles in near-infrared (NIR) light, which is significantly less sensitive to dust extinction. An SN Ia spectral energy distribution (SED) model capable of fitting rest-frame NIR observations is necessary to fully leverage current and future SN Ia data sets from ground- and space-based telescopes including HST, LSST, JWST, and RST. We construct a hierarchical Bayesian model for SN Ia SEDs, continuous over time and wavelength, from the optical to NIR (B through H, or $0.35{-}1.8\, \mu$m). We model the SED as a combination of physically distinct host galaxy dust and intrinsic spectral components. The distribution of intrinsic SEDs over time and wavelength is modelled with probabilistic functional principal components and the covariance of residual functions. We train the model on a nearby sample of 79 SNe Ia with joint optical and NIR light curves by sampling the global posterior distribution over dust and intrinsic latent variables, SED components and population hyperparameters. Photometric distances of SNe Ia with NIR data near maximum obtain a total RMS error of 0.10 mag with our BayeSN model, compared to 0.13–0.14 mag with SALT2 and SNooPy for the same sample. Jointly fitting the optical and NIR data of the full sample up to moderate reddening (host E(B − V) < 0.4) for a global host dust law, we find RV = 2.9 ± 0.2, consistent with the Milky Way average.

     
    more » « less
  4. 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
  5. 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