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 derive upper and lower limits for the mass-radius ratio of spin-fluid spheres in Einstein-Cartan theory, with matter satisfying a linear barotropic equation of state, and in the presence of a cosmological constant. Adopting a spherically symmetric interior geometry, we obtain the generalized continuity and Tolman-Oppenheimer-Volkoff equations for a Weyssenhoff spin-fluid in hydrostatic equilibrium, expressed in terms of the effective mass, density and pressure, all of which contain additional contributions from the spin. The generalized Buchdahl inequality, which remains valid at any point in the interior, is obtained, and general theoretical limits for the maximum and minimum mass-radius ratios are derived. As an application of our results we obtain gravitational red shift bounds for compact spin-fluid objects, which may (in principle) be used for observational tests of Einstein-Cartan theory in an astrophysical context. We also briefly consider applications of the torsion-induced minimum mass to the spin-generalized strong gravity model for baryons/mesons, and show that the existence of quantum spin imposes a lower bound for spinning particles, which almost exactly reproduces the electron mass.

New classes of modified teleparallel theories of gravity are introduced. The action of this theory is constructed to be a function of the irreducible parts of torsion $f(T_{\rm ax},T_{\rm ten},T_{\rm vec})$, where $T_{\rm ax},T_{\rm ten}$ and $T_{\rm vec}$ are squares of the axial, tensor and vector components of torsion, respectively. This is the most general (well-motivated) second order teleparallel theory of gravity that can be constructed from the torsion tensor. Different particular second order theories can be recovered from this theory such as new general relativity, conformal teleparallel gravity or $f(T)$ gravity. Additionally, the boundary term $B$ which connects the Ricci scalar with the torsion scalar via $R=-T+B$ can also be incorporated into the action. By performing a conformal transformation, it is shown that the two unique theories which have an Einstein frame are either the teleparallel equivalent of general relativity or $f(-T+B)=f(R)$ gravity, as expected.

We investigate modified theories of gravity in the context of teleparallel geometries. It is well known that modified gravity models based on the torsion scalar are not invariant under local Lorentz transformations while modifications based on the Ricci scalar are. This motivates the study of a model depending on the torsion scalar and the divergence of the torsion vector. We derive the teleparallel equivalent of $f(R)$ gravity as a particular subset of these models and also show that this is the unique theory in this class that is invariant under local Lorentz transformation. Furthermore one can show that $f(T)$ gravity is the unique theory admitting second order field equations.

A new Lagrangian framework has recently been proposed to describe interactions between relativistic perfect fluids and scalar fields. In this paper we investigate the Einstein static universe in this new class of theories, which have been named Scalar-Fluid theories. The stability of the static solutions to both homogeneous and inhomogeneous perturbations is analysed deriving the relevant cosmological perturbation equations at the linear order. We can find several configurations corresponding to an Einstein static universes which are stable against inhomogeneous perturbations, but unstable against homogeneous perturbations. This shows the possible applications of Scalar-Fluid theories to the inflationary emergent universe scenario.

We present a new approach to build models of quintessence interacting with dark or baryonic matter. We use a variational approach for relativistic fluids to realize an effective description of matter fields at the Lagrangian level. The coupling is introduced directly in the action by considering a single function mixing the dynamical degrees of freedom of the theory. The resulting gravitational field equations are derived by variations with respect to the independent variables. New interesting phenomenology can be obtained at both small scales, where new screening mechanisms for scalar fields can be realized, and large scales, where one finds an original and rich class of interacting quintessence models. The background cosmology of two of these models is studied in detail using dynamical system techniques. We find a variety of interesting results: for instance, these models contain dark energy dominated late time attractors and scaling solutions, both with early time matter dominated epochs and a possible inflationary origin. In general this new approach provides the starting point for future in depth studies on new interacting quintessence models.

We consider an original variational approach for building new models of quintessence interacting with dark or baryonic matter. The coupling is introduced at the Lagrangian level using a variational formulation for relativistic fluids, where the interacting term generally depends on both the dynamical degrees of freedom of the theory and their spacetime derivatives. After deriving the field equations from the action, we consider applications in the context of cosmology. Two simple models are studied using dynamical system techniques showing the interesting phenomenology arising in this framework. We find that these models contain dark energy dominated late time attractors with early time matter dominated epochs and also obtain a possible dynamical crossing of the phantom barrier. The formulation and results presented here complete and expand the analysis exposed in the first part of this work, where only algebraic couplings, without spacetime derivatives, were considered.

We propose a new point of view for interpreting Newton's and Einstein's theories of gravity. By taking inspiration from Continuum Mechanics and its treatment of anisotropies, we formulate new gravitational actions for modified theories of gravity. These models are simple and natural generalisations with many interesting properties. Above all, their precise form can, in principle, be determined experimentally.

We analyse dark energy models where self-interacting three-forms or phantom fields drive the accelerated expansion of the Universe. The dynamics of such models is often studied by rewriting the cosmological field equations in the form of a system of autonomous differential equations, or simply a dynamical system. Properties of these systems are usually studied via linear stability theory. In situations where this method fails, for instance due to the presence of zero eigenvalues in the Jacobian, centre manifold theory can be applied. We present a concise introduction and show explicitly how to use this theory in two concrete examples.

We examine the existence of relativistic stars in f(T) modified gravity and explicitly construct several classes of static perfect fluid solutions. We derive the conservation equation from the complete f(T) gravity field equations and present the differences with its teleparallel counterpart. Firstly, we choose the tetrad field in the diagonal gauge and study the resulting field equations. Some exact solutions are explicitly constructed and it is noted that these solutions have to give a constant torsion scalar. Next, we choose a non diagonal tetrad field which results in field equations similar to those of general relativity. For specific models we are able to construct exact solutions of these field equations. Among those new classes of solutions, we find negative pressure solutions, and an interesting class of polynomial solutions.

We investigate the evolution of sub-horizon cold dark matter perturbation in the dark energy dominated era of the Universe. By generalising the Meszaros equation to be valid for an arbitrary equation of state parameter, we derive the $w$-Meszaros equation. Its solutions determine the evolution of the cold dark matter perturbation by neglecting dark energy perturbations. Our analytical results provide a qualitative understanding of this evolution.

We introduce and carefully define an entire class of field theories based on non-standard spinors. Their dominant interaction is via the gravitational field which makes them naturally dark; we refer to them as Dark Spinors. We provide a critical analysis of previous proposals for dark spinors noting that they violate Lorentz invariance. As a working assumption we restrict our analysis to non-standard spinors which preserve Lorentz invariance, whilst being non-local and explicitly construct such a theory. We construct the complete energy-momentum tensor and derive its components explicitly by assuming a specific projection operator. It is natural to next consider dark spinors in a cosmological setting. We find various interesting solutions where the spinor field leads to slow roll and fast roll de Sitter solutions. We also analyse models where the spinor is coupled conformally to gravity, and consider the perturbations and stability of the spinor.

The apparent alignment of the cosmic microwave background multipoles on large scales challenges the standard cosmological model. Scalar field inflation is isotropic and cannot account for the observed alignment. We explore the imprints, a non-standard spinor driven inflation would leave on the cosmic microwave background anisotropies. We show it is natural to expect an anisotropic inflationary expansion of the Universe which has the effect of suppressing the low multipole amplitude of the primordial power spectrum, while at the same time to provide the usual inflationary features.

We introduce a new form of coupling between dark energy and dark matter that is quadratic in their energy densities. Then we investigate the background dynamics when dark energy is in the form of exponential quintessence. The three types of quadratic coupling all admit late-time accelerating critical points, but these are not scaling solutions. We also show that two types of coupling allow for a suitable matter era at early times and acceleration at late times, while the third type of coupling does not admit a suitable matter era.

Ever since the first observations that we are living in an accelerating universe, it has been asked what dark energy is. There are various explanations all of which with have various draw backs or inconsistencies. Here we show that using a dark spinor field it is possible to have an equation of state that crosses the phantom divide, becoming a dark phantom spinor which evolves into dark energy. This type of equation of state has been mildly favored by experimental data, however, in the past there were hardly any candidate theories that satisfied this crossing without creating ghosts or causing a singularity which results in the universe essentially ripping. The dark spinor model converges to dark energy in a reasonable time frame avoiding the big rip and without attaining negative kinetic energy as it crosses the phantom divide.

We analyze the stability of the Einstein static universe by considering homogeneous perturbations in the context of f(G) modified Gauss-Bonnet theories of gravity. By considering a generic form of f(G), the stability region of the Einstein static universe is parameterized by the linear equation of state parameter w=p/rho and the second derivative f"(G) of the Gauss-Bonnet term.

We study Einstein static universes in the context of generic f(R) models. It is shown that Einstein static solutions exist for a wide variety of modified gravity models sourced by a barotropic perfect fluid with equation of state w=p/rho, but these solutions are always unstable to either homogeneous or inhomogeneous perturbations. Our general results are in agreement with specific models investigated in that past. We also discuss how our techniques can be applied to other scenarios in f(R) gravity.

We solve one of the open problems in Einstein-Cartan theory, namely we find a natural matter source whose spin angular momentum tensor is compatible with the cosmological principle. We analyze the resulting evolution equations and find that an epoch of accelerated expansion is an attractor. The torsion field quickly decays in that period. Our results are interpreted in the context of the standard model of cosmology.

Inflation driven by a single dark spinor field is discussed. We define the notion of a dark spinor field and derive the cosmological field equations for such a matter source. The conditions for inflation are determined and an exactly solvable model is presented. We find the power spectrum of the quantum fluctuation of this field and compare the results with scalar field inflation.

We consider the behavior of the tangential velocity of test particles moving in stable circular orbits in f(R) modified theories of gravity. A large number of observations at the galactic scale have shown that the rotational velocities of massive test particles (hydrogen clouds) tend towards constant values at large distances from the galactic center. We analyze the vacuum gravitational field equations in f(R) models in the constant velocity region, and the general form of the metric tensor is derived in a closed form. The resulting modification of the Einstein-Hilbert Lagrangian is of the form R^{1+n}, with the parameter n expressed in terms of the tangential velocity. Therefore we find that to explain the motion of test particles around galaxies requires only very mild deviations from classical general relativity, and that modified gravity can explain the galactic dynamics without the need of introducing dark matter.

We generalize the virial theorem in f(R) modified gravity using the collisionless Boltzmann equation. We find supplementary geometric terms in the modified Einstein equation providing an effective contribution to the gravitational energy. The total virial mass is proportional to the effective mass associated with the new geometrical term, which may account for the well-known virial theorem mass discrepancy in clusters of galaxies. The model predicts that the geometric mass and its effects extend beyond the virial radius of the clusters. We also consider the behavior of the galaxy cluster velocity dispersion in f(R) models. Thus, the f(R) virial theorem can be an efficient tool in observationally testing the viability of this class of generalized gravity models.

The classical tests of general relativity (perihelion precession, deflection of light, and the radar echo delay) are considered for the Dadhich, Maartens, Papadopoulos and Rezania (DMPR) solution of the spherically symmetric static vacuum field equations in brane world models. For this solution the metric in the vacuum exterior to a brane world star is similar to the Reissner-Nordstrom form of classical general relativity, with the role of the charge played by the tidal effects arising from projections of the fifth dimension. The existing observational solar system data on the perihelion shift of Mercury, on the light bending around the Sun (obtained using long-baseline radio interferometry), and ranging to Mars using the Viking lander, constrain the numerical values of the bulk tidal parameter and of the brane tension.

Dark energy and dark matter are the dominant sources in the evolution of the late universe. They are currently only indirectly detected via their gravitational effects, and there could be a coupling between them without violating observational constraints. We investigate the background dynamics when dark energy is modelled as exponential quintessence, and is coupled to dark matter via simple models of energy exchange. We introduce a new form of dark sector coupling, which leads to a more complicated dynamical phase space and has a better physical motivation than previous mathematically similar couplings.

We consider the astrophysical and cosmological implications of the existence of a minimum density and mass due to the presence of the cosmological constant. If there is a minimum length in nature, then there is an absolute minimum mass corresponding to a hypothetical particle with radius of the order of the Planck length. On the other hand, quantum mechanical considerations suggest a different minimum mass. These particles associated with the dark energy can be interpreted as the ``quanta'' of the cosmological constant. We study the possibility that these particles can form stable stellar-type configurations through gravitational condensation, and their Jeans and Chandrasekhar masses are estimated. From the requirement of the energetic stability of the minimum density configuration on a macroscopic scale one obtains a mass of the order of 10^55 g, of the same order of magnitude as the mass of the universe. This mass can also be interpreted as the Jeans mass of the dark energy fluid. Furthermore we present a representation of the cosmological constant and of the total mass of the universe in terms of `classical' fundamental constants.

We examine the Schwarzschild interior of a black hole, incorporating quantum gravitational modifications due to loop quantum gravity. We consider an improved loop quantization using techniques that have proven successful in loop quantum cosmology. The central Schwarzschild singularity is resolved and the implications for the fate of an in-falling test particle in the interior region is discussed. The singularity is replaced by a Nariai type Universe. We discuss the resulting conformal diagram, providing a clear geometrical interpretation of the quantum effects.

The cosmological constant sets certain scales important in cosmology. We show that Lambda in conjunction with other parameters like the Schwarzschild radius leads to scales relevant not only for cosmological but also for astrophysical applications. Of special interest is the extension of orbits and velocity of test particles traveling over Mpc distances. We will show that there exists a lower and an upper cut-off on the possible velocities of test particles. For a test body moving in a central gravitational field Lambda enforces a maximal value of the angular momentum if we insist on bound orbits of the test body which move at a distance larger than the Schwarzschild radius.

Astrophysical bounds on the cosmological constant are examined for spherically symmetric bodies. Similar limits emerge from hydrostatical and gravitational equilibrium and the validity of the Newtonian limit. It is argued that the bound coming from the analysis of exact solutions is not a feature of the relativistic theory since it also emerges from the Newtonian theory. The methods in use seem to be disjoint from the basic principles, however they have the same implication regarding the upper bounds. Therefore we will compare different inequalities and comment on the possible relationship between them.

Exact solutions of traversable wormholes are found under the assumption of spherical symmetry and the existence of a {\it non-static} conformal symmetry, which presents a more systematic approach in searching for exact wormhole solutions. In this work, a wide variety of solutions are deduced by considering choices for the form function, a specific linear equation of state relating the energy density and the pressure anisotropy, and various phantom wormhole geometries are explored. A large class of solutions impose that the spatial distribution of the exotic matter is restricted to the throat neighborhood, with a cut-off of the stress-energy tensor at a finite junction interface, although asymptotically flat exact solutions are also found. Using the ``volume integral quantifier,'' it is found that the conformally symmetric phantom wormhole geometries may, in principle, be constructed by infinitesimally small amounts of averaged null energy condition violating matter. Considering the tidal acceleration traversability conditions for the phantom wormhole geometry, specific wormhole dimensions and the traversal velocity are also deduced.

We consider the possibility that the dark matter, which is required to explain the dynamics of the neutral hydrogen clouds at large distances from the galactic center, could be in the form of a Bose-Einstein condensate. To study the condensate we use the non-relativistic Gross-Pitaevskii equation. By introducing the Madelung representation of the wave function, we formulate the dynamics of the system in terms of the continuity equation and of the hydrodynamic Euler equations. Hence dark matter can be described as a non-relativistic, Newtonian Bose-Einstein gravitational condensate gas, whose density and pressure are related by a barotropic equation of state. In the case of a condensate with quartic non-linearity, the equation of state is polytropic with index $n=1$. To test the validity of the model we fit the Newtonian tangential velocity equation of the model with a sample of rotation curves of low surface brightness and dwarf galaxies, respectively. We find a very good agreement between the theoretical rotation curves and the observational data for the low surface brightness galaxies. The deflection of photons passing through the dark matter halos is also analyzed, and the bending angle of light is computed. The bending angle obtained for the Bose-Einstein condensate is larger than that predicted by standard general relativistic and dark matter models. Therefore the study of the light deflection by galaxies and the gravitational lensing could discriminate between the Bose-Einstein condensate dark matter model and other dark matter models.

We analyze the stability of the Einstein static universe by considering homogeneous scalar perturbations in the context of f(R) modified theories of gravity. By considering specific forms of f(R), the stability regions of the solutions are parameterized by a linear equation of state parameter w=p/rho. Contrary to classical general relativity, it is found that in f(R) gravity a stable Einstein cosmos with a positive cosmological constant does indeed exist. Thus, we are lead to conclude that, in principle, modifications in f(R) gravity stabilize solutions which are unstable in general relativity.

The behavior of the angular velocity of a test particle moving in a stable circular orbit in the vacuum on the brane is considered. In the brane world scenario, the four dimensional effective Einstein equation acquire extra terms, called dark radiation and dark pressure, respectively, which arise from the embedding of the 3-brane in the bulk. A large number of independent observations have shown that the rotational velocities of test particles gravitating around galaxies tend, as a function of the distance from the galactic center, toward constant values. By assuming a constant tangential velocity, the general solution of the vacuum gravitational field equations on the brane can be obtained in an exact analytic form. This allows us to obtain the explicit form of the projections of the bulk Weyl tensor on the brane, and the equation of state of the dark pressure as a function of the dark radiation. The physical and geometrical quantities are expressed in terms of observable/measurable parameters, like the tangential velocity, the baryonic mass and the radius of the galaxy. We also analyze the dynamics of test particles by using methods from the qualitative analysis of dynamical systems, by assuming a simple linear equation of state for the dark pressure. The obtained results provide a theoretical framework for the observational testing at the extra-galactic scale of the predictions of the brane world models.

The equation of motion for test particles in $f(R)$ modified theories of gravity is derived. By considering an explicit coupling between an arbitrary function of the scalar curvature, $R$, and the Lagrangian density of matter, it is shown that an extra force arises. This extra force is orthogonal to the four-velocity and the corresponding acceleration law is obtained in the weak field limit. Connections with MOND and with the Pioneer anomaly are further discussed.

We consider global and gravitational lensing properties of the recently suggested Einstein clusters of WIMPs as galactic dark matter halos. Being tangential pressure dominated, Einstein clusters are strongly anisotropic systems which can describe any galactic rotation curve by specifying the anisotropy. Due to this property, Einstein clusters may be considered as dark matter candidates. We analyse the stability of the Einstein clusters against both radial and non-radial pulsations, and we show that the Einstein clusters are dynamically stable. With the use of the Buchdahl type inequalities for anisotropic bodies, we derive upper limits on the velocity of the particles defining the cluster. These limits are consistent with those obtained from stability considerations. The study of light deflection shows that the gravitational lensing effect is slightly smaller for the Einstein clusters, as compared to the singular isothermal density sphere model for dark matter. Therefore lensing observations may discriminate, at least in principle, between Einstein cluster and other dark matter models.

The equations describing the adiabatic, small radial oscillations of general relativistic stars are generalized to include the effects of a cosmological constant. The generalized eigenvalue equation for the normal modes is used to study the changes in the stability of the homogeneous sphere induced by the presence of the cosmological constant. The variation of the critical adiabatic index as a function of the central pressure is studied numerically for different trial functions. The presence of a large cosmological constant significantly increases the value of the critical adiabatic index. The dynamical stability condition of the homogeneous star in the Schwarzschild-de Sitter geometry is obtained and several bounds on the maximum allowable value for a cosmological constant are derived from stability considerations.

Recently, a spin one half matter field with mass dimension one was discovered, called Elko spinors. The present work shows how to introduce these fields into a curved spacetime by the standard covariantisation scheme. After formulating the coupled Einstein-Elko field equations, the spacetime is assumed to be homogeneous and isotropic in order to simplify the resulting field equations. Analytical ghost Elko solutions are constructed which have vanishing energy-momentum tensor without and with cosmological constant. The cosmological Elko theory is finally related to the standard scalar field theory with self interaction that gives rise to inflation and it is pointed out that the Elko spinors are not only prime dark matter candidates but also prime candidates for inflation.

We propose that the Universe is filled with a massive vector field, non-minimally coupled to gravitation. The field equations of the model are consistently derived and their application to cosmology is considered. The Friedmann equations acquire an extra dark-energy component, which is proportional to the mass of the vector particle. This leads to a late-time accelerated de Sitter type expansion. The free parameters of the model (gravitational coupling constants and initial value of the cosmological vector field) can be estimated by using the PPN solar system constraints. The mass of the cosmological massive vector particle, which may represent the main component of the Universe, is of the order of 10^-63 g.

In the field equations of Einstein-Cartan theory with cosmological constant a static spherically symmetric perfect fluid with spin density satisfying the Weyssenhoff restriction is considered. This serves as a rough model of space filled with (fermionic) dark matter. From this the Einstein static universe with constant torsion is constructed, generalising the Einstein Cosmos to Einstein-Cartan theory. The interplay between torsion and the cosmological constant is discussed. A possible way out of the cosmological constant's sign problem is suggested.

This paper has been withdrawn by the authors because the ratio presented in Eq. (13) does not remain frozen in time for all values of w = p/rho. As such a separate more detailed analysis containing the generalized Buchdahl inequality valid for all values of w is required. The new analysis will in addition attend to the time evolution of the ratio (13). The essential results remain unchanged, but result from the w=0 case.

We show that in the framework of the classical general relativity the presence of a positive cosmological constant implies the existence of a minimal mass and of a minimal density in nature. These results rigorously follow from the generalized Buchdahl inequality in the presence of a cosmological constant.