Abstract
We develop a novel approach to constrain the Hubble parameter $H_0$ and theprimordial power spectrum amplitude $A_\mathrm{s}$ using supernovae type Ia(SNIa) data. By considering SNIa as tracers of the peculiar velocity field, wecan model their distance and their covariance as a function of cosmologicalparameters without the need of calibrators like Cepheids; this yields a newindependent probe of the large-scale structure based on SNIa data withoutdistance anchors. Crucially, we implement a differentiable pipeline in JAX,including efficient emulators and affine sampling, reducing inference time fromyears to hours on a single GPU. We first validate our method on mock datasets,demonstrating that we can constrain $H_0$ and $\log 10^{10}A_\mathrm{s}$ within$\sim10\%$ using $\sim10^3$ SNIa. We then test our pipeline with SNIa from an$N$-body simulation, obtaining $7\%$-level unbiased constraints on $H_0$ with amoderate noise level. We finally apply our method to Pantheon+ data,constraining $H_0$ at the $10\%$ level without Cepheids when fixing$A_\mathrm{s}$ to its $\it{Planck}$ value. On the other hand, we obtain$15\%$-level constraints on $\log 10^{10}A_\mathrm{s}$ in agreement with$\it{Planck}$ when including Cepheids in the analysis. In light of upcomingobservations of low redshift SNIa from the Zwicky Transient Facility and theVera Rubin Legacy Survey of Space and Time, surveys for which our method willdevelop its full potential, we make our code publicly available.