Fast Pixel Space Convolution for CMB Surveys with Asymmetric Beams and Complex Scan Strategies: FEBeCoP

Precise measurement of the angular power spectrum of the Cosmic Microwave Background (CMB) temperature and polarization anisotropy can tightly constrain many cosmological models and parameters. However, accurate measurements can only be realized in practice provided all major systematic effects have been taken into account. Beam asymmetry, coupled with the scan strategy, is a major source of systematic error in scanning CMB experiments such as Planck, the focus of our current interest. We envision Monte Carlo methods to rigorously study and account for the systematic effect of beams in CMB analysis. Toward that goal, we have developed a fast pixel space convolution method that can simulate sky maps observed by a scanning instrument, taking into account real beam shapes and scan strategy. The essence is to pre-compute the "effective beams" using a computer code, "Fast Effective Beam Convolution in Pixel space" (FEBeCoP), that we have developed for the Planck mission. The code computes effective beams given the focal plane beam characteristics of the Planck instrument and the full history of actual satellite pointing, and performs very fast convolution of sky signals using the effective beams. In this paper, we describe the algorithm and the computational scheme that has been implemented. We also outline a few applications of the effective beams in the precision analysis of Planck data, for characterizing the CMB anisotropy and for detecting and measuring properties of point sources.Comment: 26 pages, 15 figures. New subsection on beam/PSF statistics, new and better figures, more explicit algebra for polarized beams, added explanatory text at many places following referees comments [Accepted for publication in ApJS

Estimating the tensor-to-scalar ratio and the effect of residual foreground contamination

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.Comment: 19 pages, 12 figures, accepted for publication in JCA

Planck 2013 results. XXVIII. The Planck Catalogue of Compact Sources

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 \u3e80% 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

Planck 2013 results. I. Overview of products and scientific results

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

Planck 2013 results. XXII. Constraints on inflation

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_: 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 V00 < 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=dln k = 0:0134 _ 0:0090. We verify these conclusions through a numerical analysis, which makes no slowroll 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 __2 e_ _ 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 _2 e_ by approximately 4 as a result of slightly lowering the theoretical prediction for the ` <_ 40 multipoles relative to the higher multipoles. Nonetheless, the data are consistent with adiabatic initial conditions

Planck intermediate results. XXII. Frequency dependence of thermal emission from Galactic dust in intensity and polarization

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

Planck intermediate results. XXI. Comparison of polarized thermal emission from Galactic dust at 353 GHz with interstellar polarization in the visible

The Planck survey provides unprecedented full-sky coverage of the submillimetre polarized emission from Galactic dust. In addition to the information on the direction of the Galactic magnetic field, this also brings new constraints on the properties of dust. The dust grains that emit the radiation seen by Planck in the submillimetre also extinguish and polarize starlight in the visible. Comparison of the polarization of the emission and of the interstellar polarization on selected lines of sight probed by stars provides unique new diagnostics of the emission and light scattering properties of dust, and therefore of the important dust model parameters, composition, size, and shape. Using ancillary catalogues of interstellar polarization and extinction of starlight, we obtain the degree of polarization, p(V), and the optical depth in the V band to the star, tau(V). Toward these stars we measure the submillimetre polarized intensity, P-S, and total intensity, I-S,I- in the Planck 353 GHz channel. We compare the column density measure in the visible, E(B - V), with that inferred from the Planck product map of the submillimetre dust optical depth and compare the polarization direction (position angle) in the visible with that in the submillimetre. For those lines of sight through the di ff use interstellar medium with comparable values of the estimated column density and polarization directions close to orthogonal, we correlate properties in the submillimetre and visible to find two ratios, R-S/V = (P-S/I-S) = (p(V)/tau(V)) and R-P/p = P-S/p(V), the latter focusing directly on the polarization properties of the aligned grain population alone. We find R-S/V = 4.2, with statistical and systematic uncertainties 0.2 and 0.3, respectively, and R-P/p = 5.4 MJy sr(-1), with uncertainties 0.2 and 0.3 MJy sr(-1), respectively. Our estimate of R-S/V is compatible with predictions based on a range of polarizing dust models that have been developed for the di ff use interstellar medium. This estimate provides new empirical validation of many of the common underlying assumptions of the models, but is not yet very discriminating among them. However, our estimate of R-P/p is not compatible with predictions, which are too low by a factor of about 2.5. This more discriminating diagnostic, R-P/p, indicates that changes to the optical properties in the models of the aligned grain population are required. These new diagnostics, together with the spectral dependence in the submillimetre from Planck, will be important for constraining and understanding the full complexity of the grain models, and for interpreting the Planck thermal dust polarization and refinement of the separation of this contamination of the cosmic microwave background.Peer reviewe

Planck intermediate results. XX. Comparison of polarized thermal emission from Galactic dust with simulations of MHD turbulence

Polarized emission observed by Planck HFI at 353GHz towards a sample of nearby fields is presented, focusing on the statistics of polarization fractions p and angles psi. The polarization fractions and column densities in these nearby fields are representative of the range of values obtained over the whole sky. We find that: (i) the largest polarization fractions are reached in the most diffuse fields; (ii) the maximum polarization fraction p(max) decreases with column density N-H in the more opaque fields with N-H > 10(21) cm(-2); and (iii) the polarization fraction along a given line of sight is correlated with the local spatial coherence of the polarization angle. These observations are compared to polarized emission maps computed in simulations of anisotropic magnetohydrodynamical turbulence in which we assume a uniform intrinsic polarization fraction of the dust grains. We find that an estimate of this parameter may be recovered from the maximum polarization fraction p(max) in diffuse regions where the magnetic field is ordered on large scales and perpendicular to the line of sight. This emphasizes the impact of anisotropies of the magnetic field on the emerging polarization signal. The decrease of the maximum polarization fraction with column density in nearby molecular clouds is well reproduced in the simulations, indicating that it is essentially due to the turbulent structure of the magnetic field: an accumulation of variously polarized structures along the line of sight leads to such an anti-correlation. In the simulations, polarization fractions are also found to anti-correlate with the angle dispersion function S. However, the dispersion of the polarization angle for a given polarization fraction is found to be larger in the simulations than in the observations, suggesting a shortcoming in the physical content of these numerical models. In summary, we find that the turbulent structure of the magnetic field is able to reproduce the main statistical properties of the dust polarization as observed in a variety of nearby clouds, dense cores excluded, and that the large-scale field orientation with respect to the line of sight plays a major role in the quantitative analysis of these statistical properties.Peer reviewe

Planck 2013 results. XIII. Galactic CO 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

Planck 2013 results : XXXII. The updated Planck catalogue of Sunyaev-Zeldovich sources

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 (similar to 80.6%) are spectroscopic, and associated mass estimates derived from the Y-z 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.Peer reviewe