The trouble with Hubble: Local versus global expansion rates in inhomogeneous cosmological simulations with numerical relativity

In a fully inhomogeneous, anisotropic cosmological simulation performed by solving Einstein’s equations with numerical relativity, we find a local measurement of the effective Hubble parameter differs by less than 1% compared to the global value. This variance is consistent with predictions from Newtonian gravity. We analyse the averaged local expansion rate on scales comparable to Type 1a supernova surveys, and find that local variance cannot resolve the tension between the \citet{riess2018b} and \citet{planck2018a} measurements.