Planck constraints on the tensor-to-scalar ratio

We present constraints on the tensor-to-scalar ratio r using Planck data. We use the latest release of Planck maps, processed with the NPIPE code, which produces calibrated frequency maps in temperature and polarisation for all Planck channels from 30 GHz to 857 GHz using the same pipeline. We computed constraints on r using the BB angular power spectrum, and we also discuss constraints coming from the TT spectrum. Given Planck?s noise level, the TT spectrum gives constraints on r that are cosmic-variance limited (with ?r?=?0.093), but we show that the marginalised posterior peaks towards negative values of r at about the 1.2? level. We derived Planck constraints using the BB power spectrum at both large angular scales (the ?reionisation bump?) and intermediate angular scales (the ?recombination bump?) from ? = 2 to 150 and find a stronger constraint than that from TT, with ?r?=?0.069. The Planck BB spectrum shows no systematic bias and is compatible with zero, given both the statistical noise and the systematic uncertainties. The likelihood analysis using B modes yields the constraint r<?0.158 at 95% confidence using more than 50% of the sky. This upper limit tightens to r<?0.069 when Planck EE, BB, and EB power spectra are combined consistently, and it tightens further to r<?0.056 when the Planck TT power spectrum is included in the combination. Finally, combining Planck with BICEP2/Keck 2015 data yields an upper limit of r<?0.044.Planck is a project of the European Space Agency (ESA) with instruments provided by two scientific consortia funded by ESA member states and led by Principal Investigators from France and Italy, telescope reflectors provided through a collaboration between ESA and a scientific consortium led and funded by Denmark, and additional contributions from NASA (USA). Some of the results in this paper have been derived using the HEALPix package. This research used resources of the National Energy Research Scientific Computing Center (NERSC), a U.S. Department of Energy Office of Science User Facility operated under Contract No. DE-AC02-05CH11231. We gratefully acknowledge support from the CNRS/IN2P3 Computing Center for providing computing and data-processing resources needed for this work

Multiresolution internal template cleaning: an application to the Wilkinson Microwave Anisotropy Probe 7-yr polarization data

The cosmic microwave background (CMB) radiation data obtained by different experiments contain, besides the desired signal, a superposition of microwave sky contributions. Using a wavelet decomposition on the sphere, we present a fast and robust method to recover the CMB signal from microwave maps. We present an application to the Wilkinson Microwave Anisotropy Probe (WMAP) polarization data, which shows its good performance, particularly in very polluted regions of the sky. The applied wavelet has the advantages that it requires little computational time in its calculations, it is adapted to the HEALPIX pixelization scheme and it offers the possibility of multiresolution analysis. The decomposition is implemented as part of a fully internal template fitting method, minimizing the variance of the resulting map at each scale. Using a χ2 characterization of the noise, we find that the residuals of the cleaned maps are compatible with those expected from the instrumental noise. The maps are also comparable to those obtained from the WMAP team, but in our case we do not make use of external data sets. In addition, at low resolution, our cleaned maps present a lower level of noise. The E-mode power spectrum [fórmula] is computed at high and low resolutions, and a cross-power spectrum [fórmula] is also calculated from the foreground reduced maps of temperature given by WMAP and our cleaned maps of polarization at high resolution. These spectra are consistent with the power spectra supplied by the WMAP team. We detect the E-mode acoustic peak at ℓ∼ 400, as predicted by the standard ΛCDM model. The B-mode power spectrum [fórmula] is compatible with zero.We acknowledge partial financial support from the Spanish Ministerio de Ciencia e Innovación through projects AYA2010-21766-C03-01 and Consolider-Ingenio 2010 CSD2010-00064. RFC is grateful for financial support from Spanish CSIC for a JAE-predoc fellowship. PV is grateful for financial support from the Ramón y Cajal program. We are grateful for the computer resources, technical expertise and assistance provided by the Spanish Supercomputing Network (RES) node at the University of Cantabria.We are grateful for the use of LAMBDA and the assistance provided by Benjamin Gold by e-mail. The HEALPIX package was used throughout the data analysis (G´orski et al. 2005)

Planck 2018 results: II. Low Frequency Instrument data processing

We present a final description of the data-processing pipeline for the Planck Low Frequency Instrument (LFI), implemented for the 2018 data release. Several improvements have been made with respect to the previous release, especially in the calibration process and in the correction of instrumental features such as the effects of nonlinearity in the response of the analogue-to-digital converters. We provide a brief pedagogical introduction to the complete pipeline, as well as a detailed description of the important changes implemented. Self-consistency of the pipeline is demonstrated using dedicated simulations and null tests. We present the final version of the LFI full sky maps at 30, 44, and 70 GHz, both in temperature and polarization, together with a refined estimate of the solar dipole and a final assessment of the main LFI instrumental parameters.Acknowledgements. The Planck Collaboration acknowledges the support of: ESA; CNES and CNRS/INSU-IN2P3-INP (France); ASI, CNR, and INAF (Italy); NASA and DoE (USA); STFC and UKSA (UK); CSIC, MINECO, JA, and RES (Spain); Tekes, AoF, and CSC (Finland); DLR and MPG (Germany); CSA (Canada); DTU Space (Denmark); SER/SSO (Switzerland); RCN (Norway); SFI (Ireland); FCT/MCTES (Portugal); and ERC and PRACE (EU). A description of the Planck Collaboration and a list of its members, indicating which technical or scientific activities they have been involved in, can be found at http://www.cosmos.esa.int/web/planck/. The simulations for systematics assessment used the Extreme Science and Engineering Discovery Environment (XSEDE, Towns et al. 2014), supported by National Science Foundation grant number ACI-1548562, in particular the Comet Supercomputer at the San Diego Supercomputer Center through allocation AST160021: “Monte Carlo simulations for calibration uncertainty of the Planck mission”, PI A. Zonca, CoPI P. Meinhold

Planck 2018 results: XII. Galactic astrophysics using polarized dust emission

Observations of the submillimetre emission from Galactic dust, in both total intensity I and polarization, have received tremendous interest thanks to the Planck full-sky maps. In this paper we make use of such full-sky maps of dust polarized emission produced from the third public release of Planck data. As the basis for expanding on astrophysical studies of the polarized thermal emission from Galactic dust, we present full-sky maps of the dust polarization fraction p, polarization angle ?, and dispersion function of polarization angles ??. The joint distribution (one-point statistics) of p and NH confirms that the mean and maximum polarization fractions decrease with increasing NH. The uncertainty on the maximum observed polarization fraction, pmax = 22.0?1.4+3.5% at 353 GHz and 80? resolution, is dominated by the uncertainty on the Galactic emission zero level in total intensity, in particular towards diffuse lines of sight at high Galactic latitudes. Furthermore, the inverse behaviour between p and ?? found earlier is seen to be present at high latitudes. This follows the ?????p?1 relationship expected from models of the polarized sky (including numerical simulations of magnetohydrodynamical turbulence) that include effects from only the topology of the turbulent magnetic field, but otherwise have uniform alignment and dust properties. Thus, the statistical properties of p, ?, and ?? for the most part reflect the structure of the Galactic magnetic field. Nevertheless, we search for potential signatures of varying grain alignment and dust properties. First, we analyse the product map ???×?p, looking for residual trends. While the polarization fraction p decreases by a factor of 3?4 between NH?=?1020?cm?2 and NH?=?2?×?1022?cm?2, out of the Galactic plane, this product ???×?p only decreases by about 25%. Because ?? is independent of the grain alignment efficiency, this demonstrates that the systematic decrease in p with NH is determined mostly by the magnetic-field structure and not by a drop in grain alignment. This systematic trend is observed both in the diffuse interstellar medium (ISM) and in molecular clouds of the Gould Belt. Second, we look for a dependence of polarization properties on the dust temperature, as we would expect from the radiative alignment torque (RAT) theory. We find no systematic trend of ???×?p with the dust temperature Td, whether in the diffuse ISM or in the molecular clouds of the Gould Belt. In the diffuse ISM, lines of sight with high polarization fraction p and low polarization angle dispersion ?? tend, on the contrary, to have colder dust than lines of sight with low p and high ??. We also compare the Planck thermal dust polarization with starlight polarization data in the visible at high Galactic latitudes. The agreement in polarization angles is remarkable, and is consistent with what we expect from the noise and the observed dispersion of polarization angles in the visible on the scale of the Planck beam. The two polarization emission-to-extinction ratios, RP/p and RS/V, which primarily characterize dust optical properties, have only a weak dependence on the column density, and converge towards the values previously determined for translucent lines of sight. We also determine an upper limit for the polarization fraction in extinction, pV/E(B???V), of 13% at high Galactic latitude, compatible with the polarization fraction p???20% observed at 353 GHz. Taken together, these results provide strong constraints for models of Galactic dust in diffuse gas

Planck 2018 results: IV. Diffuse component separation

We present full-sky maps of the cosmic microwave background (CMB) and polarized synchrotron and thermal dust emission, derived from the third set of Planck frequency maps. These products have significantly lower contamination from instrumental systematic effects than previous versions. The methodologies used to derive these maps follow closely those described in earlier papers, adopting four methods (Commander, NILC, SEVEM, and SMICA) to extract the CMB component, as well as three methods (Commander, GNILC, and SMICA) to extract astrophysical components. Our revised CMB temperature maps agree with corresponding products in the Planck 2015 delivery, whereas the polarization maps exhibit significantly lower large-scale power, reflecting the improved data processing described in companion papers; however, the noise properties of the resulting data products are complicated, and the best available end-to-end simulations exhibit relative biases with respect to the data at the few percent level. Using these maps, we are for the first time able to fit the spectral index of thermal dust independently over 3° regions. We derive a conservative estimate of the mean spectral index of polarized thermal dust emission of ?d?=?1.55? ± ?0.05, where the uncertainty marginalizes both over all known systematic uncertainties and different estimation techniques. For polarized synchrotron emission, we find a mean spectral index of ?s?=??3.1? ± ?0.1, consistent with previously reported measurements. We note that the current data processing does not allow for construction of unbiased single-bolometer maps, and this limits our ability to extract CO emission and correlated components. The foreground results for intensity derived in this paper therefore do not supersede corresponding Planck 2015 products. For polarization the new results supersede the corresponding 2015 products in all respects.The Planck Collaboration acknowledges the support of: ESA; CNES, and CNRS/INSU-IN2P3-INP (France); ASI, CNR, and INAF (Italy); NASA and DoE (USA); STFC and UKSA (UK); CSIC, MINECO, JA, and RES (Spain); Tekes, AoF, and CSC (Finland); DLR and MPG (Germany); CSA (Canada); DTU Space (Denmark); SER/SSO (Switzerland); RCN (Norway); SFI (Ireland); FCT/MCTES (Portugal); ERC and PRACE (EU). A description of the Planck Collaboration and a list of its members, indicating which technical or scientific activities they have been involved in, can be found at http: //www.cosmos.esa.int/web/planck/planck-collaboration. This work has received funding from the European Union’s Horizon 2020 research and innovation programme under grant agreement numbers 687312, 776282 and 772253

Planck intermediate results III.The relation between galaxy cluster mass and Sunyaev-Zeldovich signal

We examine the relation between the galaxy cluster mass M and Sunyaev-Zeldovich (SZ) effect signal 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 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 R500 are on average ~20 percent larger than the corresponding weak lensing masses, which is contrary to 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, that the mass discrepancy and difference in mass concentration are 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

Planck intermediate results VI. The dynamical structure of PLCKG214.6+37.0, a Planck discovered triple system of galaxy clusters

The survey of galaxy clusters performed by Planck through the Sunyaev-Zeldovich effect has already discovered many interesting objects, thanks to its full sky coverage. One of the SZcandidates detected in the early months of the mission near to the signal-to-noise threshold, PLCKG214.6+37.0, was later revealed by XMM-Newton to be a triple system of galaxy clusters. We present the results from a deep XMM-Newton re-observation of PLCKG214.6+37.0, part of a multi-wavelength programme to investigate Planck discovered superclusters. The characterisation of the physical properties of the three components has allowed us to build a template model to extract the total SZ signal of this system with Planck data. We have partly reconciled the discrepancy between the expected SZ signal derived from X-rays and the observed one, which are now consistent within 1.2?. We measured the redshift of the three components with the iron lines in the X-ray spectrum, and confirm that the three clumps are likely part of the same supercluster structure. The analysis of the dynamical state of the three components, as well as the absence of detectable excess X-ray emission, suggests that we are witnessing the formation of a massive cluster at an early phase of interaction

Planck intermediate results XI. The gas content of dark matter halos: The sunyaev-zeldovich-stellar mass relation for locally brightest galaxies

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 M? ~ 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 × 1013 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

Planck intermediate results VII. Statistical properties of infrared and radio extragalactic sources from the Planck Early Release Compact Source Catalogue at frequencies between 100 and 857 GHz

We make use of the Planck all-sky survey to derive number counts and spectral indices of extragalactic sources ? infrared and radio sources ? from the Planck Early Release Compact Source Catalogue (ERCSC) at 100 to 857 GHz (3mm to 350 ?m). Three zones (deep, medium and shallow) of approximately homogeneous coverage are used to permit a clean and controlled correction for incompleteness, which was explicitly not done for the ERCSC, as it was aimed at providing lists of sources to be followed up. Our sample, prior to the 80% completeness cut, contains between 217 sources at 100 GHz and 1058 sources at 857 GHz over about 12800 to 16550 deg2 (31 to 40% of the sky). After the 80% completeness cut, between 122 and 452 and sources remain, with flux densities above 0.3 and 1.9 Jy at 100 and 857 GHz. The sample so defined can be used for statistical analysis. Using the multi-frequency coverage of the Planck High Frequency Instrument, all the sources have been classified as either dust-dominated (infrared galaxies) or synchrotron-dominated (radio galaxies) on the basis of their spectral energy distributions (SED). Our sample is thus complete, flux-limited and color-selected to differentiate between the two populations. We find an approximately equal number of synchrotron and dusty sources between 217 and 353 GHz; at 353 GHz or higher (or 217 GHz and lower) frequencies, the number is dominated by dusty (synchrotron) sources, as expected. For most of the sources, the spectral indices are also derived. We provide for the first time counts of bright sources from 353 to 857 GHz and the contributions from dusty and synchrotron sources at all HFI frequencies in the key spectral range where these spectra are crossing. The observed counts are in the Euclidean regime. The number counts are compared to previously published data (from earlier Planck results, Herschel, BLAST, SCUBA, LABOCA, SPT, and ACT) and models taking into account both radio or infrared galaxies, and covering a large range of flux densities. We derive the multi-frequency Euclidean level ? the plateau in the normalised differential counts at high flux-density ? and compare it to WMAP, Spitzer and IRAS results. The submillimetre number counts are not well reproduced by current evolution models of dusty galaxies, whereas the millimetre part appears reasonably well fitted by the most recent model for synchrotron-dominated sources. Finally we provide estimates of the local luminosity density of dusty galaxies, providing the first such measurements at 545 and 857 GHz

Planck intermediate results IV. The XMM-Newton validation programme for new Planck galaxy clusters

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. Fourteen new clusters were detected by XMM-Newton, ten single clusters and two double systems. Redshifts from X-ray spectroscopy lie in the range 0.2 to 0.9, with six clusters at z > 0.5. Estimated masses (M500) range from 2.5 × 1014 to 8 × 1014 M?. We discuss our results in the context of the full XMM-Newton validation programme, in which 51 new clusters have been detected. This includes four double and two triple systems, some of which are chance projections on the sky of clusters at different redshifts. We find thatassociation with a source from the RASS-Bright Source Catalogue is a robust indicator of the reliability of a candidate, whereas association with a source from the RASS-Faint Source Catalogue 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? being a good indication that the candidate is a real cluster. Candidate validation from association with SDSS galaxy overdensity at z > 0.5 is also discussed. The full sample gives a Planck sensitivity threshold of Y500 ~ 4 × 10-4 arcmin2, with indication for Malmquist bias in the YX?Y500 relation below this threshold. The corresponding mass threshold depends on redshift. Systems with M500 > 5 × 1014 M? 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. In particular, there is no significant evolution of the YX/Y500 ratio