We propose a new class of gravity theories which are characterized by a nontrivial coupling between the gravitational metric and matter mediated by an auxiliary rank-2 tensor. The actions generating the field equations are constructed so that these theories are equivalent to general relativity in a vacuum, and only differ from general relativity theory within a matter distribution. We analyze in detail one of the simplest realizations of these generalized coupling theories. We show that in this case the propagation speed of gravitational radiation in matter is different from its value in vacuum and that this can be used to weakly constrain the (single) additional parameter of the theory. An analysis of the evolution of homogeneous and isotropic spacetimes in the same framework shows that there exist cosmic histories with both an inflationary phase and a dark era characterized by a different expansion rate.

We compare theoretical predictions with observations for a class of cosmological models in which the dark energy component is modeled as a fermionic condensate, non-minimally coupled with the gravitational field and characterized by some specific self-interaction potentials. Our analysis is based on the Markov Chain Monte Carlo Method (MCMC) and employs different data sets. It turns out that with an appropriate choice of parameters our models are fully compatible with several observed data. We combine these parameter values with phase space analysis to deduce the features of the entire cosmic history of the considered models.

The grand challenges of contemporary fundamental physics---dark matter, dark energy, vacuum energy, inflation and early universe cosmology, singularities and the hierarchy problem---all involve gravity as a key component. And of all gravitational phenomena, black holes stand out in their elegant simplicity, while harbouring some of the most remarkable predictions of General Relativity: event horizons, singularities and ergoregions. The hitherto invisible landscape of the gravitational Universe is being unveiled before our eyes: the historical direct detection of gravitational waves by the LIGO-Virgo collaboration marks the dawn of a new era of scientific exploration. Gravitational-wave astronomy will allow us to test models of black hole formation, growth and evolution, as well as models of gravitational-wave generation and propagation. It will provide evidence for event horizons and ergoregions, test the theory of General Relativity itself, and may reveal the existence of new fundamental fields. The synthesis of these results has the potential to radically reshape our understanding of the cosmos and of the laws of Nature. The purpose of this work is to present a concise, yet comprehensive overview of the state of the art in the relevant fields of research, summarize important open problems, and lay out a roadmap for future progress.

The Nobel Prize winning confirmation in 1998 of the accelerated expansion of our Universe put into sharp focus the need of a consistent theoretical model to explain the origin of this acceleration. As a result over the past two decades there has been a huge theoretical and observational effort into improving our understanding of the Universe. The cosmological equations describing the dynamics of a homogeneous and isotropic Universe are systems of ordinary differential equations, and one of the most elegant ways these can be investigated is by casting them into the form of dynamical systems. This allows the use of powerful analytical and numerical methods to gain a quantitative understanding of the cosmological dynamics derived by the models under study. In this review we apply these techniques to cosmology. We begin with a brief introduction to dynamical systems, fixed points, linear stability theory, Lyapunov stability, centre manifold theory and more advanced topics relating to the global structure of the solutions. Using this machinery we then analyse a large number of cosmological models and show how the stability conditions allow them to be tightly constrained and even ruled out on purely theoretical grounds. We are also able to identify those models which deserve further in depth investigation through comparison with observational data. This review is a comprehensive and detailed study of dynamical systems applications to cosmological models focusing on the late-time behaviour of our Universe, and in particular on its accelerated expansion. In self contained sections we present a large number of models ranging from canonical and non-canonical scalar fields, interacting models and non-scalar field models through to modified gravity scenarios. Selected models are discussed in detail and interpreted in the context of late-time cosmology.

We generalise the covariant Tolman-Oppenheimer-Volkoff equations proposed in arXiv:1709.02818 [gr-qc] to the case of static and spherically symmetric spacetimes with anisotropic sources. The extended equations allow a detailed analysis of the role of the anisotropic terms in the interior solution of relativistic stars and lead to the generalisation of some well known solutions of this type. We show that, like in the isotropic case, one can define generating theorems for the anisotropic Tolman-Oppenheimer-Volkoff equations. We also find that it is possible to define a reconstruction algorithm able to generate a double infinity of interior solutions. Among these, we derive a class of solutions that can represent "quasi-isotropic" stars.

We construct a covariant version of the Tolman-Oppenheimer-Volkoff equations in the case of isotropic sources. The new equations make evident the mathematical problems in the determination of interior solutions of relativistic stellar objects. Using a reconstruction algorithm we find two physically interesting generalisations of previously known stellar interior solutions. The variables that we use also allow an easier formulation of known generating theorems for solutions associated to relativistic stellar objects.

We construct exact solutions representing a Friedmann-Lema\^itre-Robsertson-Walker (FLRW) universe in a generalized hybrid metric-Palatini theory. By writing the gravitational action in a scalar-tensor representation, the new solutions are obtained by either making an ansatz on the scale factor or on the effective potential. Among other relevant results, we show that it is possible to obtain exponentially expanding solutions for flat universes even when the cosmology is not purely vacuum. We then derive the classes of actions for the original theory which generate these solutions.

We explore the idea that the coupling between matter and spacetime is more complex than the one originally envisioned by Einstein. We propose that such coupling takes the form of a new fundamental tensor in the Einstein field equations. We then show that the introduction of this tensor can account for dark phenomenology in General Relativity, maintaining a weak field limit compatible with standard Newtonian gravitation. The same paradigm can be applied any other theory of gravitation. We show, as an example, that in the context of conformal gravity a generalised coupling is able to solve compatibility issues between the matter and the gravitational sector.

We use the 1+1+2 covariant approach to clarify a number of aspects of spherically symmetric solutions of non-minimally coupled scalar tensor theories. Particular attention is focused on the extension of Birkhoff's theorem and the nature of quasi-local horizons in this context.

We discuss dynamical systems approaches and methods applied to flat Robertson-Walker models in $f(R)$-gravity. We argue that a complete description of the solution space of a model requires a global state space analysis that motivates globally covering state space adapted variables. This is shown explicitly by an illustrative example, $f(R) = R + \alpha R^2$, $\alpha > 0$, for which we introduce new regular dynamical systems on global compactly extended state spaces for the Jordan and Einstein frames. This example also allows us to illustrate several local and global dynamical systems techniques involving, e.g., blow ups of nilpotent fixed points, center manifold analysis, averaging, and use of monotone functions. As a result of applying dynamical systems methods to globally state space adapted dynamical systems formulations, we obtain pictures of the entire solution spaces in both the Jordan and the Einstein frames. This shows, e.g., that due to the domain of the conformal transformation between the Jordan and Einstein frames, not all the solutions in the Jordan frame are completely contained in the Einstein frame. We also make comparisons with previous dynamical systems approaches to $f(R)$ cosmology and discuss their advantages and disadvantages.

In this work, we perform a detailed dynamical analysis for the cosmological applications of a nonminimal torsion-matter coupled gravity. Two alternative formalisms are proposed, which enable one to choose between the easier approach for a given problem, and furthermore, we analyze six specific models. In general, we extract fixed points corresponding either to dark-matter dominated, scaling decelerated solutions, or to dark-energy dominated accelerated solutions. Additionally, we find that there is a small parameter region in which the model can experience the transition from the matter epoch to a dark-energy era. These features are in agreement with the observed universe evolution, and make the theory a successful candidate for the description of Nature.

We propose a new dynamical system formalism for the analysis of f(R) cosmologies. The new approach eliminates the need for cumbersome inversions to close the dynamical system and allows the analysis of the phase space of f(R)-gravity models which cannot be investigated using the standard technique. Differently form previously proposed similar techniques, the new method is constructed in such a way to associate to the fixed points scale factors, which contain four integration constants (i.e. solutions of fourth order differential equations). In this way a new light is shed on the physical meaning of the fixed points. We apply this technique to some f(R) Lagrangians relevant for inflationary and dark energy models.

The so called $f(X)$ hybrid metric-Palatini gravity presents a unique viable generalisation of the $f(R)$ theories within the metric-affine formalism. Here the cosmology of the $f(X)$ theories is studied using the dynamical system approach. The method consists of formulating the propagation equation in terms of suitable (expansion-normalised) variables as an autonomous system. The fixed points of the system then represent exact cosmological solutions described by power-law or de Sitter expansion. The formalism is applied to two classes of $f(X)$ models, revealing both standard cosmological fixed points and new accelerating solutions that can be attractors in the phase space. In addition, the fixed point with vanishing expansion rate are considered with special care in order to characterise the stability of Einstein static spaces and bouncing solutions.

We investigate some cosmological models arising from a non-minimal coupling of a fermionic field to gravity in the geometrical setting of Einstein-Cartan-Sciama-Kibble gravity. The role played by the non-minimal coupling together with torsion in facing problems such as cosmological singularity, inflation and dark energy is discussed.

We propose a new reconstruction method for scalar--tensor gravity based on the use of conformal transformations. The new method allows the derivation of a set of interesting exact cosmological solutions in brans Dicke gravity as well as other extensions of General Relativity.

We revisit the relativistic restricted two-body problem with spin employing a perturbation scheme based on Lie series. Starting from a post-Newtonian expansion of the field equations, we develop a first-order secular theory that reproduces well-known relativistic effects such as the precession of the pericentre and the Lense-Thirring and geodetic effects. Additionally, our theory takes into full account the complex interplay between the various relativistic effects, and provides a new explicit solution of the averaged equations of motion in terms of elliptic functions. Our analysis reveals the presence of particular configurations for which non-periodical behaviour can arise. The application of our results to real astrodynamical systems (such as Mercury-like and pulsar planets) highlights the contribution of relativistic effects to the long-term evolution of the spin and orbit of the secondary body.

We propose the most general modified first-order Ho\v{r}ava-Lifshitz (HL) gravity, whose action does not contain time derivatives higher than the second order. The Hamiltonian structure of this theory is studied in all the details in the case of the spatially-flat FRW space-time, demonstrating many of the features of the general theory. It is shown that, with some plausible assumptions, including the projectability of the lapse function, this model is consistent. As a large class of such theories, the modified HL $F(R)$ gravity is introduced. The study of its ultraviolet properties shows that its $z=3$ version seems to be renormalizable in the same way as the original HL proposal. The Hamiltonian analysis of the modified HL $F(R)$ gravity shows that it is in general a consistent theory. The $F(R)$ gravity action is also studied in the fixed-gauge form, where the appearance of a scalar field is particularly illustrative. Then the spatially-flat FRW cosmology for this $F(R)$ gravity is investigated. It is shown that a special choice of parameters for this theory leads to the same equations of motion as in the case of traditional $F(R)$ gravity. Nevertheless, the cosmological structure of the modified HL $F(R)$ gravity turns out to be much richer than for its traditional counterpart. The emergence of multiple de Sitter solutions indicates to the possibility of unification of early-time inflation with late-time acceleration within the same model. Power-law $F(R)$ theories are investigated in detail. It is analytically shown that they have a quite rich cosmological structure: early/late-time cosmic acceleration of quintessence, as well as of phantom types. Also it is demonstrated that all the four known types of finite-time future singularities may occur in the power-law HL $F(R)$ gravity. Finally, a covariant proposal for (renormalizable) $F(R)$ gravity within the HL spirit is presented.

We discuss a mechanism that induces a time-dependent vacuum energy on cosmological scales. It is based on the instability induced renormalization triggered by the low energy quantum fluctuations in a Universe with a positive cosmological constant. We employ the dynamical systems approach to study the qualitative behavior of Friedmann-Robertson-Walker cosmologies where the cosmological constant is dynamically evolving according with this nonperturbative scaling at low energies. It will be shown that it is possible to realize a "two regimes" dark energy phases, where an unstable early phase of power-law evolution of the scale factor is followed by an accelerated expansion era at late times.

The standard model of cosmology is based on homogeneous-isotropic solutions of Einstein's equations. These solutions are known to be gravitationally unstable to local inhomogeneous perturbations, commonly described as evolving on a background given by the same solutions. In this picture, the FLRW backgrounds are taken to describe the average over inhomogeneous perturbations for all times. We study in the present article the (in)stability of FLRW dust backgrounds within a class of averaged inhomogeneous cosmologies. We examine the phase portraits of the latter, discuss their fixed points and orbital structure and provide detailed illustrations. We show that FLRW cosmologies are unstable in some relevant cases: averaged models are driven away from them through structure formation and accelerated expansion. We find support for the proposal that the dark components of the FLRW framework may be associated to these instability sectors. Our conclusion is that FLRW cosmologies have to be considered critically as for their role to serve as reliable models for the physical background.

We discuss the classical tests of general relativity in the presence of Rindler acceleration. Among these tests the perihelion shifts give the tightest constraints and indicate that the Pioneer anomaly cannot be caused by a universal solar system Rindler acceleration. We address potential caveats for massive test-objects. Our tightest bound on Rindler acceleration that comes with no caveats is derived from radar echo delay and yields |a|<3nm/s^2.

We present a new approach of the reconstruction method based on the use of the cosmic parameters instead of a time law for the scale factor. This allows the derivation and analysis of a set of new non-trivial cosmological solutions for $f(R)$-gravity. A number of simple examples are given.

We develop a new covariant formalism to treat spherically symmetric spacetimes in metric} f(R) theories of gravity. Using this formalism we derive the general equations for a static and spherically symmetric metric in a general f(R)-gravity. These equations are used to determine the conditions for which the Schwarzschild metric is the only vacuum solution with vanishing Ricci scalar. We also show that our general framework provides a clear way of showing that the Schwarzschild solution is not a unique static spherically symmetric solution, providing some incite on how the current form of Birkhoff's theorem breaks down for these theories.

We present a framework for the study of lensing in spherically symmetric spacetimes within the context of f(R) gravity. Equations for the propagation of null geodesics, together with an expression for the bending angle are derived for any f(R) theory and then applied to an exact spherically symmetric solution of R^n gravity. We find that for this case more bending is expected for R^n gravity theories in comparison to GR and is dependent on the value of n and the value of distance of closest approach of the incident null geodesic.

The 1+3 covariant approach and the covariant gauge-invariant approach to perturbations are used to analyze in depth conformal transformations in cosmology. Such techniques allow us to obtain very interesting insights on the physical content of these transformations, when applied to non-standard gravity. The results obtained lead to a number of general conclusions on the change of some key quantities describing any two conformally related cosmological models. In particular, it is shown that the physics in the Einstein frame has characteristics which are completely different from those in the Jordan frame. Even if some of the geometrical properties of the cosmology are preserved (homogeneous and isotropic Universes are mapped into homogeneous and isotropic universes), it can happen that decelerating cosmologies are mapped into accelerated ones. Differences become even more pronounced when first-order perturbations are considered: from the 1+3 equations it is seen that first-order vector and tensor perturbations are left unchanged in their structure by the conformal transformation, but this cannot be said of the scalar perturbations, which include the matter density fluctuations. Behavior in the two frames of the growth rate, as well as other evolutionary features, like the presence or absence of oscillations, etc., appear to be different too. The results obtained are then explicitly interpreted and verified with the help of some clarifying examples based on $f(R)$-gravity cosmologies.

Using the dynamical system approach, properties of cosmological models based on the Horava-Lifshitz gravity are systematically studied. In particular, the cosmological phase space of the Horava-Lifshitz model is characterized. The analysis allows to compare some key physical consequences of the imposition (or not) of detailed balance. A result of the investigation is that in the detailed balance case one of the attractors in the theory corresponds to an oscillatory behavior. Such oscillations can be associated to a bouncing universe, as previously described by Brandenberger, and will prevent a possible evolution towards a de Sitter universe. Other results obtained show that the cosmological models generated by Horava-Lifshitz gravity without the detailed balance assumption have indeed the potential to describe the transition between the Friedmann and the dark energy eras. The whole analysis leads to the plausible conclusion that a cosmology compatible with the present observations of the universe can be achieved only if the detailed balance condition is broken.

In this paper we show how the covariant gauge invariant equations for the evolution of scalar, vector and tensor perturbations for a generic $f(R)$-gravity theory can be recast in order to exploit the power of dynamical system methodology. In this way, recent results describing the dynamics of the background FRW model can be easily combined with these equations to reveal important details pertaining to the evolution of cosmological models in fourth order gravity.

We present for the first time the complete matter power spectrum for $R^n$ gravity which has been derived from the fourth order scalar perturbation equations. This leads to the discovery of a characteristic signature of fourth order gravity in the matter power spectrum, the details of which have not seen before in other studies in this area and therefore provides a crucial test for fourth order gravity on cosmological scales.

We investigate the connection between dark energy and fourth order gravity by analyzing the behavior of scalar perturbations around a Friedmann-Robertson-Walker background. The evolution equations for scalar perturbation are derived using the covariant and gauge invariant approach and applied to two widely studied $f(R)$ gravity models. The structure of the general fourth order perturbation equations and the analysis of scalar perturbations lead to the discovery of a characteristic signature of fourth order gravity in the matter power spectrum, the details of which have not seen before in other works in this area. This could provide a crucial test for fourth order gravity on cosmological scales.

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.

We give a rigorous and mathematically clear presentation of the Covariant and Gauge Invariant theory of gravitational waves in a perturbed Friedmann-Lemaitre-Robertson-Walker universe for Fourth Order Gravity, where the matter is described by a perfect fluid with a barotropic equation of state. As an example of a consistent analysis of tensor perturbations in Fourth Order Gravity, we apply the formalism to a simple background solution of R^n gravity. We obtain the exact solutions of the perturbation equations for scales much bigger than and smaller than the Hubble radius. It is shown that the evolution of tensor modes is highly sensitive to the choice of n and an interesting new feature arises. During the radiation dominated era, their exist a growing tensor perturbation for nearly all choices of n. This occurs even when the background model is undergoing accelerated expansion as opposed to the case of General Relativity. Consequently, cosmological gravitational wave modes can in principle provide a strong constraint on the theory of gravity independent of other cosmological data sets.

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 present a detailed investigation of the cosmological dynamics based on $\exp (-R/{\Lambda})$ gravity. We apply the dynamical system approach to both the vacuum and matter cases and obtain exact solutions and their stability in the finite and asymptotic regimes. The results show that cosmic histories exist which admit a double de-Sitter phase which could be useful for describing the early and the late-time accelerating universe.

The evolution equations for tensor perturbations in a generic scalar tensor theory of gravity are presented. Exact solution are given for a specific class of theories and Friedmann-Lema\^{i}tre-Robertson-Walker backgrounds. In these cases it is shown that, although the evolution of tensor models depends on the choice of parameters of the theory, no amplification is possible if the gravitational interaction is attractive.

The covariant gauge invariant perturbation theory of scalar cosmological perturbations is developed for a general Scalar-Tensor Friedmann-Lemaitre-Robertson-Walker cosmology in a vacuum. The perturbation equations are then solved exactly in the long wavelength limit for a specific coupling, potential and background. Differences with the minimally coupled case are briefly discussed.

We give the equations governing the shear evolution in Bianchi spacetimes for general f(R)-theories of gravity. We consider the case of R^n-gravity and perform a detailed analysis of the dynamics in Bianchi I cosmologies which exhibit local rotational symmetry. We find exact solutions and study their behaviour and stability in terms of the values of the parameter n. In particular, we found a set of cosmic histories in which the universe is initially isotropic, then develops shear anisotropies which approaches a constant value.

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 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 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.

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.