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.

The Schwarzschild-de Sitter (SdS) solution exists in the large majority of modified gravity theories, as expected, and in particular the effective cosmological constant is determined by the specific parameters of the given theory. We explore the possibility to use future extended radio-tracking data from the currently ongoing New Horizons mission in the outskirts peripheries of the Solar System, at about 40 au, in order to constrain this effective cosmological constant, and thus to impose constrain on each scenario's parameters. We investigate some of the recently most studied modified gravities, namely $f(R)$ and $f(T)$ theories, dRGT massive gravity, and Ho\v{r}ava-Lifshitz gravity, and we show that New Horizons mission may bring an improvement of one-two orders of magnitude with respect to the present bounds from planetary orbital dynamics.

In the framework of $f(T)$ theories of gravity, we solve the field equations for $f(T)=T+\alpha T^{n}$, in the weak-field approximation and for spherical symmetry spacetime. Since $f(T)=T$ corresponds to Teleparallel Gravity, which is equivalent to General Relativity, the non linearity of the Lagrangian are expected to produce perturbations of the general relativistic solutions, parameterized by $\alpha$. Hence, we use the $f(T)$ solutions to model the gravitational field of the Sun, and exploit data from accurate tracking of spacecrafts orbiting Mercury and Saturn to infer preliminary insights on what could be obtained about the model parameter $\alpha$ and the cosmological constant $\Lambda$. It turns out that improvements of about one-three orders with respect to the present-day constraints in the literature of magnitude seem possible.

We study weak-field solutions having spherical symmetry in $f(T)$ gravity; to this end, we solve the field equations for a non diagonal tetrad, starting from Lagrangian in the form $f(T)=T+\alpha T^{n}$, where $\alpha$ is a small constant, parameterizing the departure of the theory from GR. We show that the classical spherically symmetric solutions of GR, i.e. the Schwarzschild and Schwarzschild-de Sitter solutions, are perturbed by terms in the form $\propto r^{2-2n}$ and discuss the impact of these perturbations in observational tests.

The Kantowski-Sachs cosmological model sourced by a Skyrme field and a cosmological constant is considered in the framework of General Relativity. Assuming a constant radial profile function for the hedgehog ansatz, the Skyrme contribution to Einstein equations is shown to be equivalent to an anisotropic fluid. Using dynamical system techniques, a qualitative analysis of the cosmological equations is presented. Physically interesting features of the model such as isotropization, bounce and recollapse are discussed.

Alterations of the gravity Lagrangian introduced in modified torsion gravity theories---also referred to as $f(T)$ gravity---allows for an accelerated expansion in a matter dominated Universe. In this framework, the cosmic speed up is driven by an effective `torsion fluid'. Besides the background evolution of the Universe, structure formation is also modified because of a time dependent effective gravitational constant. Here, we investigate the imprints of $f(T)$ gravity on galaxy clustering and weak gravitational lensing to the aim of understanding whether future galaxy surveys could constrain torsion gravity and discriminate amongst it and standard general relativity. Specifically, we compute Fisher matrix forecasts for two viable $f(T)$ models to both infer the accuracy on the measurement of the model parameters and evaluate the power that a combined clustering and shear analysis will have as a tool for model selection. We find that with such a combination of probes it will indeed be possible to tightly constrain $f(T)$ model parameters. Moreover, the Occam's razor provided by the Bayes factor will allow us to confirm an $f(T)$ power-law extension of the concordance $\Lambda{\rm CDM}$ model, were a value larger than 0.02 of its power-law slope measured, whereas in $\Lambda{\rm CDM}$ it is exactly 0.

The origin of the accelerated expansion of the universe is still unclear and new physics is needed on cosmological scales. We propose and test a novel interpretation of dark energy as originated by an elastic strain due to a cosmic defect in an otherwise Euclidean space-time. The strain modifies the expansion history of the universe. This new effective contribution tracks radiation at early times and mimics a cosmological constant at late times. The theory is tested against observations, from nucleosynthesis to the cosmic microwave background and formation and evolution of large scale structure to supernovae. Data are very well reproduced with Lam\'e parameters of the order of 10^{-52} m^{-2}.

We consider static cosmological solutions along with their stability properties in the framework of a recently proposed theory of massive gravity. We show that the modifcation introduced in the cosmological equations leads to several new solutions, only sourced by a perfect fluid, generalizing the Einstein Static Universe found in General Relativity. Using dynamical system techniques and numerical analysis, we show that the found solutions can be either neutrally stable or unstable against spatially homogeneous and isotropic perturbations.

A covariant formulation of a theory with a massive graviton and no negative energy state has been recently proposed as an alternative to the usual General Relativity framework. For a spatially flat homogenous and isotropic universe, the theory introduces modified Friedmann equations where the standard matter term is supplemented by four effective fluids mimicking dust, cosmological constant, quintessence and stiff matter, respectively. We test the viability of this massive gravity formulation by contrasting its theoretical prediction to the Hubble diagram as traced by Type Ia Supernovae (SNeIa) and Gamma Ray Bursts (GRBs), the $H(z)$ measurements from passively evolving galaxies, Baryon Acoustic Oscillations (BAOs) from galaxy surveys and the distance priors from the Cosmic Microwave Background Radiation (CMBR) anisotropy spectrum. It turns out that the model is indeed able to very well fit this large dataset thus offering a viable alternative to the usual dark energy framework. We finally set stringent constraints on its parameters also narrowing down the allowed range for the graviton mass.

Generalised Teleparallel gravity, also referred to as f(T) gravity, has been recently proposed as an extended theory of gravitation able to give rise to an accelerated expansion in a matter only universe. The cosmic speed up is driven by an effective torsion fluid whose equation of state depend on the f(T) function entering the modified gravity Lagrangian. We focus on two particular choices for f(T) which share the nice property to emulate a phantom divide crossing as suggested by some recent data. We check their viability contrasting the predicted background dynamics to the Hubble diagram as traced by both Type Ia Supernovae (SNeIa) and Gamma Ray Bursts (GRBs), the measurement of the rate expansion H(z), the Baryon Acoustic Oscillations (BAOs) at different redshifts, and the Cosmic Microwave Background Radiation (CMBR) distance priors. Both f(T) models turn out to be in very good agreement with this large dataset so that we also investigate whether it is possible to discriminate among them relying on the different growth factors.

The properties of a spherically symmetric static space-time permeated of dark energy are worked out. Dark energy is viewed as the strain energy of an elastically deformable four dimensional manifold. The metric is worked out in the vacuum region around a central spherical mass/defect in the linear approximation. We discuss analogies and differences with the analogue in the de Sitter space time and how these competing scenarios could be differentiated on an observational ground. The comparison with the tests at the solar system scale puts upper limits to the parameters of the theory, consistent with the values obtained applying the classical cosmological tests.

We focus on a modified version of Horava - Lifschitz theory and, in particular, we consider the impact of its weak - field static spherically symmetric limit on the galaxy dynamics. In a previous paper, we used the modified gravitational potential obtained in this theory to evaluate the Milky Way rotation curve using a spheroidal truncated power - law bulge and a double exponential disc as the only sources of the gravitational field and showed that the modified rotation curved is not in agreement with the data. Making a step forward, we here include also the contribution from a dark matter halo in order to see whether this helps fitting the rotation curve data. As a test case, we consider a sample of spiral galaxies with smooth baryon matter distribution and well measured circular velocity profiles. It turns out that, although a marginal agreement with the data can be found, this can only be obtained if the dark matter halo has an unrealistically small virial mass and incredibly large concentration. Such results can be interpreted as a strong evidence against the reliability of the gravitational potential obtained in the modified version of Horava -Lifschitz theory that we consider.

It is argued that the discovery of cosmic acceleration could have been anticipated on thermodynamic grounds, namely, the generalized second law and the approach to equilibrium at large scale factor. Therefore, the existence of dark energy -or equivalently, some modified gravity theory- should have been expected. In general, cosmological models that satisfy the above criteria show compatibility with observational data.

The Cosmic Defect theory has been confronted with four observational constraints: primordial nuclear species abundances emerging from the big bang nucleosynthesis; large scale structure formation in the universe; cosmic microwave background acoustic scale; luminosity distances of type Ia supernovae. The test has been based on a statistical analysis of the a posteriori probabilities for three parameters of the theory. The result has been quite satisfactory and such that the performance of the theory is not distinguishable from the one of the Lambda-CDM theory. The use of the optimal values of the parameters for the calculation of the Hubble constant and the age of the universe confirms the compatibility of the Cosmic Defect approach with observations.

The Horava - Lifshitz (HL) theory has recently attracted a lot of interest as a viable solution to some quantum gravity related problems and the presence of an effective cosmological constant able to drive the cosmic speed up. We show here that, in the weak field limit, the HL proposal leads to a modification of the gravitational potential because of two additive terms (scaling respectively as $r^2$ and $r^{-4}$) to the Newtonian $1/r$ potential. We then derive a general expression to compute the rotation curve of an extended system under the assumption that the mass density only depends on the cylindrical coordinates $(R, z)$ showing that the HL modification induces a dependence of the circular velocity on the mass function which is a new feature of the theory. As a first exploratory analysis, we then try fitting the Milky Way rotation curve using its visible components only in order to see whether the HL modified potential can be an alternative to the dark matter framework. This turns out not to be the case so that we argue that dark matter is still needed, but the amount of dark matter and the dark halo density profile have to be revised according to the new HL potential.

Generalized Einstein - Aether vector field models have been shown to provide, in the weak field regime, modifications to gravity which can be reconciled with the successfull MOND proposal. Very little is known, however, on the function F(K) defining the vector field Lagrangian so that an analysis of the viability of such theories at the cosmological scales has never been performed. As a first step along this route, we rely on the relation between F(K) and the MOND interpolating function $\mu(a/a_0)$ to assign the vector field Lagrangian thus obtaining what we refer to as "MONDian vector models". Since they are able by construction to recover the MOND successes on galaxy scales, we investigate whether they can also drive the observed accelerated expansion by fitting the models to the Type Ia Supernovae data. Should be this the case, we have a unified framework where both dark energy and dark matter can be seen as different manifestations of a single vector field. It turns out that both MONDian vector models are able to well fit the low redshift data on Type Ia Supernovae, while some tension could be present in the high z regime.

The cosmological expansion of the universe affects the behaviour of all physical systems and, in the case of gravitationally bound ones, could correspond to or mimic a time dependent Newton's constant. Here we discuss the case of a locally spherical mass distribution embedded in a generic Robertson Walker universe. Choosing the most appropriate metric tensor for the problem and assuming that the local time scale is much much lower than the cosmic one, we show that G is practically unaffected thus leaving the absolute magnitude of type Ia supernovae unaltered at all epochs.