We present constraints on cosmological parameters from the Pantheon+ analysis of 1701 light curves of 1550 distinct Type Ia supernovae (SNe Ia) ranging in redshift from z = 0.001 to 2.26. This work features an increased sample size from the addition of multiple cross-calibrated photometric systems of SNe covering an increased redshift span, and improved treatments of systematic uncertainties in comparison to the original Pantheon analysis, which together result in a factor of 2 improvement in cosmological constraining power. For a flat ΛCDM model, we find Ω M = 0.334 ± 0.018 from SNe Ia alone. For a flat w 0CDM model, we measure w 0 = −0.90 ± 0.14 from SNe Ia alone, H 0 = 73.5 ± 1.1 km s−1 Mpc−1 when including the Cepheid host distances and covariance (SH0ES), and w 0 = − 0.978 − 0.031 + 0.024 when combining the SN likelihood with Planck constraints from the cosmic microwave background (CMB) and baryon acoustic oscillations (BAO); both w 0 values are consistent with a cosmological constant. We also present the most precise measurements to date on the evolution of dark energy in a flat w 0 w a CDM universe, and measure w a = − 0.1 − 2.0 + 0.9 from Pantheon+ SNe Ia alone, H 0 = 73.3 ± 1.1 km s−1 Mpc−1 when including SH0ES Cepheid distances, and w a = − 0.65 − 0.32 + 0.28 when combining Pantheon+ SNe Ia with CMB and BAO data. Finally, we find that systematic uncertainties in the use of SNe Ia along the distance ladder comprise less than one-third of the total uncertainty in the measurement of H 0 and cannot explain the present “Hubble tension” between local measurements and early universe predictions from the cosmological model.

Here we present 1701 light curves of 1550 unique, spectroscopically confirmed Type Ia supernovae (SNe Ia) that will be used to infer cosmological parameters as part of the Pantheon+ SN analysis and the Supernovae and H 0 for the Equation of State of dark energy distance-ladder analysis. This effort is one part of a series of works that perform an extensive review of redshifts, peculiar velocities, photometric calibration, and intrinsic-scatter models of SNe Ia. The total number of light curves, which are compiled across 18 different surveys, is a significant increase from the first Pantheon analysis (1048 SNe), particularly at low redshift (z). Furthermore, unlike in the Pantheon analysis, we include light curves for SNe with z < 0.01 such that SN systematic covariance can be included in a joint measurement of the Hubble constant (H 0) and the dark energy equation-of-state parameter (w). We use the large sample to compare properties of 151 SNe Ia observed by multiple surveys and 12 pairs/triplets of “SN siblings”—SNe found in the same host galaxy. Distance measurements, application of bias corrections, and inference of cosmological parameters are discussed in the companion paper by Brout et al., and the determination of H 0 is discussed by Riess et al. These analyses will measure w with ∼3% precision and H 0 with ∼1 km s−1 Mpc−1 precision.

Improving the use of Type Ia supernovae (SNe Ia) as standard candles requires a better approach to incorporate the relationship between SNe Ia and the properties of their host galaxies. Using a spectroscopically confirmed sample of ∼1600 SNe Ia, we develop the first empirical model of underlying populations for SNe Ia light-curve properties that includes their dependence on host-galaxy stellar mass; we find a significant correlation between stretch population and stellar mass (99.9% confidence) and a weaker correlation between color and stellar mass (90% confidence). These populations are important inputs to simulations that are used to model selection effects and correct distance biases within the BEAMS with Bias Correction (BBC) framework. Here we improve BBC to also account for SNe Ia-host correlations, and we validate this technique on simulated data samples. We recover the input relationship between SNe Ia luminosity and host-galaxy stellar mass (the mass step, γ) with a bias of 0.004 ±0.001 mag, which is a factor of 5 improvement over previous methods that have a γ bias of ∼0.02 ± 0.001 mag. We adapt BBC for a novel dust-based model of intrinsic brightness variations, which results in a greatly reduced mass step for data (γ = 0.017 ± 0.008) and for simulations (γ = 0.006 ± 0.007). Analyzing simulated SNe Ia, the biases on the dark energy equation of state, w, vary from Δw = 0.006(5) to 0.010(5) with our new BBC method; these biases are significantly smaller than the 0.02(5) w bias using previous BBC methods that ignore SNe Ia-host correlations.

We present constraints on cosmological parameters from the Pantheon+ analysis of 1701 light curves of 1550 distinct Type Ia supernovae (SNe Ia) ranging in redshift from z = 0.001 to 2.26. This work features an increased sample size, increased redshift span, and improved treatment of systematic uncertainties in comparison to the original Pantheon analysis and results in a factor of 2 improvement in cosmological constraining power. For a FlatΛCDM model, we find Ω M = 0.338 ± 0.018 from SNe Ia alone. For a Flatw 0 CDM model, we measure w 0 = −0.89 ± 0.13 from SNe Ia alone, H 0 = 72.86 +0.94 −1.06 km s −1 Mpc −1 when including the Cepheid host distances and covariance (SH0ES), and w 0 = −0.978 +0.024 −0.031 when combining the SN likelihood with constraints from the cosmic microwave background (CMB) and baryon acoustic oscillations (BAO); both w 0 values are consistent with a cosmological constant. We also present the most precise measurements to date on the evolution of dark energy in a Flatw 0 w a CDM universe, and measure w a = −0.4 +1.0 −1.8 from Pantheon+ alone, H 0 = 73.40 +0.99 −1.22 km s −1 Mpc −1 when including SH0ES, and w a = −0.65 +0.28 −0.32 when combining Pan-theon+ with CMB and BAO data. Finally, we find that systematic uncertainties in the use of SNe Ia along the distance ladder comprise less than one third of the total uncertainty in the measurement of H 0 and cannot explain the present "Hubble tension" between local measurements and early-Universe predictions from the cosmological model.

Here we present 1701 light curves of spectroscopically confirmed Type Ia supernovae (SNe Ia) that will be used to infer cosmological parameters as part of the Pantheon+ SN analysis and the SH0ES (Supernovae and H 0 for the Equation of State of dark energy) distance-ladder analysis. This effort is one part of a series of works that perform an extensive review of redshifts, peculiar velocities, photometric calibration, and intrinsic-scatter models of SNe Ia. The total number of light curves, which are compiled across 18 different surveys, is a significant increase from the first Pantheon analysis (1048 SNe), particularly at low redshift (z). Furthermore, unlike in the Pantheon analysis, we include light curves for SNe with z < 0.01 such that SN systematic covariance can be included in a joint measurement of the Hubble constant (H 0 ) and the dark energy equation-of-state parameter (w). We use the large sample to compare properties of 170 SNe Ia observed by multiple surveys and 12 pairs/triplets of "SN siblings" -SNe found in the same host galaxy. Distance measurements, application of bias corrections, and inference of cosmological parameters are discussed in the companion paper by Brout et al. (2022b), and the determination of H 0 is discussed by Riess et al. (2022). These analyses will measure w with ∼ 3% precision and H 0 with ∼ 1 km/s/Mpc precision.

Recent studies have shown that the observed colour distributions of Type Ia SNe (SNIa) are welldescribed by a combination of distributions from dust and intrinsic colour. Here we present a new forward-modeling fitting method (Dust2Dust) to measure the parent dust and colour distributions, including their dependence on host-galaxy mass. At each fit step, the SNIa selection efficiency is determined from a large simulated sample that is re-weighted to reflect the proposed distributions. We use five separate metrics to constrain the Dust2Dust parameters: distribution of fitted lightcurve colour c, cosmological residual trends with c, cosmological residual scatter with c, fitted colourluminosity relationship β SALT2 , and intrinsic scatter σ int . Using the Pantheon+ data sample, we present results for a Dust2Dust fit that includes 4 parameters describing intrinsic colour variations and 8 parameters describing dust. Furthermore, we propagate the Dust2Dust parameter uncertainties and covariance to the dark energy equation-of-state w and Hubble constant H 0 : we find σw = 0.005 and σH 0 = 0.145 km/s/Mpc. The Dust2Dust code is publically available.

While many studies have shown a correlation between properties of the light curves of SNe Ia and properties of their host galaxies, it remains unclear what is driving these correlations. We introduce a new direct method to study these correlations by analyzing "parent" galaxies that host multiple SNe Ia "siblings." Here, we search the Dark Energy Survey SN sample, one of the largest samples of discovered SNe, and find eight galaxies that hosted two likely SNe Ia. Comparing the light-curve properties of these SNe and recovered distances from the light curves, we find no better agreement between properties of SNe in the same galaxy as any random pair of galaxies, with the exception of the SN light-curve stretch. We show at 2.8σ significance that at least one-half of the intrinsic scatter of SNe Ia distance modulus residuals is not from common host properties. We also discuss the robustness with which we could make this evaluation with LSST, which will find 100× more pairs of galaxies, and pave a new line of study on the consistency of SNe Ia in the same parent galaxies. Finally, we argue that it is unlikely that some of these SNe are actually single, lensed SN with multiple images.

We use a sample of 809 photometrically classified type Ia supernovae (SNe Ia) discovered by the Dark Energy Survey (DES) along with 40415 field galaxies to calculate the rate of SNe Ia per galaxy in the redshift range 0.2 &lt; z &lt; 0.6. We recover the known correlation between SN Ia rate and galaxy stellar mass across a broad range of scales 8.5 ≤ log (M*/M⊙) ≤ 11.25. We find that the SN Ia rate increases with stellar mass as a power-law with index 0.63 ± 0.02, which is consistent with previous work. We use an empirical model of stellar mass assembly to estimate the average star-formation histories (SFHs) of galaxies across the stellar mass range of our measurement. Combining the modelled SFHs with the SN Ia rates to estimate constraints on the SN Ia delay time distribution (DTD), we find the data are fit well by a power-law DTD with slope index β = −1.13 ± 0.05 and normalisation A = 2.11 ± 0.05 × 10−13SNeM⊙−1yr−1,, which corresponds to an overall SN Ia production efficiency $N_{\mathrm{Ia}}/M_* = 0.9 _{-0.7}^{+4.0} \times 10^{-3} \mathrm{SNe} \mathrm{M}_{\odot }^{-1}$,. Upon splitting the SN sample by properties of the light curves, we find a strong dependence on DTD slope with the SN decline rate, with slower-declining SNe exhibiting a steeper DTD slope. We interpret this as a result of a relationship between intrinsic luminosity and progenitor age, and explore the implications of the result in the context of SN Ia progenitors.