Relativistic cosmological perturbation scheme on a general background: scalar perturbations for irrotational dust

In standard perturbation approaches and N-body simulations, inhomogeneities are described to evolve on a predefined background cosmology, commonly taken as the homogeneous-isotropic solutions of Einstein's field equations (Friedmann-Lemaitre-Robertson-Walker (FLRW) cosmologies). In order to make physical sense, this background cosmology must provide a reasonable description of the effective, i.e. spatially averaged, evolution of structure inhomogeneities also in the nonlinear regime. Guided by the insights that (i) the average over an inhomogeneous distribution of matter and geometry is in general not given by a homogeneous solution of general relativity, and that (ii) the class of FLRW cosmologies is not only locally but also globally gravitationally unstable in relevant cases, we here develop a perturbation approach that describes the evolution of inhomogeneities on a general background being defined by the spatially averaged evolution equations. This physical background interacts with the formation of structures. We derive and discuss the resulting perturbation scheme for the matter model 'irrotational dust' in the Lagrangian picture, restricting our attention to scalar perturbations.