Sample Variance for Supernovae Distance Measurements and the Hubble tension

TL;DR

This paper evaluates the impact of sample variance on local supernovae measurements of the Hubble constant, finding it contributes minimally to the observed tension with CMB-based estimates, and explores hemispherical asymmetries.

Contribution

It introduces and compares four methods for calculating sample variance in supernova measurements, confirming its limited effect on H0 and analyzing hemispherical asymmetries.

Findings

Sample variance can only cause fluctuations of about ±1 km s$^{-1}$Mpc$^{-1}$ in H0.

Different methods for estimating sample variance produce consistent results.

Hemispherical analysis shows a ~4 km s$^{-1}$Mpc$^{-1}$ higher H0 in the CMB dipole direction.

Abstract

Recent local measurements of the Hubble constant made using supernovae have delivered a value that differs by $\sim$5$\sigma$ (statistical error) from predictions using the Cosmic Microwave Background (CMB), or using Baryon Acoustic Oscillations (BAO) and Big-Bang Nucleosynthesis (BBN) constraints, which are themselves consistent. The effective volume covered by the supernovae is small compared to the other probes, and it is therefore interesting to consider whether sample variance (often also called cosmic variance) is a significant contributor to the offset. We consider four ways of calculating the sample variance: (i) perturbation theory applied to the luminosity distance, which is the most common method considered in the literature; (ii) perturbation of cosmological parameters, as is commonly used to alleviate super-sample covariance in sets of N-body simulations; (iii) a new method based on the variance between perturbed spherical top-hat regions; (iv) using numerical N-body simulations. All give consistent results showing that, for the Pantheon supernova sample, sample variance can only lead to fluctuations in $H_0$ of order $\pm1$ km s$^{-1}$Mpc$^{-1}$ or less. While this is not in itself a new result, the agreement between the methods used adds to its robustness. Furthermore, it is instructive to see how the different methods fit together. We also investigate the internal variance of the $H_{0}$ measurement using SH0ES and Pantheon data. By searching for an offset between measurements in opposite hemispheres, we find that the direction coincident with the CMB dipole has a higher $H_{0}$ measurement than the opposite hemisphere by roughly 4 km s$^{-1}$Mpc$^{-1}$. We compare this with a large number of simulations and find that the size of this asymmetry is statistically likely, but the preference of direction may indicate that further calibration is needed.