skip to main content

Title: SALT3-NIR: Taking the Open-source Type Ia Supernova Model to Longer Wavelengths for Next-generation Cosmological Measurements

A large fraction of Type Ia supernova (SN Ia) observations over the next decade will be in the near-infrared (NIR), at wavelengths beyond the reach of the current standard light-curve model for SN Ia cosmology, SALT3 (∼2800–8700 Å central filter wavelength). To harness this new SN Ia sample and reduce future light-curve standardization systematic uncertainties, we train SALT3 at NIR wavelengths (SALT3-NIR) up to 2μm with the open-source model-training softwareSALTshaker, which can easily accommodate future observations. Using simulated data, we show that the training process constrains the NIR model to ∼2%–3% across the phase range (−20 to 50 days). We find that Hubble residual (HR) scatter is smaller using the NIR alone or optical+NIR compared to optical alone, by up to ∼30% depending on filter choice (95% confidence). There is significant correlation between NIR light-curve stretch measurements and luminosity, with stretch and color corrections often improving HR scatter by up to ∼20%. For SN Ia observations expected from the Roman Space Telescope, SALT3-NIR increases the amount of usable data in the SALT framework by ∼20% at redshiftz≲ 0.4 and by ∼50% atz≲ 0.15. The SALT3-NIR model is part of the open-sourceSNCosmoandSNANASN Ia cosmology packages.

more » « less
Award ID(s):
1815935 1911225 1911151
Author(s) / Creator(s):
; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; more » ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; « less
Publisher / Repository:
DOI PREFIX: 10.3847
Date Published:
Journal Name:
The Astrophysical Journal
Medium: X Size: Article No. 11
["Article No. 11"]
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract

    We present the Keck Infrared Transient Survey, a NASA Key Strategic Mission Support program to obtain near-infrared (NIR) spectra of astrophysical transients of all types, and its first data release, consisting of 105 NIR spectra of 50 transients. Such a data set is essential as we enter a new era of IR astronomy with the James Webb Space Telescope (JWST) and the upcoming Nancy Grace Roman Space Telescope (Roman). NIR spectral templates will be essential to search JWST images for stellar explosions of the first stars and to plan an effective Roman SN Ia cosmology survey, both key science objectives for mission success. Between 2022 February and 2023 July, we systematically obtained 274 NIR spectra of 146 astronomical transients, representing a significant increase in the number of available NIR spectra in the literature. Here, we describe the first release of data from the 2022A semester. We systematically observed three samples: a flux-limited sample that includes all transients <17 mag in a red optical band (usually ZTFror ATLASobands); a volume-limited sample including all transients within redshiftz< 0.01 (D≈ 50 Mpc); and an SN Ia sample targeting objects at phases and light-curve parameters that had scant existing NIR data in the literature. The flux-limited sample is 39% complete (60% excluding SNe Ia), while the volume-limited sample is 54% complete and is 79% complete toz= 0.005. Transient classes observed include common Type Ia and core-collapse supernovae, tidal disruption events, luminous red novae, and the newly categorized hydrogen-free/helium-poor interacting Type Icn supernovae. We describe our observing procedures and data reduction usingPypeIt, which requires minimal human interaction to ensure reproducibility.

    more » « less

    Type Ia supernova (SN Ia) cosmology analyses include a luminosity step function in their distance standardization process to account for an observed yet unexplained difference in the post-standardization luminosities of SNe Ia originating from different host galaxy populations [e.g. high-mass ($M \gtrsim 10^{10} \, {\rm M}_{\odot }$) versus low-mass galaxies]. We present a novel method for including host-mass correlations in the SALT3 (Spectral Adaptive Light curve Template 3) light curve model used for standardizing SN Ia distances. We split the SALT3 training sample according to host-mass, training independent models for the low- and high-host-mass samples. Our models indicate that there are different average Si ii spectral feature strengths between the two populations, and that the average spectral energy distribution of SNe from low-mass galaxies is bluer than the high-mass counterpart. We then use our trained models to perform an SN cosmology analysis on the 3-yr spectroscopically confirmed Dark Energy Survey SN sample, treating SNe from low- and high-mass host galaxies as separate populations throughout. We find that our mass-split models reduce the Hubble residual scatter in the sample, albeit at a low statistical significance. We do find a reduction in the mass-correlated luminosity step but conclude that this arises from the model-dependent re-definition of the fiducial SN absolute magnitude rather than the models themselves. Our results stress the importance of adopting a standard definition of the SN parameters (x0, x1, c) in order to extract the most value out of the light curve modelling tools that are currently available and to correctly interpret results that are fit with different models.

    more » « less

    For the past decade, SALT2 has been the most common model used to fit Type Ia supernova (SN Ia) light curves for dark energy analyses. Recently, the SALT3 model was released, which upgraded a number of model features but has not yet been used for measurements of dark energy. Here, we evaluate the impact of switching from SALT2 to SALT3 for a SN cosmology analysis. We train SALT2 and SALT3 on an identical training sample of 1083 well-calibrated Type Ia supernovae, ensuring that any differences found come from the underlying model framework. We publicly release the results of this training (the SALT ‘surfaces’). We then run a cosmology analysis on the public Dark Energy Survey 3-Yr Supernova data sample (DES-SN3YR), and on realistic simulations of those data. We provide the first estimate of the SN + CMB systematic uncertainty arising from the choice of SALT model framework (i.e. SALT2 versus SALT3), Δw  = + 0.001 ± 0.005 – a negligible effect at the current level of dark energy analyses. We also find that the updated surfaces are less sensitive to photometric calibration uncertainties than previous SALT2 surfaces, with the average spectral energy density dispersion reduced by a factor of two over optical wavelengths. This offers an opportunity to reduce the contribution of calibration errors to SN cosmology uncertainty budgets.

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

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