Nonparametric late-time expansion history reconstruction and implications for the Hubble tension in light of recent DESI and type Ia supernovae data

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 ΛCDM model in the z 2 unnormalized expansion rate E(z) are constrained to be 10%, but also consistently identify two features in E(z): a bump at z∼0.5, and a depression at z∼0.9, which cannot be simultaneously captured by a w0waCDM 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/(rdH0)=29.90±0.33, with rd the sound horizon at baryon drag and H0 the Hubble constant. Breaking the rd-H0 degeneracy with the SH0ES prior on H0, the significance of the tension between our nonparametric determination of rd=136.20-2.40+2.20 Mpc and the Planck ΛCDM-based determination is at the 5σ level, slightly lower than the 6σ 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 2, of potential interest for dark energy model-building.