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 study flat FLRW $\alpha$-attractor $\mathrm{E}$- and $\mathrm{T}$-models by introducing a dynamical systems framework that yields regularized unconstrained field equations on two-dimensional compact state spaces. This results in both illustrative figures and a complete description of the entire solution spaces of these models, including asymptotics. In particular, it is shown that observational viability, which requires a sufficient number of $e$-folds, is associated with a solution given by a one-dimensional center manifold of a past asymptotic de Sitter state, where the center manifold structure also explains why nearby solutions are attracted to this `inflationary attractor solution.' A center manifold expansion yields a description of the inflationary regime with arbitrary analytic accuracy, where the slow-roll approximation asymptotically describes the tangency condition of the center manifold at the asymptotic de Sitter state.

We discuss dynamical systems approaches and methods applied to flat Robertson-Walker models in $f(R)$-gravity. We argue that a complete description of the solution space of a model requires a global state space analysis that motivates globally covering state space adapted variables. This is shown explicitly by an illustrative example, $f(R) = R + \alpha R^2$, $\alpha > 0$, for which we introduce new regular dynamical systems on global compactly extended state spaces for the Jordan and Einstein frames. This example also allows us to illustrate several local and global dynamical systems techniques involving, e.g., blow ups of nilpotent fixed points, center manifold analysis, averaging, and use of monotone functions. As a result of applying dynamical systems methods to globally state space adapted dynamical systems formulations, we obtain pictures of the entire solution spaces in both the Jordan and the Einstein frames. This shows, e.g., that due to the domain of the conformal transformation between the Jordan and Einstein frames, not all the solutions in the Jordan frame are completely contained in the Einstein frame. We also make comparisons with previous dynamical systems approaches to $f(R)$ cosmology and discuss their advantages and disadvantages.

This paper treats nonrelativistic matter and a scalar field $\phi$ with a monotonically decreasing potential minimally coupled to gravity in flat Friedmann-Lema\^{i}tre-Robertson-Walker cosmology. The field equations are reformulated as a three-dimensional dynamical system on an extended compact state space, complemented with cosmographic diagrams. A dynamical systems analysis provides global dynamical results describing possible asymptotic behavior. It is shown that one should impose \emph{global and asymptotic} bounds on $\lambda=-V^{-1}\,dV/d\phi$ to obtain viable cosmological models that continuously deform $\Lambda$CDM cosmology. In particular we introduce a regularized inverse power-law potential as a simple specific example.

We consider a minimally coupled scalar field with a monomial potential and a perfect fluid in flat FLRW cosmology. We apply local and global dynamical systems techniques to a new three-dimensional dynamical systems reformulation of the field equations on a compact state space. This leads to a visual global description of the solution space and asymptotic behavior. At late times we employ averaging techniques to prove statements about how the relationship between the equation of state of the fluid and the monomial exponent of the scalar field affects asymptotic source dominance and asymptotic manifest self-similarity breaking. We also situate the `attractor' solution in the three-dimensional state space and show that it corresponds to the one-dimensional unstable center manifold of a de Sitter fixed point, located on an unphysical boundary associated with the dynamics at early times. By deriving a center manifold expansion we obtain approximate expressions for the attractor solution. We subsequently improve the accuracy and range of the approximation by means of Pad\'e approximants and compare with the slow-roll approximation.

We consider the familiar problem of a minimally coupled scalar field with quadratic potential in flat Friedmann-Lema\^itre-Robertson-Walker cosmology to illustrate a number of techniques and tools, which can be applied to a wide range of scalar field potentials and problems in e.g. modified gravity. We present a global and regular dynamical systems description that yields a global understanding of the solution space, including asymptotic features. We introduce dynamical systems techniques such as center manifold expansions and use Pad\'e approximants to obtain improved approximations for the `attractor solution' at early times. We also show that future asymptotic behavior is associated with a limit cycle, which shows that manifest self-similarity is asymptotically broken toward the future, and give approximate expressions for this behavior. We then combine these results to obtain global approximations for the attractor solution, which, e.g., might be used in the context of global measures. In addition we elucidate the connection between slow-roll based approximations and the attractor solution, and compare these approximations with the center manifold based approximants.

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.

In this paper we investigate the global dynamics for the minimally coupled scalar field representation of the modified Chaplygin gas in the context of flat FLRW cosmology. The tool for doing this is a new set of bounded variables that lead to a regular dynamical system. It is shown that the exact modified Chaplygin gas perfect fluid solution appears as a straight line in the associated phase plane. It is also shown that no other solutions stay close to this solution during their entire temporal evolution, but that there exists an open subset of solutions that stay arbitrarily close during an intermediate time interval, and into the future in the case the scalar field potential exhibits a global minimum.

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.

We investigate spherically symmetric equilibrium states of the Vlasov-Poisson system, relevant in galactic dynamics. We recast the equations into a regular three-dimensional system of autonomous first order ordinary differential equations on a region with compact closure.Based on a dynamical systems analysis we derive theorems that guarantee that the steady state solutions have finite mass and compact support.

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 investigate relativistic spherically symmetric static perfect fluid models in the framework of the theory of dynamical systems. The field equations are recast into a regular dynamical system on a 3-dimensional compact state space, thereby avoiding the non-regularity problems associated with the Tolman-Oppenheimer-Volkoff equation. The global picture of the solution space thus obtained is used to derive qualitative features and to prove theorems about mass-radius properties. The perfect fluids we discuss are described by barotropic equations of state that are asymptotically polytropic at low pressures and, for certain applications, asymptotically linear at high pressures. We employ dimensionless variables that are asymptotically homology invariant in the low pressure regime, and thus we generalize standard work on Newtonian polytropes to a relativistic setting and to a much larger class of equations of state. Our dynamical systems framework is particularly suited for numerical computations, as illustrated by several numerical examples, e.g., the ideal neutron gas and examples that involve phase transitions.

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.

We investigate static spherically symmetric perfect fluid models in Newtonian gravity for barotropic equations of state that are asymptotically polytropic at low and high pressures. This is done by casting the equations into a 3-dimensional regular dynamical system with bounded dependent variables. The low and high central pressure limits correspond to two 2-dimensional boundary subsets, described by homology invariant equations for exact polytropes. Thus the formulation naturally places work about polytropes in a more general context. The introduced framework yields a visual aid for obtaining qualitative information about the solution space and is also suitable for numerical investigations. Last, but not least, it makes a host of mathematical tools from dynamical systems theory available. This allows us to prove a number of theorems about the relationship between the equation of state and properties concerning total masses and radii.

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.

We revisit the issue of integrability conditions for the irrotational silent cosmological models. We formulate the problem both in 1+3 covariant and 1+3 orthonormal frame notation, and show there exists a series of constraint equations that need to be satisfied. These conditions hold identically for FLRW-linearised silent models, but not in the general exact non-linear case. Thus there is a linearisation instability, and it is highly unlikely that there is a large class of silent models. We conjecture that there are no spatially inhomogeneous solutions with Weyl curvature of Petrov type I, and indicate further issues that await clarification.

After a brief overview of the so-called silent models and their present status, we consider the subclass of Bianchi Type--I models with a magnetic field source. Due to the presence of the magnetic field, the initial singularity shows ``oscillatory'' features reminiscent of the Bianchi Type--IX case. The Bianchi Type--I models with a magnetic field are therefore a counterexample to the folklore that matter fields can be neglected in the vicinity of the singularity.