We examine the spatial extrema (local maxima, minima and saddle points) of the covariant scalars (density, Hubble expansion, spatial curvature and eigenvalues of the shear and electric Weyl tensors) of the quasi-spherical Szekeres dust models. Sufficient conditions are obtained for the existence of distributions of multiple extrema in spatial comoving locations that can be prescribed through initial conditions. These distributions evolve without shell crossing singularities at least for ever expanding models (with or without cosmological constant) in the full evolution range where the models are valid. By considering the local maxima and minima of the density, our results allow for setting up elaborated networks of "pancake" shaped evolving cold dark matter over-densities and density voids whose spatial distribution and amplitudes can be controlled from initial data compatible with standard early Universe initial conditions. We believe that these results have an enormous range of potential application by providing a fully relativistic non-perturbative coarse grained modelling of cosmic structure at all scales.

We show that the full dynamical freedom of the well known Szekeres models allows for the description of elaborated 3-dimensional networks of cold dark matter structures (over-densities and/or density voids) undergoing "pancake" collapse. By reducing Einstein's field equations to a set of evolution equations, which themselves reduce in the linear limit to evolution equations for linear perturbations, we determine the dynamics of such structures, with the spatial comoving location of each structure uniquely specified by standard early Universe initial conditions. By means of a representative example we examine in detail the density contrast, the Hubble flow and peculiar velocities of structures that evolved, from linear initial data at the last scattering surface, to fully non-linear 10-20 Mpc scale configurations today. To motivate further research, we provide a qualitative discussion on the connection of Szekeres models with linear perturbations and the pancake collapse of the Zeldovich approximation. This type of structure modelling provides a coarse grained -but fully relativistic non-linear and non-perturbative -description of evolving large scale cosmic structures before their virialisation, and as such it has an enormous potential for applications in cosmological research.

We study the differences and equivalences between the non-perturbative description of the evolution of cosmic structure furnished by the Szekeres dust models (a non-spherical exact solution of Einstein's equations) and the dynamics of Cosmological Perturbation Theory (CPT) for dust sources in a ΛCDM background. We show how the dynamics of Szekeres models can be described by evolution equations given in terms of "exact fluctuations" that identically reduce (at all orders) to evolution equations of CPT in the comoving isochronous gauge. We explicitly show how Szekeres linearised exact fluctuations are specific (deterministic) realisations of standard linear perturbations of CPT given as random fields but, as opposed to the latter perturbations, they can be evolved exactly into the full non-linear regime. We prove two important results: (i) the conservation of the curvature perturbation (at all scales) also holds for the appropriate linear approximation of the exact Szekeres fluctuations in a ΛCDM background, and (ii) the different collapse morphologies of Szekeres models yields, at nonlinear order, different functional forms for the growth factor that follows from the study of redshift space distortions. The metric based potentials used in linear CPT are computed in terms of the parameters of the linearised Szekeres models, thus allowing us to relate our results to linear CPT results in other gauges. We believe that these results provide a solid starting stage to examine the role of non-perturbative General Relativity in current cosmological research.PACS numbers: 98.80.-k, 04.20.-q, 95.36.+x, 95.35.+d

We examine the fully relativistic evolution of cosmic voids constituted by baryons and cold dark matter (CDM), represented by two non-comoving dust sources in a ΛCDM background. For this purpose, we consider numerical solutions of Einstein's field equations in a fluid-flow representation adapted to spherical symmetry and multiple components. We present a simple example that explores the frame-dependence of the local expansion and the Hubble flow for this mixture of two dusts, revealing that the relative velocity between the sources yields a significantly different evolution in comparison with that of the two sources in a common 4-velocity (which reduces to a Lemaître-Tolman-Bondi model). In particular, significant modifications arise for the density contrast depth and void size, as well as in the amplitude of the surrounding over-densities. We show that an adequate model of a frame-dependent evolution that incorporates initial conditions from peculiar velocities and large-scale density contrast observations may contribute to understand the discrepancy between the local value of H 0 and that inferred from the CMB.

We examine the relation between the Szekeres models and relativistic Lagrangian perturbation schemes, in particular the relativistic Zel'dovich approximation (RZA). We show that the second class of the Szekeres solutions is exactly contained within the RZA when the latter is restricted to an irrotational dust source with a flow-orthogonal foliation of spacetime. In such a case, the solution is governed by the first principal scalar invariant of the deformation field, proving a direct connection with a class of Newtonian three-dimensional solutions without symmetry. For the second class, a necessary and sufficient condition for the vanishing of cosmological backreaction on a scale of homogeneity is expressed through integral constraints. Domains with no backreaction can be smoothly matched, forming a lattice model, where exact deviations average out at a given scale of homogeneity, and the homogeneous and isotropic background is recovered as an average property of the model. Although the connection with the first class of Szekeres solutions is not straightforward, this class allows for the interpretation in terms of a spatial superposition of nonintersecting fluid lines, where each world line evolves independently and under the RZA model equations, but with different associated "local backgrounds". This points to the possibility of generalizing the Lagrangian perturbation schemes to structure formation models on evolving backgrounds, including global cosmological backreaction.