We present an analysis of Planck satellite data on the Coma Cluster observed via the Sunyaev-Zeldovich effect. Planck is able, for the first time, to detect SZ emission up to r ~ 3 X R_500. We test previously proposed models for the pressure distribution in clusters against the azimuthally averaged data. We find that the Arnaud et al. universal pressure profile does not fit Coma, and that their pressure profile for merging systems provides a good fit of the data only at r<R_500: by r=2XR_500 it underestimates the observed y profile by a factor of ~2. This may indicate that at these larger radii either i) the cluster SZ emission is contaminated by unresolved SZ sources along the line of sight or ii) the pressure profile of Coma is higher at r>R_500 than the mean pressure profile predicted by the simulations. The Planck image shows significant local steepening of the y profile in two regions about half a degree to the west and to the south-east of the cluster centre. These features are consistent with the presence of shock fronts at these radii, and indeed the western feature was previously noticed in the ROSAT PSPC mosaic as well as in the radio. Using Planck y profiles extracted from corresponding sectors we find pressure jumps of 4.5+0.4-0.2 and 5.0+1.3-0.1 in the west and southeast, respectively. Assuming Rankine-Hugoniot pressure jump conditions, we deduce that the shock waves should propagate with Mach number M_w=2.03+0.09-0.04 and M_se=2.05+0.25-0.02 in the West and Southeast, respectively. Finally, we find that the y and radio-synchrotron signals are quasi-linearly correlated on Mpc scales with small intrinsic scatter. This implies either that the energy density of cosmic-ray electrons is relatively constant throughout the cluster, or that the magnetic fields fall off much more slowly with radius than previously thought.

We examine the relation between the galaxy cluster mass M and Sunyaev-Zeldovich (SZ) effect signal D_A^2 Y for a sample of 19 objects for which weak lensing (WL) mass measurements obtained from Subaru Telescope data are available in the literature. Hydrostatic X-ray masses are derived from XMM-Newton archive data and the SZ effect signal is measured from Planck all-sky survey data. We find an M_WL-D_A^2 Y relation that is consistent in slope and normalisation with previous determinations using weak lensing masses; however, there is a normalisation offset with respect to previous measures based on hydrostatic X-ray mass-proxy relations. We verify that our SZ effect measurements are in excellent agreement with previous determinations from Planck data. For the present sample, the hydrostatic X-ray masses at R_500 are on average ~ 20 per cent larger than the corresponding weak lensing masses, at odds with expectations. We show that the mass discrepancy is driven by a difference in mass concentration as measured by the two methods, and, for the present sample, the mass discrepancy and difference in mass concentration is especially large for disturbed systems. The mass discrepancy is also linked to the offset in centres used by the X-ray and weak lensing analyses, which again is most important in disturbed systems. We outline several approaches that are needed to help achieve convergence in cluster mass measurement with X-ray and weak lensing observations.

Using precise full-sky observations from Planck, and applying several methods of component separation, we identify and characterize the emission from the Galactic "haze" at microwave wavelengths. The haze is a distinct component of diffuse Galactic emission, roughly centered on the Galactic centre, and extends to |b| ~35 deg in Galactic latitude and |l| ~15 deg in longitude. By combining the Planck data with observations from the WMAP we are able to determine the spectrum of this emission to high accuracy, unhindered by the large systematic biases present in previous analyses. The derived spectrum is consistent with power-law emission with a spectral index of -2.55 +/- 0.05, thus excluding free-free emission as the source and instead favouring hard-spectrum synchrotron radiation from an electron population with a spectrum (number density per energy) dN/dE ~ E^-2.1. At Galactic latitudes |b|<30 deg, the microwave haze morphology is consistent with that of the Fermi gamma-ray "haze" or "bubbles," indicating that we have a multi-wavelength view of a distinct component of our Galaxy. Given both the very hard spectrum and the extended nature of the emission, it is highly unlikely that the haze electrons result from supernova shocks in the Galactic disk. Instead, a new mechanism for cosmic-ray acceleration in the centre of our Galaxy is implied.

We present the final results from the XMM-Newton validation follow-up of new Planck galaxy cluster candidates. We observed 15 new candidates, detected with signal-to-noise ratios between 4.0 and 6.1 in the 15.5-month nominal Planck survey. The candidates were selected using ancillary data flags derived from the ROSAT All Sky Survey (RASS) and Digitized Sky Survey all-sky maps, with the aim of pushing into the low SZ flux, high-z regime and testing RASS flags as indicators of candidate reliability. 14 new clusters were detected by XMM, including 2 double systems. Redshifts lie in the range 0.2 to 0.9, with 6 clusters at z>0.5. Estimated M500 range from 2.5 10^14 to 8 10^14 Msun. We discuss our results in the context of the full XMM validation programme, in which 51 new clusters have been detected. This includes 4 double and 2 triple systems, some of which are chance projections on the sky of clusters at different z. We find that association with a RASS-BSC source is a robust indicator of the reliability of a candidate, whereas association with a FSC source does not guarantee that the SZ candidate is a bona fide cluster. Nevertheless, most Planck clusters appear in RASS maps, with a significance greater than 2 sigma being a good indication that the candidate is a real cluster. The full sample gives a Planck sensitivity threshold of Y500 ~ 4 10^-4 arcmin^2, with indication for Malmquist bias in the YX-Y500 relation below this level. The corresponding mass threshold depends on z. Systems with M500 > 5 10^14 Msun at z > 0.5 are easily detectable with Planck. The newly-detected clusters follow the YX-Y500 relation derived from X-ray selected samples. Compared to X-ray selected clusters, the new SZ clusters have a lower X-ray luminosity on average for their mass. There is no indication of departure from standard self-similar evolution in the X-ray versus SZ scaling properties. (abridged)

A comparison is presented of Sunyaev-Zeldovich measurements for 11 galaxy clusters as obtained by Planck and by the ground-based interferometer, the Arcminute Microkelvin Imager. Assuming a universal spherically-symmetric Generalised Navarro, Frenk & White (GNFW) model for the cluster gas pressure profile, we jointly constrain the integrated Compton-Y parameter (Y_500) and the scale radius (theta_500) of each cluster. Our resulting constraints in the Y_500-theta_500 2D parameter space derived from the two instruments overlap significantly for eight of the clusters, although, overall, there is a tendency for AMI to find the Sunyaev-Zeldovich signal to be smaller in angular size and fainter than Planck. Significant discrepancies exist for the three remaining clusters in the sample, namely A1413, A1914, and the newly-discovered Planck cluster PLCKESZ G139.59+24.18. The robustness of the analysis of both the Planck and AMI data is demonstrated through the use of detailed simulations, which also discount confusion from residual point (radio) sources and from diffuse astrophysical foregrounds as possible explanations for the discrepancies found. For a subset of our cluster sample, we have investigated the dependence of our results on the assumed pressure profile by repeating the analysis adopting the best-fitting GNFW profile shape which best matches X-ray observations. Adopting the best-fitting profile shape from the X-ray data does not, in general, resolve the discrepancies found in this subset of five clusters. Though based on a small sample, our results suggest that the adopted GNFW model may not be sufficiently flexible to describe clusters universally.

We present further results from the ongoing XMM-Newton validation follow-up of Planck cluster candidates, detailing X-ray observations of eleven candidates detected at a signal-to-noise ratio of 4.5<S/N<5.3 in the same 10-month survey maps used in the construction of the Early SZ sample. The sample was selected in order to test internal SZ quality flags, and the pertinence of these flags is discussed in light of the validation results. Ten of the candidates are found to be bona fide clusters lying below the RASS flux limit. Redshift estimates are available for all confirmed systems via X-ray Fe-line spectroscopy. They lie in the redshift range 0.19<z<0.94, demonstrating Planck's capability to detect clusters up to high z. The X-ray properties of the new clusters appear to be similar to previous new detections by Planck at lower z and higher SZ flux: the majority are X-ray underluminous for their mass, estimated using Y_X as mass proxy, and many have a disturbed morphology. We find tentative indication for Malmquist bias in the Y_SZ-Y_X relation, with a turnover at Y_SZ \sim 4 e-4 arcmin^2. We present additional new optical redshift determinations with ENO and ESO telescopes of candidates previously confirmed with XMM-Newton. The X-ray and optical redshifts for a total of 20 clusters are found to be in excellent agreement. We also show that useful lower limits can be put on cluster redshifts using X-ray data only via the use of the Y_X vs. Y_SZ and X-ray flux F_X vs. Y_SZ relations.

The performance of the Planck instruments in space is enabled by their low operating temperatures, 20K for LFI and 0.1K for HFI, achieved through a combination of passive radiative cooling and three active mechanical coolers. The scientific requirement for very broad frequency coverage led to two detector technologies with widely different temperature and cooling needs. Active coolers could satisfy these needs; a helium cryostat, as used by previous cryogenic space missions (IRAS, COBE, ISO, Spitzer, AKARI), could not. Radiative cooling is provided by three V-groove radiators and a large telescope baffle. The active coolers are a hydrogen sorption cooler (<20K), a 4He Joule-Thomson cooler (4.7K), and a 3He-4He dilution cooler (1.4K and 0.1K). The flight system was at ambient temperature at launch and cooled in space to operating conditions. The HFI bolometer plate reached 93mK on 3 July 2009, 50 days after launch. The solar panel always faces the Sun, shadowing the rest of Planck, andoperates at a mean temperature of 384K. At the other end of the spacecraft, the telescope baffle operates at 42.3K and the telescope primary mirror operates at 35.9K. The temperatures of key parts of the instruments are stabilized by both active and passive methods. Temperature fluctuations are driven by changes in the distance from the Sun, sorption cooler cycling and fluctuations in gas-liquid flow, and fluctuations in cosmic ray flux on the dilution and bolometer plates. These fluctuations do not compromise the science data.

Planck's all sky surveys at 30-857 GHz provide an unprecedented opportunity to follow the radio spectra of a large sample of extragalactic sources to frequencies 2-20 times higher than allowed by past, large area, ground-based surveys. We combine the results of the Planck Early Release Compact Source Catalog (ERCSC) with quasi-simultaneous ground-based observations, as well as archival data, at frequencies below or overlapping Planck frequency bands, to validate the astrometry and photometry of the ERCSC radio sources and study the spectral features shown in this new frequency window opened by Planck. The ERCSC source positions and flux density scales are found to be consistent with the ground-based observations. We present and discuss the spectral energy distributions (SEDs) of a sample of "extreme" radio sources to illustrate the richness of the ERCSC for the study of extragalactic radio sources. Variability is found to play a role in the unusual spectral features of some of these sources.

Spectral energy distributions (SEDs) and radio continuum spectra are presented for a northern sample of 104 extragalactic radio sources, based on the Planck Early Release Compact Source Catalogue (ERCSC) and simultaneous multifrequency data. The nine Planck frequencies, from 30 to 857 GHz, are complemented by a set of simultaneous observations ranging from radio to gamma-rays. This is the first extensive frequency coverage in the radio and millimetre domains for an essentially complete sample of extragalactic radio sources, and it shows how the individual shocks, each in their own phase of development, shape the radio spectra as they move in the relativistic jet. The SEDs presented in this paper were fitted with second and third degree polynomials to estimate the frequencies of the synchrotron and inverse Compton (IC) peaks, and the spectral indices of low and high frequency radio data, including the Planck ERCSC data, were calculated. SED modelling methods are discussed, with an emphasis on proper, physical modelling of the synchrotron bump using multiplecomponents. Planck ERCSC data also suggest that the original accelerated electron energy spectrum could be much harder than commonly thought, with power-law index around 1.5 instead of the canonical 2.5. The implications of this are discussed for the acceleration mechanisms effective in blazar shock. Furthermore in many cases the Planck data indicate that gamma-ray emission must originate in the same shocks that produce the radio emission.

All-sky data from the Planck survey and the Meta-Catalogue of X-ray detected Clusters of galaxies (MCXC) are combined to investigate the relationship between the thermal Sunyaev-Zeldovich (SZ) signal and X-ray luminosity. The sample comprises ~ 1600 X-ray clusters with redshifts up to ~ 1 and spans a wide range in X-ray luminosity. The SZ signal is extracted for each object individually, and the statistical significance of the measurement is maximised by averaging the SZ signal in bins of X-ray luminosity, total mass, or redshift. The SZ signal is detected at very high significance over more than two decades in X-ray luminosity (10^43 erg/s < L_500 E(z)^-7/3 < 2 X 10^45 erg/s). The relation between intrinsic SZ signal and X-ray luminosity is investigated and the measured SZ signal is compared to values predicted from X-ray data. Planck measurements and X-ray based predictions are found to be in excellent agreement over the whole explored luminosity range. No significant deviation from standard evolution of the scaling relations is detected. For the first time the intrinsic scatter in the scaling relation between SZ signal and X-ray luminosity is measured and found to be consistent with the one in the luminosity -- mass relation from X-ray studies. There is no evidence of any deficit in SZ signal strength in Planck data relative to expectations from the X-ray properties of clusters, underlining the robustness and consistency of our overall view of intra-cluster medium properties.

(abridged) Planck has observed the entire sky from 30 GHz to 857GHz. The observed foreground emission contains contributions from different phases of the interstellar medium (ISM). We have separated the observed Galactic emission into the different gaseous components (atomic, molecular and ionised) in each of a number of Galactocentric rings. Templates are created for various Galactocentric radii using velocity information from atomic (neutral hydrogen) and molecular (12CO) observations. The ionised template is assumed to be traced by free-free emission as observed by WMAP, while 408 MHz emission is used to trace the synchrotron component. Gas emission not traced by the above templates, namely "ark gas", as evidenced using Planck data, is included as an additional template, the first time such a component has been used in this way. These templates are then correlated with each of the Planck frequency bands, as well as other ancillary data. The emission per column density of the gas templates allows us to create distinct spectral energy distributions (SEDs) per Galactocentric ring and in each of the gaseous tracers from 1.4 GHz to 25 THz (12\mu m). Apart from the thermal dust and free-free emission, we have probed the Galaxy for anomalous (e.g., spinning) dust as well as synchrotron emission. We show that anomalous dust emission is present in the atomic, molecular and dark gas phases throughout the Galactic disk. The derived dust propeties associated with the dark gas phase are derived but do not allow us to reveal the nature of this phase. For all environments, the anomalous emission is consistent with rotation from polycyclic aromatic hydrocarbons (PAHs) and, according to our simple model, accounts for $(25\pm5)%$ (statistical) of the total emission at 30 GHz.

We describe the processing of data from the Low Frequency Instrument (LFI) used in production of the Planck Early Release Compact Source Catalogue (ERCSC). In particular, we discuss the steps involved in reducing the data from telemetry packets to cleaned, calibrated, time-ordered data (TOD) and frequency maps. Data are continuously calibrated using the modulation of the temperature of the cosmic microwave background radiation induced by the motion of the spacecraft. Noise properties are estimated from TOD from which the sky signal has been removed using a generalized least square map-making algorithm. Measured 1/f noise knee-frequencies range from 100mHz at 30GHz to a few tens of mHz at 70GHz. 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 approx -10dB level using Jupiter transits, which are also used for geometrical calibration of the focal plane.

Anomalous microwave emission (AME) has been observed by numerous experiments in the frequency range ~10-60 GHz. Using Planck maps and multi-frequency ancillary data, we have constructed spectra for two known AME regions: the Perseus and Rho Ophiuchi molecular clouds. The spectra are well fitted by a combination of free-free radiation, cosmic microwave background, thermal dust, and electric dipole radiation from small spinning dust grains. The spinning dust spectra are the most precisely measured to date, and show the high frequency side clearly for the first time. The spectra have a peak in the range 20-40 GHz and are detected at high significances of 17.1 sigma for Perseus and 8.4 sigma for Rho Ophiuchi. In Perseus, spinning dust in the dense molecular gas can account for most of the AME; the low density atomic gas appears to play a minor role. In Rho Ophiuchi, the ~30 GHz peak is dominated by dense molecular gas, but there is an indication of an extended tail at frequencies 50-100 GHz, which can be accounted for by irradiated low density atomic gas. The dust parameters are consistent with those derived from other measurements. We have also searched the Planck map at 28.5 GHz for candidate AME regions, by subtracting a simple model of the synchrotron, free-free, and thermal dust. We present spectra for two of the candidates; S140 and S235 are bright HII regions that show evidence for AME, and are well fitted by spinning dust models.

This paper presents the first results of comparison of Planck along with IRAS data with Green Bank Telescope 21-cm observations in 14 fields covering more than 800 square degrees at high Galactic latitude. Galactic dust emission for fields with average HI column density lower than 2 x 10^20 cm^-2 is well correlated with 21-cm emission. The residual emission in these fields, once the HI-correlated emission is removed, is consistent with the expected statistical properties of the cosmic infrared background fluctuations. Fields with larger column densities show significant excess dust emission compared to the HI column density. Regions of excess lie in organized structures that suggest the presence of hydrogen in molecular form, though they are not always correlated with CO emission. Dust emission from intermediate-velocity clouds is detected with high significance. Its spectral properties are consistent with, compared to the local ISM values, significantly hotter dust (T~20 K), lower sub-millimeter dust opacity, and a relative abundance of very small grains to large grains about four times higher. These results are compatible with expectations for clouds that are part of the Galactic fountain in which there is dust shattering and fragmentation. Correlated dust emission in HVCs is not detected; the average of the 99.9% confidence upper limits to the emissivity is 0.15 times the local ISM value at 857 and 3000 GHz, in accordance with gas phase evidence for lower metallicity and depletion in these clouds. Unexpected anti-correlated variations of the dust temperature and emission cross-section per H atom are identified in the local ISM and IVCs, a trend that continues into molecular environments. This suggests that dust growth through aggregation, seen in molecular clouds, is active much earlier in the cloud condensation and star formation processes.

We present the XMM-Newton follow-up for confirmation of Planck cluster candidates. Twenty-five candidates have been observed to date using snapshot (~10 ksec) exposures, ten as part of a pilot programme to sample a low range of signal-to-noise ratios (4<S/N<6), and a further 15 in a programme to observe a sample of S/N>5 candidates. The sensitivity and spatial resolution of XMM-Newton allows unambiguous discrimination between clusters and false candidates. The 4 false candidates have S/N <= 4.1. A total of 21 candidates are confirmed as extended X-ray sources. Seventeen are single clusters, the majority of which are found to have highly irregular and disturbed morphologies (about ~70%). The remaining four sources are multiple systems, including the unexpected discovery of a supercluster at z=0.45. For 20 sources we are able to derive a redshift estimate from the X-ray Fe K line (albeit of variable quality). The new clusters span the redshift range 0.09 <= z <= 0.54, with a median redshift of z~0.37. A first determination is made of their X-ray properties including the characteristic size, which is used to improve the estimate of the SZ Compton parameter, Y_SZ. The follow-up validation programme has helped to optimise the Planck candidate selection process. It has also provided a preview of the X-ray properties of these newly-discovered clusters, allowing comparison with their SZ properties, and to the X-ray and SZ properties of known clusters observed in the Planck survey. Our results suggest that Planck may have started to reveal a non-negligible population of massive dynamically perturbed objects that is under-represented in X-ray surveys. However, despite their particular properties, these new clusters appear to follow the Y_SZ-Y_X relation established for X-ray selected objects, where Y_X is the product of the gas mass and temperature.

We present precise Sunyaev-Zeldovich (SZ) effect measurements in the direction of 62 nearby galaxy clusters (z <0.5) detected at high signal-to-noise in the first Planck all-sky dataset. The sample spans approximately a decade in total mass, 10^14 < M_500 < 10^15, where M_500 is the mass corresponding to a total density contrast of 500. Combining these high quality Planck measurements with deep XMM-Newton X-ray data, we investigate the relations between D_A^2 Y_500, the integrated Compton parameter due to the SZ effect, and the X-ray-derived gas mass M_g,500, temperature T_X, luminosity L_X, SZ signal analogue Y_X,500 = M_g,500 * T_X, and total mass M_500. After correction for the effect of selection bias on the scaling relations, we find results that are in excellent agreement with both X-ray predictions and recently-published ground-based data derived from smaller samples. The present data yield an exceptionally robust, high-quality local reference, and illustrate Planck's unique capabilities for all-sky statistical studies of galaxy clusters.

We present a formalism for performance forecasting and optimization of future cosmic microwave background (CMB) experiments. We implement it in the context of nearly full sky, multifrequency, B-mode polarization observations, incorporating statistical uncertainties due to the CMB sky statistics, instrumental noise, as well as the presence of the foreground signals. We model the effects of a subtraction of these using a parametric maximum likelihood technique and optimize the instrumental configuration with predefined or arbitrary observational frequency channels, constraining either a total number of detectors or a focal plane area. We showcase the proposed formalism by applying it to two cases of experimental setups based on the CMBpol and COrE mission concepts looked at as dedicated B-mode experiments. We find that, if the models of the foregrounds available at the time of the optimization are sufficiently precise, the procedure can help to either improve the potential scientific outcome of the experiment by a factor of a few, while allowing one to avoid excessive hardware complexity, or simplify the instrument design without compromising its science goals. However, our analysis also shows that even if the available foreground models are not considered to be sufficiently reliable, the proposed procedure can guide a design of more robust experimental setups. While better suited to cope with a plausible, greater complexity of the foregrounds than that foreseen by the models, these setups could ensure science results close to the best achievable, should the models be found to be correct.

Planck allows unbiased mapping of Galactic sub-millimetre and millimetre emission from the most diffuse regions to the densest parts of molecular clouds. We present an early analysis of the Taurus molecular complex, on line-of-sight-averaged data and without component separation. The emission spectrum measured by Planck and IRAS can be fitted pixel by pixel using a single modified blackbody. Some systematic residuals are detected at 353 GHz and 143 GHz, with amplitudes around -7 % and +13 %, respectively, indicating that the measured spectra are likely more complex than a simple modified blackbody. Significant positive residuals are also detected in the molecular regions and in the 217 GHz and 100 GHz bands, mainly caused by to the contribution of the J=2-1 and J=1-0 12CO and 13CO emission lines. We derive maps of the dust temperature T, the dust spectral emissivity index beta, and the dust optical depth at 250 microns tau. The temperature map illustrates the cooling of the dust particles in thermal equilibrium with the incident radiation field, from 16-17 K in the diffuse regions to 13-14 K in the dense parts. The distribution of spectral indices is centred at 1.78, with a standard deviation of 0.08 and a systematic error of 0.07. We detect a significant T-beta anti-correlation. The dust optical depth map reveals the spatial distribution of the column density of the molecular complex from the densest molecular regions to the faint diffuse regions. We use near-infrared extinction and HI data at 21-cm to perform a quantitative analysis of the spatial variations of the measured dust optical depth at 250 microns per hydrogen atom tau/NH. We report an increase of tau/NH by a factor of about 2 between the atomic phase and the molecular phase, which has a strong impact on the equilibrium temperature of the dust particles.

Using Planck maps of six regions of low Galactic dust emission with a total area of about 140 square degrees, we determine the angular power spectra of cosmic infrared background (CIB) anisotropies from multipole l = 200 to l = 2000 at 217, 353, 545 and 857 GHz. We use 21-cm observations of HI as a tracer of thermal dust emission to reduce the already low level of Galactic dust emission and use the 143 GHz Planck maps in these fields to clean out cosmic microwave background anisotropies. Both of these cleaning processes are necessary to avoid significant contamination of the CIB signal. We measure correlated CIB structure across frequencies. As expected, the correlation decreases with increasing frequency separation, because the contribution of high-redshift galaxies to CIB anisotropies increases with wavelengths. We find no significant difference between the frequency spectrum of the CIB anisotropies and the CIB mean, with Delta I/I=15% from 217 to 857 GHz. In terms of clustering properties, the Planck data alone rule out the linear scale- and redshift-independent bias model. Non-linear corrections are significant. Consequently, we develop an alternative model that couples a dusty galaxy, parametric evolution model with a simple halo-model approach. It provides an excellent fit to the measured anisotropy angular power spectra and suggests that a different halo occupation distribution is required at each frequency, which is consistent with our expectation that each frequency is dominated by contributions from different redshifts. In our best-fit model, half of the anisotropy power at l=2000 comes from redshifts z<0.8 at 857 GHz and z<1.5 at 545 GHz, while about 90% come from redshifts z>2 at 353 and 217 GHz, respectively.

We consider future balloon-borne and ground-based suborbital experiments designed to search for inflationary gravitational waves, and investigate the impact of residual foregrounds that remain in the estimated cosmic microwave background maps. This is achieved by propagating foreground modelling uncertainties from the component separation, under the assumption of a spatially uniform foreground frequency scaling, through to the power spectrum estimates, and up to measurement of the tensor to scalar ratio in the parameter estimation step. We characterize the error covariance due to subtracted foregrounds, and find it to be subdominant compared to instrumental noise and sample variance in our simulated data analysis. We model the unsubtracted residual foreground contribution using a two-parameter power law and show that marginalization over these foreground parameters is effective in accounting for a bias due to excess foreground power at low $\ell$. We conclude that, at least in the suborbital experimental setups we have simulated, foreground errors may be modeled and propagated up to parameter estimation with only a slight degradation of the target sensitivity of these experiments derived neglecting the presence of the foregrounds.

(Abridged) A brief description of the methodology of construction, contents and usage of the Planck Early Release Compact Source Catalogue (ERCSC), including the Early Cold Cores (ECC) and the Early Sunyaev-Zeldovich (ESZ) cluster catalogue is provided. The catalogue is based on data that consist of mapping the entire sky once and 60% of the sky a second time by Planck, thereby comprising the first high sensitivity radio/submillimetre observations of the entire sky. A Monte-Carlo algorithm based on the injection and extraction of artificial sources into the Planck maps was implemented to select reliable sources among all extracted candidates such that the cumulative reliability of the catalogue is >=90%. The 10sigma photometric flux density limit of the catalogue at |b|>30 deg is 0.49, 1.0, 0.67, 0.5, 0.33, 0.28, 0.25, 0.47 and 0.82 Jy at each of the nine frequencies between 30 and 857 GHz. Sources which are up to a factor of ~2 fainter than this limit, and which are present in "clean" regions of the Galaxy where the sky background due to emission from the interstellar medium is low, are included in the ERCSC if they meet the high reliability criterion. The Planck ERCSC sources have known associations to stars with dust shells, stellar cores, radio galaxies, blazars, infrared luminous galaxies and Galactic interstellar medium features. A significant fraction of unclassified sources are also present in the catalogs. In addition, two early release catalogs that contain 915 cold molecular cloud core candidates and 189 SZ cluster candidates that have been generated using multi-frequency algorithms are presented. The entire source list, with more than 15000 unique sources, is ripe for follow-up characterisation with Herschel, ATCA, VLA, SOFIA, ALMA and other ground-based observing facilities.

We present the first all-sky sample of galaxy clusters detected blindly by the Planck satellite through the Sunyaev-Zeldovich (SZ) effect from its six highest frequencies. This early SZ (ESZ) sample is comprised of 189 candidates, which have a high signal-to-noise ratio ranging from 6 to 29. Its high reliability (purity above 95%) is further ensured by an extensive validation process based on Planck internal quality assessments and by external cross-identification and follow-up observations. Planck provides the first measured SZ signal for about 80% of the 169 previously-known ESZ clusters. Planck furthermore releases 30 new cluster candidates, amongst which 20 meet the ESZ signal-to-noise selection criterion. At the submission date, twelve of the 20 ESZ candidates were confirmed as new clusters, with eleven confirmed using XMM-Newton snapshot observations, most of them with disturbed morphologies and low luminosities. The ESZ clusters are mostly at moderate redshifts (86% with z below 0.3) and span more than a decade in mass, up to the rarest and most massive clusters with masses above 10^15 Msol.

The European Space Agency's Planck satellite was launched on 14 May 2009, and has been surveying the sky stably and continuously since 13 August 2009. Its performance is well in line with expectations, and it will continue to gather scientific data until the end of its cryogenic lifetime. We give an overview of the history of Planck in its first year of operations, and describe some of the key performance aspects of the satellite. This paper is part of a package submitted in conjunction with Planck's Early Release Compact Source Catalogue, the first data product based on Planck to be released publicly. The package describes the scientific performance of the Planck payload, and presents results on a variety of astrophysical topics related to the sources included in the Catalogue, as well as selected topics on diffuse emission.

The data reported in Planck's Early Release Compact Source Catalogue (ERCSC) are exploited to measure the number counts (dN/dS) of extragalactic radio sources at 30, 44, 70, 100, 143 and 217 GHz. Due to the full-sky nature of the catalogue, this measurement extends to the rarest and brightest sources in the sky. At lower frequencies (30, 44, and 70 GHz) our counts are in very good agreement with estimates based on WMAP data, being somewhat deeper at 30 and 70 GHz, and somewhat shallower at 44 GHz. Planck's source counts at 143 and 217 GHz join smoothly with the fainter ones provided by the SPT and ACT surveys over small fractions of the sky. An analysis of source spectra, exploiting Planck's uniquely broad spectral coverage, finds clear evidence of a steepening of the mean spectral index above about 70 GHz. This implies that, at these frequencies, the contamination of the CMB power spectrum by radio sources below the detection limit is significantly lower than previously estimated.

The all-sky coverage of the Planck Early Release Compact Source Catalogue (ERCSC) provides an unsurpassed survey of galaxies at submillimetre (submm) wavelengths, representing a major improvement in the numbers of galaxies detected, as well as the range of far-IR/submm wavelengths over which they have been observed. We here present the first results on the properties of nearby galaxies using these data. We match the ERCSC catalogue to IRAS-detected galaxies in the Imperial IRAS Faint Source Redshift Catalogue (IIFSCz), so that we can measure the spectral energy distributions (SEDs) of these objects from 60 to 850 microns. This produces a list of 1717 galaxies with reliable associations between Planck and IRAS, from which we select a subset of 468 for SED studies, namely those with strong detections in the three highest frequency Planck bands and no evidence of cirrus contamination. The SEDs are fitted using parametric dust models to determine the range of dust temperatures and emissivities. We find evidence for colder dust than has previously been found in external galaxies, with T<20K. Such cold temperatures are found using both the standard single temperature dust model with variable emissivity beta, or a two dust temperature model with beta fixed at 2. We also compare our results to studies of distant submm galaxies (SMGs) which have been claimed to contain cooler dust than their local counterparts. We find that including our sample of 468 galaxies significantly reduces the distinction between the two populations. Fits to SEDs of selected objects using more sophisticated templates derived from radiative transfer models confirm the presence of the colder dust found through parameteric fitting. We thus conclude that cold (T<20K) dust is a significant and largely unexplored component of many nearby galaxies.

We present the Sunyaev-Zeldovich (SZ) signal-to-richness scaling relation (Y500-N200) for the MaxBCG cluster catalogue. Employing a multi-frequency matched filter on the Planck sky maps, we measure the SZ signal for each cluster by adapting the filter according to weak-lensing calibrated mass-richness relations (N200-M500). We bin our individual measurements and detect the SZ signal down to the lowest richness systems (N200=10) with high significance, achieving a detection of the SZ signal in systems with mass as low as M500~5e13 Msolar. The observed Y500-N200 relation is well modeled by a power law over the full richness range. It has a lower normalisation at given N200 than predicted based on X-ray models and published mass-richness relations. An X-ray subsample, however, does conform to the predicted scaling, and model predictions do reproduce the relation between our measured bin-average SZ signal and measured bin-average X-ray luminosities. At fixed richness, we find an intrinsic dispersion in the Y500-N200 relation of 60% rising to of order 100% at low richness. Thanks to its all-sky coverage, Planck provides observations for more than 13,000 MaxBCG clusters and an unprecedented SZ/optical data set, extending the list of known cluster scaling laws to include SZ-optical properties. The data set offers essential clues for models of galaxy formation. Moreover, the lower normalisation of the SZ-mass relation implied by the observed SZ-richness scaling has important consequences for cluster physics and cosmological studies with SZ clusters.

We construct an all-sky map of the apparent temperature and optical depth of thermal dust emission using the Planck-HFI and IRAS data. The optical depth maps are correlated to tracers of the atomic and molecular gas. The correlation is linear in the lowest column density regions at high galactic latitudes. At high NH, the correlation is consistent with that of the lowest NH. In the intermediate NH range, we observe departure from linearity, with the dust optical depth in excess to the correlation. We attribute this excess emission to thermal emission by dust associated with a Dark-Gas phase, undetected in the available HI and CO measurements. We show the 2D spatial distribution of the Dark-Gas in the solar neighborhood and show that it extends around known molecular regions traced by CO. The average dust emissivity in the HI phase in the solar neighborhood follows roughly a power law distribution with beta = 1.8 all the way down to 3 mm, although the SED flattens slightly in the millimetre. The threshold for the existence of the Dark-Gas is found at NH = (8.0\pm 0.58) 10^{20} Hcm-2. Assuming the same dust emissivity at high frequencies for the dust in the atomic and molecular phases leads to an average XCO = (2.54\pm0.13) 10^{20} H2cm-2/(K km s-1). The mass of Dark-Gas is found to be 28% of the atomic gas and 118% of the CO emitting gas in the solar neighborhood. A possible origin for the Dark-Gas is the existence of a dark molecular phase, where H2 survives photodissociation but CO does not. The observed transition for the onset of this phase in the solar neighborhood (AV = 0.4 mag) appears consistent with recent theoretical predictions. We also discuss the possibility that up to half of the Dark-Gas could be in atomic form, due to optical depth effects in the HI measurements.

(abridged) We perform a detailed investigation of sources from the Cold Cores Catalogue of Planck Objects (C3PO). Our goal is to probe the reliability of the detections, validate the separation between warm and cold dust emission components, provide the first glimpse at the nature, internal morphology and physical characterictics of the Planck-detected sources. We focus on a sub-sample of ten sources from the C3PO list, selected to sample different environments, from high latitude cirrus to nearby (150pc) and remote (2kpc) molecular complexes. We present Planck surface brightness maps and derive the dust temperature, emissivity spectral index, and column densities of the fields. With the help of higher resolution Herschel and AKARI continuum observations and molecular line data, we investigate the morphology of the sources and the properties of the substructures at scales below the Planck beam size.

We present the statistical properties of the first version of the Cold Core Catalogue of Planck Objects (C3PO), in terms of their spatial distribution, temperature, distance, mass, and morphology. We also describe the statistics of the Early Cold Core Catalogue (ECC, delivered with the Early Release Compact Source Catalogue, ERCSC) that is the subset of the 915 most reliable detections of the complete catalogue. We have used the CoCoCoDeT algorithm to extract 10783 cold sources. Temperature and dust emission spectral index {\beta} values are derived using the fluxes in the IRAS 100 \mum band and the three highest frequency Planck bands. Temperature spans from 7K to 17K, and peaks around 13K. Data are not consistent with a constant value of {\beta} over the all temperature range. {\beta} ranges from 1.4 to 2.8 with a mean value around 2.1, and several possible scenarios are possible, including {\beta}(T) and the effect of multiple T components folded into the measurements. For one third of the objects the distances are obtained. Most of the detections are within 2 kpc in the Solar neighbourhood, but a few are at distances greater than 4 kpc. The cores are distributed from the deep Galactic plane, despite the confusion, to high latitudes (>30$^{\circle}$). The associated mass estimates range from 1 to $10^5$ solar masses. Using their physical properties these cold sources are shown to be cold clumps, defined as the intermediate cold sub-structures between clouds and cores. These cold clumps are not isolated but mostly organized in filaments associated with molecular clouds. The Cold Core Catalogue of Planck Objects (C3PO) is the first unbiased all-sky catalogue of cold objects. It gives an unprecedented statistical view to the properties of these potential pre-stellar clumps and offers a unique possibility for their classification in terms of their intrinsic properties and environment.

The integrated Spectral Energy Distributions of the Large and Small Magellanic Cloud appear significantly flatter than expected from dust models based on their FIR and radio emission. The origin of this millimetre excess is still unexplained, and is here investigated using the Planck data. The background CMB contribution is subtracted using an ILC method performed locally around the galaxies. The foreground emission from the Milky Way is subtracted. After subtraction, the emission of both galaxies correlates closely with the gas emission of the LMC and SMC. The millimetre excess in the LMC can be explained by CMB fluctuations, but a significant excess is still present in the SMC SED. The Planck and IRIS data at 100 micron are combined to produce thermal dust temperature and optical depth maps of the two galaxies. The LMC temperature map shows the presence of a warm inner arm already found with the Spitzer data, but also shows the existence of a previously unidentified cold outer arm. Several cold regions are found along this arm, some of which are associated with known molecular clouds. The average emissivity spectral index is found to be consistent with beta=1.5 and beta=1.2 below 500 microns for the LMC and SMC respectively, significantly flatter than the values observed in the Milky Way. Furthermore, there is evidence in the SMC for a further flattening of the SED in the sub-mm. The spatial distribution of the millimetre dust excess in the SMC follows the gas and thermal dust distribution. Different models are explored in order to fit the dust emission in the SMC. It is concluded that the millimetre excess is unlikely to be caused by very cold dust emission and that it could be due to a combination of spinning dust. emission and thermal dust emission by more amorphous dust grains than those present in our Galaxy.

EBEX is a NASA-funded balloon-borne experiment designed to measure the polarization of the cosmic microwave background (CMB). Observations will be made using 1432 transition edge sensor (TES) bolometric detectors read out with frequency multiplexed SQuIDs. EBEX will observe in three frequency bands centered at 150, 250, and 410 GHz, with 768, 384, and 280 detectors in each band, respectively. This broad frequency coverage is designed to provide valuable information about polarized foreground signals from dust. The polarized sky signals will be modulated with an achromatic half wave plate (AHWP) rotating on a superconducting magnetic bearing (SMB) and analyzed with a fixed wire grid polarizer. EBEX will observe a patch covering ~1% of the sky with 8' resolution, allowing for observation of the angular power spectrum from \ell = 20 to 1000. This will allow EBEX to search for both the primordial B-mode signal predicted by inflation and the anticipated lensing B-mode signal. Calculations to predict EBEX constraints on r using expected noise levels show that, for a likelihood centered around zero and with negligible foregrounds, 99% of the area falls below r = 0.035. This value increases by a factor of 1.6 after a process of foreground subtraction. This estimate does not include systematic uncertainties. An engineering flight was launched in June, 2009, from Ft. Sumner, NM, and the long duration science flight in Antarctica is planned for 2011. These proceedings describe the EBEX instrument and the North American engineering flight.

We investigate the performance of the parametric Maximum Likelihood component separation method in the context of the CMB B-mode signal detection and its characterization by small-scale CMB suborbital experiments. We consider high-resolution (FWHM=8') balloon-borne and ground-based observatories mapping low dust-contrast sky areas of 400 and 1000 square degrees, in three frequency channels, 150, 250, 410 GHz, and 90, 150, 220 GHz, with sensitivity of order 1 to 10 micro-K per beam-size pixel. These are chosen to be representative of some of the proposed, next-generation, bolometric experiments. We study the residual foreground contributions left in the recovered CMB maps in the pixel and harmonic domain and discuss their impact on a determination of the tensor-to-scalar ratio, r. In particular, we find that the residuals derived from the simulated data of the considered balloon-borne observatories are sufficiently low not to be relevant for the B-mode science. However, the ground-based observatories are in need of some external information to permit satisfactory cleaning. We find that if such information is indeed available in the latter case, both the ground-based and balloon-borne experiments can detect the values of r as low as ~0.04 at 95% confidence level. The contribution of the foreground residuals to these limits is found to be then subdominant and these are driven by the statistical uncertainty due to CMB, including E-to-B leakage, and noise. We emphasize that reaching such levels will require a sufficient control of the level of systematic effects present in the data.

This paper provides an overview of the Low Frequency Instrument (LFI) programme within the ESA Planck mission. The LFI instrument has been developed to produce high precision maps of the microwave sky at frequencies in the range 27-77 GHz, below the peak of the cosmic microwave background (CMB) radiation spectrum. The scientific goals are described, ranging from fundamental cosmology to Galactic and extragalactic astrophysics. The instrument design and development are outlined, together with the model philosophy and testing strategy. The instrument is presented in the context of the Planck mission. The LFI approach to ground and inflight calibration is described. We also describe the LFI ground segment. We present the results of a number of tests demonstrating the capability of the LFI data processing centre (DPC) to properly reduce and analyse LFI flight data, from telemetry information to calibrated and cleaned time ordered data, sky maps at each frequency (in temperature and polarization), component emission maps (CMB and diffuse foregrounds), catalogs for various classes of sources (the Early Release Compact Source Catalogue and the Final Compact Source Catalogue). The organization of the LFI consortium is briefly presented as well as the role of the core team in data analysis and scientific exploitation. All tests carried out on the LFI flight model demonstrate the excellent performance of the instrument and its various subunits. The data analysis pipeline has been tested and its main steps verified. In the first three months after launch, the commissioning, calibration, performance, and verification phases will be completed, after which Planck will begin its operational life, in which LFI will have an integral part.

How did the universe evolve? The fine angular scale (l>1000) temperature and polarization anisotropies in the CMB are a Rosetta stone for understanding the evolution of the universe. Through detailed measurements one may address everything from the physics of the birth of the universe to the history of star formation and the process by which galaxies formed. One may in addition track the evolution of the dark energy and discover the net neutrino mass. We are at the dawn of a new era in which hundreds of square degrees of sky can be mapped with arcminute resolution and sensitivities measured in microKelvin. Acquiring these data requires the use of special purpose telescopes such as the Atacama Cosmology Telescope (ACT), located in Chile, and the South Pole Telescope (SPT). These new telescopes are outfitted with a new generation of custom mm-wave kilo-pixel arrays. Additional instruments are in the planning stages.

We discuss an approach to the component separation of microwave, multi-frequency sky maps as those typically produced from Cosmic Microwave Background (CMB) Anisotropy data sets. The algorithm is based on the two step, parametric, likelihood-based technique recently elaborated on by Eriksen et al., (2006), where the foreground spectral parameters are estimated prior to the actual separation of the components. In contrast with the previous approaches, we accomplish the former task with help of an analytically-derived likelihood function for the spectral parameters, which, we show, yields estimates equal to the maximum likelihood values of the full multi-dimensional data problem. We then use these estimates to perform the second step via the standard, generalized-least-square-like procedure. We demonstrate that the proposed approach is equivalent to a direct maximization of the full data likelihood, which is recast in a computationally tractable form. We use the corresponding curvature matrices to characterize statistical properties of the recovered parameters. We incorporate in the formalism some of the essential features of the CMB data sets, such as inhomogeneous pixel domain noise, unknown map offsets as well as calibration errors and study their consequences for the separation. We find that the calibration is likely to have a dominant effect on the precision of the spectral parameter determination for a realistic CMB experiment. We apply the algorithm to simulated data and discuss the results. Our focus is on partial-sky, total-intensity and polarization, CMB experiments such as planned balloon-borne and ground-based efforts, however, the techniques presented here should be also applicable to the full-sky data as for instance, those produced by WMAP and anticipated from Planck.

The Planck satellite will map the full sky at nine frequencies from 30 to 857 GHz. The CMB intensity and polarization that are its prime targets are contaminated by foreground emission. The goal of this paper is to compare proposed methods for separating CMB from foregrounds based on their different spectral and spatial characteristics, and to separate the foregrounds into components of different physical origin. A component separation challenge has been organized, based on a set of realistically complex simulations of sky emission. Several methods including those based on internal template subtraction, maximum entropy method, parametric method, spatial and harmonic cross correlation methods, and independent component analysis have been tested. Different methods proved to be effective in cleaning the CMB maps from foreground contamination, in reconstructing maps of diffuse Galactic emissions, and in detecting point sources and thermal Sunyaev-Zeldovich signals. The power spectrum of the residuals is, on the largest scales, four orders of magnitude lower than that of the input Galaxy power spectrum at the foreground minimum. The CMB power spectrum was accurately recovered up to the sixth acoustic peak. The point source detection limit reaches 100 mJy, and about 2300 clusters are detected via the thermal SZ effect on two thirds of the sky. We have found that no single method performs best for all scientific objectives. We foresee that the final component separation pipeline for Planck will involve a combination of methods and iterations between processing steps targeted at different objectives such as diffuse component separation, spectral estimation and compact source extraction.

The Correlated Component Analysis (CCA) allows us to estimate how the different diffuse emissions mix in CMB experiments, exploiting also complementary information from other surveys. It is especially useful to deal with possible additional components. An application of CCA to WMAP maps assuming that only the canonical Galactic emissions are present, highlights the widespread presence of a spectrally flat "synchrotron" component, largely uncorrelated with the synchrotron template, suggesting that an additional foreground is indeed required. We have tested various spectral shapes for such component, namely a power law as expected if it is flat synchrotron, and two spectral shapes that may fit the spinning dust emission: a parabola in the logS - log(frequency) plane, and a grey body. Quality tests applied to the reconstructed CMB maps clearly disfavour two of the models. The CMB power spectra, estimated from CMB maps reconstructed exploiting the three surviving foreground models, are generally consistent with the WMAP ones, although at least one of them gives a significantly higher quadrupole moment than found by the WMAP team. Taking foreground modeling uncertainties into account, we find that the mean quadrupole amplitude for the three "good" models is less than 1 sigma below the expectation from the standard LambdaCDM model. Also the other reported deviations from model predictions are found not to be statistically significant, except for the excess power at l~40. We confirm the evidence for a marked North-South asymmetry in the large scale (l < 20) CMB anisotropies. We also present a first, albeit preliminary, all-sky map of the "anomalous" component.

In this work we study the relevance of the component separation technique based on the Independent Component Analysis (ICA) and investigate its performance in the context of a limited sky coverage observation and from the viewpoint of our ability to differentiate between cosmological models with different primordial B-mode content. We focus on the low Galactic emission sky patch, corresponding to the target of several operating and planned CMB experiments and which, in many respects, adequately represents a typical "clean" high latitude sky. We consider two fiducial observations, one operating at low (40, 90 GHz) and one at high (150, 350 GHz) frequencies and thus dominated by the synchrotron and thermal dust emission, respectively. We use a parallel version of the FASTICA code to explore a substantial parameter space including Gaussian pixel noise level, observed sky area and the amplitude of the foreground emission and employ large Monte Carlo simulations to quantify errors and biases pertinent to the reconstruction for different choices of the parameter values. We identify a large subspace of the parameter space for which the quality of the CMB reconstruction is excellent. For both the cosmological models, with and without the primordial gravitational waves, we find that FASTICA performs extremely well even in the cases when the B mode CMB signal is up to a few times weaker than the foreground contamination and the noise amplitude is comparable with the total CMB polarized emission. In addition we discuss limiting cases of the noise and foreground amplitudes, for which the ICA approach fails.

We present an application of the fast Independent Component Analysis (FastICA) to the WMAP 3yr data with the goal of extracting the CMB signal. We evaluate the confidence of our results by means of Monte Carlo simulations including CMB, foreground contaminations and instrumental noise specific of each WMAP frequency band. We perform a complete analysis involving all or a subset of the WMAP channels in order to select the optimal combination for CMB extraction, using the frequency scaling of the reconstructed component as a figure of merit. We found that the combination KQVW provides the best CMB frequency scaling, indicating that the low frequency foreground contamination in Q, V and W bands is better traced by the emission in the K band. The CMB angular power spectrum is recovered up to the degree scale, it is consistent within errors for all WMAP channel combination considered, and in close agreement with the WMAP 3yr results. We perform a statistical analysis of the recovered CMB pattern, and confirm the sky asymmetry reported in several previous works with independent techniques.