Search

We have used the Planck all-sky submillimetre and millimetre maps to search for rare sources distinguished by extreme brightness, a few hundred millijanskies, 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-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. With a 4.′5 beam at the four highest frequencies, our sample is expected to include overdensities of galaxies in groups or clusters, lensed galaxies, and chance line-of-sight projections. We perform a dedicated Herschel-SPIRE follow-up of 234 such Planck targets, finding a significant excess of red 350 and 500 μm 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 350 μm, with 3% peaking at 500 μm, and none peaking at 250 μm. About 3% are candidate lensed systems, all 12 of which have secure spectroscopic confirmations, placing them at redshifts z > 2.2. Only four targets are Galactic cirrus, yielding a success rate in our search strategy for identifying extragalactic sources within the Planck beam of better than 98%. The galaxy overdensities are detected with high significance, half of the sample showing statistical significance above 10σ. The SPIRE photometric redshifts of galaxies in overdensities suggest a peak at z ≃ 2, assuming a single common dust temperature for the sources of Td = 35 K. Under this assumption, we derive an infrared (IR) luminosity for each SPIRE source of about 4 × 1012 L⊙, yielding star formation rates of typically 700 M⊙ yr-1. If the observed overdensities are actual gravitationally-bound structures, the total IR luminosity of all their SPIRE-detected sources peaks at 4 × 1013 L⊙, leading to total star formation rates of perhaps 7 × 103 M⊙yr-1 per overdensity. Taken together, these sources show the signatures of high-z (z > 2) protoclusters of intensively star-forming galaxies. All these observations confirm the uniqueness of our sample compared to reference samples and demonstrate the ability of the all-sky Planck-HFI cold sources to select populations of cosmological and astrophysical interest for structure formation studies.

We update the all-sky Planck catalogue of 1227 clusters and cluster candidates (PSZ1) published in March 2013, derived from detections of the Sunyaev-Zeldovich (SZ) effect using the first 15.5 months of Planck satellite observations. As an 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 Yz mass proxy. We also provide a new SZ quality flag for the remaining 280 candidates. This flag was derived from a novel artificial neural-network classification of the SZ signal. 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.

We update the all-sky Planck catalogue of 1227 clusters and cluster candidates (PSZ1) published in March 2013, derived from detections of the Sunyaev-Zeldovich (SZ) effect using the first 15.5 months of Planck satellite observations. As an 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 Yz mass proxy. We also provide a new SZ quality flag for the remaining 280 candidates. This flag was derived from a novel artificial neural-network classification of the SZ signal. 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.

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° → 60° wherestar-formation is highest and the emission is strong enough to make meaningful component separation. The latitude widths in this longitude range lie between 1 and 2, which correspond to FWHM z-widths of 100-200 pc at a typical distance of 6 kpc. 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.3 GHz plus WMAP data at 23-94 GHz, along with far-infrared (FIR) data from COBE-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-40 GHz with a width in latitude similar to that of the thermal dust; it comprises 45 ± 1% of the total 28.4 GHz emission in the longitude range l = 300° → 0° → 60°. 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 100 GHz. RRLs from this ionized gas are used to assess its distance, leading to a free-free z-width of FWHM ≈ 100 pc. The narrow synchrotron component has a low-frequency brightness spectral index βsynch ≈ -2.7 that is similar to the broad synchrotron component indicating that they are both populated by the cosmic ray electrons of the same spectral index. The width of this narrow synchrotron component is significantly larger than that of the other three components, suggesting that it is generated in an assembly of older supernova remnants that have expanded to sizes of order 150 pc in 3 × 105 yr; pulsars of a similar age have a similar spread in latitude. The thermal dust is identified in the SEDs with average parameters of Tdust = 20.4 ± 0.4 K, βFIR = 1.94 ± 0.03 (> 353 GHz), and βmm = 1.67 ± 0.02 (< 353 GHz). The latitude distributions of gamma-rays, CO, and the emission in high-frequency Planck bands have similar widths, showing that they are all indicators of the total gaseous matter on the plane in the inner Galaxy.

Any variation in the fundamental physical constants, more particularly in the fine structure constant, α, or in the mass of the electron, me, affects the recombination history of the Universe and cause an imprint on the cosmic microwave background angular power spectra. We show that the Planck data allow one to improve the constraint on the time variation of the fine structure constant at redshift z ∼ 103 by about a factor of 5 compared to WMAP data, as well as to break the degeneracy with the Hubble constant, H0. In addition to α, we can set a constraint on the variation in the mass of the electron, me, and in the simultaneous variation of the two constants. We examine in detail the degeneracies between fundamental constants and the cosmological parameters, in order to compare the limits obtained from Planck and WMAP and to determine the constraining power gained by including other cosmological probes. We conclude that independent time variations of the fine structure constant and of the mass of the electron are constrained by Planck to Δα/α = (3.6 ± 3.7) × 10-3 and Δme/me = (4 ± 11) × 10-3 at the 68% confidence level. We also investigate the possibility of a spatial variation of the fine structure constant. The relative amplitude of a dipolar spatial variation in α (corresponding to a gradient across our Hubble volume) is constrained to be δα/α = (-2.4 ± 3.7) × 10-2.

This paper presents an overview of the polarized sky as seen by Planck HFI at 353 GHz, which is the most sensitive Planck channel for dust polarization. We construct and analyse maps of dust polarization fraction and polarization angle at 1° resolution, taking into account noise bias and possible systematic effects. The sensitivity of the Planck HFI polarization measurements allows for the first time a mapping of Galactic dust polarized emission on large scales, including low column density regions. We find that the maximum observed dust polarization fraction is high (pmax = 19.8%), in particular in some regions of moderate hydrogen column density (NH < 2 × 1021 cm-2). The polarization fraction displays a large scatter at NH below a few 1021 cm-2. There is a general decrease in the dust polarization fraction with increasing column density above NH ∼ 1 × 1021 cm-2 and in particular a sharp drop above NH ∼ 1.5 × 1022 cm-2. We characterize the spatial structure of the polarization angle using the angle dispersion function. We find that the polarization angle is ordered over extended areas of several square degrees, separated by filamentary structures of high angle dispersion function. These appear as interfaces where the sky projection of the magnetic field changes abruptly without variations in the column density. The polarization fraction is found to be anti-correlated with the dispersion of polarization angles. These results suggest that, at the resolution of 1°, depolarization is due mainly to fluctuations in the magnetic field orientation along the line of sight, rather than to the loss of grain alignment in shielded regions. We also compare the polarization of thermal dust emission with that of synchrotron measured with Planck, low-frequency radio data, and Faraday rotation measurements toward extragalactic sources. These components bear resemblance along the Galactic plane and in some regions such as the Fan and North Polar Spur regions. The poor match observed in other regions shows, however, that dust, cosmic-ray electrons, and thermal electrons generally sample different parts of the line of sight.

The Planck survey provides unprecedented full-sky coverage of the submillimetre polarized emission from Galactic dust. In addition to the information on the direction of the Galactic magnetic field, this also brings 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, and therefore of the important dust model parameters, composition, size, and shape. Using ancillary catalogues of interstellar polarization and extinction of starlight, we obtain the degree of polarization, pV, and the optical depth in the V band to the star, τV. Toward these stars we measure the submillimetre polarized intensity, PS, and total intensity, IS, in the Planck 353 GHz channel. We compare the column density measure in the visible, E(B - V), with that inferred from the Planck product map of the submillimetre dust optical depth and compare the polarization direction (position angle) in the visible with that in the submillimetre. 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, RS/V = (PS/IS)/(pV/τV) and RP/p = PS/pV, the latter focusing directly on the polarization properties of the aligned grain population alone. We find RS/V = 4.2, with statistical and systematic uncertainties 0.2 and 0.3, respectively, and RP/p = 5.4 MJy sr-1, with uncertainties 0.2 and 0.3 MJy sr-1, respectively. Our estimate of RS/V is compatible with predictions based on a range of polarizing dust models that have been developed for the diffuse interstellar medium. This estimate provides new empirical validation of many of the common underlying assumptions of the models, but is not yet very discriminating among them. However, our estimate of RP/p is not compatible with predictions, which are too low by a factor of about 2.5. This more discriminating diagnostic, RP/p, indicates that changes to the optical properties in the models of the aligned grain population are required. These new diagnostics, together with the spectral dependence in the submillimetre from Planck, will be important for constraining and understanding the full complexity of the grain models, and for interpreting the Planck thermal dust polarization and refinement of the separation of this contamination of the cosmic microwave background.

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 ψ. 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 pmax decreases with column density NH in the more opaque fields with NH> 1021 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 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 pmax 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 maximum 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 S. However, the dispersion of the polarization angle for a given polarization fraction is found to be larger in the simulations than in the observations, suggesting a shortcoming in the physical content of these numerical models. In summary, we find that the turbulent structure of the magnetic field is able to reproduce the main statistical properties of the dust polarization as observed in a variety of nearby clouds, dense cores excluded, and that the large-scale field orientation with respect to the line of sight plays a major role in the quantitative analysis of these statistical properties.

Planck has mapped the intensity and polarization of the sky at microwave frequencies with unprecedented sensitivity. We use these data to characterize the frequency dependence of dust emission. 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° radius. The cross-correlation analysis is performed for both intensity and polarization data in a consistent manner. The results are corrected for the chance correlation between the templates and the anisotropies of the cosmic microwave background. 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. Both indices are found to be remarkably constant over the sky. The mean values, 1.59 ± 0.02 for polarization and 1.51 ± 0.01 for intensity, for a mean dust temperature of 19.6 K, are close, but significantly different (3.6σ). 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 ν< 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 ± 6)% from 353 to 70 GHz. We discuss this result within the context of existing dust models. The decrease in 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.

Planck has mapped the intensity and polarization of the sky at microwave frequencies with unprecedented sensitivity. We use these data to characterize the frequency dependence of dust emission. 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° radius. The cross-correlation analysis is performed for both intensity and polarization data in a consistent manner. The results are corrected for the chance correlation between the templates and the anisotropies of the cosmic microwave background. 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. Both indices are found to be remarkably constant over the sky. The mean values, 1.59 ± 0.02 for polarization and 1.51 ± 0.01 for intensity, for a mean dust temperature of 19.6 K, are close, but significantly different (3.6σ). 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 ν< 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 ± 6)% from 353 to 70 GHz. We discuss this result within the context of existing dust models. The decrease in 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.

We report the results of a joint analysis of data from BICEP2/Keck Array and Planck. BICEP2 and Keck Array have observed the same approximately 400  deg^{2} patch of sky centered on RA 0 h, Dec. -57.5°. The combined maps reach a depth of 57 nK deg in Stokes Q and U in a band centered at 150 GHz. Planck has observed the full sky in polarization at seven frequencies from 30 to 353 GHz, but much less deeply in any given region (1.2  μK deg in Q and U at 143 GHz). We detect 150×353 cross-correlation in B modes at high significance. We fit the single- and cross-frequency power spectra at frequencies ≥150  GHz to a lensed-ΛCDM model that includes dust and a possible contribution from inflationary gravitational waves (as parametrized by the tensor-to-scalar ratio r), using a prior on the frequency spectral behavior of polarized dust emission from previous Planck analysis of other regions of the sky. We find strong evidence for dust and no statistically significant evidence for tensor modes. We probe various model variations and extensions, including adding a synchrotron component in combination with lower frequency data, and find that these make little difference to the r constraint. Finally, we present an alternative analysis which is similar to a map-based cleaning of the dust contribution, and show that this gives similar constraints. The final result is expressed as a likelihood curve for r, and yields an upper limit r_{0.05}

The multi-frequency capability of the Planck satellite provides information both on the integrated history of star formation (via the cosmic infrared background, or CIB) and on the distribution of dark matter (via the lensing effect on the cosmic microwave background, or CMB). The conjunction of these two unique probes allows us to measure directly the connection between dark and luminous matter in the high redshift (1 < z 1. We measure directly the SFR density with around 2 sigma significance for three redshift bins between z=1 and 7, thus opening a new window into the study of the formation of stars at early times.

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 70 GHz (amplifier) and 100 GHz (bolometer) 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 143 GHz 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, and 143 GHz 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 8 to 15' 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 (measured by deviation of the ratio from unity) between 70 and 100 GHz power spectra averaged over 70 ≤ l ≤390 at the 0.8% level, and agreement between 143 and 100 GHz 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, and 217 GHz power spectra. We analyse calibration procedures and beams to determine what fraction of these differences can be accounted for by known approximations or systematicerrors that could be controlled even better in the future, reducing uncertainties still further. Several possible small improvements are described. Subsequent analysis of the beams quantifies the importance of asymmetry in the near sidelobes, which was not fully accounted for initially, affecting the 70/100 ratio. Correcting for this, the 70, 100, and 143 GHz power spectra agree to 0.4% over the first two acoustic peaks. The likelihood analysis that produced the 2013 cosmological parameters incorporated uncertainties larger than this. We show explicitly that correction of the missing near sidelobe power in the HFI channels would result in shifts in the posterior distributions of parameters of less than 0.3σ except for As, the amplitude of the primordial curvature perturbations at 0.05 Mpc-1, which changes by about 1σ. We extend these comparisons to include the sky maps from the complete nine-year mission of the Wilkinson Microwave Anisotropy Probe (WMAP), and find a roughly 2% difference between the Planck and WMAP power spectra in the region of the first acoustic peak.

This paper presents an all-sky model of dust emission from the Planck 353, 545, and 857 GHz, and IRAS 100 μm data. Using a modified blackbody 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 good representation of the IRAS and Planck data at 5′ between 353 and 3000 GHz (850 and 100 μm). 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, down to an angular resolution of 5′, providing an improved estimate of the dust optical depth compared to previous all-sky dust model, especially in high-contrast molecular regions where the dust temperature varies strongly at small scales in response to dust evolution, extinction, and/or local production of heating photons. An increase of the dust opacity at 353 GHz, τ353/NH, from the diffuse to the denser interstellar medium (ISM) is reported. It is associated with a decrease in the observed dust temperature, Tobs, that could be due at least in part to the increased dust opacity. We also report an excess of dust emission at H i column densities lower than 1020 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 anticorrelation between τ353/NH and Tobs while the dust specific luminosity, i.e., the total dust emission integrated over frequency (the radiance) per hydrogen atom, stays about constant, confirming one of the Planck Early Results obtained on selected fields. This effect is compatible with the view that, in the diffuse ISM, Tobs responds to spatial variations of the dust opacity, due to variations of dust properties, in addition to (small) variations of the radiation field strength. The implication is that in the diffuse high-latitude ISM τ353 is not as reliable a tracer of dust column density as we conclude it is in molecular clouds where the correlation of τ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.

© 2014 ESO. We describe the processing of the 531 billion raw data samples from the High Frequency Instrument (HFI), which we performed to produce six temperature maps from the first 473 days of Planck-HFI survey data. These maps provide an accurate rendition of the sky emission at 100, 143, 217, 353, 545, and 857 GHz with an angular resolution ranging from 9.I 7 to 4.I 6. The detector noise per (effective) beam solid angle is respectively, 10, 6, 12, and 39 μK in the four lowest HFI frequency channels (100-353 GHz) and 13 and 14 kJy sr-1 in the 545 and 857 GHz channels. Relative to the 143 GHz channel, these two high frequency channels are calibrated to within 5% and the 353 GHz channel to the percent level. The 100 and 217 GHz channels, which together with the 143 GHz channel determine the high-multipole part of the CMB power spectrum (50., 2500), are calibrated relative to 143 GHz to better than 0.2%.

© 2014 ESO. Planck has produced detailed all-sky observations over nine frequency bands between 30 and 857 GHz. These observations allow robust reconstruction of the primordial cosmic microwave background (CMB) temperature fluctuations over nearly the full sky, as well as new constraints on Galactic foregrounds, including thermal dust and line emission from molecular carbon monoxide (CO). This paper describes the component separation framework adopted by Planck for many cosmological analyses, including CMB power spectrum determination and likelihood construction on large angular scales, studies of primordial non-Gaussianity and statistical isotropy, the integrated Sachs-Wolfe effect, gravitational lensing, and searches for topological defects. We test four foreground-cleaned CMB maps derived using qualitatively different component separation algorithms. The quality of our reconstructions is evaluated through detailed simulations and internal comparisons, and shown through various tests to be internally consistent and robust for CMB power spectrum and cosmological parameter estimation up to â.," = 2000. The parameter constraints on ΛCDM cosmologies derived from these maps are consistent with those presented in the cross-spectrum based Planck likelihood analysis. We choose two of the CMB maps for specific scientific goals. We also present maps and frequency spectra of the Galactic low-frequency, CO, and thermal dust emission. The component maps are found to provide a faithful representation of the sky, as evaluated by simulations, with the largest bias seen in the CO component at 3%. For the low-frequency component, the spectral index varies widely over the sky, ranging from about β =-4 to-2. Considering both morphology and prior knowledge of the low frequencycomponents, the index map allows us to associate a steep spectral index (β-2.3 with strong free-free emission, and intermediate values with synchrotron emission.

Rotational transition lines of CO play a major role in molecular radio astronomy as a mass tracer and in particular in the study of star formation and Galactic structure. Although a wealth of data exists for the Galactic plane and some well-known molecular clouds, there is no available high sensitivity all-sky survey of CO emission to date. Such all-sky surveys can be constructed using the Planck HFI data because the three lowest CO rotational transition lines at 115, 230 and 345 GHz significantly contribute to the signal of the 100, 217 and 353 GHz HFI channels, respectively. Two different component separation methods are used to extract the CO maps from Planck HFI data. The maps obtained are then compared to one another and to existing external CO surveys. From these quality checks the best CO maps, in terms of signal to noise ratio and/or residual contamination by other emission, are selected. Three different sets of velocity-integrated CO emission maps are produced with different trade-offs between signal-to-noise, angular resolution, and reliability. Maps for the CO J = 1 → 0, J = 2 → 1, and J = 3 → 2 rotational transitions are presented and described in detail. They are shown to be fully compatible with previous surveys of parts of the Galactic plane as well as with undersampled surveys of the high latitude sky. The Planck HFI velocity-integrated CO maps for the J = 1 → 0, J = 2 → 1, and J = 3 →2 rotational transitions provide an unprecedented all-sky CO view of the Galaxy. These maps are also of great interest to monitor potential CO contamination of the Planck studies of the cosmological microwave background.

This paper presents the Planck 2013 likelihood, a complete statistical description of the two-point correlation function of the CMB temperature fluctuations that accounts for all known relevant uncertainties, both instrumental and astrophysical in nature. We use this likelihood to derive our best estimate of the CMB angular power spectrum from Planck over three decades in multipole moment, covering 22500. The main source of uncertainty at 1500 is cosmic variance. Uncertainties in small-scale foreground modelling and instrumental noise dominate the error budget at higher s. For < 50, our likelihood exploits all Planck frequency channels from 30 to 353 GHz, separating the cosmological CMB signal from diffuse Galactic foregrounds through a physically motivated Bayesian component separation technique. At 50, we employ a correlated Gaussian likelihood approximation based on a fine-grained set of angular cross-spectra derived from multiple detector combinations between the 100, 143, and 217 GHz frequency channels, marginalising over power spectrum foreground templates. We validate our likelihood through an extensive suite of consistency tests, and assess the impact of residual foreground and instrumental uncertainties on the final cosmological parameters. We find good internal agreement among the high-cross-spectra with residuals below a few K2at 1000, in agreement with estimated calibration uncertainties. We compare our results with foreground-cleaned CMB maps derived from all Planck frequencies, as well as with cross-spectra derived from the 70 GHz Planck map, and find broad agreement in terms of spectrum residuals and cosmological parameters. We further show that the best-fit CDM cosmology is in excellent agreement with preliminary PlanckEE and TE polarisation spectra. We find that the standard CDM cosmology is well constrained by Planck from the measurements at 1500. One specific example is the spectral index of scalar perturbations, for which we report a 5.4 deviation from scale invariance, n= 1. Increasingthe multipole range beyond 1500 does not increase our accuracy for the CDM parameters, but instead allows us to study extensions beyond the standard model. We find no indication of significant departures from the CDM framework. Finally, we report a tension between the Planck best-fit CDM model and the low-spectrum in the form of a power deficit of 510% at 40, with a statistical significance of 2.53. Without a theoretically motivated model for this power deficit, we do not elaborate further on its cosmological implications, but note that this is our most puzzling finding in an otherwise remarkably consistent data set.

© 2014 ESO. Based on cosmic microwave background (CMB) maps from the 2013 Planck Mission data release, this paper presents the detection of the integrated Sachs-Wolfe (ISW) effect, that is, the correlation between the CMB and large-scale evolving gravitational potentials. The significance of detection ranges from 2 to 4σ, depending on which method is used. We investigated three separate approaches, which essentially cover all previous studies, and also break new ground. (i) We correlated the CMB with the Planck reconstructed gravitational lensing potential (for the first time). This detection was made using the lensing-induced bispectrum between the low-â.," and high-â.," temperature anisotropies; the correlation between lensing and the ISW effect has a significance close to 2.5σ. (ii) We cross-correlated with tracers of large-scale structure, which yielded a significance of about 3σ, based on a combination of radio (NVSS) and optical (SDSS) data. (iii) We used aperture photometry on stacked CMB fields at the locations of known large-scale structures, which yielded and confirms a 4σ signal, over a broader spectral range, when using a previously explored catalogue, but shows strong discrepancies in amplitude and scale when compared with expectations. More recent catalogues give more moderate results that range from negligible to 2.5σ at most, but have a more consistent scale and amplitude, the latter being still slightly higher than what is expected from numerical simulations within ΛCMD. Where they can be compared, these measurements are compatible with previous work using data from WMAP, where these scales have been mapped to the limits of cosmic variance. Planck's broader frequency coverage allows for better foreground cleaning and confirms that the signal is achromatic, which makes it preferable for ISW detection. As a final step we used tracers of large-scale structure to filter the CMB data, from which we present maps of the ISW temperature perturbation. These results provide complementary and independent evidence for the existence of a dark energy component that governs the currently accelerated expansion of the Universe.

We describe the detection, interpretation, and removal of the signal resulting from interactions of high energy particles with the Planck High Frequency Instrument (HFI). There are two types of interactions: heating of the 0.1 K bolometer plate; and glitches in each detector time stream. The transientresponses to detector glitch shapes are not simple single-pole exponential decays and fall into three families. The glitch shape for each family has been characterized empirically in flight data and these shapes have been used to remove glitches from the detector time streams. The spectrum of the count rate per unit energy is computed for each family and a correspondence is made to the location on the detector of the particle hit. Most of the detected glitches are from Galactic protons incident on the die frame supporting the micro-machined bolometric detectors. In the Planck orbit at L2, the particle flux is around 5 cm-2 s-1 and is dominated by protons incident on the spacecraft with energy >39 MeV, at a rate of typically one event per second per detector. Different categories of glitches have different signatures in the time stream. Two of the glitch types have a low amplitude component that decays over nearly 1 s. This component produces excess noise if not properly removed from the time-ordered data. We have used a glitch detection and subtraction method based on the joint fit of population templates. The application of this novel glitch subtraction method removes excess noise from the time streams. Using realistic simulations, we find that this method does not introduce signal bias into the Planck data.

The new cosmic microwave background (CMB) temperature maps from Planck provide the highest-quality full-sky view of the surface of last scattering available to date. This allows us to detect possible departures from the standard model of a globally homogeneous and isotropic cosmology on the largest scales. We search for correlations induced by a possible non-trivial topology with a fundamental domain intersecting, or nearly intersecting, the last scattering surface (at comoving distance χrec), both via a direct search for matched circular patterns at the intersections and by an optimal likelihood search for specific topologies. For the latter we consider flat spaces with cubic toroidal (T3), equal-sided chimney (T2) and slab (T1) topologies, three multi-connected spaces of constant positive curvature (dodecahedral, truncated cube and octahedral) and two compact negative-curvature spaces. These searches yield no detection of the compact topology with the scale below the diameter of the last scattering surface. For most compact topologies studied the likelihood maximized over the orientation of the space relative to the observed map shows some preference for multi-connected models just larger than the diameter of the last scattering surface. Since this effect is also present in simulated realizations of isotropic maps, we interpret it as the inevitable alignment of mild anisotropic correlations with chance features in a single sky realization; such a feature can also be present, in milder form, when the likelihood is marginalized over orientations. Thus marginalized, the limits on the radius i of the largest sphere inscribed in topological domain (at log-likelihood-ratio Δln ' >-5 relative to a simply-connected flat Planck best-fit model) are: in a flat Universe, i> 0.92χrec for the T3 cubic torus; i> 0.71χrec for the T2 chimney; i> 0.50χrec for the T1 slab; and in a positively curved Universe, i> 1.03χrec for the dodecahedral space; i> 1.0χrec for the truncated cube; and i> 0.89χrec for the octahedral space. The limit for a wider class of topologies, i.e., those predicting matching pairs of back-to-back circles, among them tori and the three spherical cases listed above, coming from the matched-circles search, is i> 0.94χrec at 99% confidence level. Similar limits apply to a wide, although not exhaustive, range of topologies. We also perform a Bayesian search for an anisotropic global Bianchi VIIh geometry. In the non-physical setting where the Bianchi cosmology is decoupled from the standard cosmology, Planck data favour the inclusion of a Bianchi component with a Bayes factor of at least 1.5 units of log-evidence. Indeed, the Bianchi pattern is quite efficient at accounting for some of the large-scale anomalies found in Planck data. However, the cosmological parameters that generate this pattern are in strong disagreement with those found from CMB anisotropy data alone. In the physically motivated setting where the Bianchi parameters are coupled and fitted simultaneously with the standard cosmological parameters, we find no evidence for a Bianchi VIIh cosmology and constrain the vorticity of such models to (ω/H)0< 8.1 × 10-10 (95% confidence level).

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. This map shows 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 (â.," < 60), the major foreground contaminant is the diffuse thermal dust emission. At small angular scales (â.," > 500) the clustered cosmic infrared background and residual point sources are the major contaminants. These foregrounds are carefully modelled and subtracted. We thus measure the tSZ power spectrum over angular scales 0.17° ≠θ ≠3.0° that were previously unexplored. The measured tSZ power spectrum is consistent with that expected from the Planck catalogue of SZ sources, with clear evidence of additional signal from unresolved clusters and, potentially, diffuse warm baryons. 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. The measured tSZ power spectrum and high order statistics are used to place constraints on σ8.

We describe the data processing pipeline of the Planck Low Frequency Instrument (LFI) data processing centre (DPC) to create and characterize full-sky maps based on the first 15.5 months of operations at 30, 44, and 70 GHz. In particular, we discuss the various steps involved in reducing the data, from telemetry packets through to the production of cleaned, calibrated timelines and calibrated frequency maps. Data are continuously calibrated using the modulation induced on the mean temperature of the cosmic microwave background radiation by the proper motion of the spacecraft. Sky signals other than the dipole are removed by an iterative procedure based on simultaneous fitting of calibration parameters and sky maps. Noise properties are estimated from time-ordered data after the sky signal has been removed, using a generalized least squares map-making algorithm. A destriping code (Madam) is employed to combine radiometric data and pointing information into sky maps, minimizing the variance of correlated noise. Noise covariance matrices, required to compute statistical uncertainties on LFI and Planck products, are also produced. Main beams are estimated down to the â‰-20 dB level using Jupiter transits, which are also used for the geometrical calibration of the focal plane.

This paper characterizes the effective beams, the effective beam window functions and the associated errors for the Planck High Frequency Instrument (HFI) detectors. The effective beam is theangular response including the effect of the optics, detectors, data processing and the scan strategy. The window function is the representation of this beam in the harmonic domain which is required to recover an unbiased measurement of the cosmic microwave background angular power spectrum. The HFI is a scanning instrument and its effective beams are the convolution of: a) the optical response of the telescope and feeds; b) the processing of the time-ordered data and deconvolution of the bolometric and electronic transfer function; and c) the merging of several surveys to produce maps. The time response transfer functions are measured using observations of Jupiter and Saturn and by minimizing survey difference residuals. The scanning beam is the post-deconvolution angular response of the instrument, and is characterized with observations of Mars. The main beam solid angles are determined to better than 0.5% at each HFI frequency band. Observations of Jupiter and Saturn limit near sidelobes (within 5°) to about 0.1% of the total solid angle. Time response residuals remain as long tails in the scanning beams, but contribute less than 0.1% of the total solid angle. The bias and uncertainty in the beam products are estimated using ensembles of simulated planet observations that include the impact of instrumental noise and known systematic effects. The correlation structure of these ensembles is well-described by five error eigenmodes that are sub-dominant to sample variance and instrumental noise in the harmonic domain. A suite of consistency tests provide confidence that the error model represents a sufficient description of the data. The total error in the effective beam window functions is below 1% at 100 GHz up to multipole â.," ~ 1500, and below 0.5% at 143 and 217 GHz up to ~ 2000.

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 = 1.23 × 10-3, or v = 369. 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 10-3 effect applied to fluctuations which are already one part in roughly 105, 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°), of 384 km s-1 ± 78 km s-1 (stat.) ± 115 km s-1 (syst.). This is a significant confirmation of the expected velocity.

The two fundamental assumptions of the standard cosmological model-that the initial fluctuations are statistically isotropic and Gaussian-are rigorously tested using maps of the cosmic microwave background (CMB) anisotropy from the Planck satellite. The detailed results are based on studies of four independent estimates of the CMB that are compared to simulations using a fiducial ΛCDM model and incorporating essential aspects of the Planck measurement process. Deviations from isotropy have been found and demonstrated to be robust against component separation algorithm, mask choice, and frequency dependence. Many of these anomalies were previously observed in the WMAP data, and are now confirmed at similar levels of significance (about 3σ). However, we find little evidence of non-Gaussianity, with the exception of a few statistical signatures that seem to be associated with specific anomalies. In particular, we find that the quadrupole-octopole alignment is also connected to a low observed variance in the CMB signal. A power asymmetry is now found to persist on scales corresponding to about â.," = 600 and can be described in the low-â.," regime by a phenomenological dipole modulation model. However, any primordial power asymmetry is strongly scale-dependent and does not extend to arbitrarily small angular scales. Finally, it is plausible that some of these features may be reflected in the angular power spectrum of the data, which shows a deficit of power on similar scales. Indeed, when the power spectra of two hemispheres defined by a preferred direction are considered separately, one shows evidence of a deficit in power, while its opposite contains oscillations between odd and even modes that may be related to the parity violation and phase correlations also detected in the data. Although these analyses represent a step forward in building an understanding of the anomalies, a satisfactory explanation based on physically motivated models is still lacking.

This paper presents the characterization of the in-flight beams, the beam window functions, and the associated uncertainties for the Planck Low Frequency Instrument (LFI). Knowledge of the beam profiles is necessary for determining the transfer function to go from the observed to the actual sky anisotropy power spectrum. The main beam distortions affect the beam window function, complicating the reconstruction of the anisotropy power spectrum at high multipoles, whereas the sidelobes affect the low and intermediate multipoles. The in-flight assessment of the LFI main beams relies on the measurements performed during Jupiter observations. By stacking the datafrom multiple Jupiter transits, the main beam profiles are measured down to-20 dB at 30 and 44 GHz, and down to-25 dB at 70 GHz. The main beam solid angles are determined to better than 0.2% at each LFI frequency band. The Planck pre-launch optical model is conveniently tuned to characterize the main beams independently of any noise effects. This approach provides an optical model whose beams fully reproduce the measurements in the main beam region, but also allows a description of the beams at power levels lower than can be achieved by the Jupiter measurements themselves. The agreement between the simulated beams and the measured beams is better than 1% at each LFI frequency band. The simulated beams are used for the computation of the window functions for the effective beams. The error budget for the window functions is estimated from both main beam and sidelobe contributions, and accounts for the radiometer bandshapes. The total uncertainties in the effective beam window functions are: 2% and 1.2% at 30 and 44 GHz, respectively (at 600), and 0.7% at 70 GHz (at 1000).

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 be ns = 0.9603 ± 0.0073, ruling out exact scale invariance at over 5 s. 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 dns/dlnk =-0.0134 ± 0.0090. We verify these conclusions through a numerical analysis, which makes no slow-roll approximation, and carry out a Bayesian parameter estimation and model-selection analysis for a number of inflationary models including monomial, natural, and hilltop potentials. For each model, we present the Planck constraints on the parameters of the potential and explore several possibilities for the post-inflationary entropy generation epoch, thus obtaining nontrivial data-driven constraints. We also present a direct reconstruction of the observable range of the inflaton potential. Unless a quartic term is allowed in the potential, we find results consistent with second-order slow-roll predictions. We also investigate whether the primordial power spectrum contains any features. We find that models with a parameterized oscillatory feature improve the fit by 2eff 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 Planck bounds on fNL. Planck constrains with unprecedented accuracy the amplitude and possible correlation (with the adiabatic mode) of non-decaying isocurvature fluctuations. The fractional primordial contributions of cold dark matter (CDM) isocurvature modes of the types expected in the curvaton and axion scenarios have upper bounds of 0.25% and 3.9% (95% CL), respectively. In models with arbitrarily correlated CDM or neutrino isocurvature modes, an anticorrelated isocurvature component can improve the 2eff by approximately 4 as a result of slightly lowering the theoretical prediction for the l40 multipoles relative to the higher multipoles. Nonetheless, the data are consistent with adiabatic initial conditions.

This paper describes the methods used to produce photometrically calibrated maps from the Planck High Frequency Instrument (HFI) cleaned, time-ordered information. HFI observes the sky over a broad range of frequencies, from 100 to 857 GHz. To obtain the best calibration accuracy over such a large range, two different photometric calibration schemes have to be used. The 545 and 857 GHz data are calibrated by comparing flux-density measurements of Uranus and Neptune with models of their atmospheric emission. The lower frequencies (below 353 GHz) are calibrated using the solar dipole. A component of this anisotropy is time-variable, owing to the orbital motion of the satellite in the solar system. Photometric calibration is thus tightly linked to mapmaking, which also addresses low-frequency noise removal. By comparing observations taken more than one year apart in the same configuration, we have identified apparent gain variations with time. These variations are induced by non-linearities in the read-out electronics chain. We have developed an effective correction to limit their effect on calibration. We present several methods to estimate the precision of the photometric calibration. We distinguish relative uncertainties (between detectors, or between frequencies) and absolute uncertainties. Absolute uncertainties lie in the range from 0.54% to 10% from 100 to 857 GHz. We describe the pipeline used to produce the maps from the HFI timelines, based on the photometric calibration parameters, and the scheme used to set the zero level of the maps a posteriori. We also discuss the cross-calibration between HFI and the SPIRE instrument on board Herschel. Finally we summarize the basic characteristics of the set of HFI maps included in the 2013 Planck data release.

Planck data have been used to provide stringent new constraints on cosmic strings and other defects. We describe forecasts of the CMB power spectrum induced by cosmic strings, calculating these from network models and simulations using line-of-sight Boltzmann solvers. We have studied Nambu-Goto cosmic strings, as well as field theory strings for which radiative effects are important, thus spanning the range of theoretical uncertainty in the underlying strings models. We have added the angular power spectrum from strings to that for a simple adiabatic model, with the extra fraction defined as f10 at multipole = 10. This parameter has been added to the standard six parameter fit using COSMOMC with flat priors. For the Nambu-Goto string model, we have obtained a constraint on the string tension of Gμ/c2 < 1.5 × 10-7 and f10 < 0.015 at 95% confidence that can be improved to Gμ/c2 < 1.3 × 10-7 and f10 < 0.010 on inclusion of high-CMB data. For the Abelian-Higgs field theory model we find, GμAH/c2< 3.2 × 10-7 and f10 < 0.028. The marginalised likelihoods for f10 and in the f10-Ωbh2 plane are also presented. We have additionally obtained comparable constraints on f10 for models with semilocal strings and global textures. In terms of the effective defect energy scale these are somewhat weaker at Gμ/c2 < 1.1 × 10-6. We have made complementarity searches for the specific non-Gaussian signatures of cosmic strings, calibrating with all-sky Planck resolution CMB maps generated from networks of post-recombination strings. We have validated our non-Gaussian searches using these simulated maps in a Planck-realistic context, estimating sensitivities of up to ΔGμ/c2 ≠4 × 10-7. We have obtained upper limits on the string tension at 95% confidence of Gμ/c2 < 9.0 × 10-7 with modal bispectrum estimation and Gμ/c2 < 7.8 × 10-7 for real space searches with Minkowski functionals. These are conservative upper bounds because only post-recombination string contributions have been included in the non-Gaussian analysis.

The Planck Catalogue of Compact Sources (PCCS) is the catalogue of sources detected in the first 15 months of Planck operations, the "nominal" mission. It consists of nine single-frequency catalogues of compact sources, both Galactic and extragalactic, detected over the entire sky. The PCCS covers the frequency range 30-857 GHz with higher sensitivity (it is 90% complete at 180 mJy in the best channel) and better angular resolution (from 32.88′ to 4.33′) than previous all-sky surveys in this frequency band. By construction its reliability is >80% and more than 65% of the sources have been detected in at least two contiguous Planck channels. In this paper we present the construction and validation of the PCCS, its contents and its statistical characterization.

The Planck High Frequency Instrument (HFI) spectral response was determined through a series of ground based tests conducted with the HFI focal plane in a cryogenic environment prior to launch. The main goal of the spectral transmission tests was to measure the relative spectral response (includingthe level of out-of-band signal rejection) of all HFI detectors to a known source of electromagnetic radiation individually. This was determined by measuring the interferometric output of a continuously scanned Fourier transform spectrometer with all HFI detectors. As there is no on-board spectrometer within HFI, the ground-based spectral response experiments provide the definitive data set for the relative spectral calibration of the HFI. Knowledge of the relative variations in the spectral response between HFI detectors allows for a more thorough analysis of the HFI data. The spectral response of the HFI is used in Planck data analysis and component separation, this includes extraction of CO emission observed within Planck bands, dust emission, Sunyaev-Zeldovich sources, and intensity to polarization leakage. The HFI spectral response data have also been used to provide unit conversion and colour correction analysis tools. While previous papers describe the pre-flight experiments conducted on the Planck HFI, this paper focusses on the analysis of the pre-flight spectral response measurements and the derivation of data products, e.g. band-average spectra, unit conversion coefficients, and colour correction coefficients, all with related uncertainties. Verifications of the HFI spectral response data are provided through comparisons with photometric HFI flight data. This validation includes use of HFI zodiacal emission observations to demonstrate out-of-band spectral signal rejection better than 108. The accuracy of the HFI relative spectral response data is verified through comparison with complementary flight-data based unit conversion coefficients and colour correction coefficients. These coefficients include those based upon HFI observations of CO, dust, and Sunyaev-Zeldovich emission. General agreement is observed between the ground-based spectral characterization of HFI and corresponding in-flight observations, within the quoted uncertainty of each; explanations are provided for any discrepancies.

The Planck satellite provides a set of all-sky maps at nine frequencies from 30GHz to 857GHz. Planets, minor bodies, and diffuse interplanetary dust emission (IPD) are all observed. The IPD can be separated from Galactic and other emissions because Planck views a given point on the celestial sphere multiple times, through different columns of IPD. We use the Planck data to investigate the behaviour of zodiacal emission over the whole sky at sub-millimetre and millimetre wavelengths. We fit the Planck data to find the emissivities of the various components of the COBE zodiacal model-a diffuse cloud, three asteroidal dust bands, a circumsolar ring, and an Earth-trailing feature. The emissivity of the diffuse cloud decreases with increasing wavelength, as expected from earlier analyses. The emissivities of the dust bands, however, decrease less rapidly, indicating that the properties of the grains in the bands are different from those in the diffuse cloud. We fit the small amount of Galactic emission seen through the telescope's far sidelobes, and place limits on possible contamination of the cosmic microwave background (CMB) results from both zodiacal and far-sidelobe emission. When necessary, the results are used in the Planck pipeline to make maps with zodiacal emission and far sidelobes removed. We show that the zodiacal correction to the CMB maps is small compared to the Planck CMB temperature power spectrum and give a list of flux densities for small solar system bodies.

We describe the all-sky Planck catalogue of clusters and cluster candidates derived from Sunyaev-Zeldovich (SZ) effect detections using the first 15.5 months of Planck satellite observations. The catalogue contains 1227 entries, making it over six times the size of the Planck Early SZ (ESZ) sample and the largest SZ-selected catalogue to date. It contains 861 confirmed clusters, of which 178 have been confirmed as clusters, mostly through follow-up observations, and a further 683 are previously-known clusters. The remaining 366 have the status of cluster candidates, and we divide them into three classes according to the quality of evidence that they are likely to be true clusters. The Planck SZ catalogue is the deepest all-sky cluster catalogue, with redshifts up to about one, and spans the broadest cluster mass range from (0.1 to 1.6) × 1015 M⊙. Confirmation of cluster candidates through comparison with existing surveys or cluster catalogues is extensively described, as is the statistical characterization of the catalogue in terms of completeness and statistical reliability. The outputs of the validation process are provided as additional information. This gives, in particular, an ensemble of 813 cluster redshifts, and for all these Planck clusters we also include a mass estimated from a newly-proposed SZ-mass proxy. A refined measure of the SZ Compton parameter for the clusters with X-ray counter-parts is provided, as is an X-ray flux for all the Planck clusters not previously detected in X-ray surveys.

The European Space Agency's Planck satellite, dedicated to studying the early Universe and its subsequent evolution, was launched 14 May 2009 and has been scanning the microwave and submillimetre sky continuously since 12 August 2009. In March 2013, ESA and the Planck Collaboration released the initial cosmology products based on the first 15.5 months of Planck data, along with a set of scientific and technical papers and a web-based explanatory supplement. This paper gives an overview of the mission and its performance, the processing, analysis, and characteristics of the data, the scientific results, and the science data products and papers in the release. The science products include maps of the cosmic microwave background (CMB) and diffuse extragalactic foregrounds, a catalogue of compact Galactic and extragalactic sources, and a list of sources detected through the Sunyaev-Zeldovich effect. The likelihood code used to assess cosmological models against the Planck data and a lensing likelihood are described. Scientific results include robust support for the standard six-parameter ΛCDM model of cosmology and improved measurements of its parameters, including a highly significant deviation from scale invariance of the primordial power spectrum. The Planck values for these parameters and others derived from them are significantly different from those previously determined. Several large-scale anomalies in the temperature distribution of the CMB, first detected by WMAP, are confirmed with higher confidence. Planck sets new limits on the number and mass of neutrinos, and has measured gravitational lensing of CMB anisotropies at greater than 25σ. Planck finds no evidence for non-Gaussianity in the CMB. Planck's results agree well with results from the measurements of baryon acoustic oscillations. Planck finds a lower Hubble constant than found in some more local measures. Some tension is also present between the amplitude of matter fluctuations (σ8) derived from CMB data and that derived from Sunyaev-Zeldovich data. The Planck and WMAP power spectra are offset from each other by an average level of about 2% around the first acoustic peak. Analysis of Planck polarization data is not yet mature, therefore polarization results are not released, although the robust detection of E-mode polarization around CMB hot and cold spots is shown graphically.

We present the current estimate of instrumental and systematic effect uncertainties for the Planck-Low Frequency Instrument relevant to the first release of the Planck cosmological results. We give an overview of the main effects and of the tools and methods applied to assess residuals in maps and power spectra. We also present an overall budget of known systematic effect uncertainties, which are dominated by sidelobe straylight pick-up and imperfect calibration. However, even these two effects are at least two orders of magnitude weaker than the cosmic microwave background fluctuations as measured in terms of the angular temperature power spectrum. A residual signal above the noise level is present in the multipole range â.," < 20, most notably at 30 GHz, and is probably caused by residual Galactic straylight contamination. Current analysis aims to further reduce the level of spurious signals in the data and to improve the systematic effects modelling, in particular with respect to straylight and calibration uncertainties.

We present new measurements of cosmic infrared background (CIB) anisotropies using Planck. Combining HFI data with IRAS, the angular auto-and cross-frequency power spectrum is measured from 143 to 3000 GHz, and the auto-bispectrum from 217 to 545 GHz. The total areas used to compute the CIB power spectrum and bispectrum are about 2240 and 4400 deg2, respectively. After careful removal of the contaminants (cosmic microwave background anisotropies, Galactic dust, and Sunyaev-Zeldovich emission), and a complete study of systematics, the CIB power spectrum is measured with unprecedented signal to noise ratio from angular multipoles â.," ~ 150 to 2500. The bispectrum due to the clustering of dusty, star-forming galaxies is measured from â.," ~ 130 to 1100, with a total signal to noise ratio of around 6, 19, and 29 at 217, 353, and 545 GHz, respectively. Two approaches are developed for modelling CIB power spectrum anisotropies. The first approach takes advantage of the unique measurements by Planck at large angular scales, and models only the linear part of the power spectrum, with a mean bias of dark matter haloes hosting dusty galaxies at a given redshift weighted by their contribution to the emissivities. The second approach is based on a model that associates star-forming galaxies with dark matter haloes and their subhaloes, using a parametrized relation between the dust-processed infrared luminosity and (sub-)halo mass. The two approaches simultaneously fit all auto-and cross-power spectra very well. We find that the star formation history is well constrained up to redshifts around 2, and agrees with recent estimates of the obscured star-formation density using Spitzer and Herschel. However, at higher redshift, the accuracy of the star formation history measurement is strongly degraded by the uncertainty in the spectral energy distribution of CIB galaxies. We also find that the mean halo mass which is most efficient at hosting star formation is log (Meff/M) = 12.6 and that CIB galaxies have warmer temperatures as redshift increases. The CIB bispectrum is steeper than that expected from the power spectrum, although well fitted by a power law; this gives some information about the contribution of massive haloes to the CIB bispectrum. Finally, we show that the same halo occupation distribution can fit all power spectra simultaneously. The precise measurements enabled by Planck pose new challenges for the modelling of CIB anisotropies, indicating the power of using CIB anisotropies to understand the process of galaxy formation.

The Planck nominal mission cosmic microwave background (CMB) maps yield unprecedented constraints on primordial non-Gaussianity (NG). Using three optimal bispectrum estimators, separable template-fitting (KSW), binned, and modal, we obtain consistent values for the primordial local, equilateral, and orthogonal bispectrum amplitudes, quoting as our final result fNLlocal = 2.7 ± 5.8, fNLequil = -42 ± 75, and fNLorth = -25 ± 39 (68% CL statistical). Non-Gaussianity is detected in the data; using skew-Cl statistics we find a nonzero bispectrum from residual point sources, and the integrated-Sachs-Wolfe-lensing bispectrum at a level expected in the ?CDM scenario. The results are based on comprehensive cross-validation of these estimators on Gaussian and non-Gaussian simulations, are stable across component separation techniques, pass an extensive suite of tests, and are confirmed by skew-Cl, wavelet bispectrum and Minkowski functional estimators. Beyond estimates of individual shape amplitudes, we present model-independent, three-dimensional reconstructions of the Planck CMB bispectrum and thus derive constraints on early-Universe scenarios that generate primordial NG, including general single-field models of inflation, excited initial states (non-Bunch-Davies vacua), and directionally-dependent vector models. We provide an initial survey of scale-dependent feature and resonance models. These results bound both general single-field and multi-field model parameter ranges, such as the speed of sound, cs ≥ 0.02 (95% CL), in an effective field theory parametrization, and the curvaton decay fraction rD ≥ 0.15 (95% CL). The Planck data significantly limit the viable parameter space of the ekpyrotic/cyclic scenarios. The amplitude of the four-point function in the local model τNL< 2800 (95% CL). Taken together, these constraints represent the highest precision tests to date of physical mechanisms for the origin of cosmic structure.

We present constraints on cosmological parameters using number counts as a function of redshift for a sub-sample of 189 galaxy clusters from the Planck SZ (PSZ) catalogue. The PSZ is selected through the signature of the Sunyaev-Zeldovich (SZ) effect, and the sub-sample used here has a signal-to-noise threshold of seven, with each object confirmed as a cluster and all but one with a redshift estimate. We discuss the completeness of the sample and our construction of a likelihood analysis. Using a relation between mass M and SZ signal Y calibrated to X-ray measurements, we derive constraints on the power spectrum amplitude σ8 and matter density parameter Ωm in a flat ΛCDM model. We test the robustness of our estimates and find that possible biases in the Y-M relation and the halo mass function are larger than the statistical uncertainties from the cluster sample. Assuming the X-ray determined mass to be biased low relative to the true mass by between zero and 30%, motivated by comparison of the observed mass scaling relations to those from a set of numerical simulations, we find that σ8 = 0.75 ± 0.03, Ωm = 0.29 ± 0.02, and σ8(Ωm/0.27)0.3 = 0.764 ± 0.025. The value of σ8 is degenerate with the mass bias; if the latter is fixed to a value of 20% (the central value from numerical simulations) we find σ8(Ωm/0.27)0.3 = 0.78 ± 0.01 and a tighter one-dimensional range σ8 = 0.77 ± 0.02. We find that the larger values of σ8 and Ωm preferred by Planck's measurements of the primary CMB anisotropies can be accommodated by a mass bias of about 40%. Alternatively, consistency with the primary CMB constraints can be achieved by inclusion of processes that suppress power on small scales relative to the ΛCDM model, such as a component of massive neutrinos. We place our results in the context of other determinations of cosmologicalparameters, and discuss issues that need to be resolved in order to make further progress in this field.

On the arcminute angular scales probed by Planck, the cosmic microwave background (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 217 GHz frequency bands with an overall significance of greater than 25σ. We use thetemperature-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 ΛCDM 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 root-mean-squared amplitude of matter fluctuations at z ~ 2.

Planck has produced detailed all-sky observations over nine frequency bands between 30 and 857 GHz. These observations allow robust reconstruction of the primordial cosmic microwave background (CMB) temperature fluctuations over nearly the full sky, as well as new constraints on Galactic foregrounds, including thermal dust and line emission from molecular carbon monoxide (CO). This paper describes the component separation framework adopted by Planck for many cosmological analyses, including CMB power spectrum determination and likelihood construction on large angular scales, studies of primordial non-Gaussianity and statistical isotropy, the integrated Sachs-Wolfe effect, gravitational lensing, and searches for topological defects. We test four foreground-cleaned CMB maps derived using qualitatively different component separation algorithms. The quality of our reconstructions is evaluated through detailed simulations and internal comparisons, and shown through various tests to be internally consistent and robust for CMB power spectrum and cosmological parameter estimation up to â.," = 2000. The parameter constraints on ΛCDM cosmologies derived from these maps are consistent with those presented in the cross-spectrum based Planck likelihood analysis. We choose two of the CMB maps for specific scientific goals. We also present maps and frequency spectra of the Galactic low-frequency, CO, and thermal dust emission. The component maps are found to provide a faithful representation of the sky, as evaluated by simulations, with the largest bias seen in the CO component at 3%. For the low-frequency component, the spectral index varies widely over the sky, ranging from about β =-4 to-2. Considering both morphology and prior knowledge of the low frequencycomponents, the index map allows us to associate a steep spectral index (β-2.3 with strong free-free emission, and intermediate values with synchrotron emission.

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 cosmic microwave background (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 releasewe rely on the characterization of the solar dipole as measured by WMAP. We also present preliminary results (at 44 GHz 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 30 GHz radiometers (the ones with the largest sidelobes) using the signal of an internal reference load at 4 K instead of the CMB dipole. We have estimated the accuracy of the LFI calibration following two strategies: (1) 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 (2) we have performed a number of internal consistency checks on the data and on the brightness temperature of Jupiter. Errors in the calibration of this Planck/LFI data release are expected to be about 0.6% at 44 and 70 GHz, and 0.8% at 30 GHz. Both these preliminary results at low and high â.," are consistent with WMAP results within uncertainties and comparison of power spectra indicates good consistency in the absolute calibration with HFI (0.3%) and a 1.4σ discrepancy with WMAP (0.9%).

Based on cosmic microwave background (CMB) maps from the 2013 Planck Mission data release, this paper presents the detection of the integrated Sachs-Wolfe (ISW) effect, that is, the correlation between the CMB and large-scale evolving gravitational potentials. The significance of detection ranges from 2 to 4σ, depending on which method is used. We investigated three separate approaches, which essentially cover all previous studies, and also break new ground. (i) We correlated the CMB with the Planck reconstructed gravitational lensing potential (for the first time). This detection was made using the lensing-induced bispectrum between the low-â.," and high-â.," temperature anisotropies; the correlation between lensing and the ISW effect has a significance close to 2.5σ. (ii) We cross-correlated with tracers of large-scale structure, which yielded a significance of about 3σ, based on a combination of radio (NVSS) and optical (SDSS) data. (iii) We used aperture photometry on stacked CMB fields at the locations of known large-scale structures, which yielded and confirms a 4σ signal, over a broader spectral range, when using a previously explored catalogue, but shows strong discrepancies in amplitude and scale when compared with expectations. More recent catalogues give more moderate results that range from negligible to 2.5σ at most, but have a more consistent scale and amplitude, the latter being still slightly higher than what is expected from numerical simulations within ΛCMD. Where they can be compared, these measurements are compatible with previous work using data from WMAP, where these scales have been mapped to the limits of cosmic variance. Planck's broader frequency coverage allows for better foreground cleaning and confirms that the signal is achromatic, which makes it preferable for ISW detection. As a final step we used tracers of large-scale structure to filter the CMB data, from which we present maps of the ISW temperature perturbation. These results provide complementary and independent evidence for the existence of a dark energy component that governs the currently accelerated expansion of the Universe.

This paper presents the first cosmological results based on Planck measurements of the cosmic microwave background (CMB) temperature and lensing-potential power spectra. We find that the Planck spectra at high multipoles (40) are extremely well described by the standard spatially-flat six-parameter λCDM cosmology with a power-law spectrum of adiabatic scalar perturbations. Within the context of this cosmology, the Planck data determine the cosmological parameters to high precision: the angular size of the sound horizon at recombination, the physical densities of baryons and cold dark matter, and the scalar spectral index are estimated to be θ- = (1.04147 ± 0.00062) × 10-2, Ωbh2 = 0.02205 ± 0.00028, Ωch2 = 0.1199 ± 0.0027, and ns = 0.9603 ± 0.0073, respectively(note that in this abstract we quote 68% errors on measured parameters and 95% upper limits on other parameters). For this cosmology, we find a low value of the Hubble constant, H0 = (67.3 ± 1.2) km s-1 Mpc-1, and a high value of the matter density parameter, Ωm = 0.315 ± 0.017. These values are in tension with recent direct measurements of H0 and the magnitude-redshift relation for Type Ia supernovae, but are in excellent agreement with geometrical constraints from baryon acoustic oscillation (BAO) surveys. Including curvature, we find that the Universe is consistent with spatial flatness to percent level precision using Planck CMB data alone. We use high-resolution CMB data together with Planck to provide greater control on extragalactic foreground components in an investigation of extensions to the six-parameter λCDM model. We present selected results from a large grid of cosmological models, using a range of additional astrophysical data sets in addition to Planck and high-resolution CMB data. None of these models are favoured over the standard six-parameter λCDM cosmology. The deviation of the scalar spectral index from unity isinsensitive to the addition of tensor modes and to changes in the matter content of the Universe. We find an upper limit of r0.002 < 0.11 on the tensor-to-scalar ratio. There is no evidence for additional neutrino-like relativistic particles beyond the three families of neutrinos in the standard model. Using BAO and CMB data, we find Neff = 3.30 ± 0.27 for the effective number of relativistic degrees of freedom, and an upper limit of 0.23 eV for the sum of neutrino masses. Our results are in excellent agreement with big bang nucleosynthesis and the standard value of Neff = 3.046. We find no evidence for dynamical dark energy; using BAO and CMB data, the dark energy equation of state parameter is constrained to be w =-1.13-0.10+0.13. We also use the Planck data to set limits on a possible variation of the fine-structure constant, dark matter annihilation and primordial magnetic fields. Despite the success of the six-parameter λCDM model in describing the Planck data at high multipoles, we note that this cosmology does not provide a good fit to the temperature power spectrum at low multipoles. The unusual shape of the spectrum in the multipole range 20 40 was seen previously in the WMAP data and is a real feature of the primordial CMB anisotropies. The poor fit to the spectrum at low multipoles is not of decisive significance, but is an "anomaly" in an otherwise self-consistent analysis of the Planck temperature data.

© 2014 ESO. We describe the processing of the 531 billion raw data samples from the High Frequency Instrument (HFI), which we performed to produce six temperature maps from the first 473 days of Planck-HFI survey data. These maps provide an accurate rendition of the sky emission at 100, 143, 217, 353, 545, and 857 GHz with an angular resolution ranging from 9.I 7 to 4.I 6. The detector noise per (effective) beam solid angle is respectively, 10, 6, 12, and 39 μK in the four lowest HFI frequency channels (100-353 GHz) and 13 and 14 kJy sr-1 in the 545 and 857 GHz channels. Relative to the 143 GHz channel, these two high frequency channels are calibrated to within 5% and the 353 GHz channel to the percent level. The 100 and 217 GHz channels, which together with the 143 GHz channel determine the high-multipole part of the CMB power spectrum (50., 2500), are calibrated relative to 143 GHz to better than 0.2%.

The dust-Hi correlation is used to characterize the emission properties of dust in the diffuse interstellar medium (ISM) from far infrared wavelengths to microwave frequencies. The field of this investigation encompasses the part of the southern sky best suited to study the cosmic infrared and microwave backgrounds. We cross-correlate sky maps from Planck, the Wilkinson Microwave Anisotropy Probe (WMAP), and the diffuse infrared background experiment (DIRBE), at 17 frequencies from 23 to 3000 GHz, with the Parkes survey of the 21 cm line emission of neutral atomic hydrogen, over a contiguous area of 7500 deg2 centred on the southern Galactic pole. We present a general methodology to study the dust-H i correlation over the sky, including simulations to quantify uncertainties. Our analysis yields four specific results. (1) We map the temperature, submillimetre emissivity, and opacity of the dust per H-atom. The dust temperature is observed to be anti-correlated with the dust emissivity and opacity. We interpret this result as evidence of dust evolution within the diffuse ISM. The mean dust opacity is measured to be (7.1 ± 0.6) × 10-27 cm2 H-1 × (v/353 GHz) 1.53 ± 0.03for 100 ≤ v ≤ 353 GHz. This is a reference value to estimate hydrogen column densities from dust emission at submillimetre and millimetre wavelengths. (2) We map the spectral index βmm of dust emission at millimetre wavelengths (defined here as v ≤ 353GHz), and find it to be remarkably constant at βmm = 1.51 ± 0.13. We compare it with the far infrared spectral index βFIR derived from greybody fits at higher frequencies, and find a systematic difference, βmm -βFIR = -0.15, which suggests that the dust spectral energy distribution (SED) flattens at v ≤ 353 GHz. (3) We present spectral fits of the microwave emission correlated with Hi from 23 to 353 GHz, which separate dust and anomalous microwave emission (AME). We show that the flattening of the dust SED can be accounted for with an additional component with a blackbody spectrum. This additional component, which accounts for (26 ± 6)% of the dust emission at 100GHz, could represent magnetic dipole emission. Alternatively, it could account for an increasing contribution of carbon dust, or a flattening of the emissivity of amorphous silicates, at millimetre wavelengths. These interpretations make different predictions for the dust polarization SED. (4) We analyse the residuals of the dust-Hi correlation. We identify a Galactic contribution to these residuals, which we model with variations of the dust emissivity on angular scales smaller than that of our correlation analysis. This model of the residuals is used to quantify uncertainties of the CIB power spectrum in a companion Planck paper.© ESO 2014.

We explore the 2013 Planck likelihood function with a high-precision multi-dimensional minimizer (Minuit). This allows a refinement of the CDM best-fit solution with respect to previously-released results, and the construction of frequentist confidence intervals using profile likelihoods. The agreement with the cosmological results from the Bayesian framework is excellent, demonstrating the robustness of the Planck results to the statistical methodology. We investigate the inclusion of neutrino masses, where more significant differences may appear due to the non-Gaussian nature of the posterior mass distribution. By applying the Feldman-Cousins prescription, we again obtain results very similar to those of the Bayesian methodology. However, the profile-likelihood analysis of the cosmic microwave background (CMB) combination (Planck+WP+highL) reveals a minimum well within the unphysical negative-mass region. We show that inclusion of the Planck CMB-lensing information regularizes this issue, and provide a robust frequentist upper limit σmv ≤0:26 eV (95% confidence) from the CMB+lensing+BAO data combination. © ESO 2014.

The dust-Hi correlation is used to characterize the emission properties of dust in the diffuse interstellar medium (ISM) from far infrared wavelengths to microwave frequencies. The field of this investigation encompasses the part of the southern sky best suited to study the cosmic infrared and microwave backgrounds. We cross-correlate sky maps from Planck, the Wilkinson Microwave Anisotropy Probe (WMAP), and the diffuse infrared background experiment (DIRBE), at 17 frequencies from 23 to 3000 GHz, with the Parkes survey of the 21 cm line emission of neutral atomic hydrogen, over a contiguous area of 7500 deg2 centred on the southern Galactic pole. We present a general methodology to study the dust-H i correlation over the sky, including simulations to quantify uncertainties. Our analysis yields four specific results. (1) We map the temperature, submillimetre emissivity, and opacity of the dust per H-atom. The dust temperature is observed to be anti-correlated with the dust emissivity and opacity. We interpret this result as evidence of dust evolution within the diffuse ISM. The mean dust opacity is measured to be (7.1 ± 0.6) × 10-27 cm2 H-1 × (v/353 GHz) 1.53 ± 0.03for 100 ≤ v ≤ 353 GHz. This is a reference value to estimate hydrogen column densities from dust emission at submillimetre and millimetre wavelengths. (2) We map the spectral index βmm of dust emission at millimetre wavelengths (defined here as v ≤ 353GHz), and find it to be remarkably constant at βmm = 1.51 ± 0.13. We compare it with the far infrared spectral index βFIR derived from greybody fits at higher frequencies, and find a systematic difference, βmm -βFIR = -0.15, which suggests that the dust spectral energy distribution (SED) flattens at v ≤ 353 GHz. (3) We present spectral fits of the microwave emission correlated with Hi from 23 to 353 GHz, which separate dust and anomalous microwave emission (AME). We show that the flattening of the dust SED can be accounted for with an additional component with a blackbody spectrum. This additional component, which accounts for (26 ± 6)% of the dust emission at 100GHz, could represent magnetic dipole emission. Alternatively, it could account for an increasing contribution of carbon dust, or a flattening of the emissivity of amorphous silicates, at millimetre wavelengths. These interpretations make different predictions for the dust polarization SED. (4) We analyse the residuals of the dust-Hi correlation. We identify a Galactic contribution to these residuals, which we model with variations of the dust emissivity on angular scales smaller than that of our correlation analysis. This model of the residuals is used to quantify uncertainties of the CIB power spectrum in a companion Planck paper.© ESO 2014.

We explore the 2013 Planck likelihood function with a high-precision multi-dimensional minimizer (Minuit). This allows a refinement of the CDM best-fit solution with respect to previously-released results, and the construction of frequentist confidence intervals using profile likelihoods. The agreement with the cosmological results from the Bayesian framework is excellent, demonstrating the robustness of the Planck results to the statistical methodology. We investigate the inclusion of neutrino masses, where more significant differences may appear due to the non-Gaussian nature of the posterior mass distribution. By applying the Feldman-Cousins prescription, we again obtain results very similar to those of the Bayesian methodology. However, the profile-likelihood analysis of the cosmic microwave background (CMB) combination (Planck+WP+highL) reveals a minimum well within the unphysical negative-mass region. We show that inclusion of the Planck CMB-lensing information regularizes this issue, and provide a robust frequentist upper limit σmv ≤0:26 eV (95% confidence) from the CMB+lensing+BAO data combination. © ESO 2014.