No Arabic abstract
Using the Planck 2015 data release (PR2) temperature maps, we separate Galactic thermal dust emission from cosmic infrared background (CIB) anisotropies. For this purpose, we implement a specifically tailored component-separation method, the so-called generalized needlet internal linear combination (GNILC) method, which uses spatial information (the angular power spectra) to disentangle the Galactic dust emission and CIB anisotropies. We produce significantly improved all-sky maps of Planck thermal dust emission, with reduced CIB contamination, at 353, 545, and 857 GHz. By reducing the CIB contamination of the thermal dust maps, we provide more accurate estimates of the local dust temperature and dust spectral index over the sky with reduced dispersion, especially at high Galactic latitudes above $b = pm 20{deg}$. We find that the dust temperature is $T = (19.4 pm 1.3)$ K and the dust spectral index is $beta = 1.6 pm 0.1$ averaged over the whole sky, while $T = (19.4 pm 1.5)$ K and $beta = 1.6 pm 0.2$ on 21 % of the sky at high latitudes. Moreover, subtracting the new CIB-removed thermal dust maps from the CMB-removed Planck maps gives access to the CIB anisotropies over 60 % of the sky at Galactic latitudes $|b| > 20{deg}$. Because they are a significant improvement over previous Planck products, the GNILC maps are recommended for thermal dust science. The new CIB maps can be regarded as indirect tracers of the dark matter and they are recommended for exploring cross-correlations with lensing and large-scale structure optical surveys. The reconstructed GNILC thermal dust and CIB maps are delivered as Planck products.
Using Planck maps of six regions of low Galactic dust emission with a total area of about 140 square degrees, we determine the angular power spectra of cosmic infrared background (CIB) anisotropies from multipole l = 200 to l = 2000 at 217, 353, 545 and 857 GHz. We use 21-cm observations of HI as a tracer of thermal dust emission to reduce the already low level of Galactic dust emission and use the 143 GHz Planck maps in these fields to clean out cosmic microwave background anisotropies. Both of these cleaning processes are necessary to avoid significant contamination of the CIB signal. We measure correlated CIB structure across frequencies. As expected, the correlation decreases with increasing frequency separation, because the contribution of high-redshift galaxies to CIB anisotropies increases with wavelengths. We find no significant difference between the frequency spectrum of the CIB anisotropies and the CIB mean, with Delta I/I=15% from 217 to 857 GHz. In terms of clustering properties, the Planck data alone rule out the linear scale- and redshift-independent bias model. Non-linear corrections are significant. Consequently, we develop an alternative model that couples a dusty galaxy, parametric evolution model with a simple halo-model approach. It provides an excellent fit to the measured anisotropy angular power spectra and suggests that a different halo occupation distribution is required at each frequency, which is consistent with our expectation that each frequency is dominated by contributions from different redshifts. In our best-fit model, half of the anisotropy power at l=2000 comes from redshifts z<0.8 at 857 GHz and z<1.5 at 545 GHz, while about 90% come from redshifts z>2 at 353 and 217 GHz, respectively.
We use analytic computations to predict the power spectrum as well as the bispectrum of Cosmic Infrared Background (CIB) anisotropies. Our approach is based on the halo model and takes into account the mean luminosity-mass relation. The model is used to forecast the possibility to simultaneously constrain cosmological, CIB and halo occupation distribution (HOD) parameters in the presence of foregrounds. For the analysis we use wavelengths in eight frequency channels between 200 and 900$;mathrm{GHz}$ with survey specifications given by Planck and LiteBird. We explore the sensitivity to the model parameters up to multipoles of $ell =1000$ using auto- and cross-correlations between the different frequency bands. With this setting, cosmological, HOD and CIB parameters can be constrained to a few percent. Galactic dust is modeled by a power law and the shot noise contribution as a frequency dependent amplitude which are marginalized over. We find that dust residuals in the CIB maps only marginally influence constraints on standard cosmological parameters. Furthermore, the bispectrum yields tighter constraints (by a factor four in $1sigma$ errors) on almost all model parameters while the degeneracy directions are very similar to the ones of the power spectrum. The increase in sensitivity is most pronounced for the sum of the neutrino masses. Due to the similarity of degeneracies a combination of both analysis is not needed for most parameters. This, however, might be due to the simplified bias description generally adopted in such halo model approaches.
The polarized thermal emission from Galactic dust is the main foreground present in measurements of the polarization of the cosmic microwave background (CMB) at frequencies above 100GHz. We exploit the Planck HFI polarization data from 100 to 353GHz to measure the dust angular power spectra $C_ell^{EE,BB}$ over the range $40<ell<600$ well away from the Galactic plane. These will bring new insights into interstellar dust physics and a precise determination of the level of contamination for CMB polarization experiments. We show that statistical properties of the emission can be characterized over large fractions of the sky using $C_ell$. For the dust, they are well described by power laws in $ell$ with exponents $alpha^{EE,BB}=-2.42pm0.02$. The amplitudes of the polarization $C_ell$ vary with the average brightness in a way similar to the intensity ones. The dust polarization frequency dependence is consistent with modified blackbody emission with $beta_d=1.59$ and $T_d=19.6$K. We find a systematic ratio between the amplitudes of the Galactic $B$- and $E$-modes of 0.5. We show that even in the faintest dust-emitting regions there are no clean windows where primordial CMB $B$-mode polarization could be measured without subtraction of dust emission. Finally, we investigate the level of dust polarization in the BICEP2 experiment field. Extrapolation of the Planck 353GHz data to 150GHz gives a dust power $ell(ell+1)C_ell^{BB}/(2pi)$ of $1.32times10^{-2}mu$K$_{CMB}^2$ over the $40<ell<120$ range; the statistical uncertainty is $pm0.29$ and there is an additional uncertainty (+0.28,-0.24) from the extrapolation, both in the same units. This is the same magnitude as reported by BICEP2 over this $ell$ range, which highlights the need for assessment of the polarized dust signal even in the cleanest windows of the sky.
We use Planck HFI data combined with ancillary radio data to study the emissivity index of the interstellar dust emission in the frequency range 10 - 353 GHz, or 3 - 0.8 mm, in the Galactic plane. We analyse the region l=20 degr - 44 degr and |b| leq 4 degr where the free-free emission can be estimated from radio recombination line data. We fit the spectra at each sky pixel with a modified blackbody model and two spectral indices, beta_mm and beta_FIR, below and above 353 GHz respectively. We find that beta_mm is smaller than beta_FIR and we detect a correlation between this low frequency power-law index and the dust optical depth at 353 GHz, tau_353. The opacity spectral index beta_mm increases from about 1.54 in the more diffuse regions of the Galactic disk, |b| = 3 degr - 4 degr and tau_353 ~ 5 x 10^{-5}, to about 1.66 in the densest regions with an optical depth of more than one order of magnitude higher. We associate this correlation with an evolution of the dust emissivity related to the fraction of molecular gas along the line of sight. This translates into beta_mm ~ 1.54 for a medium that is mostly atomic and beta_mm ~ 1.66 when the medium is dominated by molecular gas. We find that both the Two-Level System model and the emission by ferromagnetic particles can explain the results. The results improve our understanding of the physics of interstellar dust and lead towards a complete model of the dust spectrum of the Milky Way from far-infrared to millimetre wavelengths.
We present a linear clustering model of cosmic infrared background (CIB) anisotropies at large scales that is used to measure the cosmic star formation rate density up to redshift 6, the effective bias of the CIB and the mass of dark-matter halos hosting dusty star-forming galaxies. This is achieved using the Planck CIB auto- and cross-power spectra (between different frequencies) and CIBxCMB lensing cross-spectra measurements, as well as external constraints (e.g. on the CIB mean brightness). We recovered an obscured star formation history which agrees well with the values derived from infrared deep surveys and we confirm that the obscured star formation dominates the unobscured one up to at least z=4. The obscured and unobscured star formation rate densities are compatible at $1sigma$ at z=5. We also determined the evolution of the effective bias of the galaxies emitting the CIB and found a rapid increase from $sim$0.8 at z$=$0 to $sim$8 at z$=$4. At 2$<$z$<$4, this effective bias is similar to that of galaxies at the knee of the mass functions and submillimeter galaxies. This effective bias is the weighted average of the true bias with the corresponding emissivity of the galaxies. The halo mass corresponding to this bias is thus not exactly the mass contributing the most to the star formation density. Correcting for this, we obtained a value of log(M$_h$/M$_{odot}$)=12.77$_{-0.125}^{+0.128}$ for the mass of the typical dark matter halo contributing to the CIB at z=2. Finally, we also computed using a Fisher matrix analysis how the uncertainties on the cosmological parameters affect the recovered CIB model parameters and find that the effect is negligible.