Backreaction in Numerical Relativity: Averaging on Newtonian gauge-like hypersurfaces in Einstein Toolkit cosmological simulations

TL;DR

This paper introduces a spatial averaging scheme in cosmological simulations to study backreaction effects, finding negligible overall backreaction but significant local curvature fluctuations and fluid flow impacts.

Contribution

It presents a novel averaging method on Newtonian gauge-like hypersurfaces in Einstein Toolkit simulations, revealing detailed backreaction and curvature fluctuation insights.

Findings

Negligible cosmic backreaction on initial hypersurfaces.

Significant curvature fluctuations up to 10% in small sub-volumes.

Fluid flow causes up to 5% density change in 200 Mpc spheres.

Abstract

We introduce a spatial averaging scheme and use it to study the evolution of spatial averages in large-scale simulations of cosmological structure formation performed with the Einstein Toolkit. The averages are performed on the spatial hypersurfaces of the simulation setup which, at least initially, represent the hypersurfaces of statistical homogeneity and isotropy. We find only negligible cosmic backreaction on these hypersurfaces even on very small scales, but find significant curvature fluctuations of up to $10\%$ in $\Omega_R$ for sub-volumes with radius $\sim200$ Mpc and even larger fluctuations in smaller sub-volumes. In addition, we quantify fluid flow in and out of these sub-volumes. We find this to be significant, up to a $5\%$ change in the density between redshift $z=1$ and $z=0$ of a single sphere of radius $\sim200$ Mpc (and larger for smaller spheres). We suggest this may be important for studies basing averages on volumes co-moving with the simulation hypersurfaces.