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 h
igher $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.078pm0.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.
Full-sky CMB maps from the 2015 Planck release allow us to detect departures from global isotropy on the largest scales. We present the first searches using CMB polarization for correlations induced by a non-trivial topology with a fundamental domain
intersecting, or nearly intersecting, the last scattering surface (at comoving distance $chi_{rec}$). We specialize to flat spaces with toroidal and slab topologies, finding that explicit searches for the latter are sensitive to other topologies with antipodal symmetry. These searches yield no detection of a compact topology at a scale below the diameter of the last scattering surface. The limits on the radius $R_i$ of the largest sphere inscribed in the topological domain (at log-likelihood-ratio $Deltaln{L}>-5$ relative to a simply-connected flat Planck best-fit model) are $R_i>0.97chi_{rec}$ for the cubic torus and $R_i>0.56chi_{rec}$ for the slab. The limit for the cubic torus from the matched-circles search is numerically equivalent, $R_i>0.97chi_{rec}$ (99% CL) from polarisation data alone. We also perform a Bayesian search for a Bianchi VII$_h$ geometry. In the non-physical setting where the Bianchi cosmology is decoupled from the standard cosmology, Planck temperature data favour the inclusion of a Bianchi component. However, the cosmological parameters generating this pattern are in strong disagreement with those found from CMB anisotropy data alone. Fitting the induced polarization pattern for this model to Planck data requires an amplitude of $-0.1pm0.04$ compared to +1 if the model were to be correct. In the physical setting where the Bianchi parameters are fit simultaneously with the standard cosmological parameters, we find no evidence for a Bianchi VII$_h$ cosmology and constrain the vorticity of such models to $(omega/H)_0<7.6times10^{-10}$ (95% CL). [Abridged]
Planck CMB temperature maps allow detection of large-scale departures from homogeneity and isotropy. We search for topology with a fundamental domain nearly intersecting the last scattering surface (comoving distance $chi_r$). For most topologies stu
died the likelihood maximized over orientation shows some preference for multi-connected models just larger than $chi_r$. This effect is also present in simulated realizations of isotropic maps and we interpret it as the alignment of mild anisotropic correlations with chance features in a single realization; such a feature can also exist, in milder form, when the likelihood is marginalized over orientations. Thus marginalized, the limits on the radius $R_i$ of the largest sphere inscribed in a topological domain (at log-likelihood-ratio -5) are: in a flat Universe, $R_i>0.9chi_r$ for the cubic torus (cf. $R_i>0.9chi_r$ at 99% CL for a matched-circles search); $R_i>0.7chi_r$ for the chimney; $R_i>0.5chi_r$ for the slab; in a positively curved Universe, $R_i>1.0chi_r$ for the dodecahedron; $R_i>1.0chi_r$ for the truncated cube; $R_i>0.9chi_r$ for the octahedron. Similar limits apply to alternate topologies. We perform a Bayesian search for an anisotropic Bianchi VII$_h$ geometry. In a non-physical setting where the Bianchi parameters are decoupled from cosmology, Planck data favour a Bianchi component with a Bayes factor of at least 1.5 units of log-evidence: a Bianchi pattern is efficient at accounting for some large-scale anomalies in Planck data. However, the cosmological parameters are in strong disagreement with those found from CMB anisotropy data alone. In the physically motivated setting where the Bianchi parameters are fitted simultaneously with standard cosmological parameters, we find no evidence for a Bianchi VII$_h$ cosmology and constrain the vorticity of such models: $(omega/H)_0<8times10^{-10}$ (95% CL). [Abridged]
abridged] A method to rapidly estimate the Fourier power spectrum of a point distribution is presented. This method relies on a Taylor expansion of the trigonometric functions. It yields the Fourier modes from a number of FFTs, which is controlled by
the order N of the expansion and by the dimension D of the system. In three dimensions, for the practical value N=3, the number of FFTs required is 20. We apply the method to the measurement of the power spectrum of a periodic point distribution that is a local Poisson realization of an underlying stationary field. We derive explicit analytic expression for the spectrum, which allows us to quantify--and correct for--the biases induced by discreteness and by the truncation of the Taylor expansion, and to bound the unknown effects of aliasing of the power spectrum. We show that these aliasing effects decrease rapidly with the order N. The only remaining significant source of errors is reduced to the unavoidable cosmic/sample variance due to the finite size of the sample. The analytical calculations are successfully checked against a cosmological N-body experiment. We also consider the initial conditions of this simulation, which correspond to a perturbed grid. This allows us to test a case where the local Poisson assumption is incorrect. Even in that extreme situation, the third-order Fourier-Taylor estimator behaves well. We also show how to reach arbitrarily large dynamic range in Fourier space (i.e., high wavenumber), while keeping statistical errors in control, by appropriately folding the particle distribution.
