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 857 GHz. 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 Lambda CDM 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 Lambda CDM 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 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 High Frequency Instrument (HFI) data processing procedures for the Planck 2018 release. Major improvements in mapmaking have been achieved since the previous Planck 2015 release, many of which were used and described already in an intermediate paper dedicated to the Planck polarized data at low multipoles. These improvements enabled the first significant measurement of the reionization optical depth parameter using Planck-HFI data. This paper presents an extensive analysis of systematic effects, including the use of end-to-end 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, especially the leakage from intensity to polarization. Calibration, based on the cosmic microwave background (CMB) dipole, is now extremely accurate and in the frequency range 100-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 Planck HFI for future CMB experiments. The removal of bandpass leakage has been improved for the main high-frequency foregrounds 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. End-to-end simulations have been shown to 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 (analogue-to-digital convertor non-linearity residuals). Using these simulations, we have been able to 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.

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 degrees 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.

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., l less than or similar to 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 report on the implications for cosmic inflation of the 2018 release of the Planck cosmic microwave background (CMB) anisotropy measurements. The results are fully consistent with those reported using the data from 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(s)=0.9649 +/- 0.0042 at 68% CL. We find no evidence for a scale dependence of n(s), either as a running or as a running of the running. The Universe is found to be consistent with spatial flatness with a precision of 0.4% at 95% CL by combining Planck with a compilation of baryon acoustic oscillation 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 standard single-field inflationary models with Einstein gravity, these results imply that: (a) the predictions of slow-roll models with a concave potential, V<double prime>(phi) < 0, are increasingly favoured by the data; and (b) based on two different methods for reconstructing the inflaton potential, we find no evidence for dynamics beyond slow roll. Three different methods for the non-parametric reconstruction of the primordial power spectrum consistently confirm a pure power law in the range of comoving scales 0.005 Mpc(-1)<less than or similar to>k less than or similar to 0.2 Mpc(-1). A complementary analysis also finds no evidence for theoretically motivated parameterized features in the Planck power spectra. For the case of oscillatory features that are logarithmic or linear in k, this result is further strengthened by a new combined analysis including the Planck bispectrum data. The new Planck polarization data provide a stringent test of the adiabaticity of the initial conditions for the cosmological fluctuations. In correlated, mixed adiabatic and isocurvature models, the non-adiabatic contribution to the observed CMB temperature variance is constrained to 1.3%, 1.7%, and 1.7% at 95% CL for cold dark matter, neutrino density, and neutrino velocity, respectively. Planck power spectra plus lensing set constraints on the amplitude of compensated cold dark matter-baryon isocurvature perturbations that are consistent with current complementary measurements. The polarization data also provide improved constraints on inflationary models that predict a small statistically anisotropic quadupolar modulation of the primordial fluctuations. However, the polarization data do not support physical models for a scale-dependent dipolar modulation. All these findings support the key predictions of the standard single-field inflationary models, which will be further tested by future cosmological observations.

Observations of the submillimetre emission from Galactic dust, in both total intensity I and polarization, have received tremendous interest thanks to the Planck full-sky maps. In this paper we make use of such full-sky maps of dust polarized emission produced from the third public release of Planck data. As the basis for expanding on astrophysical studies of the polarized thermal emission from Galactic dust, we present full-sky maps of the dust polarization fraction p, polarization angle psi, and dispersion function of polarization angles ?. The joint distribution (one-point statistics) of p and N-H confirms that the mean and maximum polarization fractions decrease with increasing N-H. The uncertainty on the maximum observed polarization fraction, (max) = 22.0(-1.4)(+3.5) p max = 22 . 0 - 1.4 + 3.5 % at 353 GHz and 80 ' resolution, is dominated by the uncertainty on the Galactic emission zero level in total intensity, in particular towards diffuse lines of sight at high Galactic latitudes. Furthermore, the inverse behaviour between p and ? found earlier is seen to be present at high latitudes. This follows the ?proportional to p(-1) relationship expected from models of the polarized sky (including numerical simulations of magnetohydrodynamical turbulence) that include effects from only the topology of the turbulent magnetic field, but otherwise have uniform alignment and dust properties. Thus, the statistical properties of p, psi, and ? for the most part reflect the structure of the Galactic magnetic field. Nevertheless, we search for potential signatures of varying grain alignment and dust properties. First, we analyse the product map ?xp, looking for residual trends. While the polarization fraction p decreases by a factor of 3-4 between N-H=10(20) cm(-2) and N-H=2x10(22) cm(-2), out of the Galactic plane, this product ?xp only decreases by about 25%. Because ? is independent of the grain alignment efficiency, this demonstrates that the systematic decrease in p with N-H is determined mostly by the magnetic-field structure and not by a drop in grain alignment. This systematic trend is observed both in the diffuse interstellar medium (ISM) and in molecular clouds of the Gould Belt. Second, we look for a dependence of polarization properties on the dust temperature, as we would expect from the radiative alignment torque (RAT) theory. We find no systematic trend of ?xp with the dust temperature T-d, whether in the diffuse ISM or in the molecular clouds of the Gould Belt. In the diffuse ISM, lines of sight with high polarization fraction p and low polarization angle dispersion ? tend, on the contrary, to have colder dust than lines of sight with low p and high ?. We also compare the Planck thermal dust polarization with starlight polarization data in the visible at high Galactic latitudes. The agreement in polarization angles is remarkable, and is consistent with what we expect from the noise and the observed dispersion of polarization angles in the visible on the scale of the Planck beam. The two polarization emission-to-extinction ratios, R-P/p and R-S/V, which primarily characterize dust optical properties, have only a weak dependence on the column density, and converge towards the values previously determined for translucent lines of sight. We also determine an upper limit for the polarization fraction in extinction, p(V)/E(B-V), of 13% at high Galactic latitude, compatible with the polarization fraction p approximate to 20% observed at 353 GHz. Taken together, these results provide strong constraints for models of Galactic dust in diffuse gas.

It is a well known fact that galaxies are biased tracers of the distribution of matter in the Universe. The galaxy bias is usually factored as a function of redshift and scale, and approximated as being scale-independent on large, linear scales. In cosmologies with massive neutrinos, the galaxy bias defined with respect to the total matter field (cold dark matter, baryons, and non-relativistic neutrinos) also depends on the sum of the neutrino masses M-nu, and becomes scale-dependent even on large scales. This effect has been usually neglected given the sensitivity of current surveys. However, it becomes a severe systematic for future surveys aiming to provide the first detection of non-zero M-nu. The effect can be corrected for by defining the bias with respect to the density field of cold dark matter and baryons, rather than the total matter field. In this work, we provide a simple prescription for correctly mitigating the neutrino-induced scale-dependent bias effect in a practical way. We clarify a number of subtleties regarding how to properly implement this correction in the presence of redshift-space distortions and non-linear evolution of perturbations. We perform a Markov Chain Monte Carlo analysis on simulated galaxy clustering data that match the expected sensitivity of the Euclid survey. We find that the neutrino-induced scale-dependent bias can lead to important shifts in both the inferred mean value of M-nu, as well as its uncertainty, and provide an analytical explanation for the magnitude of the shifts. We show how these shifts propagate to the inferred values of other cosmological parameters correlated with M-nu, such as the cold dark matter physical density Omega(cdm)h(2) and the scalar spectral index n(s). In conclusion, we find that correctly accounting for the neutrino-induced scale-dependent bias will be of crucial importance for future galaxy clustering analyses. We encourage the cosmology community to correctly account for this effect using the simple prescription we present in our work. The tools necessary to easily correct for the neutrino-induced scale-dependent bias will be made publicly available in an upcoming release of the Boltzmann solver CLASS.

We explore cosmological constraints on the sum of the three active neutrino masses M-v in the context of dynamical dark energy (DDE) models with equation of state (EoS) parametrized as a function of redshift z by w(z) = w(0) + w(a)z/ (1 + z), and satisfying w(z) >= -1 for all z. We make use of cosmic microwave background data from the Planck satellite, baryon acoustic oscillation measurements, and supernovae la luminosity distance measurements, and perform a Bayesian analysis. We show that, within these models, the bounds on M-v do not degrade with respect to those obtained in the Lambda CDM case; in fact, the bounds arc slightly tighter, despite the enlarged parameter space. We explain our results based on the observation that, for fixed choices of w(0), w(a) such that w(z) >= -1 (but not w = -1 for all z), the upper limit on M-v is tighter than the Lambda CDM limit because of the well-known degeneracy between w and M-v. The Bayesian analysis we have carried out then integrates over the possible values of w(0)-w(a) such that w(z) >= -1, all of which correspond to tighter limits on M-v than the Lambda CDM limit. We find a 95% credible interval (C.I.) upper bound of M-v < 0.13 eV. This bound can be compared with the 95% C.I. upper bounds of M-v < 0.16 eV, obtained within the Lambda CDM model, and M-v < 0.41 eV, obtained in a DDE model with arbitrary EoS (which allows values of w < -1). Contrary to the results derived for DDE models with arbitrary EoS, we find that a dark energy component with w(z) >= -1 is unable to alleviate the tension between high-redshift observables and direct measurements of the Hubble constant H o . Finally, in light of the results of this analysis, we also discuss the implications for DDE models of a possible determination of the neutrino mass ordering by laboratory searches.

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 greater than or similar to 5 x 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 = (5 +/- 0.4) x 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 x 10(-3) after foreground cleaning and 60% de lensing. 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, l similar to 50. Finally, we report two sources of potential bias for the detection of the primordial B-modes by future CMB experiments: (i) the use of incorrect foreground models, e.g. a modelling error of Delta beta(s) = 0.02 on the synchrotron spectral indices may result in an excess in the recovered reionization peak corresponding to an effective Delta r > 10(-3); (ii) the average of the foreground line-of-sight spectral indices by the combined effects of pixelization and beam convolution, which adds an effective curvature to the foreground spectral energy distribution and may cause spectral degeneracies with the CMB in the frequency range probed by the experiment.