We place constraints on sterile neutrino resonant production (RP) and scalar decay production (SDP) mechanisms assuming that sterile neutrino represents a fraction from the total Cold Dark Matter energy density. For the cosmological analysis, we complement the CMB anisotropies measurements with CMB lensing gravitational potential measurements, that are sensitive to the DM distribution out to high redshifts and with the cosmic shear data, that constraints the gravitational potential at lower redshifts than CMB. We show that our datasets have enough sensitivity to constrain the sterile neutrino mass and mass fraction inside the co-moving free-streaming horizon in both RP and SDP scenarios. For RP case we find that the best fit values of sterile neutrino mass and mixing angle are in the parameter space of interest for sterile neutrino DM decay interpretation of the 3.5 keV X-ray line with a DM mass fraction f_S=0.28 \pm 0.3 (at 68% CL) that excludes the assumption of sterile neutrinos as being all of the DM. For SDP case we find f_S=0.86 \pm 0.07 (at 68% CL), in agreement with the upper limit constraint on f_S from the X-ray non-detection and Ly-alpha forest measurements. The sterile neutrino mass predicted by both RP and SDP models are consistent within 0.3 \sigma.

We present cluster counts and corresponding cosmological constraints from the Planck full mission data set. Our catalogue consists of 439 clusters detected via their Sunyaev-Zeldovich (SZ) signal down to a signal-to-noise 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 two-dimensional distribution in redshift and signal-to-noise. 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, $(1-b)$. 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 present-day 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 $(1-b)$ 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 non-minimal 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.

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 tensor-to-scalar ratio is $r_{0.002} < 0.11$ (95% CL), consistent with the B-mode 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 tensor-to-scalar 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.

This paper presents a study of the ISW effect from the Planck 2015 temperature and polarization data release. The CMB is cross-correlated with different LSS tracers: the NVSS, SDSS and WISE catalogues, and the Planck 2015 lensing map. This cross-correlation yields a detection at $4\,\sigma$, where most of the signal-to-noise is due to the Planck lensing and NVSS. In fact, the ISW effect is detected only from the Planck data (through the ISW-lensing bispectrum) at $\approx 3\,\sigma$, which is similar to the detection level achieved by combining the cross-correlation signal coming from all the catalogues. The ISW signal allow us to detect $\Omega_\Lambda$ at more than $3\,\sigma$. This cross-correlation 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 CMB-LSS cross-correlation 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 cross-correlation with these large-scale 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.

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 full-mission release. Using a polarization-only estimator we detect lensing at a significance of 5 sigma. We cross-check 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 best-fitting LCDM model based on the Planck temperature and polarization power spectra. Using the lensing likelihood alone we obtain a percent-level 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 E-mode polarization also measured by Planck to generate an estimate of the lensing B-mode. We show that this lensing B-mode estimate is correlated with the B-modes 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 large-scale temperature anisotropies, detecting a cross-correlation at the 3 sigma level, as expected due to dark energy in the concordance LCDM model.

The Planck full mission cosmic microwave background(CMB) temperature and E-mode polarization maps are analysed to obtain constraints on primordial non-Gaussianity(NG). Using three classes of optimal bispectrum estimators - separable template-fitting (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 cross-validation 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 time-domain de-glitching 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 model-independent reconstructions of the CMB bispectrum and derive constraints on early universe scenarios that generate NG, including general single-field and axion inflation, initial state modifications, parity-violating tensor bispectra, and directionally dependent vector models. We also present a wide survey of scale-dependent 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]

We present results based on full-mission Planck observations of temperature and polarization anisotropies of the CMB. These data are consistent with the six-parameter 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 tensor-to-scalar ratio of r <0.11 consistent with the B-mode 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 four-epoch Planck data. The nine Planck frequencies, from 30 to 857 GHz, are complemented by a set of simultaneous ground-based radio observations between 1.1 and 37 GHz. The single-survey Planck data confirm that the flattest high-frequency radio spectral indices are close to zero, indicating that the original accelerated electron energy spectrum is much harder than commonly thought, with power-law 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 low-z 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.

Full-sky CMB maps from the 2015 Planck release allow us to detect departures from global isotropy on the largest scales. We present the first searches using CMB polarization for correlations induced by a non-trivial topology with a fundamental domain intersecting, or nearly intersecting, the last scattering surface (at comoving distance $\chi_{rec}$). We specialize to flat spaces with toroidal and slab topologies, finding that explicit searches for the latter are sensitive to other topologies with antipodal symmetry. These searches yield no detection of a compact topology at a scale below the diameter of the last scattering surface. The limits on the radius $R_i$ of the largest sphere inscribed in the topological domain (at log-likelihood-ratio $\Delta\ln{L}>-5$ relative to a simply-connected flat Planck best-fit 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 matched-circles 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 non-physical setting where the Bianchi cosmology is decoupled from the standard cosmology, Planck temperature data favour the inclusion of a Bianchi component. However, the cosmological parameters generating this pattern are in strong disagreement with those found from CMB anisotropy data alone. Fitting the induced polarization pattern for this model to Planck data requires an amplitude of $-0.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]

We study the implications of Planck data for models of dark energy (DE) and modified gravity (MG), beyond the cosmological constant scenario. We start with cases where the DE only directly affects the background evolution, considering Taylor expansions of the equation of state, principal component analysis and parameterizations related to the potential of a minimally coupled DE scalar field. When estimating the density of DE at early times, we significantly improve present constraints. We then move to general parameterizations of the DE or MG perturbations that encompass both effective field theories and the phenomenology of gravitational potentials in MG models. Lastly, we test a range of specific models, such as k-essence, f(R) theories and coupled DE. In addition to the latest Planck data, for our main analyses we use baryonic acoustic oscillations, type-Ia supernovae and local measurements of the Hubble constant. We further show the impact of measurements of the cosmological perturbations, such as redshift-space distortions and weak gravitational lensing. These additional probes are important tools for testing MG models and for breaking degeneracies that are still present in the combination of Planck and background data sets. All results that include only background parameterizations are in agreement with LCDM. When testing models that also change perturbations (even when the background is fixed to LCDM), some tensions appear in a few scenarios: the maximum one found is \sim 2 sigma for Planck TT+lowP when parameterizing observables related to the gravitational potentials with a chosen time dependence; the tension increases to at most 3 sigma when external data sets are included. It however disappears when including CMB lensing.

We compute and investigate four types of imprint of a stochastic background of primordial magnetic fields (PMFs) on the cosmic microwave background (CMB) anisotropies: the impact of PMFs on the CMB spectra; the effect on CMB polarization induced by Faraday rotation; the impact of PMFs on the ionization history; magnetically-induced non-Gaussianities; and the magnetically-induced breaking of statistical isotropy. Overall, Planck data constrain the amplitude of PMFs to less than a few nanogauss. In particular, individual limits coming from the analysis of the CMB angular power spectra, using the Planck likelihood, are $B_{1\,\mathrm{Mpc}}< 4.4$ nG (where $B_{1\,\mathrm{Mpc}}$ is the comoving field amplitude at a scale of 1 Mpc) at 95% confidence level, assuming zero helicity, and $B_{1\,\mathrm{Mpc}}< 5.6$ nG for a maximally helical field.For nearly scale-invariant PMFs we obtain $B_{1\,\mathrm{Mpc}}<2.0$ nG and $B_{1\,\mathrm{Mpc}}<0.9$ nG if the impact of PMFs on the ionization history of the Universe is included. From the analysis of magnetically-induced non-Gaussianity we obtain three different values, corresponding to three applied methods, all below 5 nG. The constraint from the magnetically-induced passive-tensor bispectrum is $B_{1\,\mathrm{Mpc}}< 2.8$ nG. A search for preferred directions in the magnetically-induced passive bispectrum yields $B_{1\,\mathrm{Mpc}}< 4.5$ nG, whereas the the compensated-scalar bispectrum gives $B_{1\,\mathrm{Mpc}}< 3$ nG. The analysis of the Faraday rotation of CMB polarization by PMFs uses the Planck power spectra in $EE$ and $BB$ at 70 GHz and gives $B_{1\,\mathrm{Mpc}}< 1380$ nG. In our final analysis, we consider the harmonic-space correlations produced by Alfv\'en waves, finding no significant evidence for the presence of these waves. Together, these results comprise a comprehensive set of constraints on possible PMFs with Planck data.

We present all-sky modelling of the high resolution Planck, IRAS, and WISE infrared (IR) observations using the physical dust model presented by Draine and Li in 2007 (DL). We study the performance and results of this model, and discuss implications for future dust modelling. The present work extends the DL dust modelling carried out on nearby galaxies using Herschel and Spitzer data to Galactic dust emission. We employ the DL dust model to generate maps of the dust mass surface density, the optical extinction Av, and the starlight intensity parametrized by Umin. The DL model reproduces the observed spectral energy distribution (SED) satisfactorily over most of the sky, with small deviations in the inner Galactic disk and in low ecliptic latitude areas. We compare the DL optical extinction Av for the diffuse interstellar medium with optical estimates for 2 10^5 quasi-stellar objects (QSOs) observed in the Sloan digital sky survey. The DL Av estimates are larger than those determined towards QSOs by a factor of about 2, which depends on Umin. The DL fitting parameter Umin, effectively determined by the wavelength where the SED peaks, appears to trace variations in the far-IR opacity of the dust grains per unit Av, and not only in the starlight intensity. To circumvent the model deficiency, we propose an empirical renormalization of the DL Av estimate, dependent of Umin, which compensates for the systematic differences found with QSO observations. This renormalization also brings into agreement the DL Av estimates with those derived for molecular clouds from the near-IR colours of stars in the 2 micron all sky survey. The DL model and the QSOs data are used to compress the spectral information in the Planck and IRAS observations for the diffuse ISM to a family of 20 SEDs normalized per Av, parameterized by Umin, which may be used to test and empirically calibrate dust models.

The all-sky Planck survey in 9 frequency bands was used to search for emission from all 274 known Galactic supernova remnants. Of these, 16 were detected in at least two Planck frequencies. The radio-through-microwave spectral energy distributions were compiled to determine the emission mechanism for microwave emission. In only one case, IC 443, is there high-frequency emission clearly from dust associated with the supernova remnant.In all cases, the low-frequency emission is from synchrotron radiation. A single power law, as predicted for a population of relativistic particles with energy distribution that extends continuously to high energies, is evident for many sources, including the Crab and PKS 1209-51/52. A decrease in flux density relative to the extrapolation of radio emission is evident in several sources. Their spectral energy distributions can be approximated as broken power laws, $S_\nu\propto\nu^{-\alpha}$, with the spectral index, alpha, increasing by 0.5-1 above a break frequency in the range 10-60 GHz. The break could be due to synchrotron losses.

We test the statistical isotropy and Gaussianity of the cosmic microwave background (CMB) anisotropies using observations made by the Planck satellite. Our results are based mainly on the full Planck mission for temperature, but also include some polarization measurements. In particular, we consider the CMB anisotropy maps derived from the multi-frequency Planck data by several component-separation methods. For the temperature anisotropies, we find excellent agreement between results based on these sky maps over both a very large fraction of the sky and a broad range of angular scales, establishing that potential foreground residuals do not affect our studies. Tests of skewness, kurtosis, multi-normality, N-point functions, and Minkowski functionals indicate consistency with Gaussianity, while a power deficit at large angular scales is manifested in several ways, for example low map variance. The results of a peak statistics analysis are consistent with the expectations of a Gaussian random field. The "Cold Spot" is detected with several methods, including map kurtosis, peak statistics, and mean temperature profile. We thoroughly probe the large-scale dipolar power asymmetry, detecting it with several independent tests, and address the subject of a posteriori correction. Tests of directionality suggest the presence of angular clustering from large to small scales, but at a significance that is dependent on the details of the approach. We perform the first examination of polarization data, finding the morphology of stacked peaks to be consistent with the expectations of statistically isotropic simulations. Where they overlap, these results are consistent with the Planck 2013 analysis based on the nominal mission data and provide our most thorough view of the statistics of the CMB fluctuations to date.

The Andromeda Galaxy (M31) is one of a few galaxies that has sufficient angular size on the sky to be resolved by the Planck satellite. Planck has detected M31 in all of its frequency bands, and has mapped out the dust emission with the High Frequency Instrument, clearly resolving multiple spiral arms and sub-features. We examine the morphology of this long-wavelength dust emission as seen by Planck, including a study of its outermost spiral arms, and investigate the dust heating mechanism across M31. We find that dust dominating the longer wavelength emission ($\gtrsim 0.3\,$mm) is heated by the diffuse stellar population (as traced by 3.6$\,\mu$m emission), with the dust dominating the shorter wavelength emission heated by a mix of the old stellar population and star-forming regions (as traced by 24$\,\mu$m emission). We also fit spectral energy distributions (SEDs) for individual 5' pixels and quantify the dust properties across the galaxy, taking into account these different heating mechanisms, finding that there is a linear decrease in temperature with galactocentric distance for dust heated by the old stellar population, as would be expected, with temperatures ranging from around 22$\,$K in the nucleus to 14$\,$K outside of the 10$\,$kpc ring. Finally, we measure the integrated spectrum of the whole galaxy, which we find to be well-fitted with a global dust temperature of ($18.2\pm1.0$)$\,$K with a spectral index of $1.62\pm0.11$ (assuming a single modified blackbody), and a significant amount of free-free emission at intermediate frequencies of 20-60$\,$GHz, which corresponds to a star formation rate of around $0.12$M$_\odot\,$yr$^{-1}$. We find a $2.3\,\sigma$ detection of the presence of spinning dust emission, with a 30$\,$GHz amplitude of $0.7\pm0.3\,$Jy, which is in line with expectations from our Galaxy.

Few independent detections of a weak X-ray line at an energy of $\sim$3.5 keV seen toward a number of astrophysical sites have been reported. If this signal will be confirmed to be the signature of decaying DM sterile neutrino with a mass of ~7.1 keV, then the cosmological observables should be consistent with its properties. In this paper we make a coupled treatment of the weak decoupling, primordial nucleosynthesis and photon decoupling epochs in the sterile neutrino resonant production scenario, including the extra radiation energy density via N_eff. We compute the radiation and matter perturbations including the full resonance sweep solution for active-sterile flavor conversion in the expanding Universe. We show that the cosmological measurements are in agreement with sub-dominant DM sterile neutrino resonant production with following parameters (errors at 95 CL): sterile neutrino mass 6.08 \pm 3.22 keV, mixing angle \sin^2 2 \theta < 5.61 x 10^{-10}, lepton number per flavor L_4 = 1.23 \pm 0.04 and sterile neutrino mass fraction f_\nu < 0.078. Our results are in good agreement with the sterile neutrino resonant production parameters inferred in Ref. [62] from the linear large scale structure constraints to produce full Dark Matter density.

The Planck design and scanning strategy provide many levels of redundancy that can be exploited to provide tests of internal consistency. One of the most important is the comparison of the 70GHz and 100GHz channels. Based on different instrument technologies, with feeds located differently in the focal plane, analysed independently by different teams using different software, and near the minimum of diffuse foreground emission, these channels are in effect two different experiments. The 143GHz channel has the lowest noise level on Planck, and is near the minimum of unresolved foreground emission. In this paper, we analyse the level of consistency achieved in the 2013 Planck data. We concentrate on comparisons between the 70/100/143GHz channel maps and power spectra, particularly over the angular scales of the first and second acoustic peaks, on maps masked for diffuse Galactic emission and for strong unresolved sources. Difference maps covering angular scales from 8deg-15arcmin are consistent with noise, and show no evidence of cosmic microwave background structure. Including small but important corrections for unresolved-source residuals, we demonstrate agreement between 70 and 100GHz power spectra averaged over 70<l<390 at the 0.8% level, and agreement between 143 and 100GHz power spectra of 0.4% over the same l range. These values are within and consistent with the overall uncertainties in calibration given in the Planck 2013 results. We also present results based on the 2013 likelihood analysis showing consistency at the 0.35% between the 100/143/217GHz power spectra. We analyse calibration procedures and beams to determine what fraction of these differences can be accounted for by known approximations or systematic errors that could be controlled even better in the future, reducing uncertainties still further. Several possible small improvements are described...(abridged)

The European Space Agency's Planck satellite, dedicated to studying the early Universe and its subsequent evolution, was launched 14~May 2009 and scanned the microwave and submillimetre sky continuously between 12~August 2009 and 23~October 2013. In February~2015, ESA and the Planck Collaboration released the second set of cosmology products based on data from the entire Planck mission, including both temperature and polarization, along with a set of scientific and technical papers and a web-based explanatory supplement. This paper gives an overview of the main characteristics of the data and the data products in the release, as well as the associated cosmological and astrophysical science results and papers. The science products include maps of the cosmic microwave background (CMB), the thermal Sunyaev-Zeldovich effect, and diffuse foregrounds in temperature and polarization, catalogues of compact Galactic and extragalactic sources (including separate catalogues of Sunyaev-Zeldovich clusters and Galactic cold clumps), and extensive simulations of signals and noise used in assessing the performance of the analysis methods and assessment of uncertainties. The likelihood code used to assess cosmological models against the Planck data are described, as well as a CMB lensing likelihood. Scientific results include cosmological parameters deriving from CMB power spectra, gravitational lensing, and cluster counts, as well as constraints on inflation, non-Gaussianity, primordial magnetic fields, dark energy, and modified gravity.

Our velocity relative to the rest frame of the cosmic microwave background (CMB) generates a dipole temperature anisotropy on the sky which has been well measured for more than 30 years, and has an accepted amplitude of v/c = 0.00123, or v = 369km/s. In addition to this signal generated by Doppler boosting of the CMB monopole, our motion also modulates and aberrates the CMB temperature fluctuations (as well as every other source of radiation at cosmological distances). This is an order 0.1% effect applied to fluctuations which are already one part in roughly one hundred thousand, so it is quite small. Nevertheless, it becomes detectable with the all-sky coverage, high angular resolution, and low noise levels of the Planck satellite. Here we report a first measurement of this velocity signature using the aberration and modulation effects on the CMB temperature anisotropies, finding a component in the known dipole direction, (l,b)=(264, 48) [deg], of 384km/s +- 78km/s (stat.) +- 115km/s (syst.). This is a significant confirmation of the expected velocity.

Planck data when combined with ancillary data provide a unique opportunity to separate the diffuse emission components of the inner Galaxy. The purpose of the paper is to elucidate the morphology of the various emission components in the strong star-formation region lying inside the solar radius and to clarify the relationship between the various components. The region of the Galactic plane covered is l=300-0-60deg where star-formation is highest and the emission is strong enough to make meaningful component separation. The latitude widths in this longitude range lie between 1deg and 2deg, which correspond to FWHM z-widths of 100-200pc at a typical distance of 6kpc. The four emission components studied here are synchrotron, free-free, anomalous microwave emission (AME), and thermal (vibrational) dust emission. These components are identified by constructing spectral energy distributions (SEDs) at positions along the Galactic plane using the wide frequency coverage of Planck (28.4-857GHz) in combination with low-frequency radio data at 0.408-2.3GHz plus WMAP data at 23-94GHz, along with far-infrared (FIR) data from DIRBE and IRAS. The free-free component is determined from radio recombination line (RRL) data. AME is found to be comparable in brightness to the free-free emission on the Galactic plane in the frequency range 20-40GHz with a width in latitude similar to that of the thermal dust; it comprises 45+/-1% of the total 28.4GHz emission in the longitude range l=300-0-60deg. The free-free component is the narrowest, reflecting the fact that it is produced by current star-formation as traced by the narrow distribution of OB stars. It is the dominant emission on the plane between 60 and 100GHz. RRLs from this ionized gas are used to assess its distance, leading to a free-free z-width of FWHM ~100pc...(abridged)

Planck has mapped the intensity and polarization of the sky at microwave frequencies with unprecedented sensitivity. We make use of the Planck 353 GHz I, Q, and U Stokes maps as dust templates, and cross-correlate them with the Planck and WMAP data at 12 frequencies from 23 to 353 GHz, over circular patches with 10 degree radius. The cross-correlation analysis is performed for both intensity and polarization data in a consistent manner. We use a mask that focuses our analysis on the diffuse interstellar medium at intermediate Galactic latitudes. We determine the spectral indices of dust emission in intensity and polarization between 100 and 353 GHz, for each sky-patch. The mean values, $1.59\pm0.02$ for polarization and $1.51\pm0.01$ for intensity, for a mean dust temperature of 19.6 K, are close, but significantly different ($3.6\,\sigma$). We determine the mean spectral energy distribution (SED) of the microwave emission, correlated with the 353 GHz dust templates, by averaging the results of the correlation over all sky-patches. We find that the mean SED increases for decreasing frequencies at $\nu < 60$ GHz, for both intensity and polarization. The rise of the polarization SED towards low frequencies may be accounted for by a synchrotron component correlated with dust, with no need for any polarization of the anomalous microwave emission. We use a spectral model to separate the synchrotron and dust polarization and to characterize the spectral dependence of the dust polarization fraction. The polarization fraction ($p$) of the dust emission decreases by $(21\pm6)$ % from 353 to 70 GHz. The decrease of $p$ could indicate differences in polarization efficiency among components of interstellar dust (e.g., carbon versus silicate grains). Our observational results provide inputs to quantify and optimize the separation between Galactic and cosmological polarization.

[Abridged] We use the Planck all-sky submm and mm maps to search for rare sources distinguished by extreme brightness, a few hundreds of mJy, and their potential for being situated at high redshift. These "cold" Planck sources, selected using the High Frequency Instrument (HFI) directly from the maps and from the Planck Catalogue of Compact Sources (PCCS), all satisfy the criterion of having their rest-frame far-infrared peak redshifted to the frequency range 353 and 857 GHz. This colour-selection favours galaxies in the redshift range z=2-4, which we consider as cold peaks in the cosmic infrared background (CIB). We perform a dedicated Herschel-SPIRE follow-up of 234 such Planck targets, finding a significant excess of red 350 and 500um sources, in comparison to reference SPIRE fields. About 94% of the SPIRE sources in the Planck fields are consistent with being overdensities of galaxies peaking at 350um. About 3% are candidate lensed systems, all 12 of which have secure spectroscopic confirmations, placing them at redshifts z>2.2. The galaxy overdensities are detected with high significance, half of the sample showing statistical significance above 10sigma. The SPIRE photometric redshifts of galaxies in overdensities suggest a peak at z~2. Under the Td=35K assumption, we derive an infrared (IR) luminosity for each SPIRE source of about 4x10^12 Lsun, yielding star formation rates of typically 700 Msun.yr^-1. If the observed overdensities are actual gravitationally-bound structures, the total total star formation rates reaches 7x10^3 Msun.yr^-1. Taken together, these sources show the signatures of high-z (z>$) protoclusters of intensively star-forming galaxies. All these observations confirm the uniqueness of our sample and demonstrate the ability of the all-sky Planck-HFI cold sources to select populations of cosmological and astrophysical interest for structure formation studies.

The role of the magnetic field in the formation of the filamentary structures observed in the interstellar medium (ISM) is a debated topic. The Planck all-sky maps of linearly polarized emission from dust at 353GHz provide the required combination of imaging and statistics to study the correlation between the structures of the Galactic magnetic field and of interstellar matter, both in the diffuse ISM and in molecular clouds. The data reveal structures, or ridges, in the intensity map with counterparts in the Stokes Q and/or U maps. We focus on structures at intermediate and high Galactic latitudes with column density from $10^{20}$ to $10^{22}$ cm$^{-2}$. We measure the magnetic field orientation on the plane of the sky from the polarization data, and present an algorithm to estimate the orientation of the ridges from the dust intensity map. We use analytical models to account for projection effects. Comparing polarization angles on and off the structures, we estimate the mean ratio between the strengths of the turbulent and mean components of the magnetic field to be between 0.6 and 1.0, with a preferred value of 0.8. We find that the ridges are preferentially aligned with the magnetic field measured on the structures. This trend becomes more striking for increasing polarization fraction and decreasing column density. We interpret the increase of alignment with polarization fraction as a consequence of projections effects. The decrease of alignment for high column density is not due to a loss of correlation between the structures and the geometry of the magnetic field. In molecular complexes, we observe structures perpendicular to the magnetic field, which cannot be accounted for by projection effects. We discuss our results in the context of models and MHD simulations, which describe the formation of structures in the magnetized ISM.

Shortened abstract: Observations of the nearby Chamaeleon clouds in gamma rays with the Fermi Large Area Telescope and in thermal dust emission with Planck and IRAS have been used with the HI and CO radio data to (i) map the gas column densities in the different phases and at the dark neutral medium (DNM) transition between the HI-bright and CO-bright media; (ii) constrain the CO-to-$H_2$ conversion factor, $X_{CO}$; (iii) probe the dust properties per gas nucleon in each gas phase and spatially across the clouds. We have separated clouds in velocity in HI and CO emission and modelled the 0.4-100 GeV intensity, the dust optical depth at 353 GHz, the thermal radiance of the large grains, and an estimate of the dust extinction empirically corrected for the starlight intensity, $A_{VQ}$. The gamma-ray emissivity spectra confirm that the GeV-TeV cosmic rays uniformly permeate all gas phases up to the CO cores. The dust and cosmic rays reveal large amounts of DNM gas, with comparable spatial distributions and twice as much mass as in the CO-bright clouds. We give constraints on the HI-DNM-CO transitions and CO-dark $H_2$ fractions for separate clouds. The corrected extinction provides the best fit to the total gas traced by the gamma rays, but we find evidence for a rise in $A_{VQ}/N_H$ and a steep rise in opacity, with increasing $N_H$ and $H_2$ fraction, and with decreasing dust temperature. We observe less variations for the specific power of the grains, except for a decline by half in the CO cores. This combined information suggests grain evolution. The gamma rays and dust radiance yield consistent $X_{CO}$ estimates near $0.7\times10^{20}$ cm$^{-2}$ (K km/s)$^{-1}$. The other dust tracers yield biased values because of the grain opacity rise in the CO clouds. These results also confirm a factor of 2 difference between $X_{CO}$ estimates at pc and kpc scales.

The Planck survey provides unprecedented full-sky coverage of the submillimetre polarized emission from Galactic dust, bringing new constraints on the properties of dust. The dust grains that emit the radiation seen by Planck in the submillimetre also extinguish and polarize starlight in the visible. Comparison of the polarization of the emission and of the interstellar polarization on selected lines of sight probed by stars provides unique new diagnostics of the emission and light scattering properties of dust. Using ancillary catalogues of interstellar polarization and extinction of starlight, we obtain the degree of polarization, $p_V$ , and the optical depth in the V band to the star, $\tau_V$. Toward these stars we measure the submillimetre polarized intensity, $P_S$, and total intensity, $I_S$, in the Planck 353 GHz channel. For those lines of sight through the diffuse interstellar medium with comparable values of the estimated column density and polarization directions close to orthogonal, we correlate properties in the submillimetre and visible to find two ratios, $R_{S/V} = (P_S/I_S)/(p_V/\tau_V)$ and $R_{P/p} = P_S/p_V$ , the latter focusing directly on the polarization properties of the aligned grain population alone. We find $R_{S/V}$ = 4.2, with statistical and systematic uncertainties 0.2 and 0.3, respectively, and $R_{P/p}$ = 5.4 MJy sr$^{-1}$, with uncertainties 0.2 and 0.3 MJy sr$^{-1}$, respectively. Our estimate of $R_{S/V}$ is compatible with predictions based on a range of polarizing dust models that have been developed for the diffuse interstellar medium. However, our estimate of $R_{P/p}$ is not compatible with predictions, which are too low by a factor of about 2.5. This more discriminating diagnostic, $R_{P/p}$, indicates that changes to the optical properties in the models of the aligned grain population are required.

We analyse the implications of the Planck data for cosmic inflation. The Planck nominal mission temperature anisotropy measurements, combined with the WMAP large-angle polarization, constrain the scalar spectral index to $n_s = 0.9603 \pm 0.0073$, ruling out exact scale invariance at over 5 $\sigma$. Planck establishes an upper bound on the tensor-to-scalar ratio of r < 0.11 (95% CL). The Planck data thus shrink the space of allowed standard inflationary models, preferring potentials with V" < 0. Exponential potential models, the simplest hybrid inflationary models, and monomial potential models of degree n > 2 do not provide a good fit to the data. Planck does not find statistically significant running of the scalar spectral index, obtaining $d n_s/d ln k = -0.0134 \pm 0.0090$. Several analyses dropping the slow-roll approximation are carried out, including detailed model comparison and inflationary potential reconstruction. We also investigate whether the primordial power spectrum contains any features. We find that models with a parameterized oscillatory feature improve the fit $\chi^2$ by ~ 10; however, Bayesian evidence does not prefer these models. We constrain several single-field inflation models with generalized Lagrangians by combining power spectrum data with bounds on $f_\mathrm{NL}$ measured by Planck. The fractional primordial contribution of CDM isocurvature modes in the curvaton and axion scenarios has upper bounds of 0.25% or 3.9% (95% CL), respectively. In models with arbitrarily correlated CDM or neutrino isocurvature modes, an anticorrelation can improve $\chi^2$ by approximatively 4 as a result of slightly lowering the theoretical prediction for the $\ell<40$ multipoles relative to the higher multipoles. Nonetheless, the data are consistent with adiabatic initial conditions.

We update the all-sky Planck catalogue of 1227 clusters and cluster candidates (PSZ1) published in March 2013, derived from Sunyaev-Zeldovich (SZ) effect detections using the first 15.5 months of Planck satellite observations. Addendum. We deliver an updated version of the PSZ1 catalogue, reporting the further confirmation of 86 Planck-discovered clusters. In total, the PSZ1 now contains 947 confirmed clusters, of which 214 were confirmed as newly discovered clusters through follow-up observations undertaken by the Planck Collaboration. The updated PSZ1 contains redshifts for 913 systems, of which 736 (~80.6%) are spectroscopic, and associated mass estimates derived from the Y_z mass proxy. We also provide a new SZ quality flag, derived from a novel artificial neural network classification of the SZ signal, for the remaining 280 candidates. Based on this assessment, the purity of the updated PSZ1 catalogue is estimated to be 94%. In this release, we provide the full updated catalogue and an additional readme file with further information on the Planck SZ detections.

The polarized thermal emission from Galactic dust is the main foreground present in measurements of the polarization of the cosmic microwave background (CMB) at frequencies above 100GHz. We exploit the Planck HFI polarization data from 100 to 353GHz to measure the dust angular power spectra $C_\ell^{EE,BB}$ over the range $40<\ell<600$ well away from the Galactic plane. These will bring new insights into interstellar dust physics and a precise determination of the level of contamination for CMB polarization experiments. We show that statistical properties of the emission can be characterized over large fractions of the sky using $C_\ell$. For the dust, they are well described by power laws in $\ell$ with exponents $\alpha^{EE,BB}=-2.42\pm0.02$. The amplitudes of the polarization $C_\ell$ vary with the average brightness in a way similar to the intensity ones. The dust polarization frequency dependence is consistent with modified blackbody emission with $\beta_d=1.59$ and $T_d=19.6$K. We find a systematic ratio between the amplitudes of the Galactic $B$- and $E$-modes of 0.5. We show that even in the faintest dust-emitting regions there are no "clean" windows where primordial CMB $B$-mode polarization could be measured without subtraction of dust emission. Finally, we investigate the level of dust polarization in the BICEP2 experiment field. Extrapolation of the Planck 353GHz data to 150GHz gives a dust power $\ell(\ell+1)C_\ell^{BB}/(2\pi)$ of $1.32\times10^{-2}\mu$K$_{CMB}^2$ over the $40<\ell<120$ range; the statistical uncertainty is $\pm0.29$ and there is an additional uncertainty (+0.28,-0.24) from the extrapolation, both in the same units. This is the same magnitude as reported by BICEP2 over this $\ell$ range, which highlights the need for assessment of the polarized dust signal even in the cleanest windows of the sky.

We have constructed the first all-sky map of the thermal Sunyaev-Zeldovich (tSZ) effect by applying specifically tailored component separation algorithms to the 100 to 857 GHz frequency channel maps from the Planck survey. These maps show an obvious galaxy cluster tSZ signal that is well matched with blindly detected clusters in the Planck SZ catalogue. To characterize the signal in the tSZ map we have computed its angular power spectrum. At large angular scales ($\ell < 60$), the major foreground contaminant is the diffuse thermal dust emission. At small angular scales ($\ell > 500$) the clustered Cosmic Infrared Background (CIB) and residual point sources are the major contaminants. These foregrounds are carefully modelled and subtracted. We measure the tSZ power spectrum in angular scales, $0.17^{\circ} \lesssim \theta \lesssim 3.0^{\circ}$, that were previously unexplored. The measured tSZ power spectrum is consistent with that expected from the Planck catalogue of SZ sources, with additional clear evidence of signal from unresolved clusters and, potentially, diffuse warm baryons. We use the tSZ power spectrum to obtain the following cosmological constraints: $\sigma_8(\Omega_{\mathrm{m}}/0.28)^{3.2/8.1}=0.784 \pm 0.016 (68% C.L.). Marginalized band-powers of the Planck tSZ power spectrum and the best-fit model are given. The non-Gaussianity of the Compton parameter map is further characterized by computing its 1D probability distribution function and its bispectrum. These are used to place additional independent constraints on $\sigma_{8}$.

Late stages of stellar evolution are characterized by copious mass-loss events whose signature is the formation of circumstellar envelopes (CSE). Planck multi-frequency measurements have provided relevant information on a sample of Galactic planetary nebulae (PNe) in the important and relatively unexplored observational band between 30 and 857GHz. Planck enables the assembly of comprehensive PNe spectral energy distributions (SEDs) from radio {\bf to} far-infrared frequencies. Modelling of the derived SEDs provides us with information on physical properties of CSEs and the mass content of both main components: ionised gas, traced by the free-free emission at cm--mm waves; and thermal dust, traced by the millimetre and far-IR emission. In particular, the amount of ionised gas and dust has been derived here. Such quantities have also been estimated for the very young PN CRL618, where the strong variability observed in its radio and millimetre emission has previously prevented the construction of its SED. A morphological study of the Helix Nebula has also been performed. Planck maps reveal, for the first time, the spatial distribution of the dust inside the envelope, allowing us to identify different components, the most interesting of which is a very extended component (up to 1pc) that may be related to a region where the slow expanding envelope is interacting with the surrounding interstellar medium.

Polarized emission observed by Planck HFI at 353 GHz towards a sample of nearby fields is presented, focusing on the statistics of polarization fractions $p$ and angles $\psi$. The polarization fractions and column densities in these nearby fields are representative of the range of values obtained over the whole sky. We find that: (i) the largest polarization fractions are reached in the most diffuse fields; (ii) the maximum polarization fraction $p_\mathrm{max}$ decreases with column density $N_\mathrm{H}$ in the more opaque fields with $N_\mathrm{H} > 10^{21}\,\mathrm{cm}^{-2}$; and (iii) the polarization fraction along a given line of sight is correlated with the local spatial coherence of the polarization angle. These observations are compared to polarized emission maps computed in simulations of anisotropic magnetohydrodynamical (MHD) turbulence in which we assume a uniform intrinsic polarization fraction of the dust grains. We find that an estimate of this parameter may be recovered from the maximum polarization fraction $p_\mathrm{max}$ in diffuse regions where the magnetic field is ordered on large scales and perpendicular to the line of sight. This emphasizes the impact of anisotropies of the magnetic field on the emerging polarization signal. The decrease of the polarization fraction with column density in nearby molecular clouds is well reproduced in the simulations, indicating that it is essentially due to the turbulent structure of the magnetic field: an accumulation of variously polarized structures along the line of sight leads to such an anti-correlation. In the simulations, polarization fractions are also found to anti-correlate with the angle dispersion function $\mathcal{S}$. [abridged]

This paper presents an all-sky model of dust emission from the Planck 857, 545 and 353 GHz, and IRAS 100 micron data. Using a modified black-body fit to the data we present all-sky maps of the dust optical depth, temperature, and spectral index over the 353-3000 GHz range. This model is a tight representation of the data at 5 arc min. It shows variations of the order of 30 % compared with the widely-used model of Finkbeiner, Davis, and Schlegel. The Planck data allow us to estimate the dust temperature uniformly over the whole sky, providing an improved estimate of the dust optical depth compared to previous all-sky dust model, especially in high-contrast molecular regions. An increase of the dust opacity at 353 GHz, tau_353/N_H, from the diffuse to the denser interstellar medium (ISM) is reported. It is associated with a decrease in the observed dust temperature, T_obs, that could be due at least in part to the increased dust opacity. We also report an excess of dust emission at HI column densities lower than 10^20 cm^-2 that could be the signature of dust in the warm ionized medium. In the diffuse ISM at high Galactic latitude, we report an anti-correlation between tau_353/N_H and T_obs while the dust specific luminosity, i.e., the total dust emission integrated over frequency (the radiance) per hydrogen atom, stays about constant. The implication is that in the diffuse high-latitude ISM tau_353 is not as reliable a tracer of dust column density as we conclude it is in molecular clouds where the correlation of tau_353 with dust extinction estimated using colour excess measurements on stars is strong. To estimate Galactic E(B-V) in extragalactic fields at high latitude we develop a new method based on the thermal dust radiance, instead of the dust optical depth, calibrated to E(B-V) using reddening measurements of quasars deduced from Sloan Digital Sky Survey data.

We discuss the methods employed to photometrically calibrate the data acquired by the Low Frequency Instrument on Planck. Our calibration is based on a combination of the Orbital Dipole plus the Solar Dipole, caused respectively by the motion of the Planck spacecraft with respect to the Sun and by motion of the Solar System with respect to the CMB rest frame. The latter provides a signal of a few mK with the same spectrum as the CMB anisotropies and is visible throughout the mission. In this data release we rely on the characterization of the Solar Dipole as measured by WMAP. We also present preliminary results (at 44GHz only) on the study of the Orbital Dipole, which agree with the WMAP value of the Solar System speed within our uncertainties. We compute the calibration constant for each radiometer roughly once per hour, in order to keep track of changes in the detectors' gain. Since non-idealities in the optical response of the beams proved to be important, we implemented a fast convolution algorithm which considers the full beam response in estimating the signal generated by the dipole. Moreover, in order to further reduce the impact of residual systematics due to sidelobes, we estimated time variations in the calibration constant of the 30GHz radiometers (the ones with the largest sidelobes) using the signal of a reference load. We have estimated the calibration accuracy in two ways: we have run a set of simulations to assess the impact of statistical errors and systematic effects in the instrument and in the calibration procedure, and we have performed a number of consistency checks on the data and on the brightness temperature of Jupiter. Calibration errors for this data release are expected to be about 0.6% at 44 and 70 GHz, and 0.8% at 30 GHz. (Abriged.)

On the arcminute angular scales probed by Planck, the CMB anisotropies are gently perturbed by gravitational lensing. Here we present a detailed study of this effect, detecting lensing independently in the 100, 143, and 217GHz frequency bands with an overall significance of greater than 25sigma. We use the temperature-gradient correlations induced by lensing to reconstruct a (noisy) map of the CMB lensing potential, which provides an integrated measure of the mass distribution back to the CMB last-scattering surface. Our lensing potential map is significantly correlated with other tracers of mass, a fact which we demonstrate using several representative tracers of large-scale structure. We estimate the power spectrum of the lensing potential, finding generally good agreement with expectations from the best-fitting LCDM model for the Planck temperature power spectrum, showing that this measurement at z=1100 correctly predicts the properties of the lower-redshift, later-time structures which source the lensing potential. When combined with the temperature power spectrum, our measurement provides degeneracy-breaking power for parameter constraints; it improves CMB-alone constraints on curvature by a factor of two and also partly breaks the degeneracy between the amplitude of the primordial perturbation power spectrum and the optical depth to reionization, allowing a measurement of the optical depth to reionization which is independent of large-scale polarization data. Discarding scale information, our measurement corresponds to a 4% constraint on the amplitude of the lensing potential power spectrum, or a 2% constraint on the RMS amplitude of matter fluctuations at z~2.

We present the results of approximately three years of observations of Planck Sunyaev-Zeldovich (SZ) sources with the Russian-Turkish 1.5-m telescope (RTT150), as a part of the optical follow-up programme undertaken by the Planck collaboration. During this time period approximately 20% of all dark and grey clear time available at the telescope was devoted to observations of Planck objects. Some observations of distant clusters were also done at the 6-m Bolshoy Telescope Azimutal'ny (BTA) of the Special Astrophysical Observatory of the Russian Academy of Sciences. In total, deep, direct images of more than one hundred fields were obtained in multiple filters. We identified 47 previously unknown galaxy clusters, 41 of which are included in the Planck catalogue of SZ sources. The redshifts of 65 Planck clusters were measured spectroscopically and 14 more were measured photometrically. We discuss the details of cluster optical identifications and redshift measurements. We also present new spectroscopic redhifts for 39 Planck clusters that were not included in the Planck SZ source catalogue and are published here for the first time.