skip to main content

Title: A hierarchical Bayesian SED model for Type Ia supernovae in the optical to near-infrared

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 more » 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.

« less
 ;  ;  ;  ;  
Publication Date:
Journal Name:
Monthly Notices of the Royal Astronomical Society
Page Range or eLocation-ID:
p. 3939-3966
Oxford University Press
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract

    Type Ia supernovae (SNe Ia) are more precise standardizable candles when measured in the near-infrared (NIR) than in the optical. With this motivation, from 2012 to 2017 we embarked on the RAISIN program with the Hubble Space Telescope (HST) to obtain rest-frame NIR light curves for a cosmologically distant sample of 37 SNe Ia (0.2 ≲z≲ 0.6) discovered by Pan-STARRS and the Dark Energy Survey. By comparing higher-zHST data with 42 SNe Ia atz< 0.1 observed in the NIR by the Carnegie Supernova Project, we construct a Hubble diagram from NIR observations (with only time of maximum light and some selection cuts from optical photometry) to pursue a unique avenue to constrain the dark energy equation-of-state parameter,w. We analyze the dependence of the full set of Hubble residuals on the SN Ia host galaxy mass and find Hubble residual steps of size ∼0.06-0.1 mag with 1.5σ−2.5σsignificance depending on the method and step location used. Combining our NIR sample with cosmic microwave background constraints, we find 1 +w= −0.17 ± 0.12 (statistical + systematic errors). The largest systematic errors are the redshift-dependent SN selection biases and the properties of the NIR mass step. We also use these data to measureH0=more »75.9 ± 2.2 km s−1Mpc−1from stars with geometric distance calibration in the hosts of eight SNe Ia observed in the NIR versusH0= 71.2 ± 3.8 km s−1Mpc−1using an inverse distance ladder approach tied to Planck. Using optical data, we find 1 +w= −0.10 ± 0.09, and with optical and NIR data combined, we find 1 +w= −0.06 ± 0.07; these shifts of up to ∼0.11 inwcould point to inconsistency in the optical versus NIR SN models. There will be many opportunities to improve this NIR measurement and better understand systematic uncertainties through larger low-zsamples, new light-curve models, calibration improvements, and eventually by building high-zsamples from the Roman Space Telescope.

    « less
  2. Aims . We present a comprehensive dataset of optical and near-infrared photometry and spectroscopy of type Ia supernova (SN) 2016hnk, combined with integral field spectroscopy (IFS) of its host galaxy, MCG -01-06-070, and nearby environment. Our goal with this complete dataset is to understand the nature of this peculiar object. Methods . Properties of the SN local environment are characterized by means of single stellar population synthesis applied to IFS observations taken two years after the SN exploded. We performed detailed analyses of SN photometric data by studying its peculiar light and color curves. SN 2016hnk spectra were compared to other 1991bg-like SNe Ia, 2002es-like SNe Ia, and Ca-rich transients. In addition, we used abundance stratification modeling to identify the various spectral features in the early phase spectral sequence and also compared the dataset to a modified non-LTE model previously produced for the sublumnious SN 1999by. Results . SN 2016hnk is consistent with being a subluminous ( M B  = −16.7 mag, s B V =0.43 ± 0.03), highly reddened object. The IFS of its host galaxy reveals both a significant amount of dust at the SN location, residual star formation, and a high proportion of old stellar populations in themore »local environment compared to other locations in the galaxy, which favors an old progenitor for SN 2016hnk. Inspection of a nebular spectrum obtained one year after maximum contains two narrow emission lines attributed to the forbidden [Ca  II ] λ λ 7291,7324 doublet with a Doppler shift of 700 km s −1 . Based on various observational diagnostics, we argue that the progenitor of SN 2016hnk was likely a near Chandrasekhar-mass ( M Ch ) carbon-oxygen white dwarf that produced 0.108 M ⊙ of 56 Ni. Our modeling suggests that the narrow [Ca  II ] features observed in the nebular spectrum are associated with 48 Ca from electron capture during the explosion, which is expected to occur only in white dwarfs that explode near or at the M Ch limit.« less
  3. 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 thismore »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.« less

    Recent analyses have found intriguing correlations between the colour (c) of type Ia supernovae (SNe Ia) and the size of their ‘mass-step’, the relationship between SN Ia host galaxy stellar mass (Mstellar) and SN Ia Hubble residual, and suggest that the cause of this relationship is dust. Using 675 photometrically classified SNe Ia from the Dark Energy Survey 5-yr sample, we study the differences in Hubble residual for a variety of global host galaxy and local environmental properties for SN Ia subsamples split by their colour. We find a 3σ difference in the mass-step when comparing blue (c < 0) and red (c > 0) SNe. We observe the lowest r.m.s. scatter (∼0.14 mag) in the Hubble residual for blue SNe in low mass/blue environments, suggesting that this is the most homogeneous sample for cosmological analyses. By fitting for c-dependent relationships between Hubble residuals and Mstellar, approximating existing dust models, we remove the mass-step from the data and find tentative ∼2σ residual steps in rest-frame galaxy U − R colour. This indicates that dust modelling based on Mstellar may not fully explain the remaining dispersion in SN Ia luminosity. Instead, accounting for a c-dependent relationship between Hubble residuals and globalmore »U − R, results in ≤1σ residual steps in Mstellar and local U − R, suggesting that U − R provides different information about the environment of SNe Ia compared to Mstellar, and motivating the inclusion of galaxy U − R colour in SN Ia distance bias correction.

    « less
  5. ABSTRACT We present a map of the total intrinsic reddening across ≃34 deg2 of the Small Magellanic Cloud (SMC) derived using optical (ugriz) and near-infrared (IR; YJKs) spectral energy distributions (SEDs) of background galaxies. The reddening map is created using a subsample of 29 274 galaxies with low levels of intrinsic reddening based on the lephare χ2 minimization SED-fitting routine. We find statistically significant enhanced levels of reddening associated with the main body of the SMC compared with regions in the outskirts [ΔE(B − V) ≃ 0.3 mag]. A comparison with literature reddening maps of the SMC shows that, after correcting for differences in the volume of the SMC sampled, there is good agreement between our results and maps created using young stars. In contrast, we find significant discrepancies between our results and maps created using old stars or based on longer wavelength far-IR dust emission that could stem from biased samples in the former and uncertainties in the far-IR emissivity and the optical properties of the dust grains in the latter. This study represents one of the first large-scale categorizations of extragalactic sources behind the SMC and as such we provide the lephare outputs for our full sample of ∼500 000 sources.