
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.

Correlations between cosmic microwave background (CMB) temperature,
polarization and spectral distortion anisotropies can be used as a probe of
primordial nonGaussianity. Here, we perform a reconstruction of
$\mu$distortion anisotropies in the presence of Galactic and extragalactic
foregrounds, applying the socalled Constrained ILC component separation method
to simulations of proposed CMB space missions (PIXIE, LiteBIRD, CORE, PICO).
Our sky simulations include Galactic dust, Galactic synchrotron, Galactic
freefree, thermal SunyaevZeldovich effect, as well as primary CMB temperature
and $\mu$distortion anisotropies, the latter being added as correlated field.
The Constrained ILC method allows us to null the CMB temperature anisotropies
in the reconstructed $\mu$map (and vice versa), in addition to mitigating the
contaminations from astrophysical foregrounds and instrumental noise. We
compute the crosspower spectrum between the reconstructed (CMBfree)
$\mu$distortion map and the ($\mu$free) CMB temperature map, after foreground
removal and component separation. Since the crosspower spectrum is
proportional to the primordial nonGaussianity parameter, $f_{\rm NL}$, on
scales $k\simeq 740$ Mpc$^{1}$, this allows us to derive $f_{\rm
NL}$detection limits for the aforementioned future CMB experiments. Our
analysis shows that foregrounds degrade the theoretical detection limits (based
mostly on instrumental noise) by more than one order of magnitude, with PICO
standing the best chance at placing upper limits on scaledependent
nonGaussianity. We also discuss the dependence of the constraints on the
channel sensitivities and chosen bands. Like for $B$mode polarization
measurements, extended coverage at frequencies $\nu\lesssim 40\,{\rm GHz}$ and
$\nu\gtrsim 400\,{\rm GHz}$ provides more leverage than increased channel
sensitivity.

Inflation is the leading theory of the first instant of the universe.
Inflation, which postulates that the universe underwent a period of rapid
expansion an instant after its birth, provides convincing explanation for
cosmological observations. Recent advancements in detector technology have
opened opportunities to explore primordial gravitational waves generated by the
inflation through Bmode (divergentfree) polarization pattern embedded in the
Cosmic Microwave Background anisotropies. If detected, these signals would
provide strong evidence for inflation, point to the correct model for
inflation, and open a window to physics at ultrahigh energies.
LiteBIRD is a satellite mission with a goal of detecting
degreeandlargerangularscale Bmode polarization. LiteBIRD will observe at
the second Lagrange point with a 400 mm diameter telescope and 2,622 detectors.
It will survey the entire sky with 15 frequency bands from 40 to 400 GHz to
measure and subtract foregrounds.
The U.S. LiteBIRD team is proposing to deliver subKelvin instruments that
include detectors and readout electronics. A lensletcoupled sinuous antenna
array will cover lowfrequency bands (40 GHz to 235 GHz) with four frequency
arrangements of trichroic pixels. An orthomodetransducercoupled corrugated
horn array will cover highfrequency bands (280 GHz to 402 GHz) with three
types of single frequency detectors. The detectors will be made with Transition
Edge Sensor (TES) bolometers cooled to a 100 milliKelvin base temperature by
an adiabatic demagnetization refrigerator.The TES bolometers will be read out
using digital frequency multiplexing with Superconducting QUantum Interference
Device (SQUID) amplifiers. Up to 78 bolometers will be multiplexed with a
single SQUID amplidier.
We report on the subKelvin instrument design and ongoing developments for
the LiteBIRD mission.

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.

HI intensity mapping is a new observational technique to survey the
largescale structure of matter using the 21 cm emission line of atomic
hydrogen (HI). In this work, we simulate BINGO (BAO from Integrated Neutral Gas
Observations) and SKA (Square Kilometre Array) phase1 dish array operating in
autocorrelation mode. For the optimal case of BINGO with no foregrounds, the
combination of the HI angular power spectra with Planck results allows $w$ to
be measured with a precision of $4\%$, while the combination of the BAO
acoustic scale with Planck gives a precision of $7\%$. We consider a number of
potentially complicating effects, including foregrounds and redshift dependent
bias, which increase the uncertainty on $w$ but not dramatically; in all cases
the final uncertainty is found to be $\Delta w < 8\%$ for BINGO. For the
combination of SKAMID in autocorrelation mode with Planck, we find that, in
ideal conditions, $w$ can be measured with a precision of $4\%$ for the
redshift range $0.35 < z < 3$ (i.e., for the bandwidth of $\Delta \nu = [350,
1050]$ MHz) and $2\%$ for $0 < z < 0.49$ (i.e., $\Delta \nu = [950, 1421]$
MHz). Extending the model to include the sum of neutrino masses yields a $95\%$
upper limit of $\sum m_\nu < 0.24$ eV for BINGO and $\sum m_\nu < 0.08$ eV for
SKA phase 1, competitive with the current best constraints in the case of BINGO
and significantly better than them in the case of SKA.

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 investigate the performance of a simple Bayesian fitting approach to
correct the cosmic microwave background (CMB) Bmode polarization for
gravitational lensing effects in the recovered probability distribution of the
tensortoscalar ratio. We perform a twodimensional power spectrum fit of the
amplitude of the primordial Bmodes (tensortoscalar ratio, $r$) and the
amplitude of the lensing Bmodes (parameter $A_{lens}$), jointly with the
estimation of the astrophysical foregrounds including both synchrotron and
thermal dust emissions. Using this Bayesian framework, we forecast the ability
of the proposed CMB space mission LiteBIRD to constrain $r$ in the presence of
realistic lensing and foreground contributions. We compute the joint posterior
distribution of $r$ and $A_{lens}$, which we improve by adopting a prior on
$A_{lens}$ taken from the South Pole Telescope (SPT) measurement. As it applies
to the power spectrum, this approach cannot mitigate the uncertainty on $r$
that is due to Emode cosmic variance transferred to Bmodes by lensing, unlike
standard delensing techniques that are performed on maps. However, the method
allows to correct for the bias on $r$ induced by lensing, at the expense of a
larger uncertainty due to the increased volume of the parameter space. We
quantify, for different values of the tensortoscalar ratio, the tradeoff
between bias correction and increase of uncertainty on $r$. For LiteBIRD
simulations, which include foregrounds and lensing contamination, we find that
correcting the foregroundcleaned CMB Bmode power spectrum for the lensing
bias, not the lensing cosmic variance, still guarantees a $3\sigma$ detection
of $r=5\times 10^{3}$. The significance of the detection is increased to
$6\sigma$ when the current SPT prior on $A_{lens}$ is adopted.

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 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]

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.

This paper presents a study of the ISW effect from the Planck 2015
temperature and polarization data release. The CMB is crosscorrelated with
different LSS tracers: the NVSS, SDSS and WISE catalogues, and the Planck 2015
lensing map. This crosscorrelation yields a detection at $4\,\sigma$, where
most of the signaltonoise is due to the Planck lensing and NVSS. In fact, the
ISW effect is detected only from the Planck data (through the ISWlensing
bispectrum) at $\approx 3\,\sigma$, which is similar to the detection level
achieved by combining the crosscorrelation signal coming from all the
catalogues. The ISW signal allow us to detect $\Omega_\Lambda$ at more than
$3\,\sigma$. This crosscorrelation analysis is performed only with the Planck
temperature data, since the polarization scales available in the 2015 release
do not permit significant improvement of the CMBLSS crosscorrelation
detectability. Nevertheless, polarization data is used to study the anomalously
large ISW signal previously reported through the stacking of CMB features at
the locations of known superstructures. We find that the current Planck
polarization data do not exclude that this signal could be caused by the ISW
effect. In addition, the stacking of the Planck lensing map on the locations of
superstructures exhibits a positive crosscorrelation with these largescale
structures. Finally, we have improved our previous reconstruction of the ISW
temperature fluctuations by combining the information encoded in all the
previously mentioned LSS tracers. In particular, we construct a map of the ISW
secondary anisotropies and the corresponding uncertainties map, obtained from
simulations. We also explore the reconstruction of the ISW anisotropies caused
by the LSS traced by the 2MPZ survey by directly inverting the density field
into the gravitational potential field.

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]