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.

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.

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.

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.

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

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 potential for Planck to detect clusters of dusty, star-forming galaxies at z > 1 is tested by examining the Herschel-SPIRE images of Planck Early Release Compact Source Catalog sources lying in fields observed by the Herschel Multitiered Extragalactic Survey. Of the 16 Planck sources that lie in the ~90 sq. deg. examined, we find that 12 are associated with single bright Herschel sources. The remaining four are associated with overdensities of Herschel sources, making them candidate clusters of dusty, star-forming galaxies. We use complementary optical/near-IR data for these 'clumps' to test this idea, and find evidence for the presence of galaxy clusters in all four cases. We use photometric redshiftsand red sequence galaxies to estimate the redshifts of these clusters, finding that they range from 0.8 to 2.3. These redshifts imply that the Herschel sources in these clusters, which contribute to the detected Planck flux, are forming stars very rapidly, with typical total cluster star formation rates > 1000M⊙ yr-1. The high-redshift clusters discovered inthese observations are used to constrain the epoch of cluster galaxy formation, finding that the galaxies in our clusters are 1-1.5 Gyr old at z ~ 1-2. Prospects for the discovery of further clusters of dusty galaxies are discussed, using not only all sky Planck surveys, but also deeper, smaller area, Herschel surveys. © 2014 The Authors Published by Oxford University Press on behalf of the Royal Astronomical Society.

We present a study of the infrared properties for a sample of seven spectroscopically confirmed submillimeter galaxies (SMGs) at z > 4.0. By combining ground-based near-infrared, Spitzer IRAC and MIPS, Herschel SPIRE, and ground-based submillimeter/millimeter photometry, we construct their spectral energy distributions (SEDs) and a composite model to fit the SEDs. The model includes a stellar emission component at λrest < 3.5 μm, a hot dust component peaking at λrest ∼ 5 μm, and cold dust component which becomes significant for λrest > 50 μm. Six objects in the sample are detected at 250 and 350 μm. The dust temperatures for the sources in this sample are in the range of 40-80 K, and their LFIR ∼ 1013 L⊙ qualifies them as hyper-luminous infrared galaxies. The mean FIR-radio index for this sample is around 〈q〉 = 2.2 indicating no radio excess in their radio emission. Most sources in the sample have 24 μm detections corresponding to a rest-frame 4.5 μm luminosity of Log10(L4.5/L ⊙) = 11 ∼ 11.5. Their L4.5/LFIR ratios are very similar to those of starburst-dominated SMGs at z ∼ 2. The L CO-LFIR relation for this sample is consistent with that determined for local ULIRGs and SMGs at z ∼ 2. We conclude that SMGs at z > 4 are hotter and more luminous in the FIR but otherwise very similar to those at z ∼ 2. None of these sources show any sign of the strong QSO phase being triggered. © 2014. The American Astronomical Society. All rights reserved.

We present a method for selecting z > 4 dusty, star-forming galaxies (DSFGs) using Herschel/Spectral and Photometric Imaging Receiver 250/350/500 μm flux densities to search for red sources. We apply this method to 21 deg2 of data from the HerMES survey to produce a catalog of 38 high-z candidates. Follow-up of the first five of these sources confirms that this method is efficient at selecting high-z DSFGs, with 4/5 at z = 4.3-6.3 (and the remaining source at z = 3.4), and that they are some of the most luminous dusty sources known. Comparison with previous DSFG samples, mostly selected at longer wavelengths (e.g., 850 μm) and in single-band surveys, shows that our method is much more efficient at selecting high-z DSFGs, in the sense that a much larger fraction are at z > 3. Correcting for the selection completeness and purity, we find that the number of bright (S 500 μm ≥ 30 mJy), red Herschel sources is 3.3 ± 0.8 deg-2. This is much higher than the number predicted by current models, suggesting that the DSFG population extends to higher redshifts than previously believed. If the shape of the luminosity function for high-z DSFGs is similar to that at z ∼ 2, rest-frame UV based studies may be missing a significant component of the star formation density at z = 4-6, even after correction for extinction. © 2014. The American Astronomical Society. All rights reserved.

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.

Strong gravitational lenses are now being routinely discovered in wide-field surveys at (sub-)millimeter wavelengths. We present Submillimeter Array (SMA) high-spatial resolution imaging and Gemini-South and Multiple Mirror Telescope optical spectroscopy of strong lens candidates discovered in the two widest extragalactic surveys conducted by the Herschel Space Observatory: the Herschel-Astrophysical Terahertz Large Area Survey (H-ATLAS) and the Herschel Multi-tiered Extragalactic Survey (HerMES). From a sample of 30 Herschel sources with S 500 > 100 mJy, 21 are strongly lensed (i.e., multiply imaged), 4 are moderately lensed (i.e., singly imaged), and the remainder require additional data to determine their lensing status. We apply a visibility-plane lens modeling technique to the SMA data to recover information about the masses of the lenses as well as the intrinsic (i.e., unlensed) sizes (r half) and far-infrared luminosities (L FIR) of the lensed submillimeter galaxies (SMGs). The sample of lenses comprises primarily isolated massive galaxies, but includes some groups and clusters as well. Several of the lenses are located at z lens > 0.7, a redshift regime that is inaccessible to lens searches based on Sloan Digital Sky Survey spectroscopy. The lensed SMGs are amplified by factors that are significantly below statistical model predictions given the 500 μm flux densities of our sample. We speculate that this may reflect a deficiency in our understanding of the intrinsic sizes and luminosities of the brightest SMGs. The lensed SMGs span nearly one decade in L FIR (median L FIR = 7.9 × 10 12 L ) and two decades in FIR luminosity surface density (median ΣFIR = 6.0 × 1011 L kpc-2). The strong lenses in this sample and others identified via (sub-)mm surveys will provide a wealth of information regarding the astrophysics of galaxy formation and evolution over a wide range in redshift. © 2013. The American Astronomical Society. All rights reserved..

Strong gravitational lenses are now being routinely discovered in wide-field surveys at (sub-)millimeter wavelengths. We present Submillimeter Array (SMA) high-spatial resolution imaging and Gemini-South and Multiple Mirror Telescope optical spectroscopy of strong lens candidates discovered in the two widest extragalactic surveys conducted by the Herschel Space Observatory: the Herschel-Astrophysical Terahertz Large Area Survey (H-ATLAS) and the Herschel Multi-tiered Extragalactic Survey (HerMES). From a sample of 30 Herschel sources with S 500 > 100 mJy, 21 are strongly lensed (i.e., multiply imaged), 4 are moderately lensed (i.e., singly imaged), and the remainder require additional data to determine their lensing status. We apply a visibility-plane lens modeling technique to the SMA data to recover information about the masses of the lenses as well as the intrinsic (i.e., unlensed) sizes (r half) and far-infrared luminosities (L FIR) of the lensed submillimeter galaxies (SMGs). The sample of lenses comprises primarily isolated massive galaxies, but includes some groups and clusters as well. Several of the lenses are located at z lens > 0.7, a redshift regime that is inaccessible to lens searches based on Sloan Digital Sky Survey spectroscopy. The lensed SMGs are amplified by factors that are significantly below statistical model predictions given the 500 μm flux densities of our sample. We speculate that this may reflect a deficiency in our understanding of the intrinsic sizes and luminosities of the brightest SMGs. The lensed SMGs span nearly one decade in L FIR (median L FIR = 7.9 × 10 12 L ) and two decades in FIR luminosity surface density (median ΣFIR = 6.0 × 1011 L kpc-2). The strong lenses in this sample and others identified via (sub-)mm surveys will provide a wealth of information regarding the astrophysics of galaxy formation and evolution over a wide range in redshift. © 2013. The American Astronomical Society. All rights reserved..

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 (\nu < 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,.

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,.