Late times dark energy transitions at redshifts $z \ll 0.1$ can raise the predicted value of the Hubble constant to the SH0ES value, $74.03\pm 1.42$ (km s$^{-1}$ Mpc$^{-1})$ or more, while providing an equally good fit as $\Lambda$CDM at $67.73 \pm 0.41$ to higher redshift data, in particular from the cosmic microwave background and baryon acoustic oscillations. These models however do not fully resolve the true source of tension between the distance ladder and high redshift observations: the local calibration of supernovae luminosities well out into the Hubble flow. When tested in this manner by transferring the SH0ES calibration to the Pantheon supernovae dataset, the ability of such transitions to raise the Hubble constant is reduced to $69.17 \pm 1.09$. Such an analysis should also be used when testing any dynamical dark energy model which can produce similarly fine features in redshift or local void models.

We apply a tension metric $Q_\textrm{UDM}$, the update difference in mean parameters, to understand the source of the difference in the measured Hubble constant $H_0$ inferred with cosmic microwave background lensing measurements from the Planck satellite ($H_0=67.9^{+1.1}_{-1.3}\, \mathrm{km/s/Mpc}$) and from the South Pole Telescope ($H_0=72.0^{+2.1}_{-2.5}\, \mathrm{km/s/Mpc}$) when both are combined with baryon acoustic oscillation (BAO) measurements with priors on the baryon density (BBN). $Q_\textrm{UDM}$ isolates the relevant parameter directions for tension or concordance where the two data sets are both informative, and aids in the identification of subsets of data that source the observed tension. With $Q_\textrm{UDM}$, we uncover that the difference in $H_0$ is driven by the tension between Planck lensing and BAO+BBN, at probability-to-exceed of 6.6%. Most of this mild tension comes from the galaxy BAO measurements parallel to the line of sight. The redshift dependence of the parallel BAOs pulls both the matter density $\Omega_m$ and $H_0$ high in $\Lambda$CDM, but these parameter anomalies are usually hidden when the BAO measurements are combined with other cosmological data sets with much stronger $\Omega_m$ constraints.

Cosmological analyses of galaxy surveys rely on knowledge of the redshift distribution of their galaxy sample. This is usually derived from a spectroscopic and/or many-band photometric calibrator survey of a small patch of sky. The uncertainties in the redshift distribution of the calibrator sample include a contribution from shot noise, or Poisson sampling errors, but, given the small volume they probe, they are dominated by sample variance introduced by large-scale structures. Redshift uncertainties have been shown to constitute one of the leading contributions to systematic uncertainties in cosmological inferences from weak lensing and galaxy clustering, and hence they must be propagated through the analyses. In this work, we study the effects of sample variance on small-area redshift surveys, from theory to simulations to the COSMOS2015 data set. We present a three-step Dirichlet method of resampling a given survey-based redshift calibration distribution to enable the propagation of both shot noise and sample variance uncertainties. The method can accommodate different levels of prior confidence on different redshift sources. This method can be applied to any calibration sample with known redshifts and phenotypes (i.e. cells in a self-organizing map, or some other way of discretizing photometric space), and provides a simple way of propagating prior redshift uncertainties into cosmological analyses. As a worked example, we apply the full scheme to the COSMOS2015 data set, for which we also present a new, principled SOM algorithm designed to handle noisy photometric data. We make available a catalog of the resulting resamplings of the COSMOS2015 galaxies.

We use supernovae measurements, calibrated by the local determination of the Hubble constant $H_0$ by SH0ES, to interpolate the distance-redshift relation using Gaussian process regression. We then predict, independent of the cosmological model, the distances that are measured with strong lensing time delays. We find excellent agreement between these predictions and the measurements. The agreement holds when we consider only the redshift dependence of the distance-redshift relation, independent of the value of $H_0$. Our results disfavor the possibility that lens mass modeling contributes a 10\% bias or uncertainty in the strong lensing analysis, as suggested recently in the literature. In general our analysis strengthens the case that residual systematic errors in both measurements are below the level of the current discrepancy with the CMB determination of $H_0$, and supports the possibility of new physical phenomena on cosmological scales. With additional data our methodology can provide more stringent tests of unaccounted for systematics in the determinations of the distance-redshift relation in the late universe.

We develop estimators of agreement and disagreement between correlated cosmological data sets. These account for data correlations when computing the significance of both tensions and excess confirmation while remaining statistically optimal. We discuss and thoroughly characterize different approaches commenting on the ones that have the best behavior in practical applications. We complement the calculation of their statistical distribution within the Gaussian model with one estimator that takes non-Gaussianities fully into account. To illustrate the use of our techniques, we apply these estimators to supernovae measurements of the distance-redshift relation, absolutely calibrated by the local distance ladder. The suite of best estimators that we discuss finds results that are in excellent agreement between estimators and find no indications of significant internal inconsistencies in this data set above the $1\%$ probability threshold. This shows the robustness of local determinations of the Hubble constant to features in the distance-redshift relation.

We present constraints on extensions of the minimal cosmological models dominated by dark matter and dark energy, $\Lambda$CDM and $w$CDM, by using a combined analysis of galaxy clustering and weak gravitational lensing from the first-year data of the Dark Energy Survey (DES Y1) in combination with external data. We consider four extensions of the minimal dark energy-dominated scenarios: 1) nonzero curvature $\Omega_k$, 2) number of relativistic species $N_{\rm eff}$ different from the standard value of 3.046, 3) time-varying equation-of-state of dark energy described by the parameters $w_0$ and $w_a$ (alternatively quoted by the values at the pivot redshift, $w_p$, and $w_a$), and 4) modified gravity described by the parameters $\mu_0$ and $\Sigma_0$ that modify the metric potentials. We also consider external information from Planck CMB measurements; BAO measurements from SDSS, 6dF, and BOSS; RSD measurements from BOSS; and SNIa information from the Pantheon compilation. Constraints on curvature and the number of relativistic species are dominated by the external data; when these are combined with DES Y1, we find $\Omega_k=0.0020^{+0.0037}_{-0.0032}$ at the 68% confidence level, and $N_{\rm eff}<3.28\, (3.55)$ at 68% (95%) confidence. For the time-varying equation-of-state, we find the pivot value $(w_p, w_a)=(-0.91^{+0.19}_{-0.23}, -0.57^{+0.93}_{-1.11})$ at pivot redshift $z_p=0.27$ from DES alone, and $(w_p, w_a)=(-1.01^{+0.04}_{-0.04}, -0.28^{+0.37}_{-0.48})$ at $z_p=0.20$ from DES Y1 combined with external data; in either case we find no evidence for the temporal variation of the equation of state. For modified gravity, we find the present-day value of the relevant parameters to be $\Sigma_0= 0.43^{+0.28}_{-0.29}$ from DES Y1 alone, and $(\Sigma_0, \mu_0)=(0.06^{+0.08}_{-0.07}, -0.11^{+0.42}_{-0.46})$ from DES Y1 combined with external data, consistent with predictions from GR.

We discuss the ability of a dark fluid becoming relevant around the time of matter radiation equality to significantly relieve the tension between local measurements of the Hubble constant and CMB inference, within the $\Lambda$CDM model. We show that the gravitational impact of acoustic oscillations in the dark fluid balance the effects on the CMB and result in an improved fit to CMB measurements themselves while simultaneously raising the Hubble constant. The required balance favors a model where the fluid is a scalar field that converts its potential to kinetic energy around matter radiation equality which then quickly redshifts away. We derive the requirements on the potential for this conversion mechanism and find that a simple canonical scalar with two free parameters for its local slope and amplitude robustly improves the fit to the combined data by $\Delta\chi^2 \approx 12.7$ over $\Lambda$CDM. We uncover the CMB polarization signatures that can definitively test this scenario with future data.

We provide an overview of the science case, instrument configuration and project plan for the next-generation ground-based cosmic microwave background experiment CMB-S4, for consideration by the 2020 Decadal Survey.

We present the science case, reference design, and project plan for the Stage-4 ground-based cosmic microwave background experiment CMB-S4.

While Type Ia Supernovae (SNe Ia) are one of the most mature cosmological probes, the next era promises to be extremely exciting in the number of different ways SNe Ia are used to measure various cosmological parameters. Here we review the experiments in the 2020s that will yield orders of magnitudes more SNe Ia, and the new understandings and capabilities to constrain systematic uncertainties at a level to match these statistics. We then discuss five different cosmological probes with SNe Ia: the conventional Hubble diagram for measuring dark energy properties, the distance ladder for measuring the Hubble constant, peculiar velocities and weak lensing for measuring sigma8 and strong-lens measurements of H0 and other cosmological parameters. For each of these probes, we discuss the experiments that will provide the best measurements and also the SN Ia-related systematics that affect each one.

Our current understanding of the Universe is established through the pristine measurements of structure in the cosmic microwave background (CMB) and the distribution and shapes of galaxies tracing the large scale structure (LSS) of the Universe. One key ingredient that underlies cosmological observables is that the field that sources the observed structure is assumed to be initially Gaussian with high precision. Nevertheless, a minimal deviation from Gaussianityis perhaps the most robust theoretical prediction of models that explain the observed Universe; itis necessarily present even in the simplest scenarios. In addition, most inflationary models produce far higher levels of non-Gaussianity. Since non-Gaussianity directly probes the dynamics in the early Universe, a detection would present a monumental discovery in cosmology, providing clues about physics at energy scales as high as the GUT scale.

The hot dense environment of the early universe is known to have produced large numbers of baryons, photons, and neutrinos. These extreme conditions may have also produced other long-lived species, including new light particles (such as axions or sterile neutrinos) or gravitational waves. The gravitational effects of any such light relics can be observed through their unique imprint in the cosmic microwave background (CMB), the large-scale structure, and the primordial light element abundances, and are important in determining the initial conditions of the universe. We argue that future cosmological observations, in particular improved maps of the CMB on small angular scales, can be orders of magnitude more sensitive for probing the thermal history of the early universe than current experiments. These observations offer a unique and broad discovery space for new physics in the dark sector and beyond, even when its effects would not be visible in terrestrial experiments or in astrophysical environments. A detection of an excess light relic abundance would be a clear indication of new physics and would provide the first direct information about the universe between the times of reheating and neutrino decoupling one second later.

We present the data-driven reconstruction of gravitational theories and Dark Energy models on cosmological scales. We showcase the power of present cosmological probes at constraining these models and quantify the knowledge of their properties that can be acquired through state of the art data. This reconstruction exploits the power of the Effective Field Theory approach to Dark Energy and Modified Gravity phenomenology, which compresses the freedom in defining such models into a finite set of functions that can be reconstructed across cosmic times using cosmological data. We consider several model classes described within this framework and thoroughly discuss their phenomenology and data implications. We find that some models can alleviate the present discrepancy in the determination of the Hubble constant as inferred from the cosmic microwave background and as directly measured. This results in a statistically significant preference for the reconstructed theories over the standard cosmological model.

Ongoing and upcoming cosmological surveys will significantly improve our ability to probe the equation of state of dark energy, $w_{\rm DE}$, and the phenomenology of Large Scale Structure. They will allow us to constrain deviations from the $\Lambda$CDM predictions for the relations between the matter density contrast and the weak lensing and the Newtonian potential, described by the functions $\Sigma$ and $\mu$, respectively. In this work, we derive the theoretical prior for the joint covariance of $w_{\rm DE}$, $\Sigma$ and $\mu$, expected in general scalar-tensor theories with second order equations of motion (Horndeski gravity), focusing on their time-dependence at certain representative scales. We employ Monte-Carlo methods to generate large ensembles of statistically independent Horndeski models, focusing on those that are physically viable and in broad agreement with local tests of gravity, the observed cosmic expansion history and the measurement of the speed of gravitational waves from a binary neutron star merger. We identify several interesting features and trends in the distribution functions of $w_{\rm DE}$, $\Sigma$ and $\mu$, as well as in their covariances; we confirm the high degree of correlation between $\Sigma$ and $\mu$ in scalar-tensor theories. The derived prior covariance matrices will allow us to reconstruct jointly $w_{\rm DE}$, $\Sigma$ and $\mu$ in a non-parametric way.

We discuss the cosmological implications of the string swampland conjectures for late-time cosmology, and test them against a wide range of state of the art cosmological observations. The refined de Sitter conjecture constrains either the minimal slope or the curvature of the scalar potential, and depends on two dimensionless constants. For constants of size one or larger, tension exists between observations, especially the Hubble constant, and the slope and curvature conjectures at a level of 4.5 sigma and 2.3 sigma, respectively. Smaller values of the constants are permitted by observations, and we determine upper bounds at varying confidence levels. We also derive and constrain the relationship between cosmological observables and the scalar field excursion during the acceleration epoch, thereby testing the distance conjecture.

We discuss the phenomenological imprints of modifications to gravity in the early universe with a specific focus on the time of recombination. We derive several interesting results regarding the effect that such modifications have on cosmological observables, especially on the driving and phasing of acoustic oscillations, observed in the CMB and BAO, as well as the weak gravitational lensing of the CMB and of galaxy shapes. This widens the pool of measurements that can be used to test gravity with present and future surveys, in particular realizing the full constraining power of the structure of the primary peaks of the CMB spectrum. We investigate whether such a phenomenology can relax tensions between cosmological measurements and find that a modification of the gravitational constant at recombination would help in reconciling measurements of the CMB with local measurements of the Hubble constant.

We investigate cosmological constraints on K-mouflage models of modified gravity. We consider two scenarios: one where the background evolution is free to deviate from $\Lambda$CDM (K-mouflage) and another one which reproduces a $\Lambda$CDM expansion (K-mimic), implementing both of them into the EFTCAMB code. We discuss the main observational signatures of these models and we compare their cosmological predictions to different datasets, including CMB, CMB lensing, SNIa and different galaxy catalogues. We argue about the possibility of relieving the $H_0$ and weak lensing tensions within these models, finding that K-mouflage scenarios effectively ease the tension on the Hubble Constant. Our final 95\% C.L. bounds on the $\epsilon_{2,0}$ parameter that measures the overall departure from $\Lambda$CDM (corresponding to $\epsilon_{2,0}=0$) are $-0.04\leq \epsilon_{2,0} <0$ for K-mouflage and $0< \epsilon_{2,0} <0.002$ for K-mimic. In the former case the main constraining power comes from changes in the background expansion history, while in the latter case the model is strongly constrained by measurements of the amplitude of matter perturbations. The sensitivity of these cosmological constraints closely matches that of solar system probes. We show that these constraints could be significantly tightened with future ideal probes like CORE.

The success of present and future cosmological studies is tied to the ability to detect discrepancies in complex data sets within the framework of a cosmological model. Tensions caused by the presence of unknown systematic effects need to be isolated and corrected to increase the overall accuracy of parameter constraints, while discrepancies due to new physical phenomena need to be promptly identified. We develop a full set of estimators of internal and mutual agreement and disagreement, whose strengths complement each other. These allow to take into account the effect of prior information and compute the statistical significance of both tensions and confirmatory biases. We apply them to a wide range of state of the art cosmological probes and show that these estimators can be easily used, regardless of model and data complexity. We derive a series of results that show that discrepancies indeed arise within the standard LCDM model. Several of them exceed the probability threshold of 95% and deserve a dedicated effort to understand their origin.

Cosmological measurements of structure are placing increasingly strong constraints on the sum of the neutrino masses, $\Sigma m_\nu$, through Bayesian inference. Because these constraints depend on the choice for the prior probability $\pi(\Sigma m_\nu)$, we argue that this prior should be motivated by fundamental physical principles rather than the ad hoc choices that are common in the literature. The first step in this direction is to specify the prior directly at the level of the neutrino mass matrix $M_\nu$, since this is the parameter appearing in the Lagrangian of the particle physics theory. Thus by specifying a probability distribution over $M_\nu$, and by including the known squared mass splittings, we predict a theoretical probability distribution over $\Sigma m_\nu$ that we interpret as a Bayesian prior probability $\pi(\Sigma m_\nu)$. We find that $\pi(\Sigma m_\nu)$ peaks close to the smallest $\Sigma m_\nu$ allowed by the measured mass splittings, roughly $0.06 \, {\rm eV}$ ($0.1 \, {\rm eV}$) for normal (inverted) ordering, due to the phenomenon of eigenvalue repulsion in random matrices. We consider three models for neutrino mass generation: Dirac, Majorana, and Majorana via the seesaw mechanism; differences in the predicted priors $\pi(\Sigma m_\nu)$ allow for the possibility of having indications about the physical origin of neutrino masses once sufficient experimental sensitivity is achieved. We present fitting functions for $\pi(\Sigma m_\nu)$, which provide a simple means for applying these priors to cosmological constraints on the neutrino masses or marginalizing over their impact on other cosmological parameters.

We revisit the cosmology of Covariant Galileon gravity in view of the most recent cosmological data sets, including weak lensing. As a higher derivative theory, Covariant Galileon models do not have a $\Lambda$CDM limit and predict a very different structure formation pattern compared with the standard $\Lambda$CDM scenario. Previous cosmological analyses suggest that this model is marginally disfavoured, yet can not be completely ruled out. In this work we use a more recent and extended combination of data, and we allow for more freedom in the cosmology, by including a massive neutrino sector with three different mass hierarchies. We use the Planck measurements of Cosmic Microwave Background temperature and polarization; Baryonic Acoustic Oscillations measurements by BOSS DR12; local measurements of $H_0$; the joint light-curve analysis supernovae sample; and, for the first time, weak gravitational lensing from the KiDS collaboration. We find, that in order to provide a reasonable fit, a non-zero neutrino mass is indeed necessary, but we do not report any sizable difference among the three neutrino hierarchies. Finally, the comparison of the Bayesian Evidence to the $\Lambda$CDM one shows that in all the cases considered, Covariant Galileon models are statistically ruled out by cosmological data.

Phenomenological functions $\Sigma$ and $\mu$, also known as $G_{\rm light}/G$ and $G_{\rm matter}/G$, are commonly used to parameterize modifications of the growth of large-scale structure in alternative theories of gravity. We study the values these functions can take in Horndeski theories, i.e. the class of scalar-tensor theories with second order equations of motion. We restrict our attention to models that are in a broad agreement with tests of gravity and the observed cosmic expansion history. In particular, we require the speed of gravity to be equal to the speed of light today, as required by the recent detection of gravitational waves and electromagnetic emission from a binary neutron star merger. We examine the correlations between the values of $\Sigma$ and $\mu$ analytically within the quasi-static approximation, and numerically, by sampling the space of allowed solutions. We confirm that the conjecture made in [Pogosian:2016pwr], that $(\Sigma-1)(\mu -1) \ge 0$ in viable Horndeski theories, holds very well. Along with that, we check the validity of the quasi-static approximation within different corners of Horndeski theory. Our results show that, even with the tight bound on the present day speed of gravitational waves, there is room within Horndeski theories for non-trivial signatures of modified gravity at the level of linear perturbations.

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 investigate the impact of general conditions of theoretical stability and cosmological viability on dynamical dark energy models. As a powerful example, we study whether minimally coupled, single field Quintessence models that are safe from ghost instabilities, can source the CPL expansion history recently shown to be mildly favored by a combination of CMB (Planck) and Weak Lensing (KiDS) data. We find that in their most conservative form, the theoretical conditions impact the analysis in such a way that smooth single field Quintessence becomes significantly disfavored with respect to the standard LCDM cosmological model. This is due to the fact that these conditions cut a significant portion of the (w0;wa) parameter space for CPL, in particular eliminating the region that would be favored by weak lensing data. Within the scenario of a smooth dynamical dark energy parametrized with CPL, weak lensing data favors a region that would require multiple fields to ensure gravitational stability.

EFTCAMB/EFTCosmoMC are publicly available patches to the CAMB/CosmoMC codes implementing the effective field theory approach to single scalar field dark energy and modified gravity models. With the present numerical notes we provide a guide to the technical details of the code. Moreover we reproduce, as they appear in the code, the complete set of the modified equations and the expressions for all the other relevant quantities used to construct these patches. We submit these notes to the arXiv to grant full and permanent access to this material which provides very useful guidance to the numerical implementation of the EFT framework. We will update this set of notes when relevant modifications to the EFTCAMB/EFTCosmoMC codes will be released. The present version is based on the version of EFTCAMB/EFTCosmoMC Sep17.

Observations of the cosmic microwave background (CMB) together with weak lensing measurements of the clustering of large scale cosmological structures and local measurements of the Hubble constant pose a challenge to the standard $\Lambda$CDM cosmological model. On one side CMB observations imply a Hubble constant that is lower than local measurements and an amplitude of the lensing signal that is higher than direct measurements from weak lensing surveys. We investigate a way of relieving these tensions by adding dark radiation tightly coupled to an acoustic part of the dark matter sector and compare it to massive neutrino solutions. While these models offer a way of separately relieving the Hubble and weak lensing tensions they are prevented from fully accommodating both at the same time since the CMB requires additional cold dark matter when adding acoustic dark matter or massive neutrinos to preserve the same sharpness of the acoustic peaks which counteracts the desired growth suppression.

A flat Friedman-Roberson-Walker universe dominated by a cosmological constant ($\Lambda$) and cold dark matter (CDM) has been the working model preferred by cosmologists since the discovery of cosmic acceleration. However, tensions of various degrees of significance are known to be present among existing datasets within the $\Lambda$CDM framework. In particular, the Lyman-$\alpha$ forest measurement of the Baryon Acoustic Oscillations (BAO) by the Baryon Oscillation Spectroscopic Survey (BOSS) prefers a smaller value of the matter density fraction $\Omega_{\rm M}$ compared to the value preferred by cosmic microwave background (CMB). Also, the recently measured value of the Hubble constant, $H_0=73.24\pm1.74 \ {\rm km}\ {\rm s}^{-1} \ {\rm Mpc}^{-1}$, is $3.4\sigma$ higher than $66.93\pm0.62 \ {\rm km}\ {\rm s}^{-1} \ {\rm Mpc}^{-1}$ inferred from the Planck CMB data. In this work, we investigate if these tensions can be interpreted as evidence for a non-constant dynamical dark energy (DE). Using the Kullback-Leibler (KL) divergence to quantify the tension between datasets, we find that the tensions are relieved by an evolving DE, with the dynamical DE model preferred at a $3.5\sigma$ significance level based on the improvement in the fit alone. While, at present, the Bayesian evidence for the dynamical DE is insufficient to favour it over $\Lambda$CDM, we show that, if the current best fit DE happened to be the true model, it would be decisively detected by the upcoming DESI survey.

Constraining the Dark Energy (DE) equation of state, w, is one of the primary science goals of ongoing and future cosmological surveys. In practice, with imperfect data and incomplete redshift coverage, this requires making assumptions about the evolution of w with redshift z. These assumptions can be manifested in a choice of a specific parametric form, which can potentially bias the outcome, or else one can reconstruct w(z) non-parametrically, by specifying a prior covariance matrix that correlates values of w at different redshifts. In this work, we derive the theoretical prior covariance for the effective DE equation of state predicted by general scalar-tensor theories with second order equations of motion (Horndeski theories). This is achieved by generating a large ensemble of possible scalar-tensor theories using a Monte Carlo methodology, including the application of physical viability conditions. We also separately consider the special sub-case of the minimally coupled scalar field, or quintessence. The prior shows a preference for tracking behaviors in the most general case. Given the covariance matrix, theoretical priors on parameters of any specific parametrization of w(z) can also be readily derived by projection.

This book lays out the scientific goals to be addressed by the next-generation ground-based cosmic microwave background experiment, CMB-S4, envisioned to consist of dedicated telescopes at the South Pole, the high Chilean Atacama plateau and possibly a northern hemisphere site, all equipped with new superconducting cameras. CMB-S4 will dramatically advance cosmological studies by crossing critical thresholds in the search for the B-mode polarization signature of primordial gravitational waves, in the determination of the number and masses of the neutrinos, in the search for evidence of new light relics, in constraining the nature of dark energy, and in testing general relativity on large scales.

Clusters of galaxies have the potential of providing powerful constraints on possible deviations from General Relativity. We use the catalogue of Sunyaev-Zel'dovich sources detected by Planck and consider a correction to the halo mass function for a f(R) class of modified gravity models, which has been recently found to reproduce well results from N-body simulations, to place constraints on the scalaron field amplitude at the present time, $f_{R}^0$. We find that applying this correction to different calibrations of the halo mass function produces upper bounds on $f_{R}^0$ tighter by more than an order of magnitude, ranging from $\log_{10}(-f_{R}^0) < -5.81$ to $\log_{10}(-f_{R}^0) < -4.40$ (95 % confidence level). This sensitivity is due to the different shape of the halo mass function, which is degenerate with the parameters used to calibrate the scaling relations between SZ observables and cluster masses. Any claim of constraints more stringent that the weaker limit above, based on cluster number counts, appear to be premature and must be supported by a careful calibration of the halo mass function and by a robust calibration of the mass scaling relations.

These notes show and comment the examples that have been used to validate the CosmicFish code. We compare the results obtained with the code to several other results available in literature finding an overall good level of agreement. We will update this set of notes when relevant modifications to the CosmicFish code will be released or other validation examples are worked out. The CosmicFish code and the package to produce all the validation results presented here are publicly available at http://cosmicfish.github.io. The present version is based on CosmicFish Jun16.

CosmicFish is a publicly available library to perform Fisher matrix forecast for several cosmological observations. With the present implementation notes we provide a guide to the physical and technical details of the library. We reproduce here the details and all the relevant equations, as they appear in the code. We submit these notes to the arXiv to grant full and permanent access to this material which provides a useful guidance to forecasting and the use of CosmicFish code. We will update this set of notes when relevant modifications to the CosmicFish code will be released. The present version is based on CosmicFish Jun16.

Facing the advent of the next generation cosmological surveys we present a method to forecast knowledge gain on cosmological models. We propose this as a well defined and general tool to quantify the performance of different experiments in relation to different theoretical models. In particular, the assessment of experimental performance will benefit enormously from the fact that this method is invariant under re-parametrization of the model. We apply this to future surveys and compare expected knowledge advancements to the most relevant experiments performed over the history of modern cosmology. When considering the standard cosmological model, we show that it will rapidly reach knowledge saturation in the near future and forthcoming improvements will not match the past ones. On the contrary, we find that new observations have the potential for unprecedented knowledge jumps when extensions of the standard scenario are considered.

We consider Horava gravity within the framework of the effective field theory (EFT) of dark energy and modified gravity. We work out a complete mapping of the theory into the EFT language for an action including all the operators which are relevant for linear perturbations with up to sixth order spatial derivatives. We then employ an updated version of the EFTCAMB/EFTCosmoMC package to study the cosmology of the low-energy limit of Horava gravity and place constraints on its parameters using several cosmological data sets. In particular we use cosmic microwave background (CMB) temperature-temperature and lensing power spectra by Planck 2013, WMAP low-l polarization spectra, WiggleZ galaxy power spectrum, local Hubble measurements, Supernovae data from SNLS, SDSS and HST and the baryon acoustic oscillations measurements from BOSS, SDSS and 6dFGS. We get improved upper bounds, with respect to those from Big Bang Nucleosynthesis, on the deviation of the cosmological gravitational constant from the local Newtonian one. At the level of the background phenomenology, we find a relevant rescaling of the Hubble rate at all epoch, which has a strong impact on the cosmological observables; at the level of perturbations, we discuss in details all the relevant effects on the observables and find that in general the quasi-static approximation is not safe to describe the evolution of perturbations. Overall we find that the effects of the modifications induced by the low-energy Horava gravity action are quite dramatic and current data place tight bounds on the theory parameters.

We show how to fully map a specific model of modified gravity into the Einstein-Boltzmann solver EFTCAMB. This approach consists in few steps and allows to obtain the cosmological phenomenology of a model with minimal effort. We discuss all these steps, from the solution of the dynamical equations for the cosmological background of the model to the use of the mapping relations to cast the model into the effective field theory language and use the latter to solve for perturbations. We choose the Hu-Sawicki f(R) model of gravity as our working example. After solving the background and performing the mapping, we interface the algorithm with EFTCAMB and take advantage of the effective field theory framework to integrate the full dynamics of linear perturbations, returning all quantities needed to accurately compare the model with observations. We discuss some observational signatures of this model, focusing on the linear growth of cosmic structures. In particular we present the behavior of $f\sigma_8$ and $E_G$ that, unlike the $\Lambda$CDM scenario, are generally scale dependent in addition to redshift dependent. Finally, we study the observational implications of the model by comparing its cosmological predictions to the Planck 2015 data, including CMB lensing, the WiggleZ galaxy survey and the CFHTLenS weak lensing survey measurements. We find that while WiggleZ data favor a non-vanishing value of the Hu-Sawicki model parameter, $\log_{10}(-f^0_{R})$, and consequently a large value of $\sigma_8$, CFHTLenS drags the estimate of $\log_{10}(-f^0_{R})$ back to the $\Lambda$CDM limit.

We use a complete and rigorous statistical indicator to measure the level of concordance between cosmological data sets, without relying on the inspection of the marginal posterior distribution of some selected parameters. We apply this test to state of the art cosmological data sets, to assess their agreement within the $\Lambda$CDM model. We find that there is a good level of concordance between all the experiments with one noticeable exception. There is substantial evidence of tension between the CMB, temperature and polarization, measurements of the Planck satellite and the data from the CFHTLenS weak lensing survey even when applying ultra conservative cuts. These results robustly point toward the possibility of having unaccounted systematic effects in the data, an incomplete modelling of the cosmological predictions or hints toward new physical phenomena.

Planck-2015 data seem to favour a large value of the lensing amplitude parameter, $A_{\rm L}=1.22\pm0.10$, in CMB spectra. This result is in $2\sigma$ tension with the lensing reconstruction result, $A_{\rm L}^{\phi\phi}=0.95\pm0.04$. In this paper, we simulate several CMB anisotropy and CMB lensing spectra based on Planck-2015 best-fit cosmological parameter values and Planck blue book beam and noise specifications. We analyse several modified gravity models within the effective field theory framework against these simulations and find that models whose effective Newton constant is enhanced can modulate the CMB anisotropy spectra in a way similar to that of the $A_{\rm L}$ parameter. However, in order to lens the CMB anisotropies sufficiently, like in the Planck-2015 results, the growth of matter perturbations is substantially enhanced and gives a high $\sigma_8$ value. This in turn proves to be problematic when combining these data to other probes, like weak lensing from CFHTLenS, that favour a smaller amplitude of matter fluctuations.

We revisit the degeneracy between massive neutrinos and generalized theories of gravity in the framework of effective field theory of cosmic acceleration. In particular we consider f(R) theories and a class of non-minimally coupled models parametrized via a coupling to gravity which is linear in the scale factor. In the former case, we find a slightly lower degeneracy with respect to what found so far in the literature, due to the fact that we implement exact designer f(R) models and evolve the full linear dynamics of perturbations. As a consequence, our bounds are slightly tighter on the f(R) parameter but looser on the summed neutrino mass. We also set a new upper bound on the Compton wavelength parameter ${\rm Log}_{10}B_0<-4.1$ at 95% C.L. with fixed summed neutrino mass ($\sum m_{\nu}=0.06$ eV) in f(R) gravity with the combined data sets from cosmic microwave background temperature and lensing power spectra of Planck collaboration as well as galaxy power spectrum from WiggleZ dark energy survey. We do not observe a sizable degeneracy between massive neutrinos and modified gravity in the linear parametrization of non-minimally gravitational coupling model. The analysis is performed with an updated version of the EFTCAMB/EFTCosmoMC package, which is now publicly available and extends the first version of the code with the consistent inclusion of massive neutrinos, tensor modes, several alternative background histories and designer quintessence models.

We introduce EFTCAMB/EFTCosmoMC as publicly available patches to the commonly used CAMB/CosmoMC codes. We briefly describe the structure of the codes, their applicability and main features. To illustrate the use of these patches, we obtain constraints on parametrized pure EFT and designer $f(R)$ models, both on $\Lambda$CDM and $w$CDM background expansion histories, using data from Planck temperature and lensing potential spectra, WMAP low-$\ell$ polarization spectra (WP), and baryon acoustic oscillation (BAO). Upon inspecting theoretical stability of the models on the given background, we find non-trivial parameter spaces that we translate into viability priors. We use different combinations of data sets to show their individual effects on cosmological and model parameters. Our data analysis results show that, depending on the adopted data sets, in the $w$CDM background case this viability priors could dominate the marginalized posterior distributions. Interestingly, with Planck+WP+BAO+lensing data, in $f(R)$ gravity models, we get very strong constraints on the constant dark energy equation of state, $w_0\in(-1,-0.9997)\,\,(95% {\rm C.L.})$.

In general relativity gravitational waves propagate at the speed of light, however in alternative theories of gravity that might not be the case. We study the effects of a modified speed of gravity, $c_T^2$, on the B-modes of the Cosmic Microwave Background (CMB) anisotropy in polarisation. We find that a departure from the light speed value would leave a characteristic imprint on the BB spectrum part induced by tensors, manifesting as a shift in the angular scale of its peaks. We derive constraints by using the available {\it Planck} and BICEP2 datasets showing how $c_T^2$ can be measured, albeit obtaining weak constraints due to the overall poor accuracy of the current BB power spectrum measurements. The present constraint corresponds to $c_T^2 = 1.30 \pm 0.79$ and $c_T^2< 2.85$ at $95%$ C.L. by assuming a power law primordial tensor power spectrum and $c_T^2<2.33$ at $95%$ C.L. if the running of the spectral index is allowed. We derive forecasts for the next generation CMB satellites, which we find capable of tightly constraining $c_T^2$ at percent level, comparable with bounds from binary pulsar measurements, largely due to the absence of degeneracy with other cosmological parameters.

We implement the effective field theory (EFT) approach to dark energy and modified gravity in the public Einstein-Boltzmann solver CAMB. The resulting code, which we dub EFTCAMB, is a powerful and versatile tool that can be used for several objectives. It can be employed to evolve the full dynamics of linear scalar perturbations in any given single field dark energy or modified gravity model, once the latter is mapped into the EFT formalism. It offers a numerical implementation of EFT as a model-independent framework to test gravity on cosmological scales. EFTCAMB has a built-in check for the fulfillment of general stability conditions such as the absence of ghost and superluminal propagation of perturbations. It handles phantom-divide crossing models and does not contain any quasi-static approximation, but rather evolves the full dynamics of perturbations on linear scales. As we will show, the latter is an important feature in view of the accuracy and scale range of upcoming surveys. We show the reliability and applicability of our code by evolving the dynamics of linear perturbations and extracting predictions for power spectra in several models. In particular we perform a thorough analysis of f(R) theories, comparing our outputs with those of an existing code for LCDM backgrounds, and finding an agreement that can reach 0.1% for models with a Compton wavelength consistent with current cosmological data. We then showcase the flexibility of our code studying two different scenarios. First we produce new results for designer f(R) models with a time-varying dark energy equation of state. Second, we extract predictions for linear observables in some parametrized EFT models with a phantom-divide crossing equation of state for dark energy.

The effective field theory (EFT) of dark energy relies on three functions of time to describe the background dynamics. The viability of these functions is investigated here by means of a thorough dynamical analysis. While the system is underdetermined, and one can always find a set of functions reproducing any expansion history, we are able to determine general compatibility conditions for these functions by requiring a viable background cosmology. In particular, we identify a set of variables that allows us to transform the non-autonomous system of equations into an infinite-dimensional one characterized by a significant recursive structure. We then analyze several autonomous sub-systems, obtained truncating the original one at increasingly higher dimension, that correspond to increasingly general models of dark energy and modified gravity. Furthermore, we exploit the recursive nature of the system to draw some general conclusions on the different cosmologies that can be recovered within the EFT formalism and the corresponding compatibility requirements for the EFT functions. The machinery that we set up serves different purposes. It offers a general scheme for performing dynamical analysis of dark energy and modified gravity models within the model independent framework of EFT; the general results, obtained with this technique, can be projected into specific models, as we show in one example. It also can be used to determine appropriate ansatze for the three EFT background functions when studying the dynamics of cosmological perturbations in the context of large scale structure tests of gravity.