We investigate generalized interacting dark matter-dark energy scenarios with a time-dependent coupling parameter, allowing also for freedom in the neutrino sector. The models are tested in the phantom and quintessence regimes, characterized by an equation of state $w_x<-1$ and $w_x>-1$, respectively. Our analyses show that for some of the scenarios the existing tensions on the Hubble constant $H_0$ and on the clustering parameter $S_8$ can be significantly alleviated. The relief is either due to \textit{(a)} a dark energy component which lies within the phantom region; or \textit{(b)} the presence of a dynamical coupling in quintessence scenarios. The inclusion of massive neutrinos into the interaction schemes does not affect neither the constraints on the cosmological parameters nor the bounds on the total number or relativistic degrees of freedom $N_{\rm eff}$, which are found to be extremely robust and, in general, strongly consistent with the canonical prediction $N_{\rm eff}=3.045$. The most stringent bound on the total neutrino mass $M_{\nu}$ is $M_{\nu}<0.116$ eV and it is obtained within a quintessence scenario in which the matter mass-energy density is only mildly affected by the presence of a dynamical dark sector coupling.

We present the NPIPE processing pipeline, which produces calibrated frequency maps in temperature and polarization from data from the Planck Low Frequency Instrument (LFI) and High Frequency Instrument (HFI) using high-performance computers. NPIPE represents a natural evolution of previous Planck analysis efforts, and combines some of the most powerful features of the separate LFI and HFI analysis pipelines. The net effect of the improvements is lower levels of noise and systematics in both frequency and component maps at essentially all angular scales, as well as notably improved internal consistency between the various frequency channels. Based on the NPIPE maps, we present the first estimate of the Solar dipole determined through component separation across all nine Planck frequencies. The amplitude is ($3366.6 \pm 2.7$)$\mu$K, consistent with, albeit slightly higher than, earlier estimates. From the large-scale polarization data, we derive an updated estimate of the optical depth of reionization of $\tau = 0.051 \pm 0.006$, which appears robust with respect to data and sky cuts. There are 600 complete signal, noise and systematics simulations of the full-frequency and detector-set maps. As a Planck first, these simulations include full time-domain processing of the beam-convolved CMB anisotropies. The release of NPIPE maps and simulations is accompanied with a complete suite of raw and processed time-ordered data and the software, scripts, auxiliary data, and parameter files needed to improve further on the analysis and to run matching simulations.

The Phenomenologically Emergent Dark Energy model, a dark energy model with the same number of free parameters as the flat $\Lambda$CDM, has been proposed as a working example of a minimal model which can avoid the current cosmological tensions. A straightforward question is whether or not the inclusion of massive neutrinos and extra relativistic species may spoil such an appealing phenomenological alternative. We present the bounds on $M_{\nu}$ and $N_{\rm eff}$ and comment on the long standing $H_0$ and $\sigma_8$ tensions within this cosmological framework with a wealth of cosmological observations. Interestingly, we find, at $95\%$ confidence level, and with the most complete set of cosmological observations, $M_{\nu}\sim 0.21^{+0.15}_{-0.14}$ eV and $N_{\rm eff}= 3.03\pm 0.32$ i.e. an indication for a non-zero neutrino mass with a significance above $2\sigma$. The well known Hubble constant tension is considerably easened, with a significance always below the $2\sigma$ level.

We examine interactions between dark matter and dark energy in light of the latest cosmological observations, focusing on a specific model with coupling proportional to the dark energy density. Our data includes Cosmic Microwave Background (CMB) measurements from the Planck 2018 legacy data release, late-time measurements of the expansion history from Baryon Acoustic Oscillations (BAO) and Supernovae Type Ia (SNeIa), galaxy clustering and cosmic shear measurements from the Dark Energy Survey Year 1 results, and the 2019 local distance ladder measurement of the Hubble constant $H_0$ from the Hubble Space Telescope. Considering Planck data both in combination with BAO or SNeIa data reduces the $H_0$ tension to a level which could possibly be compatible with a statistical fluctuation. The very same model also significantly reduces the $\Omega_{\rm m}-\sigma_8$ tension between CMB and cosmic shear measurements. Interactions between the dark sectors of our Universe remain therefore a promising joint solution to these persisting cosmological tensions.

We do not solve tensions with concordance cosmology; we do obtain $H_0\approx 74\,$km/s/Mpc from CMB+BAO+SN data in our model, but that is not the point. Discrepancies in Hubble constant values obtained by various astrophysical probes should not be viewed in isolation. While one can resolve at least some of the differences through either an early time transition or late time transition in the expansion rate, these introduce other changes. We advocate a holistic approach, using a wide variety of cosmic data, rather than focusing on one number, $H_0$. Vacuum metamorphosis, a late time transition physically motivated by quantum gravitational effects and with the same number of parameters as \lcdm, can successfully give a high $H_0$ value from cosmic microwave background data but fails when combined with multiple distance probes. We also explore the influence of spatial curvature, and of a conjoined analysis of cosmic expansion and growth.

Type Ia Supernovae (SNeIa) used as standardizable candles have been instrumental in the discovery of cosmic acceleration, usually attributed to some form of dark energy (DE). Recent studies have raised the issue of whether intrinsic SNeIa luminosities might evolve with redshift. While the evidence for cosmic acceleration is robust to this possible systematic, the question remains of how much the latter can affect the inferred properties of the DE component responsible for cosmic acceleration. This is the question we address in this work. We use SNeIa distance moduli measurements from the Pantheon and JLA samples. We consider models where the DE equation of state is a free parameter, either constant or time-varying, as well as models where DE and dark matter interact, and finally a model-agnostic parametrization of effects due to modified gravity (MG). When SNeIa data are combined with Cosmic Microwave Background (CMB) temperature and polarization anisotropy measurements, we find strong degeneracies between parameters governing the SNeIa systematics, the DE parameters, and the Hubble constant $H_0$. These degeneracies significantly broaden the DE parameter uncertainties, in some cases leading to ${\cal O}(\sigma)$ shifts in the central values. However, including low-redshift Baryon Acoustic Oscillation and Cosmic Chronometer measurements, as well as CMB lensing measurements, considerably improves the previous constraints, and the only remaining effect of the examined systematic is a $\lesssim 40\%$ broadening of the uncertainties on the DE parameters. The constraints we derive on the MG parameters are instead basically unaffected by the systematic in question. We therefore confirm the overall soundness of dark energy properties.

We revisit our previous work [Phys. Rev. D 95, 096014 (2017)] where neutrino oscillation and nonoscillation data were analyzed in the standard framework with three neutrino families, in order to constrain their absolute masses and to probe their ordering (either normal, NO, or inverted, IO). We include updated oscillation results to discuss best fits and allowed ranges for the two squared mass differences $\delta m^2$ and $\Delta m^2$, the three mixing angles $\theta_{12}$, $\theta_{23}$ and $\theta_{13}$, as well as constraints on the CP-violating phase $\delta$, plus significant indications in favor of NO vs IO at the level of $\Delta\chi^2=10.0$. We then consider nonoscillation data from beta decay, from neutrinoless double beta decay (if neutrinos are Majorana), and from various cosmological input variants (in the data or the model) leading to results dubbed as default, aggressive, and conservative. In the default option, we obtain from nonoscillation data an extra contribution $\Delta\chi^2 = 2.2$ in favor of NO, and an upper bound on the sum of neutrino masses $\Sigma < 0.15$ eV at $2\sigma$; both results - dominated by cosmology - can be strengthened or weakened by using more aggressive or conservative options, respectively. Taking into account such variations, we find that the combination of all (oscillation and nonoscillation) neutrino data favors NO at the level of $3.2-3.7\sigma$, and that $\Sigma$ is constrained at the $2\sigma$ level within $\Sigma < 0.12-0.69$ eV. The upper edge of this allowed range corresponds to an effective $\beta$-decay neutrino mass $m_\beta = \Sigma/3 = 0.23$ eV, at the sensitivity frontier of the KATRIN experiment.

Consistent observations indicate that some of the important cosmological parameters measured through the local observations are in huge tension with their measurements from the global observations (within the minimal $\Lambda$CDM cosmology). The tensions in those cosmological parameters have been found to be either weakened or reconciled with the introduction of new degrees of freedom that effectively increases the underlying parameter space compared to the minimal $\Lambda$CDM cosmology. It might be interesting to investigate the above tensions within the context of an emergent dark energy scenario proposed recently by Li and Shafieloo~\cite{Li:2019yem}. We find that the tension on $H_0$ is clearly alleviated within 68\% confidence level with an improvement of the $\chi^2$ for CMB, for the above emergent dark energy model having only six free parameters similar to the spatially flat $\Lambda$CDM model. The tension on $H_0$ is still alleviated for every combined datasets considered in the work, however, such alleviation occurs by worsening the $\chi^2$ compared to the $\chi^2$ for $\Lambda$CDM model obtained for the same combined dataset.

We investigate the possibility of phantom crossing in the dark energy sector and solution for the Hubble tension between early and late Universe observations. We use robust combinations of different cosmological observations, namely the CMB, local measurement of Hubble constant ($H_0$), BAO and SnIa for this purpose. For a combination of CMB+BAO data which is related to early Universe physics, phantom crossing in the dark energy sector is confirmed at $95$\% confidence level and we obtain the constraint $H_0=71.0^{+2.9}_{-3.8}$ km/s/Mpc at 68\% confidence level which is in perfect agreement with the local measurement by Riess et al. We show that constraints from different combination of data are consistent with each other and all of them are consistent with phantom crossing in the dark energy sector. For the combination of all data considered, we obtain the constraint $H_0=70.25\pm 0.78$ km/s/Mpc at 68\% confidence level and the phantom crossing happening at the scale factor $a_m=0.851^{+0.048}_{-0.031}$ at 68\% confidence level.

The physics of the dark sector has remained one of the controversial areas of modern cosmology at present and hence it naturally attracts massive attention to the scientific community. With the developments of the astronomical data, the physics of the dark sector is becoming much more transparent than it was some twenty years back. The detection of gravitational waves (GWs) has now opened a cluster of possibilities in the cosmological regime. Being motivated by the detection of GWs and its possible impact on the physics of dark matter and dark energy, in this work we focus on the interacting dark energy models. Assuming the simplest possibility in which the vacuum energy with equation-of-state $w_x =-1$ is allowed to interact with the pressureless dark matter, we have extracted the constraints of the cosmological parameters. We find that the addition of the GWs data to the CMB measurements significantly improves up to a factor 4 of the parameters space, and up to a factor 2 for the full combination of current cosmological datasets, namely CMB+BAO+Pantheon+RSD+R16+CC+WL. The most affected parameters by the inclusion of the GWs are $\Omega_ch^2$, $\theta_{MC}$, $\xi$, and the derived parameters $\Omega_{m0}$, $\sigma_8$ and $H_0$.

In this paper we show how the string landscape can be constrained using observational data. We illustrate this idea by focusing on Fibre Inflation which is a promising class of string inflationary models in type IIB flux compactifications. We determine the values of the microscopic flux-dependent parameters which yield the best fit to the most recent cosmological datasets.

The largest temperature anisotropy in the cosmic microwave background (CMB) is the dipole, which has been measured with increasing accuracy for more than three decades, particularly with the Planck satellite. The simplest interpretation of the dipole is that it is due to our motion with respect to the rest frame of the CMB. Since current CMB experiments infer temperature anisotropies from angular intensity variations, the dipole modulates the temperature anisotropies with the same frequency dependence as the thermal Sunyaev-Zeldovich (tSZ) effect. We present the first (and significant) detection of this signal in the tSZ maps, and find that it is consistent with direct measurements of the CMB dipole, as expected. The signal contributes power in the tSZ maps, modulated in a quadrupolar pattern, and we estimate its contribution to the tSZ bispectrum, noting that it contributes negligible noise to the bispectrum at relevant scales.

In this paper we explore possible extensions of Interacting Dark Energy cosmologies, where Dark Energy and Dark Matter interact non-gravitationally with one another. In particular, we focus on the neutrino sector, analyzing the effect of both neutrino masses and the effective number of neutrino species. We consider the Planck 2018 legacy release data combined with several other cosmological probes, finding no evidence for new physics in the dark radiation sector. The current neutrino constraints from cosmology should be therefore regarded as robust, as they are not strongly dependent on the dark sector physics, once all the available observations are combined. Namely, we find a total neutrino mass $M_\nu<0.15$ eV and a number of effective relativistic degrees of freedom of $N_{\rm eff}=3.03^{+0.33}_{-0.33}$, both at 95\% CL, which are close to those obtained within the $\Lambda$CDM cosmology, $M_\nu<0.12$ eV and $N_{\rm eff}=3.00^{+0.36}_{-0.35}$ for the same data combination.

We show that a combined analysis of CMB anisotropy power spectra obtained by the Planck satellite and luminosity distance data simultaneously excludes a flat universe and a cosmological constant at $99 \%$ C.L. These results hold separately when combining Planck with three different datasets: the two determinations of the Hubble constant from Riess et al. 2019 and Freedman et al. 2020, and the Pantheon catalog of high redshift supernovae type-Ia. We conclude that either LCDM needs to be replaced by a drastically different model, or else there are significant but still undetected systematics. Our result calls for new observations and stimulates the investigation of alternative theoretical models and solutions.

We explore whether non-standard dark sector physics might be required to solve the existing cosmological tensions. The properties we consider in combination are an interaction between the dark matter and dark energy components, and a dark energy equation of state $w$ different from that of the canonical cosmological constant $w=-1$. In principle, these two parameters are independent. In practice, to avoid early-time, superhorizon instabilities, their allowed parameter spaces are correlated. We analyze three classes of extended interacting dark energy models in light of the 2019 Planck CMB results and Cepheid-calibrated local distance ladder $H_0$ measurements of Riess et al. (R19), as well as recent BAO and SNeIa distance data. We find that in quintessence coupled dark energy models, where $w > -1$, the evidence for a non-zero coupling between the two dark sectors can surpass the $5\sigma$ significance. On the other hand, in phantom coupled dark energy models, there is no such preference for a non-zero dark sector coupling. All the models we consider significantly raise the value of the Hubble constant easing the $H_0$ tension. The addition of low-redshift BAO and SNeIa measurements leaves some residual tension with R19 but at a level that could be justified by a statistical fluctuation. Bayesian evidence considerations mildly disfavour both the coupled quintessence and phantom models, while mildly favouring a coupled vacuum scenario, even when late-time datasets are considered. We conclude that non-minimal dark energy cosmologies, such as coupled quintessence, phantom, or vacuum models, are still an interesting route towards softening existing cosmological tensions, even when low-redshift datasets and Bayesian evidence considerations are taken into account. (abstract severely abridged)

We investigate the recently introduced metastable dark energy (DE) models after the final Planck 2018 legacy release. The essence of the present work is to analyze their evolution at the level of perturbations. Our analyses show that both the metastable dark energy models considered in this article, are excellent candidates to alleviate the $H_0$ tension. In particular, for the present models, Planck 2018 alone can alleviate the $H_0$ tension within 68\% CL. Along with the final cosmic microwave background data from the Planck 2018 legacy release, we also include external cosmological datasets in order to asses the robustness of our findings.

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.

We present new constraints on extended cosmological scenarios using the recent data from the Planck 2018 Legacy release. In addition to the six parameters of the standard LCDM model, we also simultaneously vary the dark energy equation of state, the neutrino mass, the neutrino effective number, the running of the spectral index and the lensing amplitude $A_L$. We confirm that a resolution of the Hubble tension is given by a dark energy equation of state with w<-1, ruling out quintessence models at high statistical significance. This solution is, however, not supported by BAO and Pantheon data. We find no evidence for evolving dark energy, i.e. $w_a \neq 0$. The neutrino effective number is always in agreement with the expectations of the standard model based on three active neutrinos. The running of the spectral index also is always consistent with zero. Despite the increase in the number of parameters, the $A_L$ lensing anomaly is still present at more than two standard deviations. The $A_L$ anomaly significantly affects the bounds on the neutrino mass that can be larger by a factor four with respect to those derived under standard $\Lambda$CDM. While the lensing data reduces the evidence for $A_L>1$, the inclusion of BAO and Pantheon increase its statistical significance.

While the origin and composition of dark matter and dark energy remains unknown, it is possible that they might represent two manifestations of a single entity, as occurring in unified dark sector models. On the other hand, advances in our understanding of the dark sector of the Universe might arise from Cosmic Dawn, the epoch when the first stars formed. In particular, the first detection of the global 21-cm absorption signal at Cosmic Dawn from the EDGES experiment opens up a new arena wherein to test models of dark matter and dark energy. Here, we consider generalized and modified Chaplygin gas models as candidate unified dark sector models. We first constrain these models against Cosmic Microwave Background data from the \textit{Planck} satellite, before exploring how the inclusion of the global 21-cm signal measured by EDGES can improve limits on the model parameters, finding that the uncertainties on the parameters of the Chaplygin gas models can be reduced by a factor between $1.5$ and $10$. We also find that within the generalized Chaplygin gas model, the tension between the CMB and local determinations of the Hubble constant $H_0$ is reduced from $\approx 4\sigma$ to $\approx 1.3\sigma$. In conclusion, we find that the global 21-cm signal at Cosmic Dawn can provide an extraordinary window onto the physics of unified dark sectors.

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

The recent Planck Legacy 2018 release has confirmed the presence of an enhanced lensing amplitude in CMB power spectra compared to that predicted in the standard $\Lambda$CDM model. A closed universe can provide a physical explanation for this effect, with the Planck CMB spectra now preferring a positive curvature at more than $99 \%$ C.L. Here we further investigate the evidence for a closed universe from Planck, showing that positive curvature naturally explains the anomalous lensing amplitude and demonstrating that it also removes a well-known tension within the Planck data set concerning the values of cosmological parameters derived at different angular scales. We show that since the Planck power spectra prefer a closed universe, discordances higher than generally estimated arise for most of the local cosmological observables, including BAO. The assumption of a flat universe could, therefore, mask a cosmological crisis where disparate observed properties of the Universe appear to be mutually inconsistent. Future measurements are needed to clarify whether the observed discordances are due to undetected systematics, or to new physics, or simply are a statistical fluctuation.

We investigate interacting scenarios which belong to a wider class, since they include a dynamical dark energy component whose equation of state follows various one-parameter parametrizations. We confront them with the latest observational data from Cosmic Microwave Background (CMB), Joint light-curve (JLA) sample from Supernovae Type Ia, Baryon Acoustic Oscillations (BAO), Hubble parameter measurements from Cosmic Chronometers (CC) and a gaussian prior on the Hubble parameter $H_0$. In all examined scenarios we find a non-zero interaction, nevertheless the non-interacting case is allowed within 2$\sigma$. Concerning the current value of the dark energy equation of state for all combination of datasets it always lies in the phantom regime at more than two/three standard deviations. Finally, for all interacting models, independently of the combination of datasets considered, the estimated values of the present Hubble parameter $H_0$ are greater compared to the $\Lambda$CDM-based Planck's estimation and close to the local measurements, thus alleviating the $H_0$ tension.

In a spatially flat Friedmann-Lema\^{i}tre-Robertson-Walker universe, we investigate a unified cosmic fluid scenario endowed with bulk viscosity in which the coefficient of the bulk viscosity has a power law evolution. The power law in the bulk viscous coefficient is a general case in this study which naturally includes several choices as special cases. Considering such a general bulk viscous scenario, in the present work we have extracted the observational constraints using the latest cosmological datasets and examine their behaviour at the level of both background and perturbations. From the observational analyses, we find that a non-zero bulk viscous coefficient is always favored and some of the models in this series are able to weaken the current tension on $H_0$ for some dataset. However, from the Bayesian evidence analysis, $\Lambda$CDM is favored over the bulk viscous model.

Coupled dark matter-dark energy scenarios are modeled via a dimensionless parameter $\xi$, which controls the strength of their interaction. While this coupling is commonly assumed to be constant, there is no underlying physical law or symmetry that forbids a time-dependent $\xi$ parameter. The most general and complete interacting scenarios between the two dark sectors should therefore allow for such a possibility, and it is the main purpose of this study to constrain two possible and well-motivated coupled cosmologies by means of the most recent and accurate early and late-time universe observations. We find that CMB data alone prefers $\xi(z) >0$ and therefore a smaller amount of dark matter, alleviating some crucial and well-known cosmological data tensions. An objective assessment of the Bayesian evidence for the coupled models explored here shows no particular preference for the presence of a dynamical dark sector coupling.

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 consider a quintessence field which transitions from a matter-like to a cosmological constant behavior between recombination and the present time. We aim at easing the tension in the measurement of the present Hubble rate, and we assess the $\Lambda$CDM model properly enlarged to include our quintessence field against cosmological observations. The model does not address the scope we proposed. This result allows us to exclude a class of quintessential models as a solution to the tension in the Hubble constant measurements.

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.

The Probe of Inflation and Cosmic Origins (PICO) is a proposed probe-scale space mission consisting of an imaging polarimeter operating in frequency bands between 20 and 800 GHz. We describe the science achievable by PICO, which has sensitivity equivalent to more than 3300 Planck missions, the technical implementation, the schedule and cost.

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.

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

This paper describes the 2018 Planck CMB likelihoods, following a hybrid approach similar to the 2015 one, with different approximations at low and high multipoles, and implementing several methodological and analysis refinements. With more realistic simulations, and better correction and modelling of systematics, we can now make full use of the High Frequency Instrument polarization data. The low-multipole 100x143 GHz EE cross-spectrum constrains the reionization optical-depth parameter $\tau$ to better than 15% (in combination with with the other low- and high-$\ell$ likelihoods). We also update the 2015 baseline low-$\ell$ joint TEB likelihood based on the Low Frequency Instrument data, which provides a weaker $\tau$ constraint. At high multipoles, a better model of the temperature-to-polarization leakage and corrections for the effective calibrations of the polarization channels (polarization efficiency or PE) allow us to fully use the polarization spectra, improving the constraints on the $\Lambda$CDM parameters by 20 to 30% compared to TT-only constraints. Tests on the modelling of the polarization demonstrate good consistency, with some residual modelling uncertainties, the accuracy of the PE modelling being the main limitation. Using our various tests, simulations, and comparison between different high-$\ell$ implementations, we estimate the consistency of the results to be better than the 0.5$\sigma$ level. Minor curiosities already present before (differences between $\ell$<800 and $\ell$>800 parameters or the preference for more smoothing of the $C_\ell$ peaks) are shown to be driven by the TT power spectrum and are not significantly modified by the inclusion of polarization. Overall, the legacy Planck CMB likelihoods provide a robust tool for constraining the cosmological model and represent a reference for future CMB observations. (Abridged)

We present measurements of the cosmic microwave background (CMB) lensing potential using the final $\textit{Planck}$ 2018 temperature and polarization data. We increase the significance of the detection of lensing in the polarization maps from $5\,\sigma$ to $9\,\sigma$. Combined with temperature, lensing is detected at $40\,\sigma$. We present an extensive set of tests of the robustness of the lensing-potential power spectrum, and construct a minimum-variance estimator likelihood over lensing multipoles $8 \le L \le 400$. We find good consistency between lensing constraints and the results from the $\textit{Planck}$ CMB power spectra within the $\rm{\Lambda CDM}$ model. Combined with baryon density and other weak priors, the lensing analysis alone constrains $\sigma_8 \Omega_{\rm m}^{0.25}=0.589\pm 0.020$ ($1\,\sigma$ errors). Also combining with baryon acoustic oscillation (BAO) data, we find tight individual parameter constraints, $\sigma_8=0.811\pm0.019$, $H_0=67.9_{-1.3}^{+1.2}\,\text{km}\,\text{s}^{-1}\,\rm{Mpc}^{-1}$, and $\Omega_{\rm m}=0.303^{+0.016}_{-0.018}$. Combining with $\textit{Planck}$ CMB power spectrum data, we measure $\sigma_8$ to better than $1\,\%$ precision, finding $\sigma_8=0.811\pm 0.006$. We find consistency with the lensing results from the Dark Energy Survey, and give combined lensing-only parameter constraints that are tighter than joint results using galaxy clustering. Using $\textit{Planck}$ cosmic infrared background (CIB) maps we make a combined estimate of the lensing potential over $60\,\%$ of the sky with considerably more small-scale signal. We demonstrate delensing of the $\textit{Planck}$ power spectra, detecting a maximum removal of $40\,\%$ of the lensing-induced power in all spectra. The improvement in the sharpening of the acoustic peaks by including both CIB and the quadratic lensing reconstruction is detected at high significance (abridged).

We consider two stable Interacting Dark Matter -- Dark Energy models and confront them against current Cosmic Microwave Background data from the \textit{Planck} satellite. We then generate luminosity distance measurements from ${\cal O}(10^3)$ mock Gravitational Wave events matching the expected sensitivity of the proposed Einstein Telescope. We use these to forecast how the addition of Gravitational Wave standard sirens data can improve current limits on the Dark Matter -- Dark Energy coupling strength ($\xi$). We find that the addition of Gravitational Waves data can reduce the current uncertainty by a factor of $5$. Moreover, if the underlying cosmological model truly features Dark Matter -- Dark Energy interactions with a value of $\xi$ within the currently allowed $1\sigma$ upper limit, the addition of Gravitational Wave data would help disentangle such an interaction from the standard case of no interaction at a significance of more than $3\sigma$.

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

Analysis of the Planck 2018 data set indicates that the statistical properties of the cosmic microwave background (CMB) temperature anisotropies are in excellent agreement with previous studies using the 2013 and 2015 data releases. In particular, they are consistent with the Gaussian predictions of the $\Lambda$CDM cosmological model, yet also confirm the presence of several so-called "anomalies" on large angular scales. The novelty of the current study, however, lies in being a first attempt at a comprehensive analysis of the statistics of the polarization signal over all angular scales, using either maps of the Stokes parameters, $Q$ and $U$, or the $E$-mode signal derived from these using a new methodology (which we describe in an appendix). Although remarkable progress has been made in reducing the systematic effects that contaminated the 2015 polarization maps on large angular scales, it is still the case that residual systematics (and our ability to simulate them) can limit some tests of non-Gaussianity and isotropy. However, a detailed set of null tests applied to the maps indicates that these issues do not dominate the analysis on intermediate and large angular scales (i.e., $\ell \lesssim 400$). In this regime, no unambiguous detections of cosmological non-Gaussianity, or of anomalies corresponding to those seen in temperature, are claimed. Notably, the stacking of CMB polarization signals centred on the positions of temperature hot and cold spots exhibits excellent agreement with the $\Lambda$CDM cosmological model, and also gives a clear indication of how Planck provides state-of-the-art measurements of CMB temperature and polarization on degree scales.

We analyse the Planck full-mission cosmic microwave background (CMB) temperature and E-mode polarization maps to obtain constraints on primordial non-Gaussianity (NG). We compare estimates obtained from separable template-fitting, binned, and modal bispectrum estimators, finding consistent values for the local, equilateral, and orthogonal bispectrum amplitudes. Our combined temperature and polarization analysis produces the following results: f_NL^local = -0.9 +\- 5.1; f_NL^equil = -26 +\- 47; and f_NL^ortho = - 38 +\- 24 (68%CL, statistical). These results include the low-multipole (4 <= l < 40) polarization data, not included in our previous analysis, pass an extensive battery of tests, and are stable with respect to our 2015 measurements. Polarization bispectra display a significant improvement in robustness; they can now be used independently to set NG constraints. We consider a large number of additional cases, e.g. scale-dependent feature and resonance bispectra, isocurvature primordial NG, and parity-breaking models, where we also place tight constraints but do not detect any signal. The non-primordial lensing bispectrum is detected with an improved significance compared to 2015, excluding the null hypothesis at 3.5 sigma. We present model-independent reconstructions and analyses of the CMB bispectrum. Our final constraint on the local trispectrum shape is g_NLl^local = (-5.8 +\-6.5) x 10^4 (68%CL, statistical), while constraints for other trispectra are also determined. We constrain the parameter space of different early-Universe scenarios, including general single-field models of inflation, multi-field and axion field parity-breaking models. Our results provide a high-precision test for structure-formation scenarios, in complete agreement with the basic picture of the LambdaCDM cosmology regarding the statistics of the initial conditions (abridged).

Following the pioneering observations with COBE in the early 1990s, studies of the cosmic microwave background (CMB) have focused on temperature and polarization anisotropies. CMB spectral distortions - tiny departures of the CMB energy spectrum from that of a perfect blackbody - provide a second, independent probe of fundamental physics, with a reach deep into the primordial Universe. The theoretical foundation of spectral distortions has seen major advances in recent years, which highlight the immense potential of this emerging field. Spectral distortions probe a fundamental property of the Universe - its thermal history - thereby providing additional insight into processes within the cosmological standard model (CSM) as well as new physics beyond. Spectral distortions are an important tool for understanding inflation and the nature of dark matter. They shed new light on the physics of recombination and reionization, both prominent stages in the evolution of our Universe, and furnish critical information on baryonic feedback processes, in addition to probing primordial correlation functions at scales inaccessible to other tracers. In principle the range of signals is vast: many orders of magnitude of discovery space could be explored by detailed observations of the CMB energy spectrum. Several CSM signals are predicted and provide clear experimental targets, some of which are already observable with present-day technology. Confirmation of these signals would extend the reach of the CSM by orders of magnitude in physical scale as the Universe evolves from the initial stages to its present form. The absence of these signals would pose a huge theoretical challenge, immediately pointing to new physics.

Here we test the predictions of the theory of the origin of the universe from the landscape multiverse, against the 2015 Planck data, for the case of the Hilltop class of inflationary models, for $p=4$ and $p=6$. By considering the quantum entanglement correction of the multiverse, we can place just a lower limit on the local 'SUSY-breaking' scale, respectively $b>8.7\times10^6 GeV$ at $95 \%$ c.l. and $b>1.3\times10^8 GeV$ at $95 \%$ c.l. from Planck TT+lowP, so the case with multiverse correction is statistically indistinguishable from the case with an unmodified inflation. We find that the series of anomalies predicted by the quantum landscape multiverse for the allowed range of $b$, is consistent with Planck's tests of the anomalies. In addition, the friction between the two cosmological probes of the Hubble parameter and with the weak lensing experiments goes away for a particular subset, the $p=6$ case of Hilltop models.

The expansion of the Universe is understood to have accelerated during two epochs: in its very first moments during a period of Inflation and much more recently, at $z < 1$, when Dark Energy is hypothesized to drive cosmic acceleration. The undiscovered mechanisms behind these two epochs represent some of the most important open problems in fundamental physics. The large cosmological volume at $2 < z < 5$, together with the ability to efficiently target high-$z$ galaxies with known techniques, enables large gains in the study of Inflation and Dark Energy. A future spectroscopic survey can test the Gaussianity of the initial conditions up to a factor of ~50 better than our current bounds, crossing the crucial theoretical threshold of $\sigma(f_{NL}^{\rm local})$ of order unity that separates single field and multi-field models. Simultaneously, it can measure the fraction of Dark Energy at the percent level up to $z = 5$, thus serving as an unprecedented test of the standard model and opening up a tremendous discovery space.

During the next decade, gravitational waves will be observed from hundreds of binary inspiral events. When the redshifts of the host galaxies are known, these events can be used as `standard sirens', sensitive to the expansion rate of the Universe. Measurements of the Hubble constant $H_0$ from standard sirens can be done independently of other cosmological probes, and events occurring at $z<0.1$ will allow us to infer $H_0$ independently of cosmological models. The next generation of spectroscopic galaxy surveys will play a crucial role in reducing systematic uncertainties in $H_0$ from standard sirens, particularly for the numerous `dark sirens' which do not have an electromagnetic counterpart. In combination with large spectroscopic data sets, standard sirens with an EM counterpart are expected to constrain $H_0$ to $\sim 1-2\%$ precision within the next decade. This is competitive with the best estimates of $H_0$ obtained to date and will help illuminate the current tension between existing measurements. Information on the galaxies that host the gravitational wave events will also shed light on the origin and evolution of compact object binaries.

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

Astrophysical observations currently provide the only robust, empirical measurements of dark matter. In the coming decade, astrophysical observations will guide other experimental efforts, while simultaneously probing unique regions of dark matter parameter space. This white paper summarizes astrophysical observations that can constrain the fundamental physics of dark matter in the era of LSST. We describe how astrophysical observations will inform our understanding of the fundamental properties of dark matter, such as particle mass, self-interaction strength, non-gravitational interactions with the Standard Model, and compact object abundances. Additionally, we highlight theoretical work and experimental/observational facilities that will complement LSST to strengthen our understanding of the fundamental characteristics of dark matter.

We present 353 GHz full-sky maps of the polarization fraction $p$, angle $\psi$, and dispersion of angles $S$ of Galactic dust thermal emission produced from the 2018 release of Planck data. We confirm that the mean and maximum of $p$ decrease with increasing $N_H$. The uncertainty on the maximum polarization fraction, $p_\mathrm{max}=22.0$% at 80 arcmin resolution, is dominated by the uncertainty on the zero level in total intensity. The observed inverse behaviour between $p$ and $S$ is interpreted with models of the polarized sky that include effects from only the topology of the turbulent Galactic magnetic field. Thus, the statistical properties of $p$, $\psi$, and $S$ mostly reflect the structure of the magnetic field. Nevertheless, we search for potential signatures of varying grain alignment and dust properties. First, we analyse the product map $S \times p$, looking for residual trends. While $p$ decreases by a factor of 3--4 between $N_H=10^{20}$ cm$^{-2}$ and $N_H=2\times 10^{22}$ cm$^{-2}$, $S \times p$ decreases by only about 25%, a systematic trend observed in both the diffuse ISM and molecular clouds. Second, we find no systematic trend of $S \times p$ with the dust temperature, even though in the diffuse ISM lines of sight with high $p$ and low $S$ tend to have colder dust. We also compare Planck data with starlight polarization in the visible at high latitudes. The agreement in polarization angles is remarkable. Two polarization emission-to-extinction ratios that characterize dust optical properties depend only weakly on $N_H$ and converge towards the values previously determined for translucent lines of sight. We determine an upper limit for the polarization fraction in extinction of 13%, compatible with the $p_\mathrm{max}$ observed in emission. These results provide strong constraints for models of Galactic dust in diffuse gas.

The Probe of Inflation and Cosmic Origins (PICO) is an imaging polarimeter that will scan the sky for 5 years in 21 frequency bands spread between 21 and 799 GHz. It will produce full-sky surveys of intensity and polarization with a final combined-map noise level of 0.87 $\mu$K arcmin for the required specifications, equivalent to 3300 Planck missions, and with our current best-estimate would have a noise level of 0.61 $\mu$K arcmin (6400 Planck missions). PICO will either determine the energy scale of inflation by detecting the tensor to scalar ratio at a level $r=5\times 10^{-4}~(5\sigma)$, or will rule out with more than $5\sigma$ all inflation models for which the characteristic scale in the potential is the Planck scale. With LSST's data it could rule out all models of slow-roll inflation. PICO will detect the sum of neutrino masses at $>4\sigma$, constrain the effective number of light particle species with $\Delta N_{\rm eff}<0.06~(2\sigma)$, and elucidate processes affecting the evolution of cosmic structures by measuring the optical depth to reionization with errors limited by cosmic variance and by constraining the evolution of the amplitude of linear fluctuations $\sigma_{8}(z)$ with sub-percent accuracy. Cross-correlating PICO's map of the thermal Sunyaev-Zeldovich effect with LSST's gold sample of galaxies will precisely trace the evolution of thermal pressure with $z$. PICO's maps of the Milky Way will be used to determine the make up of galactic dust and the role of magnetic fields in star formation efficiency. With 21 full sky legacy maps in intensity and polarization, which cannot be obtained in any other way, the mission will enrich many areas of astrophysics. PICO is the only single-platform instrument with the combination of sensitivity, angular resolution, frequency bands, and control of systematic effects that can deliver this compelling, timely, and broad science.

The phenomenological parametrizations of dark-energy (DE) equation of state can be very helpful, since they allow for the investigation of its cosmological behavior despite the fact that its underlying theory is unknown. However, although there has been a large amount of research on DE parametrizations which involve two or more free parameters, the one-parameter parametrizations seem to be underestimated. We perform a detailed observational confrontation of five one-parameter DE models, with observational data from cosmic microwave background (CMB), Joint light-curve analysis sample from Supernovae Type Ia observations (JLA), baryon acoustic oscillations (BAO) distance measurements, and cosmic chronometers (CC). We find that all models favor a phantom DE equation of state at present time, while they lead to $H_0$ values in perfect agreement with its direct measurements and therefore they offer an alleviation to the $H_0$-tension. Finally, performing a Bayesian analysis we show that although $\Lambda$CDM cosmology is still favored, one-parameter DE models have similar or better efficiency in fitting the data comparing to two-parameter DE parametrizations, and thus they deserve a thorough investigation.

Almost in all interacting dark energy models present in the literature, the stability of the model becomes potentially sensitive to the dark energy equation of state parameter $w_x$, and a singularity arises at `$w_x = -1$'. Thus, it becomes mandatory to test the stability of the model into two separate regions, namely, for quintessence and phantom. This essentially brings in a discontinuity into the parameters space for $w_x$. Such discontinuity can be removed with some specific choices of the interaction or coupling function. In the present work we choose one particular coupling between dark matter and dark energy which can successfully remove such instability and we allow a dynamical dark energy equation of state parameter instead of the constant one. In particular, considering a dynamical dark energy equation of state with only one free parameter $w_0$, representing the current value of the dark energy equation of state, we confront the interacting scenario with several observational datasets. The results show that the present cosmological data allow an interaction in the dark sector, in agreement with some latest claims by several authors, and additionally, a phantom behaviour in the dark energy equation of state is suggested at present. Moreover, for this case the tension on $H_0$ is clearly released. As a final remark, we mention that according to the Bayesian analysis, $\Lambda$-cold dark matter ($\Lambda$CDM) is always favored over this interacting dark energy model.

We investigate the generalized Chevallier-Polarski-Linder (CPL) parametrization, which contains the pivoting redshift $z_p$ as an extra free parameter. We use various data combinations from cosmic microwave background (CMB), baryon acoustic oscillations (BAO), redshift space distortion (RSD), weak lensing (WL), joint light curve analysis (JLA), cosmic chronometers (CC), and we include a Gaussian prior on the Hubble constant value, in order to extract the observational constraints on various quantities. For the case of free $z_p$ we find that for all data combinations it always remains unconstrained, and there is a degeneracy with the current value of the dark energy equation of state $w_0$. For the case where $z_p$ is fixed to specific values, and for the full data combination, we find that with increasing $z_p$ the mean value of $w_0$ slowly moves into the phantom regime, however the cosmological constant is always allowed within 1$\sigma$ confidence-level. However, the significant effect is that with increasing $z_p$ the correlations between $w_0$ and $w_a$ change from negative to positive, with the case $z_p =0.35$ corresponding to no correlation. This feature indeed justifies why a non-zero pivoting redshift should be taken into account.

The study of polarized dust emission has become entwined with the analysis of the cosmic microwave background (CMB) polarization. We use new Planck maps to characterize Galactic dust emission as a foreground to the CMB polarization. We present Planck EE, BB, and TE power spectra of dust polarization at 353 GHz for six nested sky regions covering from 24 to 71 % of the sky. We present power-law fits to the angular power spectra, yielding evidence for statistically significant variations of the exponents over sky regions and a difference between the values for the EE and BB spectra. The TE correlation and E/B power asymmetry extend to low multipoles that were not included in earlier Planck polarization papers. We also report evidence for a positive TB dust signal. Combining data from Planck and WMAP, we determine the amplitudes and spectral energy distributions (SEDs) of polarized foregrounds, including the correlation between dust and synchrotron polarized emission, for the six sky regions as a function of multipole. This quantifies the challenge of the component separation procedure required for detecting the reionization and recombination peaks of primordial CMB B modes. The SED of polarized dust emission is fit well by a single-temperature modified blackbody emission law from 353 GHz to below 70 GHz. For a dust temperature of 19.6 K, the mean spectral index for dust polarization is $\beta_{\rm d}^{P} = 1.53\pm0.02 $. By fitting multi-frequency cross-spectra, we examine the correlation of the dust polarization maps across frequency. We find no evidence for decorrelation. If the Planck limit for the largest sky region applies to the smaller sky regions observed by sub-orbital experiments, then decorrelation might not be a problem for CMB experiments aiming at a primordial B-mode detection limit on the tensor-to-scalar ratio $r\simeq0.01$ at the recombination peak.

Gravitational-wave standard sirens present a novel approach for the determination of the Hubble constant. After the recent spectacular confirmation of the method thanks to GW170817 and its optical counterpart, additional standard siren measurements from future gravitational-wave sources are expected to constrain the Hubble constant to high accuracy. At the same time, improved constraints are expected from observations of cosmic microwave background (CMB) polarization and from baryon acoustic oscillations (BAO) surveys. We explore the role of future standard siren constraints on $H_0$ in light of expected CMB+BAO data. Considering a $10$-parameters cosmological model, in which curvature, the dark energy equation of state, and the Hubble constant are unbounded by CMB observations, we find that a combination of future CMB+BAO data will constrain the Hubble parameter to $\sim 1.5 \%$. Further extending the parameter space to a time-varying dark energy equation of state, we find that future CMB+BAO constraints on $H_0$ are relaxed to $\sim 3.0 \%$. These accuracies are within reach of future standard siren measurements from the Hanford-Livingston-Virgo and the Hanford-Livingston-Virgo-Japan-India networks of interferometers, showing the cosmological relevance of these sources. If future gravitational-wave standard siren measurements reach $1\%$ on $H_0$, as expected, they would significantly improve future CMB+BAO constraints on curvature and on the dark energy equation of state by up to a factor $\sim 3$. We also show that the inclusion of $H_0$ constraints from gravitational-wave standard sirens could result in a reduction of the dark energy figure-of-merit (i.e., the cosmological parameter volume) by up to a factor of $\sim 400$.

We present a final description of the data-processing pipeline for the Planck, Low Frequency Instrument (LFI), implemented for the 2018 data release. Several improvements have been made with respect to the previous release, especially in the calibration process and in the correction of instrumental features such as the effects of nonlinearity in the response of the analogue-to-digital converters. We provide a brief pedagogical introduction to the complete pipeline, as well as a detailed description of the important changes implemented. Self-consistency of the pipeline is demonstrated using dedicated simulations and null tests. We present the final version of the LFI full sky maps at 30, 44, and 70 GHz, both in temperature and polarization, together with a refined estimate of the Solar dipole and a final assessment of the main LFI instrumental parameters.

This paper presents the Planck Multi-frequency Catalogue of Non-thermal (i.e. synchrotron-dominated) Sources (PCNT) observed between 30 and 857 GHz by the ESA Planck mission. This catalogue was constructed by selecting objects detected in the full mission all-sky temperature maps at 30 and 143 GHz, with a signal-to-noise ratio (S/N)>3 in at least one of the two channels after filtering with a particular Mexican hat wavelet. As a result, 29400 source candidates were selected. Then, a multi-frequency analysis was performed using the Matrix Filters methodology at the position of these objects, and flux densities and errors were calculated for all of them in the nine Planck channels. The present catalogue is the first unbiased, full-sky catalogue of synchrotron-dominated sources published at millimetre and submillimetre wavelengths and constitutes a powerful database for statistical studies of non-thermal extragalactic sources, whose emission is dominated by the central active galactic nucleus. Together with the full multi-frequency catalogue, we also define the Bright Planck Multi-frequency Catalogue of Non-thermal Sources PCNTb, where only those objects with a S/N>4 at both 30 and 143 GHz were selected. In this catalogue 1146 compact sources are detected outside the adopted Planck GAL070 mask; thus, these sources constitute a highly reliable sample of extragalactic radio sources. We also flag the high-significance subsample PCNThs, a subset of 151 sources that are detected with S/N>4 in all nine Planck channels, 75 of which are found outside the Planck mask adopted here. The remaining 76 sources inside the Galactic mask are very likely Galactic objects.

We investigate the observational consequences of a novel class of stable interacting dark energy (IDE) models, featuring interactions between dark matter (DM) and dark energy (DE). In the first part of our work, we start by considering two IDE models which are known to present early-time linear perturbation instabilities. Applying a transformation depending on the dark energy equation of state (EoS) to the DM-DE coupling, we then obtain two novel stable IDE models. Subsequently, we derive robust and accurate constraints on the parameters of these models, assuming a constant EoS $w_x$ for the DE fluid, in light of some of the most recent publicly available cosmological data. These include Cosmic Microwave Background (CMB) temperature and polarization anisotropy measurements from the \textit{Planck} satellite, a selection of Baryon Acoustic Oscillation measurements, Supernovae Type-Ia luminosity distance measurements from the JLA sample, and measurements of the Hubble parameter up to redshift $2$ from cosmic chronometers. Our analysis displays a mild preference for the DE fluid residing in the phantom region ($w_x<-1$), with significance up to 95\% confidence level, while we obtain new upper limits on the coupling parameter between the dark components. The preference for a phantom DE suggests a coupling function $Q<0$, thus a scenario where energy flows from the DE to the DM. We also examine the possibility of addressing the $H_0$ and $\sigma_8$ tensions, finding that only the former can be partially alleviated. Finally, we perform a Bayesian model comparison analysis to quantify the possible preference for the two IDE models against the standard concordance $\Lambda$CDM model, finding that the latter is always preferred with the strength of the evidence ranging from positive to very strong.

Current cosmological constraints on the scalar spectral index of primordial fluctuations $n_{\rm s}$ in the $\Lambda$CDM model have excluded the minimal scale-invariant Harrison-Zel'dovich model ($n_{\rm s}=1$; hereafter HZ) at high significance, providing support for inflation. In recent years, however, some tensions have emerged between different cosmological datasets that, if not due to systematics, could indicate the presence of new physics beyond the $\Lambda$CDM model. In the light of these developments, we evaluate the Bayesian evidence against HZ in different data combinations and model extensions. Considering only the Planck temperature data, we find inconclusive evidence against HZ when including variations in the neutrino number $N_{\rm eff}$ and/or the Helium abundance $Y_{\rm He}$. Adding the Planck polarization data, on the other hand, yields strong evidence against HZ in the extensions we considered. Perhaps most interestingly, Planck temperature data combined with local measurements of the Hubble constant give as the most probable model an HZ spectrum, with additional neutrinos. However, with the inclusion of polarisation, standard $\Lambda$CDM is once again preferred, but the HZ model with extra neutrinos is not strongly disfavored. The possibility of fully ruling out the HZ spectrum is therefore ultimately connected with the solution to current tensions between cosmological datasets. If these tensions are confirmed by future data, then new physical mechanisms could be at work and an HZ spectrum could still offer a valid alternative.

We discuss whether an unaccounted contribution to the Cosmic Microwave Background polarization $B$-mode by primordial magnetic fields (PMFs) can bias future constraints on inflationary gravitational waves. As a case-study, we consider a scale-invariant PMF spectrum with amplitude of $\approx{1}\,{\mathrm{nG}}$ on $1\,\mathrm{Mpc}$ scales, compatible with current cosmological bounds. We find a degeneracy in the $B$-mode spectra between PMFs and inflationary gravitational waves. If PMFs of this amplitude are not accounted for, future CMB experiments could claim a false detection of a tensor-to-scalar ratio $r\approx 0.007$, close to the predictions of Starobinsky and $\alpha$-attractor models. The degeneracy can be broken if $B$-modes are measured also at multipoles $\ell\gtrsim 900$: more precisely experiments like CMB-S4 or CORE-M5 would be able to discriminate PMFs from primordial GWs at high statistical significance. Experiments like LiteBIRD or PIXIE will not be able to break the degeneracy and will need complementary bounds coming, for example, from measurements of anisotropies in the Faraday rotation angle of CMB polarization. This reinforces the importance of future experimental constraints on PMFs.

Using the ${\it Planck}$ full-mission data, we present a detection of the temperature (and therefore velocity) dispersion due to the kinetic Sunyaev-Zeldovich (kSZ) effect from clusters of galaxies. To suppress the primary CMB and instrumental noise we derive a matched filter and then convolve it with the ${\it Planck}$ foreground-cleaned `${\tt 2D-ILC\,}$' maps. By using the Meta Catalogue of X-ray detected Clusters of galaxies (MCXC), we determine the normalized ${\it rms}$ dispersion of the temperature fluctuations at the positions of clusters, finding that this shows excess variance compared with the noise expectation. We then build an unbiased statistical estimator of the signal, determining that the normalized mean temperature dispersion of $1526$ clusters is $\langle \left(\Delta T/T \right)^{2} \rangle = (1.64 \pm 0.48) \times 10^{-11}$. However, comparison with analytic calculations and simulations suggest that around $0.7\,\sigma$ of this result is due to cluster lensing rather than the kSZ effect. By correcting this, the temperature dispersion is measured to be $\langle \left(\Delta T/T \right)^{2} \rangle = (1.35 \pm 0.48) \times 10^{-11}$, which gives a detection at the $2.8\,\sigma$ level. We further convert uniform-weight temperature dispersion into a measurement of the line-of-sight velocity dispersion, by using estimates of the optical depth of each cluster (which introduces additional uncertainty into the estimate). We find that the velocity dispersion is $\langle v^{2} \rangle =(123\,000 \pm 71\,000)\,({\rm km}\,{\rm s}^{-1})^{2}$, which is consistent with findings from other large-scale structure studies, and provides direct evidence of statistical homogeneity on scales of $600\,h^{-1}{\rm Mpc}$. Our study shows the promise of using cross-correlations of the kSZ effect with large-scale structure in order to constrain the growth of structure.

We present full-sky maps of the cosmic microwave background (CMB) and polarized synchrotron and thermal dust emission, derived from the third set of Planck frequency maps. These products have significantly lower contamination from instrumental systematic effects than previous versions. The methodologies used to derive these maps follow closely those described in earlier papers, adopting four methods (Commander, NILC, SEVEM, and SMICA) to extract the CMB component, as well as three methods (Commander, GNILC, and SMICA) to extract astrophysical components. Our revised CMB temperature maps agree with corresponding products in the Planck 2015 delivery, whereas the polarization maps exhibit significantly lower large-scale power, reflecting the improved data processing described in companion papers; however, the noise properties of the resulting data products are complicated, and the best available end-to-end simulations exhibit relative biases with respect to the data at the few percent level. Using these maps, we are for the first time able to fit the spectral index of thermal dust independently over 3 degree regions. We derive a conservative estimate of the mean spectral index of polarized thermal dust emission of beta_d = 1.55 +/- 0.05, where the uncertainty marginalizes both over all known systematic uncertainties and different estimation techniques. For polarized synchrotron emission, we find a mean spectral index of beta_s = -3.1 +/- 0.1, consistent with previously reported measurements. We note that the current data processing does not allow for construction of unbiased single-bolometer maps, and this limits our ability to extract CO emission and correlated components. The foreground results for intensity derived in this paper therefore do not supersede corresponding Planck 2015 products. For polarization the new results supersede the corresponding 2015 products in all respects.

This paper presents the High Frequency Instrument (HFI) data processing procedures for the Planck 2018 release. Major improvements in mapmaking have been achieved since the previous 2015 release. They enabled the first significant measurement of the reionization optical depth parameter using HFI data. This paper presents an extensive analysis of systematic effects, including the use of simulations to facilitate their removal and characterize the residuals. The polarized data, which presented a number of known problems in the 2015 Planck release, are very significantly improved. Calibration, based on the CMB dipole, is now extremely accurate and in the frequency range 100 to 353 GHz reduces intensity-to-polarization leakage caused by calibration mismatch. The Solar dipole direction has been determined in the three lowest HFI frequency channels to within one arc minute, and its amplitude has an absolute uncertainty smaller than $0.35\mu$K, an accuracy of order $10^{-4}$. This is a major legacy from the HFI for future CMB experiments. The removal of bandpass leakage has been improved by extracting the bandpass-mismatch coefficients for each detector as part of the mapmaking process; these values in turn improve the intensity maps. This is a major change in the philosophy of "frequency maps", which are now computed from single detector data, all adjusted to the same average bandpass response for the main foregrounds. Simulations reproduce very well the relative gain calibration of detectors, as well as drifts within a frequency induced by the residuals of the main systematic effect. Using these simulations, we measure and correct the small frequency calibration bias induced by this systematic effect at the $10^{-4}$ level. There is no detectable sign of a residual calibration bias between the first and second acoustic peaks in the CMB channels, at the $10^{-3}$ level.

The precise measurements of Cosmic Microwave Background Anisotropy angular power spectra made by the Planck satellite show an anomalous value for the lensing amplitude, defined by the parameter $A_{lens}$, at more than $2$ standard deviations. In this paper, after discussing the current status of the anomaly, we quantify the potential of future CMB measurements in confirming/falsifying the $A_{lens}$ tension. We find that a space-based experiment as LiteBIRD could falsify the current $A_{lens}$ tension at the level of $5$ standard deviations. Similar constraints can be achieved by a Stage-III experiment assuming an external prior on the reionization optical depth of $\tau=0.055\pm0.010$ as already provided by the Planck satellite. A Stage-IV experiment could further test the $A_{lens}$ tension at the level of $10$ standard deviations. A comparison between temperature and polarization measurements made at different frequencies could further identify possible systematics responsible for $A_{lens}>1$. We show that, in the case of the CMB-S4 experiment, polarization data alone will have the potential of falsifying the current $A_{lens}$ anomaly at more than five standard deviation and to strongly bound its frequency dependence. We also evaluate the future constraints on a possible scale dependence for $A_{lens}$.

One of the major goals of future Cosmic Microwave Background measurements is the accurate determination of the effective number of neutrinos $N_{\rm eff}$. Reaching an experimental sensitivity of $\Delta N_{\rm eff} = 0.013$ could indeed falsify the presence of any non-standard relativistic particles at $95 \%$ c.l.. In this paper, we test how this future constraint can be affected by the removal of two common assumptions: a negligible running of the inflationary spectral index $n_{\rm run}$ and a precise determination of the neutron lifetime $\tau_n$. We first show that the constraints on $N_{\rm eff}$ could be significantly biased by the unaccounted presence of a running of the spectral index. Considering the Stage-IV experiment, a negative running of ${\rm d}n/{\rm d}\ln k= - 0.002$ could mimic a positive variation of $\Delta N_{\rm eff} = 0.03$. Moreover, given the current discrepancies between experimental measurements of the neutron lifetime $\tau_n$, we show that the assumption of a conservative error of $\Delta\tau_n \sim 10$s could bring to a systematic error of $\Delta N_{\rm eff} = 0.02$. Complementary cosmological constraints on the running of the spectral index and a solution to the neutron lifetime discrepancy are therefore needed for an accurate and reliable future CMB bound of $N_{\rm eff}$ at percent level.

The recent observations of gravitational-wave and electromagnetic emission produced by the merger of the binary neutron-star system GW170817 have opened the possibility of using standard sirens to constrain the value of the Hubble constant. While the reported bound of $H_0=70_{-8}^{+12}$ at $68 \%$ C.L. is significantly weaker than those recently derived by observations of Cepheid variables, it does not require any form of cosmic distance ladder and can be considered as complementary and, in principle, more conservative. Here we combine, for the first time, the new measurement with the Planck Cosmic Microwave Background observations in a $12$ parameters extended $\Lambda$CDM scenario, where the Hubble constant is weakly constrained from CMB data alone and bound to a low value $H_0=55^{+7}_{-20}$ km/s/Mpc at $68 \%$ C.L. We point out that the non-Gaussian shape of the GW170817 bound makes lower values of the Hubble constant in worst agreement with observations than what expected from a Gaussian form. The inclusion of the new GW170817 Hubble constant measurement therefore significantly reduces the allowed parameter space, improving the cosmological bounds on several parameters as the neutrino mass, curvature and the dark energy equation of state.

The introduction of Dark Matter-neutrino interactions modifies the Cosmic Microwave Background (CMB) angular power spectrum at all scales, thus affecting the reconstruction of the cosmological parameters. Such interactions can lead to a slight increase of the value of $H_0$ and a slight decrease of $\sigma_8$, which can help reduce somewhat the tension between the CMB and lensing or Cepheids datasets. Here we show that it is impossible to solve both tensions simultaneously. While the 2015 Planck temperature and low multipole polarisation data combined with the Cepheids datasets prefer large values of the Hubble rate (up to $H_0 = 72.1^{+1.5}_{-1.7} \rm{km/s/Mpc}$, when $N_{\rm{eff}}$ is free to vary), the $\sigma_8$ parameter remains too large to reduce the $\sigma_8$ tension. Adding high multipole Planck polarization data does not help since this data shows a strong preference for low values of $H_0$, thus worsening current tensions, even though they also prefer smaller value of $\sigma_8$.

Taking the Planck cosmic microwave background data and the more direct Hubble constant measurement data as unaffected by systematic offsets, the values of the Hubble constant $H_0$ interpreted within the $\Lambda$CDM cosmological constant and cold dark matter cosmological model are in $\sim 3.3 \sigma$ tension. We show that the Parker vacuum metamorphosis model, physically motivated by quantum gravitational effects and with the same number of parameters as $\Lambda$CDM, can remove the $H_0$ tension, and can give an improved fit to data (up to $\Delta\chi^2=-7.5$). It also ameliorates tensions with weak lensing data and the high redshift Lyman alpha forest data. We separately consider a scale dependent scaling of the gravitational lensing amplitude, such as provided by modified gravity, neutrino mass, or cold dark energy, motivated by the somewhat different cosmological parameter estimates for low and high CMB multipoles. We find that no such scale dependence is preferred.

A time-dependent dark energy component of the Universe may be able to explain tensions between local and primordial measurements of cosmological parameters, shaking current confidence in the concept of a cosmological constant.

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]

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.

The answer is Yes! We indeed find that interacting dark energy can alleviate the current tension on the value of the Hubble constant $H_0$ between the Cosmic Microwave Background anisotropies constraints obtained from the Planck satellite and the recent direct measurements reported by Riess et al. 2016. The combination of these two datasets points towards an evidence for a non-zero dark matter-dark energy coupling $\xi$ at more than two standard deviations, with $\xi=-0.26_{-0.12}^{+0.16}$ at $95\%$ CL. However the $H_0$ tension is better solved when the equation of state of the interacting dark energy component is allowed to freely vary, with a phantom-like equation of state $w=-1.184\pm0.064$ (at $68 \%$ CL), ruling out the pure cosmological constant case, $w=-1$, again at more than two standard deviations. When Planck data are combined with external datasets, as BAO, JLA Supernovae Ia luminosity distances, cosmic shear or lensing data, we find good consistency with the cosmological constant scenario and no compelling evidence for a dark matter-dark energy coupling.

The six parameters of the standard $\Lambda$CDM model have best-fit values derived from the Planck temperature power spectrum that are shifted somewhat from the best-fit values derived from WMAP data. These shifts are driven by features in the Planck temperature power spectrum at angular scales that had never before been measured to cosmic-variance level precision. We investigate these shifts to determine whether they are within the range of expectation and to understand their origin in the data. Taking our parameter set to be the optical depth of the reionized intergalactic medium $\tau$, the baryon density $\omega_{\rm b}$, the matter density $\omega_{\rm m}$, the angular size of the sound horizon $\theta_*$, the spectral index of the primordial power spectrum, $n_{\rm s}$, and $A_{\rm s}e^{-2\tau}$ (where $A_{\rm s}$ is the amplitude of the primordial power spectrum), we examine the change in best-fit values between a WMAP-like large angular-scale data set (with multipole moment $\ell<800$ in the Planck temperature power spectrum) and an all angular-scale data set ($\ell<2500$ Planck temperature power spectrum), each with a prior on $\tau$ of $0.07\pm0.02$. We find that the shifts, in units of the 1$\sigma$ expected dispersion for each parameter, are $\{\Delta \tau, \Delta A_{\rm s} e^{-2\tau}, \Delta n_{\rm s}, \Delta \omega_{\rm m}, \Delta \omega_{\rm b}, \Delta \theta_*\} = \{-1.7, -2.2, 1.2, -2.0, 1.1, 0.9\}$, with a $\chi^2$ value of 8.0. We find that this $\chi^2$ value is exceeded in 15% of our simulated data sets, and that a parameter deviates by more than 2.2$\sigma$ in 9% of simulated data sets, meaning that the shifts are not unusually large. Comparing $\ell<800$ instead to $\ell>800$, or splitting at a different multipole, yields similar results. We examine the $\ell<800$ model residuals in the $\ell>800$ power spectrum data and find that the features there... [abridged]

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)

Dynamical dark energy has been recently suggested as a promising and physical way to solve the 3.4 sigma tension on the value of the Hubble constant $H_0$ between the direct measurement of Riess et al. (2016) (R16, hereafter) and the indirect constraint from Cosmic Microwave Anisotropies obtained by the Planck satellite under the assumption of a $\Lambda$CDM model. In this paper, by parameterizing dark energy evolution using the $w_0$-$w_a$ approach, and considering a $12$ parameter extended scenario, we find that: a) the tension on the Hubble constant can indeed be solved with dynamical dark energy, b) a cosmological constant is ruled out at more than $95 \%$ c.l. by the Planck+R16 dataset, and c) all of the standard quintessence and half of the "downward going" dark energy model space (characterized by an equation of state that decreases with time) is also excluded at more than $95 \%$ c.l. These results are further confirmed when cosmic shear, CMB lensing, or SN~Ia luminosity distance data are also included. However, tension remains with the BAO dataset. A cosmological constant and small portion of the freezing quintessence models are still in agreement with the Planck+R16+BAO dataset at between 68\% and 95\% c.l. Conversely, for Planck plus a phenomenological $H_0$ prior, both thawing and freezing quintessence models prefer a Hubble constant of less than 70 km/s/Mpc. The general conclusions hold also when considering models with non-zero spatial curvature.

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]

Within the standard three-neutrino framework, the absolute neutrino masses and their ordering (either normal, NO, or inverted, IO) are currently unknown. However, the combination of current data coming from oscillation experiments, neutrinoless double beta decay searches, and cosmological surveys, can provide interesting constraints for such unknowns in the sub-eV mass range, down to O(0.1) eV in some cases. We discuss current limits on absolute neutrino mass observables by performing a global data analysis, that includes the latest results from oscillation experiments, neutrinoless double beta decay bounds from the KamLAND-Zen experiment, and constraints from representative combinations of Planck measurements and other cosmological data sets. In general, NO appears to be somewhat favored with respect to IO at the level of ~2 sigma, mainly by neutrino oscillation data (especially atmospheric), corroborated by cosmological data in some cases. Detailed constraints are obtained via the chi^2 method, by expanding the parameter space either around separate minima in NO and IO, or around the absolute minimum in any ordering. Implications for upcoming oscillation and non-oscillation neutrino experiments, including beta-decay searches, are also discussed.

The 2015 Planck data release has placed tight constraints on the allowed class of inflationary models. The current data favors concave downwards inflationary potentials while offering interesting hints on possible deviations from the standard picture of CMB perturbations. We here test the predictions of the theory of the origin of the universe from the landscape multiverse, against the most recent Planck data, for the case of concave downwards inflationary potentials, such as the Starobinsky model of inflation. By considering the quantum entanglement correction of the multiverse, we can place a lower limit on the local 'SUSY breaking' scale $b>1.2\times10^7 GeV$ at $95 \%$ c.l. from Planck TT+lowTEB. We find that this limit is consistent with the range for $b$ that allows the landscape multiverse to explain a serie of anomalies present in the current data.

The 2015 Planck data release tightened the region of the allowed inflationary models. Inflationary models with convex potentials have now been ruled out since they produce a large tensor to scalar ratio. Meanwhile the same data offers interesting hints on possible deviations from the standard picture of CMB perturbations. Here we revisit the predictions of the theory of the origin of the universe from the landscape multiverse for the case of exponential inflation, for two reasons: firstly to check the status of the anomalies associated with this theory, in the light of the recent Planck data; secondly, to search for a counterexample whereby new physics modifications may bring convex inflationary potentials, thought to have been ruled out, back into the region of potentials allowed by data. Using the exponential inflation as an example of convex potentials, we find that the answer to both tests is positive: modifications to the perturbation spectrum and to the Newtonian potential of the universe originating from the quantum entanglement, bring the exponential potential, back within the allowed region of current data; and, the series of anomalies previously predicted in this theory, is still in good agreement with current data. Hence our finding for this convex potential comes at the price of allowing for additional thermal relic particles, equivalently dark radiation, in the early universe.

Measurements of flux density are described for five planets, Mars, Jupiter, Saturn, Uranus, and Neptune, across the six Planck High Frequency Instrument frequency bands (100-857 GHz) and these are then compared with models and existing data. In our analysis, we have also included estimates of the brightness of Jupiter and Saturn at the three frequencies of the Planck Low Frequency Instrument (30, 44, and 70 GHz). The results provide constraints on the intrinsic brightness and the brightness time-variability of these planets. The majority of the planet flux density estimates are limited by systematic errors, but still yield better than 1% measurements in many cases. Applying data from Planck HFI, the Wilkinson Microwave Anisotropy Probe (WMAP), and the Atacama Cosmology Telescope (ACT) to a model that incorporates contributions from Saturn's rings to the planet's total flux density suggests a best fit value for the spectral index of Saturn's ring system of $\beta _\mathrm{ring} = 2.30\pm0.03$ over the 30-1000 GHz frequency range. The average ratio between the Planck-HFI measurements and the adopted model predictions for all five planets (excluding Jupiter observations for 353 GHz) is 0.997, 0.997, 1.018, and 1.032 for 100, 143, 217, and 353 GHz, respectively. Model predictions for planet thermodynamic temperatures are therefore consistent with the absolute calibration of Planck-HFI detectors at about the three-percent-level. We compare our measurements with published results from recent cosmic microwave background experiments. In particular, we observe that the flux densities measured by Planck HFI and WMAP agree to within 2%. These results allow experiments operating in the mm-wavelength range to cross-calibrate against Planck and improve models of radiative transport used in planetary science.

Recent Cosmic Microwave Background (CMB) temperature and polarization anisotropy measurements from the Planck mission have significantly improved previous constraints on the neutrino masses as well as the bounds on extended models with massless or massive sterile neutrino states. However, due to parameter degeneracies, additional low redshift priors are mandatory in order to sharpen the CMB neutrino bounds. We explore here the role of different priors on low redshift quantities, such as the Hubble constant, the cluster mass bias, and the reionization optical depth $\tau$. Concerning current priors on the Hubble constant and the cluster mass bias, the bounds on the neutrino parameters may differ appreciably depending on the choices adopted in the analyses. With regard to future improvements in the priors on the reionization optical depth, a value of $\tau=0.05\pm 0.01$, motivated by astrophysical estimates of the reionization redshift, would lead to $\sum m_\nu<0.0926$~eV at $90\%$~CL, when combining the full \textit{Planck} measurements, Baryon Acoustic Oscillation and Planck clusters data, thereby opening the window to unravel the neutrino mass hierarchy with existing cosmological probes.

We investigate constraints on cosmic reionization extracted from the Planck cosmic microwave background (CMB) data. We combine the Planck CMB anisotropy data in temperature with the low-multipole polarization data to fit LCDM models with various parameterizations of the reionization history. We obtain a Thomson optical depth tau=0.058 +/- 0.012 for the commonly adopted instantaneous reionization model. This confirms, with only data from CMB anisotropies, the low value suggested by combining Planck 2015 results with other data sets and also reduces the uncertainties. We reconstruct the history of the ionization fraction using either a symmetric or an asymmetric model for the transition between the neutral and ionized phases. To determine better constraints on the duration of the reionization process, we also make use of measurements of the amplitude of the kinetic Sunyaev-Zeldovich (kSZ) effect using additional information from the high resolution Atacama Cosmology Telescope and South Pole Telescope experiments. The average redshift at which reionization occurs is found to lie between z=7.8 and 8.8, depending on the model of reionization adopted. Using kSZ constraints and a redshift-symmetric reionization model, we find an upper limit to the width of the reionization period of Dz < 2.8. In all cases, we find that the Universe is ionized at less than the 10% level at redshifts above z~10. This suggests that an early onset of reionization is strongly disfavoured by the Planck data. We show that this result also reduces the tension between CMB-based analyses and constraints from other astrophysical sources.

Tram et al. 2016 recently pointed out that power-law inflation in presence of a dark radiation component may relieve the 3.3 sigma tension which exists within standard LCDM between the determination of the local value of the Hubble constant by Riess et al. (2016) and the value derived from CMB anisotropy data by the Planck collaboration. In this comment, we simply point out that this interesting proposal does not help in solving the $\sigma_8$ tension between the Planck data and, e.g., the weak lensing measurements. Moreover, when the latest constraints on the reionization optical depth obtained from Planck HFI data are included in the analysis, the $H_0$ tension reappears and this scenario looses appeal.

The recent determination of the local value of the Hubble constant by Riess et al, 2016 (hereafter R16) is now 3.3 sigma higher than the value derived from the most recent CMB anisotropy data provided by the Planck satellite in a LCDM model. Here we perform a combined analysis of the Planck and R16 results in an extended parameter space, varying simultaneously 12 cosmological parameters instead of the usual 6. We find that a phantom-like dark energy component, with effective equation of state $w=-1.29_{-0.12}^{+0.15}$ at 68 % c.l. can solve the current tension between the Planck dataset and the R16 prior in an extended $\Lambda$CDM scenario. On the other hand, the neutrino effective number is fully compatible with standard expectations. This result is confirmed when including cosmic shear data from the CFHTLenS survey and CMB lensing constraints from Planck. However, when BAO measurements are included we find that some of the tension with R16 remains, as also is the case when we include the supernova type Ia luminosity distances from the JLA catalog.

Using the Planck 2015 data release (PR2) temperature maps, we separate Galactic thermal dust emission from cosmic infrared background (CIB) anisotropies. For this purpose, we implement a specifically tailored component-separation method, the so-called generalized needlet internal linear combination (GNILC) method, which uses spatial information (the angular power spectra) to disentangle the Galactic dust emission and CIB anisotropies. We produce significantly improved all-sky maps of Planck thermal dust emission, with reduced CIB contamination, at 353, 545, and 857 GHz. By reducing the CIB contamination of the thermal dust maps, we provide more accurate estimates of the local dust temperature and dust spectral index over the sky with reduced dispersion, especially at high Galactic latitudes above $b = \pm 20{\deg}$. We find that the dust temperature is $T = (19.4 \pm 1.3)$ K and the dust spectral index is $\beta = 1.6 \pm 0.1$ averaged over the whole sky, while $T = (19.4 \pm 1.5)$ K and $\beta = 1.6 \pm 0.2$ on 21 % of the sky at high latitudes. Moreover, subtracting the new CIB-removed thermal dust maps from the CMB-removed Planck maps gives access to the CIB anisotropies over 60 % of the sky at Galactic latitudes $|b| > 20{\deg}$. Because they are a significant improvement over previous Planck products, the GNILC maps are recommended for thermal dust science. The new CIB maps can be regarded as indirect tracers of the dark matter and they are recommended for exploring cross-correlations with lensing and large-scale structure optical surveys. The reconstructed GNILC thermal dust and CIB maps are delivered as Planck products.

Parity violating extensions of the standard electromagnetic theory cause in vacuo rotation of the plane of polarization of propagating photons. This effect, also known as cosmic birefringence, impacts the cosmic microwave background (CMB) anisotropy angular power spectra, producing non-vanishing $T$--$B$ and $E$--$B$ correlations that are otherwise null when parity is a symmetry. Here we present new constraints on an isotropic rotation, parametrized by the angle $\alpha$, derived from Planck 2015 CMB polarization data. To increase the robustness of our analyses, we employ two complementary approaches, in harmonic space and in map space, the latter based on a peak stacking technique. The two approaches provide estimates for $\alpha$ that are in agreement within statistical uncertainties and very stable against several consistency tests. Considering the $T$--$B$ and $E$--$B$ information jointly, we find $\alpha = 0.31^{\circ} \pm 0.05^{\circ} \, ({\rm stat.})\, \pm 0.28^{\circ} \, ({\rm syst.})$ from the harmonic analysis and $\alpha = 0.35^{\circ} \pm 0.05^{\circ} \, ({\rm stat.})\, \pm 0.28^{\circ} \, ({\rm syst.})$ from the stacking approach. These constraints are compatible with no parity violation and are dominated by the systematic uncertainty in the orientation of Planck's polarization-sensitive bolometers.

We use the recent observations of Cosmic Microwave Background temperature and polarization anisotropies provided by the Planck satellite experiment to place constraints on the running $\alpha_\mathrm{s} = \mathrm{d}n_{\mathrm{s}} / \mathrm{d}\log k$ and the running of the running $\beta_{\mathrm{s}} = \mathrm{d}\alpha_{\mathrm{s}} / \mathrm{d}\log k$ of the spectral index $n_{\mathrm{s}}$ of primordial scalar fluctuations. We find $\alpha_\mathrm{s}=0.011\pm0.010$ and $\beta_\mathrm{s}=0.027\pm0.013$ at $68\%\,\mathrm{CL}$, suggesting the presence of a running of the running at the level of two standard deviations. We find no significant correlation between $\beta_{\mathrm{s}}$ and foregrounds parameters, with the exception of the point sources amplitude at $143\,\mathrm{GHz}$, $A^{PS}_{143}$, which shifts by half sigma when the running of the running is considered. We further study the cosmological implications of such preference for $\alpha_\mathrm{s},\beta_\mathrm{s}\sim0.01$ by including in the analysis the lensing amplitude $A_L$, the curvature parameter $\Omega_k$, and the sum of neutrino masses $\sum m_{\nu}$. We find that when the running of the running is considered, Planck data are more compatible with the standard expectations of $A_L = 1$ and $\Omega_k = 0$ but still hint at possible deviations. The indication for $\beta_\mathrm{s} > 0$ survives at two standard deviations when external datasets such as BAO and CFHTLenS are included in the analysis, and persists at $\sim 1.7$ standard deviations when CMB lensing is considered. We discuss the possibility of constraining $\beta_\mathrm{s}$ with current and future measurements of CMB spectral distortions, showing that an experiment like PIXIE could provide strong constraints on $\alpha_\mathrm{s}$ and $\beta_\mathrm{s}$.

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.

The characterization of the Galactic foregrounds has been shown to be the main obstacle in the challenging quest to detect primordial B-modes in the polarized microwave sky. We make use of the Planck-HFI 2015 data release at high frequencies to place new constraints on the properties of the polarized thermal dust emission at high Galactic latitudes. Here, we specifically study the spatial variability of the dust polarized spectral energy distribution, and its potential impact on the determination of the tensor-to-scalar ratio. We use the correlation ratio of the $C_\ell^{BB}$ angular power spectra between the 217- and 353-GHz channels as a tracer of these potential variations, computed on different high Galactic latitude regions, ranging from 80% to 20% of the sky. The new insight from Planck data is a departure of the correlation ratio from unity that cannot be attributed to a spurious decorrelation due to the cosmic microwave background, instrumental noise, or instrumental systematics. The effect is marginally detected on each region, but the statistical combination of all the regions gives more than 99% confidence for this variation in polarized dust properties. In addition, we show that the decorrelation increases when there is a decrease in the mean column density of the region of the sky being considered, and we propose a simple power-law empirical model for this dependence, which matches what is seen in the Planck data. We explore the effect that this measured decorrelation has on simulations of the BICEP2-Keck Array/Planck analysis and show that the 2015 constraints from those data still allow a decorrelation between the dust at 150 and 353GHz of the order of the one we measure. Finally we show that either spatial variation of the dust SED or of the dust polarization angle could produce decorrelations between 217- and 353-GHz data similar to those we observe in the data.

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.