Nonparametric late-time expansion history reconstruction and implications for the Hubble tension in light of recent DESI and type Ia supernovae data
Jun-Qian Jiang, Davide Pedrotti, S. Santos da Costa, Sunny Vagnozzi
Abstract
We nonparametrically reconstruct the late-time expansion history in light of the latest baryon acoustic oscillation (BAO) measurements from DESI combined with various type Ia supernovae (SNeIa) catalogs, using interpolation through piece-wise natural cubic splines, and a reconstruction procedure based on Gaussian processes (GPs). Applied to DESI BAO and PantheonPlus SNeIa data, both methods indicate that deviations from a reference $\mathrm{\ensuremath{\Lambda}}\mathrm{CDM}$ model in the $z\ensuremath{\lesssim}2$ unnormalized expansion rate $E(z)$ are constrained to be $\ensuremath{\lesssim}10%$, but also consistently identify two features in $E(z)$: a bump at $z\ensuremath{\sim}0.5$, and a depression at $z\ensuremath{\sim}0.9$, which cannot be simultaneously captured by a ${w}_{0}{w}_{a}\mathrm{CDM}$ fit. These features, which are stable against assumptions regarding spatial curvature, interpolation knots, and GP kernel, disappear if one adopts the older SDSS BAO measurements in place of DESI, and decrease in significance when replacing the PantheonPlus catalog with the Union3 and DESY5 ones. We infer $c/({r}_{d}{H}_{0})=29.90\ifmmode\pm\else\textpm\fi{}0.33$, with ${r}_{d}$ the sound horizon at baryon drag and ${H}_{0}$ the Hubble constant. Breaking the ${r}_{d}\ensuremath{-}{H}_{0}$ degeneracy with the SH0ES prior on ${H}_{0}$, the significance of the tension between our nonparametric determination of ${r}_{d}=136.2{0}_{\ensuremath{-}2.40}^{+2.20}\text{ }\text{ }\mathrm{Mpc}$ and the Planck $\mathrm{\ensuremath{\Lambda}}\mathrm{CDM}$-based determination is at the $5\ensuremath{\sigma}$ level, slightly lower than the $6\ensuremath{\sigma}$ obtained when adopting the older SDSS dataset in place of DESI. This indicates the persistence at very high significance of the ``sound horizon tension,'' reinforcing the need for prerecombination new physics. If substantiated in forthcoming data releases, our results tentatively point to oscillatory/nonmonotonic features in the shape of the expansion rate at $z\ensuremath{\lesssim}2$, of potential interest for dark energy model-building.