We outline the experimental concept and key scientific capabilities of AION (Atom Interferometer Observatory and Network), a proposed UK-based experimental programme using cold strontium atoms to search for ultra-light dark matter, to explore gravitational waves in the mid-frequency range between the peak sensitivities of the LISA and LIGO/Virgo/ KAGRA/INDIGO/Einstein Telescope/Cosmic Explorer experiments, and to probe other frontiers in fundamental physics. AION would complement other planned searches for dark matter, as well as probe mergers involving intermediate mass black holes and explore early universe cosmology. AION would share many technical features with the MAGIS experimental programme in the US, and synergies would flow from operating AION in a network with this experiment, as well as with other atom interferometer experiments such as MIGA, ZAIGA and ELGAR. Operating AION in a network with other gravitational wave detectors such as LIGO, Virgo and LISA would also offer many synergies.

Dark matter (DM) can consist of very light bosons behaving as a classical scalar field that experiences coherent oscillations. The presence of this DM field would perturb the dynamics of celestial bodies, either because the (oscillating) DM stress tensor modifies the gravitational potentials of the galaxy, or if DM is directly coupled to the constituents of the body. We study secular variations of the orbital parameters of binary systems induced by such perturbations. Two classes of effects are identified. Effects of the first class appear if the frequency of DM oscillations is in resonance with the orbital motion; these exist for general DM couplings including the case of purely gravitational interaction. Effects of the second class arise if DM is coupled quadratically to the masses of the binary system members and do not require any resonant condition. The exquisite precision of binary pulsar timing can be used to constrain these effects. Current observations are not sensitive to oscillations in the galactic gravitational field, though a discovery of pulsars in regions of high DM density may improve the situation. For DM with direct coupling to ordinary matter, the current timing data are already competitive with other existing constraints in the range of DM masses $\sim 10^{-22}-10^{-18}\,{\rm eV}$. Future observations are expected to increase the sensitivity and probe new regions of parameters.

New light states thermally coupled to the Standard Model plasma alter the expansion history of the Universe and impact the synthesis of the primordial light elements. In this work, we carry out an exhaustive and precise analysis of the implications of MeV-scale BSM particles in Big Bang Nucleosynthesis (BBN) and for Cosmic Microwave Background (CMB) observations. We find that, BBN observations set a lower bound on the thermal dark matter mass of $m_\chi > 0.4\,\text{MeV}$ at $2\sigma$. This bound is independent of the spin and number of internal degrees of freedom of the particle, of the annihilation being s-wave or p-wave, and of the annihilation final state. Furthermore, we show that current BBN plus CMB observations constrain purely electrophilic and neutrinophilic BSM species to have a mass, $m_\chi > 3.7\,\text{MeV}$ at $2\sigma$. We explore the reach of future BBN measurements and show that upcoming CMB missions should improve the bounds on light BSM thermal states to $m_\chi > (10-15)\,\text{MeV}$. Finally, we demonstrate that very light BSM species thermally coupled to the SM plasma are highly disfavoured by current cosmological observations.

We study the propagation of light in the presence of a parity-violating coupling between photons and axion-like particles (ALPs). Naively, this interaction could lead to a split of light rays into two separate beams of different polarization chirality and with different refraction angles. However, by using the eikonal method we explicitly show that this is not the case and that ALP clumps do not produce any spatial birefringence. This happens due to non-trivial variations of the photon's frequency and wavevector, which absorb time-derivatives and gradients of the ALP field. We argue that these variations represent a new way to probe the ALP-photon couping with precision frequency measurements.

We propose in this White Paper a concept for a space experiment using cold atoms to search for ultra-light dark matter, and to detect gravitational waves in the frequency range between the most sensitive ranges of LISA and the terrestrial LIGO/Virgo/KAGRA/INDIGO experiments. This interdisciplinary experiment, called Atomic Experiment for Dark Matter and Gravity Exploration (AEDGE), will also complement other planned searches for dark matter, and exploit synergies with other gravitational wave detectors. We give examples of the extended range of sensitivity to ultra-light dark matter offered by AEDGE, and how its gravitational-wave measurements could explore the assembly of super-massive black holes, first-order phase transitions in the early universe and cosmic strings. AEDGE will be based upon technologies now being developed for terrestrial experiments using cold atoms, and will benefit from the space experience obtained with, e.g., LISA and cold atom experiments in microgravity. This paper is based on a submission (v1) in response to the Call for White Papers for the Voyage 2050 long-term plan in the ESA Science Programme. ESA limited the number of White Paper authors to 30. However, in this version (v2) we have welcomed as supporting authors participants in the Workshop on Atomic Experiments for Dark Matter and Gravity Exploration held at CERN: ({\tt https://indico.cern.ch/event/830432/}), as well as other interested scientists, and have incorporated additional material.

Particle dark matter could have a mass anywhere from that of ultralight candidates, $m_\chi\sim 10^{-21}\,$eV, to scales well above the GeV. Conventional laboratory searches are sensitive to a range of masses close to the weak scale, while new techniques are required to explore candidates outside this realm. In particular lighter candidates are difficult to detect due to their small momentum. Here we study two experimental set-ups which {\it do not require transfer of momentum} to detect dark matter: atomic clocks and co-magnetometers. These experiments probe dark matter that couples to the spin of matter via the very precise measurement of the energy difference between atomic states of different angular momenta. This coupling is possible (even natural) in most dark matter models, and we translate the current experimental sensitivity into implications for different dark matter models. It is found that the constraints from current atomic clocks and co-magnetometers can be competitive in the mass range $m_\chi\sim 10^{-21}-10^3\,$eV, depending on the model. We also comment on the (negligible) effect of different astrophysical neutrino backgrounds.

We derive novel constraints on millicharged dark matter and ultralight axion-like particles using pulsar timing and fast radio burst observations. Millicharged dark matter affects the dispersion measure of the time of arrival of radio pulses in a way analogous to free electrons. Light pseudo-scalar dark matter, on the other hand, causes the polarization angle of radio signals to oscillate. We show that current and future data can set strong constraints in both cases. For dark matter particles of charge $\epsilon e$, these constraints are ${\epsilon}/{m_{\rm milli}} \lesssim 10^{-8}{\rm eV}^{-1}$, for masses $m_{\rm milli}\gtrsim 10^{-6}\,$eV. For axion-like particles, the analysis of signals from pulsars yields constraints in the axial coupling of the order of $g/m_a\lesssim 10^{-13} {\rm GeV}^{-1}/(10^{-22}{\rm eV})$. Both bounds scale as $(\rho/\rho_{\rm dm})^{1/2}$ if the energy density $\rho$ of the components is a fraction of the total dark matter energy density $\rho_{\rm dm}$. We do a detailed study of both effects using data from two samples of pulsars in the galaxy and in globular clusters, as well as data from FRB 121102 and PSR J0437$-$4715. We show that in both cases actual pulsar data constrain a new region of the parameter space for these models, and will improve with future pulsar-timing observations.

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.

Bosonic ultra-light dark matter (ULDM) would form cored density distributions at the center of galaxies. These cores, seen in numerical simulations, admit analytic description as the lowest energy bound state solution ("soliton") of the Schroedinger-Poisson equations. Numerical simulations of ULDM galactic halos found empirical scaling relations between the mass of the large-scale host halo and the mass of the central soliton. We discuss how the simulation results of different groups can be understood in terms of the basic properties of the soliton. Importantly, simulations imply that the energy per unit mass in the soliton and in the virialised host halo should be approximately equal. This relation lends itself to observational tests, because it predicts that the peak circular velocity, measured for the host halo in the outskirts of the galaxy, should approximately repeat itself in the central region. Contrasting this prediction to the measured rotation curves of well-resolved near-by galaxies, we show that ULDM in the mass range $m\sim (10^{-22}\div 10^{-21})$ eV, which has been invoked as a possible solution to the small-scale puzzles of $\Lambda$CDM, is in tension with the data. We suggest that a dedicated analysis of the Milky Way inner gravitational potential could probe ULDM up to $m\lesssim 10^{-19}$ eV.

We present a detailed analysis of the effect of light Dark Matter (DM) on atomic clocks, for the case where DM mass and density are such that occupation numbers are low and DM must be considered as particles scattering off the atoms, rather than a classical field. We show that the resulting atomic clock frequency shifts are first order in the scattering amplitudes, and particularly suited to constrain DM models in the regime where the DM mass $m_\chi \ll$ GeV. We provide some rough order of magnitude estimates of sensitivity that can be confronted to any DM model that allows for non zero differential scattering amplitudes of the two atomic states involved in the clock.

We consider the scenario where dark matter (DM) is represented by an ultra-light classical scalar field performing coherent periodic oscillations. We point out that such DM perturbs the dynamics of binary systems either through its gravitational field or via direct coupling to ordinary matter. This perturbation gets resonantly amplified if the frequency of DM oscillations is close to a (half-)integer multiple of the orbital frequency of the system and leads to a secular variation of the orbital period. We suggest to use binary pulsars as probes of this scenario and estimate their sensitivity. While the current accuracy of observations is not yet sufficient to probe the purely gravitational effect of DM, it already yields constraints on direct coupling that are competitive with other bounds. The sensitivity will increase with the upcoming radio observatories such as Square Kilometer Array.

As a binary pulsar moves through a wind of dark matter particles, the resulting dynamical friction modifies the binary's orbit. We study this effect for the double disk dark matter (DDDM) scenario, where a fraction of the dark matter is dissipative and settles into a thin disk. For binaries within the dark disk, this effect is enhanced due to the higher dark matter density and lower velocity dispersion of the dark disk, and due to its co-rotation with the baryonic disk.We estimate the effect and compare it with observations for two different limits in the Knudsen number ($Kn$). First, in the case where DDDM is effectively collisionless within the characteristic scale of the binary ($Kn\gg1$) and ignoring the possible interaction between the pair of dark matter wakes. Second, in the fully collisional case ($Kn\ll1$), where a fluid description can be adopted and the interaction of the pair of wakes is taken into account. We find that the change in the orbital period is of the same order of magnitude in both limits. A comparison with observations reveals good prospects to probe currently allowed DDDM models with timing data from binary pulsars in the near future. We finally comment on the possibility of extending the analysis to the intermediate (rarefied gas) case with $Kn\sim1$.

Binary pulsars are ideal to test the foundations of General Relativity, such as Lorentz symmetry, which requires that experiments produce the same results in all free-falling (i.e.inertial) frames. We here break this symmetry in the gravitational sector by specifying a preferred time direction, and thus a preferred frame, at each spacetime point. We then examine the consequences of this gravitational Lorentz symmetry breaking in the orbital evolution of binary pulsars, focusing on the dissipative effects. We find that Lorentz symmetry breaking modifies these effects, and thus the orbital dynamics, in two different ways. First, it generically causes the emission of dipolar radiation, which makes the orbital separation decrease faster than in General Relativity. Second, the quadrupole component of the emission is also modified. The orbital evolution depends critically on the sensitivities of the stars, which measure how their binding energies depend on the motion relative to the preferred frame. We calculate the sensitivities numerically and compute the predicted orbital decay rate of binary pulsars in Lorentz-violating gravity. By testing these predictions against observations, we place very stringent constraints on gravitational Lorentz violation.

We develop the framework for testing Lorentz invariance in the dark matter sector using galactic dynamics. We consider a Lorentz violating (LV) vector field acting on the dark matter component of a satellite galaxy orbiting in a host halo. We introduce a numerical model for the dynamics of satellites in a galactic halo and for a galaxy in a rich cluster to explore observational consequences of such an LV field. The orbital motion of a satellite excites a time dependent LV force which greatly affects its internal dynamics. Our analysis points out key observational signatures which serve as probes of LV forces. These include modifications to the line of sight velocity dispersion, mass profiles and shapes of satellites. With future data and a more detailed modeling these signatures can be exploited to constrain a new region of the parameter space describing the LV in the dark matter sector.

We find multi-scalar effective field theories (EFTs) that can achieve a slow inflationary roll despite having a scalar potential that does not satisfy the usual slow-roll condition (d V)^2 << V^2/Mp^2. They evade the usual slow-roll conditions on $V$ because their kinetic energies are dominated by single-derivative terms rather than the usual two-derivative terms. Single derivatives dominate during slow roll and so do not require a breakdown of the usual derivative expansion that underpins calculational control in much of cosmology. The presence of such terms requires some sort of UV Lorentz-symmetry breaking during inflation (besides the usual cosmological breaking). Chromo-natural inflation provides an example of a UV theory that can generate the multi-field single-derivative terms we consider, and we argue that the EFT we find indeed captures the slow-roll conditions for the background evolution for Chromo-natural inflation. We also show that our EFT can be understood as a multi-field generalization of the single-field Cuscuton models. The multi-field case introduces a new feature, however: the scalar kinetic terms define a target-space 2-form, F_{ab}, whose antisymmetry gives new ways for slow roll to be achieved.

We present a new analytic approach to describe large scale structure formation in the mildly non-linear regime. The central object of the method is the time-dependent probability distribution function generating correlators of the cosmological observables at a given moment of time. Expanding the distribution function around the Gaussian weight we formulate a perturbative technique to calculate non-linear corrections to cosmological correlators, similar to the diagrammatic expansion in a three-dimensional Euclidean quantum field theory, with time playing the role of an external parameter. For the physically relevant case of cold dark matter in an Einstein--de Sitter universe, the time evolution of the distribution function can be found exactly and is encapsulated by a time-dependent coupling constant controlling the perturbative expansion. We show that all building blocks of the expansion are free from spurious infrared enhanced contributions that plague the standard cosmological perturbation theory. This paves the way towards the systematic resummation of infrared effects in large scale structure formation. We also argue that the approach proposed here provides a natural framework to account for the influence of short-scale dynamics on larger scales along the lines of effective field theory.

We use time-sliced perturbation theory (TSPT) to give an accurate description of the infrared non-linear effects affecting the baryonic acoustic oscillations (BAO) present in the distribution of matter at very large scales. In TSPT this can be done via a systematic resummation that has a simple diagrammatic representation and does not involve uncontrollable approximations. We discuss the power counting rules and derive explicit expressions for the resummed matter power spectrum up to next-to leading order and the bispectrum at the leading order. The two-point correlation function agrees well with N-body data at BAO scales. The systematic approach also allows to reliably assess the shift of the baryon acoustic peak due to non-linear effects.

We point out that the observed time delay between the detection of the signal at the Hanford and Livingston LIGO sites from the gravitational wave event GW150914 places an upper bound on the speed of propagation of gravitational waves, $c_{gw}\lesssim 1.7$ in the units of speed of light. Combined with the lower bound from the absence of gravitational Cherenkov losses by cosmic rays that rules out most of subluminal velocities, this gives a model-independent double-sided constraint $1\lesssim c_{gw}\lesssim 1.7$. We compare this result to model-specific constraints from pulsar timing and cosmology.

Cosmological perturbations of sufficiently long wavelength admit a fluid dynamic description. We consider modes with wavevectors below a scale $k_m$ for which the dynamics is only mildly non-linear. The leading effect of modes above that scale can be accounted for by effective non-equilibrium viscosity and pressure terms. For mildly non-linear scales, these mainly arise from momentum transport within the ideal and cold but inhomogeneous fluid, while momentum transport due to more microscopic degrees of freedom is suppressed. As a consequence, concrete expressions with no free parameters, except the matching scale $k_m$, can be derived from matching evolution equations to standard cosmological perturbation theory. Two-loop calculations of the matter power spectrum in the viscous theory lead to excellent agreement with $N$-body simulations up to scales $k=0.2 \, h/$Mpc. The convergence properties in the ultraviolet are better than for standard perturbation theory and the results are robust with respect to variations of the matching scale.

Theories of gravitation without Lorentz invariance are candidates of low-energy descriptions of quantum gravity. In this review we will describe the phenomenological consequences of the candidates associated to the existence of a preferred time direction

We consider a scenario where local Lorentz invariance is violated by the existence of a preferred time direction at every space-time point. This scenario can arise in the context of quantum gravity and its description at low energies contains a unit time-like vector field which parameterizes the preferred direction. The particle physics tests of Lorentz invariance preclude a direct coupling of this vector to the fields of the Standard Model, but do not bear implications for dark matter. We discuss how the presence of this vector and its possible coupling to dark matter affect the evolution of the Universe. At the level of homogeneous cosmology the only effect of Lorentz invariance violation is a rescaling of the expansion rate. The physics is richer at the level of perturbations. We identify three effects crucial for observations: the rescaling of the matter contribution to the Poisson equation, the appearance of an extra contribution to the anisotropic stress and the scale-dependent enhancement of dark matter clustering. These effects result in distinctive features in the power spectra of the CMB and density fluctuations. Making use of the data from Planck and WiggleZ we obtain the most stringent cosmological constraints to date on departures from Lorentz symmetry. Our analysis provides the first direct bounds on deviations from Lorentz invariance in the dark matter sector.

Theories with massive gravitons are interesting for a variety of physical applications, ranging from cosmological phenomena to holographic modeling of condensed matter systems. To date, they have been formulated as effective field theories with a cutoff proportional to a positive power of the graviton mass m_g and much smaller than that of the massless theory (M_P ~ 10^19 GeV in the case of general relativity). In this paper we present an ultraviolet completion for massive gravity valid up to a high energy scale independent of the graviton mass. The construction is based on the existence of a preferred time foliation combined with spontaneous condensation of vector fields. The perturbations of these fields are massive and below their mass the theory reduces to a model of Lorentz violating massive gravity. The latter theory possesses instantaneous modes whose consistent quantization we discuss in detail. We briefly study some modifications to gravitational phenomenology at low-energies. The homogeneous cosmological solutions are the same as in the standard cosmology. The gravitational potential of point sources agrees with the Newtonian one at distances small with respect to m_g^(-1). Interestingly, it becomes repulsive at larger distances.

We compute non-linear corrections to the matter power spectrum taking the time- and scale-dependent free-streaming length of neutrinos into account. We adopt a hybrid scheme that matches the full Boltzmann hierarchy to an effective two-fluid description at an intermediate redshift. The non-linearities in the neutrino component are taken into account by using an extension of the time-flow framework. We point out that this remedies a spurious behaviour that occurs when neglecting non-linear terms for neutrinos. This behaviour is related to how efficiently short modes decouple from long modes and can be traced back to the violation of momentum conservation if neutrinos are treated linearly. Furthermore, we compare our results at next to leading order to various other methods and quantify the accuracy of the fluid description. Due to the correct decoupling behaviour of short modes, the two-fluid scheme is a suitable starting point to compute higher orders in perturbations or for resummation methods.

Binary pulsars are excellent laboratories to test the building blocks of Einstein's theory of General Relativity. One of these is Lorentz symmetry which states that physical phenomena appear the same for all inertially moving observers. We study the effect of violations of Lorentz symmetry in the orbital evolution of binary pulsars and find that it induces a much more rapid decay of the binary's orbital period due to the emission of dipolar radiation. The absence of such behavior in recent observations allows us to place the most stringent constraints on Lorentz violation in gravity, thus verifying one of the cornerstones of Einstein's theory much more accurately than any previous gravitational observation.

We analyze the dark matter power spectrum at three-loop order in standard perturbation theory of large scale structure. We observe that at late times the loop expansion does not converge even for large scales (small momenta) well within the linear regime, but exhibits properties compatible with an asymptotic series. We propose a technique to restore the convergence in the limit of small momentum, and use it to obtain a perturbative expansion with improved convergence for momenta in the range where baryonic acoustic oscillations are present. Our results are compared with data from N-body simulations at different redshifts, and we find good agreement within this range.

We discuss the convergence of cosmological perturbation theory. We prove that the polynomial enhancement of the non-linear corrections expected from the effects of soft modes is absent in equal-time correlators like the power or bispectrum. We first show this at leading order by resumming the most important corrections of soft modes to an arbitrary skeleton of hard fluctuations. We derive the same result in the eikonal approximation, which also allows us to show the absence of enhancement at any order. We complement the proof by an explicit calculation of the power spectrum at two-loop order, and by further numerical checks at higher orders. Using these insights, we argue that the modification of the power spectrum from soft modes corresponds at most to logarithmic corrections. Finally, we discuss the asymptotic behavior in the large and small momentum regimes and identify the expansion parameter pertinent to non-linear corrections.

The role of Lorentz invariance as a fundamental symmetry of nature has been lately reconsidered in different approaches to quantum gravity. It is thus natural to study whether other puzzles of physics may be solved within these proposals. This may be the case for the cosmological constant problem. Indeed, it has been shown that breaking Lorentz invariance provides Lagrangians that can drive the current acceleration of the universe without experiencing large corrections from ultraviolet physics. In this work, we focus on the simplest model of this type, called ThetaCDM, and study its cosmological implications in detail. At the background level, this model cannot be distinguished from LambdaCDM. The differences appear at the level of perturbations. We show that in ThetaCDM, the spectrum of CMB anisotropies and matter fluctuations may be affected by a rescaling of the gravitational constant in the Poisson equation, by the presence of extra contributions to the anisotropic stress, and finally by the existence of extra clustering degrees of freedom. To explore these modifications accurately, we modify the Boltzmann code CLASS. We then use the parameter inference code Monte Python to confront ThetaCDM with data from WMAP-7, SPT and WiggleZ. We obtain strong bounds on the parameters accounting for deviations from LambdaCDM. In particular, we find that the discrepancy between the gravitational constants appearing in the Poisson and Friedmann equations is constrained at the level 1.8%.

We construct a model of dark energy with a technically natural small contribution to cosmic acceleration, i.e. this contribution does not receive corrections from other scales in the theory. The proposed acceleration mechanism appears generically in the low-energy limit of gravity theories with violation of Lorentz invariance that contain a derivatively coupled scalar field Theta. The latter may be the Goldstone field of a broken global symmetry. The model, that we call Theta-CDM, is a valid effective field theory up to a high cutoff just a few orders of magnitude below the Planck scale. Furthermore, it can be ultraviolet-completed in the context of Horava gravity. We discuss the observational predictions of the model. Even in the absence of a cosmological constant term, the expansion history of the Universe is essentially indistinguishable from that of Lambda-CDM. The difference between the two theories appears at the level of cosmological perturbations. We find that in Theta-CDM the matter power spectrum is enhanced at subhorizon scales compared to Lambda-CDM. This property can be used to discriminate the model from Lambda-CDM with current cosmological data.

We study the possibility to constrain deviations from Lorentz invariance in dark matter (DM) with cosmological observations. Breaking of Lorentz invariance generically introduces new light gravitational degrees of freedom, which we represent through a dynamical timelike vector field. If DM does not obey Lorentz invariance, it couples to this vector field. We find that this coupling affects the inertial mass of small DM halos which no longer satisfy the equivalence principle. For large enough lumps of DM we identify a (chameleon) mechanism that restores the inertial mass to its standard value. As a consequence, the dynamics of gravitational clustering are modified. Two prominent effects are a scale dependent enhancement in the growth of large scale structure and a scale dependent bias between DM and baryon density perturbations. The comparison with the measured linear matter power spectrum in principle allows to bound the departure from Lorentz invariance of DM at the per cent level.

We study the radiation of gravitational waves by self-gravitating binary systems in the low-energy limit of Horava gravity. We find that the predictions for the energy-loss formula of General Relativity are modified already for Newtonian sources: the quadrupole contribution is altered, in part due to the different speed of propagation of the tensor modes; furthermore, there is a monopole contribution stemming from an extra scalar degree of freedom. A dipole contribution only appears at higher post-Newtonian order. We use these findings to constrain the low-energy action of Horava gravity by comparing them with the radiation damping observed for binary pulsars. Even if this comparison is not completely appropriate - since compact objects cannot be described within the post-Newtonian approximation - it represents an order of magnitude estimate. In the limit where the post-Newtonian metric coincides with that of General Relativity, our energy-loss formula provides the strongest constraints for Horava gravity at low-energies.

Boltzmann codes are used extensively by several groups for constraining cosmological parameters with Cosmic Microwave Background and Large Scale Structure data. This activity is computationally expensive, since a typical project requires from 10'000 to 100'000 Boltzmann code executions. The newly released code CLASS (Cosmic Linear Anisotropy Solving System) incorporates improved approximation schemes leading to a simultaneous gain in speed and precision. We describe here the three approximations used by CLASS for basic LambdaCDM models, namely: a baryon-photon tight-coupling approximation which can be set to first order, second order or to a compromise between the two; an ultra-relativistic fluid approximation which had not been implemented in public distributions before; and finally a radiation streaming approximation taking reionisation into account.

Horava's proposal for non-relativistic quantum gravity introduces a preferred time foliation of space-time which violates the local Lorentz invariance. The foliation is encoded in a dynamical scalar field which we call `khronon'. The dynamics of the khronon field is sensitive to the symmetries and other details of the particular implementations of the proposal. In this paper we examine several consistency issues present in three non-relativistic gravity theories: Horava's projectable theory, the healthy non-projectable extension, and a new extension related to ghost condensation. We find that the only model which is free from instabilities and strong coupling is the non-projectable one. We elaborate on the phenomenology of the latter model including a discussion of the couplings of the khronon to matter. In particular, we obtain the parameters of the post-Newtonian expansion in this model and show that they are compatible with current observations.

We address the consistency of Horava's proposal for a theory of quantum gravity from the low-energy perspective. We uncover the additional scalar degree of freedom arising from the explicit breaking of the general covariance and study its properties. The analysis is performed both in the original formulation of the theory and in the Stueckelberg picture. A peculiarity of the new mode is that it satisfies an equation of motion that is of first order in time derivatives. At linear level the mode is manifest only around spatially inhomogeneous and time-dependent backgrounds. We find two serious problems associated with this mode. First, the mode develops very fast exponential instabilities at short distances. Second, it becomes strongly coupled at an extremely low cutoff scale. We also discuss the "projectable" version of Horava's proposal and argue that this version can be understood as a certain limit of the ghost condensate model. The theory is still problematic since the additional field generically forms caustics and, again, has a very low strong coupling scale. We clarify some subtleties that arise in the application of the Stueckelberg formalism to Horava's model due to its non-relativistic nature.

We propose a natural extension of Horava's model for quantum gravity, which is free from the notorious pathologies of the original proposal. The new model endows the scalar graviton mode with a regular quadratic action and remains power-counting renormalizable. At low energies, it reduces to a Lorentz-violating scalar-tensor gravity theory. The deviations with respect to general relativity can be made weak by an appropriate choice of parameters.

A systematic study of the different phases of Lorentz-breaking massive gravity in a curved background is performed. For tensor and vector modes, the analysis is very close to that of Minkowski space. The most interesting results are in the scalar sector where, generically, there are two propagating degrees of freedom (DOF). While in maximally symmetric spaces ghost-like instabilities are inevitable, they can be avoided in a FRW background. The phases with less than two DOF in the scalar sector are also studied. Curvature allows an interesting interplay with the mass parameters; in particular, we have extended the Higuchi bound of dS to FRW and Lorentz breaking masses. As in dS, when the bound is saturated there is no propagating DOF in the scalar sector. In a number of phases the smallness of the kinetic terms gives rise to strongly coupled scalar modes at low energies. Finally, we have computed the gravitational potentials for point-like sources. In the general case we recover the GR predictions at small distances, whereas the modifications appear at distances of the order of the characteristic mass scale. In contrast with Minkowski space, these corrections may not spoil the linear approximation at large distances.