We propose a relativistic gravitational theory leading to Modified Newtonian Dynamics, a paradigm that explains the observed universal acceleration and associated phenomenology in galaxies. We discuss phenomenological requirements leading to its construction and demonstrate its agreement with the observed Cosmic Microwave Background and matter power spectra on linear cosmological scales. We show that its action expanded to 2nd order is free of ghost instabilities and discuss its possible embedding in a more fundamental theory.

We perform the first test of Dark Matter (DM) stress-energy evolution through cosmic history, using Cosmic Microwave Background measurements supplemented with Baryon Acoustic Oscillation data and the Hubble Space Telescope key project data. We constrain the DM equation of state (EoS) in 8 redshift bins, and its sound speed and (shear) viscosity in 9 redshift bins, finding no convincing evidence for non-$\Lambda$CDM values in any of the redshift bins. Despite this enlarged parameter space, the sound speed and viscosity are constrained relatively well at late times (due to the inclusion of CMB lensing), whereas the EoS is most strongly constrained around recombination. These results constrain for the first time the level of ''coldness'' required of DM across various cosmological epochs at both the background and perturbative levels. We show that simultaneously allowing time dependence for both the EoS and sound speed parameters shifts the posterior of the DM abundance before recombination to a higher value, whilst keeping the present day DM abundance similar to the $\Lambda$CDM value. This shifts the posterior for the present day Hubble constant compared to $\Lambda$CDM, suggesting that DM with time-dependent parameters is well-suited to explore possible solutions to persistent tensions within the $\Lambda$CDM model.

A number of theories of gravity have been proposed as proxies for dark matter in the regime of galaxies and cosmology. The recent observations of gravitational waves (GW170817) from the merger of two neutron stars, followed by an electromagnetic counterpart (GW170817a) have placed stringent constraints on the difference of the speed of gravity to the speed of light, severely restricting the phenomenological viability of such theories. We revisit the impact of these observations on the Tensor-Vector-Scalar (TeVeS) paradigm of relativistic Modified Newtonian Dynamics (MOND) and demonstrate the existence of a previously unknown class of this paradigm where the speed of gravity always equals the speed of light. We show that this holds without altering the usual (bimetric) MOND phenomenology in galaxies.

Cold Dark Matter (CDM) is a crucial constituent of the current concordance cosmological model. Having a vanishing equation of state (EoS), its energy density scales with the inverse cosmic volume and is thus uniquely described by a single number, its present abundance. We test the inverse cosmic volume law for Dark Matter (DM) by allowing its EoS to vary independently in eight redshift bins in the range $z=10^5$ and $z=0$. We use the latest measurements of the Cosmic Microwave Background radiation from the Planck satellite and supplement them with Baryon Acoustic Oscillation (BAO) data from the 6dF and SDSS-III BOSS surveys, and with the Hubble Space Telescope (HST) key project data. We find no evidence for nonzero EoS in any of the eight redshift bins. With Planck data alone, the DM abundance is most strongly constrained around matter-radiation equality $\omega^{\rm eq}_g = 0.1193^{+0.0036}_{-0.0035}$ (95% c.l.), whereas its present day value is more weakly constrained $\omega^{(0)}_g = 0.16^{+0.12}_{-0.10}$ (95% c.l.). Adding BAO or HST data does not significantly change the $\omega^{\rm eq}_g$ constraint, while $\omega^{(0)}_g$ tightens to $0.160^{+0.069}_{-0.065} $ (95% c.l.) and $0.124^{+0.081}_{-0.067}$ (95% c.l.) respectively. Our results constrain for the first time the level of "coldness" required of the DM across various cosmological epochs and show that the DM abundance is strictly positive at all times.

We demonstrate that the Vlasov equation describing collisionless self-gravitating matter may be solved with the so-called Schr\"odinger method (ScM). With the ScM, one solves the Schr\"odinger-Poisson system of equations for a complex wave function in d dimensions, rather than the Vlasov equation for a 2d-dimensional phase space density. The ScM also allows calculating the d-dimensional cumulants directly through quasi-local manipulations of the wave function, avoiding the complexity of 2d-dimensional phase space. We perform for the first time a quantitive comparison of the ScM and a conventional Vlasov solver in d=2 dimensions. Our numerical tests were carried out using two types of cold cosmological initial conditions: the classic collapse of a sine wave and those of a gaussian random field as commonly used in cosmological cold dark matter N-body simulations. We compare the first three cumulants, that is, the density, velocity and velocity dispersion, to those obtained by solving the Vlasov equation using the publicly available code ColDICE. We find excellent qualitative and quantitative agreement between these codes, demonstrating the feasibility and advantages of the ScM as an alternative to N-body simulations. We discuss, the emergence of effective vorticity in the ScM through the winding number around the points where the wave function vanishes. As an application we evaluate the background pressure induced by the non-linearity of large scale structure formation, thereby estimating the magnitude of cosmological backreaction. We find that it is negligibly small and has time dependence and magnitude compatible with expectations from the effective field theory of large scale structure.

We compare Einstein-Boltzmann solvers that include modifications to General Relativity and find that, for a wide range of models and parameters, they agree to a high level of precision. We look at three general purpose codes that primarily model general scalar-tensor theories, three codes that model Jordan-Brans-Dicke (JBD) gravity, a code that models f(R) gravity, a code that models covariant Galileons, a code that models Ho\v{r}ava-Lifschitz gravity and two codes that model non-local models of gravity. Comparing predictions of the angular power spectrum of the cosmic microwave background and the power spectrum of dark matter for a suite of different models, we find agreement at the sub-percent level. This means that this suite of Einstein-Boltzmann solvers is now sufficiently accurate for precision constraints on cosmological and gravitational parameters.

The maximum size of a cosmic structure is given by the maximum turnaround radius -- the scale where the attraction due to its mass is balanced by the repulsion due to dark energy. We derive generic formulae for the estimation of the maximum turnaround radius in any theory of gravity obeying the Einstein equivalence principle, in two situations: on a spherically symmetric spacetime and on a perturbed Friedman-Robertson-Walker spacetime. We show that the two formulae agree. As an application of our formula, we calculate the maximum turnaround radius in the case of the Brans-Dicke theory of gravity. We find that for this theory, such maximum sizes always lie above the \LCDM value, by a factor $1 + \frac{1}{3\omega}$, where $\omega\gg 1$ is the Brans-Dicke parameter, implying consistency of the theory with current data.

A generalization of Bekenstein's Tensor-Vector-Scalar (TeVeS) model of modified gravity has recently been proposed as an alternative to dark matter. This model -- which we will refer to as g-TeVeS -- utilizes a Galileon-induced Vainshtein mechanism to suppress modifications to General Relativity in strong gravity regimes and so avoids the need to introduce the baroque kinetic terms that typically exist in relativistic models of Modified Newtonian Dynamics (MOND). We explore the behavior of this model in spacetimes with exact Friedmann-Robertson-Walker (FRW) symmetry. The ability of the theory to recover MOND phenomenology places restrictions on the theory's parameter space and it is found that within an estimate of this area of parameter space the theory cannot successfully approximate the Friedmannian cosmological behavior of cold dark matter. It is found that much closer agreement may be recovered in other regions of the theory's parameter space and the reasons for this are discussed.

We examine how the properties of dark matter, parameterised by an equation of state parameter $w$ and two perturbative Generalised Dark Matter (GDM) parameters $c^2_s$ (the sound speed) and $c^2_\text{vis}$ (the viscosity), are constrained by existing cosmological data, particularly the Planck 2015 data release. We find that the GDM parameters are consistent with zero, and are strongly constrained, showing no evidence for extending the dark matter model beyond the Cold Dark Matter (CDM) paradigm. The dark matter equation of state is constrained to be within $-0.000896<w<0.00238$ at the $99.7\%$ confidence level, which is several times stronger than constraints found previously using WMAP data. The parameters $c^2_s$ and $c^2_\text{vis}$ are constrained to be less than $3.21\times10^{-6}$ and $6.06\times10^{-6}$ respectively at the $99.7\%$ confidence level. The inclusion of the GDM parameters does significantly affect the error bars on several $\Lambda$CDM parameters, notably the dimensionless dark matter density $\omega_g$ and the derived parameters $\sigma_8$ and $H_0$. This can be partially alleviated with the inclusion of data constraining the expansion history of the universe.

The Cold Dark Matter (CDM) model, wherein the dark matter is treated as a pressureless perfect fluid, provides a good fit to galactic and cosmological data. With the advent of precision cosmology, it should be asked whether this simplest model needs to be extended, and whether doing so could improve our understanding of the properties of dark matter. One established parameterisation for generalising the CDM fluid is the Generalised Dark Matter (GDM) model, in which dark matter is an imperfect fluid with pressure and shear viscosity that fulfill certain closure equations. We investigate these closure equations and the three new parametric functions they contain: the background equation of state w, the speed of sound c_s^2 and the viscosity c_{vis}^2. Taking these functions to be constant parameters, we analyse an exact solution of the perturbed Einstein equations in a GDM-dominated universe and discuss the main effects of the three parameters on the Cosmic Microwave Background (CMB). Our analysis suggests that the CMB alone is not able to distinguish between the GDM sound speed and viscosity parameters, but that other observables, such as the matter power spectrum, are required to break this degeneracy. We also consider other descriptions of imperfect fluids that have a non-perturbative definition and relate these to the GDM model. In particular, we consider scalar fields, an effective field theory (EFT) of fluids, an EFT of Large Scale Structure, non-equilibrium thermodynamics and tightly-coupled fluids. These descriptions could be used to extend the GDM model into the nonlinear regime of structure formation, which is necessary if the wealth of data available on those scales is to be employed in constraining the model. We also derive the initial conditions for adiabatic and isocurvature perturbations and provide the result in a form ready for implementation in Einstein-Boltzmann solvers.

Euclid is a European Space Agency medium class mission selected for launch in 2020 within the Cosmic Vision 2015 2025 program. The main goal of Euclid is to understand the origin of the accelerated expansion of the universe. Euclid will explore the expansion history of the universe and the evolution of cosmic structures by measuring shapes and redshifts of galaxies as well as the distribution of clusters of galaxies over a large fraction of the sky. Although the main driver for Euclid is the nature of dark energy, Euclid science covers a vast range of topics, from cosmology to galaxy evolution to planetary research. In this review we focus on cosmology and fundamental physics, with a strong emphasis on science beyond the current standard models. We discuss five broad topics: dark energy and modified gravity, dark matter, initial conditions, basic assumptions and questions of methodology in the data analysis. This review has been planned and carried out within Euclid's Theory Working Group and is meant to provide a guide to the scientific themes that will underlie the activity of the group during the preparation of the Euclid mission.

Euclid is a European Space Agency medium class mission selected for launch in 2019 within the Cosmic Vision 2015-2025 programme. The main goal of Euclid is to understand the origin of the accelerated expansion of the Universe. Euclid will explore the expansion history of the Universe and the evolution of cosmic structures by measuring shapes and redshifts of galaxies as well as the distribution of clusters of galaxies over a large fraction of the sky. Although the main driver for Euclid is the nature of dark energy, Euclid science covers a vast range of topics, from cosmology to galaxy evolution to planetary research. In this review we focus on cosmology and fundamental physics, with a strong emphasis on science beyond the current standard models. We discuss five broad topics: dark energy and modified gravity, dark matter, initial conditions, basic assumptions and questions of methodology in the data analysis. This review has been planned and carried out within Euclid's Theory Working Group and is meant to provide a guide to the scientific themes that will underlie the activity of the group during the preparation of the Euclid mission.

Light degrees of freedom that modify gravity on cosmological scales must be "screened" on solar system scales in order to be compatible with data. The Vainshtein mechanism achieves this through a breakdown of classical perturbation theory, as large interactions involving new degrees of freedom become important below the so-called Vainshtein radius. We begin to develop an extension of the Parameterized Post-Newtonian (PPN) formalism that is able to handle Vainshteinian corrections. We argue that theories with a unique Vainshtein scale must be expanded using two small parameters. In this Parameterized Post-Newtonian-Vainshteinian (PPNV) expansion, the primary expansion parameter that controls the PPN order is, as usual, the velocity $v$. The secondary expansion parameter, $\alpha$, controls the strength of the Vainshteinian correction and is a theory-specific combination of the Schwarzschild radius and the Vainshtein radius of the source that is independent of its mass. We present the general framework and apply it to Cubic Galileon theory both inside and outside the Vainshtein radius. The PPNV framework can be used to determine the compatibility of such theories with solar system and other strong-field data.

The current effort to test General Relativity employs multiple disparate formalisms for different observables, obscuring the relations between laboratory, astrophysical and cosmological constraints. To remedy this situation, we develop a parameter space for comparing tests of gravity on all scales in the universe. In particular, we present new methods for linking cosmological large-scale structure, the Cosmic Microwave Background and gravitational waves with classic PPN tests of gravity. Diagrams of this gravitational parameter space reveal a noticeable untested regime. The untested window, which separates small-scale systems from the troubled cosmological regime, could potentially hide the onset of corrections to General Relativity.

We present the most general parametrisation of models of dark energy in the form of a scalar field which is explicitly coupled to dark matter. We follow and extend the Parameterized Post-Friedmannian approach, previously applied to modified gravity theories, in order to include interacting dark energy. We demonstrate its use through a number of worked examples and show how the initially large parameter space of free functions can be significantly reduced and constrained to include only a few non-zero coefficients. This paves the way for a model-independent approach to classify and test interacting dark energy theories.

The growth rate of the large-scale structure of the universe has been advocated as the observable par excellence for testing gravity on cosmological scales. By considering linear-order deviations from General Relativity, we show that corrections to the growth rate, f, can be expressed as an integral over a `source' term, weighted by a theory-independent `response kernel'. This leads to an efficient and accurate `plug-and-play' expression for generating growth rates in alternative gravity theories, bypassing lengthy theory-specific computations. We use this approach to explicitly show that f is sensitive to a degenerate combination of modified expansion and modified clustering effects. Hence the growth rate, when used in isolation, is not a straightforward diagnostic of modified gravity.

We report strong cosmological constraints on the Brans-Dicke (BD) theory of gravity using Cosmic Microwave Background data from Planck.We consider two types of models. First, the initial condition of the scalar field is fixed to give the same effective gravitational strength $G_{eff}$ today as the one measured on the Earth, $G_N$. In this case the BD parameter $\omega$ is constrained to $\omega > 692$ at the $99\%$ confidence level, an order of magnitude improvement over previous constraints.In the second type the initial condition for the scalar is a free parameter leading to a somewhat stronger constraint of $\omega > 890$ while $G_{eff}$ is constrained to $0.981 <\frac{G_{eff}}{G_N} <1.285$ at the same confidence level. We argue that these constraints have greater validity than for the BD theory and are valid for any Horndeski theory, the most general second-order scalar-tensor theory, which approximates BD on cosmological scales. In this sense, our constraints place strong limits on possible modifications of gravity that might explain cosmic acceleration.

We propose a completely new parametrisation of the dark energy equation of state, which uses the dark energy density, $\Omega_e$ as a cosmic clock. We expand the equation of state in a series of orthogonal polynomials, with $\Omega_e$ as the expansion parameter and determine the expansion coefficients by fitting to SNIa and $H(z)$ data. Assuming that $\Omega_e$ is a monotonic function of time, we show that our parametrisation performs better than the popular Chevallier--Polarski--Linder (CPL) and Gerke and Efstathiou (GE) parametrisations, and we demonstrate that it is robust to the choice of prior. Expanding in orthogonal polynomials allows us to relate models of dark energy directly to our parametrisation, which we illustrate by placing constraints on the expansion coefficients extracted from two popular quintessence models. Finally, we comment on how this parametrisation could be modified to accommodate high redshift data, where any non--monotonicity of $\Omega_e$ would need to be accounted for.

We present three distinct types of models of dark energy in the form of a scalar field which is explicitly coupled to dark matter. Our construction draws from the pull-back formalism for fluids and generalises the fluid action to involve couplings to the scalar field. We investigate the cosmology of each class of model both at the background and linearly perturbed level. We choose a potential for the scalar field and a specific coupling function for each class of models and we compute the Cosmic Microwave Background and matter power spectra.

A unified framework for theories of modified gravity will be an essential tool for interpreting the forthcoming deluge of cosmological data. We present such a formalism, the Parameterized Post-Friedmann framework (PPF), which parameterizes the cosmological perturbation theory of a wide variety of modified gravity models. PPF is able to handle spin-0 degrees of freedom from new scalar, vector and tensor fields, meaning that it is not restricted to simple models based solely on cosmological scalar fields. A direct correspondence is maintained between the parameterization and the underlying space of theories, which allows us to build up a `dictionary' of modified gravity theories and their PPF correspondences. In this paper we describe the construction of the parameterization and demonstrate its use through a number of worked examples relevant to the current literature. We indicate how the formalism will be implemented numerically, so that the dictionary of modified gravity can be pitted against forthcoming observations.

There are a number of approaches to testing General Relativity (GR) on linear scales using parameterized frameworks for modifying cosmological perturbation theory. It is sometimes assumed that the details of any given parameterization are unimportant if one uses it as a diagnostic for deviations from GR. In this brief report we argue that this is not necessarily so. First we show that adopting alternative combinations of modifications to the field equations significantly changes the constraints that one obtains. In addition, we show that using a parameterization with insufficient freedom significantly tightens the apparent theoretical constraints. Fundamentally we argue that it is almost never appropriate to consider modifications to the perturbed Einstein equations as being constraints on the effective gravitational constant, for example, in the same sense that solar system constraints are. The only consistent modifications are either those that grant near-total freedom, as in decomposition methods, or ones which map directly to a particular part of theory space.

In this review we present a thoroughly comprehensive survey of recent work on modified theories of gravity and their cosmological consequences. Amongst other things, we cover General Relativity, Scalar-Tensor, Einstein-Aether, and Bimetric theories, as well as TeVeS, f(R), general higher-order theories, Horava-Lifschitz gravity, Galileons, Ghost Condensates, and models of extra dimensions including Kaluza-Klein, Randall-Sundrum, DGP, and higher co-dimension braneworlds. We also review attempts to construct a Parameterised Post-Friedmannian formalism, that can be used to constrain deviations from General Relativity in cosmology, and that is suitable for comparison with data on the largest scales. These subjects have been intensively studied over the past decade, largely motivated by rapid progress in the field of observational cosmology that now allows, for the first time, precision tests of fundamental physics on the scale of the observable Universe. The purpose of this review is to provide a reference tool for researchers and students in cosmology and gravitational physics, as well as a self-contained, comprehensive and up-to-date introduction to the subject as a whole.

There is a distinct possibility that current and future cosmological data can be used to constrain Einstein's theory of gravity on the very largest scales. To be able to do this in a model-independent way, it makes sense to work with a general parameterization of modified gravity. Such an approach would be analogous to the Parameterized Post-Newtonian (PPN) approach which is used on the scale of the Solar System. A few such parameterizations have been proposed and preliminary constraints have been obtained. We show that the majority of such parameterizations are only exactly applicable in the quasistatic regime. On larger scales they fail to encapsulate the full behaviour of typical models currently under consideration. We suggest that it may be possible to capture the additions to the `Parameterized Post-Friedmann' (PPF) formalism by treating them akin to fluid perturbations.

Euclid is a space-based survey mission from the European Space Agency designed to understand the origin of the Universe's accelerating expansion. It will use cosmological probes to investigate the nature of dark energy, dark matter and gravity by tracking their observational signatures on the geometry of the universe and on the cosmic history of structure formation. The mission is optimised for two independent primary cosmological probes: Weak gravitational Lensing (WL) and Baryonic Acoustic Oscillations (BAO). The Euclid payload consists of a 1.2 m Korsch telescope designed to provide a large field of view. It carries two instruments with a common field-of-view of ~0.54 deg2: the visual imager (VIS) and the near infrared instrument (NISP) which contains a slitless spectrometer and a three bands photometer. The Euclid wide survey will cover 15,000 deg2 of the extragalactic sky and is complemented by two 20 deg2 deep fields. For WL, Euclid measures the shapes of 30-40 resolved galaxies per arcmin2 in one broad visible R+I+Z band (550-920 nm). The photometric redshifts for these galaxies reach a precision of dz/(1+z) < 0.05. They are derived from three additional Euclid NIR bands (Y, J, H in the range 0.92-2.0 micron), complemented by ground based photometry in visible bands derived from public data or through engaged collaborations. The BAO are determined from a spectroscopic survey with a redshift accuracy dz/(1+z) =0.001. The slitless spectrometer, with spectral resolution ~250, predominantly detects Ha emission line galaxies. Euclid is a Medium Class mission of the ESA Cosmic Vision 2015-2025 programme, with a foreseen launch date in 2019. This report (also known as the Euclid Red Book) describes the outcome of the Phase A study.

We consider the post-Newtonian limit of a general class of bimetric theories of gravity, in which both metrics are dynamical. The established parameterised post-Newtonian approach is followed as closely as possible, although new potentials are found that do not exist within the standard frame-work. It is found that these theories can evade solar system tests of post-Newtonian gravity remarkably well. We show that perturbations about Minkowski space in these theories contain both massless and massive degrees of freedom, and that in general there are two different types of massive mode, each with a different mass parameter. If both of these masses are sufficiently large then the predictions of the most general class of theories we consider are indistinguishable from those of general relativity, up to post-Newtonian order in a weak field, low velocity expansion. In the limit that the massive modes become massless, we find that these general theories do not exhibit a van Dam-Veltman-Zakharov-like discontinuity in their gamma parameter, although there are discontinuities in other post-Newtonian parameters as the massless limit is approached. This smooth behaviour in gamma is due to the discontinuities from each of the two different massive modes cancelling each other out. Such cancellations cannot occur in special cases with only one massive mode, such as the Isham-Salam-Strathdee theory.

A possible solution to the dark energy problem is that Einstein's theory of general relativity is modified. A suite of models have been proposed that, in general, are unable to predict the correct amount of large scale structure in the distribution of galaxies or anisotropies in the Cosmic Microwave Background. It has been argued, however, that it should be possible to constrain a general class of theories of modified gravity by focusing on properties such as the growing mode, gravitational slip and the effective, time varying Newton's constant. We show that assuming certain physical requirements such as stability, metricity and gauge invariance, it is possible to come up with consistency conditions between these various parameters. In this paper we focus on theories which have, at most, 2nd derivatives in the metric variables and find restrictions that shed light on current and future experimental constraints without having to resort to a (as yet unknown) complete theory of modified gravity. We claim that future measurements of the growth of structure on small scales (i.e. from 1-200 h^{-1} Mpc) may lead to tight constraints on both dark energy and modified theories of gravity.

It has been argued that a Universe governed by Eddington-Born-Infeld gravity can be compatible with current cosmological constraints. The extra fields introduced in this theory can behave both as dark matter and dark energy, unifying the dark sector in one coherent framework. We show the various roles the extra fields can play in the expansion of the Universe and study the evolution of linear perturbations in the various regimes. We find that, as a unified theory of the dark sector, Eddington-Born-Infeld gravity will lead to excessive fluctuations in the Cosmic Microwave Background on large scales. In the presence of a cosmological constant, however, the extra fields can behave as a form of non-particulate dark matter and can lead to a cosmology which is entirely compatible with current observations of large scale structure. We discuss the interpretation of this form of dark matter and how it can differ from standard, particulate dark matter.

We show that a class of bi-gravity theories contain solutions describing dark matter. A particular member of this class is also shown to be equivalent to the Eddington-Born-Infeld gravity, recently proposed as a candidate for dark matter. Bigravity theories also have cosmological de Sitter backgrounds and we find solutions interpolating between matter and acceleration eras.

Over the last few decades, astronomers and cosmologists have accumulated vast amounts of data clearly demonstrating that our current theories of fundamental particles and of gravity are inadequate to explain the observed discrepancy between the dynamics and the distribution of the visible matter in the Universe. The Modified Newtonian Dynamics (MOND) proposal aims at solving the problem by postulating that Newton's second law of motion is modified for accelerations smaller than ~10^{-10}m/s^2. This simple amendment, has had tremendous success in explaining galactic rotation curves. However, being non-relativistic, it cannot make firm predictions for cosmology. A relativistic theory called Tensor-Vector-Scalar (TeVeS) has been proposed by Bekenstein building on earlier work of Sanders which has a MOND limit for non-relativistic systems. In this article I give a short introduction to TeVeS theory and focus on its predictions for cosmology as well as some non-cosmological studies.

We show that models of new particle physics containing massless pseudoscalar fields super-weakly coupled to photons can be very efficiently probed with CMB polarization anisotropies. The stochastic pseudoscalar fluctuations generated during inflation provide a mechanism for converting E-mode polarization to B-mode during photon propagation from the surface of last scattering. The efficiency of this conversion process is controlled by the dimensionless ratio H/(2\pi f_a), where H is the Hubble scale during inflation, and f_a^{-1} is the strength of the pseudoscalar coupling to photons. The current observational limits on the B-mode constrain this ratio to be less than 0.07, which in many models of inflation translates to a sensitivity to values of f_a in excess of 10^{14} GeV, surpassing the sensitivity of other tests.

General Relativity (GR) is a phenomenologically successful theory that rests on firm foundations, but has not been tested on cosmological scales. The advent of dark energy (and possibly even the requirement of cold dark matter), has increased the need for testing modifications to GR, as the inference of such otherwise undetected fluids, depends crucially on the theory of gravity. In this work I outline a general scheme for constructing consistent and covariant modifications to the Einstein equations. This framework is such that there is a clear connection between the modification and the underlying field content that produces it. I conclude by a simple metric based modification of the fluctuation equations for which the background is exact Lambda-CDM and present its impact on observations of the cosmic microwave background radiation.

It has been suggested that Einstein's theory of General Relativity can be modified to accomodate mismatches between the gravitational field and luminous matter on a wide range of scales. Covariant theories of modified gravity generically predict the existence of extra degrees of freedom which may be interpreted as dark matter. We study a subclass of these theories where the overall energy density in these extra degrees of freedom is subdominant relative to the baryon density and show that they favour the presence of massive neutrinos. For some specific cases (such as a flat Universes with a cosmological constant) one finds a conservative lower bound on the neutrinos mass of $m_\nu>0.31$ eV.

I consider an extented version of Bekenstein's Tensor-Vector-Scalar theory where the action of the vector field is of a general Einstein-Ether form. This work presents the cosmological equations of this theory, both at the background and perturbed level, for scalar, vector and tensor perturbation modes. By solving the background equations in the radiation era analytically, to an excellent approximation, I construct the primordial adiabatic perturbation for a general family of scalar field kinetic functions.

We study the background cosmology governed by the Tensor-Vector-Scalar theory of gravity proposed by Bekenstein. We consider a broad family of potentials that lead to modified gravity and calculate the evolution of the field variables both numerically and analytically. We find a range of possible behaviors, from scaling to the late time domination of either the additional gravitational degrees of freedom or the background fluid.

A relativistic theory of modified gravity has been recently proposed by Bekenstein. The tensor field in Einstein's theory of gravity is replaced by a scalar, a vector, and a tensor field which interact in such a way to give Modified Newtonian Dynamics (MOND) in the weak-field non-relativistic limit. We study the evolution of the universe in such a theory, identifying its key properties and comparing it with the standard cosmology obtained in Einstein gravity. The evolution of the scalar field is akin to that of tracker quintessence fields. We expand the theory to linear order to find the evolution of perturbations on large scales. The impact on galaxy distributions and the cosmic microwave background is calculated in detail. We show that it may be possible to reproduce observations of the cosmic microwave background and galaxy distributions with Bekenstein's theory of MOND.

A relativistic theory of gravity has recently been proposed by Bekenstein, where gravity is mediated by a tensor, a vector and a scalar field, thus called TeVeS. The theory aims at modifying gravity in such a way as to reproduce Milgrom's MOdified Newtonian Dynamics (MOND) in the weak field, non-relativistic limit, which provides a framework to solve the missing mass problem in galaxies without invoking dark matter. In this paper I apply a covariant approach to formulate the cosmological equations for this theory, for both the background and linear perturbations. I derive the necessary perturbed equations for scalar, vector and tensor modes without adhering to a particular gauge. Special gauges are considered in the appendix.

The Cosmic Microwave Background (CMB) anisotropy constrains the geometry of the Universe because the positions of the acoustic peaks of the angular power spectrum depend strongly on the curvature of underlying three-dimensional space. In this Letter we exploit current observations to determine the spatial geometry of the Universe in the presence of isocurvature modes. Previous analyses have always assumed that the cosmological perturbations were initially adiabatic. A priori one might expect that allowing additional isocurvature modes would substantially degrade the constraints on the curvature of the Universe. We find, however, that if one considers additional data sets, the geometry remains well constrained. When the most general isocurvature perturbation is allowed, the CMB alone can only poorly constrain the geometry to Omega_0=1.6+-0.3. Including large-scale structure (LSS) data one obtains Omega_0=1.07+-0.03, and Omega_0=1.06+-0.02 when supplemented by the Hubble Space Telescope (HST) Key Project determination of H_0 and SNIa data.

We investigate the constraints imposed by the first-year WMAP CMB data extended to higher multipole by data from ACBAR, BOOMERANG, CBI and the VSA and by the LSS data from the 2dF galaxy redshift survey on the possible amplitude of primordial isocurvature modes. A flat universe with CDM and Lambda is assumed, and the baryon, CDM (CI), and neutrino density (NID) and velocity (NIV) isocurvature modes are considered. Constraints on the allowed isocurvature contributions are established from the data for various combinations of the adiabatic mode and one, two, and three isocurvature modes, with intermode cross-correlations allowed. Since baryon and CDM isocurvature are observationally virtually indistinguishable, these modes are not considered separately. We find that when just a single isocurvature mode is added, the present data allows an isocurvature fraction as large as 13+-6, 7+-4, and 13+-7 percent for adiabatic plus the CI, NID, and NIV modes, respectively. When two isocurvature modes plus the adiabatic mode and cross-correlations are allowed, these percentages rise to 47+-16, 34+-12, and 44+-12 for the combinations CI+NID, CI+NIV, and NID+NIV, respectively. Finally, when all three isocurvature modes and cross-correlations are allowed, the admissible isocurvature fraction rises to 57+-9 per cent. The sensitivity of the results to the choice of prior probability distribution is examined.

Markov Chain Monte Carlo (MCMC) techniques are now widely used for cosmological parameter estimation. Chains are generated to sample the posterior probability distribution obtained following the Bayesian approach. An important issue is how to optimize the efficiency of such sampling and how to diagnose whether a finite-length chain has adequately sampled the underlying posterior probability distribution. We show how the power spectrum of a single such finite chain may be used as a convergence diagnostic by means of a fitting function, and discuss strategies for optimizing the distribution for the proposed steps. The methods developed are applied to current CMB and LSS data interpreted using both a pure adiabatic cosmological model and a mixed adiabatic/isocurvature cosmological model including possible correlations between modes. For the latter application, because of the increased dimensionality and the presence of degeneracies, the need for tuning MCMC methods for maximum efficiency becomes particularly acute.

We present a new method for measuring the CMB temperature quadrupole, using large scale CMB polarization. The method exploits the fact that CMB polarization is partially sourced by the local temperature quadrupole. We link the temperature with the polarization spectrum directly by relating the local quadrupole at the onset of reionization to both of them. The dominant contribution is at l<30 and since we use many l values, we can reduce the error significantly below cosmic variance. In particular, for our fiducial model, the error on the temperature quadrupole is reduced to 24%. This has the potential of reducing the probability of a low quadrupole by two orders of magnitude.

We investigate the constraints imposed by the current data on correlated mixtures of adiabatic and non-adiabatic primordial perturbations. We discover subtle flat directions in parameter space that tolerate large (~60%) contributions of non-adiabatic fluctuations. In particular, larger values of the baryon density and a spectral tilt are allowed. The cancellations in the degenerate directions are explored and the role of priors elucidated.

The recent discovery of a Gunn--Peterson (GP) trough in the spectrum of the redshift 6.28 SDSS quasar has raised the tantalizing possibility that we have detected the reionization of the universe. However, a neutral fraction (of hydrogen) as small as 0.1% is sufficient to cause the GP trough, hence its detection alone cannot rule out reionization at a much earlier epoch. The Cosmic Microwave Background (CMB) polarization anisotropy offers an alternative way to explore the dark age of the universe. We show that for most models constrained by the current CMB data and by the discovery of a GP trough (showing that reionization occurred at z > 6.3), MAP can detect the reionization signature in the polarization power spectrum. The expected 1-sigma error on the measurement of the electron optical depth is around 0.03 with a weak dependence on the value of that optical depth. Such a constraint on the optical depth will allow MAP to achieve a 1-sigma error on the amplitude of the primordial power spectrum of 6%. MAP with two years (Planck with one year) of observation can distinguish a model with 50% (6%) partial ionization between redshifts of 6.3 and 20 from a model in which hydrogen was completely neutral at redshifts greater than 6.3. Planck will be able to distinguish between different reionization histories even when they imply the same optical depth to electron scattering for the CMB photons.

We have developed a fast method for predicting the angular power spectrum, C_l, of the cosmic microwave background given cosmological parameters and a primordial power spectrum of perturbations. After pre--computing the radiation temperature and gravitational potential transfer functions over a small sub--space of the total model parameter space, the rest of the model space (six or more cosmological parameters and arbitrarily many primordial power spectrum parameters) is reached via rapid analytic and semi--analytic approximations which are highly accurate on all angular scales for which linear perturbation theory applies. A single power spectrum can be calculated in ~ 1 second on a desktop computer. We discuss applications to cosmological parameter estimation.

If Omega_tot = 1 and structure formed from adiabatic initial conditions then the age of the Universe, as constrained by measurements of the cosmic microwave background (CMB), is t_0=14.0 +/- 0.5 Gyr. The uncertainty is surprisingly small given that CMB data alone constrain neither h nor Omega_Lambda significantly. It is due to the tight (and accidental) correlation, in these models, of the age with the angle subtended by the sound horizon on the last--scattering surface and thus with the well-determined acoustic peak locations. If we assume either the HST Key Project result h = 0.72 \pm .08 or simply that h > 0.55, we find Omega_Lambda > 0.4 at 95% confidence--another argument for dark energy, independent of supernovae observations. Our analysis is greatly simplified by the Monte Carlo Markov chain approach to Bayesian inference combined with a fast method for calculating angular power spectra.

We show how the presence of a very light scalar with a cubic self-interaction in six dimensions can stabilize the extra dimensions at radii which are naturally exponentially large, $r \sim \ell \exp [(4\pi)^3/g^2]$, where $\ell$ is a microscopic physics scale and $g$ is the (dimensionless) cubic coupling constant. The resulting radion mode of the metric becomes a very light degree of freedom whose mass, $m \sim 1/(M_p r^2)$ is stable under radiative corrections. For $1/r \sim 10^{-3}$ eV the radion is extremely light, $m \sim 10^{-33}$ eV. Its couplings cause important deviations from General Relativity in the very early universe, but naturally evolve to phenomenologically acceptable values at present.

We examine the late-time (nucleosynthesis and later) cosmological implications of brane-world scenarios having large (millimeter sized) extra dimensions. In particular, recent proposals for understanding why the extra dimensions are so large in these models indicate that moduli like the radion appear (to four-dimensional observers) to be extremely light, with a mass of order 10^{-33} eV, allowing them to play the role of the light scalar of quintessence models. The radion-as-quintessence solves a long-standing problem since its small mass is technically natural, in that it is stable against radiative corrections. Its challenges are to explain why such a light particle has not been seen in precision tests of gravity, and why Newton's constant has not appreciably evolved since nucleosynthesis. We find the couplings suggested by stabilization models can provide explanations for both of these questions. We identify the features which must be required of any earlier epochs of cosmology in order for these explanations to hold.

In a recent paper we considered the possibility of a scalar field providing an explanation for the cosmic acceleration. Our model had the interesting properties of attractor-like behavior and having its parameters of O(1) in Planck units. Here we discuss the effect of the field on large scale structure and CMB anisotropies. We show how some versions of our model inspired by "brane" physics have novel features due to the fact that the scalar field has a significant role over a wider range of redshifts than for typical "dark energy" models. One of these features is the additional suppression of the formation of large scale structure, as compared with cosmological constant models. In light of the new pressures being placed on cosmological parameters (in particular H_0) by CMB data, this added suppression allows our "brane" models to give excellent fits to both CMB and large scale structure data.

Modern data is showing increasing evidence that the Universe is accelerating. So far, all attempts to account for the acceleration have required some fundamental dimensionless quantities to be extremely small. We show how a class of scalar field models (which may emerge naturally from superstring theory) can account for acceleration which starts in the present epoch with all the potential parameters O(1) in Planck units.