The introduction of a kind of exponential damping in the correction to General Relativity gives promising results in the construction of viable dark energy model in the context of modified gravity. The model behaves as $R-2\Lambda$ at large curvature and tends to zero at $R\rightarrow 0$, containing flat spacetime solution and implying that the curvature effect that induces the accelerated expansion is unrelated to quantum vacuum energy in flat space-time. The cosmic evolution of the main density parameters in this model is consistent with current observations with an equation of state very close to $-1$ and the effective equation of state showing the transition deceleration-acceleration at $z\sim 0.5$. Apart from viable cosmology, satisfying cosmological and local gravity restrictions, the model may also show measurable differences with $\Lambda$CDM at recent times. Some features of the present model that highlight differences with the Hu-Sawicki model are discussed.

The role of an exponential function of the scalar curvature in the modified gravity is analyzed. Two models are proposed. A toy model that complies with local and cosmological constraints and gives appropriate qualitative description of the cosmic evolution. This model contains a saddle matter-dominant critical point that can evolve towards a late time de Sitter attractor. Initial conditions have been proposed, showing that this toy model has an acceptable matter era and gives an approximate qualitative behavior of cosmic evolution. A second viable model, behaves very close to $\Lambda$CDM at early times and can satisfy local and cosmological constraints. It behaves as $R-2\Lambda$ at $R\rightarrow \infty$ and tends to zero at $R\rightarrow 0$, containing flat spacetime solution. The model gives viable cosmological trajectories that, as the first model, connect the matter dominated point with a late time de Sitter attractor. The cosmic evolution of the main density parameters in this model is consistent with current observations with an equation of state very close to $-1$.

We consider new models of dark energy with finite time future singularities, by introducing the pressure density as a function of the scale factor. This approach gives acceptable phenomenological models of dark energy, practically indistinguishable from the cosmological constant up to the present, which face future singularities at finite time and finite scale factor. Exact scalar field model representation was found for quintessence, Big Rip and type III singularity models. The simple form of the equation of state allows to establish a relationship between its current value, $w_0$, and the time or redshift at which the singularity takes place. The effect on the growth of matter perturbations was calculated.

The main characteristic of the dark energy is its negative pressure. In a homogeneous and isotropic FRW background, we consider several models for the dark energy fluid, which lead to finite time future singularities of the type I-IV, by introducing the pressure density as a function of the scale factor. This approach gives acceptable behavior of the dark energy equation of state. We give various numerical examples of models with type I-IV singularities, that show very similar late time behavior, making it difficult to determine the type of singularity that would take place in the future.

We consider a scalar-tensor model of dark energy with kinetic and Gauss Bonnet couplings. We study the conditions for the existence of quintessential and phantom power-law expansion, and also analyze these conditions in absence of potential (closely related to string theory). A mechanism to avoid the Big Rip singularity in various asymptotic limits of the model has been studied. It was found that the kinetic and Gauss-Bonnet couplings might prevent the Big Rip singularity in a phantom scenario. The autonomous system for the model has been used to study the stability properties of the power-law solution, and the centre manifold analysis was used to treat zero eigenvalues.