The growth of structure in the Szekeres inhomogeneous cosmological models and the matter-dominated era

Abstract: This study belongs to a series devoted to using Szekeres inhomogeneous models to develop a theoretical framework where observations can be investigated with a wider range of possible interpretations. We look here into the growth of large-scale structure in the models. The Szekeres models are exact solutions to Einstein's equations that were originally derived with no symmetries. We use a formulation of the models that is due to Goode and Wainwright, who considered the models as exact perturbations of an FLRW background. Using the Raychaudhuri equation, we write for the two classes of the models, exact growth equations in terms of the under/overdensity and measurable cosmological parameters. The new equations in the overdensity split into two informative parts. The first part, while exact, is identical to the growth equation in the usual linearly perturbed FLRW models, while the second part constitutes exact non-linear perturbations. We integrate numerically the full exact growth rate equations for the flat and curved cases. We find that for the matter-dominated era, the Szekeres growth rate is up to a factor of three to five stronger than the usual linearly perturbed FLRW cases, reflecting the effect of exact Szekeres non-linear perturbations. The growth is also stronger than that of the non-linear spherical collapse model, and the difference between the two increases with time. This highlights the distinction when we use general inhomogeneous models where shear and a tidal gravitational field are present and contribute to the gravitational clustering. Additionally, it is worth observing that the enhancement of the growth found in the Szekeres models during the matter-dominated era could suggest a substitute to the argument that dark matter is needed when using FLRW models to explain the enhanced growth and resulting large-scale structures that we observe today (abridged)