Abstract
In this work, we investigate the impact of cosmic flows and density perturbations on Hubble constant H0 measurements using non-linear phase–space reconstructions of the Local Universe (LU). In particular, we rely on a set of 25 precise constrained N-body simulations based on Bayesian initial conditions reconstructions of the LU using the Two-Micron Redshift Survey galaxy sample within distances of about 90 h−1 Mpc. These have been randomly extended up to volumes enclosing distances of 360 h−1 Mpc with augmented Lagrangian perturbation theory (750 simulations in total), accounting in this way for gravitational mode coupling from larger scales, correcting for periodic boundary effects, and estimating systematics of missing attractors (σlarge = 134 s−1 km). We report on Local Group (LG) speed reconstructions, which for the first time are compatible with those derived from cosmic microwave background-dipole measurements: |vLG| = 685 ± 137 s−1 km. The direction (l, b) = (260$_{.}^{\circ}$5 ± 13$_{.}^{\circ}$3, 39$_{.}^{\circ}$1 ± 10$_{.}^{\circ}$4) is found to be compatible with the observations after considering the variance of large scales. Considering this effect of large scales, our local bulk flow estimations assuming a Λ cold dark matter model are compatible with the most recent estimates based on velocity data derived from the Tully–Fisher relation. We focus on low-redshift supernova measurements out to 0.01 < z < 0.025, which have been found to disagree with probes at larger distances. Our analysis indicates that there are two effects related to cosmic variance contributing to this tension. The first one is caused by the anisotropic distribution of supernovae, which aligns with the velocity dipole and hence induces a systematic boost in H0. The second one is due to the inhomogeneous matter fluctuations in the LU. In particular, a divergent region surrounding the Virgo Supercluster is responsible for an additional positive bias in H0. Taking these effects into account yields a correction of ΔH0 = -1.76 ± 0.21 s− 1 km Mpc− 1, thereby reducing the tension between local probes and more distant probes. Effectively H0 is lower by about 2 per cent.
The multiwavelength Tully–Fisher relation with spatially resolved H i kinematics
