Dispersion in the Hubble-Lema\^{i}tre constant measurements from gravitational clustering

TL;DR

This study uses N-body simulations to quantify how peculiar velocities from gravitational clustering affect local measurements of the Hubble-Lemaître constant, highlighting scale-dependent uncertainties and potential biases.

Contribution

It provides a detailed analysis of the impact of peculiar velocities on $H_0$ estimates, emphasizing the significance of scale and local density environments.

Findings

Dispersion in $H_0$ estimates decreases with scale, reaching statistical error levels beyond 135-220 Mpc/h.

Local over-density correlates negatively with deviations in $H_0$, affecting measurements up to 40 Mpc/h.

Errors from peculiar velocities can cause deviations of several percent in local $H_0$ measurements.

Abstract

Measurements of the Hubble-Lema\^{i}tre constant ($H_0$) require us to estimate the distance and recession velocity of galaxies independently. Gravitational clustering that leads to the formation of galaxies and the large scale structure leaves its imprints in the form of peculiar velocities of galaxies. In general, it is not possible to disentangle the peculiar velocity component from the recession velocities of galaxies, and this introduces an uncertainty in the determination of $H_0$. Using N-body simulations, we quantify the impact of peculiar velocities on the $H_0$ estimation. We consider observers to be located in dark matter halos and compute the distribution of the estimated value of $H_0$ across all such observers. We find that the dispersion of this distribution is large at small scales, and it diminishes as we go to large separations, reaching the level of the quoted statistical error in Planck and SH0ES measurements well beyond $\sim$135 Mpc/h and $\sim$220 Mpc/h, respectively. Measurements at smaller scales are susceptible to errors arising from peculiar motions, and this error can propagate to measurements at larger scales in the distance ladder. Notably, we observe a negative correlation between the local over-density around an observer and the deviation of the local and the global value of $H_0$. We show that deviations more significant than 5% of the global values can be encountered frequently at scales of up to 40 Mpc/h, and this is considerably larger than the statistical errors on local estimates. We also analyse the cumulative effect of such errors on mock measurements of $H_0$ as measured from Milky Way-sized halos. We find that this error is sensitive to the lowest distance at which we use measurements. The distribution of $H_0$ in mock measurements has a large tail, and deviations of a few percent from the global value cannot be ruled out.