We present a conjugate analysis of two different dark energy models investigating both their agreement with recent data and their thermodynamical properties. The successful match with the data allows to both constrain the model parameters and characterize their kinematical properties. As a novel step, we exploit the strong connection between gravity and thermodynamics to further check models viability by investigating their thermodynamical quantities. In particular, we study whether the cosmological scenario fulfills the generalized second law of thermodynamics and, moreover, we contrast the two model asking whether the evolution of the total entropy is in agreement with the expectation for a closed system. As a general result, we discuss whether thermodynamic constraints can be a valid complementary way to both constrain dark energy models and differentiate among rival scenarios.

We formally discuss the post-Minkowskian limit of $f(R)$-gravity without adopting conformal transformations but developing all the calculations in the original Jordan frame. It is shown that such an approach gives rise, in general, together with the standard massless graviton, to massive scalar modes whose masses are directly related to the analytic parameters of the theory. In this sense, the presence of massless gravitons only is a peculiar feature of General Relativity. This fact is never stressed enough and could have dramatic consequences in detection of gravitational waves. Finally the role of curvature stress-energy tensor of $f(R)$-gravity is discussed showing that it generalizes the so called Landau-Lifshitz tensor of General Relativity. The further degrees of freedom, giving rise to the massive modes, are directly related to the structure of such a tensor.

We show that any analytic $f(R)$-gravity model, in the metric approach, presents a weak field limit where the standard Newtonian potential is corrected by a Yukawa-like term. This general result has never been pointed out but often derived for some particular theories. This means that only $f(R)=R$ allows to recover the standard Newton potential while this is not the case for other relativistic theories of gravity. Some considerations on the physical consequences of such a general solution are addressed.

Spherical symmetry in $f(R)$ gravity is discussed in details considering also the relations with the weak field limit. Exact solutions are obtained for constant Ricci curvature scalar and for Ricci scalar depending on the radial coordinate. In particular, we discuss how to obtain results which can be consistently compared with General Relativity giving the well known post-Newtonian and post-Minkowskian limits. Furthermore, we implement a perturbation approach to obtain solutions up to the first order starting from spherically symmetric backgrounds. Exact solutions are given for several classes of $f(R)$ theories in both $R =$ constant and $R = R(r)$.

Building on earlier work, we discuss a general framework for exploring the cosmological dynamics of Higher Order Theories of Gravity. We show that once the theory of gravity has been specified, the cosmological equations can be written as a first-order autonomous system and we give several examples which illustrate the utility of our method. We also discuss a number of results which have appeared recently in the literature.

A general analytic procedure is developed to deal with the Newtonian limit of $f(R)$ gravity. A discussion comparing the Newtonian and the post-Newtonian limit of these models is proposed in order to point out the differences between the two approaches. We calculate the post-Newtonian parameters of such theories without any redefinition of the degrees of freedom, in particular, without adopting some scalar fields and without any change from Jordan to Einstein frame. Considering the Taylor expansion of a generic $f(R)$ theory, it is possible to obtain general solutions in term of the metric coefficients up to the third order of approximation. In particular, the solution relative to the $g_{tt}$ component gives a gravitational potential always corrected with respect to the Newtonian one of the linear theory $f(R)=R$. Furthermore, we show that the Birkhoff theorem is not a general result for $f(R)$-gravity since time-dependent evolution for spherically symmetric solutions can be achieved depending on the order of perturbations. Finally, we discuss the post-Minkowskian limit and the emergence of massive gravitational wave solutions.

We give a rigorous and mathematically well defined presentation of the Covariant and Gauge Invariant theory of scalar perturbations of a Friedmann-Lemaitre-Robertson-Walker universe for Fourth Order Gravity, where the matter is described by a perfect fluid with a barotropic equation of state. The general perturbations equations are applied to a simple background solution of R^n gravity. We obtain exact solutions of the perturbations equations for scales much bigger than the Hubble radius. These solutions have a number of interesting features. In particular, we find that for all values of n there is always a growing mode for the density contrast, even if the universe undergoes an accelerated expansion. Such a behaviour does not occur in standard General Relativity, where as soon as Dark Energy dominates, the density contrast experiences an unrelenting decay. This peculiarity is sufficiently novel to warrant further investigation on fourth order gravity models.

We investigate the possibility that the observed flatness of the rotation curves of spiral galaxies is not an evidence for the existence of dark matter haloes, but rather a signal of the breakdown of General Relativity. To this aim, we consider power - law fourth order theories of gravity obtained by replacing the scalar curvature $R$ with $f(R) = f_0 R^n$ in the gravity Lagrangian. We show that, in the low energy limit, the gravitational potential generated by a pointlike source may be written as $\Phi(r) \propto r^{-1} \left [ 1 + (r/r_c)^{\beta} \right ]$ with $\beta$ a function of the slope $n$ of the gravity Lagrangian and $r_c$ a scalelength depending on the gravitating system properties. In order to apply the model to realistic systems, we compute the modified potential and the rotation curve for spherically symmetric and for thin disk mass distributions. It turns out that the potential is still asymptotically decreasing, but the corrected rotation curve, although not flat, is higher than the Newtonian one thus offering the possibility to fit rotation curves without dark matter. To test the viability of the model, we consider a sample of 15 low surface brightness (LSB) galaxies with combined HI and H$\alpha$ measurements of the rotation curve extending in the putative dark matter dominated region. We find a very good agreement between the theoretical rotation curve and the data using only stellar disk and interstellar gas.

PPN-limit of higher order theories of gravity represents a still controversial matter of debate and no definitive answer has been provided, up to now, about this issue. By exploiting the analogy between scalar-tensor and fourth-order theories of gravity, one can generalize the PPN-limit formulation. By using the definition of the PPN-parameters $\gamma$ and $\beta$ in term of the $f(R)$ derivatives, we show that a family of third-order polynomial theories, in the Ricci scalar $R$, turns out to be compatible with the PPN-limit and the deviation from General Relativity theoretically predicted agree with experimental data.

The inflessence model has recently been proposed in an attempt to explain both early inflation and present day accelerated expansion within a single mechanism. The model has been successfully tested against the Hubble diagram of Type Ia Supernovae, the shift parameter, and the acoustic peak parameter. As a further mandatory test, we investigate here structure formation in the inflessence model determining the evolution of matter density contrast $\delta \equiv \delta \rho_M/\rho_M$ in the linear regime. We compare the growth factor $D(a) \equiv \delta/a$ and the growth index $f(z) \equiv d\ln{\delta}/d\ln{a}$ to these same quantities for the successful concordance $\Lambda$CDM model with a particular emphasis on the role of the inflessence parameters $(\gamma, z_Q)$. We also evaluate the anisotropy spectrum of the cosmic microwave background radiation (CMBR) to check whether the inflessence model may be in agreement with the observations. We find that, for large values of $(\gamma, z_Q)$, structure formation proceeds in a similar way to that in the $\Lambda$CDM scenario, and it is also possible to nicely fit the CMBR spectrum.

We show that f(R)-gravity can, in general, give rise to cosmological viable models compatible with a matter-dominated epoch evolving into a late accelerated phase. We discuss the various representations of f(R)-gravity as an ideal fluid or a scalar-tensor gravity theory, taking into account conformal transformations. We point out that mathematical equivalence does not correspond, in several cases, to the physical equivalence of Jordan frame and Einstein frame. Finally, we show that wide classes of f(R)-gravity models, including matter and accelerated phases, can be phenomenologically reconstructed by means of observational data. In principle, any popular quintessence models could be "reframed" as an f(R)-gravity model.

Gravitational lensing is investigated in the weak field limit of fourth order gravity in which the Lagrangian of the gravitational field is modified by replacing the Ricci scalar curvature R with an analytical expression $f(R)$. Considering the case of a pointlike lens, we study the behaviour of the deflection angle in the case of power law Lagrangians, i.e. with f(R) = f_0 R^n. In order to investigate possible detectable signatures, the position of the Einstein ring and the solutions of the lens equation are evaluated considering the change with respect to the standard case. Effects on the amplification of the images and the Paczynski curve in microlensing experiments are also estimated.

Astrophysical observations are pointing out huge amounts of dark matter and dark energy needed to explain the observed large scale structures and cosmic accelerating expansion. Up to now, no experimental evidence has been found, at fundamental level, to explain such mysterious components. The problem could be completely reversed considering dark matter and dark energy as shortcomings of General Relativity and claiming for the correct theory of gravity as that derived by matching the largest number of observational data. As a result, accelerating behavior of cosmic fluid and rotation curves of spiral galaxies are reproduced by means of curvature effects.

We investigate the possibility that part of the dark matter is not made out of the usual cold dark matter (CDM) dustlike particles, but is under the form of a fluid of strings with barotropic factor $w_s = -1/3$ of cosmic origin. To this aim, we split the dark matter density parameter in two terms and investigate the dynamics of a spatially flat universe filled with baryons, CDM, fluid of strings and dark energy, modeling this latter as a cosmological constant or a negative pressure fluid with a constant equation of state $w < 0$. To test the viability of the models and to constrain their parameters, we use the Type Ia Supernovae Hubble diagram and the data on the gas mass fraction in galaxy clusters. We also discuss the weak field limit of a model comprising a significant fraction of dark matter in the form of a fluid of strings and show that this mechanism makes it possible to reduce the need for the elusive and up to now undetected CDM. We finally find that a model comprising both a cosmological constant and a fluid of strings fits very well the data and eliminates the need of phantom dark energy thus representing a viable candidate to alleviate some of the problems plaguing the dark side of the universe.

Abandoning the perfect fluid hypothesis, we investigate here the possibility that the dark energy equation of state (EoS) $w$ is a nonlinear function of the energy density $\rho$. To this aim, we consider four different EoS describing classical fluids near thermodynamical critical points and discuss the main features of cosmological models made out of dust matter and a dark energy term with the given EoS. Each model is tested against the data on the dimensionless coordinate distance to Type Ia Supernovae and radio galaxies, the shift and the acoustic peak parameters and the positions of the first three peaks in the anisotropy spectrum of the comic microwave background radation. We propose a possible interpretation of each model in the framework of scalar field quintessence determining the shape of the self interaction potential $V(\phi)$ that gives rise to each one of the considered thermodynamical EoS. As a general result, we demonstrate that replacing the perfect fluid EoS with more generar expressions gives both the possibility of successfully solving the problem of cosmic acceleration escaping the resort to phantom models.

Based on the {\it dynamical} equivalence between higher order gravity and scalar-tensor gravity the PPN-limit of fourth order gravity is discussed. We exploit this analogy developing a fourth order gravity version of the Eddington PPN-parameters. As a result, Solar System experiments can be reconciled with higher order gravity, if physical constraints descending from experiments are fulfilled.

A detailed analysis of dynamics of cosmological models based on $R^{n}$ gravity is presented. We show that the cosmological equations can be written as a first order autonomous system and analyzed using the standard techniques of dynamical system theory. In absence of perfect fluid matter, we find exact solutions whose behavior and stability are analyzed in terms of the values of the parameter $n$. When matter is introduced, the nature of the (non-minimal) coupling between matter and higher order gravity induces restrictions on the allowed values of $n$. Selecting such intervals of values and following the same procedure used in the vacuum case, we present exact solutions and analyze their stability for a generic value of the parameter $n$. From this analysis emerges the result that for a large set of initial conditions an accelerated expansion is an attractor for the evolution of the $R^n$ cosmology. When matter is present a transient almost-Friedman phase can also be present before the transition to an accelerated expansion.

We present a phenomenologically motivated model which is able to give rise both to an inflationary epoch and to the present day cosmic acceleration. We introduce an approach where the energy density depends on the scale factor $a$ in such a way that a smooth transition from the inflation to the matter dominated era and then to the nowadays accelerated expansion is achieved. We use the dimensionless coordinate distance $y(z)$ to the Gold Type Ia Supernovae dataset and to a sample comprising 20 radio galaxies, the shift parameter ${\cal{R}}$ and the acoustic peak parameter ${\cal{A}}$ to test whether the model is in agreement with observations and to constrain its main parameters. As an independent cross check, we also compare model predictions with the lookback time to galaxy clusters and the age of the universe. Although phenomenologically inspired, the model may be theoretically motivated either resorting to a scenario with a scalar and phantom fields (eventually interacting) or as an effective description arising from higher order gravity theories.

We discuss an alternative approach to quintessence modifying the usual equation of state of the cosmological fluid in order to see if going further than the approximation of perfect fluid allows to better reproduce the available data. We consider a cosmological model comprising only two fluids, namely baryons (modelled as dust) and dark matter with a Van der Waals equation of state. First, the general features of the model are presented and then the evolution of the energy density, the Hubble parameter and the scale factor are determined showing that it is possible to obtain accelerated expansion choosing suitably the model parameters. We use the estimated age of the universe and the data on the dimensionless coordinate distances to Type Ia supernovae and distant radio galaxies to see whether Van der Waals quintessence is viable to explain dark energy and to constrain its parameters.

Higher order theories of gravity have recently attracted a lot of interest as alternative candidates to explain the observed cosmic acceleration without the need of introducing any scalar field. A critical ingredient is the choice of the function f(R) of the Ricci scalar curvature entering the gravity Lagrangian and determining the dynamics of the universe. We describe an efficient procedure to reconstruct f(R) from the Hubble parameter $H$ depending on the redshift z. Using the metric formulation of f(R) theories, we derive a third order linear differential equation for f(R(z)) which can be numerically solved after setting the boundary conditions on the basis of physical considerations. Since H(z) can be reconstructed from the astrophysical data, the method we present makes it possible to determine, in principle, what is the f(R) theory which best reproduces the observed cosmological dynamics. Moreover, the method allows to reconcile dark energy models with f(R) theories finding out what is the expression of f(R) which leads to the same H(z) of the given quintessence model. As interesting examples, we consider "quiessence" (dark energy with constant equation of state) and the Chaplygin gas.

Higher order curvature gravity has recently received a lot of attention due to the fact that it gives rise to cosmological models which seem capable of solving dark energy and quintessence issues without using "ad hoc" scalar fields. Such an approach is naturally related to fundamental theories of quantum gravity which predict higher order terms for loop expansions of quantum fields in curved spacetimes. In this framework, we obtain a class of cosmological solutions which are fitted against cosmological data. We reproduce reliable models able to fit high redshift supernovae and WMAP observations. The age of the universe and other cosmological parameters are recovered in this context. Furthermore, in the weak field limit, we obtain gravitational potentials which differ from the Newtonian one because of repulsive corrections increasing with distance. We evaluate the rotation curve of our Galaxy and compare it with the observed data in order to test the viability of these theories and to estimate the scale-length of the correction. It is remarkable that the Milky Way rotation curve is well fitted without the need of any dark matter halo and similar results hold also for other galaxies.

A phenomenological approach is proposed to the problem of universe accelerated expansion and of the dark energy nature. A general class of models is introduced whose energy density depends on the redshift $z$ in such a way that a smooth transition among the three main phases of the universe evolution (radiation era, matter domination, asymptotical de Sitter state) is naturally achieved. We use the estimated age of the universe, the Hubble diagram of Type Ia Supernovae and the angular size - redshift relation for compact and ultracompact radio structures to test whether the model is in agreement with astrophysical observation and to constrain its main parameters. Although phenomenologically motivated, the model may be straightforwardly interpreted as a two fluids scenario in which the quintessence is generated by a suitably chosen scalar field potential. On the other hand, the same model may also be read in the context of unified dark energy models or in the framework of modified Friedmann equation theories.

Astrophysical observations provide a picture of the universe as a 4-dim homogeneous and isotropic flat space-time dominated by an unknown form of dark energy. To achieve such a cosmology one has to consider in the early universe an inflationary era able to overcome problems of standard cosmological models. Here an inhomogeneous model is proposed which allows to obtain a Friedmann-Robertson-Walker behaviour far away from the inhomogeneities and it naturally describes structures formation. We also obtain that the cosmological term does not prevent structure formation, avoiding a fine tuning problem in initial conditions. The asymptotic exact solution have been calculated. A simple test with universe age prediction has been performed. A relation between the inhomogeneity, the breaking of time reversal, parity and the matter-antimatter asymmetry is briefly discussed.

Quintessence issues can be achieved by taking into account higher order curvature invariants into the effective action of gravitational field. Such an approach is naturally related to fundamental theories of quantum gravity which predict higher order terms in loop expansion of quantum fields in curved space-times. In this framework, we obtain a class of cosmological solutions which are fitted against cosmological data. We reproduce encouraging results able to fit high redshift supernovae and WMAP observations. The age of the universe and other cosmological parameters are discussed in this context.

The issues of quintessence and cosmic acceleration can be discussed in the framework of theories which do not include necessarily scalar fields. It is possible to define pressure and energy density for new components considering effective theories derived from fundamental physics like the extended theories of gravity or simply generalizing the state equation of matter. Exact accelerated expanding solutions can be achieved in several schemes: either in models containing higher order curvature and torsion terms or in models where the state equation of matter is corrected by a second order Van der Waals terms. In this review, we present such new approaches and compare them with observations.

Taking into account a torsion field gives rise to a negative pressure contribution in cosmological dynamics and then to an accelerated behaviour of Hubble fluid. The presence of torsion has the same effect of a Lambda - term. We obtain a general exact solution which well fits the data coming from high - redshift supernovae and Sunyaev - Zeldovich/X - ray method for the determination of the cosmological parameters. On the other hand, it is possible to obtain observational constraints on the amount of torsion density. A dust dominated Friedmann behaviour is recovered as soon as torsion effects are not relevant.