Within the Friedmann–Lemaître–Robertson– Walker (FLRW) framework, the Hubble constant H0 is an integration constant. Thus, consistency of the model demands observational constancy of H0. We demonstrate redshift evolution of best fit CDM parameters (H0, m) in Pantheon+ supernove (SNe). Redshift evolution of best fit cosmological parameters is a prerequisite to finding a statistically significant evolution as well as identifying alternative models that are competitive with CDM in a Bayesian model comparison. To assess statistical significance, we employ three different methods: (i) Bayesian model comparison, (ii) mock simulations and (iii) profile distributions. The first shows a marginal preference for the vanilla CDM model over an ad hoc model with 3 additional parameters and an unphysical jump in cosmological parameters at z = 1. From mock simulations, we estimate the statistical significance of redshift evolution of best fit parameters and negative dark energy density (m > 1) to be in the 1 − 2σ range, depending on the criteria employed. Importantly, in direct comparison to the same analysis with the earlier Pantheon sample we find that statistical significance of redshift evolution of best fit parameters has increased, as expected for a physical effect. Our profile distribution analysis demonstrates a shift in (H0, m) in excess of 95% confidence level for SNe with redshifts z > 1 and also shows that a degeneracy in MCMC posteriors is not equivalent to a curve of constant χ2. Our findings can be interpreted as a statistical fluctuation or unexplored systematics in Pantheon+ or CDM model breakdown. The first two possibilities are disfavoured by similar trends in independent probes.
