
Using the ${\it Planck}$ fullmission data, we present a detection of the
temperature (and therefore velocity) dispersion due to the kinetic
SunyaevZeldovich (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}$ foregroundcleaned `${\tt 2DILC\,}$' maps. By using
the Meta Catalogue of Xray 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 uniformweight temperature dispersion into a
measurement of the lineofsight 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 largescale structure studies, and provides direct
evidence of statistical homogeneity on scales of $600\,h^{1}{\rm Mpc}$. Our
study shows the promise of using crosscorrelations of the kSZ effect with
largescale structure in order to constrain the growth of structure.

The "Square Kilometre Array" (SKA) is a large international radio telescope
project characterised, as suggested by its name, by a total collecting area of
approximately one square kilometre, and consisting of several interferometric
arrays to observe at metric and centimetric wavelengths. The deployment of the
SKA will take place in two sites, in South Africa and Australia, and in two
successive phases. From its Phase 1, the SKA will be one of the most formidable
scientific machines ever deployed by mankind, and by far the most impressive in
terms of data throughput and required computing power. With the participation
of almost 200 authors from forty research institutes and six private companies,
the publication of this French SKA white paper illustrates the strong
involvement in the SKA project of the French astronomical community and of a
rapidly growing number of major scientific and technological players in the
fields of Big Data, high performance computing, energy production and storage,
as well as system integration.

This paper presents the Planck Multifrequency Catalogue of Nonthermal
(i.e., synchrotrondominated) Sources (PMNT) observed between 30 and 857 GHz by
ESA's Planck mission. This catalogue has been constructed by selecting objects
detected in the full mission allsky temperature maps at 30 and 143 GHz, with a
signaltonoise ratio (S/N)>3 in at least one of the two channels after
filtering with a particular Mexicanhat wavelet. As a result, 29400 source
candidates have been selected. Then, a multifrequency analysis has been
performed using the Matrix Filters methodology at the position of these
objects, and flux densities and errors have been calculated for all of them in
the nine Planck channels. This catalogue has been built using a different
methodology than the one adopted for the previous Planck Catalogues of Compact
Sources (PCCS and PCCS2), although the initial detection is done with the same
pipeline that was used to produce them. The present catalogue is the first
unbiased, fullsky catalogue of synchrotrondominated sources published at
millimetre and submillimetre wavelengths and constitutes a powerful database
for statistical studies of nonthermal extragalactic sources, whose emission is
dominated by the central active galactic nucleus. Together with the full
multifrequency catalogue, we also define the socalled Bright Planck
Multifrequency Catalogue of Nonthermal Sources (BPMNT), where only those
objects with a S/N>4 at both 30 and 143 GHz have been selected. In this
catalogue 1146 compact sources are detected outside the adopted Planck GAL070
mask. We also flag the socalled highsignificance subsample (HSPMNT), a
subset of 151 sources that are detected with S/N>4 in all nine Planck channels,
out of which 75 are found outside the Planck mask adopted here. The remaining
76 sources inside the Galactic mask are very likely Galactic objects.

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

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
powerlaw 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 singletemperature 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 multifrequency crossspectra, 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 suborbital experiments, then decorrelation might not be a problem
for CMB experiments aiming at a primordial Bmode detection limit on the
tensortoscalar ratio $r\simeq0.01$ at the recombination peak.

The dust properties inferred from the analysis of Planck observations in
total and polarized emission challenge current dust models. We propose new dust
models compatible with polarized and unpolarized data in extinction and
emission for translucent lines of sight ($0.5 < A_V < 2.5$). We amended the
DustEM tool to model polarized extinction and emission. We fit the spectral
dependence of the mean extinction, polarized extinction, SED and polarized SED
with PAHs, astrosilicates and amorphous carbon (aC). The astrosilicate
population is aligned along the magnetic field lines, while the aC population
may be aligned or not. With their current optical properties, oblate
astrosilicate grains are not emissive enough to reproduce the emission to
extinction polarization ratio $P_{353}/p_V$ derived with Planck data. Models
using prolate astrosilicate grains with an elongation $a/b=3$ and an inclusion
of 20% of porosity succeed. The spectral dependence of the polarized SED is
steeper in our models than in the data. Models perform slightly better when aC
grains are aligned. A small (6%) volume inclusion of aC in the astrosilicate
matrix removes the need for porosity and perfect grain alignment, and improves
the fit to the polarized SED. Dust models based on astrosilicates can be
reconciled with Planck data by adapting the shape of grains and adding
inclusions of porosity or aC in the astrosilicate matrix.

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

The six parameters of the standard $\Lambda$CDM model have bestfit values
derived from the Planck temperature power spectrum that are shifted somewhat
from the bestfit 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 cosmicvariance 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 bestfit values
between a WMAPlike large angularscale data set (with multipole moment
$\ell<800$ in the Planck temperature power spectrum) and an all angularscale
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]

The Draco nebula is a high Galactic latitude interstellar cloud likely to
have been formed by the collision of a Galactic halo cloud entering the disk of
the Milky Way. Such conditions are ideal to study the formation of cold and
dense gas in colliding flows of warm gas. We present HerschelSPIRE
observations that reveal the fragmented structure of the interface between the
infalling cloud and the Galactic layer. This front is characterized by a
RayleighTaylor instability structure. From the determination of the typical
length of the periodic structure (2.2 pc) we estimated the gas kinematic
viscosity and the turbulence dissipation scale (0.1 pc) that is compatible with
that expected if ambipolar diffusion is the main mechanism of energy
dissipation in the WNM. The smallscale structures of the nebula are typical of
that seen in some molecular clouds. The gas density has a lognormal
distribution with an average value of $10^3$ cm$^{3}$. The size of the
structures is 0.10.2 pc but this estimate is limited by the resolution of the
observations. The mass ranges from 0.2 to 20 M$_{\odot}$ and the distribution
of the more massive clumps follows a power law $dN/d\log(M) \sim M^{1.4}$. We
identify a masssize relation with the same exponent as that found in GMCs
($M\sim L^{2.3}$) but only 15% of the mass of the cloud is in gravitationally
bound structures. We conclude that the increase of pressure in the collision is
strong enough to trigger the WNMCNM transition caused by the interplay between
turbulence and thermal instability as selfgravity is not dominating the
dynamics.

Measurements of flux density are described for five planets, Mars, Jupiter,
Saturn, Uranus, and Neptune, across the six Planck High Frequency Instrument
frequency bands (100857 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 timevariability 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 301000 GHz frequency range. The average ratio between
the PlanckHFI 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 PlanckHFI detectors at about the threepercentlevel. 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 mmwavelength range to crosscalibrate against
Planck and improve models of radiative transport used in planetary science.

Observations of Galactic H I gas for seven targeted regions at intermediate
Galactic latitude are presented at 1' angular resolution using data from the
DRAO Synthesis Telescope (ST) and the Green Bank Telescope (GBT). The DHIGLS
data are the most extensive arcminute resolution measurements of the diffuse
atomic interstellar medium beyond those in the Galactic plane. The acquisition,
reduction, calibration, and mosaicking of the DRAO ST data and the cross
calibration and incorporation of the shortspacing information from the GBT are
described. The high quality of the resulting DHIGLS products enables a variety
of new studies in directions of low Galactic column density. We analyze the
angular power spectra of maps of the integrated H I emission (column density)
from the data cubes for several distinct velocity ranges. Fitting power
spectrum models based on a power law, but including the effects of the
synthesized beam and noise at high spatial frequencies, we find exponents
ranging from 2.5 to 3.0.. Power spectra of maps of the centroid velocity for
these components give similar results. These exponents are interpreted as being
representative of the 3D density and 3D velocity fields of the atomic gas,
respectively. We find evidence for dramatic changes in the H I structures in
channel maps over even small changes in velocity. This narrow line emission has
counterparts in absorption spectra against bright background radio sources,
quantifying that the gas is cold and dense and can be identified as the cold
neutral medium phase. Fully reduced DHIGLS H I data cubes and other data
products are available at www.cita.utoronto.ca/DHIGLS.

The secondary cosmic microwave background (CMB) $B$modes stem from the
postdecoupling distortion of the polarization $E$modes due to the
gravitational lensing effect of largescale structures. These lensinginduced
$B$modes constitute both a valuable probe of the dark matter distribution and
an important contaminant for the extraction of the primary CMB $B$modes from
inflation. Planck provides accurate nearly allsky measurements of both the
polarization $E$modes and the integrated mass distribution via the
reconstruction of the CMB lensing potential. By combining these two data
products, we have produced an allsky template map of the lensinginduced
$B$modes using a realspace algorithm that minimizes the impact of sky masks.
The crosscorrelation of this template with an observed (primordial and
secondary) $B$mode map can be used to measure the lensing $B$mode power
spectrum at multipoles up to $2000$. In particular, when crosscorrelating with
the $B$mode contribution directly derived from the Planck polarization maps,
we obtain lensinginduced $B$mode power spectrum measurement at a significance
level of $12\,\sigma$, which agrees with the theoretical expectation derived
from the Planck bestfit $\Lambda$CDM model. This unique nearly allsky
secondary $B$mode template, which includes the lensinginduced information
from intermediate to small ($10\lesssim \ell\lesssim 1000$) angular scales, is
delivered as part of the Planck 2015 public data release. It will be
particularly useful for experiments searching for primordial $B$modes, such as
BICEP2/Keck Array or LiteBIRD, since it will enable an estimate to be made of
the lensinginduced contribution to the measured total CMB $B$modes.

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

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 lowmultipole 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 SunyaevZeldovich (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 redshiftsymmetric 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 CMBbased analyses and constraints from other astrophysical
sources.

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
componentseparation method, the socalled 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 allsky 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 CIBremoved thermal dust maps from the
CMBremoved 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 crosscorrelations with
lensing and largescale 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 nonvanishing
$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 polarizationsensitive bolometers.

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

This paper presents the Planck 2015 likelihoods, statistical descriptions of
the 2point correlations of CMB data, using the hybrid approach employed
previously: pixelbased 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 endtoend
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 CMBforeground separation. Within the baseline model
this requires a reionization optical depth $\tau=0.078\pm0.019$, significantly
lower than without highfrequency 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.

Diffuse Galactic light has been observed in the optical since the 1930s. We
propose that, when observed in the optical with deep imaging surveys, it can be
used as a tracer of the turbulent cascade in the diffuse interstellar medium
(ISM), down to scales of about 1 arcsec. Here we present a power spectrum
analysis of the dust column density of a diffuse cirrus at high Galactic
latitude (l ~ 198 deg, b ~ 32 deg) as derived from the combination of a MegaCam
gband image, obtained as part of the MATLAS large programme at the CFHT, with
Planck radiance and Wise 12 micron data. The combination of these three
datasets have allowed us to compute the density power spectrum of the HI over
scales of more than three orders of magnitude. We found that the density field
is well described by a single power law over scales ranging from 0.01 to 50 pc.
The exponent of the power spectrum, gamma=2.9 + 0.1, is compatible with what
is expected for thermally bistable and turbulent HI. We did not find any
steepening of the power spectrum at small scales indicating that the typical
scale at which turbulent energy is dissipated in this medium is smaller than
0.01pc. The ambipolar diffusion scenario that is usually proposed as the main
dissipative agent, is consistent with our data only if the density of the cloud
observed is higher than the typical values assumed for the cold neutral medium
gas. We discuss the new avenue offered by deep optical imaging surveys for the
study of the low density ISM structure and turbulence.

The characterization of the Galactic foregrounds has been shown to be the
main obstacle in the challenging quest to detect primordial Bmodes in the
polarized microwave sky. We make use of the PlanckHFI 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 tensortoscalar ratio. We use
the correlation ratio of the $C_\ell^{BB}$ angular power spectra between the
217 and 353GHz 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 powerlaw
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 BICEP2Keck 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 353GHz data similar to those we observe in the
data.

We present results based on fullmission Planck observations of temperature
and polarization anisotropies of the CMB. These data are consistent with the
sixparameter 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 tensortoscalar ratio of r <0.11 consistent
with the Bmode 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.

Continuum spectra covering centimetre to submillimetre wavelengths are
presented for a northern sample of 104 extragalactic radio sources, mainly
active galactic nuclei, based on fourepoch Planck data. The nine Planck
frequencies, from 30 to 857 GHz, are complemented by a set of simultaneous
groundbased radio observations between 1.1 and 37 GHz. The singlesurvey
Planck data confirm that the flattest highfrequency radio spectral indices are
close to zero, indicating that the original accelerated electron energy
spectrum is much harder than commonly thought, with powerlaw index around 1.5
instead of the canonical 2.5. The radio spectra peak at high frequencies and
exhibit a variety of shapes. For a small set of lowz sources, we find a
spectral upturn at high frequencies, indicating the presence of intrinsic cold
dust. Variability can generally be approximated by achromatic variations, while
sources with clear signatures of evolving shocks appear to be limited to the
strongest outbursts.

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

Fullsky 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 nontrivial 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 loglikelihoodratio $\Delta\ln{L}>5$
relative to a simplyconnected flat Planck bestfit model) are
$R_i>0.97\chi_{rec}$ for the cubic torus and $R_i>0.56\chi_{rec}$ for the slab.
The limit for the cubic torus from the matchedcircles search is numerically
equivalent, $R_i>0.97\chi_{rec}$ (99% CL) from polarisation data alone. We also
perform a Bayesian search for a Bianchi VII$_h$ geometry. In the nonphysical
setting where the Bianchi cosmology is decoupled from the standard cosmology,
Planck temperature data favour the inclusion of a Bianchi component. However,
the cosmological parameters generating this pattern are in strong disagreement
with those found from CMB anisotropy data alone. Fitting the induced
polarization pattern for this model to Planck data requires an amplitude of
$0.1\pm0.04$ compared to +1 if the model were to be correct. In the physical
setting where the Bianchi parameters are fit simultaneously with the standard
cosmological parameters, we find no evidence for a Bianchi VII$_h$ cosmology
and constrain the vorticity of such models to $(\omega/H)_0<7.6\times10^{10}$
(95% CL). [Abridged]

Within four nearby (d < 160 pc) molecular clouds, we statistically evaluate
the structure of the interstellar magnetic field, projected on the plane of the
sky and integrated along the line of sight, as inferred from the polarized
thermal emission of Galactic dust observed by Planck at 353 GHz and from the
optical and NIR polarization of background starlight. We compare the dispersion
of the field orientation directly in vicinities with an area equivalent to that
subtended by the Planck effective beam at 353 GHz (10') and using the
secondorder structure functions of the field orientation angles. We find that
the average dispersion of the starlightinferred field orientations within
10'diameter vicinities is less than 20 deg, and that at these scales the mean
field orientation is on average within 5 deg of that inferred from the
submillimetre polarization observations in the considered regions. We also find
that the dispersion of starlight polarization orientations and the polarization
fractions within these vicinities are well reproduced by a Gaussian model of
the turbulent structure of the magnetic field, in agreement with the findings
reported by the Planck collaboration at scales greater than 10' and for
comparable column densities. At scales greater than 10', we find differences of
up to 14.7 deg between the secondorder structure functions obtained from
starlight and submillimetre polarization observations in the same positions in
the plane of the sky, but comparison with a Gaussian model of the turbulent
structure of the magnetic field indicates that these differences are small and
are consistent with the difference in angular resolution between both
techniques.