The observational success and simplicity of the $\Lambda$CDM model, and the explicit analytic perturbations thereof, set the standard for any alternative cosmology. It therefore serves as a comparison ground and as a test case for methods which can be extended and applied to other cosmological models. In this paper we introduce dynamical systems and methods to describe linear scalar and tensor perturbations of the $\Lambda$CDM model, which serve as pedagogical examples that show the global illustrative powers of dynamical systems in the context of cosmological perturbations. We also study the asymptotic properties of the shear and Weyl tensors and discuss the validity of the perturbations as approximations to the Einstein field equations. Furthermore, we give a new approximation for the linear growth rate, $f(z) = \frac{d\ln \delta}{d\ln a} = \Omega^{\frac{6}{11}}_m - \frac{1}{70}(1-\Omega_m)^{\frac{5}{2}}$, where $z$ is the cosmological redshift, $\Omega_m=\Omega_m(z)$, while $a$ is the background scale factor, and show that it is much more accurate than the previous ones in the literature.

In this paper we consider second order perturbations of a flat Friedmann-Lema\^{i}tre universe whose stress-energy content is a single minimally coupled scalar field with an arbitrary potential. We derive the general solution of the perturbed Einstein equations in explicit form for this class of models when the perturbations are in the super-horizon regime. As a by-product we obtain a new conserved quantity for long wavelength perturbations of a single scalar field at second order.

The study of long wavelength scalar perturbations, in particular the existence of conserved quantities when the perturbations are adiabatic, plays an important role in e.g. inflationary cosmology. In this paper we present some new conserved quantities at second order and relate them to the curvature perturbation in the uniform density gauge, $\zeta$, and the comoving curvature perturbation, ${\cal R}$. We also, for the first time, derive the general solution of the perturbed Einstein equations at second order, which thereby contains both growing and decaying modes, for adiabatic long wavelength perturbations for a stress-energy tensor with zero anisotropic stresses and zero heat flux. The derivation uses the total matter gauge, but results are subsequently translated to the uniform curvature and Poisson (longitudinal, zero shear) gauges.

In this paper we give five gauge-invariant systems of governing equations for first and second order scalar perturbations of flat Friedmann-Lema\^{i}tre universes that are minimal in the sense that they contain no redundant equations or variables. We normalize the variables so that they are dimensionless, which leads to systems of equations that are simple and ready-to-use. We compare the properties and utility of the different systems. For example, they serve as a starting point for finding explicit solutions for two benchmark problems in cosmological perturbation theory at second order: adiabatic perturbations in the super-horizon regime (the long wavelength limit) and perturbations of $\Lambda$CDM universes. However, our framework has much wider applicability and serves as a reference for future work in the field.

We present simple expressions for the relativistic first and second order fractional density perturbations for FL cosmologies with dust, in four different gauges: the Poisson, uniform curvature, total matter and synchronous gauges. We include a cosmological constant and arbitrary spatial curvature in the background. A distinctive feature of our approach is our description of the spatial dependence of the perturbations using a canonical set of quadratic differential expressions involving an arbitrary spatial function that arises as a conserved quantity. This enables us to unify, simplify and extend previous seemingly disparate results. We use the primordial matter and metric perturbations that emerge at the end of the inflationary epoch to determine the additional arbitrary spatial function that arises when integrating the second order perturbation equations. This introduces a non-Gaussianity parameter into the expressions for the second order density perturbation. In the special case of zero spatial curvature we show that the time evolution simplifies significantly, and requires the use of only two non-elementary functions, the so-called growth supression factor at the linear level, and one new function at the second order level. We expect that the results will be useful in applications, for example, studying the effects of primordial non-Gaussianity on the large scale structure of the universe.

In this paper we present four simple expressions for the relativistic first and second order fractional density perturbations for $\Lambda$CDM cosmologies in different gauges: the Poisson, uniform curvature, total matter and synchronous gauges. A distinctive feature of our approach is the use of a canonical set of quadratic differential expressions involving an arbitrary spatial function, the so-called comoving curvature perturbation, to describe the spatial dependence, which enables us to unify, simplify and extend previous seemingly disparate results. The simple structure of the expressions makes the evolution of the density perturbations completely transparent and clearly displays the effect of the cosmological constant on the dynamics, namely that it stabilizes the perturbations. We expect that the results will be useful in applications, for example, studying the effects of primordial non-Gaussianity on the large scale structure of the universe.

Increasingly accurate observations of the cosmic microwave background and the large scale distribution of galaxies necessitate the study of nonlinear perturbations of Friedmann-Lemaitre cosmologies, whose equations are notoriously complicated. In this paper we present a new derivation of the governing equations for second order perturbations within the framework of the metric-based approach that is minimal, as regards amount of calculation and length of expressions, and flexible, as regards choice of gauge and stress-energy tensor. Because of their generality and the simplicity of their structure our equations provide a convenient starting point for determining the behaviour of nonlinear perturbations of FL cosmologies with any given stress-energy content, using either the Poisson gauge or the uniform curvature gauge.

Scalar perturbations of Friedmann-Lemaitre cosmologies can be analyzed in a variety of ways using Einstein's field equations, the Ricci and Bianchi identities, or the conservation equations for the stress-energy tensor, and possibly introducing a timelike reference congruence. The common ground is the use of gauge invariants derived from the metric tensor, the stress-energy tensor, or from vectors associated with a reference congruence, as basic variables. Although there is a complication in that there is no unique choice of gauge invariants, we will show that this can be used to advantage. With this in mind our first goal is to present an efficient way of constructing dimensionless gauge invariants associated with the tensors that are involved, and of determining their inter-relationships. Our second goal is to give a unified treatment of the various ways of writing the governing equations in dimensionless form using gauge-invariant variables, showing how simplicity can be achieved by a suitable choice of variables and normalization factors. Our third goal is to elucidate the connection between the metric-based approach and the so-called 1+3 gauge-invariant approach to cosmological perturbations. We restrict our considerations to linear perturbations, but our intent is to set the stage for the extension to second order perturbations.

In this paper we investigate asymptotic isotropization. We derive the asymptotic dynamics of spatially inhomogeneous cosmological models with a perfect fluid matter source and a positive cosmological constant near the de Sitter equilibrium state at late times, and near the flat FL equilibrium state at early times. Our results show that there exists an open set of solutions approaching the de Sitter state at late times, consistent with the cosmic no-hair conjecture. On the other hand, solutions that approach the flat FL state at early times are special and admit a so-called isotropic initial singularity. For both classes of models the asymptotic expansion of the line element contains an arbitrary spatial metric at leading order, indicating asymptotic spatial inhomogeneity. We show, however, that in the asymptotic regimes this spatial inhomogeneity is significant only at super-horizon scales.

We present a general framework for analyzing spatially inhomogeneous cosmological dynamics. It employs Hubble-normalized scale-invariant variables which are defined within the orthonormal frame formalism, and leads to the formulation of Einstein's field equations with a perfect fluid matter source as an autonomous system of evolution equations and constraints. This framework incorporates spatially homogeneous dynamics in a natural way as a special case, thereby placing earlier work on spatially homogeneous cosmology in a broader context, and allows us to draw on experience gained in that field using dynamical systems methods. One of our goals is to provide a precise formulation of the approach to the spacelike initial singularity in cosmological models, described heuristically by Belinski\v{\i}, Khalatnikov and Lifshitz. Specifically, we construct an invariant set which we conjecture forms the local past attractor for the evolution equations. We anticipate that this new formulation will provide the basis for proving rigorous theorems concerning the asymptotic behavior of spatially inhomogeneous cosmological models.

In this letter we will show that, contrary to what is widely believed, an almost isotropic cosmic microwave background (CMB) temperature does not imply that the universe is ``close to a Friedmann-Lemaitre universe''. There are two important manifestations of anisotropy in the geometry of the universe, (i) the anisotropy in the overall expansion, and (ii) the intrinsic anisotropy of the gravitational field, described by the Weyl curvature tensor, although the former usually receives more attention than the latter in the astrophysical literature. Here we consider a class of spatially homogeneous models for which the anisotropy of the CMB temperature is within the current observational limits but whose Weyl curvature is not negligible, i.e. these models are not close to isotropy even though the CMB temperature is almost isotropic.