We discuss in detail the weak field limit of f(R,G) gravity taking into account analytic functions of the Ricci scalar R and the Gauss-Bonnet invariant G. Specifically, we develop, in metric formalism, the Newtonian, Post Newtonian and Parameterized Post Newtonian limits starting from general f(R, G) Lagrangian. The special cases of f(R) and f(G) gravities are considered. In the case of the Newtonian limit of f(R, G) gravity, a general solution in terms of Green's functions is achieved.

A generic feature of viable exponential $F(R)$-gravity is investigated. An additional modification to stabilize the effective dark energy oscillations during matter era is proposed and applied to two viable models. An analysis on the future evolution of the universe is performed. Furthermore, a unified model for early and late-time acceleration is proposed and studied.

A generic feature of viable $F(R)$ gravity is investigated: It is demonstrated that during the matter dominated era the large frequency oscillations of the effective dark energy may influence the behavior of higher derivatives of the Hubble parameter with the risk to produce some singular unphysical solutions at high redshift. This behavior is explicitly analyzed for realistic $F(R)$ models, in particular, exponential gravity and a power form model. To stabilize such oscillations, we consider the additional modification of the models via a correction term which does not destroy the viability properties. A detailed analysis on the future evolution of the universe and the evolution history of the growth index of the matter density perturbations are performed. Furthermore, we explore two applications of exponential gravity to the inflationary scenario. We show how it is possible to obtain different numbers of $e$-folds during the early-time acceleration by making different choices of the model parameters in the presence of ultrarelativistic matter, which destabilizes inflation and eventually leads to the exit from the inflationary stage. We execute the numerical analysis of inflation in two viable exponential gravity models. It is proved that at the end of the inflation, the effective energy density and curvature of the universe decrease and thus a unified description between inflation and the $\Lambda$CDM-like dark energy dominated era can be realized.

Cosmological models with an SU(2) Yang-Mills field are studied. For a specific model with a minimally coupled Yang-Mills Lagrangian, which includes an arbitrary function of the second-order term and a fourth-order term, a corresponding reconstruction program is proposed. It is shown that the model with minimal coupling has no de Sitter solutions, for any nontrivial function of the second-order term. To get de Sitter solutions, a gravitational model with nonminimally coupled Yang-Mills fields is then investigated. It is shown that the model with non-minimal coupling has in fact a de Sitter solution, even in absence of the cosmological constant term.

A variant of the accelerating cosmology reconstruction program is developed for $f(R)$ gravity and for a modified Yang-Mills/Maxwell theory. Reconstruction schemes in terms of e-foldings and by using an auxiliary scalar field are developed and carefully compared, for the case of $f(R)$ gravity. An example of a model with a transient phantom behavior without real matter is explicitly discussed in both schemes. Further, the two reconstruction schemes are applied to the more physically interesting case of a Yang-Mills/Maxwell theory, again with explicit examples. Detailed comparison of the two schemes of reconstruction is presented also for this theory. It seems to support, as well, physical non-equivalence of the two frames.