Search

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.

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.

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.

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.

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.

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 K2 at 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.

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.

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.

We perform an analysis of the diffuse low-frequency Galactic components inthe Southern part of the Gould Belt system (130^\circ\leq l\leq 230^\circ and-50^\circ\leq b\leq -10^\circ). Strong ultra-violet (UV) flux coming from theGould Belt super-association is responsible for bright diffuse foregrounds thatwe observe from our position inside the system and that can help us improve ourknowledge of the Galactic emission. Free-free emission and anomalous microwaveemission (AME) are the dominant components at low frequencies (u < 40 GHz),while synchrotron emission is very smooth and faint. We separate diffusefree-free emission and AME from synchrotron emission and thermal dust emissionby using Planck data, complemented by ancillary data, using the "CorrelatedComponent Analysis" (CCA) component separation method and we compare with theresults of cross-correlation of foreground templates with the frequency maps.We estimate the electron temperature T_e from H$\alpha$ and free-free emissionusing two methods (temperature-temperature plot and cross-correlation) and weobtain T_e ranging from 3100 to 5200 K, for an effective fraction of absorbingdust along the line of sight of 30% (f_d=0.3). We estimate the frequencyspectrum of the diffuse AME and we recover a peak frequency (in flux densityunits) of 25.5 \pm 1.5 GHz. We verify the reliability of this result withrealistic simulations that include the presence of biases in the spectral modelfor the AME and in the free-free template. By combining physical models forvibrational and rotational dust emission and adding the constraints from thethermal dust spectrum from Planck and IRAS we are able to get a gooddescription of the frequency spectrum of the AME for plausible values of thelocal density and radiation field.

We present the scaling relation between Sunyaev-Zeldovich (SZ) signal and stellar mass for almost 260,000 locally brightest galaxies (LBGs) selected from the Sloan Digital Sky Survey (SDSS). These are predominantly the central galaxies of their dark matter halos. We calibrate the stellar-to-halo mass conversion using realistic mock catalogues based on the Millennium Simulation. Applying a multi-frequency matched filter to the Planck data for each LBG, and averaging the results in bins of stellar mass, we measure the mean SZ signal down to M500 ∼ 2 × 1011 M⊙, with a clear indication of signal at even lower stellar mass. We derive the scaling relation between SZ signal and halo mass by assigning halo properties from our mock catalogues to the real LBGs and simulating the Planck observation process. This relation shows no evidence for deviation from a power law over a halo mass range extending from rich clusters down to M500 ~ 2 × 10 13 M⊙, and there is a clear indication of signal down to M500 ~ 4 × 1012 M ⊙. Planck's SZdetections in such low-mass halos imply that about a quarter of all baryons have now been seen in the form of hot halo gas, and that this gas must be less concentrated than the dark matter in such halos in order to remain consistent with X-ray observations. At the high-mass end, the measured SZ signal is 20% lower than found from observations of X-ray clusters, a difference consistent with the magnitude of Malmquist bias effects that were previously estimated for the X-ray sample. © 2013 ESO,.