We study if the discrepancy between the local and cosmological measurements of the Hubble constant $H_0$ can be reformulated as a tension in the cosmic microwave background (CMB) monopole temperature $T_0$. The latter is customarily fixed to the FIRAS best-fit value in CMB data analyses. Although this value was confirmed by several independent experiments, it is interesting to see how much parameter constraints depend on this prior. We first provide a detailed pedagogical description of the $T_0$ effects on cosmological observables. We show that the recombination history and transfer functions do not depend on $T_0$, provided they are parametrized by the energy scale rather than redshift, and at fixed dark matter and baryon densities per CMB photon. Thus, $T_0$ is only a property of the observer, quantifying the amount of expansion between key cosmological events and today. As a consequence, the sole effect of $T_0$ on small-scale primary CMB anisotropies is through the angular diameter distance to the epoch of last scattering, resulting in a near-perfect degeneracy between $T_0$ and $H_0$. This geometric degeneracy is partially lifted by the late-time integrated Sachs-Wolfe effect and CMB lensing. Still, Planck data alone is consistent with a broad region in the $H_0-T_0$ plane, implying that removing the FIRAS prior on $T_0$ can make Planck and SH0ES less discrepant, without introducing new physics beyond $\Lambda$CDM. One may break the degeneracy by combining Planck with SH0ES, yielding an independent measurement of $T_0$, which happens to be in $3\sigma$ tension with FIRAS. Therefore, the Hubble tension indeed can be recast into the $T_0$ tension. The agreement with FIRAS is restored when combining Planck with the baryon acoustic oscillation data instead of SH0ES. Thus, the tension between SH0ES and cosmological measurements of $H_0$ persists even if we discard the FIRAS $T_0$ prior.

Measurements of the cosmic microwave background (CMB) spectral distortions (SDs) will open a new window on the very early universe, providing new information complementary to that gathered from CMB temperature and polarization anisotropies. In this paper, we study their synergy as a function of the characteristics of the considered experiments. In particular, we examine a wide range of sensitivities for possible SD measurements, spanning from FIRAS up to noise levels 1000 times better than PIXIE, and study their constraining power when combined with current or future CMB anisotropy experiments such as Planck or LiteBIRD plus CMB-S4. We consider a number of different cosmological models such as the $\Lambda$CDM, as well as its extensions with the running of the scalar spectral index, the decay or the annihilation of dark matter (DM) particles. While upcoming CMB anisotropy experiments will be able to decrease the uncertainties on inflationary parameters such as $A_s$ and $n_s$ by about a factor 2 in the $\Lambda$CDM case, we find that an SD experiment 10 times more sensitive than PIXIE (comparable to the proposed PRISM satellite) could potentially further contribute to constrain these parameters. This is even more significant in the case of the running of the scalar spectral index. Furthermore, as expected, constraints on DM particles decaying at redshifts probed by SDs will improve by orders of magnitude even with an experiment 10 times worse than PIXIE as compared to CMB anisotropies or Big Bang Nucleosynthesis bounds. On the contrary, DM annihilation constraints will not significantly improve over CMB anisotropy measurements. For some of the cases considered, we study the impact of marginalizing over the contribution from reionization and structure formation. Finally, we forecast the constraints obtainable with sensitivities achievable either from the ground or from a balloon.

Within cosmological perturbation theory, the cosmic microwave background anisotropies are usually computed from a Boltzmann hierarchy coupled to the perturbed Einstein equations. In this setup, one set of multipoles describes the temperature anisotropies, while two other sets, of electric and magnetic types, describe the polarization anisotropies. In order to reduce the number of multipoles types needed for polarization, and thus to speed up the numerical resolution, an optimal hierarchy has been proposed in the literature for Einstein-Boltzmann codes. However, it has been recently shown that the separability between directional and orbital eigenfunctions employed in the optimal hierarchy is not correct in the presence of spatial curvature. We investigate how the assumption of separability affects the optimal hierarchy, and show that it introduces relative errors of order $\Omega_K$ with respect to the full hierarchy. Despite of that, we show that the optimal hierarchy still gives extremely good results for temperature and polarization angular spectra, with relative errors that are much smaller than cosmic variance even for curvatures as large as $|\Omega_K|=0.1$. Still, we find that the polarization angular spectra from tensor perturbations are significantly altered when using the optimal hierarchy, leading to errors that are typically of order $50 |\Omega_K| \%$ on that component.

The Ly-$\alpha$ forest 1D flux power spectrum is a powerful probe of several cosmological parameters. Assuming a $\Lambda$CDM cosmology including massive neutrinos, we find that the latest SDSS DR14 BOSS and eBOSS Ly-$\alpha$ forest data is in very good agreement with current weak lensing constraints on $(\Omega_m, \sigma_8)$ and has the same small level of tension with Planck. We did not identify a systematic effect in the data analysis that could explain this small tension, but we show that it can be reduced in extended cosmological models where the spectral index is not the same on the very different times and scales probed by CMB and Ly-$\alpha$ data. A particular case is that of a $\Lambda$CDM model including a running of the spectral index on top of massive neutrinos. With combined Ly-$\alpha$ and Planck data, we find a slight (3$\sigma$) preference for negative running, $\alpha_s= -0.010 \pm 0.004$ (68% CL). Neutrino mass bounds are found to be robust against different assumptions. In the $\Lambda$CDM model with running, we find $\sum m_\nu <0.11$ eV at the 95% confidence level for combined Ly-$\alpha$ and Planck (temperature and polarisation) data, or $\sum m_\nu < 0.09$ eV when adding CMB lensing and BAO data. We further provide strong and nearly model-independent bounds on the mass of thermal warm dark matter. For a conservative configuration consisting of SDSS data restricted to $z<4.5$ combined with XQ-100 \lya data, we find $m_X > 5.3\;\mathrm{keV}$ (95\%CL).

We study the impact of assumptions made about the neutrino mass ordering on cosmological parameter estimation with the purpose of understanding whether in the future it will be possible to infer the specific neutrino mass distribution from cosmological data. We find that although the commonly used assumption of a degenerate neutrino hierarchy is manifestly wrong and leads to changes in cosmological observables such as the cosmic microwave background and large scale structure compared to the correct (normal or inverted) neutrino hierarchy, the induced changes are so small that even with extremely optimistic assumptions about future data they will remain undetectable. We are thus able to conclude that while cosmology can probe the neutrino contribution to the cosmic energy density extremely precisely (and hence provide a detection of a non-zero total neutrino mass at high significance), it will not be possible to directly measure the individual neutrino masses.

Spectral distortions and anisotropies of the CMB provide independent and complementary probes to study energy injection processes in the early universe. Here we discuss the synergy between these observables, and show the promising future of spectral distortion missions to constrain both exotic and non-exotic energy injections. We show that conventional probes such as Big Bang Nucleosynthesis and CMB anisotropies can benefit from and even be surpassed by future spectral distortion experiments. For this, we have implemented a unified framework within the Boltzmann code CLASS to consistently treat the thermal evolution of photons and baryons. Furthermore, we give an extensive and pedagogical introduction into the topic of spectral distortions and energy injections throughout the thermal history of the universe, highlighting some of their unique features and potential as a novel probe for cosmology and particle physics.

The European Space Agency's Planck satellite, which was dedicated to studying the early Universe and its subsequent evolution, was launched on 14 May 2009. It scanned the microwave and submillimetre sky continuously between 12 August 2009 and 23 October 2013, producing deep, high-resolution, all-sky maps in nine frequency bands from 30 to 857GHz. This paper presents the cosmological legacy of Planck, which currently provides our strongest constraints on the parameters of the standard cosmological model and some of the tightest limits available on deviations from that model. The 6-parameter LCDM model continues to provide an excellent fit to the cosmic microwave background data at high and low redshift, describing the cosmological information in over a billion map pixels with just six parameters. With 18 peaks in the temperature and polarization angular power spectra constrained well, Planck measures five of the six parameters to better than 1% (simultaneously), with the best-determined parameter (theta_*) now known to 0.03%. We describe the multi-component sky as seen by Planck, the success of the LCDM model, and the connection to lower-redshift probes of structure formation. We also give a comprehensive summary of the major changes introduced in this 2018 release. The Planck data, alone and in combination with other probes, provide stringent constraints on our models of the early Universe and the large-scale structure within which all astrophysical objects form and evolve. We discuss some lessons learned from the Planck mission, and highlight areas ripe for further experimental advances.

Reionisation blurring is a non-linear correction to the cosmic microwave background that acts similar to weak gravitational lensing and that can be computed from linear perturbations through a blurring potential. Its impact on the cosmic microwave background is roughly two order of magnitude smaller than that of lensing, in isolation. But the blurring potential is strongly correlated with the lensing potential thereby generating a potentially observable cross-correlation. We compute for the first time the inclusive impact of reionisation blurring on the temperature angular power spectrum and discuss how much it could induce lensing anomalies

A number of anomalous results in short-baseline oscillation may hint at the existence of one or more light sterile neutrino states in the eV mass range and have triggered a wave of new experimental efforts to search for a definite signature of oscillations between active and sterile neutrino states. The present paper aims to provide a comprehensive review on the status of light sterile neutrino searches in mid-2019: we discuss not only the basic experimental approaches and sensitivities of reactor, source, atmospheric, and accelerator neutrino oscillation experiments but also the complementary bounds arising from direct neutrino mass experiments and cosmological observations. Moreover, we review current results from global oscillation analyses that include the constraints set by running reactor and atmospheric neutrino experiments. They permit to set tighter bounds on the active-sterile oscillation parameters but as yet are not able to provide a definite conclusion on the existence of eV-scale sterile neutrinos.

The standard cosmological model predicts the existence of a Cosmic Neutrino Background, which has not yet been observed directly. Some experiments aiming at its detection are currently under development, despite the tiny kinetic energy of the cosmological relic neutrinos, which makes this task incredibly challenging. Since massive neutrinos are attracted by the gravitational potential of our Galaxy, they can cluster locally. Neutrinos should be more abundant at the Earth position than at an average point in the Universe. This fact may enhance the expected event rate in any future experiment. Past calculations of the local neutrino clustering factor only considered a spherical distribution of matter in the Milky Way and neglected the influence of other nearby objects like the Virgo cluster, although recent $N$-body simulations suggest that the latter may actually be important. In this paper, we adopt a back-tracking technique, well established in the calculation of cosmic rays fluxes, to perform the first three-dimensional calculation of the number density of relic neutrinos at the Solar System, taking into account not only the matter composition of the Milky Way, but also the contribution of the Andromeda galaxy and the Virgo cluster. The effect of Virgo is indeed found to be relevant and to depend non-trivially on the value of the neutrino mass. Our results show that the local neutrino density is enhanced by 0.53% for a neutrino mass of 10 meV, 12% for 50 meV, 50% for 100 meV or 500% for 300 meV.

We present constraints on Horndeski gravity from a combined analysis of cosmic shear, galaxy-galaxy lensing and galaxy clustering from $450\,\mathrm{deg}^2$ of the Kilo-Degree Survey (KiDS) and the Galaxy And Mass Assembly (GAMA) survey. The Horndeski class of dark energy/modified gravity models includes the majority of universally coupled extensions to $\Lambda$CDM with one scalar field in addition to the metric. We study the functions of time that fully describe the evolution of linear perturbations in Horndeski gravity. Our results are compatible throughout with a $\Lambda$CDM model. By imposing gravitational wave constraints, we fix the tensor speed excess to zero and consider a subset of models including e.g. quintessence and $f(R)$ theories. Assuming proportionality of the Horndeski functions $\alpha_B$ and $\alpha_M$ (kinetic braiding and the Planck mass run rate, respectively) to the dark energy density fraction $\Omega_{\mathrm{DE}}(a) = 1 - \Omega_{\mathrm{m}}(a)$, we find for the proportionality coefficients $\hat{\alpha}_B = 0.20_{-0.33}^{+0.20} \,$ and $\, \hat{\alpha}_M = 0.25_{-0.29}^{+0.19}$. Our value of $S_8 \equiv \sigma_8 \sqrt{\Omega_{\mathrm{m}}/0.3}$ is in better agreement with the $Planck$ estimate when measured in the enlarged Horndeski parameter space than in a pure $\Lambda$CDM scenario. In our joint three-probe analysis we report a downward shift of the $S_8$ best fit value from the $Planck$ measurement of $\Delta S_8 = 0.016_{-0.046}^{+0.048}$ in Horndeski gravity, compared to $\Delta S_8 = 0.059_{-0.039}^{+0.040}$ in $\Lambda$CDM. Our constraints are robust to the modelling uncertainty of the non-linear matter power spectrum in Horndeski gravity. Our likelihood code for multi-probe analysis in both $\Lambda$CDM and Horndeski gravity is publicly available at http://github.com/alessiospuriomancini/KiDSHorndeski .

We present cosmological parameter results from the final full-mission Planck measurements of the CMB anisotropies. We find good consistency with the standard spatially-flat 6-parameter $\Lambda$CDM cosmology having a power-law spectrum of adiabatic scalar perturbations (denoted "base $\Lambda$CDM" in this paper), from polarization, temperature, and lensing, separately and in combination. A combined analysis gives dark matter density $\Omega_c h^2 = 0.120\pm 0.001$, baryon density $\Omega_b h^2 = 0.0224\pm 0.0001$, scalar spectral index $n_s = 0.965\pm 0.004$, and optical depth $\tau = 0.054\pm 0.007$ (in this abstract we quote $68\,\%$ confidence regions on measured parameters and $95\,\%$ on upper limits). The angular acoustic scale is measured to $0.03\,\%$ precision, with $100\theta_*=1.0411\pm 0.0003$. These results are only weakly dependent on the cosmological model and remain stable, with somewhat increased errors, in many commonly considered extensions. Assuming the base-$\Lambda$CDM cosmology, the inferred late-Universe parameters are: Hubble constant $H_0 = (67.4\pm 0.5)$km/s/Mpc; matter density parameter $\Omega_m = 0.315\pm 0.007$; and matter fluctuation amplitude $\sigma_8 = 0.811\pm 0.006$. We find no compelling evidence for extensions to the base-$\Lambda$CDM model. Combining with BAO we constrain the effective extra relativistic degrees of freedom to be $N_{\rm eff} = 2.99\pm 0.17$, and the neutrino mass is tightly constrained to $\sum m_\nu< 0.12$eV. The CMB spectra continue to prefer higher lensing amplitudes than predicted in base -$\Lambda$CDM at over $2\,\sigma$, which pulls some parameters that affect the lensing amplitude away from the base-$\Lambda$CDM model; however, this is not supported by the lensing reconstruction or (in models that also change the background geometry) BAO data. (Abridged)

We extend the formalism of time-sliced perturbation theory (TSPT) for cosmological large-scale structure to include non-Gaussian initial conditions. We show that in such a case the TSPT interaction vertices acquire new contributions whose time-dependence factorizes for the Einstein-de Sitter cosmology. The new formulation is free from spurious infrared (IR) enhancements and reveals a clear IR structure of non-Gaussian vertices. We use the new technique to study the evolution of oscillating features in primordial statistics and show that they are damped due to non-linear effects of large bulk flows. We derive the damping factors for the oscillating primordial power spectrum and bispectrum by means of a systematic IR resummation of relevant Feynman diagrams.

Voyage 2050 White Paper highlighting the unique science opportunities using spectral distortions of the cosmic microwave background (CMB). CMB spectral distortions probe many processes throughout the history of the Universe. Precision spectroscopy, possible with existing technology, would provide key tests for processes expected within the cosmological standard model and open an enormous discovery space to new physics. This offers unique scientific opportunities for furthering our understanding of inflation, recombination, reionization and structure formation as well as dark matter and particle physics. A dedicated experimental approach could open this new window to the early Universe in the decades to come, allowing us to turn the long-standing upper distortion limits obtained with COBE/FIRAS some 25 years ago into clear detections of the expected standard distortion signals.

This paper discusses the science case for a sensitive spectro-polarimetric survey of the microwave sky. Such a survey would provide a tomographic and dynamic census of the three-dimensional distribution of hot gas, velocity flows, early metals, dust, and mass distribution in the entire Hubble volume, exploit CMB temperature and polarisation anisotropies down to fundamental limits, and track energy injection and absorption into the radiation background across cosmic times by measuring spectral distortions of the CMB blackbody emission. In addition to its exceptional capability for cosmology and fundamental physics, such a survey would provide an unprecedented view of microwave emissions at sub-arcminute to few-arcminute angular resolution in hundreds of frequency channels, a data set that would be of immense legacy value for many branches of astrophysics. We propose that this survey be carried-out with a large space mission featuring a broad-band polarised imager and a moderate resolution spectro-imager at the focus of a 3.5m aperture telescope actively cooled to about 8K, complemented with absolutely-calibrated Fourier Transform Spectrometer modules observing at degree-scale angular resolution in the 10-2000 GHz frequency range. We propose two observing modes: a survey mode to map the entire sky as well as a few selected wide fields, and an observatory mode for deeper observations of regions of specific interest.

We present constraints on neutrino masses, the primordial fluctuation spectrum from inflation, and other parameters of the $\Lambda$CDM model, using the one-dimensional Ly$\alpha$-forest power spectrum measured by Palanque-Delabrouille et al. (2013) from SDSS-III/BOSS, complemented by Planck 2015 cosmic microwave background (CMB) data and other cosmological probes. This paper improves on the previous analysis by Palanque-Delabrouille et al. (2015) by using a more powerful set of calibrating hydrodynamical simulations that reduces uncertainties associated with resolution and box size, by adopting a more flexible set of nuisance parameters for describing the evolution of the intergalactic medium, by including additional freedom to account for systematic uncertainties, and by using Planck 2015 constraints in place of Planck 2013. Fitting Ly$\alpha$ data alone leads to cosmological parameters in excellent agreement with the values derived independently from CMB data, except for a weak tension on the scalar index $n_s$. Combining BOSS Ly$\alpha$ with Planck CMB constrains the sum of neutrino masses to $\sum m_\nu < 0.12$ eV (95\% C.L.) including all identified systematic uncertainties, tighter than our previous limit (0.15 eV) and more robust. Adding Ly$\alpha$ data to CMB data reduces the uncertainties on the optical depth to reionization $\tau$, through the correlation of $\tau$ with $\sigma_8$. Similarly, correlations between cosmological parameters help in constraining the tensor-to-scalar ratio of primordial fluctuations $r$. The tension on $n_s$ can be accommodated by allowing for a running ${\mathrm d}n_s/{\mathrm d}\ln k$. Allowing running as a free parameter in the fits does not change the limit on $\sum m_\nu$. We discuss possible interpretations of these results in the context of slow-roll inflation.

Many attempts to solve the Hubble tension with extended cosmological models combine an enhanced relic radiation density, acting at the level of background cosmology, with new physical ingredients affecting the evolution of cosmological perturbations. Several authors have pointed out the ability of combined Baryon Acoustic Oscillation (BAO) and Big Bang Nucleosynthesis (BBN) data to probe the background cosmological history independently of both CMB maps and supernovae data. Using state-of-the-art assumptions on BBN, we confirm that combined BAO, deuterium, and helium data are in tension with the SH0ES measurements under the $\Lambda$CDM assumption at the 3.2$\sigma$ level, while being in close agreement with the CMB value. We subsequently show that floating the radiation density parameter $N_\mathrm{eff}$ only reduces the tension down to the 2.6$\sigma$ level. This conclusion, totally independent of any CMB data, shows that a high $N_\mathrm{eff}$ accounting for extra relics (either free-streaming or self-interacting) does not provide an obvious solution to the crisis, not even at the level of background cosmology. To circumvent this strong bound, (i) the extra radiation has to be generated after BBN to avoid helium bounds, and (ii) additional ingredients have to be invoked at the level of perturbations to reconcile this extra radiation with CMB and LSS data.

We report on the implications for cosmic inflation of the 2018 Release of the Planck CMB anisotropy measurements. The results are fully consistent with the two previous Planck cosmological releases, but have smaller uncertainties thanks to improvements in the characterization of polarization at low and high multipoles. Planck temperature, polarization, and lensing data determine the spectral index of scalar perturbations to be $n_\mathrm{s}=0.9649\pm 0.0042$ at 68% CL and show no evidence for a scale dependence of $n_\mathrm{s}.$ Spatial flatness is confirmed at a precision of 0.4% at 95% CL with the combination with BAO data. The Planck 95% CL upper limit on the tensor-to-scalar ratio, $r_{0.002}<0.10$, is further tightened by combining with the BICEP2/Keck Array BK15 data to obtain $r_{0.002}<0.056$. In the framework of single-field inflationary models with Einstein gravity, these results imply that: (a) slow-roll models with a concave potential, $V" (\phi) < 0,$ are increasingly favoured by the data; and (b) two different methods for reconstructing the inflaton potential find no evidence for dynamics beyond slow roll. Non-parametric reconstructions of the primordial power spectrum consistently confirm a pure power law. A complementary analysis also finds no evidence for theoretically motivated parameterized features in the Planck power spectrum, a result further strengthened for certain oscillatory models by a new combined analysis that includes Planck bispectrum data. The new Planck polarization data provide a stringent test of the adiabaticity of the initial conditions. The polarization data also provide improved constraints on inflationary models that predict a small statistically anisotropic quadrupolar modulation of the primordial fluctuations. However, the polarization data do not confirm physical models for a scale-dependent dipolar modulation.

Einstein-Boltzmann Solvers (EBSs) are run on a massive scale by the cosmology community when fitting cosmological models to data. We present a new concept for speeding up such codes with neural networks. The originality of our approach stems from not substituting the whole EBS by a machine learning algorithm, but only its most problematic and least parallelizable step: the integration of perturbation equations over time. This approach offers two significant advantages: the task depends only on a subset of cosmological parameters, and it is blind to the characteristics of the experiment for which the output must be computed (for instance, redshift bins). These allow us to construct a fast and highly re-usable network. In this proof-of-concept paper, we focus on the prediction of CMB source functions, and design our networks according to physical considerations and analytical approximations. This allows us to reduce the depth and training time of the networks compared to a brute-force approach. Indeed, the calculation of the source functions using the networks is fast enough so that it is not a bottleneck in the EBS anymore. Finally, we show that their accuracy is more than sufficient for accurate MCMC parameter inference from Planck data. This paves the way for a new project, CosmicNet, aimed at gradually extending the use and the range of validity of neural networks within EBSs, and saving massive computational time in the context of cosmological parameter extraction.

Several interesting Dark Matter (DM) models invoke a dark sector leading to two types of relic particles, possibly interacting with each other: non-relativistic DM, and relativistic Dark Radiation (DR). These models have interesting consequences for cosmological observables, and could in principle solve problems like the small-scale cold DM crisis, Hubble tension, and/or low $\sigma_8$ value. Their cosmological behaviour is captured by the ETHOS parametrisation, which includes a DR-DM scattering rate scaling like a power-law of the temperature, $T^n$. Scenarios with $n=0$, $2$, or $4$ can easily be realised in concrete dark sector set-ups. Here we update constraints on these three scenarios using recent CMB, BAO, and high-resolution Lyman-$\alpha$ data. We introduce a new Lyman-$\alpha$ likelihood that is applicable to a wide range of cosmological models with a suppression of the matter power spectrum on small scales. For $n=2$ and $4$, we find that Lyman-$\alpha$ data strengthen the CMB+BAO bounds on the DM-DR interaction rate by many orders of magnitude. However, models offering a possible solution to the missing satellite problem are still compatible with our new bounds. For $n=0$, high-resolution Lyman-$\alpha$ data bring no stronger constraints on the interaction rate than CMB+BAO data, except for extremely small values of the DR density. Using CMB+BAO data and a theory-motivated prior on the minimal density of DR, we find that the $n=0$ model can reduce the Hubble tension from $4.1\sigma$ to $2.7\sigma$, while simultaneously accommodating smaller values of the $\sigma_8$ and $S_8$ parameters hinted by cosmic shear data.

We consistently include the effect of massive neutrinos in the thermal Sunyaev Zeldovich (SZ) power spectrum and cluster counts analyses, highlighting subtle dependencies on the parameterisation and data combination. In $\Lambda$CDM, with an X-ray mass bias corresponding the expected hydrostatic mass bias, i.e., $b\simeq0.2$, our constraints from Planck SZ data are consistent with the latest results from SPT, DES-Y1 and KiDS+VIKING-450. In $\nu\Lambda$CDM, without prior information on $b$, our joint analyses of Planck SZ with Planck 2015 primary CMB yield a small improvement on the total neutrino mass bound compared to the Planck 2015 primary CMB constraint, as well as $(1-b)=0.64\pm0.04$~(68\%~CL). For forecasts, we find that a combination of a mock cosmic variance limited SZ power spectrum with primary CMB and BAO can improve the uncertainty on the total neutrino mass by 14\% with respect to CMB combined with BAO. This requires masking the heaviest clusters and probing the small-scale SZ power spectrum up to multipoles of $\ell_\mathrm{max}=10^4$. For a future competitive measurement of the total neutrino mass using CMB and the SZ power spectrum, but excluding BAO and lensing power spectrum, we find that a 1\% precision on the mass calibration is needed. Although SZ power spectrum-based measurements of the neutrino masses are challenging, we find that the SZ power spectrum can be used to tightly constrain intra-cluster medium properties.

We forecast the sensitivity of thirty-five different combinations of future Cosmic Microwave Background and Large Scale Structure data sets to cosmological parameters and to the total neutrino mass. We work under conservative assumptions accounting for uncertainties in the modelling of systematics. In particular, for galaxy redshift surveys, we remove the information coming from non-linear scales. We use Bayesian parameter extraction from mock likelihoods to avoid Fisher matrix uncertainties. Our grid of results allows for a direct comparison between the sensitivity of different data sets. We find that future surveys will measure the neutrino mass with high significance and will not be substantially affected by potential parameter degeneracies between neutrino masses, the density of relativistic relics, and a possible time-varying equation of state of Dark Energy.

Cosmological observations offer unique and robust avenues for probing the fundamental nature of dark matter particles-they broadly test a range of compelling theoretical scenarios, often surpassing or complementing the reach of terrestrial and other experiments. We discuss observational and theoretical advancements that will play a pivotal role in realizing a strong program of cosmological searches for the identity of dark matter in the coming decade. Specifically, we focus on measurements of the cosmic-microwave-background anisotropy and spectral distortions, and tracers of structure (such as the Lyman-$\alpha$ forest, galaxies, and the cosmological 21-cm signal).

Theoretical uncertainties on non-linear scales are among the main obstacles to exploit the sensitivity of forthcoming galaxy and hydrogen surveys like Euclid or the Square Kilometre Array (SKA). Here, we devise a new method to model the theoretical error that goes beyond the usual cut-off on small scales. The advantage of this more efficient implementation of the non-linear uncertainties is tested through a Markov-Chain-Monte-Carlo (MCMC) forecast of the sensitivity of Euclid and SKA to the parameters of the standard $\Lambda$CDM model, including massive neutrinos with total mass $M_\nu$, and to 3 extended scenarios, including 1) additional relativistic degrees of freedom ($\Lambda$CDM + $M_\nu$ + $N_\mathrm{eff}$), 2) a deviation from the cosmological constant ($\Lambda$CDM + $M_\nu$ + $w_0$), and 3) a time-varying dark energy equation of state parameter ($\Lambda$CDM + $M_\nu$ + $\left(w_0,w_a \right)$). We compare the sensitivity of 14 different combinations of cosmological probes and experimental configurations. For Euclid combined with Planck, assuming a plain cosmological constant, our method gives robust predictions for a high sensitivity to the primordial spectral index $n_{\rm s}$ ($\sigma(n_s)=0.00085$), the Hubble constant $H_0$ ($\sigma(H_0)=0.141 \, {\rm km/s/Mpc}$), the total neutrino mass $M_\nu$ ($\sigma(M_\nu)=0.020 \, {\rm eV}$). Assuming dynamical dark energy we get $\sigma(M_\nu)=0.030 \, {\rm eV}$ for the mass and $(\sigma(w_0), \sigma(w_a)) = (0.0214, 0.071)$ for the equation of state parameters. The predicted sensitivity to $M_\nu$ is mostly stable against the extensions of the cosmological model considered here. Interestingly, a significant improvement of the constraints on the extended model parameters is also obtained when combining Euclid with a low redshift HI intensity mapping survey by SKA1, demonstrating the importance of the synergy of Euclid and SKA.

In cosmologies with massive neutrinos, the galaxy bias defined with respect to the total matter field (cold dark matter, baryons, and non-relativistic neutrinos) depends on the sum of the neutrino masses $M_{\nu}$, and becomes scale-dependent even on large scales. This effect has been usually neglected given the sensitivity of current surveys, but becomes a severe systematic for future surveys aiming to provide the first detection of non-zero $M_{\nu}$. The effect can be corrected for by defining the bias with respect to the density field of cold dark matter and baryons instead of the total matter field. In this work, we provide a simple prescription for correctly mitigating the neutrino-induced scale-dependent bias effect in a practical way. We clarify a number of subtleties regarding how to properly implement this correction in the presence of redshift-space distortions and non-linear evolution of perturbations. We perform a MCMC analysis on simulated galaxy clustering data that match the expected sensitivity of the \textit{Euclid} survey. We find that the neutrino-induced scale-dependent bias can lead to important shifts in both the inferred mean value of $M_{\nu}$, as well as its uncertainty. We show how these shifts propagate to other cosmological parameters correlated with $M_{\nu}$, such as the cold dark matter physical density $\Omega_{cdm} h^2$ and the scalar spectral index $n_s$. In conclusion, we find that correctly accounting for the neutrino-induced scale-dependent bias will be of crucial importance for future galaxy clustering analyses. We encourage the cosmology community to correctly account for this effect using the simple prescription we present in our work. The tools necessary to easily correct for the neutrino-induced scale-dependent bias will be made publicly available in an upcoming release of the Boltzmann solver \texttt{CLASS}.

We present a new efficient method to compute the angular power spectra of large-scale structure observables that circumvents the numerical integration over Bessel functions, expanding on a recently proposed algorithm based on FFTlog. This new approach has better convergence properties. The method is explicitly implemented in the CLASS code for the case of number count $C_\ell$'s (including redshift-space distortions, weak lensing, and all other relativistic corrections) and cosmic shear $C_\ell$'s. In both cases our approach speeds up the calculation of the exact $C_\ell$'s (without the Limber approximation) by a factor of order 400 at a fixed precision target of 0.1%.

MontePython is a parameter inference package for cosmology. We present the latest development of the code over the past couple of years. We explain, in particular, two new ingredients both contributing to improve the performance of Metropolis-Hastings sampling: an adaptation algorithm for the jumping factor, and a calculation of the inverse Fisher matrix, which can be used as a proposal density. We present several examples to show that these features speed up convergence and can save many hundreds of CPU-hours in the case of difficult runs, with a poor prior knowledge of the covariance matrix. We also summarise all the functionalities of MontePython in the current release, including new likelihoods and plotting options.

We derive updated constraints on the Integrated Sachs-Wolfe (ISW) effect through cross-correlation of the cosmic microwave background with galaxy surveys. We improve with respect to similar previous analyses in several ways. First, we use the most recent versions of extragalactic object catalogs: SDSS DR12 photometric redshift (photo-$z$) and 2MASS Photo-$z$ datasets, as well as employed earlier for ISW, SDSS QSO photo-$z$ and NVSS samples. Second, we use for the first time the WISE~$\times$~SuperCOSMOS catalog, which allows us to perform an all-sky analysis of the ISW up to $z\sim0.4$. Third, thanks to the use of photo-$z$s, we separate each dataset into different redshift bins, deriving the cross-correlation in each bin. This last step leads to a significant improvement in sensitivity. We remove cross-correlation between catalogs using masks which mutually exclude common regions of the sky. We use two methods to quantify the significance of the ISW effect. In the first one, we fix the cosmological model, derive linear galaxy biases of the catalogs, and then evaluate the significance of the ISW using a single parameter. In the second approach we perform a global fit of the ISW and of the galaxy biases varying the cosmological model. We find significances of the ISW in the range 4.7-5.0 $\sigma$ thus reaching, for the first time in such an analysis, the threshold of 5 $\sigma$. Without the redshift tomography we find a significance of $\sim$ 4.0 $\sigma$, which shows the importance of the binning method. Finally we use the ISW data to infer constraints on the Dark Energy redshift evolution and equation of state. We find that the redshift range covered by the catalogs is still not optimal to derive strong constraints, although this goal will be likely reached using future datasets such as from Euclid, LSST, and SKA.

We present cluster counts and corresponding cosmological constraints from the Planck full mission data set. Our catalogue consists of 439 clusters detected via their Sunyaev-Zeldovich (SZ) signal down to a signal-to-noise ratio of 6, and is more than a factor of 2 larger than the 2013 Planck cluster cosmology sample. The counts are consistent with those from 2013 and yield compatible constraints under the same modelling assumptions. Taking advantage of the larger catalogue, we extend our analysis to the two-dimensional distribution in redshift and signal-to-noise. We use mass estimates from two recent studies of gravitational lensing of background galaxies by Planck clusters to provide priors on the hydrostatic bias parameter, $(1-b)$. In addition, we use lensing of cosmic microwave background (CMB) temperature fluctuations by Planck clusters as an independent constraint on this parameter. These various calibrations imply constraints on the present-day amplitude of matter fluctuations in varying degrees of tension with those from the Planck analysis of primary fluctuations in the CMB; for the lowest estimated values of $(1-b)$ the tension is mild, only a little over one standard deviation, while it remains substantial ($3.7\,\sigma$) for the largest estimated value. We also examine constraints on extensions to the base flat $\Lambda$CDM model by combining the cluster and CMB constraints. The combination appears to favour non-minimal neutrino masses, but this possibility does little to relieve the overall tension because it simultaneously lowers the implied value of the Hubble parameter, thereby exacerbating the discrepancy with most current astrophysical estimates. Improving the precision of cluster mass calibrations from the current 10%-level to 1% would significantly strengthen these combined analyses and provide a stringent test of the base $\Lambda$CDM model.

We consider a recently proposed model in which dark matter interacts with a thermal background of dark radiation. Dark radiation consists of relativistic degrees of freedom which allow larger values of the expansion rate of the universe today to be consistent with CMB data ($H_0$-problem). Scattering between dark matter and radiation suppresses the matter power spectrum at small scales and can explain the apparent discrepancies between $\Lambda$CDM predictions of the matter power spectrum and direct measurements of Large Scale Structure LSS ($\sigma_8$-problem). We go beyond previous work in two ways: 1. we enlarge the parameter space of our previous model and allow for an arbitrary fraction of the dark matter to be interacting and 2. we update the data sets used in our fits, most importantly we include LSS data with full $k$-dependence to explore the sensitivity of current data to the shape of the matter power spectrum. We find that LSS data prefer models with overall suppressed matter clustering due to dark matter - dark radiation interactions over $\Lambda$CDM at 3-4 $\sigma$. However recent weak lensing measurements of the power spectrum are not yet precise enough to clearly distinguish two limits of the model with different predicted shapes for the linear matter power spectrum. In two Appendices we give a derivation of the coupled dark matter and dark radiation perturbation equations from the Boltzmann equation in order to clarify a confusion in the recent literature, and we derive analytic approximations to the solutions of the perturbation equations in the two physically interesting limits of all dark matter weakly interacting or a small fraction of dark matter strongly interacting.

We devise a new user-friendly tool interfaced with the Boltzmann code CLASS to deal with any kind of exotic electromagnetic energy injection in the universe and its impact on anisotropies of the Cosmic Microwave Background. It makes use of the results from standard electromagnetic cascade calculations develop in the context of WIMP annihilation, generalized to incorporate any injection history. We first validate it on a specific WIMP scenario, the Higgs Portal model, confirming that the standard effective on-the-spot treatment is accurate enough. We then analyze the more involved example of evaporating Primordial Black Holes (PBHs) with masses in the range $[3\times10^{13},5\times10^{16}]$g, for which the standard approximations break down. We derive robust CMB bounds on the relic density of evaporating PBHs, ruling out the possibility for PBHs with a monochromatic distribution of masses in the range $[3\times10^{13},2.5\times10^{16}]$g to represent all of the Dark Matter in our Universe. Remarkably, we confirm with an accurate study that the CMB bounds are several orders of magnitude stronger than those from the galactic gamma-ray background in the range $[3\times10^{13},3\times10^{14}]$g. A future CMB experiment like CORE+, or an experiment attempting at measuring the 21 cm signal from the Dark Ages could greatly improve the sensitivity to these models.

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 use BOSS DR9 quasars to constrain 2 cases of dark matter models: cold-plus-warm (C+WDM) where the warm component is a thermal relic, and sterile neutrinos resonantly produced in the presence of a lepton asymmetry (RPSN). We establish constraints on the relic mass m_x and its relative abundance $F=\Omega_{wdm}/\Omega_{dm}$ using a suite of hydrodynamical simulations in 28 C+WDM configurations. We find that the 3 sigma bounds approximately follow F ~ $0.35 (keV/m_x)^{-1.37}$ from BOSS alone. We also establish constraints on sterile neutrino mass and mixing angle by producing the non-linear flux power spectrum of 8 RPSN models, where the input linear power spectrum is computed directly from the particles distribution functions. We find values of lepton asymmetries for which sterile neutrinos as light as 6.5 keV (resp. 3.5 keV) are consistent with BOSS at the 2 sigma (resp. 3sigma) level. These limits tighten by close to a factor of 2 for lepton asymmetries departing from those yielding the coolest distribution functions. Our Ly-a forest bounds can be strengthened if we include higher-resolution data from XQ-100, HIRES and MIKE. At these smaller scales, the flux power spectrum exhibits a suppression that can be due to Doppler broadening, IGM pressure smoothing or free-streaming of WDM particles. In the current work, we show that if one extrapolates temperatures from lower redshifts via broken power laws in T_0 and gamma, then our 3 sigma C+WDM bounds strengthen to F~ $0.20 (keV/m_x)^{-1.37}$, and the lightest RPSN consistent with our extended data set have masses of 7.0 keV at the 3 sigma level. Using dedicated hydrodynamical simulations, we show that a 7 keV sterile neutrino produced in a lepton asymmetry $L = 8 \times 10^{-6}$ is consistent at 1.9 sigma (resp. 3.1 sigma) with BOSS (resp. BOSS + higher-resolution), for the thermal history models tested in this work.

Nonzero neutrino masses are required by the existence of flavour oscillations, with values of the order of at least 50 meV. We consider the gravitational clustering of relic neutrinos within the Milky Way, and used the $N$-one-body simulation technique to compute their density enhancement factor in the neighbourhood of the Earth with respect to the average cosmic density. Compared to previous similar studies, we pushed the simulation down to smaller neutrino masses, and included an improved treatment of the baryonic and dark matter distributions in the Milky Way. Our results are important for future experiments aiming at detecting the cosmic neutrino background, such as the Princeton Tritium Observatory for Light, Early-universe, Massive-neutrino Yield (PTOLEMY) proposal. We calculate the impact of neutrino clustering in the Milky Way on the expected event rate for a PTOLEMY-like experiment. We find that the effect of clustering remains negligible for the minimal normal hierarchy scenario, while it enhances the event rate by 10 to 20% (resp. a factor 1.7 to 2.5) for the minimal inverted hierarchy scenario (resp. a degenerate scenario with 150 meV masses). Finally we compute the impact on the event rate of a possible fourth sterile neutrino with a mass of 1.3 eV.

We present the implications for cosmic inflation of the Planck measurements of the cosmic microwave background (CMB) anisotropies in both temperature and polarization based on the full Planck survey. The Planck full mission temperature data and a first release of polarization data on large angular scales measure the spectral index of curvature perturbations to be $n_\mathrm{s} = 0.968 \pm 0.006$ and tightly constrain its scale dependence to $d n_s/d \ln k =-0.003 \pm 0.007$ when combined with the Planck lensing likelihood. When the high-$\ell$ polarization data is included, the results are consistent and uncertainties are reduced. The upper bound on the tensor-to-scalar ratio is $r_{0.002} < 0.11$ (95% CL), consistent with the B-mode polarization constraint $r< 0.12$ (95% CL) obtained from a joint BICEP2/Keck Array and Planck analysis. These results imply that $V(\phi) \propto \phi^2$ and natural inflation are now disfavoured compared to models predicting a smaller tensor-to-scalar ratio, such as $R^2$ inflation. Three independent methods reconstructing the primordial power spectrum are investigated. The Planck data are consistent with adiabatic primordial perturbations. We investigate inflationary models producing an anisotropic modulation of the primordial curvature power spectrum as well as generalized models of inflation not governed by a scalar field with a canonical kinetic term. The 2015 results are consistent with the 2013 analysis based on the nominal mission data.

We discuss the effects on the CMB, CIB, and thermal SZ effect due to the peculiar motion of an observer with respect to the CMB rest frame, which induces boosting effects. We investigate the scientific perspectives opened by future CMB space missions, focussing on the CORE proposal. The improvements in sensitivity offered by a mission like CORE, together with its high resolution over a wide frequency range, will provide a more accurate estimate of the CMB dipole. The extension of boosting effects to polarization and cross-correlations will enable a more robust determination of purely velocity-driven effects that are not degenerate with the intrinsic CMB dipole, allowing us to achieve a S/N ratio of 13; this improves on the Planck detection and essentially equals that of an ideal cosmic-variance-limited experiment up to a multipole l of 2000. Precise inter-frequency calibration will offer the opportunity to constrain or even detect CMB spectral distortions, particularly from the cosmological reionization, because of the frequency dependence of the dipole spectrum, without resorting to precise absolute calibration. The expected improvement with respect to COBE-FIRAS in the recovery of distortion parameters (in principle, a factor of several hundred for an ideal experiment with the CORE configuration) ranges from a factor of several up to about 50, depending on the quality of foreground removal and relative calibration. Even for 1% accuracy in both foreground removal and relative calibration at an angular scale of 1 deg, we find that dipole analyses for a mission like CORE will be able to improve the recovery of the CIB spectrum amplitude by a factor of 17 in comparison with current results based on FIRAS. In addition to the scientific potential of a mission like CORE for these analyses, synergies with other planned and ongoing projects are also discussed.

We present an analysis of the main systematic effects that could impact the measurement of CMB polarization with the proposed CORE space mission. We employ timeline-to-map simulations to verify that the CORE instrumental set-up and scanning strategy allow us to measure sky polarization to a level of accuracy adequate to the mission science goals. We also show how the CORE observations can be processed to mitigate the level of contamination by potentially worrying systematics, including intensity-to-polarization leakage due to bandpass mismatch, asymmetric main beams, pointing errors and correlated noise. We use analysis techniques that are well validated on data from current missions such as Planck to demonstrate how the residual contamination of the measurements by these effects can be brought to a level low enough not to hamper the scientific capability of the mission, nor significantly increase the overall error budget. We also present a prototype of the CORE photometric calibration pipeline, based on that used for Planck, and discuss its robustness to systematics, showing how CORE can achieve its calibration requirements. While a fine-grained assessment of the impact of systematics requires a level of knowledge of the system that can only be achieved in a future study phase, the analysis presented here strongly suggests that the main areas of concern for the CORE mission can be addressed using existing knowledge, techniques and algorithms.

Lensing of the CMB is now a well-developed probe of large-scale clustering over a broad range of redshifts. By exploiting the non-Gaussian imprints of lensing in the polarization of the CMB, the CORE mission can produce a clean map of the lensing deflections over nearly the full-sky. The number of high-S/N modes in this map will exceed current CMB lensing maps by a factor of 40, and the measurement will be sample-variance limited on all scales where linear theory is valid. Here, we summarise this mission product and discuss the science that it will enable. For example, the summed mass of neutrinos will be determined to an accuracy of 17 meV combining CORE lensing and CMB two-point information with contemporaneous BAO measurements, three times smaller than the minimum total mass allowed by neutrino oscillations. In the search for B-mode polarization from primordial gravitational waves with CORE, lens-induced B-modes will dominate over instrument noise, limiting constraints on the gravitational wave power spectrum amplitude. With lensing reconstructed by CORE, one can "delens" the observed polarization internally, reducing the lensing B-mode power by 60%. This improves to 70% by combining lensing and CIB measurements from CORE, reducing the error on the gravitational wave amplitude by 2.5 compared to no delensing (in the null hypothesis). Lensing measurements from CORE will allow calibration of the halo masses of the 40000 galaxy clusters that it will find, with constraints dominated by the clean polarization-based estimators. CORE can accurately remove Galactic emission from CMB maps with its 19 frequency channels. We present initial findings that show that residual Galactic foreground contamination will not be a significant source of bias for lensing power spectrum measurements with CORE. [abridged]

One of the open questions in modern cosmology is the small scale crisis of the cold dark matter paradigm. Increasing attention has recently been devoted to self-interacting dark matter models as a possible answer. However, solving the so-called "missing satellites" problem requires in addition the presence of an extra relativistic particle (dubbed dark radiation) scattering with dark matter in the early universe. Here we investigate the impact of different theoretical models devising dark matter dark radiation interactions on large scale cosmological observables. We use cosmic microwave background data to put constraints on the dark radiation component and its coupling to dark matter. We find that the values of the coupling allowed by the data imply a cut-off scale of the halo mass function consistent with the one required to match the observations of satellites in the Milky Way.

We demonstrate that, for the baseline design of the CORE satellite mission, the polarized foregrounds can be controlled at the level required to allow the detection of the primordial cosmic microwave background (CMB) $B$-mode polarization with the desired accuracy at both reionization and recombination scales, for tensor-to-scalar ratio values of ${r\gtrsim 5\times 10^{-3}}$. We consider detailed sky simulations based on state-of-the-art CMB observations that consist of CMB polarization with $\tau=0.055$ and tensor-to-scalar values ranging from $r=10^{-2}$ to $10^{-3}$, Galactic synchrotron, and thermal dust polarization with variable spectral indices over the sky, polarized anomalous microwave emission, polarized infrared and radio sources, and gravitational lensing effects. Using both parametric and blind approaches, we perform full component separation and likelihood analysis of the simulations, allowing us to quantify both uncertainties and biases on the reconstructed primordial $B$-modes. Under the assumption of perfect control of lensing effects, CORE would measure an unbiased estimate of $r=\left(5 \pm 0.4\right)\times 10^{-3}$ after foreground cleaning. In the presence of both gravitational lensing effects and astrophysical foregrounds, the significance of the detection is lowered, with CORE achieving a $4\sigma$-measurement of $r=5\times 10^{-3}$ after foreground cleaning and $60$% delensing. For lower tensor-to-scalar ratios ($r=10^{-3}$) the overall uncertainty on $r$ is dominated by foreground residuals, not by the 40% residual of lensing cosmic variance. Moreover, the residual contribution of unprocessed polarized point-sources can be the dominant foreground contamination to primordial B-modes at this $r$ level, even on relatively large angular scales, $\ell \sim 50$. Finally, we report two sources of potential bias for the detection of the primordial $B$-modes.[abridged]

Future observations of cosmic microwave background (CMB) polarisation have the potential to answer some of the most fundamental questions of modern physics and cosmology. In this paper, we list the requirements for a future CMB polarisation survey addressing these scientific objectives, and discuss the design drivers of the CORE space mission proposed to ESA in answer to the "M5" call for a medium-sized mission. The rationale and options, and the methodologies used to assess the mission's performance, are of interest to other future CMB mission design studies. CORE is designed as a near-ultimate CMB polarisation mission which, for optimal complementarity with ground-based observations, will perform the observations that are known to be essential to CMB polarisation scienceand cannot be obtained by any other means than a dedicated space mission.

We describe a space-borne, multi-band, multi-beam polarimeter aiming at a precise and accurate measurement of the polarization of the Cosmic Microwave Background. The instrument is optimized to be compatible with the strict budget requirements of a medium-size space mission within the Cosmic Vision Programme of the European Space Agency. The instrument has no moving parts, and uses arrays of diffraction-limited Kinetic Inductance Detectors to cover the frequency range from 60 GHz to 600 GHz in 19 wide bands, in the focal plane of a 1.2 m aperture telescope cooled at 40 K, allowing for an accurate extraction of the CMB signal from polarized foreground emission. The projected CMB polarization survey sensitivity of this instrument, after foregrounds removal, is 1.7 {\mu}K$\cdot$arcmin. The design is robust enough to allow, if needed, a downscoped version of the instrument covering the 100 GHz to 600 GHz range with a 0.8 m aperture telescope cooled at 85 K, with a projected CMB polarization survey sensitivity of 3.2 {\mu}K$\cdot$arcmin.

We discuss the potential of a next generation space-borne Cosmic Microwave Background (CMB) experiment for studies of extragalactic sources. Our analysis has particular bearing on the definition of the future space project, CORE, that has been submitted in response to ESA's call for a Medium-size mission opportunity as the successor of the Planck satellite. Even though the effective telescope size will be somewhat smaller than that of Planck, CORE will have a considerably better angular resolution at its highest frequencies, since, in contrast with Planck, it will be diffraction limited at all frequencies. The improved resolution implies a considerable decrease of the source confusion, i.e. substantially fainter detection limits. In particular, CORE will detect thousands of strongly lensed high-z galaxies distributed over the full sky. The extreme brightness of these galaxies will make it possible to study them, via follow-up observations, in extraordinary detail. Also, the CORE resolution matches the typical sizes of high-z galaxy proto-clusters much better than the Planck resolution, resulting in a much higher detection efficiency; these objects will be caught in an evolutionary phase beyond the reach of surveys in other wavebands. Furthermore, CORE will provide unique information on the evolution of the star formation in virialized groups and clusters of galaxies up to the highest possible redshifts. Finally, thanks to its very high sensitivity, CORE will detect the polarized emission of thousands of radio sources and, for the first time, of dusty galaxies, at mm and sub-mm wavelengths, respectively.

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.

The interplay between cosmology and earth based experiments is crucial in order to pin down neutrino physics. Indeed cosmology can provide very tight, yet model dependent, constraints on some neutrino properties. Here we focus on the neutrino mass sum, reviewing the up to date current bounds and showing the results of our forecast of the sensitivity of future experiments. Finally, we discuss the case for sterile neutrinos, explaining how non standard sterile neutrino self-interactions can reconcile the oscillation anomalies with cosmology.

We forecast the scientific capabilities to improve our understanding of cosmic inflation of CORE, a proposed CMB space satellite submitted in response to the ESA fifth call for a medium-size mission opportunity. The CORE satellite will map the CMB anisotropies in temperature and polarization in 19 frequency channels spanning the range 60-600 GHz. CORE will have an aggregate noise sensitivity of $1.7 \mu$K$\cdot \,$arcmin and an angular resolution of 5' at 200 GHz. We explore the impact of telescope size and noise sensitivity on the inflation science return by making forecasts for several instrumental configurations. This study assumes that the lower and higher frequency channels suffice to remove foreground contaminations and complements other related studies of component separation and systematic effects, which will be reported in other papers of the series "Exploring Cosmic Origins with CORE." We forecast the capability to determine key inflationary parameters, to lower the detection limit for the tensor-to-scalar ratio down to the $10^{-3}$ level, to chart the landscape of single field slow-roll inflationary models, to constrain the epoch of reheating, thus connecting inflation to the standard radiation-matter dominated Big Bang era, to reconstruct the primordial power spectrum, to constrain the contribution from isocurvature perturbations to the $10^{-3}$ level, to improve constraints on the cosmic string tension to a level below the presumptive GUT scale, and to improve the current measurements of primordial non-Gaussianities down to the $f_{NL}^{\rm local} < 1$ level. For all the models explored, CORE alone will improve significantly on the present constraints on the physics of inflation. Its capabilities will be further enhanced by combining with complementary future cosmological observations.

We forecast the main cosmological parameter constraints achievable with the CORE space mission which is dedicated to mapping the polarisation of the Cosmic Microwave Background (CMB). CORE was recently submitted in response to ESA's fifth call for medium-sized mission proposals (M5). Here we report the results from our pre-submission study of the impact of various instrumental options, in particular the telescope size and sensitivity level, and review the great, transformative potential of the mission as proposed. Specifically, we assess the impact on a broad range of fundamental parameters of our Universe as a function of the expected CMB characteristics, with other papers in the series focusing on controlling astrophysical and instrumental residual systematics. In this paper, we assume that only a few central CORE frequency channels are usable for our purpose, all others being devoted to the cleaning of astrophysical contaminants. On the theoretical side, we assume LCDM as our general framework and quantify the improvement provided by CORE over the current constraints from the Planck 2015 release. We also study the joint sensitivity of CORE and of future Baryon Acoustic Oscillation and Large Scale Structure experiments like DESI and Euclid. Specific constraints on the physics of inflation are presented in another paper of the series. In addition to the six parameters of the base LCDM, which describe the matter content of a spatially flat universe with adiabatic and scalar primordial fluctuations from inflation, we derive the precision achievable on parameters like those describing curvature, neutrino physics, extra light relics, primordial helium abundance, dark matter annihilation, recombination physics, variation of fundamental constants, dark energy, modified gravity, reionization and cosmic birefringence. (ABRIDGED)

We examine the cosmological constraints that can be achieved with a galaxy cluster survey with the future CORE space mission. Using realistic simulations of the millimeter sky, produced with the latest version of the Planck Sky Model, we characterize the CORE cluster catalogues as a function of the main mission performance parameters. We pay particular attention to telescope size, key to improved angular resolution, and discuss the comparison and the complementarity of CORE with ambitious future ground-based CMB experiments that could be deployed in the next decade. A possible CORE mission concept with a 150 cm diameter primary mirror can detect of the order of 50,000 clusters through the thermal Sunyaev-Zeldovich effect (SZE). The total yield increases (decreases) by 25% when increasing (decreasing) the mirror diameter by 30 cm. The 150 cm telescope configuration will detect the most massive clusters ($>10^{14}\, M_\odot$) at redshift $z>1.5$ over the whole sky, although the exact number above this redshift is tied to the uncertain evolution of the cluster SZE flux-mass relation; assuming self-similar evolution, CORE will detect $\sim 500$ clusters at redshift $z>1.5$. This changes to 800 (200) when increasing (decreasing) the mirror size by 30 cm. CORE will be able to measure individual cluster halo masses through lensing of the cosmic microwave background anisotropies with a 1-$\sigma$ sensitivity of $4\times10^{14} M_\odot$, for a 120 cm aperture telescope, and $10^{14} M_\odot$ for a 180 cm one. [abridged]

We compute cosmic microwave background (CMB) anisotropy constraints on exotic forms of energy injection in electromagnetic (e.m.) channels over a large range of timescales. These constraints are very powerful around or just after recombination, although CMB keeps some sensitivity e.g. to decaying species with lifetimes as long as $10^{25}\,$s. We review here complementary with CMB spectral distortions and primordial nucleosynthesis bounds, which dominate at earlier timescales. For the first time, we describe the effects of the e.m. energy injection on the CMB power spectra as a function of the injection epoch, using the lifetime of a decaying particle as proxy. We identify a suitable on-the-spot approximation. Our results are of interest not only for early universe relics constituting (a fraction of) the dark matter, but also for other exotic injection of e.m. radiation. For illustration, we apply our formalism to: i) Primordial black holes of mass $\in [10^{13.5},10^{16.8}]$ g, showing that the constraints are comparable to the ones obtained from gamma-ray background studies and even dominate below $\sim 10^{14}$g. ii) To a peculiar mass-mixing range in the sterile neutrino parameter space, complementary to other astrophysical and laboratory probes. iii) Finally, we provide a first estimate of the room for improvement left for forthcoming 21 cm experiments, comparing it with the reach of proposed CMB spectral distortion (PiXiE) and CMB angular power spectrum (CORE) missions. We show that the best and most realistic opportunity to look for this signal (or to improve over current constraints) in the 21 cm probe is to focus on the Cosmic Dawn epoch, $15\lesssim z\lesssim30$, where the qualitatively unambiguous signature of a spectrum in emission can be expected for models that evade all current constraints.

We present a comprehensive review of keV-scale sterile neutrino Dark Matter, collecting views and insights from all disciplines involved - cosmology, astrophysics, nuclear, and particle physics - in each case viewed from both theoretical and experimental/observational perspectives. After reviewing the role of active neutrinos in particle physics, astrophysics, and cosmology, we focus on sterile neutrinos in the context of the Dark Matter puzzle. Here, we first review the physics motivation for sterile neutrino Dark Matter, based on challenges and tensions in purely cold Dark Matter scenarios. We then round out the discussion by critically summarizing all known constraints on sterile neutrino Dark Matter arising from astrophysical observations, laboratory experiments, and theoretical considerations. In this context, we provide a balanced discourse on the possibly positive signal from X-ray observations. Another focus of the paper concerns the construction of particle physics models, aiming to explain how sterile neutrinos of keV-scale masses could arise in concrete settings beyond the Standard Model of elementary particle physics. The paper ends with an extensive review of current and future astrophysical and laboratory searches, highlighting new ideas and their experimental challenges, as well as future perspectives for the discovery of sterile neutrinos.

Future Cosmic Microwave Background experiments together with upcoming galaxy and 21-cm surveys will provide extremely accurate measurements of different cosmological observables located at different epochs of the cosmic history. The new data will be able to constrain the neutrino mass sum with the best precision ever. In order to exploit the complementarity of the different redshift probes, a deep understanding of the physical effects driving the impact of massive neutrinos on CMB and large scale structures is required. The goal of this work is to describe these effects, assuming a summed neutrino mass close to its minimum allowed value. We find that parameter degeneracies can be removed by appropriate combinations, leading to robust and model independent constraints. A joint forecast of the sensitivity of Euclid and DESI surveys together with a CORE-like CMB experiment leads to a $1\sigma$ uncertainty of $14$~meV on the summed neutrino mass. However this particular combination gives rise to a peculiar degeneracy between $M_\nu$ and the optical depth at reionization. Independent constraints from 21-cm surveys can break this degeneracy and decrease the $1\sigma$ uncertainty down to $12$~meV.

This paper presents a study of the ISW effect from the Planck 2015 temperature and polarization data release. The CMB is cross-correlated with different LSS tracers: the NVSS, SDSS and WISE catalogues, and the Planck 2015 lensing map. This cross-correlation yields a detection at $4\,\sigma$, where most of the signal-to-noise is due to the Planck lensing and NVSS. In fact, the ISW effect is detected only from the Planck data (through the ISW-lensing bispectrum) at $\approx 3\,\sigma$, which is similar to the detection level achieved by combining the cross-correlation signal coming from all the catalogues. The ISW signal allow us to detect $\Omega_\Lambda$ at more than $3\,\sigma$. This cross-correlation analysis is performed only with the Planck temperature data, since the polarization scales available in the 2015 release do not permit significant improvement of the CMB-LSS cross-correlation detectability. Nevertheless, polarization data is used to study the anomalously large ISW signal previously reported through the stacking of CMB features at the locations of known superstructures. We find that the current Planck polarization data do not exclude that this signal could be caused by the ISW effect. In addition, the stacking of the Planck lensing map on the locations of superstructures exhibits a positive cross-correlation with these large-scale structures. Finally, we have improved our previous reconstruction of the ISW temperature fluctuations by combining the information encoded in all the previously mentioned LSS tracers. In particular, we construct a map of the ISW secondary anisotropies and the corresponding uncertainties map, obtained from simulations. We also explore the reconstruction of the ISW anisotropies caused by the LSS traced by the 2MPZ survey by directly inverting the density field into the gravitational potential field.

We present the most significant measurement of the cosmic microwave background (CMB) lensing potential to date (at a level of 40 sigma), using temperature and polarization data from the Planck 2015 full-mission release. Using a polarization-only estimator we detect lensing at a significance of 5 sigma. We cross-check the accuracy of our measurement using the wide frequency coverage and complementarity of the temperature and polarization measurements. Public products based on this measurement include an estimate of the lensing potential over approximately 70% of the sky, an estimate of the lensing potential power spectrum in bandpowers for the multipole range 40<L<400 and an associated likelihood for cosmological parameter constraints. We find good agreement between our measurement of the lensing potential power spectrum and that found in the best-fitting LCDM model based on the Planck temperature and polarization power spectra. Using the lensing likelihood alone we obtain a percent-level measurement of the parameter combination $\sigma_8 \Omega_m^{0.25} = 0.591\pm 0.021$. We combine our determination of the lensing potential with the E-mode polarization also measured by Planck to generate an estimate of the lensing B-mode. We show that this lensing B-mode estimate is correlated with the B-modes observed directly by Planck at the expected level and with a statistical significance of 10 sigma, confirming Planck's sensitivity to this known sky signal. We also correlate our lensing potential estimate with the large-scale temperature anisotropies, detecting a cross-correlation at the 3 sigma level, as expected due to dark energy in the concordance LCDM model.

The Planck full mission cosmic microwave background(CMB) temperature and E-mode polarization maps are analysed to obtain constraints on primordial non-Gaussianity(NG). Using three classes of optimal bispectrum estimators - separable template-fitting (KSW), binned, and modal - we obtain consistent values for the local, equilateral, and orthogonal bispectrum amplitudes, quoting as our final result from temperature alone fNL^local=2.5+\-5.7, fNL^equil=-16+\-70 and fNL^ortho=-34+\-33(68%CL). Combining temperature and polarization data we obtain fNL^local=0.8+\-5.0, fNL^equil=-4+\-43 and fNL^ortho=-26+\-21 (68%CL). The results are based on cross-validation of these estimators on simulations, are stable across component separation techniques, pass an extensive suite of tests, and are consistent with Minkowski functionals based measurements. The effect of time-domain de-glitching systematics on the bispectrum is negligible. In spite of these test outcomes we conservatively label the results including polarization data as preliminary, owing to a known mismatch of the noise model in simulations and the data. Beyond fNL estimates, we present model-independent reconstructions of the CMB bispectrum and derive constraints on early universe scenarios that generate NG, including general single-field and axion inflation, initial state modifications, parity-violating tensor bispectra, and directionally dependent vector models. We also present a wide survey of scale-dependent oscillatory bispectra, and we look for isocurvature NG. Our constraint on the local primordial trispectrum amplitude is gNL^local=(-9.0+\-7.7)x10^4 (68%CL), and we perform an analysis of additional trispectrum shapes. The global picture is one of consistency with the premises of the LambdaCDM cosmology, namely that the structure we observe today was sourced by adiabatic, passive, Gaussian, and primordial seed perturbations.[abridged]

This paper presents the Planck 2015 likelihoods, statistical descriptions of the 2-point correlations of CMB data, using the hybrid approach employed previously: pixel-based at $\ell<30$ and a Gaussian approximation to the distribution of spectra at higher $\ell$. The main improvements are the use of more and better processed data and of Planck polarization data, and more detailed foreground and instrumental models, allowing further checks and enhanced immunity to systematics. Progress in foreground modelling enables a larger sky fraction. Improvements in processing and instrumental models further reduce uncertainties. For temperature, we perform an analysis of end-to-end instrumental simulations fed into the data processing pipeline; this does not reveal biases from residual instrumental systematics. The $\Lambda$CDM cosmological model continues to offer a very good fit to Planck data. The slope of primordial scalar fluctuations, $n_s$, is confirmed smaller than unity at more than 5{\sigma} from Planck alone. We further validate robustness against specific extensions to the baseline cosmology. E.g., the effective number of neutrino species remains compatible with the canonical value of 3.046. This first detailed analysis of Planck polarization concentrates on E modes. At low $\ell$ we use temperature at all frequencies and a subset of polarization. The frequency range improves CMB-foreground separation. Within the baseline model this requires a reionization optical depth $\tau=0.078\pm0.019$, significantly lower than without high-frequency data for explicit dust monitoring. At high $\ell$ we detect residual errors in E, typically O($\mu$K$^2$); we recommend temperature alone as the high-$\ell$ baseline. Nevertheless, Planck high-$\ell$ polarization allows a separate determination of $\Lambda$CDM parameters consistent with those from temperature alone.

We present results based on full-mission Planck observations of temperature and polarization anisotropies of the CMB. These data are consistent with the six-parameter inflationary LCDM cosmology. From the Planck temperature and lensing data, for this cosmology we find a Hubble constant, H0= (67.8 +/- 0.9) km/s/Mpc, a matter density parameter Omega_m = 0.308 +/- 0.012 and a scalar spectral index with n_s = 0.968 +/- 0.006. (We quote 68% errors on measured parameters and 95% limits on other parameters.) Combined with Planck temperature and lensing data, Planck LFI polarization measurements lead to a reionization optical depth of tau = 0.066 +/- 0.016. Combining Planck with other astrophysical data we find N_ eff = 3.15 +/- 0.23 for the effective number of relativistic degrees of freedom and the sum of neutrino masses is constrained to < 0.23 eV. Spatial curvature is found to be |Omega_K| < 0.005. For LCDM we find a limit on the tensor-to-scalar ratio of r <0.11 consistent with the B-mode constraints from an analysis of BICEP2, Keck Array, and Planck (BKP) data. Adding the BKP data leads to a tighter constraint of r < 0.09. We find no evidence for isocurvature perturbations or cosmic defects. The equation of state of dark energy is constrained to w = -1.006 +/- 0.045. Standard big bang nucleosynthesis predictions for the Planck LCDM cosmology are in excellent agreement with observations. We investigate annihilating dark matter and deviations from standard recombination, finding no evidence for new physics. The Planck results for base LCDM are in agreement with BAO data and with the JLA SNe sample. However the amplitude of the fluctuations is found to be higher than inferred from rich cluster counts and weak gravitational lensing. Apart from these tensions, the base LCDM cosmology provides an excellent description of the Planck CMB observations and many other astrophysical data sets.

(abridged) We discuss the Galactic foreground emission between 20 and 100GHz based on observations by Planck/WMAP. The Commander component-separation tool has been used to separate the various astrophysical processes in total intensity. Comparison with RRL templates verifies the recovery of the free-free emission along the Galactic plane. Comparison of the high-latitude Halpha emission with our free-free map shows residuals that correlate with dust optical depth, consistent with a fraction (~30%) of Halpha having been scattered by high-latitude dust. We highlight a number of diffuse spinning dust morphological features at high latitude. There is substantial spatial variation in the spinning dust spectrum, with the emission peak ranging from below 20GHz to more than 50GHz. There is a strong tendency for the spinning dust component near many prominent HII regions to have a higher peak frequency, suggesting that this increase in peak frequency is associated with dust in the photodissociation regions around the nebulae. The emissivity of spinning dust in these diffuse regions is of the same order as previous detections in the literature. Over the entire sky, the commander solution finds more anomalous microwave emission than the WMAP component maps, at the expense of synchrotron and free-free emission. This can be explained by the difficulty in separating multiple broadband components with a limited number of frequency maps. Future surveys (5-20GHz), will greatly improve the separation by constraining the synchrotron spectrum. We combine Planck/WMAP data to make the highest S/N ratio maps yet of the intensity of the all-sky polarized synchrotron emission at frequencies above a few GHz. Most of the high-latitude polarized emission is associated with distinct large-scale loops and spurs, and we re-discuss their structure...

We consider a cosmological model in which a fraction $f$ of the Dark Matter (DM) is allowed to decay in an invisible relativistic component, and compute the resulting constraints on both the decay width (or inverse lifetime) $\Gamma$ and $f$ from purely gravitational arguments. We report a full derivation of the Boltzmann hierarchy, correcting a mistake in previous literature, and compute the impact of the decay --as a function of the lifetime-- on the CMB and matter power spectra. From CMB only, we obtain that no more than 3.8 % of the DM could have decayed in the time between recombination and today (all bounds quoted at 95 % CL). We also comment on the important application of this bound to the case where primordial black holes constitute DM, a scenario notoriously difficult to constrain. For lifetimes longer than the age of the Universe, the bounds can be cast as $f\Gamma < 6.3\times10^{-3}$ Gyr$^{-1}$. For the first time, we also checked that degeneracies with massive neutrinos are broken when information from the large scale structure is used. Even secondary effects like CMB lensing suffice to this purpose. Decaying DM models have been invoked to solve a possible tension between low redshift astronomical measurements of $\sigma_8$ and $\Omega_{\rm m}$ and the ones inferred by Planck. We reassess this claim finding that with the most recent BAO, HST and $\sigma_8$ data extracted from the CFHT survey, the tension is only slightly reduced despite the two additional free parameters, loosening the bound to $f\Gamma < 15.9\times10^{-3}$ Gyr$^{-1}$. The bound however improves to $f\Gamma < 5.9\times10^{-3}$ Gyr$^{-1}$ if only data consistent with the CMB are included. This highlights the importance of establishing whether the tension is due to real physical effects or unaccounted systematics, for settling the reach of achievable constraints on decaying DM.

Full-sky CMB maps from the 2015 Planck release allow us to detect departures from global isotropy on the largest scales. We present the first searches using CMB polarization for correlations induced by a non-trivial topology with a fundamental domain intersecting, or nearly intersecting, the last scattering surface (at comoving distance $\chi_{rec}$). We specialize to flat spaces with toroidal and slab topologies, finding that explicit searches for the latter are sensitive to other topologies with antipodal symmetry. These searches yield no detection of a compact topology at a scale below the diameter of the last scattering surface. The limits on the radius $R_i$ of the largest sphere inscribed in the topological domain (at log-likelihood-ratio $\Delta\ln{L}>-5$ relative to a simply-connected flat Planck best-fit model) are $R_i>0.97\chi_{rec}$ for the cubic torus and $R_i>0.56\chi_{rec}$ for the slab. The limit for the cubic torus from the matched-circles search is numerically equivalent, $R_i>0.97\chi_{rec}$ (99% CL) from polarisation data alone. We also perform a Bayesian search for a Bianchi VII$_h$ geometry. In the non-physical setting where the Bianchi cosmology is decoupled from the standard cosmology, Planck temperature data favour the inclusion of a Bianchi component. However, the cosmological parameters generating this pattern are in strong disagreement with those found from CMB anisotropy data alone. Fitting the induced polarization pattern for this model to Planck data requires an amplitude of $-0.1\pm0.04$ compared to +1 if the model were to be correct. In the physical setting where the Bianchi parameters are fit simultaneously with the standard cosmological parameters, we find no evidence for a Bianchi VII$_h$ cosmology and constrain the vorticity of such models to $(\omega/H)_0<7.6\times10^{-10}$ (95% CL). [Abridged]

We modified the CLASS code in order to include relativistic galaxy number counts in spatially curved geometries; we present the formalism and study the effect of relativistic corrections on spatial curvature. The new version of the code is now publicly available. Using a Fisher matrix analysis, we investigate how measurements of the spatial curvature parameter $\Omega_K$ with future galaxy surveys are affected by relativistic effects, which influence observations of the large scale galaxy distribution. These effects include contributions from cosmic magnification, Doppler terms and terms involving the gravitational potential. As an application, we consider angle and redshift dependent power spectra, which are especially well suited for model independent cosmological constraints. We compute our results for a representative deep, wide and spectroscopic survey, and our results show the impact of relativistic corrections on the spatial curvature parameter estimation. We show that constraints on the curvature parameter may be strongly biased if, in particular, cosmic magnification is not included in the analysis. Other relativistic effects turn out to be subdominant in the studied configuration. We analyze how the shift in the estimated best-fit value for the curvature and other cosmological parameters depends on the magnification bias parameter, and find that significant biases are to be expected if this term is not properly considered in the analysis.

We study the implications of Planck data for models of dark energy (DE) and modified gravity (MG), beyond the cosmological constant scenario. We start with cases where the DE only directly affects the background evolution, considering Taylor expansions of the equation of state, principal component analysis and parameterizations related to the potential of a minimally coupled DE scalar field. When estimating the density of DE at early times, we significantly improve present constraints. We then move to general parameterizations of the DE or MG perturbations that encompass both effective field theories and the phenomenology of gravitational potentials in MG models. Lastly, we test a range of specific models, such as k-essence, f(R) theories and coupled DE. In addition to the latest Planck data, for our main analyses we use baryonic acoustic oscillations, type-Ia supernovae and local measurements of the Hubble constant. We further show the impact of measurements of the cosmological perturbations, such as redshift-space distortions and weak gravitational lensing. These additional probes are important tools for testing MG models and for breaking degeneracies that are still present in the combination of Planck and background data sets. All results that include only background parameterizations are in agreement with LCDM. When testing models that also change perturbations (even when the background is fixed to LCDM), some tensions appear in a few scenarios: the maximum one found is \sim 2 sigma for Planck TT+lowP when parameterizing observables related to the gravitational potentials with a chosen time dependence; the tension increases to at most 3 sigma when external data sets are included. It however disappears when including CMB lensing.

We estimate the sensitivity of future galaxy surveys to cosmological parameters, using the redshift dependent angular power spectra of galaxy number counts, $C_\ell(z_1,z_2)$, calculated with all relativistic corrections at first order in perturbation theory. We pay special attention to the redshift dependence of the non-linearity scale and present Fisher matrix forecasts for Euclid-like and DES-like galaxy surveys. We compare the standard $P(k)$ analysis with the new $C_\ell(z_1,z_2)$ method. We show that for surveys with photometric redshifts the new analysis performs significantly better than the $P(k)$ analysis. For spectroscopic redshifts, however, the large number of redshift bins which would be needed to fully profit from the redshift information, is severely limited by shot noise. We also identify surveys which can measure the lensing contribution and we study the monopole, $C_0(z_1,z_2)$.

We present some accurate and efficient computations of large scale structure observables, obtained with a modified version of the CLASS code which is made publicly available. This code includes all relativistic corrections and computes both the power spectrum Cl(z1,z2) and the corresponding correlation function xi(theta,z1,z2) in linear perturbation theory. For Gaussian initial perturbations, these quantities contain the full information encoded in the large scale matter distribution at the level of linear perturbation theory. We illustrate the usefulness of our code for cosmological parameter estimation through a few simple examples.

We study a two-parameter extension of the cosmological standard model $\Lambda$CDM in which cold dark matter interacts with a new form of dark radiation. The two parameters correspond to the energy density in the dark radiation fluid $\Delta N_\mathrm{fluid}$ and the interaction strength between dark matter and dark radiation. The interactions give rise to a very weak "dark matter drag" which damps the growth of matter density perturbations throughout radiation domination, allowing to reconcile the tension between predictions of large scale structure from the CMB and direct measurements of $\sigma_8$. We perform a precision fit to Planck CMB data, BAO, large scale structure, and direct measurements of the expansion rate of the universe today. Our model lowers the $\chi$-squared relative to $\Lambda$CDM by about 12, corresponding to a preference for non-zero dark matter drag by more than $3 \sigma$. Particle physics models which naturally produce a dark matter drag of the required form include the recently proposed non-Abelian dark matter model in which the dark radiation corresponds to massless dark gluons.

We compute and investigate four types of imprint of a stochastic background of primordial magnetic fields (PMFs) on the cosmic microwave background (CMB) anisotropies: the impact of PMFs on the CMB spectra; the effect on CMB polarization induced by Faraday rotation; the impact of PMFs on the ionization history; magnetically-induced non-Gaussianities; and the magnetically-induced breaking of statistical isotropy. Overall, Planck data constrain the amplitude of PMFs to less than a few nanogauss. In particular, individual limits coming from the analysis of the CMB angular power spectra, using the Planck likelihood, are $B_{1\,\mathrm{Mpc}}< 4.4$ nG (where $B_{1\,\mathrm{Mpc}}$ is the comoving field amplitude at a scale of 1 Mpc) at 95% confidence level, assuming zero helicity, and $B_{1\,\mathrm{Mpc}}< 5.6$ nG for a maximally helical field.For nearly scale-invariant PMFs we obtain $B_{1\,\mathrm{Mpc}}<2.0$ nG and $B_{1\,\mathrm{Mpc}}<0.9$ nG if the impact of PMFs on the ionization history of the Universe is included. From the analysis of magnetically-induced non-Gaussianity we obtain three different values, corresponding to three applied methods, all below 5 nG. The constraint from the magnetically-induced passive-tensor bispectrum is $B_{1\,\mathrm{Mpc}}< 2.8$ nG. A search for preferred directions in the magnetically-induced passive bispectrum yields $B_{1\,\mathrm{Mpc}}< 4.5$ nG, whereas the the compensated-scalar bispectrum gives $B_{1\,\mathrm{Mpc}}< 3$ nG. The analysis of the Faraday rotation of CMB polarization by PMFs uses the Planck power spectra in $EE$ and $BB$ at 70 GHz and gives $B_{1\,\mathrm{Mpc}}< 1380$ nG. In our final analysis, we consider the harmonic-space correlations produced by Alfv\'en waves, finding no significant evidence for the presence of these waves. Together, these results comprise a comprehensive set of constraints on possible PMFs with Planck data.

We present an updated description of the Planck Low Frequency Instrument (LFI) data processing pipeline, associated with the 2015 data release. We point out the places where our results and methods have remained unchanged since the 2013 paper and we highlight the changes made for the 2015 release, describing the products (especially timelines) and the ways in which they were obtained. We demonstrate that the pipeline is self-consistent (principally based on simulations) and report all null tests. For the first time, we present LFI maps in Stokes Q and U polarization. We refer to other related papers where more detailed descriptions of the LFI data processing pipeline may be found if needed.

The Second Planck Catalogue of Compact Sources is a catalogue of sources detected in single-frequency maps from the full duration of the Planck mission and supersedes previous versions of the Planck compact source catalogues. It consists of compact sources, both Galactic and extragalactic, detected over the entire sky. Compact sources detected in the lower frequency channels are assigned to the PCCS2, while at higher frequencies they are assigned to one of two sub-catalogues, the PCCS2 or PCCS2E, depending on their location on the sky. The first of these catalogues covers most of the sky and allows the user to produce subsamples at higher reliabilities than the target 80% integral reliability of the catalogue. The PCCS2E contains sources detected in sky regions where the diffuse emission makes it difficult to quantify the reliability of the detections. Both the PCCS2 and PCCS2E include polarization measurements, in the form of polarized flux densities, or upper limits, and orientation angles for all seven polarization-sensitive Planck channels. The improved data-processing of the full-mission maps and their reduced noise levels allow us to increase the number of objects in the catalogue, improving its completeness for the target 80 % reliability as compared with the previous versions, the PCCS and ERCSC catalogues.

We test the statistical isotropy and Gaussianity of the cosmic microwave background (CMB) anisotropies using observations made by the Planck satellite. Our results are based mainly on the full Planck mission for temperature, but also include some polarization measurements. In particular, we consider the CMB anisotropy maps derived from the multi-frequency Planck data by several component-separation methods. For the temperature anisotropies, we find excellent agreement between results based on these sky maps over both a very large fraction of the sky and a broad range of angular scales, establishing that potential foreground residuals do not affect our studies. Tests of skewness, kurtosis, multi-normality, N-point functions, and Minkowski functionals indicate consistency with Gaussianity, while a power deficit at large angular scales is manifested in several ways, for example low map variance. The results of a peak statistics analysis are consistent with the expectations of a Gaussian random field. The "Cold Spot" is detected with several methods, including map kurtosis, peak statistics, and mean temperature profile. We thoroughly probe the large-scale dipolar power asymmetry, detecting it with several independent tests, and address the subject of a posteriori correction. Tests of directionality suggest the presence of angular clustering from large to small scales, but at a significance that is dependent on the details of the approach. We perform the first examination of polarization data, finding the morphology of stacked peaks to be consistent with the expectations of statistically isotropic simulations. Where they overlap, these results are consistent with the Planck 2013 analysis based on the nominal mission data and provide our most thorough view of the statistics of the CMB fluctuations to date.

It is well known that annihilations in the homogeneous fluid of dark matter (DM) can leave imprints in the cosmic microwave background (CMB) anisotropy power spectrum. However, the relevance of DM annihilations in halos for cosmological observables is still subject to debate, with previous works reaching different conclusions on this point. Also, all previous studies used a single type of parameterization for the astrophysical reionization, and included no astrophysical source for the heating of the intergalactic medium. In this work, we revisit these problems. When standard approaches are adopted, we find that the ionization fraction does exhibit a very particular (and potentially constraining) pattern, but the currently measurable optical depth to reionization is left almost unchanged: In agreement with the most of the previous literature, for plausible halo models we find that the modification of the signal with respect to the one coming from annihilations in the smooth background is tiny, below cosmic variance within currently allowed parameter space. However, if different and probably more realistic treatments of the astrophysical sources of reionization and heating are adopted, a more pronounced effect of the DM annihilation in halos is possible. We thus conclude that within currently adopted baseline models the impact of the virialised DM structures cannot be uncovered by CMB power spectra measurements, but a larger impact is possible if peculiar models are invoked for the redshift evolution of the DM annihilation signal or different assumptions are made for the astrophysical contributions. A better understanding (both theoretical and observational) of the reionization and temperature history of the universe, notably via the 21 cm signal, seems the most promising way for using halo formation as a tool in DM searches, improving over the sensitivity of current cosmological probes.

We present a description of the pipeline used to calibrate the Planck Low Frequency Instrument (LFI) timelines into thermodynamic temperatures for the Planck 2015 data release, covering four years of uninterrupted operations. As in the 2013 data release, our calibrator is provided by the spin-synchronous modulation of the cosmic microwave background dipole, but we now use the orbital component, rather than adopting the Wilkinson Microwave Anisotropy Probe (WMAP) solar dipole. This allows our 2015 LFI analysis to provide an independent Solar dipole estimate, which is in excellent agreement with that of HFI and within $1\sigma$ (0.3% in amplitude) of the WMAP value. This 0.3% shift in the peak-to-peak dipole temperature from WMAP and a global overhaul of the iterative calibration code increases the overall level of the LFI maps by 0.45% (30 GHz), 0.64% (44 GHz), and 0.82% (70 GHz) in temperature with respect to the 2013 Planck data release, thus reducing the discrepancy with the power spectrum measured by WMAP. We estimate that the LFI calibration uncertainty is now at the level of 0.20% for the 70 GHz map, 0.26% for the 44 GHz map, and 0.35% for the 30 GHz map. We provide a detailed description of the impact of all the changes implemented in the calibration since the previous data release.

We present a new flexible, fast and accurate way to implement massive neutrinos, warm dark matter and any other non-cold dark matter relics in Boltzmann codes. For whatever analytical or numerical form of the phase-space distribution function, the optimal sampling in momentum space compatible with a given level of accuracy is automatically found by comparing quadrature methods. The perturbation integration is made even faster by switching to an approximate viscous fluid description inside the Hubble radius, which differs from previous approximations discussed in the literature. When adding one massive neutrino to the minimal cosmological model, CLASS becomes just 1.5 times slower, instead of about 5 times in other codes (for fixed accuracy requirements). We illustrate the flexibility of our approach by considering a few examples of standard or non-standard neutrinos, as well as warm dark matter models.

This paper presents the characterization of the in-flight beams, the beam window functions, and the associated uncertainties for the Planck Low Frequency Instrument (LFI). The structure of the paper is similar to that presented in the 2013 Planck release; the main differences concern the beam normalization and the delivery of the window functions to be used for polarization analysis. The in-flight assessment of the LFI main beams relies on measurements performed during observations of Jupiter. By stacking data from seven Jupiter transits, the main beam profiles are measured down to -25 dB at 30 and 44 GHz, and down to -30 dB at 70 GHz. It has been confirmed that the agreement between the simulated beams and the measured beams is better than 1% at each LFI frequency band (within the 20 dB contour from the peak, the rms values are 0.1% at 30 and 70 GHz; 0.2% at 44 GHz). Simulated polarized beams are used for the computation of the effective beam window functions. The error budget for the window functions is estimated from both main beam and sidelobe contributions, and accounts for the radiometer band shapes. The total uncertainties in the effective beam window functions are 0.7% and 1% at 30 and 44 GHz, respectively (at $\ell \approx 600$); and 0.5% at 70 GHz (at $\ell \approx 1000$).

We present the 8th Full Focal Plane simulation set (FFP8), deployed in support of the Planck 2015 results. FFP8 consists of 10 fiducial mission realizations reduced to 18144 maps, together with the most massive suite of Monte Carlo realizations of instrument noise and CMB ever generated, comprising $10^4$ mission realizations reduced to about $10^6$ maps. The resulting maps incorporate the dominant instrumental, scanning, and data analysis effects; remaining subdominant effects will be included in future updates. Generated at a cost of some 25 million CPU-hours spread across multiple high-performance-computing (HPC) platforms, FFP8 is used for the validation and verification of analysis algorithms, as well as their implementations, and for removing biases from and quantifying uncertainties in the results of analyses of the real data.

We present joint constraints on the number of effective neutrino species N_eff and the sum of neutrino masses M_nu, based on a technique which exploits the full information contained in the one-dimensional Lyman-Alpha forest flux power spectrum, complemented by additional cosmological probes. In particular, we obtain N_eff=2.91(+0.21)(-0.22) (95% CL) and M_nu<0.15 eV (95% CL) when we combine BOSS Lyman-Alpha forest data with CMB (Planck+ACT+SPT+WMAP polarization) measurements, and N_eff=2.88(+0.20)(-0.20) (95% CL) and M_nu<0.14 eV (95% CL) when we further add baryon acoustic oscillations. Our results provide evidence for the Cosmic Neutrino Background from N_eff~3 (N_eff=0 is rejected at more than 14 sigma), and rule out the possibility of a sterile neutrino thermalized with active neutrinos (i.e., N_eff=4) - or more generally any decoupled relativistic relic with Delta N_eff ~ 1 - at a significance of over 5 sigma, the strongest bound to date, implying that there is no need for exotic neutrino physics in the concordance LCDM model.

The Planck design and scanning strategy provide many levels of redundancy that can be exploited to provide tests of internal consistency. One of the most important is the comparison of the 70GHz and 100GHz channels. Based on different instrument technologies, with feeds located differently in the focal plane, analysed independently by different teams using different software, and near the minimum of diffuse foreground emission, these channels are in effect two different experiments. The 143GHz channel has the lowest noise level on Planck, and is near the minimum of unresolved foreground emission. In this paper, we analyse the level of consistency achieved in the 2013 Planck data. We concentrate on comparisons between the 70/100/143GHz channel maps and power spectra, particularly over the angular scales of the first and second acoustic peaks, on maps masked for diffuse Galactic emission and for strong unresolved sources. Difference maps covering angular scales from 8deg-15arcmin are consistent with noise, and show no evidence of cosmic microwave background structure. Including small but important corrections for unresolved-source residuals, we demonstrate agreement between 70 and 100GHz power spectra averaged over 70<l<390 at the 0.8% level, and agreement between 143 and 100GHz power spectra of 0.4% over the same l range. These values are within and consistent with the overall uncertainties in calibration given in the Planck 2013 results. We also present results based on the 2013 likelihood analysis showing consistency at the 0.35% between the 100/143/217GHz power spectra. We analyse calibration procedures and beams to determine what fraction of these differences can be accounted for by known approximations or systematic errors that could be controlled even better in the future, reducing uncertainties still further. Several possible small improvements are described...(abridged)

The European Space Agency's Planck satellite, dedicated to studying the early Universe and its subsequent evolution, was launched 14~May 2009 and scanned the microwave and submillimetre sky continuously between 12~August 2009 and 23~October 2013. In February~2015, ESA and the Planck Collaboration released the second set of cosmology products based on data from the entire Planck mission, including both temperature and polarization, along with a set of scientific and technical papers and a web-based explanatory supplement. This paper gives an overview of the main characteristics of the data and the data products in the release, as well as the associated cosmological and astrophysical science results and papers. The science products include maps of the cosmic microwave background (CMB), the thermal Sunyaev-Zeldovich effect, and diffuse foregrounds in temperature and polarization, catalogues of compact Galactic and extragalactic sources (including separate catalogues of Sunyaev-Zeldovich clusters and Galactic cold clumps), and extensive simulations of signals and noise used in assessing the performance of the analysis methods and assessment of uncertainties. The likelihood code used to assess cosmological models against the Planck data are described, as well as a CMB lensing likelihood. Scientific results include cosmological parameters deriving from CMB power spectra, gravitational lensing, and cluster counts, as well as constraints on inflation, non-Gaussianity, primordial magnetic fields, dark energy, and modified gravity.