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.

We perform an analysis of the diffuse low-frequency Galactic components in the Southern part of the Gould Belt system (130^\circ\leq l\leq 230^\circ and -50^\circ\leq b\leq -10^\circ). Strong ultra-violet (UV) flux coming from the Gould Belt super-association is responsible for bright diffuse foregrounds that we observe from our position inside the system and that can help us improve our knowledge of the Galactic emission. Free-free emission and anomalous microwave emission (AME) are the dominant components at low frequencies ( u < 40 GHz), while synchrotron emission is very smooth and faint. We separate diffuse free-free emission and AME from synchrotron emission and thermal dust emission by using Planck data, complemented by ancillary data, using the "Correlated Component Analysis" (CCA) component separation method and we compare with the results of cross-correlation of foreground templates with the frequency maps. We estimate the electron temperature T_e from H$\alpha$ and free-free emission using two methods (temperature-temperature plot and cross-correlation) and we obtain T_e ranging from 3100 to 5200 K, for an effective fraction of absorbing dust along the line of sight of 30% (f_d=0.3). We estimate the frequency spectrum of the diffuse AME and we recover a peak frequency (in flux density units) of 25.5 \pm 1.5 GHz. We verify the reliability of this result with realistic simulations that include the presence of biases in the spectral model for the AME and in the free-free template. By combining physical models for vibrational and rotational dust emission and adding the constraints from the thermal dust spectrum from Planck and IRAS we are able to get a good description of the frequency spectrum of the AME for plausible values of the local density and radiation field.

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 ∼ 10³ 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.

Taking advantage of the all-sky coverage and broadfrequency range of the Planck satellite, we study the Sunyaev-Zeldovich (SZ) and pressure profiles of 62 nearby massive clusters detected at high significance in the 14-month nominal survey. Careful reconstruction of the SZ signal indicates that most clusters are individually detected at least out to R500. By stacking the radial profiles, we have statistically detected the radial SZ signal out to 3 × R500, i.e., at a density contrast of about 50-100, though the dispersion about the mean profile dominates the statistical errors across the whole radial range. Our measurement is fully consistent with previous Planck results on integrated SZ fluxes, further strengthening the agreement between SZ and X-ray measurements inside R500. Correcting for the effects of the Planck beam, we have calculated the corresponding pressure profiles. This new constraint from SZ measurements is consistent with the X-ray constraints from XMM-Newton in the region in which the profiles overlap (i.e., [0.1-1] R 500), and is in fairly good agreement with theoretical predictions within the expected dispersion. At larger radii the average pressure profile is slightly flatter than most predictions from numerical simulations. Combining the SZ and X-ray observed profiles into a joint fit to a generalised pressure profile gives best-fit parameters [P0,c500,γ, α,β] = [6.41,1.81,0.31,1.33,4.13]. Using a reasonable hypothesis for the gas temperature in the cluster outskirts we reconstruct from our stacked pressure profile the gas mass fraction profile out to 3 R500. Within the temperature driven uncertainties, our Planck constraints are compatible with the cosmic baryon fraction and expected gas fraction in halos. © 2013 ESO.

The dust-Hi correlation is used to characterize the emission properties of dust in the diffuse interstellar medium (ISM) from far infrared wavelengths to microwave frequencies. The field of this investigation encompasses the part of the southern sky best suited to study the cosmic infrared and microwave backgrounds. We cross-correlate sky maps from Planck, the Wilkinson Microwave Anisotropy Probe (WMAP), and the diffuse infrared background experiment (DIRBE), at 17 frequencies from 23 to 3000 GHz, with the Parkes survey of the 21 cm line emission of neutral atomic hydrogen, over a contiguous area of 7500 deg2 centred on the southern Galactic pole. We present a general methodology to study the dust-H i correlation over the sky, including simulations to quantify uncertainties. Our analysis yields four specific results. (1) We map the temperature, submillimetre emissivity, and opacity of the dust per H-atom. The dust temperature is observed to be anti-correlated with the dust emissivity and opacity. We interpret this result as evidence of dust evolution within the diffuse ISM. The mean dust opacity is measured to be (7.1 ± 0.6) × 10-27 cm2 H-1 × (v/353 GHz) 1.53 ± 0.03for 100 ≤ v ≤ 353 GHz. This is a reference value to estimate hydrogen column densities from dust emission at submillimetre and millimetre wavelengths. (2) We map the spectral index βmm of dust emission at millimetre wavelengths (defined here as v ≤ 353GHz), and find it to be remarkably constant at βmm = 1.51 ± 0.13. We compare it with the far infrared spectral index βFIR derived from greybody fits at higher frequencies, and find a systematic difference, βmm -βFIR = -0.15, which suggests that the dust spectral energy distribution (SED) flattens at v ≤ 353 GHz. (3) We present spectral fits of the microwave emission correlated with Hi from 23 to 353 GHz, which separate dust and anomalous microwave emission (AME). We show that the flattening of the dust SED can be accounted for with an additional component with a blackbody spectrum. This additional component, which accounts for (26 ± 6)% of the dust emission at 100GHz, could represent magnetic dipole emission. Alternatively, it could account for an increasing contribution of carbon dust, or a flattening of the emissivity of amorphous silicates, at millimetre wavelengths. These interpretations make different predictions for the dust polarization SED. (4) We analyse the residuals of the dust-Hi correlation. We identify a Galactic contribution to these residuals, which we model with variations of the dust emissivity on angular scales smaller than that of our correlation analysis. This model of the residuals is used to quantify uncertainties of the CIB power spectrum in a companion Planck paper.© ESO 2014.

We describe the detection, interpretation, and removal of the signal resulting from interactions of high energy particles with the Planck High Frequency Instrument (HFI). There are two types of interactions: heating of the 0.1 K bolometer plate; and glitches in each detector time stream. The transientresponses to detector glitch shapes are not simple single-pole exponential decays and fall into three families. The glitch shape for each family has been characterized empirically in flight data and these shapes have been used to remove glitches from the detector time streams. The spectrum of the count rate per unit energy is computed for each family and a correspondence is made to the location on the detector of the particle hit. Most of the detected glitches are from Galactic protons incident on the die frame supporting the micro-machined bolometric detectors. In the Planck orbit at L2, the particle flux is around 5 cm-2 s-1 and is dominated by protons incident on the spacecraft with energy >39 MeV, at a rate of typically one event per second per detector. Different categories of glitches have different signatures in the time stream. Two of the glitch types have a low amplitude component that decays over nearly 1 s. This component produces excess noise if not properly removed from the time-ordered data. We have used a glitch detection and subtraction method based on the joint fit of population templates. The application of this novel glitch subtraction method removes excess noise from the time streams. Using realistic simulations, we find that this method does not introduce signal bias into the Planck data.

Using precise full-sky observations from Planck, and applying several methods of component separation, we identify and characterise the emission from the Galactic "haze" at microwave wavelengths. The haze is a distinct component of diffuse Galactic emission, roughly centered on the Galactic centre, and extends to | b | ∼ 35-50 in Galactic latitude and | l | ∼ 15-20 in longitude. By combining the Planck data with observations from the Wilkinson Microwave Anisotropy Probe, we were able to determine the spectrum of this emission to high accuracy, unhindered by the strong systematic biases present in previous analyses. The derived spectrum is consistent with power-law emission with a spectral index of-2.56 ± 0.05, thus excluding free-free emission as the source and instead favouring hard-spectrum synchrotron radiation from an electron population with a spectrum (number density per energy) dN/dE ∞ E-2.1. At Galactic latitudes | b | < 30, the microwave haze morphology is consistent with that of the Fermi gamma-ray "haze" or "bubbles", while at b ∼-50 we have identified an edge in the microwave haze that is spatially coincident with the edge in the gamma-ray bubbles. Taken together, this indicates that we have a multi-wavelength view of a distinct component of our Galaxy. Given both the very hard spectrum and the extended nature of the emission, it is highly unlikely that the haze electrons result from supernova shocks in the Galactic disk. Instead, a new astrophysical mechanism for cosmic-ray acceleration in the inner Galaxy is implied. © 2013 ESO.

We explore the 2013 Planck likelihood function with a high-precision multi-dimensional minimizer (Minuit). This allows a refinement of the CDM best-fit solution with respect to previously-released results, and the construction of frequentist confidence intervals using profile likelihoods. The agreement with the cosmological results from the Bayesian framework is excellent, demonstrating the robustness of the Planck results to the statistical methodology. We investigate the inclusion of neutrino masses, where more significant differences may appear due to the non-Gaussian nature of the posterior mass distribution. By applying the Feldman-Cousins prescription, we again obtain results very similar to those of the Bayesian methodology. However, the profile-likelihood analysis of the cosmic microwave background (CMB) combination (Planck+WP+highL) reveals a minimum well within the unphysical negative-mass region. We show that inclusion of the Planck CMB-lensing information regularizes this issue, and provide a robust frequentist upper limit σmv ≤0:26 eV (95% confidence) from the CMB+lensing+BAO data combination. © ESO 2014.

Maps of cosmic microwave background (CMB) temperature and polarization from the 2015 release of Planck data provide the highestquality full-sky view of the surface of last scattering available to date. This enables us to detect possible departures from a globally isotropic cosmology. We present the first searches using CMB polarization for correlations induced by a possible non-trivial topology with a fundamental domain that intersects, or nearly intersects, the last-scattering surface (at comoving distance χrec), both via a direct scan for matched circular patterns at the intersections and by an optimal likelihood calculation for specific topologies. We specialize to flat spaces with cubic toroidal (T3) and slab (T1) topologies, finding that explicit searches for the latter are sensitive to other topologies with antipodal symmetry. These searches yield no detection of a compact topology with a scale below the diameter of the last-scattering surface. The limits on the radius i of the largest sphere inscribed in the fundamental domain (at log-likelihood ratio Δ >-5 relative to a simply-connected flat Planck best-fit model) are: i > 0.97 χrec for the T3 cubic torus; and i > 0.56 χrec for the T1 slab. The limit for the T3 cubic torus from the matched-circles search is numerically equivalent, i > 0.97 χrec at 99% confidence level from polarization data alone. We also perform a Bayesian search for an anisotropic global Bianchi VIIh geometry. In the non-physical setting, where the Bianchi cosmology is decoupled from the standard cosmology, Planck temperature data favour the inclusion of a Bianchi component with a Bayes factor of at least 2.3 units of log-evidence. However, the cosmological parameters that generate this pattern are in strong disagreement with those found from CMB anisotropy data alone. Fitting the induced polarization pattern for this model to the Planck data requires an amplitude of-0.10 ± 0.04 compared to the value of + 1 if the model were to be correct. In the physically motivated setting, where the Bianchi parameters are coupled and fitted simultaneously with the standard cosmological parameters, we find no evidence for a Bianchi VIIh cosmology and constrain the vorticity of such models to (ω/H)0 < 7.6 × 10-10 (95% CL).