Constraining the local variance of H₀ from directional analyses

We evaluate the local variance of the Hubble Constant $H_0$ with low-z Type Ia Supernovae (SNe). Our analyses are performed using a hemispherical comparison method in order to test whether taking the bulk flow motion into account can reconcile the measurement of the Hubble Constant $H_0$ from standard candles ($H_0 = 73.8 \pm 2.4 \; \mathrm{km \; s}^{-1}\; \mathrm{Mpc}^{-1}$) with that of the Planck's Cosmic Microwave Background data ($H_0 = 67.8 \pm 0.9 \; \mathrm{km \; s}^{-1} \mathrm{Mpc}^{-1}$). We obtaina Hubble Constant maximal variance of $\delta H_0 = (2.30 \pm 0.86) \; \mathrm{km \; s}^{-1} \mathrm{Mpc}^{-1}$ towards the $(l,b) = (315^{\circ},27^{\circ})$ direction. Interestingly, this result agrees with the bulk flow direction estimates found in the literature, as well as previous evaluations of the $H_0$ variance due to the presence of nearby inhomogeneities. We assess the statistical significance of this result with different prescriptions of Monte Carlo simulations, obtaining moderate statistical significance, i.e., $68.7$\% confidence level (CL) for such variance. Furthermore, we test the hypothesis of a higher $H_0$ value in the presence of a bulk flow velocity dipole, finding some evidence for this result which, however, cannot be claimed to be significant due to the current large uncertainty in the SNe distance modulus. Then, we conclude that the tension between different $H_0$ determinations can plausibly be caused to the bulk flow motion of the local Universe, even though the current incompleteness of the SNe data set, both in terms of celestial coverage and distance uncertainties, does not allow a high statistical significance for these results or a definitive conclusion about this issue.