Inflation generically produces primordial gravitational waves with a red spectral tilt. In this paper we calculate the backreaction produced by these gravitational waves on the expansion of the universe. We find that in radiation domination the backreaction acts as a relativistic fluid, while in matter domination a small dark energy emerges with an equation of state w=-8/9.

We investigate the effect that the average backreaction of structure formation has on the dynamics of the cosmological expansion, within the concordance model. Our approach in the Poisson gauge is fully consistent up to second-order in a perturbative expansion about a flat Friedmann background, including a cosmological constant. We discuss the key length scales which are inherent in any averaging procedure of this kind. We identify an intrinsic homogeneity scale that arises from the averaging procedure, beyond which a residual offset remains in the expansion rate and deceleration parameter. In the case of the deceleration parameter, this can lead to a quite large increase in the value - more than 10% - and may therefore have important ramifications for dark energy measurements, even if the underlying nature of dark energy is a cosmological constant. We give the intrinsic variance that affects the value of the effective Hubble rate and deceleration parameter. These considerations serve to add extra intrinsic errors to our determination of the cosmological parameters, and, in particular, may render attempts to measure the Hubble constant to percent precision overly optimistic.

We investigate the spectrum of vector modes today which is generated at second order by density perturbations. The vector mode background that is generated by structure formation is small but in principle it contributes to the integrated Sachs-Wolfe effect, to redshift-space distortions and to weak lensing. We recover, clarify and extend previous results, and explain carefully why no vorticity is generated in the fluid at second order. The amplitude of the induced vector mode in the metric is around 1% that of the first-order scalars on small scales. We also calculate the power spectrum and the energy density of the vector part of the shear at second order.

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.

While vector modes are usually ignored in cosmology since they are not produced during inflation they are inevitably produced from the interaction of density fluctuations of differing wavelengths. This effect may be calculated via a second-order perturbative expansion. We investigate this effect during the radiation era. We discuss the generation mechanism by investigating two scalar modes interacting, and we calculate the power of vector modes generated by a power-law spectrum of density perturbations on all scales.

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 derive the Klein--Gordon equation for a single scalar field coupled to gravity at second order in perturbation theory and leading order in slow-roll. This is done in two ways: we derive the Klein--Gordon equation first using the Einstein field equations, and then directly from the action after integrating out the constraint equations. We also point out an unexpected result regarding the treatment of the field equations.

We discuss the gravitational wave background generated by primordial density perturbations evolving during the radiation era. At second-order in a perturbative expansion, density fluctuations produce gravitational waves. We calculate the power spectra of gravitational waves from this mechanism, and show that, in principle, future gravitational wave detectors could be used to constrain the primordial power spectrum on scales vastly different from those currently being probed by large-scale structure. As examples we compute the gravitational wave background generated by both a power-law spectrum on all scales, and a delta-function power spectrum on a single scale.

In general relativity, for fluids with a linear equation of state (EoS) or scalar fields, the high isotropy of the universe requires special initial conditions, and singularities are anisotropic in general. In the brane world scenario anisotropy at the singularity is suppressed by an effective quadratic equation of state. There is no reason why the effective EoS of matter should be linear at the highest energies, and a non-linear EoS may describe dark energy or unified dark matter (Paper I, astro-ph/0512224). In view of this, here we study the effects of a quadratic EoS in homogenous and inhomogeneous cosmological models in general relativity, in order to understand if in this context the quadratic EoS can isotropize the universe at early times. With respect to Paper I, here we use the simplified EoS P=alpha rho + rho^2/rho_c, which still allows for an effective cosmological constant and phantom behavior, and is general enough to analyze the dynamics at high energies. We first study anisotropic Bianchi I and V models, focusing on singularities. Using dynamical systems methods, we find the fixed points of the system and study their stability. We find that models with standard non-phantom behavior are in general asymptotic in the past to an isotropic fixed point IS, i.e. in these models even an arbitrarily large anisotropy is suppressed in the past: the singularity is matter dominated. Using covariant and gauge invariant variables, we then study linear perturbations about the homogenous and isotropic spatially flat models with a quadratic EoS. We find that, in the large scale limit, all perturbations decay asymptotically in the past, indicating that the isotropic fixed point IS is the general asymptotic past attractor for non phantom inhomogeneous models with a quadratic EoS. (Abridged)

We investigate the general relativistic dynamics of Robertson-Walker models with a non-linear equation of state (EoS), focusing on the quadratic case P = P_0 + \alpha \rho + \beta \rho^2. This may be taken to represent the Taylor expansion of any arbitrary barotropic EoS, P(\rho). With the right combination of P_0, \alpha and \beta, it serves as a simple phenomenological model for dark energy, or even unified dark matter. Indeed we show that this simple model for the EoS can produce a large variety of qualitatively different dynamical behaviors that we classify using dynamical systems theory. An almost universal feature is that accelerated expansion phases are mostly natural for these non-linear EoS's. These are often asymptotically de Sitter thanks to the appearance of an effective cosmological constant. Other interesting possibilities that arise from the quadratic EoS are closed models that can oscillate with no singularity, models that bounce between infinite contraction/expansion and models which evolve from a phantom phase, asymptotically approaching a de Sitter phase instead of evolving to a "Big Rip". In a second paper we investigate the effects of the quadratic EoS in inhomogeneous and anisotropic models, focusing in particular on singularities.