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.

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.

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.

We revise the conditions for the physical viability of a cosmological model in which dark matter has bulk viscosity and also interacts with dark energy. We have also included radiation and baryonic matter components; all matter components are represented by perfect fluids, except the dark matter, that is treated as an imperfect fluid. We impose upon the model the condition of a complete cosmological dynamics that results in an either null or negative bulk viscosity, but the latter also disagrees with the Local Second Law of Thermodynamics. The model is also compared with cosmological observations at different redshifts: type Ia supernova, the shift parameter of CMB, the acoustic peak of BAO, and the Hubble parameter H(z). In general, observations consistently point out to a negative value of the bulk viscous coefficient, and in overall the fitting procedure shows no preference for the model over the standard LCDM model.

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.

In this work we investigate the evolution of matter density perturbations for quintessence models with a self-interaction potential that is a combination of exponentials. One of the models is based on the Einstein theory of gravity, while the other is based on the Brans-Dicke scalar tensor theory. We constrain the parameter space of the models using the determinations for the growth rate of perturbations derived from data of the 2-degree Field Galaxy Redshift Survey.