The existence of static self-gravitating Newtonian elastic balls is proved under general assumptions on the constitutive equations of the elastic material. The proof uses methods from the theory of finite-dimensional dynamical systems and the Euler formulation of elasticity theory for spherically symmetric bodies introduced recently by the authors. Examples of elastic materials covered by the results of this paper are Saint Venant-Kirchhoff, John and Hadamard materials.

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.

We consider the wave equation, $\square_g\psi=0$, in fixed flat Friedmann-Lema\^itre-Robertson-Walker and Kasner spacetimes with topology $\mathbb{R}_+\times\mathbb{T}^3$. We obtain generic blow up results for solutions to the wave equation towards the Big Bang singularity in both backgrounds. In particular, we characterize open sets of initial data prescribed at a spacelike hypersurface close to the singularity, which give rise to solutions that blow up in an open set of the Big Bang hypersurface $\{t=0\}$. The initial data sets are characterized by the condition that the Neumann data should dominate, in an appropriate $L^2$-sense, up to two spatial derivatives of the Dirichlet data. For these initial configurations, the $L^2(\mathbb{T}^3)$ norms of the solutions blow up towards the Big Bang hypersurfaces of FLRW and Kasner with inverse polynomial and logarithmic rates respectively. Our method is based on deriving suitably weighted energy estimates in physical space. No symmetries of solutions are assumed.

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 study the dynamics of spatially homogeneous and isotropic spacetimes containing a fluid undergoing microscopic velocity diffusion in a cosmological scalar field. After deriving a few exact solutions of the equations, we continue by analyzing the qualitative behavior of general solutions. To this purpose we recast the equations in the form of a two dimensional dynamical system and perform a global analysis of the flow. Among the admissible behaviors, we find solutions that are asymptotically de-Sitter both in the past and future time direction and which undergo accelerated expansion at all times.

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.