Replacing dark energy by silent virialisation

Abstract

Context. Standard cosmological N-body simulations have background scale factor evolution that is decoupled from non-linear structure formation. Prior to gravitational collapse, kinematical backreaction (ℚ'¬ ℚ') justifies this approach in a Newtonian context. Aims. However, the final stages of a gravitational collapse event are sudden; a globally imposed smooth expansion rate forces at least one expanding region to suddenly and instantaneously decelerate in compensation for the virialisation event. This is relativistically unrealistic. A more conservative hypothesis is to allow non-collapsed domains to continue their volume evolution according to the ℚ'¬ ℚ' Zel'dovich approximation (QZA). We aim to study the inferred average expansion under this "silent" virialisation hypothesis. Methods. We set standard (MPGRAFIC) EdS 3-torus (T 3 ) cosmological N-body initial conditions. Using RAMSES, we partitioned the volume into domains and called the DTFE library to estimate the per-domain initial values of the three invariants of the extrinsic curvature tensor that determine the QZA. We integrated the Raychaudhuri equation in each domain using the INHOMOG library, and adopted the stable clustering hypothesis to represent virialisation (VQZA). We spatially averaged to obtain the effective global scale factor. We adopted an early-epoch-normalised EdS reference-model Hubble constant H 1 EdS = 37.7 km s -1 â• Mpc and an effective Hubble constant H eff,0 = 67.7 km s -1 â• Mpc. Results. From 2000 simulations at resolution 256 3 , we find that reaching a unity effective scale factor at 13.8 Gyr (16% above EdS), occurs for an averaging scale of L 13.8 13 = 2.5 -0.1 +0.1 Mpc/h eff . Relativistically interpreted, this corresponds to strong average negative curvature evolution, with the mean (median) curvature functional Ω â., ℚ growing from zero to about 1.5-2 by the present. Over 100 realisations, the virialisation fraction and super-EdS expansion correlate strongly at fixed cosmological time. Conculsions. Thus, starting from EdS initial conditions and averaging on a typical non-linear structure formation scale, the VQZA dark-energy-free average expansion matches ΛCDM expansion to first order. The software packages used here are free-licensed.