The late time acceleration of the Universe can be characterized in terms of an extra, time dependent, component of the universe -- dark energy. The simplest proposal for dark energy is a scalar-tensor theory -- quintessence -- which consists of a scalar field, $\phi$, whose dynamics is solely dictated by its potential, $V(\phi)$. Such a theory can be uniquely characterized by the equation of state of the scalar field energy momentum-tensor. We find the time dependence of the equation of state for a broad family of potentials and, using this information, we propose an analytic prior distribution for the most commonly used parametrization. We show that this analytic prior can be used to accurately predict the distribution of observables for the next generation of cosmological surveys. Including the theoretical priors in the comparison with observations considerably improves the constraints on the equation of state.

Cosmological datasets have great potential to elucidate the nature of dark energy and test gravity on the largest scales available to observation. Theoretical predictions can be computed with hi_class (www.hiclass-code.net), an accurate, fast and flexible code for linear cosmology, incorporating a wide range of dark energy theories and modifications to general relativity. We introduce three new functionalities into hi_class: (1) Support for models based on covariant Lagrangians, including a constraint-preserving integration scheme for the background evolution and a series of worked-out examples: Galileon, nKGB, quintessence (monomial, tracker) and Brans-Dicke. (2) Consistent initial conditions for the scalar-field perturbations in the deep radiation era, identifying the conditions under which modified-gravity isocurvature perturbations may grow faster than adiabatic modes leading to a loss of predictivity. (3) An automated quasistatic approximation scheme allowing order-of-magnitude improvement in computing performance without sacrificing accuracy for wide classes of models. These enhancements bring the treatment of dark energy and modified gravity models to the level of detail comparable to software tools restricted to standard $\Lambda$CDM cosmologies. The hi_class code is publicly available (https://github.com/miguelzuma/hi_class_public), ready to explore current data and prepare for next-generation experiments.

3D data compression techniques can be used to determine the natural basis of radial eigenmodes that encode the maximum amount of information in a tomographic large-scale structure survey. We explore the potential of the Karhunen-Lo\`eve decomposition in reducing the dimensionality of the data vector for cosmic shear measurements, and apply it to the final data from the \cfh survey. We find that practically all of the cosmological information can be encoded in one single radial eigenmode, from which we are able to reproduce compatible constraints with those found in the fiducial tomographic analysis (done with 7 redshift bins) with a factor of ~30 fewer datapoints. This simplifies the problem of computing the two-point function covariance matrix from mock catalogues by the same factor, or by a factor of ~800 for an analytical covariance. The resulting set of radial eigenfunctions is close to ell-independent, and therefore they can be used as redshift-dependent galaxy weights. This simplifies the application of the Karhunen-Lo\`eve decomposition to real-space and Fourier-space data, and allows one to explore the effective radial window function of the principal eigenmodes as well as the associated shear maps in order to identify potential systematics. We also apply the method to extended parameter spaces and verify that additional information may be gained by including a second mode to break parameter degeneracies. The data and analysis code are publicly available at https://github.com/emiliobellini/kl_sample.

The disconnected part of the power spectrum covariance matrix (also known as the "Gaussian" covariance) is the dominant contribution on large scales for galaxy clustering and weak lensing datasets. The presence of a complicated sky mask causes non-trivial correlations between different Fourier/harmonic modes, which must be accurately characterized in order to obtain reliable cosmological constraints. This is particularly relevant for galaxy survey data. Unfortunately, an exact calculation of these correlations involves $O(\ell_{\rm max}^6)$ operations that become computationally impractical very quickly. We present an implementation of approximate methods to estimate the Gaussian covariance matrix of power spectra involving spin-0 and spin-2 flat- and curved-sky fields, expanding on existing algorithms. These methods achieve an $O(\ell_{\rm max}^3)$ scaling, which makes the computation of the covariance matrix as fast as the computation of the power spectrum itself. We quantify the accuracy of these methods on large-scale structure and weak lensing data, making use of a large number of Gaussian but otherwise realistic simulations. We show that, using the approximate covariance matrix, we are able to recover the true posterior distribution of cosmological parameters to high accuracy. We also quantify the shortcomings of these methods, which become unreliable on the very largest scales, as well as for covariance matrix elements involving cosmic shear $B$ modes. The algorithms presented here are implemented in the public code NaMaster (https://github.com/LSSTDESC/NaMaster).

We study the phenomenology of the beyond Horndeski class of scalar-tensor theories of gravity, which on cosmological scales can be characterised in terms of one extra function of time, $\alpha_{\rm H}$, as well as the usual four Horndeski set of free functions. We show that $\alpha_{\rm H}$ can be directly related to the the damping of the matter power spectrum on both large and small scales. We also find that the temperature power spectrum of the cosmic microwave background (CMB) is enhanced at low multipoles and the lensing potential is decreased, as a function of $\alpha_{\rm H}$. We find constraints on $\alpha_{\rm H}$ of order ${\cal O}(1)$ using measurements of the temperature and polarisation of the CMB, as well as the lensing potential derived from it, combined with large scale structure data. We find that redshift space distortion measurements can play a significant role in constraining these theories. Finally, we comment on the recent constraints from the observation of an electromagnetic counterpart to a gravitational wave signal; we find that these constraints reduce the number of free parameters of the model but do not significantly change the constraints on the remaining parameters.

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.

We analyse cosmological perturbations around a homogeneous and isotropic background for scalar-tensor, vector-tensor and bimetric theories of gravity. Building on previous results, we propose a unified view of the effective parameters of all these theories. Based on this structure, we explore the viable space of parameters for each family of models by imposing the absence of ghosts and gradient instabilities. We then focus on the quasistatic regime and confirm that all these theories can be approximated by the phenomenological two-parameter model described by an effective Newton's constant and the gravitational slip. Within the quasistatic regime we pinpoint signatures which can distinguish between the broad classes of models (scalar-tensor, vector-tensor or bimetric). Finally, we present the equations of motion for our unified approach in such a way that they can be implemented in Einstein-Boltzmann solvers.

The detection of an electromagnetic counterpart (GRB 170817A) to the gravitational wave signal (GW170817) from the merger of two neutron stars opens a completely new arena for testing theories of gravity. We show that this measurement allows us to place stringent constraints on general scalar-tensor and vector-tensor theories, while allowing us to place an independent bound on the graviton mass in bimetric theories of gravity. These constraints severely reduce the viable range of cosmological models that have been proposed as alternatives to general relativistic cosmology.

We present the public version of hi_class (www.hiclass-code.net), an extension of the Boltzmann code CLASS to a broad ensemble of modifications to general relativity. In particular, hi_class can calculate predictions for models based on Horndeski's theory, which is the most general scalar-tensor theory described by second-order equations of motion and encompasses any perfect-fluid dark energy, quintessence, Brans-Dicke, $f(R)$ and covariant Galileon models. hi_class has been thoroughly tested and can be readily used to understand the impact of alternative theories of gravity on linear structure formation as well as for cosmological parameter extraction.

Cosmological observables show a dependence with the neutrino mass, which is partially degenerate with parameters of extended models of gravity. We study and explore this degeneracy in Horndeski generalized scalar-tensor theories of gravity. Using forecasted cosmic microwave background and galaxy power spectrum datasets, we find that a single parameter in the linear regime of the effective theory dominates the correlation with the total neutrino mass. For any given mass, a particular value of this parameter approximately cancels the power suppression due to the neutrino mass at a given redshift. The extent of the cancellation of this degeneracy depends on the cosmological large-scale structure data used at different redshifts. We constrain the parameters and functions of the effective gravity theory and determine the influence of gravity on the determination of the neutrino mass from present and future surveys.

The next generation of surveys will greatly improve our knowledge of cosmological gravity. In this paper we focus on how Stage IV photometric redshift surveys, including weak lensing and multiple tracers of the matter distribution and radio experiments combined with measurements of the cosmic microwave background will lead to precision constraints on deviations from General Relativity. We use a broad subclass of Horndeski scalar-tensor theories to forecast the accuracy with which we will be able to determine these deviations and their degeneracies with other cosmological parameters. Our analysis includes relativistic effects, does not rely on the quasi-static evolution and makes conservative assumptions about the effect of screening on small scales. We define a figure of merit for cosmological tests of gravity and show how the combination of different types of surveys, probing different length scales and redshifts, can be used to pin down constraints on the gravitational physics to better than a few percent, roughly an order of magnitude better than present probes. Future cosmological experiments will be able to constrain the Brans-Dicke parameter at a level comparable to Solar System and astrophysical tests

We investigate our knowledge of early universe cosmology by exploring how much additional energy density can be placed in different components beyond those in the $\Lambda$CDM model. To do this we use a method to separate early- and late-universe information enclosed in observational data, thus markedly reducing the model-dependency of the conclusions. We find that the 95\% credibility regions for extra energy components of the early universe at recombination are: non-accelerating additional fluid density parameter $\Omega_{\rm MR} < 0.006$ and extra radiation parameterised as extra effective neutrino species $2.3 < N_{\rm eff} < 3.2$ when imposing flatness. Our constraints thus show that even when analyzing the data in this largely model-independent way, the possibility of hiding extra energy components beyond $\Lambda$CDM in the early universe is seriously constrained by current observations. We also find that the standard ruler, the sound horizon at radiation drag, can be well determined in a way that does not depend on late-time Universe assumptions, but depends strongly on early-time physics and in particular on additional components that behave like radiation. We find that the standard ruler length determined in this way is $r_{\rm s} = 147.4 \pm 0.7$ Mpc if the radiation and neutrino components are standard, but the uncertainty increases by an order of magnitude when non-standard dark radiation components are allowed, to $r_{\rm s} = 150 \pm 5$ Mpc.

Recent anomalies found in cosmological datasets such as the low multipoles of the Cosmic Microwave Background or the low redshift amplitude and growth of clustering measured by e.g., abundance of galaxy clusters and redshift space distortions in galaxy surveys, have motivated explorations of models beyond standard $\Lambda$CDM. Of particular interest are models where general relativity (GR) is modified on large cosmological scales. Here we consider deviations from $\Lambda$CDM+GR within the context of Horndeski gravity, which is the most general theory of gravity with second derivatives in the equations of motion. We adopt a parametrization in which the four additional Horndeski functions of time $\alpha_i(t)$ are proportional to the cosmological density of dark energy $\Omega_{DE}(t)$. Constraints on this extended parameter space using a suite of state-of-the art cosmological observations are presented for the first time. Although the theory is able to accommodate the low multipoles of the Cosmic Microwave Background and the low amplitude of fluctuations from redshift space distortions, we find no significant tension with $\Lambda$CDM+GR when performing a global fit to recent cosmological data and thus there is no evidence against $\Lambda$CDM+GR from an analysis of the value of the Bayesian evidence ratio of the modified gravity models with respect to $\Lambda$CDM, despite introducing extra parameters. The posterior distribution of these extra parameters that we derive return strong constraints on any possible deviations from $\Lambda$CDM+GR in the context of Horndeski gravity. We illustrate how our results can be applied to a more general frameworks of modified gravity models.

Despite its continued observational successes, there is a persistent (and growing) interest in extending cosmology beyond the standard model, $\Lambda$CDM. This is motivated by a range of apparently serious theoretical issues, involving such questions as the cosmological constant problem, the particle nature of dark matter, the validity of general relativity on large scales, the existence of anomalies in the CMB and on small scales, and the predictivity and testability of the inflationary paradigm. In this paper, we summarize the current status of $\Lambda$CDM as a physical theory, and review investigations into possible alternatives along a number of different lines, with a particular focus on highlighting the most promising directions. While the fundamental problems are proving reluctant to yield, the study of alternative cosmologies has led to considerable progress, with much more to come if hopes about forthcoming high-precision observations and new theoretical ideas are fulfilled.

We investigate the limits of applicability of the quasi-static approximation in cosmologies featuring general models of dark energy or modified gravity. We show that, at best, the quasi-static approximation breaks down outside of the sound horizon of the dark-energy, rather than the cosmological horizon as is frequently assumed. When the sound speed of dark energy is significantly below that of light, the quasi-static limit is only valid in a limited range of observable scales and this must be taken into account when computing effects on observations in such models. As an order of magnitude estimate, in the analysis of data from today's weak-lensing and peculiar-velocity surveys, dark energy can be modelled as quasi-static only if the sound speed is larger than order 1% of that of light. In upcoming surveys, such as Euclid, it should only be used when the sound speed exceeds around 10% of the speed of light. In the analysis of the cosmic microwave background, the quasi-static limit should never be used for the integrated Sachs-Wolf effect and for lensing only when the sound speed exceeds 10% of the speed of light.

The scale of Baryon Acoustic Oscillations (BAO) imprinted in the matter power spectrum provides an almost-perfect standard ruler: it only suffers sub-percent deviations from fixed comoving length due to non-linear effects. We study the BAO shift in the large Horndeski class of gravitational theories and compute its magnitude in momentum space using second-order perturbation theory and a peak-background split. The standard prediction is affected by the modified linear growth, as well as by non-linear gravitational effects that alter the mode-coupling kernel. For covariant Galileon models, we find a $14-45\%$ enhancement of the BAO shift with respect to standard gravity and a distinct time evolution depending on the parameters. Despite the larger values, the shift remains well below the forecasted precision of next-generation galaxy surveys. Models that produce significant BAO shift would cause large redshift-space distortions or affect the bispectrum considerably. Our computation therefore validates the use of the BAO scale as a comoving standard ruler for tests of general dark energy models.

We present a detailed study of second-order matter perturbations for the general Horn- deski class of models. Being the most general scalar-tensor theory having second-order equations of motion, it includes many known gravity and dark energy theories and General Relativity with a cosmological constant as a specific case. This enables us to estimate the leading order dark matter bispectrum generated at late-times by gravitational instability. We parametrize the evolution of the first and second-order equations of motion as proposed by Bellini and Sawicki (2014), where the free functions of the theory are assumed to be proportional to the dark energy density. We show that it is unnatural to have large 10% ( 1%) deviations of the bispectrum introducing even larger ~ 30% (~ 5%) deviations in the linear growth rate. Considering that measurements of the linear growth rate have much higher signal-to-noise than bispectrum measurements, this indicates that for Horndeski models which reproduce the expansion history and the linear growth rate as predicted by GR the dark matter bispectrum kernel can be effectively modelled as the standard GR one. On the other hand, an observation of a large bispectrum deviation that can not be explained in terms of bias would imply either that the evolution of perturbations is strongly different than the evolution predicted by GR or that the theory of gravity is exotic (e.g., breaks the weak equivalence principle) and/or fine-tuned.

The existence of a cosmic neutrino background can be probed indirectly by CMB experiments, not only by measuring the background density of radiation in the universe, but also by searching for the typical signatures of the fluctuations of free-streaming species in the temperature and polarisation power spectrum. Previous studies have already proposed a rather generic parametrisation of these fluctuations, that could help to discriminate between the signature of ordinary free-streaming neutrinos, or of more exotic dark radiation models. Current data are compatible with standard values of these parameters, which seems to bring further evidence for the existence of a cosmic neutrino background. In this work, we investigate the robustness of this conclusion under various assumptions. We generalise the definition of an effective sound speed and viscosity speed to the case of massive neutrinos or other dark radiation components experiencing a non-relativistic transition. We show that current bounds on these effective parameters do not vary significantly when considering an arbitrary value of the particle mass, or extended cosmological models with a free effective neutrino number, dynamical dark energy or a running of the primordial spectrum tilt. We conclude that it is possible to make a robust statement about the detection of the cosmic neutrino background by CMB experiments.

We present a turnkey solution, ready for implementation in numerical codes, for the study of linear structure formation in general scalar-tensor models involving a single universally coupled scalar field. We show that the totality of cosmological information on the gravitational sector can be compressed - without any redundancy - into five independent and arbitrary functions of time only and one constant. These describe physical properties of the universe: the observable background expansion history, fractional matter density today, and four functions of time describing the properties of the dark energy. We show that two of those dark-energy property functions control the existence of anisotropic stress, the other two - dark-energy clustering, both of which are can be scale-dependent. All these properties can in principle be measured, but no information on the underlying theory of acceleration beyond this can be obtained. We present a translation between popular models of late-time acceleration (e.g. perfect fluids, f (R), kinetic gravity braiding, galileons), as well as the effective field theory framework, and our formulation. In this way, implementing this formulation numerically would give a single tool which could consistently test the majority of models of late-time acceleration heretofore proposed.

We use recent measurements of the expansion history of the universe to place constraints on the parameter space of cubic Galileon models, in particular we concentrate on those models which contain the simplest Galileon term plus a linear potential. This gives strong constraints on the Lagrangian of these models. Most dynamical terms in the Galileon Lagrangian are constraint to be small and the acceleration is effectively provided by a constant term in the scalar potential, thus reducing, effectively, to a LCDM model for current acceleration. The effective equation of state is indistinguishable from that of a cosmological constant w = -1 and the data constraint it to have no temporal variations of more than at the few % level. The energy density of the Galileon can contribute only to about 10% of the acceleration energy density, being the other 90% a cosmological constant term. This demonstrates how useful direct measurements of the expansion history of the universe are at constraining the dynamical nature of dark energy.

In this paper we obtain the bispectrum of dark matter density perturbations in the frame of covariant cubic Galileon theories. This result is obtained by means of a semi-analytic approach to second-order perturbations in Galileon cosmologies, assuming Gaussian initial conditions. In particular, we show that, even in the presence of large deviations of the linear growth-rate w.r.t. the {\Lambda}CDM one, at the bispectrum level such deviations are reduced to a few percent.

In this paper we study the evolution of a spherical matter overdensity in the context of the recently introduced Galileon field theory. Our analysis considers the complete covariant Lagrangian in four dimensions. This theory is composed by a potential and a standard kinetic term, a cubic kinetic term and two additional terms that include the coupling between the Galileon and the metric, to preserve the original properties of Galileons also in curved space-times. Here we extend previous studies, which considered both the quintessence and the cubic terms, by focussing on the role of the last two terms. The background evolution we consider is driven by a tracker solution. Studying scalar perturbations in the non-linear regime, we find constraints on the parameter of the model. We will show how the new terms contribute to the collapse phase and how they modify physical parameters, such as the linearized density contrast and the virial overdensity. The results show that the Galileon modifies substantially the dynamics of the collapse, thus making it possible to observationally constrain the parameters of this theory.