Search

We present foreground-reduced cosmic microwave background (CMB) maps derived from the full Planck data set in both temperature and polarization. Compared to the corresponding Planck 2013 temperature sky maps, the total data volume is larger by a factor of 3.2 for frequencies between 30 and 70 GHz, and by 1.9 for frequencies between 100 and 857 GHz. In addition, systematic errors in the forms of temperature-to-polarization leakage, analogue-to-digital conversion uncertainties, and very long time constant errors have been dramatically reduced, to the extent that the cosmological polarization signal may now be robustly recovered on angular scales â"" μ 40. On the very largest scales, instrumental systematic residuals are still non-negligible compared to the expected cosmological signal, and modes with â"" < 20 are accordingly suppressed in the current polarization maps by high-pass filtering. As in 2013, four different CMB component separation algorithms are applied to these observations, providing a measure of stability with respect to algorithmic and modelling choices. The resulting polarization maps have rms instrumental noise ranging between 0.21 and 0.27μK averaged over 55′ pixels, and between 4.5 and 6.1μK averaged over \hbox{$3\parcm4$} pixels. The cosmological parameters derived from the analysis of temperature power spectra are in agreement at the 1σ level with the Planck 2015 likelihood. Unresolved mismatches between the noise properties of the data and simulations prevent a satisfactory description of the higher-order statistical properties of the polarization maps. Thus, the primary applications of these polarization maps are those that do not require massive simulations for accurate estimation of uncertainties, for instance estimation of cross-spectra and cross-correlations, or stacking analyses. However, the amplitude of primordial non-Gaussianity is consistent with zero within 2σ for all local, equilateral, and orthogonal configurations of the bispectrum, including for polarization E-modes. Moreover, excellent agreement is found regarding the lensing B-mode power spectrum, both internally among the various component separation codes and with the best-fit Planck 2015 Λ cold dark matter model.

The Planck full mission cosmic microwave background (CMB) temperature and E-mode polarization maps are analysed to obtain constraints on primordial non-Gaussianity (NG). Using three classes of optimal bispectrum estimators-separable template-fitting (KSW), binned, and modal-we obtain consistent values for the primordial local, equilateral, and orthogonal bispectrum amplitudes, quoting as our final result from temperature alone localNL = 2.5 ± 5.7, equilNL=-16 ± 70, fNLlocal=2.5±5.7, fNLequil=-16±70, and orthoNL =-34 ± 32fNLortho=-34±33 (68% CL, statistical). Combining temperature and polarization data we obtain fNLlocal=0.8±5.0, fNLequil=-4±43, and fNLortho=-26±21localNL = 0.8 ± 5.0, equilNL=-4 ± 43, and orthoNL =-26 ± 21 (68% CL, statistical). The results are based on comprehensive cross-validation of these estimators on Gaussian and non-Gaussian simulations, are stable across component separation techniques, pass an extensive suite of tests, and are consistent with estimators based on measuring the Minkowski functionals of the CMB. The effect of time-domain de-glitching systematics on the bispectrum is negligible. In spite of these test outcomes we conservatively label the results including polarization data as preliminary, owing to a known mismatch of the noise model in simulations and the data. Beyond estimates of individual shape amplitudes, we present model-independent, three-dimensional reconstructions of the Planck CMB bispectrum and derive constraints on early universe scenarios that generate primordial NG, including general single-field models of inflation, axion inflation, initial state modifications, models producing parity-violating tensor bispectra, and directionally dependent vector models. We present a wide survey of scale-dependent feature and resonance models, accounting for the "look elsewhere" effect in estimating the statistical significance of features. We also look for isocurvature NG, and find no signal, but we obtain constraints that improve significantly with the inclusion of polarization. The primordial trispectrum amplitude in the local model is constrained to be glocalNL = (-0.9 ± 7.7) X 104(68% CL statistical), and we perform an analysis of trispectrum shapes beyond the local case. The global picture that emerges is one of consistency with the premises of the ΛCDM cosmology, namely that the structure we observe today was sourced by adiabatic, passive, Gaussian, and primordial seed perturbations.

We present the implications for cosmic inflation of the Planck measurements of the cosmic microwave background (CMB) anisotropies in both temperature and polarization based on the full Planck survey, which includes more than twice the integration time of the nominal survey used for the 2013 release papers. The Planck full mission temperature data and a first release of polarization data on large angular scales measure the spectral index of curvature perturbations to be ns = 0.968 ± 0.006 and tightly constrain its scale dependence to dns/dlnk =-0.003 ± 0.007 when combined with the Planck lensing likelihood. When the Planck high-â"" polarization data are included, the results are consistent and uncertainties are further reduced. The upper bound on the tensor-to-scalar ratio is r0.002< 0.11 (95% CL). This upper limit is consistent with the B-mode polarization constraint r< 0.12 (95% CL) obtained from a joint analysis of the BICEP2/Keck Array and Planck data. These results imply that V(φ) â φ2 and natural inflation are now disfavoured compared to models predicting a smaller tensor-to-scalar ratio, such as R2 inflation. We search for several physically motivated deviations from a simple power-law spectrum of curvature perturbations, including those motivated by a reconstruction of the inflaton potential not relying on the slow-roll approximation. We find that such models are not preferred, either according to a Bayesian model comparison or according to a frequentist simulation-based analysis. Three independent methods reconstructing the primordial power spectrum consistently recover a featureless and smooth \hbox{${\cal P}-{\cal{R}}(k)$} over the range of scales 0.008 Mpc-1 â‰2 k â‰2 0.1 Mpc-1. At large scales, each method finds deviations from a power law, connected to a deficit at multipoles â"" ≠20-40 in the temperature power spectrum, but at an uncompelling statistical significance owing to the large cosmic variance present at these multipoles. By combining power spectrum and non-Gaussianity bounds, we constrain models with generalized Lagrangians, including Galileon models and axion monodromy models. The Planck data are consistent with adiabatic primordial perturbations, and the estimated values for the parameters of the base Λ cold dark matter (ΛCDM) model are not significantly altered when more general initial conditions are admitted. In correlated mixed adiabatic and isocurvature models, the 95% CL upper bound for the non-adiabatic contribution to the observed CMB temperature variance is | αnon-adi | < 1.9%, 4.0%, and 2.9% for CDM, neutrino density, and neutrino velocity isocurvature modes, respectively. We have tested inflationary models producing an anisotropic modulation of the primordial curvature power spectrum findingthat the dipolar modulation in the CMB temperature field induced by a CDM isocurvature perturbation is not preferred at a statistically significant level. We also establish tight constraints on a possible quadrupolar modulation of the curvature perturbation. These results are consistent with the Planck 2013 analysis based on the nominal mission data and further constrain slow-roll single-field inflationary models, as expected from the increased precision of Planck data using the full set of observations.

We present a description of the pipeline used to calibrate the Planck Low Frequency Instrument (LFI) timelines into thermodynamic temperatures for the Planck 2015 data release, covering four years of uninterrupted operations. As in the 2013 data release, our calibrator is provided by the spin-synchronous modulation of the cosmic microwave background dipole, but we now use the orbital component, rather than adopting the Wilkinson Microwave Anisotropy Probe (WMAP) solar dipole. This allows our 2015 LFI analysis to provide an independent Solar dipole estimate, which is in excellent agreement with that of HFI and within 1σ (0.3% in amplitude) of the WMAP value. This 0.3% shift in the peak-to-peak dipole temperature from WMAP and a general overhaul of the iterative calibration code increases the overall level of the LFI maps by 0.45% (30 GHz), 0.64% (44 GHz), and 0.82% (70 GHz) in temperature with respect to the 2013 Planck data release, thus reducing the discrepancy with the power spectrum measured by WMAP. We estimate that the LFI calibration uncertainty is now at the level of 0.20% for the 70 GHz map, 0.26% for the 44 GHz map, and 0.35% for the 30 GHz map. We provide a detailed description of the impact of all the changes implemented in the calibration since the previous data release.

We present the 8th full focal plane simulation set (FFP8), deployed in support of the Planck 2015 results. FFP8 consists of 10 fiducial mission realizations reduced to 18 144 maps, together with the most massive suite of Monte Carlo realizations of instrument noise and CMB ever generated, comprising 104 mission realizations reduced to about 106 maps. The resulting maps incorporate the dominant instrumental, scanning, and data analysis effects, and the remaining subdominant effects will be included in future updates. Generated at a cost of some 25 million CPU-hours spread across multiple high-performance-computing (HPC) platforms, FFP8 is used to validate and verify analysis algorithms and their implementations, and to remove biases from and quantify uncertainties in the results of analyses of the real data.

The European Space Agency's Planck satellite, which is dedicated to studying the early Universe and its subsequent evolution, was launched on 14 May 2009. It scanned the microwave and submillimetre sky continuously between 12 August 2009 and 23 October 2013. In February 2015, ESA and the Planck Collaboration released the second set of cosmology products based ondata from the entire Planck mission, including both temperature and polarization, along with a set of scientific and technical papers and a web-based explanatory supplement. This paper gives an overview of the main characteristics of the data and the data products in the release, as well as the associated cosmological and astrophysical science results and papers. The data products include maps of the cosmic microwave background (CMB), the thermal Sunyaev-Zeldovich effect, diffuse foregrounds in temperature and polarization, catalogues of compact Galactic and extragalactic sources (including separate catalogues of Sunyaev-Zeldovich clusters and Galactic cold clumps), and extensive simulations of signals and noise used in assessing uncertainties and the performance of the analysis methods. The likelihood code used to assess cosmological models against the Planck data is described, along with a CMB lensing likelihood. Scientific results include cosmological parameters derived from CMB power spectra, gravitational lensing, and cluster counts, as well as constraints on inflation, non-Gaussianity, primordial magnetic fields, dark energy, and modified gravity, and new results on low-frequency Galactic foregrounds.

We present the all-sky Planck catalogue of Sunyaev-Zeldovich (SZ) sources detected from the 29 month full-mission data. The catalogue (PSZ2) is the largest SZ-selected sample of galaxy clusters yet produced and the deepest systematic all-sky surveyof galaxy clusters. It contains 1653 detections, of which 1203 are confirmed clusters with identified counterparts in external data sets, and is the first SZ-selected cluster survey containing >103 confirmed clusters. We present a detailed analysis of the survey selection function in terms of its completeness and statistical reliability, placing a lower limit of 83% on the purity. Using simulations, we find that the estimates of the SZ strength parameter Y5R500are robust to pressure-profile variation and beam systematics, but accurate conversion to Y500 requires the use of prior information on the cluster extent. We describe the multi-wavelength search for counterparts in ancillary data, which makes use of radio, microwave, infra-red, optical, and X-ray data sets, and which places emphasis on the robustness of the counterpart match. We discuss the physical properties of the new sample and identify a population of low-redshift X-ray under-luminous clusters revealed by SZ selection. These objects appear in optical and SZ surveys with consistent properties for their mass, but are almost absent from ROSAT X-ray selected samples.

This paper describes the mapmaking procedure applied to Planck Low Frequency Instrument (LFI) data. The mapmaking step takes as input the calibrated timelines and pointing information. The main products are sky maps of I, Q, and U Stokes components. For the first time, we present polarization maps at LFI frequencies. The mapmaking algorithm is based on a destriping technique, which is enhanced with a noise prior. The Galactic region is masked to reduce errors arising from bandpass mismatch and high signal gradients. We apply horn-uniform radiometer weights to reduce the effects of beam-shape mismatch. The algorithm is the same as used for the 2013 release, apart from small changes in parameter settings. We validate the procedure through simulations. Special emphasis is put on the control of systematics, which is particularly important for accurate polarization analysis. We also produce low-resolution versions of the maps and corresponding noise covariance matrices. These serve as input in later analysis steps and parameter estimation. The noise covariance matrices are validated through noise Monte Carlo simulations. The residual noise in the map products is characterized through analysis of half-ring maps, noise covariance matrices, and simulations.

This paper presents a study of the integrated Sachs-Wolfe (ISW) effect from the Planck 2015 temperature and polarization data release. This secondary cosmic microwave background (CMB) anisotropy caused by the large-scale time-evolving gravitational potential is probed from different perspectives. The CMB is cross-correlated with different large-scale structure (LSS) tracers: radio sources from the NVSS catalogue; galaxies from the optical SDSS and the infrared WISE surveys; and the Planck 2015 convergence lensing map. The joint cross-correlation of the CMB with the tracers yields a detection at 4σ where most of the signal-to-noise is due to the Planck lensing and the NVSS radio catalogue. In fact, the ISW effect is detected from the Planck data only at ≠3σ (through the ISW-lensing bispectrum), which is similar to the detection level achieved by combining the cross-correlation signal coming from all the galaxy catalogues mentioned above. We study the ability of the ISW effect to place constraints on the dark-energy parameters; in particular, we show that ΩΛ is detected at more than 3σ. This cross-correlation analysis is performed only with the Planck temperature data, since the polarization scales available in the 2015 release do not permit significant improvement of the CMB-LSS cross-correlation detectability. Nevertheless, the Planck polarization data are used to study the anomalously large ISW signal previously reported through the aperture photometry on stacked CMB features at the locations of known superclusters and supervoids, which is in conflict with ΛCDM expectations. We find that the current Planck polarization data do not exclude that this signal could be caused by the ISW effect. In addition, the stacking of the Planck lensing map on the locations of superstructures exhibits a positive cross-correlation with these large-scale structures. Finally, we have improved our previous reconstruction of the ISW temperature fluctuations by combining the information encoded in all the previously mentioned LSS tracers. In particular, we construct a map of the ISW secondary anisotropies and the corresponding uncertainties map, obtained from simulations. We also explore the reconstruction of the ISW anisotropies caused by the large-scale structure traced by the 2MASS Photometric Redshift Survey (2MPZ) by directly inverting the density field into the gravitational potential field.

We present cluster counts and corresponding cosmological constraints from the Planck full mission data set. Our catalogue consists of 439 clusters detected via their Sunyaev-Zeldovich (SZ) signal down to a signal-to-noise ratio of 6, and is more than a factor of 2 larger than the 2013 Planck cluster cosmology sample. The counts are consistent with those from 2013 and yield compatible constraints under the same modelling assumptions. Taking advantage of the larger catalogue, we extend our analysis to the two-dimensional distribution in redshift and signal-to-noise. We use mass estimates from two recent studies of gravitational lensing of background galaxies by Planck clusters to provide priors on the hydrostatic bias parameter, (1-b). In addition, we use lensing of cosmic microwave background (CMB) temperature fluctuations by Planck clusters as an independent constraint on this parameter. These various calibrations imply constraints on the present-day amplitude of matter fluctuations in varying degrees of tension with those from the Planck analysis of primary fluctuations in the CMB; for the lowest estimated values of (1-b) the tension is mild, only a little over one standard deviation, while it remains substantial (3.7σ) for the largest estimated value. We also examine constraints on extensions to the base flat ΛCDM model by combining the cluster and CMB constraints. The combination appears to favour non-minimal neutrino masses, but this possibility does little to relieve the overall tension because it simultaneously lowers the implied value of the Hubble parameter, thereby exacerbating the discrepancy with most current astrophysical estimates. Improving the precision of cluster mass calibrations from the current 10%-level to 1% would significantly strengthen these combined analyses and provide a stringent test of the base ΛCDM model.

We present an updated description of the Planck Low Frequency Instrument (LFI) data processing pipeline, associated with the 2015 data release. We point out the places where our results and methods have remained unchanged since the 2013 paper and we highlight the changes made for the 2015 release, describing the products (especially timelines) and the ways in which they were obtained. We demonstrate that the pipeline is self-consistent (principally based on simulations) and report all null tests. For the first time, we present LFI maps in Stokes Q and U polarization. We refer to other related papers where more detailed descriptions of the LFI data processing pipeline may be found if needed.

This paper presents the characterization of the in-flight beams, the beam window functions, and the associated uncertainties for the Planck Low Frequency Instrument (LFI). The structure of the paper is similar to that presented in the 2013 Planck release; the main differences concern the beam normalization and the delivery of the window functions to be used for polarization analysis. The in-flight assessment of the LFI main beams relies on measurements performed during observations of Jupiter. By stacking data from seven Jupiter transits, the main beam profiles are measured down to -25 dB at 30 and 44 GHz, and down to -30 dB at 70 GHz. It has been confirmed that the agreement between the simulated beams and the measured beams is better than 1% at each LFI frequency band (within the 20 dB contour from the peak, the rms values are 0.1% at 30 and 70 GHz; 0.2% at 44 GHz). Simulated polarized beams are used for the computation of the effective beam window functions. The error budget for the window functions is estimated from both main beam and sidelobe contributions, and accounts for the radiometer band shapes. The total uncertainties in the effective beam window functions are 0.7% and 1% at 30 and 44 GHz, respectively (at ℓ ≈ 600); and 0.5% at 70 GHz (at ℓ ≈ 1000).

The all-sky Planck survey in 9 frequency bands was used to search for emission from all 274 known Galactic supernova remnants. Of these, 16 were detected in at least two Planck frequencies. The radio-through-microwave spectral energy distributions were compiled to determine the mechanism for microwave emission. In only one case, IC 443, is there high-frequency emission clearly from dust associated with the supernova remnant. In all cases, the low-frequency emission is from synchrotron radiation. As predicted for a population of relativistic particles with energy distribution that extends continuously to high energies, a single power law is evident for many sources, including the Crab and PKS 1209-51/52. A decrease in flux density relative to the extrapolation of radio emission is evident in several sources. Their spectral energy distributions can be approximated as broken power laws, Sν ∝ ν-α, with the spectral index, α, increasing by 0.5-1 above a break frequency in the range 10-60 GHz. The break could be due to synchrotron losses.

We present all-sky modelling of the high resolution Planck, IRAS, and WISE infrared (IR) observations using the physical dust model presented by Draine & Li in 2007 (DL, ApJ, 657, 810). We study the performance and results of this model, and discuss implications for future dust modelling. The present work extends the DL dust modelling carried out on nearby galaxies using Herschel and Spitzer data to Galactic dust emission. We employ the DL dust model to generate maps of the dust mass surface density ΣMd, the dust optical extinction AV, and the starlight intensity heating the bulk of the dust, parametrized by Umin. The DL model reproduces the observed spectral energy distribution (SED) satisfactorily over most of the sky, with small deviations in the inner Galactic disk and in low ecliptic latitude areas, presumably due to zodiacal light contamination. In the Andromeda galaxy (M31), the present dust mass estimates agree remarkably well (within 10%) with DL estimates based on independent Spitzer and Herschel data. We compare the DL optical extinction AV for the diffuse interstellar medium (ISM) with optical estimates for approximately 2 × 105 quasi-stellar objects (QSOs) observed inthe Sloan Digital Sky Survey (SDSS). The DL AV estimates are larger than those determined towards QSOs by a factor of about 2, which depends on Umin. The DL fitting parameter Umin, effectively determined by the wavelength where the SED peaks, appears to trace variations in the far-IR opacity of the dust grains per unit AV, and not only in the starlight intensity. These results show that some of the physical assumptions of the DL model will need to be revised. To circumvent the model deficiency, we propose an empirical renormalization of the DL AV estimate, dependent of Umin, which compensates for the systematic differences found with QSO observations. This renormalization, made to match the AV estimates towards QSOs, also brings into agreement the DL AV estimates with those derived for molecular clouds from the near-IR colours of stars in the 2 micron all sky survey (2MASS). The DL model and the QSOs data are also used to compress the spectral information in the Planck and IRAS observations for the diffuse ISM to a family of 20 SEDs normalized per AV, parameterized by Umin, which may be used to test and empirically calibrate dust models. The family of SEDs and the maps generated with the DL model are made public in the Planck Legacy Archive.

The role of the magnetic field in the formation of the filamentary structuresobserved in the interstellar medium (ISM) is a debated topic. The Planckall-sky maps of linearly polarized emission from dust at 353GHz provide therequired combination of imaging and statistics to study the correlation betweenthe structures of the Galactic magnetic field and of interstellar matter, bothin the diffuse ISM and in molecular clouds. The data reveal structures, orridges, in the intensity map with counterparts in the Stokes Q and/or U maps.We focus on structures at intermediate and high Galactic latitudes with columndensity from $10^{20}$ to $10^{22}$ cm$^{-2}$. We measure the magnetic fieldorientation on the plane of the sky from the polarization data, and present analgorithm to estimate the orientation of the ridges from the dust intensitymap. We use analytical models to account for projection effects. Comparingpolarization angles on and off the structures, we estimate the mean ratiobetween the strengths of the turbulent and mean components of the magneticfield to be between 0.6 and 1.0, with a preferred value of 0.8. We find thatthe ridges are preferentially aligned with the magnetic field measured on thestructures. This trend becomes more striking for increasing polarizationfraction and decreasing column density. We interpret the increase of alignmentwith polarization fraction as a consequence of projections effects. Thedecrease of alignment for high column density is not due to a loss ofcorrelation between the structures and the geometry of the magnetic field. Inmolecular complexes, we observe structures perpendicular to the magnetic field,which cannot be accounted for by projection effects. We discuss our results inthe context of models and MHD simulations, which describe the formation ofstructures in the magnetized ISM.

Within ten nearby (d < 450 pc) Gould belt molecular clouds we evaluate statistically the relative orientation between the magnetic field projected on the plane of sky, inferred from the polarized thermal emission of Galactic dust observed by Planck at 353 GHz, and the gas column density structures, quantified by the gradient of the column density, NH. The selected regions, covering several degrees in size, are analysed at an effective angular resolution of 10′ FWHM, thus sampling physical scales from 0.4 to 40 pc in the nearest cloud. The column densities in the selected regions range from NH ≈ 1021 to 1023 cm-2, and hence they correspond to the bulk of the molecular clouds. The relative orientation is evaluated pixel by pixel and analysed in bins of column density using the novel statistical tool called "histogram of relative orientations". Throughout this study, we assume that the polarized emission observed by Planck at 353 GHz is representative of the projected morphology of the magnetic field in each region, i.e., we assume a constant dust grain alignment efficiency, independent of the local environment. Within most clouds we find that the relative orientation changes progressively with increasing NH, from mostly parallel or having no preferred orientation to mostly perpendicular. In simulations of magnetohydrodynamic turbulence in molecular clouds this trend in relative orientation is a signature of Alfvénic or sub-Alfvénic turbulence, implying that the magnetic field is significant for the gas dynamics at the scales probed by Planck. We compare the deduced magnetic field strength with estimates we obtain from other methods and discuss the implications of the Planck observations for the general picture of molecular cloud formation and evolution.

By looking at the kinetic Sunyaev-Zeldovich effect (kSZ) in Planck nominal mission data, we present a significant detection of baryons participating in large-scale bulk flows around central galaxies (CGs) at redshift z ≈ 0.1. We estimate the pairwise momentum of the kSZ temperature fluctuations at the positions of the Central Galaxy Catalogue (CGC) samples extracted from Sloan Digital Sky Survey (SDSS-DR7) data. For the foreground-cleaned SEVEM, SMICA, NILC, and COMMANDER maps, we find 1.8-2.5σ detections of the kSZ signal, which are consistent with the kSZ evidence found in individual Planck raw frequency maps, although lower than found in the WMAP-9yr W-band (3.3σ). We further reconstruct the peculiar velocity field from the CG density field, and compute for the first time the cross-correlation function between kSZ temperature fluctuations and estimates of CG radial peculiar velocities. This correlation function yields a 3.0-3.7σ detection of the peculiar motion of extended gas on Mpc scales in flows correlated up to distances of 80-100 h-1 Mpc. Both the pairwise momentum estimates and the kSZ temperature-velocity field correlation find evidence for kSZ signatures out to apertures of 8 arcmin and beyond, corresponding to a physical radius of >1 Mpc, more than twice the mean virial radius of halos. This is consistent with the predictions from hydrodynamical simulations that most of the baryons are outside the virialized halos. We fit a simple model, in which the temperature-velocity cross-correlation is proportional to the signal seen in a semi-analytic model built upon N-body simulations, and interpret the proportionality constant as an effective optical depth to Thomson scattering. We find τT = (1.4 ± 0.5) × 10-4; the simplest interpretation of this measurement is that much of the gas is in a diffuse phase, which contributes little signal to X-ray or thermal Sunyaev-Zeldovich observations.

We present the results of approximately three years of observations of Planck Sunyaev-Zeldovich (SZ) sources with telescopes at the Canary Islands observatories as part of the general optical follow-up programme undertaken by the Planck Collaboration. In total, 78 SZ sources are discussed. Deep-imaging observations were obtained for most of these sources; spectroscopic observations in either in long-slit or multi-object modes were obtained for many. We effectively used 37.5 clear nights. We found optical counterparts for 73 of the 78 candidates. This sample includes 53 spectroscopic redshift determinations, 20 of them obtained with a multi-object spectroscopic mode. The sample contains new redshifts for 27 Planck clusters that were not included in the first Planck SZ source catalogue (PSZ1).

The all-sky Planck survey in 9 frequency bands was used to search for emission from all 274 known Galactic supernova remnants. Of these, 16 were detected in at least two Planck frequencies. The radio-through-microwave spectral energy distributions were compiled to determine the mechanism for microwave emission. In only one case, IC 443, is there high-frequency emission clearly from dust associated with the supernova remnant. In all cases, the low-frequency emission is from synchrotron radiation. As predicted for a population of relativistic particles with energy distribution that extends continuously to high energies, a single power law is evident for many sources, including the Crab and PKS 1209-51/52. A decrease in flux density relative to the extrapolation of radio emission is evident in several sources. Their spectral energy distributions can be approximated as broken power laws, Sν ∝ ν-α, with the spectral index, α, increasing by 0.5-1 above a break frequency in the range 10-60 GHz. The break could be due to synchrotron losses.

The quest for a B-mode imprint from primordial gravity waves on the polarization of the cosmic microwave background (CMB) requires the characterization of foreground polarization from Galactic dust. We present a statistical study of the filamentary structure of the 353 GHz Planck Stokes maps at high Galactic latitude, relevant to the study of dust emission as a polarized foreground to the CMB. We filter the intensity and polarization maps to isolate filaments in the range of angular scales where the power asymmetry between E-modes and B-modes is observed. Using the Smoothed Hessian Major Axis Filament Finder (SMAFF), we identify 259 filaments at high Galactic latitude, with lengths larger or equal to 2° (corresponding to 3.5 pc in length for a typical distance of 100 pc). These filaments show a preferred orientation parallel to the magnetic field projected onto the plane of the sky, derived from their polarization angles. We present mean maps of the filaments in Stokes I, Q, U, E, and B, computed by stacking individual images rotated to align the orientations of the filaments. Combining the stacked images and the histogram of relative orientations, we estimate the mean polarization fraction of the filaments to be 11%. Furthermore, we show that the correlation between the filaments and the magnetic field orientations may account for the E and B asymmetry and the CℓTE/CℓEE ratio, reported in the power spectra analysis of the Planck353 GHz polarization maps. Future models of the dust foreground for CMB polarization studies will need to take into account the observed correlation between the dust polarization and the structure of interstellar matter.

The nearby Chamaeleon clouds have been observed in γ rays by the Fermi Large Area Telescope (LAT) and in thermal dust emission by Planck and IRAS. Cosmic rays and large dust grains, if smoothly mixed with gas, can jointly serve with the Hi and 12CO radio data to (i) map the hydrogen column densities, NH, in the different gas phases, in particular at the dark neutral medium (DNM) transition between the Hi-bright and CO-bright media; (ii) constrain the CO-to-H2 conversion factor, XCO; and (iii) probe the dust properties per gas nucleon in each phase and map their spatial variations across the clouds. We have separated clouds at local, intermediate, and Galactic velocities in Hi and 12CO line emission to model in parallel the γ-ray intensity recorded between 0.4 and 100 GeV; the dust optical depth at 353 GHz, τ353; the thermal radiance of the large grains; and an estimate of the dust extinction, AVQ, empirically corrected for the starlight intensity. The dust and γ-ray models have been coupled to account for the DNM gas. The consistent γ-ray emissivity spectra recorded in the different phases confirm that the GeV-TeV cosmic rays probed by the LAT uniformly permeate all gas phases up to the 12CO cores. The dust and cosmic rays both reveal large amounts of DNM gas, with comparable spatial distributions and twice as much mass as in the CO-bright clouds. We give constraints on the Hi-DNM-CO transitions for five separate clouds. CO-dark H2 dominates the molecular columns up to AV ≃ 0.9 and its mass often exceeds the one-third of the molecular mass expected by theory. The corrected AVQ extinction largely provides the best fit to the total gas traced by the γ rays. Nevertheless, we find evidence for a marked rise in AVQ/NH with increasing NH and molecular fraction, and with decreasing dust temperature. The rise in τ353/NH is even steeper. We observe variations of lesser amplitude and orderliness for the specific power of the grains, except for a coherent decline by half in the CO cores. This combined information suggests grain evolution. We provide average values for the dust properties per gas nucleon in the different phases. The γ rays and dust radiance yield consistent XCO estimates near 0.7 × 1020cm-2K-1km-1s. The AVQ and τ353 tracers yield biased values because of the large rise in grain opacity in the CO clouds. These results clarify a recurrent disparity in the γ-ray versus dust calibration of XCO, but they confirm the factor of 2 difference found between the XCO estimates in nearby clouds and in the neighbouring spiral arms.

We present the results of approximately three years of observations of Planck Sunyaev-Zeldovich (SZ) sources with the Russian-Turkish 1.5 m telescope (RTT150), as a part of the optical follow-up programme undertaken by the Planck collaboration. During this time period approximately 20% of all dark and grey clear time available at the telescope was devoted to observations of Planck objects. Some observations of distant clusters were also done at the 6 m Bolshoi Telescope Alt-azimutalnyi (BTA) of the Special Astrophysical Observatory of the Russian Academy of Sciences. In total, deep, direct images of more than one hundred fields were obtained in multiple filters. We identified 47 previously unknown galaxy clusters, 41 of which are included in the Planck catalogue of SZ sources. The redshifts of 65 Planck clusters were measured spectroscopically and 14 more were measured photometrically. We discuss the details of cluster optical identifications and redshift measurements. We also present new spectroscopic redshifts for 39 Planck clusters that were not included in the Planck SZ source catalogue and are published here for the first time.

We have used the Planck all-sky submillimetre and millimetre maps to search for rare sources distinguished by extreme brightness, a few hundred millijanskies, and their potential for being situated at high redshift. These "cold" Planck sources, selected using the High Frequency Instrument (HFI) directly from the maps and from the Planck Catalogue of Compact Sources (PCCS), all satisfy the criterion of having their rest-frame far-infrared peak redshifted to the frequency range 353-857 GHz. This colour-selection favours galaxies in the redshift range z = 2-4, which we consider as cold peaks in the cosmic infrared background. With a 4.′5 beam at the four highest frequencies, our sample is expected to include overdensities of galaxies in groups or clusters, lensed galaxies, and chance line-of-sight projections. We perform a dedicated Herschel-SPIRE follow-up of 234 such Planck targets, finding a significant excess of red 350 and 500 μm sources, in comparison to reference SPIRE fields. About 94% of the SPIRE sources in the Planck fields are consistent with being overdensities of galaxies peaking at 350 μm, with 3% peaking at 500 μm, and none peaking at 250 μm. About 3% are candidate lensed systems, all 12 of which have secure spectroscopic confirmations, placing them at redshifts z > 2.2. Only four targets are Galactic cirrus, yielding a success rate in our search strategy for identifying extragalactic sources within the Planck beam of better than 98%. The galaxy overdensities are detected with high significance, half of the sample showing statistical significance above 10σ. The SPIRE photometric redshifts of galaxies in overdensities suggest a peak at z ≃ 2, assuming a single common dust temperature for the sources of Td = 35 K. Under this assumption, we derive an infrared (IR) luminosity for each SPIRE source of about 4 × 1012 L⊙, yielding star formation rates of typically 700 M⊙ yr-1. If the observed overdensities are actual gravitationally-bound structures, the total IR luminosity of all their SPIRE-detected sources peaks at 4 × 1013 L⊙, leading to total star formation rates of perhaps 7 × 103 M⊙yr-1 per overdensity. Taken together, these sources show the signatures of high-z (z > 2) protoclusters of intensively star-forming galaxies. All these observations confirm the uniqueness of our sample compared to reference samples and demonstrate the ability of the all-sky Planck-HFI cold sources to select populations of cosmological and astrophysical interest for structure formation studies.

The nearby Chamaeleon clouds have been observed in γ rays by the Fermi Large Area Telescope (LAT) and in thermal dust emission by Planck and IRAS. Cosmic rays and large dust grains, if smoothly mixed with gas, can jointly serve with the Hi and 12CO radio data to (i) map the hydrogen column densities, NH, in the different gas phases, in particular at the dark neutral medium (DNM) transition between the Hi-bright and CO-bright media; (ii) constrain the CO-to-H2 conversion factor, XCO; and (iii) probe the dust properties per gas nucleon in each phase and map their spatial variations across the clouds. We have separated clouds at local, intermediate, and Galactic velocities in Hi and 12CO line emission to model in parallel the γ-ray intensity recorded between 0.4 and 100 GeV; the dust optical depth at 353 GHz, τ353; the thermal radiance of the large grains; and an estimate of the dust extinction, AVQ, empirically corrected for the starlight intensity. The dust and γ-ray models have been coupled to account for the DNM gas. The consistent γ-ray emissivity spectra recorded in the different phases confirm that the GeV-TeV cosmic rays probed by the LAT uniformly permeate all gas phases up to the 12CO cores. The dust and cosmic rays both reveal large amounts of DNM gas, with comparable spatial distributions and twice as much mass as in the CO-bright clouds. We give constraints on the Hi-DNM-CO transitions for five separate clouds. CO-dark H2 dominates the molecular columns up to AV ≃ 0.9 and its mass often exceeds the one-third of the molecular mass expected by theory. The corrected AVQ extinction largely provides the best fit to the total gas traced by the γ rays. Nevertheless, we find evidence for a marked rise in AVQ/NH with increasing NH and molecular fraction, and with decreasing dust temperature. The rise in τ353/NH is even steeper. We observe variations of lesser amplitude and orderliness for the specific power of the grains, except for a coherent decline by half in the CO cores. This combined information suggests grain evolution. We provide average values for the dust properties per gas nucleon in the different phases. The γ rays and dust radiance yield consistent XCO estimates near 0.7 × 1020cm-2K-1km-1s. The AVQ and τ353 tracers yield biased values because of the large rise in grain opacity in the CO clouds. These results clarify a recurrent disparity in the γ-ray versus dust calibration of XCO, but they confirm the factor of 2 difference found between the XCO estimates in nearby clouds and in the neighbouring spiral arms.

The Andromeda galaxy (M 31) is one of a few galaxies that has sufficient angular size on the sky to be resolved by the Planck satellite. Planck has detected M 31 in all of its frequency bands, and has mapped out the dust emission with the High Frequency Instrument, clearly resolving multiple spiralarms and sub-features. We examine the morphology of this long-wavelength dust emission as seen by Planck, including a study of its outermost spiral arms, and investigate the dust heating mechanism across M 31. We find that dust dominating the longer wavelength emission (≥ 0.3 mm) is heated by the diffuse stellar population (as traced by 3.6 μm emission), with the dust dominating the shorter wavelength emission heated by a mix of the old stellar population and star-forming regions (as traced by 24 μm emission). We also fit spectral energy distributions for individual 5′ pixels and quantify the dust properties across the galaxy, taking into account these different heating mechanisms, finding that there is a linear decrease in temperature with galactocentric distance for dust heated by the old stellar population, as would be expected, with temperatures ranging from around 22 K in the nucleus to 14 K outside of the 10 kpc ring. Finally, we measure the integrated spectrum of the whole galaxy, which we find to be well-fitted with a global dust temperature of (18.2 ± 1.0) K with a spectral index of 1.62 ± 0.11 (assuming a single modified blackbody), and a significant amount of free-free emission at intermediate frequencies of 20-60 GHz, which corresponds to a star formation rate of around 0.12 M⊙ yr-1. We find a 2.3σ detection of the presence of spinning dust emission, with a 30 GHz amplitude of 0.7 ± 0.3 Jy, which is in line with expectations from our Galaxy.

We have used the Planck all-sky submillimetre and millimetre maps to search for rare sources distinguished by extreme brightness, a few hundred millijanskies, and their potential for being situated at high redshift. These "cold" Planck sources, selected using the High Frequency Instrument (HFI) directly from the maps and from the Planck Catalogue of Compact Sources (PCCS), all satisfy the criterion of having their rest-frame far-infrared peak redshifted to the frequency range 353-857 GHz. This colour-selection favours galaxies in the redshift range z = 2-4, which we consider as cold peaks in the cosmic infrared background. With a 4.′5 beam at the four highest frequencies, our sample is expected to include overdensities of galaxies in groups or clusters, lensed galaxies, and chance line-of-sight projections. We perform a dedicated Herschel-SPIRE follow-up of 234 such Planck targets, finding a significant excess of red 350 and 500 μm sources, in comparison to reference SPIRE fields. About 94% of the SPIRE sources in the Planck fields are consistent with being overdensities of galaxies peaking at 350 μm, with 3% peaking at 500 μm, and none peaking at 250 μm. About 3% are candidate lensed systems, all 12 of which have secure spectroscopic confirmations, placing them at redshifts z > 2.2. Only four targets are Galactic cirrus, yielding a success rate in our search strategy for identifying extragalactic sources within the Planck beam of better than 98%. The galaxy overdensities are detected with high significance, half of the sample showing statistical significance above 10σ. The SPIRE photometric redshifts of galaxies in overdensities suggest a peak at z ≃ 2, assuming a single common dust temperature for the sources of Td = 35 K. Under this assumption, we derive an infrared (IR) luminosity for each SPIRE source of about 4 × 1012 L⊙, yielding star formation rates of typically 700 M⊙ yr-1. If the observed overdensities are actual gravitationally-bound structures, the total IR luminosity of all their SPIRE-detected sources peaks at 4 × 1013 L⊙, leading to total star formation rates of perhaps 7 × 103 M⊙yr-1 per overdensity. Taken together, these sources show the signatures of high-z (z > 2) protoclusters of intensively star-forming galaxies. All these observations confirm the uniqueness of our sample compared to reference samples and demonstrate the ability of the all-sky Planck-HFI cold sources to select populations of cosmological and astrophysical interest for structure formation studies.

We update the all-sky Planck catalogue of 1227 clusters and cluster candidates (PSZ1) published in March 2013, derived from detections of the Sunyaev-Zeldovich (SZ) effect using the first 15.5 months of Planck satellite observations. As an addendum, we deliver an updated version of the PSZ1 catalogue, reporting the further confirmation of 86 Planck-discovered clusters. In total, the PSZ1 now contains 947 confirmed clusters, of which 214 were confirmed as newly discovered clusters through follow-up observations undertaken by the Planck Collaboration. The updated PSZ1 contains redshifts for 913 systems, of which 736 (~ 80.6%) are spectroscopic, and associated mass estimates derived from the Yz mass proxy. We also provide a new SZ quality flag for the remaining 280 candidates. This flag was derived from a novel artificial neural-network classification of the SZ signal. Based on this assessment, the purity of the updated PSZ1 catalogue is estimated to be 94%. In this release, we provide the full updated catalogue and an additional readme file with further information on the Planck SZ detections.

We update the all-sky Planck catalogue of 1227 clusters and cluster candidates (PSZ1) published in March 2013, derived from detections of the Sunyaev-Zeldovich (SZ) effect using the first 15.5 months of Planck satellite observations. As an addendum, we deliver an updated version of the PSZ1 catalogue, reporting the further confirmation of 86 Planck-discovered clusters. In total, the PSZ1 now contains 947 confirmed clusters, of which 214 were confirmed as newly discovered clusters through follow-up observations undertaken by the Planck Collaboration. The updated PSZ1 contains redshifts for 913 systems, of which 736 (~ 80.6%) are spectroscopic, and associated mass estimates derived from the Yz mass proxy. We also provide a new SZ quality flag for the remaining 280 candidates. This flag was derived from a novel artificial neural-network classification of the SZ signal. Based on this assessment, the purity of the updated PSZ1 catalogue is estimated to be 94%. In this release, we provide the full updated catalogue and an additional readme file with further information on the Planck SZ detections.

Planck data when combined with ancillary data provide a unique opportunity to separate the diffuse emission components of the inner Galaxy. The purpose of the paper is to elucidate the morphology of the various emission components in the strong star-formation region lying inside the solar radius and to clarify the relationship between the various components. The region of the Galactic plane covered is l = 300° → 0° → 60° wherestar-formation is highest and the emission is strong enough to make meaningful component separation. The latitude widths in this longitude range lie between 1 and 2, which correspond to FWHM z-widths of 100-200 pc at a typical distance of 6 kpc. The four emission components studied here are synchrotron, free-free, anomalous microwave emission (AME), and thermal (vibrational) dust emission. These components are identified by constructing spectral energy distributions (SEDs) at positions along the Galactic plane using the wide frequency coverage of Planck (28.4-857GHz) in combination with low-frequency radio data at 0.408-2.3 GHz plus WMAP data at 23-94 GHz, along with far-infrared (FIR) data from COBE-DIRBE and IRAS. The free-free component is determined from radio recombination line (RRL) data. AME is found to be comparable in brightness to the free-free emission on the Galactic plane in the frequency range 20-40 GHz with a width in latitude similar to that of the thermal dust; it comprises 45 ± 1% of the total 28.4 GHz emission in the longitude range l = 300° → 0° → 60°. The free-free component is the narrowest, reflecting the fact that it is produced by current star-formation as traced by the narrow distribution of OB stars. It is the dominant emission on the plane between 60 and 100 GHz. RRLs from this ionized gas are used to assess its distance, leading to a free-free z-width of FWHM ≈ 100 pc. The narrow synchrotron component has a low-frequency brightness spectral index βsynch ≈ -2.7 that is similar to the broad synchrotron component indicating that they are both populated by the cosmic ray electrons of the same spectral index. The width of this narrow synchrotron component is significantly larger than that of the other three components, suggesting that it is generated in an assembly of older supernova remnants that have expanded to sizes of order 150 pc in 3 × 105 yr; pulsars of a similar age have a similar spread in latitude. The thermal dust is identified in the SEDs with average parameters of Tdust = 20.4 ± 0.4 K, βFIR = 1.94 ± 0.03 (> 353 GHz), and βmm = 1.67 ± 0.02 (< 353 GHz). The latitude distributions of gamma-rays, CO, and the emission in high-frequency Planck bands have similar widths, showing that they are all indicators of the total gaseous matter on the plane in the inner Galaxy.

Any variation in the fundamental physical constants, more particularly in the fine structure constant, α, or in the mass of the electron, me, affects the recombination history of the Universe and cause an imprint on the cosmic microwave background angular power spectra. We show that the Planck data allow one to improve the constraint on the time variation of the fine structure constant at redshift z ∼ 103 by about a factor of 5 compared to WMAP data, as well as to break the degeneracy with the Hubble constant, H0. In addition to α, we can set a constraint on the variation in the mass of the electron, me, and in the simultaneous variation of the two constants. We examine in detail the degeneracies between fundamental constants and the cosmological parameters, in order to compare the limits obtained from Planck and WMAP and to determine the constraining power gained by including other cosmological probes. We conclude that independent time variations of the fine structure constant and of the mass of the electron are constrained by Planck to Δα/α = (3.6 ± 3.7) × 10-3 and Δme/me = (4 ± 11) × 10-3 at the 68% confidence level. We also investigate the possibility of a spatial variation of the fine structure constant. The relative amplitude of a dipolar spatial variation in α (corresponding to a gradient across our Hubble volume) is constrained to be δα/α = (-2.4 ± 3.7) × 10-2.

This paper presents an overview of the polarized sky as seen by Planck HFI at 353 GHz, which is the most sensitive Planck channel for dust polarization. We construct and analyse maps of dust polarization fraction and polarization angle at 1° resolution, taking into account noise bias and possible systematic effects. The sensitivity of the Planck HFI polarization measurements allows for the first time a mapping of Galactic dust polarized emission on large scales, including low column density regions. We find that the maximum observed dust polarization fraction is high (pmax = 19.8%), in particular in some regions of moderate hydrogen column density (NH < 2 × 1021 cm-2). The polarization fraction displays a large scatter at NH below a few 1021 cm-2. There is a general decrease in the dust polarization fraction with increasing column density above NH ∼ 1 × 1021 cm-2 and in particular a sharp drop above NH ∼ 1.5 × 1022 cm-2. We characterize the spatial structure of the polarization angle using the angle dispersion function. We find that the polarization angle is ordered over extended areas of several square degrees, separated by filamentary structures of high angle dispersion function. These appear as interfaces where the sky projection of the magnetic field changes abruptly without variations in the column density. The polarization fraction is found to be anti-correlated with the dispersion of polarization angles. These results suggest that, at the resolution of 1°, depolarization is due mainly to fluctuations in the magnetic field orientation along the line of sight, rather than to the loss of grain alignment in shielded regions. We also compare the polarization of thermal dust emission with that of synchrotron measured with Planck, low-frequency radio data, and Faraday rotation measurements toward extragalactic sources. These components bear resemblance along the Galactic plane and in some regions such as the Fan and North Polar Spur regions. The poor match observed in other regions shows, however, that dust, cosmic-ray electrons, and thermal electrons generally sample different parts of the line of sight.

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, pV, and the optical depth in the V band to the star, τV. Toward these stars we measure the submillimetre polarized intensity, PS, and total intensity, IS, 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 diffuse 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, RS/V = (PS/IS)/(pV/τV) and RP/p = PS/pV, the latter focusing directly on the polarization properties of the aligned grain population alone. We find RS/V = 4.2, with statistical and systematic uncertainties 0.2 and 0.3, respectively, and RP/p = 5.4 MJy sr-1, with uncertainties 0.2 and 0.3 MJy sr-1, respectively. Our estimate of RS/V is compatible with predictions based on a range of polarizing dust models that have been developed for the diffuse 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 RP/p is not compatible with predictions, which are too low by a factor of about 2.5. This more discriminating diagnostic, RP/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.

Polarized emission observed by Planck HFI at 353? GHz towards a sample of nearby fields is presented, focusing on the statistics of polarization fractions p and angles ψ. 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 pmax decreases with column density NH in the more opaque fields with NH> 1021 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 pmax 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.