
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.

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.

QUIJOTE (QUI JOint TEnerife) is an experiment designed to achieve CMB
Bmode polarization detection and sensitive enough to detect a primordial
gravitationalwave component if the Bmode amplitude is larger than r = 0.05.
It consists in two telescopes and three instruments observing in the frequency
range 1042 GHz installed at the Teide Observatory in the Canary Islands,
Spain. The observing strategy includes three raster scan deep integration
fields for cosmology, a nominal wide survey covering the Northen Sky and
specific raster scan deep integration observations in regions of specific
interest. The main goals of the project are presented and the first scientific
results obtained with the first instrument are reviewed.

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.

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.

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
crosscorrelations will enable a more robust determination of purely
velocitydriven 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 cosmicvariancelimited experiment up
to a multipole l of 2000. Precise interfrequency 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 COBEFIRAS 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
timelinetomap simulations to verify that the CORE instrumental setup 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 intensitytopolarization 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 finegrained 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.

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 tensortoscalar ratio values of ${r\gtrsim 5\times 10^{3}}$. We
consider detailed sky simulations based on stateoftheart CMB observations
that consist of CMB polarization with $\tau=0.055$ and tensortoscalar 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 tensortoscalar
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 pointsources can be the
dominant foreground contamination to primordial Bmodes 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 mediumsized 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 nearultimate
CMB polarisation mission which, for optimal complementarity with groundbased
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 spaceborne, multiband, multibeam 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 mediumsize space mission within the Cosmic Vision Programme
of the European Space Agency. The instrument has no moving parts, and uses
arrays of diffractionlimited 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 spaceborne 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 Mediumsize 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 highz galaxies distributed over the full sky. The
extreme brightness of these galaxies will make it possible to study them, via
followup observations, in extraordinary detail. Also, the CORE resolution
matches the typical sizes of highz galaxy protoclusters 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 submm
wavelengths, respectively.

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]

We present a comprehensive study of the regularity of the covariance matrix
of a discretized field on the sphere. In a particular situation, the rank of
the matrix depends on the number of pixels, the number of spherical harmonics,
the symmetries of the pixelization scheme and the presence of a mask. Taking
into account the above mentioned components, we provide analytical expressions
that constrain the rank of the matrix. They are obtained by expanding the
determinant of the covariance matrix as a sum of determinants of matrices made
up of spherical harmonics. We investigate these constraints for five different
pixelizations that have been used in the context of Cosmic Microwave Background
(CMB) data analysis: Cube, Icosahedron, Igloo, GLESP and HEALPix, finding that,
at least in the considered cases, the HEALPix pixelization tends to provide a
covariance matrix with a rank closer to the maximum expected theoretical value
than the other pixelizations. The effect of the propagation of numerical errors
in the regularity of the covariance matrix is also studied for different
computational precisions, as well as the effect of adding a certain level of
noise in order to regularize the matrix. In addition, we investigate the
application of the previous results to a particular example that requires the
inversion of the covariance matrix: the estimation of the CMB temperature power
spectrum through the Quadratic Maximum Likelihood algorithm. Finally, some
general considerations in order to achieve a regular covariance matrix are also
presented.

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 groundbased 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 SunyaevZeldovich 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 fluxmass relation; assuming selfsimilar evolution, CORE will
detect $\sim 500$ clusters at redshift $z>1.5$. This changes to 800 (200) when
increasing (decreasing) the mirror size by 30 cm. CORE will be able to measure
individual cluster halo masses through lensing of the cosmic microwave
background anisotropies with a 1$\sigma$ sensitivity of $4\times10^{14}
M_\odot$, for a 120 cm aperture telescope, and $10^{14} M_\odot$ for a 180 cm
one. [abridged]

We study the Planck CMB temperature at different scales through its
derivatives up to second order, which allows one to characterize the local
shape and isotropy of the field. The problem of having an incomplete sky in the
calculation and statistical characterization of the derivatives is addressed in
the paper. The analysis confirms the existence of a low variance in the CMB at
large scales, which is also noticeable in the derivatives. Moreover, deviations
from the standard model in the gradient, curvature and the eccentricity tensor
are studied in terms of extreme values on the data. As it is expected, the Cold
Spot is detected as one of the most prominent peaks in terms of curvature, but
additionally, when the information of the temperature and its Laplacian are
combined, another feature with similar probability at the scale of $10^\circ$
is also observed. However, the $p$value of these two deviations increase above
the $6\%$ when they are referred to the variance calculated from the
theoretical fiducial model, indicating that these deviations can be associated
to the low variance anomaly. Finally, an estimator of the directional
anisotropy for spinorial quantities is introduced, which is applied to the
spinors derived from the field derivatives. An anisotropic direction whose
probability is $<1\%$ is detected in the eccentricity tensor.

In the present work, we study the largest structures of the CMB temperature
measured by Planck in terms of the most prominent peaks on the sky, which, in
particular, are located in the southern galactic hemisphere. Besides these
largescale features, the wellknown Cold Spot anomaly is included in the
analysis. All these peaks would contribute significantly to some of the CMB
largescale anomalies, as the parity and hemispherical asymmetries, the dipole
modulation, the alignment between the quadrupole and the octopole, or in the
case of the Cold Spot, to the nonGaussianity of the field. The analysis of the
peaks is performed by using their multipolar profiles, which characterize the
local shape of the peaks in terms of the discrete Fourier transform of the
azimuthal angle. In order to quantify the local anisotropy of the peaks, the
distribution of the phases of the multipolar profiles is studied by using the
Rayleigh random walk methodology. Finally, a direct analysis of the
2dimensional field around the peaks is performed in order to take into account
the effect of the galactic mask. The results of the analysis conclude that,
once the peak amplitude and its first and second order derivatives at the
centre are conditioned, the rest of the field is compatible with the standard
model. In particular, it is observed that the Cold Spot anomaly is caused by
the large value of curvature at the centre.

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.

We present QUI JOint TEnerife (QUIJOTE) intensity and polarisation maps at
1020 GHz covering a region along the Galactic plane 24<l<45 deg, b<8 deg.
These maps result from 210 h of data, have a sensitivity in polarisation of ~40
muK/beam and an angular resolution of ~1 deg. Our intensity data are crucial to
confirm the presence of anomalous microwave emission (AME) towards the two
molecular complexes W43 (22 sigma) and W47 (8 sigma). We also detect at high
significance (6 sigma) AME associated with W44, the first clear detection of
this emission towards a SNR. The new QUIJOTE polarisation data, in combination
with WMAP, are essential to: i) Determine the spectral index of the synchrotron
emission in W44, beta_sync =0.62 +/0.03, in good agreement with the value
inferred from the intensity spectrum once a freefree component is included in
the fit. ii) Trace the change in the polarisation angle associated with Faraday
rotation in the direction of W44 with rotation measure 404 +/ 49 rad/m2. And
iii) set upper limits on the polarisation of W43 of Pi_AME <0.39 per cent (95
per cent C.L.) from QUIJOTE 17~GHz, and <0.22 per cent from WMAP 41 GHz data,
which are the most stringent constraints ever obtained on the polarisation
fraction of the AME. For typical physical conditions (grain temperature and
magnetic field strengths), and in the case of perfect alignment between the
grains and the magnetic field, the models of electric or magnetic dipole
emissions predict higher polarisation fractions.

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]