Improved Treatment of Host-galaxy Correlations in Cosmological Analyses with Type Ia Supernovae

We apply BayeSN, our new hierarchical Bayesian model for the SEDs of Type Ia supernovae (SNe Ia), to analyse the griz light curves of 157 nearby SNe Ia (0.015 < z < 0.08) from the public Foundation DR1 dataset. We train a new version of BayeSN, continuous from 0.35–0.95 μm, which we use to model the properties of SNe Ia in the rest-frame z-band, study the properties of dust in their host galaxies, and construct a Hubble diagram of SN Ia distances determined from full griz light curves. Our griz Hubble diagram has a low total RMS of 0.13 mag using BayeSN, compared to 0.16 mag using SALT2. Additionally, we test the consistency of the dust law RV between low- and high-mass host galaxies by using our model to fit the full time- and wavelength-dependent SEDs of SNe Ia up to moderate reddening (peak apparent B − V ≲ 0.3). Splitting the population at the median host mass, we find RV = 2.84 ± 0.31 in low-mass hosts, and RV = 2.58 ± 0.23 in high-mass hosts, both consistent with the global value of RV = 2.61 ± 0.21 that we estimate for the full sample. For all choices of mass split we consider, RV is consistent across the step within ≲ 1.2σ. Modelling population distributions of dust laws in low- and high-mass hosts, we find that both subsamples are highly consistent with the full sample’s population mean μ(RV) = 2.70 ± 0.25 with a 95 per cent upper bound on the population σ(RV) < 0.61. The RV population means are consistent within ≲ 1.2σ. We find that simultaneous fitting of host-mass-dependent dust properties within our hierarchical model does not account for the conventional mass step.

Using Zwicky Transient Facility (ZTF) observations, we identify a pair of ”sibling” Type Ia supernovae (SNe Ia), i.e., hosted by the same galaxy at z = 0.0541. They exploded within 200 days from each other at a separation of $0.6^{\prime \prime}$corresponding to a projected distance of only 0.6 kpc. Performing SALT2 light curve fits to the gri ZTF photometry, we show that for these equally distant ‘standardizable candles”, there is a difference of 2 magnitudes in their rest frame B-band peaks, and the fainter SN has a significantly red SALT2 colour c = 0.57 ± 0.04, while the stretch values x1 of the two SNe are similar, suggesting that the fainter SN is attenuated by dust in the interstellar medium of the host galaxy. We use these measurements to infer the SALT2 colour standardization parameter, β = 3.5 ± 0.3, independent of the underlying cosmology and Malmquist bias. Assuming the colour excess is entirely due to dust, the result differs by 2σ from the average Milky-Way total-to-selective extinction ratio, but is in good agreement with the colour-brightness corrections empirically derived from the most recent SN Ia Hubble-Lemaitre diagram fits. Thus we suggest that SN ”siblings”, which will increasingly be discovered in the coming years, can be used to probe the validity of the colour and lightcurve shape corrections using in SN Ia cosmology while avoiding important systematic effects in their inference from global multi-parameter fits to inhomogeneous data-sets, and also help constrain the role of interstellar dust in SN Ia cosmology.

Current and future cosmological analyses with Type Ia Supernovae (SNe Ia) face three critical challenges: i) measuring redshifts from the supernova or its host galaxy; ii) classifying SNe without spectra; and iii) accounting for correlations between the properties of SNe Ia and their host galaxies. We present here a novel approach that addresses each challenge. In the context of the Dark Energy Survey (DES), we analyze a SNIa sample with host galaxies in the redMaGiC galaxy catalog, a selection of Luminous Red Galaxies. Photo-z estimates for these galaxies are expected to be accurate to σ ∆z/(1+z) ∼ 0.02. The DES-5YR photometrically classified SNIa sample contains approximately 1600 SNe and 125 of these SNe are in redMaGiC galaxies. We demonstrate that redMaGiC galaxies almost exclusively host SNe Ia, reducing concerns with classification uncertainties. With this subsample, we find similar Hubble scatter (to within ∼ 0.01 mag) using photometric redshifts in place of spectroscopic redshifts. With detailed simulations, we show the bias due to using photo-zs from redMaGiC host galaxies on the measurement of the dark energy equation-of-state w is up to ∆w ∼ 0.01 − 0.02. With real data, we measure a difference in w when using redMaGiC photometric redshifts versus spectroscopic redshifts of ∆w = 0.005. Finally, we discuss how SNe in redMaGiC galaxies appear to be a more standardizable population due to a weaker relation between color and luminosity (β) compared to the DES-3YR population by ∼ 5σ; this finding is consistent with predictions that redMaGiC galaxies exhibit lower reddening ratios (R V ) than the general population of SN host galaxies. These results establish the feasibility of performing redMaGiC SN cosmology with photometric survey data in the absence of spectroscopic data.