In this book are studied, from the perspective of the dynamical systems, several Universe models. In chapter 1 we give a bird's eye view on cosmology and cosmological problems. Chapter 2 is devoted to a brief review on some results and useful tools from the qualitative theory of dynamical systems. They provide the theoretical basis for the qualitative study of concrete cosmological models. Chapters 1 and 2 are a review of well-known results. Chapters 3, 4, 5 and 6 are devoted to our main results. In these chapters are extended and settled in a substantially different, more strict mathematical language, several results obtained by one of us in arXiv:0812.1013 [gr-qc]; arXiv:1009.0689 [gr-qc]; arXiv:0904.1577[gr-qc]; and arXiv:0909.3571 [hep-th]. In chapter 6, we provide a different approach to the subject discussed in astro-ph/0503478. Additionally, we perform a Poincar\'e compactification process allowing to construct a global phase space containing all the cosmological information in both finite and infinite regions for all the models.

We perform a detailed dynamical analysis of anisotropic scalar-field cosmologies, and in particular of the most significant Kantowski-Sachs, Locally Rotationally Symmetric (LRS) Bianchi I and LRS Bianchi III cases. We follow the new and powerful method of $f$-devisers, which allows us to perform the whole analysis for a wide range of potentials. Thus, one can just substitute the specific potential form in the final results and obtain the corresponding behavior, without the need of new calculations. We find a very rich behavior, and amongst others the universe can result in isotropized solutions with observables in agreement with observations, such as de Sitter, quintessence-like, or stiff-dark energy solutions. In particular, all expanding, accelerating, stable attractors are isotropic. Additionally, we prove that as long as matter obeys the null energy condition, bounce behavior is impossible. Finally, applying the general results to the well-studied exponential and power-law potentials, we find that some of the general stable solutions disappear. This feature may be an indication that such simple potentials might restrict the dynamics in scalar-field cosmology, opening the way to the introduction of more complicated ones.

In this work we present a phase space analysis of a quintessence field and a perfect fluid trapped in a Randall-Sundrum's Braneworld of type 2. We consider a homogeneous but anisotropic Bianchi I brane geometry. Moreover, we consider the effect of the projection of the five dimensional Weyl tensor onto the three-brane in the form of a negative Dark Radiation term. For the treatment of the potential we use the "Method of $f$-devisers" that allows investigating arbitrary potentials in a phase space. We present general conditions on the potential in order to obtain the stability of standard 4D and non-standard 5D de Sitter solutions, and we provide the stability conditions for both scalar field-matter scaling solutions, scalar field-dark radiation solutions and scalar field-dominated solutions. We find that the shear-dominated solutions are unstable (particularly, contracting shear-dominated solutions are of saddle type). As a main difference with our previous work, the traditionally ever-expanding models could potentially re-collapse due to the negativity of the dark radiation. Additionally, our system admits a large class of static solutions that are of saddle type. These kinds of solutions are important at intermediate stages in the evolution of the universe, since they allow the transition from contracting to expanding models and viceversa. New features of our scenario are the existence of a bounce and a turnaround, which lead to cyclic behavior, that are not allowed in Bianchi I branes with positive dark radiation term. Finally, as specific examples we consider the potentials $V\propto\sinh^{-\alpha}(\beta\phi)$ and $V\propto\left[\cosh\left(\xi\phi \right)-1\right]$ which have simple $f$-devisers.

In this paper we consider a cosmological model whose main components are a scalar field and a generalized Chaplygin gas. We obtain an exact solution for a flat arbitrary potential. This solution have the right dust limit when the Chaplygin parameter $A\rightarrow 0$. We use the dynamical systems approach in order to describe the cosmological evolution of the mixture for an exponential self-interacting scalar field potential. We study the scalar field with an arbitrary self-interacting potential using the "Method of $f$-devisers." Our results are illustrated for the special case of a coshlike potential. We find that usual scalar-field-dominated and scaling solutions cannot be late-time attractors in the presence of the Chaplygin gas (with $\alpha>0$). We recover the standard results at the dust limit ($A\rightarrow 0$). In particular, for the exponential potential, the late-time attractor is a pure generalized Chaplygin solution mimicking an effective cosmological constant. In the case of arbitrary potentials, the late-time attractors are de Sitter solutions in the form of a cosmological constant, a pure generalized Chaplygin solution or a continuum of solutions, when the scalar field and the Chaplygin gas densities are of the same orders of magnitude. The different situations depend on the parameter choices.

We investigate the phase-space of a flat FRW universe including both a scalar field, $\phi,$ coupled to matter, and radiation. The model is inspired in scalar-tensor theories of gravity, and thus, related with $F(R)$ theories through conformal transformation. The aim of the chapter is to extent several results to the more realistic situation when radiation is included in the cosmic budget particularly for studying the early time dynamics. Under mild conditions on the potential we prove that the equilibrium points corresponding to the non-negative local minima for $V(\phi)$ are asymptotically stable. Normal forms are employed to obtain approximated solutions associated to the inflection points and the strict degenerate local minimum of the potential. We prove for arbitrary potentials and arbitrary coupling functions $\chi(\phi),$ of appropriate differentiable class, that the scalar field almost always diverges into the past. It is designed a dynamical system adequate to studying the stability of the critical points in the limit $|\phi|\to\infty.$ We obtain there: radiation-dominated cosmological solutions; power-law scalar-field dominated inflationary cosmological solutions; matter-kinetic-potential scaling solutions and radiation-kinetic-potential scaling solutions. Using the mathematical apparatus developed here, we investigate the important examples of higher order gravity theories $F(R) = R + \alpha R^2$ (quadratic gravity) and $F(R) =R^n.$ We illustrated both analytically and numerically our principal results. In the case of quadratic gravity we prove, by an explicit computation of the center manifold, that the equilibrium point corresponding to de Sitter solution is locally asymptotically unstable (saddle point).