Recently \citet{PEDE:2019ApJ} proposed a phenomenologically emergent dark energy (PEDE) which consider that the dark energy density evolves as $\widetilde{\Omega}_{\rm{DE}}(z)\,=\,\Omega_{\rm{DE,0}}\left[ 1 - {\rm{tanh}}\left( {\log}_{10}(1+z) \right) \right]$ with the advantage that it does not have degree of freedom. Later on, \citet{PEDE:2020} proposed a generalized model by adding one degree of freedom to the PEDE model, encoded in the parameter $\Delta$. Motivated by these proposals, we constrain the parameter space ($h,\Omega_m$) and ($h,\Omega_m, \Delta$) for PEDE and Generalized Emergent Dark Energy (GEDE) respectively, by employing the most recent observational (non-) homogeneous Hubble data. Additionally, we reconstruct the deceleration and jerk parameters and estimate yield values at $z=0$ of $q_0 = -0.784^{+0.028}_{-0.027}$ and $j_0 = 1.241^{+0.164}_{-0.149}$ for PEDE and $q_0 = -0.730^{+0.059}_{-0.067}$ and $j_0 = 1.293^{+0.194}_{-0.187}$ for GEDE using the homogeneous sample. We report values on the deceleration-acceleration transition redshift with those reported in the literature within $2\sigma$ CL. Furthermore, we perform a stability analysis of the PEDE and GEDE models to study the global evolution of the Universe around their critical points. Although the PEDE and GEDE dynamics are similar to the standard model, our stability analysis indicates that in both models there is an accelerated phase at early epochs of the Universe.

In a recent paper [J. Martin and V. Vennin, Phys. Rev. Lett. 124, 080402 (2020)] it was argued that, for most natural choices, the direct application of the continuous spontaneous localization (CSL) theory to the inflationary case, as it is known to work in non-relativistic laboratory situations, is ruled out by cosmological observational data, thus casting a shadow on models based on CSL theory. We point out that such results are based on the consideration of a rather narrow set of choices for the application of the theory to the cosmological context and that the landscape of open and different possibilities is extremely vast.

Cosmic inflation, which describes an accelerated expansion of the early Universe, yields the most successful predictions regarding temperature anisotropies in the cosmic microwave background (CMB). Nevertheless, the precise origin of the primordial perturbations and their quantum-to-classical transition is still an open issue. The continuous spontaneous localization model (CSL), in the cosmological context, might be used to provide a solution to the mentioned puzzles by considering an objective reduction of the inflaton wave function. In this work, we calculate the primordial power spectrum at the next leading order in the Hubble flow functions that results from applying the CSL model to slow roll inflation within the semiclassical gravity framework. We employ the method known as uniform approximation along with a second order expansion in the Hubble flow functions. We analyze some features in the CMB temperature and primordial power spectra that could help to distinguish between the standard prediction and our approach.

In the context of Einstein-aether scalar field cosmology we solve the field equations and determine exact and analytic solutions. In particular, we consider a model proposed by Kanno and Soda where the aether and the scalar fields interact through the aether coefficient parameters, which are promoted to be functions of the scalar field. For this model, we write the field equations by using the minisuperspace approach and we determine the scalar field potentials which leads to Liouville--integrable systems. We solve the field equations for five families of scalar field potentials and, whether it is feasible, we write down the analytic solutions by using closed-form functions.

In the present article we study the cosmological evolution of a two-scalar field gravitational theory defined in the Jordan frame. Specifically, we assume one of the scalar fields to be minimally coupled to gravity, while the second field which is the Brans-Dicke scalar field is nonminimally coupled to gravity and also coupled to the other scalar field. In the Einstein frame this theory reduces to a two-scalar field theory where the two fields can interact only in the potential term, which means that the quintom theory is recovered. The cosmological evolution is studied by analyzing the equilibrium points of the field equations in the Jordan frame. We find that the theory can describe the cosmological evolution in large scales, while inflationary solutions are also provided.

We examine the phase space of Ho\v{r}ava-Lifshitz cosmology for a wide range of self-interacting potentials for the scalar field under the detailed-balance condition and without imposing it, by means of the powerful method of $f$-devisers. A compactification approach is performed for the exponential potential and for potentials beyond the exponential one, extending the previous findings in the literature. By using this approach it is possible to describe the finite region of the phase space and the region where the phase-space variables becomes infinity. Furthermore, we present several results concerning the stability of the \emph{de Sitter} solution in Ho\v{r}ava-Lifshitz cosmology using Center Manifold theory. The advantages of this procedure are unveiled immediately when it is compared with the Normal Forms Calculations presented before in the literature.

We perform for the first time a dynamical system analysis of both the background and perturbation equations, of $\Lambda$CDM cosmology and quintessence scenario with an exponential potential. In the former case the perturbations do not change the stability of the late-time attractor of the background equations, and the system still results in the dark-energy dominated, de Sitter solution, having passed from the correct dark-matter era with $\gamma\approx6/11$. However, in the case of quintessence the incorporation of perturbations changes the stability and properties of the background evolution, and the only conditionally stable points present either an exponentially increasing matter clustering not favored by observations, or Laplacian instabilities, and thus not physically interesting. This result is a severe disadvantage of quintessence cosmology comparing to $\Lambda$CDM paradigm.

The physical mechanism responsible for the emergence of primordial cosmic seeds from a perfect isotropic and homogeneous Universe has not been fully addressed in standard cosmic inflation. To handle this shortcoming, D. Sudarsky et al have developed a proposal: the self-induced collapse hypothesis. In this scheme, the objective collapse of the inflaton's wave function generates the inhomogeneity and anisotropy at all scales. In this paper we analyze the viability of a set of inflationary potentials in both the context of the collapse proposal and within the standard inflationary framework. For this, we perform a statistical analysis using recent CMB and BAO data to obtain the prediction for the scalar spectral index $n_s$ in the context of a particular collapse model: the Wigner scheme. The predicted $n_s$ and the tensor-to-scalar ratio $r$ in terms of the slow roll parameters is different between the collapse scheme and the standard inflationary model. For each potential considered we compare the prediction of $n_s$ and $r$ with the limits established by observational data in both pictures. The result of our analysis shows in most cases a difference in the inflationary potentials allowed by the observational limits in both frameworks. In particular, in the standard approach the more concave a potential is, the more is favored by the data. On the other hand, in the Wigner scheme, the data favors equally all type of concave potentials, including those at the border between convex and concave families.

In a higher-order modified teleparallel theory cosmological we present analytical cosmological solutions. In particular we determine forms of the unknown potential which drives the scalar field such that the field equations form a Liouville integrable system. For the determination of the conservation laws we apply the Cartan symmetries. Furthermore, inspired from our solutions, a toy model is studied and it is shown that it can describe the Supernova data, while at the same time introduces dark matter components in the Hubble function. When the extra matter source is a stiff fluid then we show how analytical solutions for Bianchi I universes can be constructed from our analysis. Finally, we perform a global dynamical analysis of the field equations by using variables different from that of the Hubble-normalization.

Recently, the nonlinear electrodynamics (NED) has been gaining attention to generate primordial magnetic fields in the Universe and also to resolve singularity problems. Moreover, recent works have shown the crucial role of the NED on the inflation. This paper provides a new approach based on a new model of NED as a source of gravitation to remove the cosmic singularity at the big bang and explain the cosmic acceleration during the inflation era on the background of stochastic magnetic field. Also, we found a realization of a cyclic Universe, free of initial singularity, due to the proposed NED energy density. In addition, we explore whether a NED field without or with matter can be the origin of the late-time acceleration. For this we obtain explicit equations for $H(z)$ and perform a MCMC analysis to constrain the NED parameters by using $31$ observational Hubble data (OHD) obtained from cosmic chronometers covering the redshift range $0 < z < 1.97$; and with the joint-light-analysis (JLA) SNIa compilation consisting in $740$ data points in the range $0.01<z<1.2$. All our constraints on the current magnetic field give $B_{0}\sim 10^{-31}\mathrm{cm^{-1}}$, which are larger than the upper limit $10^{-33}\mathrm{cm^{-1}}$ by the Planck satellite implying that NED cosmologies could not be suitable to explain the Universe late-time dynamics. However, the current data is able to falsify the scenario at late times. Indeed, one is able to reconstruct the deceleration parameter $q(z)$ using the best fit values of the parameters obtained from OHD and SNIa data sets. If the matter component is not included, the data sets predict an accelerated phase in the early Universe, but a non accelerated Universe is preferred in the current epoch...

We consider the eternal inflation scenario of the slow-roll/chaotic type with the additional element of an objective collapse of the wave function. The incorporation of this new agent to the traditional inflationary setting might represent a possible solution to the quantum measurement problem during inflation, a subject that has not reached a consensus among the community. Specifically, it could provide an explanation for the generation of the primordial anisotropies and inhomogeneities, starting from a perfectly symmetric background and invoking symmetric dynamics. We adopt the continuous spontaneous localization model, in the context of inflation, as the dynamical reduction mechanism that generates the primordial inhomogeneities. Furthermore, when enforcing the objective reduction mechanism, the condition for eternal inflation can be bypassed. In particular, the collapse mechanism incites the wave function, corresponding to the inflaton, to localize itself around the zero mode of the field. Then the zero mode will evolve essentially unperturbed, driving inflation to an end in any region of the Universe where inflation occurred. Also, our approach achieves a primordial spectrum with an amplitude and shape consistent with the one that best fits the observational data.

A class of generalized Galileon cosmological models, which can be described by a point-like Lagrangian, is considered in order to utilize Noether's Theorem to determine conservation laws for the field equations. In the Friedmann-Lema\^itre-Robertson-Walker universe, the existence of a nontrivial conservation law indicates the integrability of the field equations. Due to the complexity of the latter, we apply the differential invariants approach in order to construct special power-law solutions and study their stability.

We investigate Kantowski-Sachs models in Einstein-{\ae}ther theory with a perfect fluid source using the singularity analysis to prove the integrability of the field equations and dynamical system tools to study the evolution. We find an inflationary source at early times, and an inflationary sink at late times, for a wide region in the parameter space. The results by A. A. Coley, G. Leon, P. Sandin and J. Latta (JCAP 12, 010, 2015), are then re-obtained as particular cases. Additionally, we select other values for the non-GR parameters which are consistent with current constraints, getting a very rich phenomenology. In particular, we find solutions with infinite shear, zero curvature, and infinite matter energy density in comparison with the Hubble scalar. We also have stiff-like future attractors, anisotropic late-time attractors, or both, in some special cases. Such results are developed analytically, and then verified by numerics. Finally, the physical interpretation of the new critical points is discussed.

In recent years, some studies have drawn attention to the lack of large-angle correlations in the observed cosmic microwave background (CMB) temperature anisotropies with respect to that predicted within the standard $\Lambda$CDM model. Lately, it has been argued that such a lack of correlations could be explained in the framework of the so-called $R_h=ct$ model without inflation. The aim of this work is to study whether there is a mechanism to generate, through a quantum field theory, the primordial power spectrum presented by these authors. Specifically, we consider two different scenarios: first, we assume a scalar field dominating the early Universe in the $R_h=ct$ cosmological model, and second, we deal with the possibility of adding an early inflationary phase to the mentioned model. During the analysis of the consistency between the predicted and observed amplitudes of the CMB temperature anisotropies in both scenarios, we run into deep issues which indicate that it is not clear how to characterize the primordial quantum perturbations within the $R_h=ct$ model.

The Continuous Spontaneous Localization (CSL) model has been proposed as a possible solution to the quantum measurement problem by modifying the Schr\"{o}dinger equation. In this work, we apply the CSL model to two cosmological models of the early Universe: the matter bounce scenario and slow roll inflation. In particular, we focus on the generation of the classical primordial inhomogeneities and anisotropies that arise from the dynamical evolution, provided by the CSL mechanism, of the quantum state associated to the quantum fields. In each case, we obtained a prediction for the shape and the parameters characterizing the primordial spectra (scalar and tensor), i.e. the amplitude, the spectral index and the tensor-to-scalar ratio. We found that there exist CSL parameter values, allowed by other non-cosmological experiments, for which our predictions for the angular power spectrum of the CMB temperature anisotropy are consistent with the best fit canonical model to the latest data released by the Planck Collaboration.

We investigate a Jordan-Brans-Dicke (JBD) scalar field, $\Phi$, with power-law potential in the presence of a second scalar field, $\phi$, with an exponential potential, in both the Jordan and the Einstein frames. We present the relation of our model with the induced gravity model with power-law potential and the integrability of this kind of models is discussed when the quintessence field $\phi$ is massless, and has a small velocity. We prove that in JBD theory, the de Sitter solution is not a natural attractor but an intermediate accelerated solution of the form $a(t)\simeq e^{\alpha_1 t^{p_1}}$, as $t\rightarrow \infty$ where $\alpha_1>0$ and $0<p_1<1$, for a wide range of parameters. Furthermore, in the Einstein frame we get that the attractor is also an intermediate accelerated solution of the form $\mathfrak{a}(\mathfrak{t})\simeq e^{\alpha_2 \mathfrak{t}^{p_2}}$ as $\mathfrak{t}\rightarrow \infty$ where $\alpha_2>0$ and $0<p_2<1$, for the same conditions on the parameters as in the Jordan frame. In the special case of a quadratic potential in the Jordan frame, or for a constant potential in the Einstein's frame, these solutions are of saddle type. Finally, we present a specific elaboration of our extension of the induced gravity model in the Jordan frame, which corresponds to a linear potential of $\Phi$. The dynamical system is then reduced to a two dimensional one, and the late-time attractor is linked with the exact solution found for the induced gravity model. In this example the intermediate accelerated solution does not exist, and the attractor solution has an asymptotic de Sitter-like evolution law for the scale factor. Apart from some fine-tuned examples such as the linear, and quadratic potential ${U}(\Phi)$ in the Jordan frame, it is true that intermediate accelerated solutions are generic late-time attractors in a modified JBD theory.

The inflationary paradigm is the most successful model that explains the observed spectrum of primordial perturbations. However, the precise emergence of such inhomogeneities and the quantum-to-classical transition of the perturbations has not yet reached a consensus among the community. The Continuous Spontaneous Localization model (CSL), in the cosmological context, might be used to provide a solution to the mentioned issues by considering a dynamical reduction of the wave function. The CSL model has been applied to the inflationary universe before and different conclusions have been obtained. In this letter, we use a different approach to implement the CSL model during inflation. In particular, in addition to accounting for the quantum-to-classical transition, we use the CSL model to generate the primordial perturbations, that is, the dynamical evolution provided by the CSL model is responsible for the transition from a homogeneous and isotropic initial state to a final one lacking such symmetries. Our approach leads to results that can be clearly distinguished from preceding works. Specifically, the scalar and tensor power spectra are not time-dependent, and retains the amplification mechanism of the CSL model. Moreover, our framework depends only on one parameter (the CSL parameter) and its value is consistent with cosmological and laboratory observations.

The inflationary paradigm is an important cornerstone of the concordance cosmological model. However, standard inflation cannot fully address the transition from an early homogeneous and isotropic stage, to another one lacking such symmetries corresponding to our present universe. In previous works, a self-induced collapse of the wave function has been suggested as the missing ingredient of inflation. Most of the analysis regarding the collapse hypothesis has been solely focused on the characteristics of the spectrum associated to scalar perturbations, and within a semiclassical gravity framework. In this Letter, working in terms of a joint metric-matter quantization for inflation, we calculate, for the first time, the tensor power spectrum and the tensor-to-scalar ratio corresponding to the amplitude of primordial gravitational waves resulting from considering a generic self-induced collapse.

We investigate spherically symmetric cosmological models in Einstein-aether theory with a tilted (non-comoving) perfect fluid source. We use a 1+3 frame formalism and adopt the comoving aether gauge to derive the evolution equations, which form a well-posed system of first order partial differential equations in two variables. We then introduce normalized variables. The formalism is particularly well-suited for numerical computations and the study of the qualitative properties of the models, which are also solutions of Horava gravity. We study the local stability of the equilibrium points of the resulting dynamical system corresponding to physically realistic inhomogeneous cosmological models and astrophysical objects with values for the parameters which are consistent with current constraints. In particular, we consider dust models in ($\beta-$) normalized variables and derive a reduced (closed) evolution system and we obtain the general evolution equations for the spatially homogeneous Kantowski-Sachs models using appropriate bounded normalized variables. We then analyse these models, with special emphasis on the future asymptotic behaviour for different values of the parameters. Finally, we investigate static models for a mixture of a (necessarily non-tilted) perfect fluid with a barotropic equations of state and a scalar field.

We investigate the cosmological behavior in a universe governed by time asymmetric extensions of general relativity, which is a novel modified gravity based on the addition of new, time-asymmetric, terms on the Hamiltonian framework, in a way that the algebra of constraints and local physics remain unchanged. Nevertheless, at cosmological scales these new terms can have significant effects that can alter the universe evolution, both at early and late times, and the freedom in the choice of the involved modification function makes the scenario able to produce a huge class of cosmological behaviors. For basic ansatzes of modification, we perform a detailed dynamical analysis, extracting the stable late-time solutions. Amongst others, we find that the universe can result in dark-energy dominated, accelerating solutions, even in the absence of an explicit cosmological constant, in which the dark energy can be quintessence-like, phantom-like, or behave as an effective cosmological constant. Moreover, it can result to matter-domination, or to a Big Rip, or experience the sequence from matter to dark energy domination. Additionally, in the case of closed curvature, the universe may experience a cosmological bounce or turnaround, or even cyclic behavior. Finally, these scenarios can easily satisfy the observational and phenomenological requirements. Hence, time asymmetric cosmology can be a good candidate for the description of the universe.

"The self-induced collapse hypothesis" was introduced by D. Sudarsky and collaborators to explain the origin of cosmic structure from a perfect isotropic and homogeneous universe during the inflationary regime. In this paper, we calculate the power spectrum for the tensor modes, within the semiclassical gravity approximation, with the additional hypothesis of a generic self-induced collapse of the inflaton's wave function; we also compute an estimate for the tensor-to-scalar ratio. Based on this calculation, we show that the considered proposal exhibits a strong suppression of the tensor modes amplitude; nevertheless, the corresponding amplitude is still consistent with the joint BICEP/KECK and Planck Collaboration's limit on the tensor-to-scalar ratio.

The statistical properties of the primordial density perturbations has been considered in the past decade as a powerful probe of the physical processes taking place in the early universe. Within the inflationary paradigm, the properties of the bispectrum are one of the keys that serves to discriminate among competing scenarios concerning the details of the origin of cosmological perturbations. However, all of the scenarios, based on the conventional approach to the so-called "quantum-to-classical transition" during inflation, lack the ability to point out the precise physical mechanism responsible for generating the inhomogeneity and anisotropy of our universe starting from and exactly homogeneous and isotropic vacuum state associated with the early inflationary regime. In past works, we have shown that the proposals involving a spontaneous dynamical reduction of the quantum state provide plausible explanations for the birth of said primordial inhomogeneities and anisotropies. In the present manuscript we show that, when considering within the context of such proposals, the characterization of the spectrum and bispectrum turn out to be quite different from those found in the traditional approach, and in particular, some of the statistical features, must be treated in a different way leading to some rather different conclusions.

We study higher derivative terms associated with scalar field cosmology. We consider a coupling between the scalar field and the geometry inspired by the Pais-Uhlenbeck oscillator, given by $\alpha\partial_{\mu}\partial^{\mu}\phi\partial_{\nu}\partial^{\nu}\phi.$ We investigate the cosmological dynamics in a phase space. For $\alpha>0$, we provide conditions for the stability of de Sitter solutions. In this case the crossing of the phantom divide $w_{DE}=-1$ occurs once; thereafter, the equation of state parameter remains under this line, asymptotically reaching towards the de Sitter solution from below. For $\alpha<0,$ which is the portion of the parameter space where in addition to crossing the phantom divide, cyclic behavior is possible, we present regions in the parameter space where, according to Smilga's classification the ghost has benign or malicious behavior.

The primordial bispectrum has been considered in the past decade as a powerful probe of the physical processes taking place in the early Universe. Within the inflationary paradigm, the properties of the bispectrum are one of the keys that serves to discriminate among competing scenarios concerning the details of the origin of cosmological perturbations. However, all of the scenarios, based on the conventional approach to the so-called "quantum-to-classical transition" during inflation, lack the ability to point out the precise physical mechanism responsible for generating the inhomogeneity and anisotropy of our Universe starting from and exactly homogeneous and isotropic vacuum state associated with the early inflationary regime. In past works, we have shown that the proposals involving a spontaneous dynamical reduction of the quantum state provide plausible explanations for the birth of said primordial inhomogeneities and anisotropies. In the present letter, we show that, when considering single-field slow-roll inflation within the context of such proposals, the expected characteristics of the bispectrum turn out to be quite different from those found in the traditional approach. In particular, the statistical features corresponding to the primordial perturbations, which are normally associated with the bispectrum, are treated here in a novel way leading to rather different conclusions.

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.

The $f(T,T_G)$ class of gravitational modification, based on the quadratic torsion scalar $T$, as well as on the new quartic torsion scalar $T_G$ which is the teleparallel equivalent of the Gauss-Bonnet term, is a novel theory, different from both $f(T)$ and $f(R,G)$ ones. We perform a detailed dynamical analysis of a spatially flat universe governed by the simplest non-trivial model of $f(T,T_G)$ gravity which does not introduce a new mass scale. We find that the universe can result in dark-energy dominated, quintessence-like, cosmological-constant-like or phantom-like solutions, according to the parameter choices. Additionally, it may result to a dark energy - dark matter scaling solution, and thus it can alleviate the coincidence problem. Finally, the analysis "at infinity" reveals that the universe may exhibit future, past, or intermediate singularities depending on the parameters.

In this paper we investigate, from the dynamical systems perspective, the evolution of a Kantowski-Sachs metric in a generic class of $f(R)$ models. We present conditions (i. e., differentiability conditions, existence of minima, monotony intervals, etc.) for a free input function related to the $f(R)$, that guarantee the asymptotic stability of well-motivated physical solutions, specially, self-accelerated solutions, allowing to describe both inflationary- and late-time acceleration stages of the cosmic evolution. We discuss which $f(R)$ theories allows for a cosmic evolution with an acceptable matter era, in correspondence to the modern cosmological paradigm. 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 phantom solutions. Additionally, we find that a cosmological bounce and turnaround are realized in a part of the parameter-space as a consequence of the metric choice.

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.

We construct general anisotropic cosmological scenarios governed by an $f(R)=R^n$ gravitational sector. Focusing then on some specific geometries, and modelling the matter content as a perfect fluid, we perform a phase-space analysis. We analyze the possibility of accelerating expansion at late times, and additionally, we determine conditions for the parameter $n$ for the existence of phantom behavior, contracting solutions as well as of cyclic cosmology. Furthermore, we analyze if the universe evolves towards the future isotropization without relying on a cosmic no-hair theorem. Our results indicate that anisotropic geometries in modified gravitational frameworks present radically different cosmological behaviors compared to the simple isotropic scenarios.

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.

The usual account for the origin of cosmic structure during inflation is not fully satisfactory, as it lacks a physical mechanism capable of generating the inhomogeneity and anisotropy of our Universe, from an exactly homogeneous and isotropic initial state associated with the early inflationary regime. The proposal in [A. Perez, H. Sahlmann, and D. Sudarsky, Classical Quantum Gravity, 23, 2317, (2006)] considers the spontaneous dynamical collapse of the wave function, as a possible answer to that problem. In this work, we review briefly the difficulties facing the standard approach, as well as the answers provided by the above proposal and explore their relevance to the investigations concerning the characterization of the primordial spectrum and other statistical aspects of the cosmic microwave background and large-scale matter distribution. We will see that the new approach leads to novel ways of considering some of the relevant questions, and, in particular, to distinct characterizations of the non-Gaussianities that might have left imprints on the available data.

We perform a detailed dynamical analysis of various cosmological scenarios in extended (varying-mass) nonlinear massive gravity. Due to the enhanced freedom in choosing the involved free functions, this cosmological paradigm allows for a huge variety of solutions that can attract the universe at late times, comparing to scalar-field cosmology or usual nonlinear massive gravity. Amongst others, it accepts quintessence, phantom, or cosmological-constant-like late-time solutions, which moreover can alleviate the coincidence problem. These features seem to be general and non-sensitive to the imposed ansantzes and model parameters, and thus extended nonlinear massive gravity can be a good candidate for the description of nature.

We perform a detailed dynamical analysis of generalized Galileon cosmology, incorporating also the requirements of ghost and instabilities absence. We find that there are not any new stable late-time solutions apart from those of standard quintessence. Furthermore, depending on the model parameters the Galileons may survive at late times or they may completely disappear by the dynamics, however the corresponding observables are always independent of the Galileon terms, determined only by the usual action terms. Thus, although the Galileons can play an important role at inflationary or at recent times, in the future, when the universe will asymptotically reach its stable state, they will not have any effect on its evolution.

We construct general anisotropic cosmological scenarios governed by an $f(R)$ gravitational sector. Focusing then on Kantowski-Sachs geometries in the case of $R^n$-gravity, and modelling the matter content as a perfect fluid, we perform a detailed phase-space analysis. We find that at late times the universe can result to a state of accelerating expansion, and additionally, for a particular $n$-range ($2<n<3$) it exhibits phantom behavior. Furthermore, isotropization has been achieved independently of the initial anisotropy degree, showing in a natural way why the observable universe is so homogeneous and isotropic, without relying on a cosmic no-hair theorem. Moreover, contracting solutions have also a large probability to be the late-time states of the universe. Finally, we can also obtain the realization of the cosmological bounce and turnaround, as well as of cyclic cosmology. These features indicate that anisotropic geometries in modified gravitational frameworks present radically different cosmological behaviors comparing to the simple isotropic scenarios.

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).

We survey the application of specific tools to distinguish amongst the wide variety of dark energy models that are nowadays under investigation. The first class of tools is more mathematical in character: the application of the theory of dynamical systems to select the better behaved models, with appropriate attractors in the past and future. The second class of tools is rather physical: the use of astrophysical observations to crack the degeneracy of classes of dark energy models. In this last case the observations related with structure formation are emphasized both in the linear and non-linear regimes. We exemplify several studies based on our research, such as quintom and quinstant dark energy ones. Quintom dark energy paradigm is a hybrid construction of quintessence and phantom fields, which does not suffer from fine-tuning problems associated to phantom field and additionally it preserves the scaling behavior of quintessence. Quintom dark energy is motivated on theoretical grounds as an explanation for the crossing of the phantom divide, i.e. the smooth crossing of the dark energy state equation parameter below the value -1. On the other hand, quinstant dark energy is considered to be formed by quintessence and a negative cosmological constant, the inclusion of this later component allows for a viable mechanism to halt acceleration. We comment that the quinstant dark energy scenario gives good predictions for structure formation in the linear regime, but fails to do that in the non-linear one, for redshifts larger than one. We comment that there might still be some degree of arbitrariness in the selection of the best dark energy models.

We perform a detailed phase-space analysis of Horava-Lifshitz cosmology, with and without the detailed-balance condition. Under detailed-balance we find that the universe can reach a bouncing-oscillatory state at late times, in which dark-energy, behaving as a simple cosmological constant, is dominant. In the case where the detailed-balance condition is relaxed, we find that the universe reaches an eternally expanding, dark-energy-dominated solution, with the oscillatory state preserving also a small probability. Although this analysis indicates that Horava-Lifshitz cosmology can be compatible with observations, it does not enlighten the discussion about its possible conceptual and theoretical problems.

We investigate several varying-mass dark-matter particle models in the framework of phantom cosmology. We examine whether there exist late-time cosmological solutions, corresponding to an accelerating universe and possessing dark energy and dark matter densities of the same order. Imposing exponential or power-law potentials and exponential or power-law mass dependence, we conclude that the coincidence problem cannot be solved or even alleviated. Thus, if dark energy is attributed to the phantom paradigm, varying-mass dark matter models cannot fulfill the basic requirement that led to their construction.

We explore quintessence models of dark energy which exhibit non-minimal coupling between the dark matter and the dark energy components of the cosmic fluid. The kind of coupling chosen is inspired in scalar-tensor theories of gravity. We impose a suitable dynamics of the expansion allowing to derive exact Friedmann-Robertson-Walker solutions once the coupling function is given as input. Self-interaction potentials of single and double exponential types emerge as result of our choice of the coupling function. The stability and existence of the solutions is discussed in some detail. Although, in general, models with appropriated interaction between the components of the cosmic mixture are useful to handle the coincidence problem, in the present study the coincidence can not be evaded due to the choice of the solution generating ansatz.

In this paper we analyze the asymptotic behavior of Cardassian cosmological models filled with a perfect fluid and a scalar field with an exponential potential. Cardassian cosmologies arise from modifications of the Friedmann equation, and among the different proposals within that framework we will choose those of the form $3H^2-\rho\propto \rho^n$ with $n<1$. We construct two three dimensional dynamical systems arising from the evolution equations, respectively adapted for studying the high and low energy limits. Using standard dynamical systems techniques we find the fixed points and characterize the solutions they represent. We pay especial attention to the properties inherent to the modifications and compare with the (standard) unmodified scenario. Among other interesting results, we find there are no late-time tracking attractors.

We explore quintessence models of dark energy which exhibit non-minimal coupling between the dark matter and the dark energy components of the cosmic fluid. The kind of coupling chosen is inspired in scalar-tensor theories of gravity. We impose a suitable dynamics of the expansion allowing to derive exact Friedmann-Robertson-Walker solutions once the coupling function is given as input. Coupling functions that lead to self-interaction potentials of the single and double exponential types are the target of the present investigation. The stability and existence of the solutions is discussed in some detail. The models are tested against the observational evidence.