Source author record

L. Bonavera

L. Bonavera appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

90works
5topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

90 published item(s)

preprint2022arXiv

Selecting a complete sample of blazars in sub-millimetre catalogues

The \textit{Herschel} Astrophysical Terahertz Large Area Survey (H-ATLAS), that has covered about 642 sq. deg. in 5 bands from 100 to 500 $μ\rm m$, allows a blind flux-limited selection of blazars at sub-mm wavelengths. However, blazars constitute a tiny fraction of H-ATLAS sources and therefore identifying them is not a trivial task. Using the data on known blazars detected by the H-ATLAS we have defined a locus for 500\,$μ$m selected blazars and exploited it to select blazar candidates in the H-ATLAS fields. Candidates and known blazars in the H-ATLAS equatorial and South Galactic Pole fields were followed up with the Australia Telescope Compact Array (ATCA) or with the Karl G. Jansky Very Large Array (VLA), and matched with existing radio- and mm-catalogues to reconstruct the spectral behaviour over at least 6 orders of magnitude in frequency. We identified a selection approach that, combining the information in the sub-mm and radio domains, efficiently singles out genuine blazars. In this way, we identified a sample of 39 blazars brighter than $S_{500μ\rm m} = 35\,$mJy in the H-ATLAS fields. Tests made cross-matching the H-ATLAS catalogues with large catalogues of blazar candidates indicate that the sample is complete. The derived counts are compared with model predictions finding good consistency with the C2Ex model and with estimates based on ALMA data.

preprint2022arXiv

Tomography-based observational measurements of the halo mass function via the submillimeter magnification bias

Aims. The main goal of this paper is to derive observational constraints on the halo mass fuction (HMF) by performing a tomographic analysis of the magnification bias signal on a sample of background submillimeter galaxies. The results can then be compared with those from a non-tomographic study. Methods. We measure the cross-correlation function between a sample of foreground GAMA galaxies with spectroscopic redshifts in the range $0.1 < z < 0.8$ (and divided up into four bins) and a sample of background submillimeter galaxies from H-ATLAS with photometric redshifts in the range $1.2 < z < 4.0$. We model the weak lensing signal within the halo model formalism and carry out a Markov chain Monte Carlo algorithm to obtain the posterior distribution of all HMF parameters, which we assume to follow the Sheth and Tormen (ST) three-parameter and two-parameter fits. Results. While the observational constraints on the HMF from the non-tomographic analysis are not stringent, there is a remarkable improvement in terms of uncertainty reduction when tomography is adopted. Moreover, with respect to the traditional ST triple of values from numerical simulations, the results from the three-parameter fit predict a higher number density of halos at masses below $10^{12}M_{\odot}/h$ at 95% credibility. The two-parameter fit yields even more restricting results, with a larger number density of halos below $10^{13}M_{\odot}/h$ and a lower one above $10^{14}M_{\odot}/h$, this time at more than 3$σ$ credibility. Our results are therefore in disagreement with the standard N-body values for the ST fit at 2$σ$ and 3$σ$, respectively.

preprint2021arXiv

A direct and robust method to observationally constrain the halo mass function via the submillimeter magnification bias: Proof of concept

Aims. The main purpose of this work is to provide a method to derive tabulated observational constraints on the halo mass function (HMF) by studying the magnification bias effect on high-redshift submillimeter galaxies. Under the assumption of universality, we parametrize the HMF according to two traditional models, namely the Sheth and Tormen (ST) and Tinker fits and assess their performance in explaining the measured data within the Λ cold dark matter (ΛCDM) model. We also study the potential influence of the halo occupation distribution (HOD) parameters in this analysis and discuss two important aspects regarding the HMF parametrization. Methods. We measure the cross-correlation function between a foreground sample of GAMA galaxies with redshifts in the range $0.2<z<0.8$ and a background sample of H-ATLAS galaxies with redshifts in the range $1.2<z<4.0$ and carry out an MCMC algorithm to check this observable against its mathematical prediction within the halo model formalism. Results. If all HMF parameters are assumed to be positive, the ST fit only seems to fully explain the measurements by forcing the mean number of satellite galaxies in a halo to increase substantially from its prior mean value. The Tinker fit, on the other hand, provides a robust description of the data without relevant changes in the HOD parameters, but with some dependence on the prior range of two of its parameters. When the normalization condition for the HMF is dropped and we allow negative values of the $p_1$ parameter in the ST fit, all the involved parameters are better determined, unlike the previous models, thus deriving the most general HMF constraints. While all cases are in agreement with the traditional fits within the uncertainties, the last one hints at a slightly higher number of halos at intermediate and high masses, raising the important point of the allowed parameter range.

preprint2021arXiv

Galaxy cluster mass density profile derived using the submillimetre galaxies magnification bias

In this work we want to study the average mass density profile of tens to hundreds of clusters of galaxies acting as lenses that produce a magnification bias on the SMGs, and to estimate their associated masses and concentrations for different richness ranges. The background sample is composed of SMGs observed by Herschel with 1.2<z<4.0 (mean redshift at ~2.3) while the foreground sample is made up of galaxy clusters extracted from the SDSS III with photometric redshifts of 0.05< z< 0.8 (mean redshift at ~0.38). Measurements are obtained by stacking the SMG--cluster pairs to estimate the cross-correlation function using the Davis-Peebles estimator. This methodology allows us to derive the mass density profile for a wide range of angular scales, ~2-250 arcsec or ~10-1300 kpc for z=0.38, with a high radial resolution. We find that It is impossible to fit the data with a single mass density profile at all scales. As for the outer part, the estimated average masses increase from $M_{200c}=5.8$ to $51.5\times 10^{13} M_\odot$ and the concentration parameter from C=0.74 to 1.74. In the small-scale regions, the obtained average masses fluctuate around $M_{200c}=3-4 \times 10^{13}M_\odot$ with average C~4. The total average masses are in perfect agreement with the M-R relationship estimated from the cluster catalogue. While the estimated average C values of the central galactic halos are in agreement with traditional M-C relationships, we find low concentrations for the outer part. Moreover, C decrease for lower R values, probably indicating that the group of galaxies cannot be considered to be relaxed systems. Finally, we notice a systematic lack of signal at the transition between the dominance of the cluster halo and the central galactic halo (~100 kpc). This feature is also present in previous studies using different catalogues and/or methodologies, but is never discussed.

preprint2021arXiv

Multi-frequency point source detection with fully convolutional networks: Performance in realistic microwave sky simulations

Point Source (PS) detection is an important issue for future Cosmic Microwave Background (CMB) experiments since they are one of the main contaminants to the recovery of CMB signal at small scales. Improving its multifrequency detection would allow to take into account valuable information otherwise neglected when extracting PS using a channel-by-channel approach. We develop a method based on Neural Networks (NNs) to detect PS in multifrequency realistic simulations and compare its performance against one of the most popular methods, the matrix filters. The frequencies used are 143, 217 and 353 GHz and we impose a Galactic cut of 30 degrees. We produce simulations by adding contaminating signals to the PS maps as the CMB, the Cosmic Infrared Background, the Galactic thermal emission, the thermal Sunyaev-Zel'dovich effect and the instrumental noise. These simulations are used to train two NNs called Flat and Spectral MultiPoSeIDoN. The first one considers PS with a flat spectrum and the second one is more realistic because it takes into account the spectral behavior of the PS. Using a detection limit of 60 mJy, Flat MultiPoSeIDoN reachs the 90% of completeness level at 58 mJy and at 79, 71 and 60 for the spectral case at 143, 217 and 353 GHz respectively, while the matrix filters reach it at 84, 79 and 123 mJy. Using safer 4σ detection limit does not help to improve these results. In all cases, MultiPoSeIDoN obtain a much lower number of spurious sources than the filter. The NNs recover the flux density of the detections with a relative error of 10% above 100 mJy, while the filter above 150 mJy. Based on the results, NNs are the perfect candidates to substitute filters to detect multifrequency PS in future CMB experiments. Moreover, we have shown that a multifrequency approach can detect sources with higher accuracy than single-frequency approaches also based on NNs.

preprint2020arXiv

Cosmological constraints with the sub-millimetre galaxies Magnification Bias after large scale bias corrections

The study of the magnification bias produced on high-redshift sub-millimetre galaxies by foreground galaxies through the analysis of the cross-correlation function was recently demonstrated as an interesting independent alternative to the weak-lensing shear as a cosmological probe. In the case of the proposed observable, most of the cosmological constraints mainly depend on the largest angular separation measurements. Therefore, we aim to study and correct the main large-scale biases that affect foreground and background galaxy samples to produce a robust estimation of the cross-correlation function. Then we analyse the corrected signal to derive updated cosmological constraintsWe measured the large-scale, bias-corrected cross-correlation functions using a background sample of H-ATLAS galaxies with photometric redshifts > 1.2 and two different foreground samples (GAMA galaxies with spectroscopic redshifts or SDSS galaxies with photometric ones, both in the range 0.2 < z < 0.8). These measurements are modelled using the traditional halo model description that depends on both halo occupation distribution and cosmological parameters. We then estimated these parameters by performing a Markov chain Monte Carlo under multiple scenarios to study the performance of this observable and how to improve its results. After the large-scale bias corrections, we obtain only minor improvements with respect to the previous magnification bias results, mainly confirming their conclusions: a lower bound on $Ω_m > 0.22$ at $95\%$ C.L. and an upper bound $σ_8 < 0.97$ at $95\%$ C.L. (results from the $z_{spec}$ sample). However, by combining both foreground samples into a simplified tomographic analysis, we were able to obtain interesting constraints on the $Ω_m$-$σ_8$ plane as follows: $Ω_m= 0.50_{- 0.20}^{+ 0.14}$ and $σ_8= 0.75_{- 0.10}^{+ 0.07}$ at 68\% CL.

preprint2020arXiv

Cosmology with the submillimetre galaxies magnification bias: Proof of concept

Context. As recently demonstrated, high-z submillimetre galaxies (SMGs) are the perfect background sample for tracing the mass density profiles of galaxies and clusters (baryonic and dark matter) and their time-evolution through gravitational lensing. Their magnification bias, a weak gravitational lensing effect, is a powerful tool for constraining the free parameters of a halo occupation distribution (HOD) model and potentially also some of the main cosmological parameters. Aims. The aim of this work is to test the capability of the magnification bias produced on high-z SMGs as a cosmological probe. We exploit cross-correlation data to constrain not only astrophysical parameters ($M_{min}$, $M_1$, and $α$), but also some of the cosmological ones ($Ω_m$, $σ_8$, and $H_0$) for this proof of concept. Methods. The measured cross-correlation function between a foreground sample of GAMA galaxies with spectroscopic redshifts in the range 0.2 < z < 0.8 and a background sample of H-ATLAS galaxies with photometric redshifts >1.2 is modelled using the traditional halo model description that depends on HOD and cosmological parameters. These parameters are then estimated by performing a Markov chain Monte Carlo analysis using different sets of priors to test the robustness of the results and to study the performance of this novel observable with the current set of data Results. With our current results, $Ω_m$ and $H_0$ cannot be well constrained. However, we can set a lower limit of >0.24 at 95\% confidence level (CL) on $Ω_m$ and we see a slight trend towards $H_0>70$ values. For our constraints on $σ_8$ we obtain only a tentative peak around 0.75, but an interesting upper limit of $σ_8\lesssim 1$ at 95\% CL. We also study the possibility to derive better constraints by imposing more restrictive priors on the astrophysical parameters.

preprint2020arXiv

The Hyper-Kamiokande Experiment -- Snowmass LOI

Hyper-Kamiokande is the next generation underground water Cherenkov detector that builds on the highly successful Super-Kamiokande experiment. The detector which has an 8.4~times larger effective volume than its predecessor will be located along the T2K neutrino beamline and utilize an upgraded J-PARC beam with 2.6~times beam power. Hyper-K's low energy threshold combined with the very large fiducial volume make the detector unique, that is expected to acquire an unprecedented exposure of 3.8~Mton$\cdot$year over a period of 20~years of operation. Hyper-Kamiokande combines an extremely diverse science program including nucleon decays, long-baseline neutrino oscillations, atmospheric neutrinos, and neutrinos from astrophysical origins. The scientific scope of this program is highly complementary to liquid-argon detectors for example in sensitivity to nucleon decay channels or supernova detection modes. Hyper-Kamiokande construction has started in early 2020 and the experiment is expected to start operations in 2027. The Hyper-Kamiokande collaboration is presently being formed amongst groups from 19 countries including the United States, whose community has a long history of making significant contributions to the neutrino physics program in Japan. US physicists have played leading roles in the Kamiokande, Super-Kamiokande, EGADS, K2K, and T2K programs.

preprint2019arXiv

A transiting super-Earth close to the inner edge of the habitable zone of an M0 dwarf star

We present a super-Earth orbiting close to the inner edge of the habitable zone of the cool dwarf star K2-286 (EPIC 249889081), detected with data from the K2 mission in its $15^{th}$ campaign. The planet has radius of $2.1\pm0.2$ R$_{\oplus}$, near the 1.5 - 2.0 R$_{\oplus}$ gap in the radii distribution. The equilibrium temperature is $347^{+21}_{-11}$ K, cooler than most of the small planets with well measured masses, and the orbital period is $27.359\pm0.005$ days. K2-286, located at a distance of $76.3\pm0.3$ pc, is an M0V star with estimated effective temperature of $3926\pm100$ K, less active than other M dwarf stars hosting exoplanets. The expected radial velocity semi-amplitude induced by the planet on the star is $1.9^{+1.3}_{-1.2}$ m$\cdot$s$^{-1}$, and the amplitude of signals in transit transmission spectroscopy is estimated at $5.0\pm3.0$ ppm. Follow-up observations for mass measurements and transit spectroscopy should be desirable for this relatively bright target ($m_V=12.76, m_{Ks}=9.32$) hosting a transiting super-Earth within the inner edge of the habitable zone.

preprint2018arXiv

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.

preprint2016arXiv

On the recovery of ISW fluctuations using large-scale structure tracers and CMB temperature and polarization anisotropies

In this work we present a method to extract the signal induced by the integrated Sachs-Wolfe (ISW) effect in the cosmic microwave background (CMB). It makes use of the Linear Covariance-Based filter introduced by Barreiro et al., and combines CMB data with any number of large-scale structure (LSS) surveys and lensing information. It also exploits CMB polarization to reduce cosmic variance. The performance of the method has been thoroughly tested with simulations taking into account the impact of non-ideal conditions such as incomplete sky coverage or the presence of noise. In particular, three galaxy surveys are simulated, whose redshift distributions peak at low ($z \simeq 0.3$), intermediate ($z \simeq 0.6$) and high redshift ($z \simeq 0.9$). The contribution of each of the considered data sets as well as the effect of a mask and noise in the reconstructed ISW map is studied in detail. When combining all the considered data sets (CMB temperature and polarization, the three galaxy surveys and the lensing map), the proposed filter successfully reconstructs a map of the weak ISW signal, finding a perfect correlation with the input signal for the ideal case and around 80 per cent, on average, in the presence of noise and incomplete sky coverage. We find that including CMB polarization improves the correlation between input and reconstruction although only at a small level. Nonetheless, given the weakness of the ISW signal, even modest improvements can be of importance. In particular, in realistic situations, in which less information is available from the LSS tracers, the effect of including polarisation is larger. For instance, for the case in which the ISW signal is recovered from CMB plus only one survey, and taking into account the presence of noise and incomplete sky coverage, the improvement in the correlation coefficient can be as large as 10 per cent.

preprint2016arXiv

Planck 2015 results. II. Low Frequency Instrument data processing

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.

preprint2016arXiv

Planck 2015 results. III. LFI systematic uncertainties

We present the current accounting of systematic effect uncertainties for the Low Frequency Instrument (LFI) that are relevant to the 2015 release of the Planck cosmological results, showing the robustness and consistency of our data set, especially for polarization analysis. We use two complementary approaches: (i) simulations based on measured data and physical models of the known systematic effects; and (ii) analysis of difference maps containing the same sky signal ("null-maps"). The LFI temperature data are limited by instrumental noise. At large angular scales the systematic effects are below the cosmic microwave background (CMB) temperature power spectrum by several orders of magnitude. In polarization the systematic uncertainties are dominated by calibration uncertainties and compete with the CMB $E$-modes in the multipole range 10--20. Based on our model of all known systematic effects, we show that these effects introduce a slight bias of around $0.2\,σ$ on the reionization optical depth derived from the 70\,GHz $EE$ spectrum using the 30 and 353\,GHz channels as foreground templates. At 30\,GHz the systematic effects are smaller than the Galactic foreground at all scales in temperature and polarization, which allows us to consider this channel as a reliable template of synchrotron emission. We assess the residual uncertainties due to LFI effects on CMB maps and power spectra after component separation and show that these effects are smaller than the CMB amplitude at all scales. We also assess the impact on non-Gaussianity studies and find it to be negligible. Some residuals still appear in null maps from particular sky survey pairs, particularly at 30 GHz, suggesting possible straylight contamination due to an imperfect knowledge of the beam far sidelobes.

preprint2016arXiv

Planck 2015 results. XI. CMB power spectra, likelihoods, and robustness of parameters

This paper presents the Planck 2015 likelihoods, statistical descriptions of the 2-point correlations of CMB data, using the hybrid approach employed previously: pixel-based at $\ell<30$ and a Gaussian approximation to the distribution of spectra at higher $\ell$. The main improvements are the use of more and better processed data and of Planck polarization data, and more detailed foreground and instrumental models, allowing further checks and enhanced immunity to systematics. Progress in foreground modelling enables a larger sky fraction. Improvements in processing and instrumental models further reduce uncertainties. For temperature, we perform an analysis of end-to-end instrumental simulations fed into the data processing pipeline; this does not reveal biases from residual instrumental systematics. The $Λ$CDM cosmological model continues to offer a very good fit to Planck data. The slope of primordial scalar fluctuations, $n_s$, is confirmed smaller than unity at more than 5σ from Planck alone. We further validate robustness against specific extensions to the baseline cosmology. E.g., the effective number of neutrino species remains compatible with the canonical value of 3.046. This first detailed analysis of Planck polarization concentrates on E modes. At low $\ell$ we use temperature at all frequencies and a subset of polarization. The frequency range improves CMB-foreground separation. Within the baseline model this requires a reionization optical depth $τ=0.078\pm0.019$, significantly lower than without high-frequency data for explicit dust monitoring. At high $\ell$ we detect residual errors in E, typically O($μ$K$^2$); we recommend temperature alone as the high-$\ell$ baseline. Nevertheless, Planck high-$\ell$ polarization allows a separate determination of $Λ$CDM parameters consistent with those from temperature alone.

preprint2016arXiv

Planck 2015 results. XIII. Cosmological parameters

We present results based on full-mission Planck observations of temperature and polarization anisotropies of the CMB. These data are consistent with the six-parameter inflationary LCDM cosmology. From the Planck temperature and lensing data, for this cosmology we find a Hubble constant, H0= (67.8 +/- 0.9) km/s/Mpc, a matter density parameter Omega_m = 0.308 +/- 0.012 and a scalar spectral index with n_s = 0.968 +/- 0.006. (We quote 68% errors on measured parameters and 95% limits on other parameters.) Combined with Planck temperature and lensing data, Planck LFI polarization measurements lead to a reionization optical depth of tau = 0.066 +/- 0.016. Combining Planck with other astrophysical data we find N_ eff = 3.15 +/- 0.23 for the effective number of relativistic degrees of freedom and the sum of neutrino masses is constrained to < 0.23 eV. Spatial curvature is found to be |Omega_K| < 0.005. For LCDM we find a limit on the tensor-to-scalar ratio of r <0.11 consistent with the B-mode constraints from an analysis of BICEP2, Keck Array, and Planck (BKP) data. Adding the BKP data leads to a tighter constraint of r < 0.09. We find no evidence for isocurvature perturbations or cosmic defects. The equation of state of dark energy is constrained to w = -1.006 +/- 0.045. Standard big bang nucleosynthesis predictions for the Planck LCDM cosmology are in excellent agreement with observations. We investigate annihilating dark matter and deviations from standard recombination, finding no evidence for new physics. The Planck results for base LCDM are in agreement with BAO data and with the JLA SNe sample. However the amplitude of the fluctuations is found to be higher than inferred from rich cluster counts and weak gravitational lensing. Apart from these tensions, the base LCDM cosmology provides an excellent description of the Planck CMB observations and many other astrophysical data sets.

preprint2016arXiv

Planck 2015 results. XIV. Dark energy and modified gravity

We study the implications of Planck data for models of dark energy (DE) and modified gravity (MG), beyond the cosmological constant scenario. We start with cases where the DE only directly affects the background evolution, considering Taylor expansions of the equation of state, principal component analysis and parameterizations related to the potential of a minimally coupled DE scalar field. When estimating the density of DE at early times, we significantly improve present constraints. We then move to general parameterizations of the DE or MG perturbations that encompass both effective field theories and the phenomenology of gravitational potentials in MG models. Lastly, we test a range of specific models, such as k-essence, f(R) theories and coupled DE. In addition to the latest Planck data, for our main analyses we use baryonic acoustic oscillations, type-Ia supernovae and local measurements of the Hubble constant. We further show the impact of measurements of the cosmological perturbations, such as redshift-space distortions and weak gravitational lensing. These additional probes are important tools for testing MG models and for breaking degeneracies that are still present in the combination of Planck and background data sets. All results that include only background parameterizations are in agreement with LCDM. When testing models that also change perturbations (even when the background is fixed to LCDM), some tensions appear in a few scenarios: the maximum one found is \sim 2 sigma for Planck TT+lowP when parameterizing observables related to the gravitational potentials with a chosen time dependence; the tension increases to at most 3 sigma when external data sets are included. It however disappears when including CMB lensing.

preprint2016arXiv

Planck 2015 results. XIX. Constraints on primordial magnetic fields

We compute and investigate four types of imprint of a stochastic background of primordial magnetic fields (PMFs) on the cosmic microwave background (CMB) anisotropies: the impact of PMFs on the CMB spectra; the effect on CMB polarization induced by Faraday rotation; the impact of PMFs on the ionization history; magnetically-induced non-Gaussianities; and the magnetically-induced breaking of statistical isotropy. Overall, Planck data constrain the amplitude of PMFs to less than a few nanogauss. In particular, individual limits coming from the analysis of the CMB angular power spectra, using the Planck likelihood, are $B_{1\,\mathrm{Mpc}}< 4.4$ nG (where $B_{1\,\mathrm{Mpc}}$ is the comoving field amplitude at a scale of 1 Mpc) at 95% confidence level, assuming zero helicity, and $B_{1\,\mathrm{Mpc}}< 5.6$ nG for a maximally helical field.For nearly scale-invariant PMFs we obtain $B_{1\,\mathrm{Mpc}}<2.0$ nG and $B_{1\,\mathrm{Mpc}}<0.9$ nG if the impact of PMFs on the ionization history of the Universe is included. From the analysis of magnetically-induced non-Gaussianity we obtain three different values, corresponding to three applied methods, all below 5 nG. The constraint from the magnetically-induced passive-tensor bispectrum is $B_{1\,\mathrm{Mpc}}< 2.8$ nG. A search for preferred directions in the magnetically-induced passive bispectrum yields $B_{1\,\mathrm{Mpc}}< 4.5$ nG, whereas the the compensated-scalar bispectrum gives $B_{1\,\mathrm{Mpc}}< 3$ nG. The analysis of the Faraday rotation of CMB polarization by PMFs uses the Planck power spectra in $EE$ and $BB$ at 70 GHz and gives $B_{1\,\mathrm{Mpc}}< 1380$ nG. In our final analysis, we consider the harmonic-space correlations produced by Alfvén waves, finding no significant evidence for the presence of these waves. Together, these results comprise a comprehensive set of constraints on possible PMFs with Planck data.

preprint2016arXiv

Planck 2015 results. XV. Gravitational lensing

We present the most significant measurement of the cosmic microwave background (CMB) lensing potential to date (at a level of 40 sigma), using temperature and polarization data from the Planck 2015 full-mission release. Using a polarization-only estimator we detect lensing at a significance of 5 sigma. We cross-check the accuracy of our measurement using the wide frequency coverage and complementarity of the temperature and polarization measurements. Public products based on this measurement include an estimate of the lensing potential over approximately 70% of the sky, an estimate of the lensing potential power spectrum in bandpowers for the multipole range 40<L<400 and an associated likelihood for cosmological parameter constraints. We find good agreement between our measurement of the lensing potential power spectrum and that found in the best-fitting LCDM model based on the Planck temperature and polarization power spectra. Using the lensing likelihood alone we obtain a percent-level measurement of the parameter combination $σ_8 Ω_m^{0.25} = 0.591\pm 0.021$. We combine our determination of the lensing potential with the E-mode polarization also measured by Planck to generate an estimate of the lensing B-mode. We show that this lensing B-mode estimate is correlated with the B-modes observed directly by Planck at the expected level and with a statistical significance of 10 sigma, confirming Planck's sensitivity to this known sky signal. We also correlate our lensing potential estimate with the large-scale temperature anisotropies, detecting a cross-correlation at the 3 sigma level, as expected due to dark energy in the concordance LCDM model.

preprint2016arXiv

Planck 2015 results. XVI. Isotropy and statistics of the CMB

We test the statistical isotropy and Gaussianity of the cosmic microwave background (CMB) anisotropies using observations made by the Planck satellite. Our results are based mainly on the full Planck mission for temperature, but also include some polarization measurements. In particular, we consider the CMB anisotropy maps derived from the multi-frequency Planck data by several component-separation methods. For the temperature anisotropies, we find excellent agreement between results based on these sky maps over both a very large fraction of the sky and a broad range of angular scales, establishing that potential foreground residuals do not affect our studies. Tests of skewness, kurtosis, multi-normality, N-point functions, and Minkowski functionals indicate consistency with Gaussianity, while a power deficit at large angular scales is manifested in several ways, for example low map variance. The results of a peak statistics analysis are consistent with the expectations of a Gaussian random field. The "Cold Spot" is detected with several methods, including map kurtosis, peak statistics, and mean temperature profile. We thoroughly probe the large-scale dipolar power asymmetry, detecting it with several independent tests, and address the subject of a posteriori correction. Tests of directionality suggest the presence of angular clustering from large to small scales, but at a significance that is dependent on the details of the approach. We perform the first examination of polarization data, finding the morphology of stacked peaks to be consistent with the expectations of statistically isotropic simulations. Where they overlap, these results are consistent with the Planck 2013 analysis based on the nominal mission data and provide our most thorough view of the statistics of the CMB fluctuations to date.

preprint2016arXiv

Planck 2015 results. XVII. Constraints on primordial non-Gaussianity

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 local, equilateral, and orthogonal bispectrum amplitudes, quoting as our final result from temperature alone fNL^local=2.5+\-5.7, fNL^equil=-16+\-70 and fNL^ortho=-34+\-33(68%CL). Combining temperature and polarization data we obtain fNL^local=0.8+\-5.0, fNL^equil=-4+\-43 and fNL^ortho=-26+\-21 (68%CL). The results are based on cross-validation of these estimators on simulations, are stable across component separation techniques, pass an extensive suite of tests, and are consistent with Minkowski functionals based measurements. 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 fNL estimates, we present model-independent reconstructions of the CMB bispectrum and derive constraints on early universe scenarios that generate NG, including general single-field and axion inflation, initial state modifications, parity-violating tensor bispectra, and directionally dependent vector models. We also present a wide survey of scale-dependent oscillatory bispectra, and we look for isocurvature NG. Our constraint on the local primordial trispectrum amplitude is gNL^local=(-9.0+\-7.7)x10^4 (68%CL), and we perform an analysis of additional trispectrum shapes. The global picture is one of consistency with the premises of the LambdaCDM cosmology, namely that the structure we observe today was sourced by adiabatic, passive, Gaussian, and primordial seed perturbations.[abridged]

preprint2016arXiv

Planck 2015 results. XVIII. Background geometry & topology

Full-sky CMB maps from the 2015 Planck release allow us to detect departures from global isotropy on the largest scales. We present the first searches using CMB polarization for correlations induced by a non-trivial topology with a fundamental domain intersecting, or nearly intersecting, the last scattering surface (at comoving distance $χ_{rec}$). We specialize to flat spaces with toroidal and slab 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 at a scale below the diameter of the last scattering surface. The limits on the radius $R_i$ of the largest sphere inscribed in the topological domain (at log-likelihood-ratio $Δ\ln{L}>-5$ relative to a simply-connected flat Planck best-fit model) are $R_i>0.97χ_{rec}$ for the cubic torus and $R_i>0.56χ_{rec}$ for the slab. The limit for the cubic torus from the matched-circles search is numerically equivalent, $R_i>0.97χ_{rec}$ (99% CL) from polarisation data alone. We also perform a Bayesian search for a Bianchi VII$_h$ 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. However, the cosmological parameters generating this pattern are in strong disagreement with those found from CMB anisotropy data alone. Fitting the induced polarization pattern for this model to Planck data requires an amplitude of $-0.1\pm0.04$ compared to +1 if the model were to be correct. In the physical setting where the Bianchi parameters are fit simultaneously with the standard cosmological parameters, we find no evidence for a Bianchi VII$_h$ cosmology and constrain the vorticity of such models to $(ω/H)_0<7.6\times10^{-10}$ (95% CL). [Abridged]

preprint2016arXiv

Planck 2015 results. XXI. The integrated Sachs-Wolfe effect

This paper presents a study of the ISW effect from the Planck 2015 temperature and polarization data release. The CMB is cross-correlated with different LSS tracers: the NVSS, SDSS and WISE catalogues, and the Planck 2015 lensing map. This cross-correlation yields a detection at $4\,σ$, where most of the signal-to-noise is due to the Planck lensing and NVSS. In fact, the ISW effect is detected only from the Planck data (through the ISW-lensing bispectrum) at $\approx 3\,σ$, which is similar to the detection level achieved by combining the cross-correlation signal coming from all the catalogues. The ISW signal allow us to detect $Ω_Λ$ 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, polarization data is used to study the anomalously large ISW signal previously reported through the stacking of CMB features at the locations of known superstructures. 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 LSS traced by the 2MPZ survey by directly inverting the density field into the gravitational potential field.

preprint2016arXiv

Planck 2015 results. XXV. Diffuse low-frequency Galactic foregrounds

(abridged) We discuss the Galactic foreground emission between 20 and 100GHz based on observations by Planck/WMAP. The Commander component-separation tool has been used to separate the various astrophysical processes in total intensity. Comparison with RRL templates verifies the recovery of the free-free emission along the Galactic plane. Comparison of the high-latitude Halpha emission with our free-free map shows residuals that correlate with dust optical depth, consistent with a fraction (~30%) of Halpha having been scattered by high-latitude dust. We highlight a number of diffuse spinning dust morphological features at high latitude. There is substantial spatial variation in the spinning dust spectrum, with the emission peak ranging from below 20GHz to more than 50GHz. There is a strong tendency for the spinning dust component near many prominent HII regions to have a higher peak frequency, suggesting that this increase in peak frequency is associated with dust in the photodissociation regions around the nebulae. The emissivity of spinning dust in these diffuse regions is of the same order as previous detections in the literature. Over the entire sky, the commander solution finds more anomalous microwave emission than the WMAP component maps, at the expense of synchrotron and free-free emission. This can be explained by the difficulty in separating multiple broadband components with a limited number of frequency maps. Future surveys (5-20GHz), will greatly improve the separation by constraining the synchrotron spectrum. We combine Planck/WMAP data to make the highest S/N ratio maps yet of the intensity of the all-sky polarized synchrotron emission at frequencies above a few GHz. Most of the high-latitude polarized emission is associated with distinct large-scale loops and spurs, and we re-discuss their structure...

preprint2016arXiv

Planck 2015 results. XXVI. The Second Planck Catalogue of Compact Sources

The Second Planck Catalogue of Compact Sources is a catalogue of sources detected in single-frequency maps from the full duration of the Planck mission and supersedes previous versions of the Planck compact source catalogues. It consists of compact sources, both Galactic and extragalactic, detected over the entire sky. Compact sources detected in the lower frequency channels are assigned to the PCCS2, while at higher frequencies they are assigned to one of two sub-catalogues, the PCCS2 or PCCS2E, depending on their location on the sky. The first of these catalogues covers most of the sky and allows the user to produce subsamples at higher reliabilities than the target 80% integral reliability of the catalogue. The PCCS2E contains sources detected in sky regions where the diffuse emission makes it difficult to quantify the reliability of the detections. Both the PCCS2 and PCCS2E include polarization measurements, in the form of polarized flux densities, or upper limits, and orientation angles for all seven polarization-sensitive Planck channels. The improved data-processing of the full-mission maps and their reduced noise levels allow us to increase the number of objects in the catalogue, improving its completeness for the target 80 % reliability as compared with the previous versions, the PCCS and ERCSC catalogues.

preprint2016arXiv

Planck intermediate results. XLI. A map of lensing-induced B-modes

The secondary cosmic microwave background (CMB) $B$-modes stem from the post-decoupling distortion of the polarization $E$-modes due to the gravitational lensing effect of large-scale structures. These lensing-induced $B$-modes constitute both a valuable probe of the dark matter distribution and an important contaminant for the extraction of the primary CMB $B$-modes from inflation. Planck provides accurate nearly all-sky measurements of both the polarization $E$-modes and the integrated mass distribution via the reconstruction of the CMB lensing potential. By combining these two data products, we have produced an all-sky template map of the lensing-induced $B$-modes using a real-space algorithm that minimizes the impact of sky masks. The cross-correlation of this template with an observed (primordial and secondary) $B$-mode map can be used to measure the lensing $B$-mode power spectrum at multipoles up to $2000$. In particular, when cross-correlating with the $B$-mode contribution directly derived from the Planck polarization maps, we obtain lensing-induced $B$-mode power spectrum measurement at a significance level of $12\,σ$, which agrees with the theoretical expectation derived from the Planck best-fit $Λ$CDM model. This unique nearly all-sky secondary $B$-mode template, which includes the lensing-induced information from intermediate to small ($10\lesssim \ell\lesssim 1000$) angular scales, is delivered as part of the Planck 2015 public data release. It will be particularly useful for experiments searching for primordial $B$-modes, such as BICEP2/Keck Array or LiteBIRD, since it will enable an estimate to be made of the lensing-induced contribution to the measured total CMB $B$-modes.

preprint2016arXiv

Planck intermediate results. XLII. Large-scale Galactic magnetic fields

Recent models for the large-scale Galactic magnetic fields in the literature have been largely constrained by synchrotron emission and Faraday rotation measures. We use three different but representative models to compare their predicted polarized synchrotron and dust emission with that measured by the Planck satellite. We first update these models to match the Planck synchrotron products using a common model for the cosmic-ray leptons. We discuss the impact on this analysis of the ongoing problems of component separation in the Planck microwave bands and of the uncertain cosmic-ray spectrum. In particular, the inferred degree of ordering in the magnetic fields is sensitive to these systematic uncertainties, and we further show the importance of considering the expected variations in the observables in addition to their mean morphology. We then compare the resulting simulated emission to the observed dust polarization and find that the dust predictions do not match the morphology in the Planck data but underpredict the dust polarization away from the plane. We modify one of the models to roughly match both observables at high latitudes by increasing the field ordering in the thin disc near the observer. Though this specific analysis is dependent on the component separation issues, we present the improved model as a proof of concept for how these studies can be advanced in future using complementary information from ongoing and planned observational projects.

preprint2016arXiv

Planck intermediate results. XLIV. The structure of the Galactic magnetic field from dust polarization maps of the southern Galactic cap

We study the statistical properties of interstellar dust polarization at high Galactic latitude, using the Stokes parameter Planck maps at 353 GHz. Our aim is to advance the understanding of the magnetized interstellar medium (ISM), and to provide a model of the polarized dust foreground for cosmic microwave background component-separation procedures. Focusing on the southern Galactic cap, we examine the statistical distributions of the polarization fraction ($p$) and angle ($ψ$) to characterize the ordered and turbulent components of the Galactic magnetic field (GMF) in the solar neighbourhood. We relate patterns at large angular scales in polarization to the orientation of the mean (ordered) GMF towards Galactic coordinates $(l_0,b_0)=(70^\circ \pm 5^\circ,24^\circ \pm 5^\circ)$. The histogram of $p$ shows a wide dispersion up to 25 %. The histogram of $ψ$ has a standard deviation of $12^\circ$ about the regular pattern expected from the ordered GMF. We use these histograms to build a phenomenological model of the turbulent component of the GMF, assuming a uniform effective polarization fraction ($p_0$) of dust emission. To model the Stokes parameters, we approximate the integration along the line of sight (LOS) as a sum over a set of $N$ independent polarization layers, in each of which the turbulent component of the GMF is obtained from Gaussian realizations of a power-law power spectrum. We are able to reproduce the observed $p$ and $ψ$ distributions using: a $p_0$ value of (26 $\pm$ 3)%; a ratio of 0.9 $\pm$ 0.1 between the strengths of the turbulent and mean components of the GMF; and a small value of $N$. We relate the polarization layers to the density structure and to the correlation length of the GMF along the LOS.

preprint2016arXiv

Planck intermediate results. XLIX. Parity-violation constraints from polarization data

Parity violating extensions of the standard electromagnetic theory cause in vacuo rotation of the plane of polarization of propagating photons. This effect, also known as cosmic birefringence, impacts the cosmic microwave background (CMB) anisotropy angular power spectra, producing non-vanishing $T$--$B$ and $E$--$B$ correlations that are otherwise null when parity is a symmetry. Here we present new constraints on an isotropic rotation, parametrized by the angle $α$, derived from Planck 2015 CMB polarization data. To increase the robustness of our analyses, we employ two complementary approaches, in harmonic space and in map space, the latter based on a peak stacking technique. The two approaches provide estimates for $α$ that are in agreement within statistical uncertainties and very stable against several consistency tests. Considering the $T$--$B$ and $E$--$B$ information jointly, we find $α= 0.31^{\circ} \pm 0.05^{\circ} \, ({\rm stat.})\, \pm 0.28^{\circ} \, ({\rm syst.})$ from the harmonic analysis and $α= 0.35^{\circ} \pm 0.05^{\circ} \, ({\rm stat.})\, \pm 0.28^{\circ} \, ({\rm syst.})$ from the stacking approach. These constraints are compatible with no parity violation and are dominated by the systematic uncertainty in the orientation of Planck's polarization-sensitive bolometers.

preprint2016arXiv

Planck intermediate results. XLV. Radio spectra of northern extragalactic radio sources

Continuum spectra covering centimetre to submillimetre wavelengths are presented for a northern sample of 104 extragalactic radio sources, mainly active galactic nuclei, based on four-epoch Planck data. The nine Planck frequencies, from 30 to 857 GHz, are complemented by a set of simultaneous ground-based radio observations between 1.1 and 37 GHz. The single-survey Planck data confirm that the flattest high-frequency radio spectral indices are close to zero, indicating that the original accelerated electron energy spectrum is much harder than commonly thought, with power-law index around 1.5 instead of the canonical 2.5. The radio spectra peak at high frequencies and exhibit a variety of shapes. For a small set of low-z sources, we find a spectral upturn at high frequencies, indicating the presence of intrinsic cold dust. Variability can generally be approximated by achromatic variations, while sources with clear signatures of evolving shocks appear to be limited to the strongest outbursts.

preprint2016arXiv

Planck intermediate results. XLVI. Reduction of large-scale systematic effects in HFI polarization maps and estimation of the reionization optical depth

This paper describes the identification, modelling, and removal of previously unexplained systematic effects in the polarization data of the Planck High Frequency Instrument (HFI) on large angular scales, including new mapmaking and calibration procedures, new and more complete end-to-end simulations, and a set of robust internal consistency checks on the resulting maps. These maps, at 100, 143, 217, and 353 GHz, are early versions of those that will be released in final form later in 2016. The improvements allow us to determine the cosmic reionization optical depth $τ$ using, for the first time, the low-multipole $EE$ data from HFI, reducing significantly the central value and uncertainty, and hence the upper limit. Two different likelihood procedures are used to constrain $τ$ from two estimators of the CMB $E$- and $B$-mode angular power spectra at 100 and 143 GHz, after debiasing the spectra from a small remaining systematic contamination. These all give fully consistent results. A further consistency test is performed using cross-correlations derived from the Low Frequency Instrument maps of the Planck 2015 data release and the new HFI data. For this purpose, end-to-end analyses of systematic effects from the two instruments are used to demonstrate the near independence of their dominant systematic error residuals. The tightest result comes from the HFI-based $τ$ posterior distribution using the maximum likelihood power spectrum estimator from $EE$ data only, giving a value $0.055\pm 0.009$. In a companion paper these results are discussed in the context of the best-fit Planck $Λ$CDM cosmological model and recent models of reionization.

preprint2016arXiv

Planck intermediate results. XLVII. Planck constraints on reionization history

We investigate constraints on cosmic reionization extracted from the Planck cosmic microwave background (CMB) data. We combine the Planck CMB anisotropy data in temperature with the low-multipole polarization data to fit LCDM models with various parameterizations of the reionization history. We obtain a Thomson optical depth tau=0.058 +/- 0.012 for the commonly adopted instantaneous reionization model. This confirms, with only data from CMB anisotropies, the low value suggested by combining Planck 2015 results with other data sets and also reduces the uncertainties. We reconstruct the history of the ionization fraction using either a symmetric or an asymmetric model for the transition between the neutral and ionized phases. To determine better constraints on the duration of the reionization process, we also make use of measurements of the amplitude of the kinetic Sunyaev-Zeldovich (kSZ) effect using additional information from the high resolution Atacama Cosmology Telescope and South Pole Telescope experiments. The average redshift at which reionization occurs is found to lie between z=7.8 and 8.8, depending on the model of reionization adopted. Using kSZ constraints and a redshift-symmetric reionization model, we find an upper limit to the width of the reionization period of Dz < 2.8. In all cases, we find that the Universe is ionized at less than the 10% level at redshifts above z~10. This suggests that an early onset of reionization is strongly disfavoured by the Planck data. We show that this result also reduces the tension between CMB-based analyses and constraints from other astrophysical sources.

preprint2016arXiv

Planck intermediate results. XLVIII. Disentangling Galactic dust emission and cosmic infrared background anisotropies

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°$. We find that the dust temperature is $T = (19.4 \pm 1.3)$ K and the dust spectral index is $β= 1.6 \pm 0.1$ averaged over the whole sky, while $T = (19.4 \pm 1.5)$ K and $β= 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°$. 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.

preprint2016arXiv

Planck intermediate results. XXIX. All-sky dust modelling with Planck, IRAS, and WISE observations

We present all-sky modelling of the high resolution Planck, IRAS, and WISE infrared (IR) observations using the physical dust model presented by Draine and Li in 2007 (DL). 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, the optical extinction Av, and the starlight intensity 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. We compare the DL optical extinction Av for the diffuse interstellar medium with optical estimates for 2 10^5 quasi-stellar objects (QSOs) observed in the Sloan digital sky survey. 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. 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 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. The DL model and the QSOs data are 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.

preprint2016arXiv

Planck Intermediate Results. XXXVI. Optical identification and redshifts of Planck SZ sources with telescopes in the Canary Islands Observatories

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).

preprint2015arXiv

A Joint Analysis of BICEP2/Keck Array and Planck Data

We report the results of a joint analysis of data from BICEP2/Keck Array and Planck. BICEP2 and Keck Array have observed the same approximately 400 deg$^2$ patch of sky centered on RA 0h, Dec. $-57.5°$. The combined maps reach a depth of 57 nK deg in Stokes $Q$ and $U$ in a band centered at 150 GHz. Planck has observed the full sky in polarization at seven frequencies from 30 to 353 GHz, but much less deeply in any given region (1.2 $μ$K deg in $Q$ and $U$ at 143 GHz). We detect 150$\times$353 cross-correlation in $B$-modes at high significance. We fit the single- and cross-frequency power spectra at frequencies $\geq 150$ GHz to a lensed-$Λ$CDM model that includes dust and a possible contribution from inflationary gravitational waves (as parameterized by the tensor-to-scalar ratio $r$), using a prior on the frequency spectral behavior of polarized dust emission from previous \planck\ analysis of other regions of the sky. We find strong evidence for dust and no statistically significant evidence for tensor modes. We probe various model variations and extensions, including adding a synchrotron component in combination with lower frequency data, and find that these make little difference to the $r$ constraint. Finally we present an alternative analysis which is similar to a map-based cleaning of the dust contribution, and show that this gives similar constraints. The final result is expressed as a likelihood curve for $r$, and yields an upper limit $r_{0.05}<0.12$ at 95% confidence. Marginalizing over dust and $r$, lensing $B$-modes are detected at $7.0\,σ$ significance.

preprint2015arXiv

Extragalactic sources in Cosmic Microwave Background maps

We discuss the potential of a next generation space-borne CMB experiment for studies of extragalactic sources with reference to COrE+, a project submitted to ESA in response to the M4 call. We consider three possible options for the telescope size: 1m, 1.5m and 2m (although the last option is probably impractical, given the M4 boundary conditions). The proposed instrument will be far more sensitive than Planck and will have a diffraction-limited angular resolution. These properties imply that even the 1m telescope option will perform substantially better than Planck for studies of extragalactic sources. The source detection limits as a function of frequency have been estimated by means of realistic simulations. The most significant improvements over Planck results are presented for each option. COrE+ will provide much larger samples of truly local star-forming galaxies, making possible analyses of the properties of galaxies (luminosity functions, dust mass functions, star formation rate functions, dust temperature distributions, etc.) across the Hubble sequence. Even more interestingly, COrE+ will detect, at |b|> 30 deg, thousands of strongly gravitationally lensed galaxies. Such large samples are of extraordinary astrophysical and cosmological value in many fields. Moreover, COrE+ high frequency maps will be optimally suited to pick up proto-clusters of dusty galaxies, i.e. to investigate the evolution of large scale structure at larger redshifts than can be reached by other means. Thanks to its high sensitivity COrE+ will also yield a spectacular advance in the blind detection of extragalactic sources in polarization. This will open a new window for studies of radio source polarization and of the global properties of magnetic fields in star forming galaxies and of their relationships with SFRs.

preprint2015arXiv

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

The European Space Agency's Planck satellite, dedicated to studying the early Universe and its subsequent evolution, was launched 14~May 2009 and 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 on data 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 science products include maps of the cosmic microwave background (CMB), the thermal Sunyaev-Zeldovich effect, and 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 the performance of the analysis methods and assessment of uncertainties. The likelihood code used to assess cosmological models against the Planck data are described, as well as a CMB lensing likelihood. Scientific results include cosmological parameters deriving 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.

preprint2015arXiv

Planck 2015 results. IV. Low Frequency Instrument beams and window functions

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 $\ell \approx 600$); and 0.5% at 70 GHz (at $\ell \approx 1000$).

preprint2015arXiv

Planck 2015 results. IX. Diffuse component separation: CMB maps

We present foreground-reduced 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 $\ell\gtrsim40$. On the very largest scales, instrumental systematic residuals are still non-negligible compared to the expected cosmological signal, and modes with $\ell < 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$\,μ\textrm{K}$ averaged over 55 arcmin pixels, and between 4.5 and 6.1$\,μ\textrm{K}$ averaged over 3.4 arcmin 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.

preprint2015arXiv

Planck 2015 results. V. LFI calibration

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 global 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.

preprint2015arXiv

Planck 2015 results. VI. LFI mapmaking

This paper describes the mapmaking procedure applied to Planck LFI (Low Frequency Instrument) 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, 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 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.

preprint2015arXiv

Planck 2015 results. VII. HFI TOI and beam processing

The Planck High Frequency Instrument (HFI) has observed the full sky at six frequencies (100, 143, 217, 353, 545, and 857 GHz) in intensity and at four frequencies in linear polarization (100, 143, 217, and 353 GHz). In order to obtain sky maps, the time-ordered information (TOI) containing the detector and pointing samples must be processed and the angular response must be assessed. The full mission TOI is included in the Planck 2015 release. This paper describes the HFI TOI and beam processing for the 2015 release. HFI calibration and map-making are described in a companion paper. The main pipeline has been modified since the last release (2013 nominal mission in intensity only), by including a correction for the non-linearity of the warm readout and by improving the model of the bolometer time response. The beam processing is an essential tool that derives the angular response used in all the Planck science papers and we report an improvement in the effective beam window function uncertainty of more than a factor 10 relative to the 2013 release. Noise correlations introduced by pipeline filtering function are assessed using dedicated simulations. Angular cross-power spectra using datasets that are decorrelated in time are immune to the main systematic effects.

preprint2015arXiv

Planck 2015 results. VIII. High Frequency Instrument data processing: Calibration and maps

This paper describes the processing applied to the Planck High Frequency Instrument (HFI) cleaned, time-ordered information to produce photometrically calibrated maps in temperature and (for the first time) in polarization. The data from the entire 2.5 year HFI mission include almost five independent full-sky surveys. HFI observes the sky over a broad range of frequencies, from 100 to 857 GHz. To obtain the best accuracy on the calibration over such a large range, two different photometric calibration schemes have been used. The 545 and 857 GHz data are calibrated using models of planetary atmospheric emission. The lower frequencies (from 100 to 353 GHz) are calibrated using the time-variable cosmological microwave background dipole, which we call the "orbital dipole". This source of calibration only depends on the satellite velocity with respect to the solar system. Using a CMB temperature of 2.7255 +/- 0.0006 K, it permits an independent measurement of the amplitude of the CMB solar dipole (3364.3 +/- 1.5 μK) which is approximatively 1σ higher than the WMAP measurement with a direction that is consistent between both experiments. We describe the pipeline used to produce the maps of intensity and linear polarization from the HFI timelines, and the scheme used to set the zero level of the maps a posteriori. We also summarize the noise characteristics of the HFI maps in the 2015 Planck data release and present some null tests to assess their quality. Finally, we discuss the major systematic effects and in particular the leakage induced by flux mismatch between the detectors that leads to spurious polarization signal.

preprint2015arXiv

Planck 2015 results. X. Diffuse component separation: Foreground maps

Planck has mapped the microwave sky in nine frequency bands between 30 and 857 GHz in temperature and seven bands between 30 and 353 GHz in polarization. In this paper we consider the problem of diffuse astrophysical component separation, and process these maps within a Bayesian framework to derive a consistent set of full-sky astrophysical component maps. For the temperature analysis, we combine the Planck observations with the 9-year WMAP sky maps and the Haslam et al. 408 MHz map to derive a joint model of CMB, synchrotron, free-free, spinning dust, CO, line emission in the 94 and 100 GHz channels, and thermal dust emission. Full-sky maps are provided with angular resolutions varying between 7.5 arcmin and 1 deg. Global parameters (monopoles, dipoles, relative calibration, and bandpass errors) are fitted jointly with the sky model, and best-fit values are tabulated. For polarization, the model includes CMB, synchrotron, and thermal dust emission. These models provide excellent fits to the observed data, with rms temperature residuals smaller than 4 uK over 93% of the sky for all Planck frequencies up to 353 GHz, and fractional errors smaller than 1% in the remaining 7% of the sky. The main limitations of the temperature model at the lower frequencies are degeneracies among the spinning dust, free-free, and synchrotron components; additional observations from external low-frequency experiments will be essential to break these. The main limitations of the temperature model at the higher frequencies are uncertainties in the 545 and 857 GHz calibration and zero-points. For polarization, the main outstanding issues are instrumental systematics in the 100-353 GHz bands on large angular scales in the form of temperature-to-polarization leakage, uncertainties in the analog-to-digital conversion, and very long time constant corrections, all of which are expected to improve in the near future.

preprint2015arXiv

Planck 2015 results. XII. Full Focal Plane simulations

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 18144 maps, together with the most massive suite of Monte Carlo realizations of instrument noise and CMB ever generated, comprising $10^4$ mission realizations reduced to about $10^6$ maps. The resulting maps incorporate the dominant instrumental, scanning, and data analysis effects; 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 for the validation and verification of analysis algorithms, as well as their implementations, and for removing biases from and quantifying uncertainties in the results of analyses of the real data.

preprint2015arXiv

Planck 2015 results. XXII. A map of the thermal Sunyaev-Zeldovich effect

We have constructed all-sky y-maps of the thermal Sunyaev-Zeldovich (tSZ) effect by applying specifically tailored component separation algorithms to the 30 to 857 GHz frequency channel maps from the Planck satellite survey. These reconstructed y-maps are delivered as part of the Planck 2015 release. The y-maps are characterised in terms of noise properties and residual foreground contamination, mainly thermal dust emission at large angular scales and CIB and extragalactic point sources at small angular scales. Specific masks are defined to minimize foreground residuals and systematics. Using these masks we compute the y-map angular power spectrum and higher order statistics. From these we conclude that the y-map is dominated by tSZ signal in the multipole range, 20-600. We compare the measured tSZ power spectrum and higher order statistics to various physically motivated models and discuss the implications of our results in terms of cluster physics and cosmology.

preprint2015arXiv

Planck 2015 results. XXIII. The thermal Sunyaev-Zeldovich effect--cosmic infrared background correlation

We use Planck data to detect the cross-correlation between the thermal Sunyaev-Zeldovich (tSZ) effect and the infrared emission from the galaxies that make up the the cosmic infrared background (CIB). We first perform a stacking analysis towards Planck-confirmed galaxy clusters. We detect infrared emission produced by dusty galaxies inside these clusters and demonstrate that the infrared emission is about 50% more extended than the tSZ effect. Modelling the emission with a Navarro--Frenk--White profile, we find that the radial profile concentration parameter is $c_{500} = 1.00^{+0.18}_{-0.15}$. This indicates that infrared galaxies in the outskirts of clusters have higher infrared flux than cluster-core galaxies. We also study the cross-correlation between tSZ and CIB anisotropies, following three alternative approaches based on power spectrum analyses: (i) using a catalogue of confirmed clusters detected in Planck data; (ii) using an all-sky tSZ map built from Planck frequency maps; and (iii) using cross-spectra between Planck frequency maps. With the three different methods, we detect the tSZ-CIB cross-power spectrum at significance levels of (i) 6 $σ$, (ii) 3 $σ$, and (iii) 4 $σ$. We model the tSZ-CIB cross-correlation signature and compare predictions with the measurements. The amplitude of the cross-correlation relative to the fiducial model is $A_{\rm tSZ-CIB}= 1.2\pm0.3$. This result is consistent with predictions for the tSZ-CIB cross-correlation assuming the best-fit cosmological model from Planck 2015 results along with the tSZ and CIB scaling relations.

preprint2015arXiv

Planck 2015 results. XXVII. The Second Planck Catalogue of Sunyaev-Zeldovich Sources

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 all-sky catalogue of 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 > $10^3$ 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 Y5R500 estimates are 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.

preprint2015arXiv

Planck 2015 results. XXVIII. The Planck Catalogue of Galactic Cold Clumps

We present the Planck Catalogue of Galactic Cold Clumps (PGCC), an all-sky catalogue of Galactic cold clump candidates detected by Planck. This catalogue is the full version of the Early Cold Core (ECC) catalogue, which was made available in 2011 with the Early Release Compact Source Catalogue (ERCSC) and contained 915 high S/N sources. It is based on the Planck 48 months mission data that are currently being released to the astronomical community. The PGCC catalogue is an observational catalogue consisting exclusively of Galactic cold sources. The three highest Planck bands (857, 545, 353 GHz) have been combined with IRAS data at 3 THz to perform a multi-frequency detection of sources colder than their local environment. After rejection of possible extragalactic contaminants, the PGCC catalogue contains 13188 Galactic sources spread across the whole sky, i.e., from the Galactic plane to high latitudes, following the spatial distribution of the main molecular cloud complexes. The median temperature of PGCC sources lies between 13 and 14.5 K, depending on the quality of the flux density measurements, with a temperature ranging from 5.8 to 20 K after removing sources with the 1% largest temperature estimates. Using seven independent methods, reliable distance estimates have been obtained for 5574 sources, which allows us to derive their physical properties such as their mass, physical size, mean density and luminosity. The PGCC sources are located mainly in the solar neighbourhood, up to a distance of 10.5 kpc towards the Galactic centre, and range from low-mass cores to large molecular clouds. Because of this diversity and because the PGCC catalogue contains sources in very different environments, the catalogue is useful to investigate the evolution from molecular clouds to cores. Finally, the catalogue also includes 54 additional sources located in the SMC and LMC.

preprint2015arXiv

Planck intermediate results. XL. The Sunyaev-Zeldovich signal from the Virgo cluster

The Virgo cluster is the largest Sunyaev-Zeldovich (SZ) source in the sky, both in terms of angular size and total integrated flux. Planck's wide angular scale and frequency coverage, together with its high sensitivity, allow a detailed study of this large object through the SZ effect. Virgo is well resolved by Planck, showing an elongated structure, which correlates well with the morphology observed from X-rays, but extends beyond the observed X-ray signal. We find a good agreement between the SZ signal (or Compton paranmeter, y_c) observed by Planck and the expected signal inferred from X-ray observations and simple analytical models. Due to its proximity to us, the gas beyond the virial radius can be studied with unprecedented sensitivity by integrating the SZ signal over tens of square degrees. We study the signal in the outskirts of Virgo and compare it with analytical models and a constrained simulation of the environment of Virgo. Planck data suggest that significant amounts of low-density plasma surround Virgo out to twice the virial radius. We find the SZ signal in the outskirts of Virgo to be consistent with a simple model that extrapolates the inferred pressure at lower radii while assuming that the temperature stays in the keV range beyond the virial radius. The observed signal is also consistent with simulations and points to a shallow pressure profile in the outskirts of the cluster. This reservoir of gas at large radii can be linked with the hottest phase of the elusive warm/hot intergalactic medium. Taking the lack of symmetry of Virgo into account, we find that a prolate model is favoured by the combination of SZ and X-ray data, in agreement with predictions.

preprint2015arXiv

Planck intermediate results. XXV. The Andromeda Galaxy as seen by Planck

The Andromeda Galaxy (M31) is one of a few galaxies that has sufficient angular size on the sky to be resolved by the Planck satellite. Planck has detected M31 in all of its frequency bands, and has mapped out the dust emission with the High Frequency Instrument, clearly resolving multiple spiral arms 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 M31. We find that dust dominating the longer wavelength emission ($\gtrsim 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 (SEDs) 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\pm1.0$)$\,$K with a spectral index of $1.62\pm0.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$_\odot\,$yr$^{-1}$. We find a $2.3\,σ$ detection of the presence of spinning dust emission, with a 30$\,$GHz amplitude of $0.7\pm0.3\,$Jy, which is in line with expectations from our Galaxy.

preprint2015arXiv

Planck intermediate results. XXVII. High-redshift infrared galaxy overdensity candidates and lensed sources discovered by Planck and confirmed by Herschel-SPIRE

[Abridged] We use the Planck all-sky submm and mm maps to search for rare sources distinguished by extreme brightness, a few hundreds of mJy, 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 and 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 (CIB). We perform a dedicated Herschel-SPIRE follow-up of 234 such Planck targets, finding a significant excess of red 350 and 500um 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 350um. About 3% are candidate lensed systems, all 12 of which have secure spectroscopic confirmations, placing them at redshifts z>2.2. The galaxy overdensities are detected with high significance, half of the sample showing statistical significance above 10sigma. The SPIRE photometric redshifts of galaxies in overdensities suggest a peak at z~2. Under the Td=35K assumption, we derive an infrared (IR) luminosity for each SPIRE source of about 4x10^12 Lsun, yielding star formation rates of typically 700 Msun.yr^-1. If the observed overdensities are actual gravitationally-bound structures, the total total star formation rates reaches 7x10^3 Msun.yr^-1. Taken together, these sources show the signatures of high-z (z>$) protoclusters of intensively star-forming galaxies. All these observations confirm the uniqueness of our sample and demonstrate the ability of the all-sky Planck-HFI cold sources to select populations of cosmological and astrophysical interest for structure formation studies.

preprint2015arXiv

Planck intermediate results. XXVIII. Interstellar gas and dust in the Chamaeleon clouds as seen by Fermi LAT and Planck

Shortened abstract: Observations of the nearby Chamaeleon clouds in gamma rays with the Fermi Large Area Telescope and in thermal dust emission with Planck and IRAS have been used with the HI and CO radio data to (i) map the gas column densities in the different phases and at the dark neutral medium (DNM) transition between the HI-bright and CO-bright media; (ii) constrain the CO-to-$H_2$ conversion factor, $X_{CO}$; (iii) probe the dust properties per gas nucleon in each gas phase and spatially across the clouds. We have separated clouds in velocity in HI and CO emission and modelled the 0.4-100 GeV intensity, the dust optical depth at 353 GHz, the thermal radiance of the large grains, and an estimate of the dust extinction empirically corrected for the starlight intensity, $A_{VQ}$. The gamma-ray emissivity spectra confirm that the GeV-TeV cosmic rays uniformly permeate all gas phases up to the CO cores. The dust and cosmic rays 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 and CO-dark $H_2$ fractions for separate clouds. The corrected extinction provides the best fit to the total gas traced by the gamma rays, but we find evidence for a rise in $A_{VQ}/N_H$ and a steep rise in opacity, with increasing $N_H$ and $H_2$ fraction, and with decreasing dust temperature. We observe less variations for the specific power of the grains, except for a decline by half in the CO cores. This combined information suggests grain evolution. The gamma rays and dust radiance yield consistent $X_{CO}$ estimates near $0.7\times10^{20}$ cm$^{-2}$ (K km/s)$^{-1}$. The other dust tracers yield biased values because of the grain opacity rise in the CO clouds. These results also confirm a factor of 2 difference between $X_{CO}$ estimates at pc and kpc scales.

preprint2015arXiv

Planck intermediate results. XXXII. The relative orientation between the magnetic field and structures traced by interstellar dust

The role of the magnetic field in the formation of the filamentary structures observed in the interstellar medium (ISM) is a debated topic. The Planck all-sky maps of linearly polarized emission from dust at 353GHz provide the required combination of imaging and statistics to study the correlation between the structures of the Galactic magnetic field and of interstellar matter, both in the diffuse ISM and in molecular clouds. The data reveal structures, or ridges, 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 column density from $10^{20}$ to $10^{22}$ cm$^{-2}$. We measure the magnetic field orientation on the plane of the sky from the polarization data, and present an algorithm to estimate the orientation of the ridges from the dust intensity map. We use analytical models to account for projection effects. Comparing polarization angles on and off the structures, we estimate the mean ratio between the strengths of the turbulent and mean components of the magnetic field to be between 0.6 and 1.0, with a preferred value of 0.8. We find that the ridges are preferentially aligned with the magnetic field measured on the structures. This trend becomes more striking for increasing polarization fraction and decreasing column density. We interpret the increase of alignment with polarization fraction as a consequence of projections effects. The decrease of alignment for high column density is not due to a loss of correlation between the structures and the geometry of the magnetic field. In molecular complexes, we observe structures perpendicular to the magnetic field, which cannot be accounted for by projection effects. We discuss our results in the context of models and MHD simulations, which describe the formation of structures in the magnetized ISM.

preprint2015arXiv

Planck intermediate results. XXXIII. Signature of the magnetic field geometry of interstellar filaments in dust polarization maps

Planck observations at 353GHz provide the first fully-sampled maps of the polarized dust emission towards interstellar filaments and their backgrounds. The polarization data provide insight on the structure of their magnetic field (B). We present the polarization maps of three nearby star forming filament of moderate column density (NH~10^22cm^-2): Musca, B211, and L1506. We use the spatial information to separate Stokes I, Q, and U of the filaments from those of their backgrounds, an essential step to measure the intrinsic polarization fraction (p) and angle (psi) of each emission component. We find that the polarization angles in the three filaments (psi_fil) are coherent along their lengths and not the same as in their backgrounds (psi_bg). The differences between psi_fil and psi_bg are 12deg, 6deg, and 54deg for Musca, B211, and L1506, respectively. These differences for Musca and L1506 are larger than the dispersions of psi, both along the filaments and in their backgrounds. The observed changes of psi are direct evidence for variations of the orientation of the plane of the sky (POS) projection of the B-field. As in previous studies, we find a decrease of several percent of p with NH. We show that the drop in p cannot be explained by random fluctuations of the orientation of B within the filaments because they are too small (sigma_psi<10deg). We recognize the degeneracy between dust alignment efficiency and the structure of B in causing variations in p, but we argue that the decrease of p from the backgrounds to the filaments results in part from depolarization associated with the 3D structure of B: both its orientation in the POS and with respect to the POS. We do not resolve the inner structure of the filaments, but at the smallest scales accessible with Planck (~0.2pc), the observed changes of psi and p hold information on the B-field structure within filaments.

preprint2015arXiv

Planck intermediate results. XXXIV. The magnetic field structure in the Rosette Nebula

Planck has mapped the polarized dust emission over the whole sky, making it possible to trace the Galactic magnetic field structure that pervades the interstellar medium (ISM). We combine polarization data from Planck with rotation measure (RM) observations towards a massive star-forming region, the Rosette Nebula in the Monoceros molecular cloud, to study its magnetic field structure and the impact of an expanding HII region on the morphology of the field. We derive an analytical solution for the magnetic field, assumed to evolve from an initially uniform configuration following the expansion of ionized gas and the formation of a shell of swept-up ISM. From the RM data we estimate a mean value of the line-of-sight component of the magnetic field of about 3microG (towards the observer) in the Rosette Nebula, for a uniform electron density of about 12cm-3. The dust shell that surrounds the Rosette HII region is clearly observed in the Planck intensity map at 353 GHz, with a polarization signal significantly different from that of the local background when considered as a whole. The Planck observations constrain the plane-of-the-sky orientation of the magnetic field in the Rosette's parent molecular cloud to be mostly aligned with the large-scale field along the Galactic plane. The Planck data are compared with the analytical model, which predicts the mean polarization properties of a spherical and uniform dust shell for a given orientation of the field. This comparison leads to an upper limit of about 45degr on the angle between the line of sight and the magnetic field in the Rosette complex, for an assumed intrinsic dust polarization fraction of 4%. This field direction can reproduce the RM values detected in the ionized region if the magnetic field strength in the Monoceros molecular cloud is in the range 6.5--9microG.

preprint2015arXiv

Planck intermediate results. XXXIX. The Planck list of high-redshift source candidates

The Planck mission, thanks to its large frequency range and all-sky coverage, has a unique potential for systematically detecting the brightest, and rarest, submillimetre sources on the sky, including distant objects in the high-redshift Universe traced by their dust emission. A novel method, based on a component-separation procedure using a combination of Planck and IRAS data, has been applied to select the most luminous cold submm sources with spectral energy distributions peaking between 353 and 857GHz at 5' resolution. A total of 2151 Planck high-z source candidates (the PHZ) have been detected in the cleanest 26% of the sky, with flux density at 545GHz above 500mJy. Embedded in the cosmic infrared background close to the confusion limit, these high-z candidates exhibit colder colours than their surroundings, consistent with redshifts z>2, assuming a dust temperature of 35K and a spectral index of 1.5. First follow-up observations obtained from optical to submm have confirmed that this list consists of two distinct populations. A small fraction (around 3%) of the sources have been identified as strongly gravitationally lensed star-forming galaxies, which are amongst the brightest submm lensed objects (with flux density at 545GHz ranging from 350mJy up to 1Jy) at redshift 2 to 4. However, the vast majority of the PHZ sources appear as overdensities of dusty star-forming galaxies, having colours consistent with z>2, and may be considered as proto-cluster candidates. The PHZ provides an original sample, complementary to the Planck Sunyaev-Zeldovich Catalogue; by extending the population of the virialized massive galaxy clusters to a population of sources at z>1.5, the PHZ may contain the progenitors of today's clusters. Hence the PHZ opens a new window on the study of the early ages of structure formation, and the understanding of the intensively star-forming phase at high-z.

preprint2015arXiv

Planck intermediate results. XXXV. Probing the role of the magnetic field in the formation of structure in molecular clouds

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, $N_H$. The selected regions, covering several degrees in size, are analyzed 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 $N_H \approx 10^{21}$ to $10^{23}$ cm$^{-2}$, and hence they correspond to the bulk of the molecular clouds. The relative orientation is evaluated pixel by pixel and analyzed 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 $N_H$, from preferentially parallel or having no preferred orientation to preferentially 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.

preprint2015arXiv

Planck intermediate results. XXXVII. Evidence of unbound gas from the kinetic Sunyaev-Zeldovich effect

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\approx 0.1$. We estimate the pairwise momentum of the kSZ temperature fluctuations at the positions of the CGC (Central Galaxy Catalogue) samples extracted from Sloan Digital Sky Survey (DR7) data. For the foreground-cleaned 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 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 hydro 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\pm0.5)\times 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 SZ observations.

preprint2015arXiv

Planck intermediate results. XXXVIII. E- and B-modes of dust polarization from the magnetized filamentary structure of the interstellar medium

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, 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_{\ell}^{TE}/C_{\ell}^{EE}$ ratio, reported in the power spectra analysis of the Planck 353 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.

preprint2015arXiv

The ASKAP/EMU Source Finding Data Challenge

The Evolutionary Map of the Universe (EMU) is a proposed radio continuum survey of the Southern Hemisphere up to declination +30 deg., with the Australian Square Kilometre Array Pathfinder (ASKAP). EMU will use an automated source identification and measurement approach that is demonstrably optimal, to maximise the reliability, utility and robustness of the resulting radio source catalogues. As part of the process of achieving this aim, a "Data Challenge" has been conducted, providing international teams the opportunity to test a variety of source finders on a set of simulated images. The aim is to quantify the accuracy of existing automated source finding and measurement approaches, and to identify potential limitations. The Challenge attracted nine independent teams, who tested eleven different source finding tools. In addition, the Challenge initiators also tested the current ASKAPsoft source-finding tool to establish how it could benefit from incorporating successful features of the other tools. Here we present the results of the Data Challenge, identifying the successes and limitations for this broad variety of the current generation of radio source finding tools. As expected, most finders demonstrate completeness levels close to 100% at 10sigma dropping to levels around 10% by 5sigma. The reliability is typically close to 100% at 10sigma, with performance to lower sensitivities varying greatly between finders. All finders demonstrate the usual trade-off between completeness and reliability, whereby maintaining a high completeness at low signal-to-noise comes at the expense of reduced reliability, and vice-versa. We conclude with a series of recommendations for improving the performance of the ASKAPsoft source-finding tool.

preprint2014arXiv

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

The ESA's Planck satellite, dedicated to studying the early Universe and its subsequent evolution, was launched 14 May 2009 and has been scanning the microwave and submillimetre sky continuously since 12 August 2009. This paper gives an overview of the mission and its performance, the processing, analysis, and characteristics of the data, the scientific results, and the science data products and papers in the release. The science products include maps of the CMB and diffuse extragalactic foregrounds, a catalogue of compact Galactic and extragalactic sources, and a list of sources detected through the SZ effect. The likelihood code used to assess cosmological models against the Planck data and a lensing likelihood are described. Scientific results include robust support for the standard six-parameter LCDM model of cosmology and improved measurements of its parameters, including a highly significant deviation from scale invariance of the primordial power spectrum. The Planck values for these parameters and others derived from them are significantly different from those previously determined. Several large-scale anomalies in the temperature distribution of the CMB, first detected by WMAP, are confirmed with higher confidence. Planck sets new limits on the number and mass of neutrinos, and has measured gravitational lensing of CMB anisotropies at greater than 25 sigma. Planck finds no evidence for non-Gaussianity in the CMB. Planck's results agree well with results from the measurements of baryon acoustic oscillations. Planck finds a lower Hubble constant than found in some more local measures. Some tension is also present between the amplitude of matter fluctuations derived from CMB data and that derived from SZ data. The Planck and WMAP power spectra are offset from each other by an average level of about 2% around the first acoustic peak.

preprint2014arXiv

Planck 2013 results. II. Low Frequency Instrument data processing

We describe the data processing pipeline of the Planck Low Frequency Instrument (LFI) data processing centre (DPC) to create and characterize full-sky maps based on the first 15.5 months of operations at 30, 44 and 70 GHz. In particular, we discuss the various steps involved in reducing the data, starting from telemetry packets through to the production of cleaned, calibrated timelines and calibrated frequency maps. Data are continuously calibrated using the modulation induced on the mean temperature of the cosmic microwave background radiation by the proper motion of the spacecraft. Sky signals other than the dipole are removed by an iterative procedure based on simultaneous fitting of calibration parameters and sky maps. Noise properties are estimated from time-ordered data after the sky signal has been removed, using a generalized least square map-making algorithm. A destriping code (Madam) is employed to combine radiometric data and pointing information into sky maps, minimizing the variance of correlated noise. Noise covariance matrices, required to compute statistical uncertainties on LFI and Planck products, are also produced. Main beams are estimated down to the -20 dB level using Jupiter transits, which are also used for the geometrical calibration of the focal plane.

preprint2014arXiv

Planck 2013 results. III. LFI systematic uncertainties

We present the current estimate of instrumental and systematic effect uncertainties for the Planck-Low Frequency Instrument relevant to the first release of the Planck cosmological results. We give an overview of the main effects and of the tools and methods applied to assess residuals in maps and power spectra. We also present an overall budget of known systematic effect uncertainties, which are dominated sidelobe straylight pick-up and imperfect calibration. However, even these two effects are at least two orders of magnitude weaker than the cosmic microwave background (CMB) fluctuations as measured in terms of the angular temperature power spectrum. A residual signal above the noise level is present in the multipole range $\ell<20$, most notably at 30 GHz, and is likely caused by residual Galactic straylight contamination. Current analysis aims to further reduce the level of spurious signals in the data and to improve the systematic effects modelling, in particular with respect to straylight and calibration uncertainties.

preprint2014arXiv

Planck 2013 results. V. LFI calibration

We discuss the methods employed to photometrically calibrate the data acquired by the Low Frequency Instrument on Planck. Our calibration is based on a combination of the Orbital Dipole plus the Solar Dipole, caused respectively by the motion of the Planck spacecraft with respect to the Sun and by motion of the Solar System with respect to the CMB rest frame. The latter provides a signal of a few mK with the same spectrum as the CMB anisotropies and is visible throughout the mission. In this data release we rely on the characterization of the Solar Dipole as measured by WMAP. We also present preliminary results (at 44GHz only) on the study of the Orbital Dipole, which agree with the WMAP value of the Solar System speed within our uncertainties. We compute the calibration constant for each radiometer roughly once per hour, in order to keep track of changes in the detectors' gain. Since non-idealities in the optical response of the beams proved to be important, we implemented a fast convolution algorithm which considers the full beam response in estimating the signal generated by the dipole. Moreover, in order to further reduce the impact of residual systematics due to sidelobes, we estimated time variations in the calibration constant of the 30GHz radiometers (the ones with the largest sidelobes) using the signal of a reference load. We have estimated the calibration accuracy in two ways: we have run a set of simulations to assess the impact of statistical errors and systematic effects in the instrument and in the calibration procedure, and we have performed a number of consistency checks on the data and on the brightness temperature of Jupiter. Calibration errors for this data release are expected to be about 0.6% at 44 and 70 GHz, and 0.8% at 30 GHz. (Abriged.)

preprint2014arXiv

Planck 2013 results. XVII. Gravitational lensing by large-scale structure

On the arcminute angular scales probed by Planck, the CMB anisotropies are gently perturbed by gravitational lensing. Here we present a detailed study of this effect, detecting lensing independently in the 100, 143, and 217GHz frequency bands with an overall significance of greater than 25sigma. We use the temperature-gradient correlations induced by lensing to reconstruct a (noisy) map of the CMB lensing potential, which provides an integrated measure of the mass distribution back to the CMB last-scattering surface. Our lensing potential map is significantly correlated with other tracers of mass, a fact which we demonstrate using several representative tracers of large-scale structure. We estimate the power spectrum of the lensing potential, finding generally good agreement with expectations from the best-fitting LCDM model for the Planck temperature power spectrum, showing that this measurement at z=1100 correctly predicts the properties of the lower-redshift, later-time structures which source the lensing potential. When combined with the temperature power spectrum, our measurement provides degeneracy-breaking power for parameter constraints; it improves CMB-alone constraints on curvature by a factor of two and also partly breaks the degeneracy between the amplitude of the primordial perturbation power spectrum and the optical depth to reionization, allowing a measurement of the optical depth to reionization which is independent of large-scale polarization data. Discarding scale information, our measurement corresponds to a 4% constraint on the amplitude of the lensing potential power spectrum, or a 2% constraint on the RMS amplitude of matter fluctuations at z~2.

preprint2014arXiv

Planck 2013 results. XXIII. Isotropy and statistics of the CMB

The two fundamental assumptions of the standard cosmological model - that the initial fluctuations are statistically isotropic and Gaussian - are rigorously tested using maps of the cosmic microwave background (CMB) anisotropy from the Planck satellite. Deviations from isotropy have been found and demonstrated to be robust against component separation algorithm, mask choice and frequency dependence. Many of these anomalies were previously observed in the WMAP data, and are now confirmed at similar levels of significance (about 3 sigma). However, we find little evidence for non-Gaussianity, with the exception of a few statistical signatures that seem to be associated with specific anomalies. In particular, we find that the quadrupole-octopole alignment is also connected to a low observed variance of the CMB signal. A power asymmetry is now found to persist to scales corresponding to about l=600, and can be described in the low-l regime by a phenomenological dipole modulation model. However, any primordial power asymmetry is strongly scale-dependent and does not extend to arbitrarily small angular scales. Finally, it is plausible that some of these features may be reflected in the angular power spectrum of the data, which shows a deficit of power on similar scales. Indeed, when the power spectra of two hemispheres defined by a preferred direction are considered separately, one shows evidence for a deficit in power, while its opposite contains oscillations between odd and even modes that may be related to the parity violation and phase correlations also detected in the data. Although these analyses represent a step forward in building an understanding of the anomalies, a satisfactory explanation based on physically motivated models is still lacking.

preprint2014arXiv

Planck intermediate results. XXVI. Optical identification and redshifts of Planck clusters with the RTT150 telescope

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 Bolshoy Telescope Azimutal'ny (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 redhifts for 39 Planck clusters that were not included in the Planck SZ source catalogue and are published here for the first time.

preprint2014arXiv

Planck intermediate results. XXX. The angular power spectrum of polarized dust emission at intermediate and high Galactic latitudes

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 $α^{EE,BB}=-2.42\pm0.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 $β_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}/(2π)$ of $1.32\times10^{-2}μ$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.

preprint2013arXiv

Dust and star-formation properties of a complete sample of local galaxies drawn from the Planck Early Release Compact Source Catalogue

We combine Planck HFI data at 857, 545, 353 & 217GHz with data from WISE, Spitzer, IRAS & Herschel to investigate the properties of a flux limited sample of local star-forming galaxies. A 545GHz flux density limit was chosen so that the sample is 80% complete at this frequency, giving a sample of 234 local galaxies. We investigate the dust emission and star formation properties of the sample via various models & calculate the local dust mass function. Although 1-component modified black bodies fit the dust emission longward of 80um very well (median beta=1.83) the degeneracy between dust temp & beta also means that the SEDs are very well described by a dust emissivity index fixed at beta=2 and 10<T<25 K. Although a second, warmer dust component is required to fit shorter wavelength data, & contributes ~1/3 of the total infrared emission, its mass is negligible. No evidence is found for a very cold (6-10 K) dust component. The temp of the cold dust component is strongly influenced by the ratio of the star formation rate to the total dust mass. This implies, contrary to what is often assumed, that a significant fraction of even the emission from ~20 K dust is powered by ongoing star formation, whether or not the dust itself is associated with star forming clouds or `cirrus'. There is statistical evidence of a free-free contribution to the 217GHz flux densities of <20%. We find a median dust-to-stellar mass ratio of 0.0046; & that this ratio is anti-correlated with galaxy mass. There is good correlation between dust mass & atomic gas mass (median M_d/M_HI = 0.022), suggesting that galaxies that have more dust have more interstellar medium in general. Our derived dust mass function implies a mean dust mass density of the local Universe (for dust within galaxies), of 7.0+-1.4 x 10^5 M_solar/Mpc, significantly greater than that found in the most recent estimate using Herschel data.

preprint2013arXiv

Planck 2013 results. XII. Component separation

Planck has produced detailed all-sky observations over nine frequency bands between 30 and 857 GHz. These observations allow robust reconstruction of the primordial cosmic microwave background (CMB) temperature fluctuations over nearly the full sky, as well as new constraints on Galactic foregrounds. This paper describes the component separation framework adopted by Planck. We test four foreground-cleaned CMB maps derived using qualitatively different component separation algorithms. The quality of our reconstructions is evaluated through detailed simulations and internal comparisons, and shown through various tests to be internally consistent and robust for CMB power spectrum and cosmological parameter estimation up to l = 2000. The parameter constraints on LambdaCDM cosmologies derived from these maps are consistent with those presented in the cross-spectrum based Planck likelihood analysis. We choose two of the CMB maps for specific scientific goals. We also present maps and frequency spectra of the Galactic low-frequency, CO, and thermal dust emission. The component maps are found to provide a faithful representation of the sky, as evaluated by simulations. For the low-frequency component, the spectral index varies widely over the sky, ranging from about beta = -4 to -2. Considering both morphology and prior knowledge of the low frequency components, the index map allows us to associate a steep spectral index (beta < -3.2) with strong anomalous microwave emission, corresponding to a spinning dust spectrum peaking below 20 GHz, a flat index of beta > -2.3 with strong free-free emission, and intermediate values with synchrotron emission.

preprint2013arXiv

Planck 2013 results. XIX. The integrated Sachs-Wolfe effect

Based on CMB maps from the 2013 Planck Mission data release, this paper presents the detection of the ISW effect, i.e., the correlation between the CMB and large-scale evolving gravitational potentials. The significance of detection ranges from 2 to 4 sigma, depending on which method is used. We investigate three separate approaches, which cover essentially all previous studies, as well as breaking new ground. (i) Correlation of the CMB with the Planck reconstructed gravitational lensing potential (for the first time). This detection is made using the lensing-induced bispectrum; the correlation between lensing and the ISW effect has a significance close to 2.5 sigma. (ii) Cross-correlation with tracers of LSS, yielding around 3 sigma significance, based on a combination of radio (NVSS) and optical (SDSS) data. (iii) Aperture photometry on stacked CMB fields at the locations of known large-scale structures, which yields a 4 sigma signal when using a previously explored catalogue, but shows strong discrepancies in amplitude and scale compared to expectations. More recent catalogues give more moderate results, ranging from negligible to 2.5 sigma at most, but with a more consistent scale and amplitude, the latter being still slightly above what is expected from numerical simulations within LCMD. Where they can be compared, these measurements are compatible with previous work using data from WMAP, which had already mapped these scales to the limits of cosmic variance. Planck's broader frequency coverage confirms that the signal is achromatic, bolstering the case for ISW detection. As a final step we use tracers of large-scale structure to filter the CMB data, presenting maps of the ISW temperature perturbation. These results provide complementary and independent evidence for the existence of a dark energy component that governs the current accelerated expansion of the Universe.

preprint2013arXiv

Planck 2013 results. XV. CMB power spectra and likelihood

We present the Planck likelihood, a complete statistical description of the two-point correlation function of the CMB temperature fluctuations. We use this likelihood to derive the Planck CMB power spectrum over three decades in l, covering 2 <= l <= 2500. The main source of error at l <= 1500 is cosmic variance. Uncertainties in small-scale foreground modelling and instrumental noise dominate the error budget at higher l's. For l < 50, our likelihood exploits all Planck frequency channels from 30 to 353 GHz through a physically motivated Bayesian component separation technique. At l >= 50, we employ a correlated Gaussian likelihood approximation based on angular cross-spectra derived from the 100, 143 and 217 GHz channels. We validate our likelihood through an extensive suite of consistency tests, and assess the impact of residual foreground and instrumental uncertainties on cosmological parameters. We find good internal agreement among the high-l cross-spectra with residuals of a few uK^2 at l <= 1000. We compare our results with foreground-cleaned CMB maps, and with cross-spectra derived from the 70 GHz Planck map, and find broad agreement in terms of spectrum residuals and cosmological parameters. The best-fit LCDM cosmology is in excellent agreement with preliminary Planck polarisation spectra. The standard LCDM cosmology is well constrained by Planck by l <= 1500. For example, we report a 5.4 sigma deviation from n_s /= 1. Considering various extensions beyond the standard model, we find no indication of significant departures from the LCDM framework. Finally, we report a tension between the best-fit LCDM model and the low-l spectrum in the form of a power deficit of 5-10% at l <~ 40, significant at 2.5-3 sigma. We do not elaborate further on its cosmological implications, but note that this is our most puzzling finding in an otherwise remarkably consistent dataset. (Abridged)

preprint2013arXiv

Planck 2013 Results. XXIV. Constraints on primordial non-Gaussianity

The Planck nominal mission cosmic microwave background (CMB) maps yield unprecedented constraints on primordial non-Gaussianity (NG). Using three 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 fNL^local= 2.7+/-5.8, fNL^equil= -42+/-75, and fNL^ortho= -25+\-39 (68% CL statistical). NG is detected in the data; using skew-C_l statistics we find a nonzero bispectrum from residual point sources, and the ISW-lensing bispectrum at a level expected in the LambdaCDM scenario. 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 confirmed by skew-C_l, wavelet bispectrum and Minkowski functional estimators. Beyond estimates of individual shape amplitudes, we present model-independent, 3-dimensional reconstructions of the Planck CMB bispectrum and thus derive constraints on early-Universe scenarios that generate primordial NG, including general single-field models of inflation, excited initial states (non-Bunch-Davies vacua), and directionally-dependent vector models. We provide an initial survey of scale-dependent feature and resonance models. These results bound both general single-field and multi-field model parameter ranges, such as the speed of sound, c_s \geq 0.02 (95% CL), in an effective field theory parametrization, and the curvaton decay fraction r_D \geq 0.15 (95% CL). The Planck data significantly limit the viable parameter space of the ekpyrotic/cyclic scenarios. The amplitude of the 4-point function in the local model tauNL < 2800 (95% CL). These constraints represent the highest precision tests to date of physical mechanisms for the origin of cosmic structure.

preprint2013arXiv

Planck 2013 results. XXV. Searches for cosmic strings and other topological defects

Planck data have been used to provide stringent new constraints on cosmic strings and other defects. We describe forecasts of the CMB power spectrum induced by cosmic strings, calculating these from network models and simulations using line-of-sight Boltzmann solvers. We have studied Nambu-Goto cosmic strings, as well as field theory strings for which radiative effects are important, thus spanning the range of theoretical uncertainty in strings models. We have added the angular power spectrum from strings to that for a simple adiabatic model, with the extra fraction defined as $f_{10}$ at multipole $\ell=10$. This parameter has been added to the standard six parameter fit using COSMOMC with flat priors. For the Nambu-Goto string model, we have obtained a constraint on the string tension of $Gμ/c^2 < 1.5 x 10^{-7}$ and $f_{10} < 0.015$ at 95% confidence that can be improved to $Gμ/c^2 < 1.3 x 10^{-7}$ and $f_{10} < 0.010$ on inclusion of high-$\ell$ CMB data. For the abelian-Higgs field theory model we find, $Gμ_{AH}/c^2 < 3.2 x 10^{-7}$ and $f_{10} < 0.028$. The marginalized likelihoods for $f_{10}$ and in the $f_{10}$--$Ω_b h^2$ plane are also presented. We have also obtained constraints on $f_{10}$ for models with semi-local strings and global textures for which $Gμ/c^2 < 1.1 x 10^{-6}$. We have made complementarity searches for the specific non-Gaussian signatures of cosmic strings, calibrating with all-sky Planck resolution CMB maps generated from networks of post-recombination strings. We have obtained upper limits on the string tension at 95% confidence of $Gμ/c^2 < 8.8 x 10^{-7}$ using modal bispectrum estimation and $Gμ/c^2 < 7.8 x 10^{-7}$ for real space searches with Minkowski functionals. These are conservative upper bounds because only post-recombination string contributions have been included in the non-Gaussian analysis.

preprint2013arXiv

Planck 2013 results. XXVI. Background geometry and topology of the Universe

Planck CMB temperature maps allow detection of large-scale departures from homogeneity and isotropy. We search for topology with a fundamental domain nearly intersecting the last scattering surface (comoving distance $χ_r$). For most topologies studied the likelihood maximized over orientation shows some preference for multi-connected models just larger than $χ_r$. This effect is also present in simulated realizations of isotropic maps and we interpret it as the alignment of mild anisotropic correlations with chance features in a single realization; such a feature can also exist, in milder form, when the likelihood is marginalized over orientations. Thus marginalized, the limits on the radius $R_i$ of the largest sphere inscribed in a topological domain (at log-likelihood-ratio -5) are: in a flat Universe, $R_i>0.9χ_r$ for the cubic torus (cf. $R_i>0.9χ_r$ at 99% CL for a matched-circles search); $R_i>0.7χ_r$ for the chimney; $R_i>0.5χ_r$ for the slab; in a positively curved Universe, $R_i>1.0χ_r$ for the dodecahedron; $R_i>1.0χ_r$ for the truncated cube; $R_i>0.9χ_r$ for the octahedron. Similar limits apply to alternate topologies. We perform a Bayesian search for an anisotropic Bianchi VII$_h$ geometry. In a non-physical setting where the Bianchi parameters are decoupled from cosmology, Planck data favour a Bianchi component with a Bayes factor of at least 1.5 units of log-evidence: a Bianchi pattern is efficient at accounting for some large-scale anomalies in Planck data. However, the cosmological parameters are in strong disagreement with those found from CMB anisotropy data alone. In the physically motivated setting where the Bianchi parameters are fitted simultaneously with standard cosmological parameters, we find no evidence for a Bianchi VII$_h$ cosmology and constrain the vorticity of such models: $(ω/H)_0<8\times10^{-10}$ (95% CL). [Abridged]

preprint2013arXiv

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

The Planck Catalogue of Compact Sources (PCCS) is the catalogue of sources detected in the first 15 months of Planck operations, the "nominal" mission. It consists of nine single-frequency catalogues of compact sources, both Galactic and extragalactic, detected over the entire sky. The PCCS covers the frequency range 30--857\,GHz with higher sensitivity (it is 90% complete at 180 mJy in the best channel) and better angular resolution (from ~33' to ~5') than previous all-sky surveys in this frequency band. By construction its reliability is >80% and more than 65% of the sources have been detected at least in two contiguous Planck channels. In this paper we present the construction and validation of the PCCS, its contents and its statistical characterization.

preprint2012arXiv

Herschel-ATLAS: Planck sources in the Phase 1 fields

We present the results of a cross-correlation of the Planck Early Release Compact Source Catalog (ERCSC) with the catalog of Herschel-ATLAS sources detected in the Phase 1 fields, covering 134.55 deg2. There are 28 ERCSC sources detected by Planck at 857 GHz in this area. As many as 16 of them are probably high Galactic latitude cirrus; 10 additional sources can be clearly identified as bright, low-z galaxies; one further source is resolved by Herschel as two relatively bright sources; and the last is resolved into an unusual condensation of low-flux, probably high-redshift point sources, around a strongly lensed Herschel-ATLAS source at z = 3.26. Our results demonstrate that the higher sensitivity and higher angular resolution H-ATLAS maps provide essential information for the interpretation of candidate sources extracted from Planck sub-mm maps.

preprint2012arXiv

Mining the Herschel-ATLAS: submillimeter-selected blazars in equatorial fields

The Herschel Astrophysical Terahertz Large Area Survey (H-ATLAS) provides an unprecedented opportunity to search for blazars at sub-mm wavelengths. We cross-matched the FIRST radio source catalogue with the 11655 sources brighter than 35 mJy at 500μm in the \sim 135 square degrees of the sky covered by the H-ATLAS equatorial fields at 9 h and 15 h, plus half of the field at 12 h. We found that 379 of the H-ATLAS sources have a FIRST counterpart within 10 arcsec, including 8 catalogued blazars (plus one known blazar that was found at the edge of one the H-ATLAS maps). To search for additional blazar candidates we have devised new diagnostic diagrams and found that known blazars occupy a region of the log(S500μm/S350μm) vs. log(S500μm/S1.4GHz) plane separated from that of the other sub-mm sources with radio counterparts. Using this diagnostic we have selected 12 further candidates that turn out to be scattered in the (r-z) vs. (u-r) plane or in the WISE colour-colour diagram proposed by Massaro et al. (2012), where known blazars are concentrated in well defined strips. This suggests that the majority of them either are not blazars or have spectral energy distributions contaminated by their host galaxies. A significant fraction of true blazars are found to be hosted by star-forming galaxies. This finding, supported by an analysis of blazars detected in Planck 545 and 857 GHz bands, is at odds with the notion that blazar hosts are passive ellipticals and indicates that the sub-mm selection is providing a novel prospect on blazar properties. Based on an inspection of the available photometric data, including the WISE all-sky survey, the unpublished VIKING survey and new radio observations, we tentatively estimate that there are 11 blazars with synchrotron flux density S500μm > 35mJy over the considered area. This result already allows us to constrain blazar evolution models.

preprint2012arXiv

Simultaneous Planck, Swift, and Fermi observations of X-ray and gamma-ray selected blazars

We present simultaneous Planck, Swift, Fermi, and ground-based data for 105 blazars belonging to three samples with flux limits in the soft X-ray, hard X-ray, and gamma-ray bands. Our unique data set has allowed us to demonstrate that the selection method strongly influences the results, producing biases that cannot be ignored. Almost all the BL Lac objects have been detected by Fermi-LAT, whereas ~40% of the flat-spectrum radio quasars (FSRQs) in the radio, soft X-ray, and hard X-ray selected samples are still below the gamma-ray detection limit even after integrating 27 months of Fermi-LAT data. The radio to sub-mm spectral slope of blazars is quite flat up to ~70GHz, above which it steepens to <α>~-0.65. BL Lacs have significantly flatter spectra than FSRQs at higher frequencies. The distribution of the rest-frame synchrotron peak frequency (\nupS) in the SED of FSRQs is the same in all the blazar samples with <\nupS>=10^13.1 Hz, while the mean inverse-Compton peak frequency, <\nupIC>, ranges from 10^21 to 10^22 Hz. The distributions of \nupS and of \nupIC of BL Lacs are much broader and are shifted to higher energies than those of FSRQs and strongly depend on the selection method. The Compton dominance of blazars ranges from ~0.2 to ~100, with only FSRQs reaching values >3. Its distribution is broad and depends strongly on the selection method, with gamma-ray selected blazars peaking at ~7 or more, and radio-selected blazars at values ~1, thus implying that the assumption that the blazar power is dominated by high-energy emission is a selection effect. Simple SSC models cannot explain the SEDs of most of the gamma-ray detected blazars in all samples. The SED of the blazars that were not detected by Fermi-LAT may instead be consistent with SSC emission. Our data challenge the correlation between bolometric luminosity and \nupS predicted by the blazar sequence.

preprint2012arXiv

The optically unbiased GRB host (TOUGH) survey. VI. Radio observations at z<1 and consistency with typical star-forming galaxies

The objective of this paper is to determine the level of obscured star formation activity and dust attenuation in a sample of gamma-ray burst (GRB) hosts; and to test the hypothesis that GRB hosts have properties consistent with those of the general star-forming galaxy populations. We present a radio continuum survey of all z<1 GRB hosts in The Optically Unbiased GRB Host (TOUGH) sample supplemented with radio data for all (mostly pre-Swift) GRB-SN hosts discovered before October 2006. We present new radio data for 22 objects and have obtained a detection for three of them (GRB 980425, 021211, 031203; none in the TOUGH sample), increasing the number of radio-detected GRB hosts from two to five. The star formation rate (SFR) for the GRB 021211 host of ~825 Mo yr^-1, the highest ever reported for a GRB host, places it in the category of ultraluminous infrared galaxies. We found that at least 63% of GRB hosts have SFR < 100 Mo yr^-1 and at most 8% can have SFR > 500 Mo yr^-1. For the undetected hosts the mean radio flux (<35 uJy 3sigma) corresponds to an average SFR < 15 Mo yr^-1. Moreover, ~88% of the z<1 GRB hosts have ultraviolet dust attenuation A_UV < 6.7 mag (visual attenuation A_V < 3 mag). Hence we did not find evidence for large dust obscuration in a majority of GRB hosts. Finally, we found that the distributions of SFRs and A_UV of GRB hosts are consistent with those of Lyman break galaxies, Halpha emitters at similar redshifts and of galaxies from cosmological simulations. The similarity of the GRB population with other star-forming galaxies is consistent with the hypothesis that GRBs, a least at z<1, trace a large fraction of all star formation, and are therefore less biased indicators than once thought.

preprint2011arXiv

Planck early results XIV: ERCSC validation and extreme radio sources

Planck's all sky surveys at 30-857 GHz provide an unprecedented opportunity to follow the radio spectra of a large sample of extragalactic sources to frequencies 2-20 times higher than allowed by past, large area, ground-based surveys. We combine the results of the Planck Early Release Compact Source Catalog (ERCSC) with quasi-simultaneous ground-based observations, as well as archival data, at frequencies below or overlapping Planck frequency bands, to validate the astrometry and photometry of the ERCSC radio sources and study the spectral features shown in this new frequency window opened by Planck. The ERCSC source positions and flux density scales are found to be consistent with the ground-based observations. We present and discuss the spectral energy distributions (SEDs) of a sample of "extreme" radio sources to illustrate the richness of the ERCSC for the study of extragalactic radio sources. Variability is found to play a role in the unusual spectral features of some of these sources.

preprint2011arXiv

Planck Early Results XVIII: The power spectrum of cosmic infrared background anisotropies

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.

preprint2011arXiv

Planck Early Results: Statistical properties of extragalactic radio sources in the Planck Early Release Compact Source Catalogue

The data reported in Planck's Early Release Compact Source Catalogue (ERCSC) are exploited to measure the number counts (dN/dS) of extragalactic radio sources at 30, 44, 70, 100, 143 and 217 GHz. Due to the full-sky nature of the catalogue, this measurement extends to the rarest and brightest sources in the sky. At lower frequencies (30, 44, and 70 GHz) our counts are in very good agreement with estimates based on WMAP data, being somewhat deeper at 30 and 70 GHz, and somewhat shallower at 44 GHz. Planck's source counts at 143 and 217 GHz join smoothly with the fainter ones provided by the SPT and ACT surveys over small fractions of the sky. An analysis of source spectra, exploiting Planck's uniquely broad spectral coverage, finds clear evidence of a steepening of the mean spectral index above about 70 GHz. This implies that, at these frequencies, the contamination of the CMB power spectrum by radio sources below the detection limit is significantly lower than previously estimated.

preprint2011arXiv

Planck early results. III. First assessment of the Low Frequency Instrument in-flight performance

The scientific performance of the Planck Low Frequency Instrument (LFI) after one year of in-orbit operation is presented. We describe the main optical parameters and discuss photometric calibration, white noise sensitivity, and noise properties. A preliminary evaluation of the impact of the main systematic effects is presented. For each of the performance parameters, we outline the methods used to obtain them from the flight data and provide a comparison with pre-launch ground assessments, which are essentially confirmed in flight.

preprint2011arXiv

Planck Early Results. V. The Low Frequency Instrument data processing

We describe the processing of data from the Low Frequency Instrument (LFI) used in production of the Planck Early Release Compact Source Catalogue (ERCSC). In particular, we discuss the steps involved in reducing the data from telemetry packets to cleaned, calibrated, time-ordered data (TOD) and frequency maps. Data are continuously calibrated using the modulation of the temperature of the cosmic microwave background radiation induced by the motion of the spacecraft. Noise properties are estimated from TOD from which the sky signal has been removed using a generalized least square map-making algorithm. Measured 1/f noise knee-frequencies range from 100mHz at 30GHz to a few tens of mHz at 70GHz. A destriping code (Madam) is employed to combine radiometric data and pointing information into sky maps, minimizing the variance of correlated noise. Noise covariance matrices required to compute statistical uncertainties on LFI and Planck products are also produced. Main beams are estimated down to the approx -10dB level using Jupiter transits, which are also used for geometrical calibration of the focal plane.

preprint2011arXiv

Planck Early Results. VII. The Early Release Compact Source Catalog

(Abridged) A brief description of the methodology of construction, contents and usage of the Planck Early Release Compact Source Catalogue (ERCSC), including the Early Cold Cores (ECC) and the Early Sunyaev-Zeldovich (ESZ) cluster catalogue is provided. The catalogue is based on data that consist of mapping the entire sky once and 60% of the sky a second time by Planck, thereby comprising the first high sensitivity radio/submillimetre observations of the entire sky. A Monte-Carlo algorithm based on the injection and extraction of artificial sources into the Planck maps was implemented to select reliable sources among all extracted candidates such that the cumulative reliability of the catalogue is >=90%. The 10sigma photometric flux density limit of the catalogue at |b|>30 deg is 0.49, 1.0, 0.67, 0.5, 0.33, 0.28, 0.25, 0.47 and 0.82 Jy at each of the nine frequencies between 30 and 857 GHz. Sources which are up to a factor of ~2 fainter than this limit, and which are present in "clean" regions of the Galaxy where the sky background due to emission from the interstellar medium is low, are included in the ERCSC if they meet the high reliability criterion. The Planck ERCSC sources have known associations to stars with dust shells, stellar cores, radio galaxies, blazars, infrared luminous galaxies and Galactic interstellar medium features. A significant fraction of unclassified sources are also present in the catalogs. In addition, two early release catalogs that contain 915 cold molecular cloud core candidates and 189 SZ cluster candidates that have been generated using multi-frequency algorithms are presented. The entire source list, with more than 15000 unique sources, is ripe for follow-up characterisation with Herschel, ATCA, VLA, SOFIA, ALMA and other ground-based observing facilities.

preprint2011arXiv

Planck early results. XV. Spectral energy distributions and radio continuum spectra of northern extragalactic radio sources

Spectral energy distributions (SEDs) and radio continuum spectra are presented for a northern sample of 104 extragalactic radio sources, based on the Planck Early Release Compact Source Catalogue (ERCSC) and simultaneous multifrequency data. The nine Planck frequencies, from 30 to 857 GHz, are complemented by a set of simultaneous observations ranging from radio to gamma-rays. This is the first extensive frequency coverage in the radio and millimetre domains for an essentially complete sample of extragalactic radio sources, and it shows how the individual shocks, each in their own phase of development, shape the radio spectra as they move in the relativistic jet. The SEDs presented in this paper were fitted with second and third degree polynomials to estimate the frequencies of the synchrotron and inverse Compton (IC) peaks, and the spectral indices of low and high frequency radio data, including the Planck ERCSC data, were calculated. SED modelling methods are discussed, with an emphasis on proper, physical modelling of the synchrotron bump using multiplecomponents. Planck ERCSC data also suggest that the original accelerated electron energy spectrum could be much harder than commonly thought, with power-law index around 1.5 instead of the canonical 2.5. The implications of this are discussed for the acceleration mechanisms effective in blazar shock. Furthermore in many cases the Planck data indicate that gamma-ray emission must originate in the same shocks that produce the radio emission.

preprint2011arXiv

The Simultaneous Medicina-Planck Experiment (SiMPlE): data acquisition, reduction and first results

The Simultaneous Medicina-{\it Planck} Experiment (SiMPlE) is aimed at observing a selected sample of 263 extragalactic and Galactic sources with the Medicina 32-m single dish radio telescope in the same epoch as the Planck satellite observations. The data acquired with a frequency coverage down to 5 GHz, also combined with Planck at frequencies above 30 GHz, will constitute a useful reference catalogue of bright sources over the whole Northern hemisphere. Furthermore, source observations performed in different epochs and comparison with other catalogues allow the investigation of source variabilities on different timescales. In this work, we describe the sample selection, the on-going data acquisition campaign, the data reduction procedures, the developed tools, and the comparison with other data-sets. We present the data at 5 and 8.3 GHz for the SiMPlE Northern Sample consisting of 79 sources with $δ\geq 45^\circ$ selected in our catalogue and observed during the first 6 months of the project. A first analysis of their spectral behaviour and long-term variability is also presented.

preprint2010arXiv

Herschel-ATLAS: blazars in the SDP field

To investigate the poorly constrained sub-mm counts and spectral properties of blazars we searched for these in the Herschel-ATLAS (H-ATLAS) science demostration phase (SDP) survey catalog. We cross-matched 500$μ$m sources brighter than 50 mJy with the FIRST radio catalogue. We found two blazars, both previously known. Our study is among the first blind blazar searches at sub-mm wavelengths, i.e., in the spectral regime where little is still known about the blazar SEDs, but where the synchrotron peak of the most luminous blazars is expected to occur. Our early results are consistent with educated extrapolations of lower frequency counts and question indications of substantial spectral curvature downwards and of spectral upturns at mm wavelengths. One of the two blazars is identified with a Fermi/LAT $γ$-ray source and a WMAP source. The physical parameters of the two blazars are briefly discussed.These observations demonstrate that the H-ATLAS survey will provide key information about the physics of blazars and their contribution to sub-mm counts.