Source author record

H. Dole

H. Dole 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

127works
4topics
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

127 published item(s)

preprint2026arXiv

Euclid preparation. Calibrated intrinsic galaxy alignments in the Euclid Flagship simulation

Intrinsic alignments of galaxies are potentially a major contaminant of cosmological analyses of weak gravitational lensing. We construct a semi-analytic model of galaxy ellipticities and alignments in the \Euclid Flagship simulation to predict this contamination in Euclid's weak lensing observations. Galaxy shapes and orientations are determined by the corresponding properties of the host haloes in the underlying $N$-body simulation, as well as the relative positions of galaxies within their halo. Alignment strengths are moderated via stochastic misalignments, separately for central and satellite galaxies and conditional on the galaxy's redshift, luminosity, and rest-frame colour. The resulting model is calibrated against galaxy ellipticity statistics from the COSMOS Survey, selected alignment measurements based on Sloan Digital Sky Survey samples, and galaxy orientations extracted from the Horizon-AGN hydrodynamic simulation at redshift $z=1$. The best-fit model has a total of 12 alignment parameters and generally reproduces the calibration data sets well within the $1σ$ statistical uncertainties of the observations and the \flagship simulation, with notable exceptions for the most luminous sub-samples on small physical scales. The statistical power of the calibration data and the volume of the single \flagship realisation are still too small to provide informative prior ranges for intrinsic alignment amplitudes in relevant galaxy samples. As a first application, we predict that \Euclid end-of-mission tomographic weak gravitational lensing two-point statistics are modified by up to order $10\,\%$ due to intrinsic alignments.

preprint2026arXiv

Euclid preparation. Galaxy 2-point correlation function modelling in redshift space

The Euclid satellite will measure spectroscopic redshifts for tens of millions of emission-line galaxies. In the context of Stage-IV surveys, the 3-dimensional clustering of galaxies plays a key role in providing cosmological constraints. In this paper, we conduct a model comparison for the multipole moments of the galaxy 2-point correlation function (2PCF) in redshift space. We test state-of-the-art models, in particular the effective field theory of large-scale structure (EFT), one based on the velocity difference generating function (VDG$_{\infty}$), and different variants of Lagrangian perturbation theory (LPT) models, such as convolutional LPT (CLPT) and its effective-field-theory extension (CLEFT). We analyse the first three even multipoles of the 2PCF in the Flagship 1 simulation, which consists of four snapshots at $z\in\{0.9,1.2,1.5,1.8\}$. We study both template-fitting and full-shape approaches and find that with the template-fitting approach, only the VDG$_{\infty}$ model is able to reach a minimum fitting scale of $s_{\rm min}=20\,h^{-1}\,{\rm Mpc}$ at $z=0.9$ without biasing the recovered parameters. Indeed, the EFT model becomes inaccurate already at $s_{\rm min}=30\,h^{-1}\,{\rm Mpc}$. Conversely, in the full-shape analysis, the CLEFT and VDG$_{\infty}$ models perform similarly well, but only the CLEFT model can reach $s_{\rm min}=20\,h^{-1}\,{\rm Mpc}$ while the VDG$_{\infty}$ model is unbiased down to $s_{\rm min}=25\,h^{-1}\,{\rm Mpc}$ at the lowest redshift. Overall, in order to achieve the accuracy required by Euclid, non-perturbative modelling such as in the VDG$_{\infty}$ or CLEFT models should be considered. At $z=1.8$, the CLPT model is sufficient to describe the data with high figure of merit. This comparison selects baseline models that perform best in ideal conditions and sets the stage for an optimal analysis of Euclid data in configuration space.

preprint2026arXiv

Euclid preparation. Testing analytic models of galaxy intrinsic alignments in the Euclid Flagship simulation

We model intrinsic alignments (IA) in Euclid's Flagship simulation to investigate its impact on Euclid's weak lensing signal. Our IA implementation in the Flagship simulation takes into account photometric properties of galaxies as well as their dark matter host halos. We compare simulations against theory predictions, determining the parameters of two of the most widely used IA models: the Non Linear Alignment (NLA) and the Tidal Alignment and Tidal Torquing (TATT) models. We measure the amplitude of the simulated IA signal as a function of galaxy magnitude and colour in the redshift range $0.1<z<2.1$. We find that both NLA and TATT can accurately describe the IA signal in the simulation down to scales of $6$-$7 \,h^{-1}\,$Mpc. We measure alignment amplitudes for red galaxies comparable to those of the observations, with samples not used in the calibration procedure. For blue galaxies, our constraints are consistent with zero alignments in our first redshift bin $0.1 < z < 0.3$, but we detect a non-negligible signal at higher redshift, which is, however, consistent with the upper limits set by observational constraints. Additionally, several hydrodynamical simulations predict alignment for spiral galaxies, in agreement with our findings. Finally, the evolution of alignment with redshift is realistic and comparable to that determined in the observations. However, we find that the commonly adopted redshift power-law for IA fails to reproduce the simulation alignments above $z=1.1$. A significantly better agreement is obtained when a luminosity dependence is included, capturing the intrinsic luminosity evolution with redshift in magnitude-limited surveys. We conclude that the Flagship IA simulation is a useful tool for translating current IA constraints into predictions for IA contamination of Euclid-like samples.

preprint2026arXiv

Euclid: Early Release Observations -- The extended stellar component of the IC10 dwarf galaxy

We present a detailed analysis of the old, extended stellar component of the Local Group dwarf galaxy IC 10 using deep resolved-star photometry in the VIS and NISP bands of the Euclid Early Release Observations. Leveraging Euclid's unique combination of wide field of view and high spatial resolution, we trace red giant branch (RGB) stars out to $\sim$8 kpc from the galaxy centre, reaching azimuthally-averaged surface brightness levels as faint as $μ_{HE}\sim$29 mag arcsec$^{-2}$. Our analysis reveals that IC 10's stellar distribution is significantly more extended than previously thought. After correcting for foreground extinction and subtracting contamination from Milky Way stars and background galaxies, we derive a radial stellar density profile from RGB star counts. The profile shows a marked flattening beyond $\sim$5 kpc, and is best fit by a two-component (Sersic + exponential) model, yielding a total stellar mass in old (age $>$1 Gyr) stars of $M_{\star}=(6.7-8.1)\times10^8 M_{\odot}$. The origin of the outer stellar component is unclear. It may be accreted, even possibly associated with the counter-rotating HI gas in the outer regions of IC 10, or it may represent an ancient in-situ stellar halo. We tentatively detect two symmetric stellar overdensities at the edge of our imagery. These roughly align with the direction of IC 10's orbit around M31, suggesting that they may be signatures of tidal stripping. As part of our analysis, we derive a new distance to IC 10 based on the RGB tip, finding $D=762\pm 20$ kpc and the distance modulus is $(m-M)_0=24.41\pm 0.05$.

preprint2026arXiv

Euclid: Galaxy SED reconstruction in the PHZ processing function: impact on the PSF and the role of medium-band filters

Weak lensing surveys require accurate correction for the point spread function (PSF) when measuring galaxy shapes. For a diffraction-limited PSF, as arises in space-based missions, this correction depends on each galaxy SED. In the Euclid mission, galaxy SED reconstruction, a tasks of the photometric-redshift processing function (PHZ PF), relies on broad- and medium-band ancillary photometry. The limited wavelength sampling of the Euclid VIS passband and signal-to-noise ratio may affect the reconstruction accuracy and translate into biases in the weak lensing measurements. In this study, we present the methodology, which is employed in the Euclid PHZ PF, for reconstructing galaxy SEDs at 55 wavelengths, sampling the VIS passband every 10 nm, and we assess whether it fulfils the accuracy requirements imposed on the Euclid PSF model. We employ both physics- and data-driven methods, focusing on a new approach of template-based flux correction and Gaussian processes, and we introduce an SED metric whose bias propagates into PSF quadrupole moment errors. Our findings demonstrate that Gaussian processes and template fitting meet the requirements only in specific, but complementary, redshift intervals. We therefore propose a hybrid approach, which leverages both methods. This solution proves to be effective in meeting the Euclid accuracy requirements for most of the redshift range of the survey. Finally, we investigate the impact on the SED reconstruction of a new set of 16 evenly-spaced medium-band filters for the Subaru telescope, providing quasi-spectroscopic coverage of the VIS passband. This study shows promising results, ensuring accurate SED reconstruction and meeting the mission PSF requirements. This work thus provides not only the methodological foundation of galaxy SED reconstruction in the Euclid PHZ PF, but also a roadmap for future improvements using a new medium-band survey.

preprint2023arXiv

Euclid preparation: XXIII. Derivation of galaxy physical properties with deep machine learning using mock fluxes and H-band images

Next generation telescopes, like Euclid, Rubin/LSST, and Roman, will open new windows on the Universe, allowing us to infer physical properties for tens of millions of galaxies. Machine learning methods are increasingly becoming the most efficient tools to handle this enormous amount of data, because they are often faster and more accurate than traditional methods. We investigate how well redshifts, stellar masses, and star-formation rates (SFR) can be measured with deep learning algorithms for observed galaxies within data mimicking the Euclid and Rubin/LSST surveys. We find that Deep Learning Neural Networks and Convolutional Neutral Networks (CNN), which are dependent on the parameter space of the training sample, perform well in measuring the properties of these galaxies and have a better accuracy than methods based on spectral energy distribution fitting. CNNs allow the processing of multi-band magnitudes together with $H_{\scriptscriptstyle\rm E}$-band images. We find that the estimates of stellar masses improve with the use of an image, but those of redshift and SFR do not. Our best results are deriving i) the redshift within a normalised error of less than 0.15 for 99.9$\%$ of the galaxies with S/N>3 in the $H_{\scriptscriptstyle\rm E}$-band; ii) the stellar mass within a factor of two ($\sim0.3 \rm dex$) for 99.5$\%$ of the considered galaxies; iii) the SFR within a factor of two ($\sim0.3 \rm dex$) for $\sim$70$\%$ of the sample. We discuss the implications of our work for application to surveys as well as how measurements of these galaxy parameters can be improved with deep learning.

preprint2023arXiv

Euclid preparation. XXX. Performance assessment of the NISP Red-Grism through spectroscopic simulations for the Wide and Deep surveys

This work focuses on the pilot run of a simulation campaign aimed at investigating the spectroscopic capabilities of the Euclid Near-Infrared Spectrometer and Photometer (NISP), in terms of continuum and emission line detection in the context of galaxy evolutionary studies. To this purpose we constructed, emulated, and analysed the spectra of 4992 star-forming galaxies at $0.3 \leq z \leq 2.5$ using the NISP pixel-level simulator. We built the spectral library starting from public multi-wavelength galaxy catalogues, with value-added information on spectral energy distribution (SED) fitting results, and from Bruzual and Charlot (2003) stellar population templates. Rest-frame optical and near-IR nebular emission lines were included using empirical and theoretical relations. We inferred the 3.5$σ$ NISP red grism spectroscopic detection limit of the continuum measured in the $H$ band for star-forming galaxies with a median disk half-light radius of \ang{;;0.4} at magnitude $H= 19.5\pm0.2\,$AB$\,$mag for the Euclid Wide Survey and at $H = 20.8\pm0.6\,$AB$\,$mag for the Euclid Deep Survey. We found a very good agreement with the red grism emission line detection limit requirement for the Wide and Deep surveys. We characterised the effect of the galaxy shape on the detection capability of the red grism and highlighted the degradation of the quality of the extracted spectra as the disk size increases. In particular, we found that the extracted emission line signal to noise ratio (SNR) drops by $\sim\,$45$\%$ when the disk size ranges from \ang{;;0.25} to \ang{;;1}. These trends lead to a correlation between the emission line SNR and the stellar mass of the galaxy and we demonstrate the effect in a stacking analysis unveiling emission lines otherwise too faint to detect.

preprint2022arXiv

Euclid preparation: XIII. Forecasts for galaxy morphology with the Euclid Survey using Deep Generative Models

We present a machine learning framework to simulate realistic galaxies for the Euclid Survey. The proposed method combines a control on galaxy shape parameters offered by analytic models with realistic surface brightness distributions learned from real Hubble Space Telescope observations by deep generative models. We simulate a galaxy field of $0.4\,\rm{deg}^2$ as it will be seen by the Euclid visible imager VIS and show that galaxy structural parameters are recovered with similar accuracy as for pure analytic Sérsic profiles. Based on these simulations, we estimate that the Euclid Wide Survey will be able to resolve the internal morphological structure of galaxies down to a surface brightness of $22.5\,\rm{mag}\,\rm{arcsec}^{-2}$, and $24.9\,\rm{mag}\,\rm{arcsec}^{-2}$ for the Euclid Deep Survey. This corresponds to approximately $250$ million galaxies at the end of the mission and a $50\,\%$ complete sample for stellar masses above $10^{10.6}\,\rm{M}_\odot$ (resp. $10^{9.6}\,\rm{M}_\odot$) at a redshift $z\sim0.5$ for the wide (resp. deep) survey. The approach presented in this work can contribute to improving the preparation of future high-precision cosmological imaging surveys by allowing simulations to incorporate more realistic galaxies.

preprint2022arXiv

Euclid preparation: XX. The Complete Calibration of the Color-Redshift Relation survey: LBT observations and data release

The Complete Calibration of the Color-Redshift Relation survey (C3R2) is a spectroscopic programme designed to empirically calibrate the galaxy color-redshift relation to the Euclid depth (I_E=24.5), a key ingredient for the success of Stage IV dark energy projects based on weak lensing cosmology. A spectroscopic calibration sample as representative as possible of the galaxies in the Euclid weak lensing sample is being collected, selecting galaxies from a self-organizing map (SOM) representation of the galaxy color space. Here, we present the results of a near-infrared H- and K-bands spectroscopic campaign carried out using the LUCI instruments at the LBT. For a total of 251 galaxies, we present new highly-reliable redshifts in the 1.3<= z <=1.7 and 2<= z<=2.7 ranges. The newly-determined redshifts populate 49 SOM cells which previously contained no spectroscopic measurements and almost double the occupation numbers of an additional 153 SOM cells. A final optical ground-based observational effort is needed to calibrate the missing cells in particular in the redshift range 1.7<= z<=2.7 that lack spectroscopic calibration. In the end, Euclid itself will deliver telluric-free NIR spectra that can complete the calibration.

preprint2022arXiv

Molecular gas properties of Planck-selected protocluster candidates at z~1.3-3

We report on IRAM-30m/EMIR observations of 38 red Herschel sources in 18 Planck-selected protocluster candidates (PHz). We detect 40 CO lines on a total of 24 bright Herschel sources in 14 of the 18 PHz fields. The measured average redshift is <z>=2.25. We measure redshifts for multiple Herschel sources in eight PHz fields. In half of those fields we detect from two to three objects at similar redshifts, supporting the idea that they contain high-z protoclusters. The detection of sources at different redshifts demonstrates that foreground and background sources also contribute to the total sub-mm emission. We compare the properties of the molecular gas, and of the star formation activity of our sources with samples of normal star-forming galaxies (SFGs), sub-mm galaxies, and CO-detected cluster and protocluster galaxies at similar redshifts. We find that the PHz-IRAM sources are mainly normal SFGs, with only ~20% undergoing a starburst phase. The PHz-IRAM sources are characterized by star formation rates and gas masses that are, on average, eight and five times higher than those typical of normal SFGs at similar redshifts. Their dust temperatures, and depletion timescales are instead consistent with those of normal SFGs. The analysis of the CO spectral line energy distribution, available for ten sources, peaks at Jup=3 in most of the cases, implying low gas excitation. These properties imply that a significant fraction of the PHz-IRAM sources contains extended, and cold molecular gas reservoirs at low excitation, and that their star-formation is driven by secular processes. Multiplicity and moderate gravitational lensing might also play a role in producing the observed properties. We find that the highest star-forming protoclusters drawn from the large volume simulations, have similar SFRs as the PHz protoclusters, but separate out into a larger number of SFGs. (Abridged)

preprint2022arXiv

Questioning Planck-selected star-forming high-redshift galaxy protoclusters and their fate

About 2100 star-forming galaxy protocluster candidates at z=1-4 were identified at sub-mm wavelengths in the Planck all sky survey. Follow-up spectroscopic observations of a few candidates have confirmed the presence of actual galaxy overdensities with large star-formation rates. In this work, we use state-of-the-art hydrodynamical simulations to investigate whether the Planck high-z sub-mm sources are progenitors of massive clusters at z=0. To match the PHz sources with simulated halos, we select the most star-forming halos from z = 3 to z =1.3 in the TNG300 simulation. The total star formation rate (SFR) of the simulated protocluster candidates is computed from the SFR of all the galaxies within an aperture corresponding to the Planck beam size, including those along the line-of-sight. The simulations reproduce the Planck derived SFRs as the sum of both, the SFR of at least one of the most SF high-z halo, and the average contribution from SF sources along the line-of-sight. Focusing on the spectroscopically confirmed PHz protoclusters, we compare the observed properties of their galaxy members with those in the most SF simulated halos. We find a good agreement in the stellar mass and SFR distributions, and in the galaxy number counts, but the SFR-stellar mass relation of the simulated galaxies tends to be shifted to lower SFRs with respect to the observed one. Based on the estimated final masses of the simulated halos, we infer that between 63% and 72% of the Planck selected protoclusters will evolve into massive galaxy clusters at z=0. Despite contamination from star-forming galaxies along the line of sight, we confirm the efficiency of Planck to select star-forming protoclusters at Cosmic Noon with the simulations, and provide a new criterion for selecting the most massive cluster progenitors at high-z, using observables like the number of galaxy members and their SFR distribution.

preprint2021arXiv

$Euclid$ preparation: XV. Forecasting cosmological constraints for the $Euclid$ and CMB joint analysis

The combination and cross-correlation of the upcoming $Euclid$ data with cosmic microwave background (CMB) measurements is a source of great expectation since it will provide the largest lever arm of epochs, ranging from recombination to structure formation across the entire past light cone. In this work, we present forecasts for the joint analysis of $Euclid$ and CMB data on the cosmological parameters of the standard cosmological model and some of its extensions. This work expands and complements the recently published forecasts based on $Euclid$-specific probes, namely galaxy clustering, weak lensing, and their cross-correlation. With some assumptions on the specifications of current and future CMB experiments, the predicted constraints are obtained from both a standard Fisher formalism and a posterior-fitting approach based on actual CMB data. Compared to a $Euclid$-only analysis, the addition of CMB data leads to a substantial impact on constraints for all cosmological parameters of the standard $Λ$-cold-dark-matter model, with improvements reaching up to a factor of ten. For the parameters of extended models, which include a redshift-dependent dark energy equation of state, non-zero curvature, and a phenomenological modification of gravity, improvements can be of the order of two to three, reaching higher than ten in some cases. The results highlight the crucial importance for cosmological constraints of the combination and cross-correlation of $Euclid$ probes with CMB data.

preprint2021arXiv

Euclid preparation: I. The Euclid Wide Survey

Euclid is an ESA mission designed to constrain the properties of dark energy and gravity via weak gravitational lensing and galaxy clustering. It will carry out a wide area imaging and spectroscopy survey (EWS) in visible and near-infrared, covering roughly 15,000 square degrees of extragalactic sky on six years. The wide-field telescope and instruments are optimized for pristine PSF and reduced straylight, producing very crisp images. This paper presents the building of the Euclid reference survey: the sequence of pointings of EWS, Deep fields, Auxiliary fields for calibrations, and spacecraft movements followed by Euclid as it operates in a step-and-stare mode from its orbit around the Lagrange point L2. Each EWS pointing has four dithered frames; we simulate the dither pattern at pixel level to analyse the effective coverage. We use up-to-date models for the sky background to define the Euclid region-of-interest (RoI). The building of the reference survey is highly constrained from calibration cadences, spacecraft constraints and background levels; synergies with ground-based coverage are also considered. Via purposely-built software optimized to prioritize best sky areas, produce a compact coverage, and ensure thermal stability, we generate a schedule for the Auxiliary and Deep fields observations and schedule the RoI with EWS transit observations. The resulting reference survey RSD_2021A fulfills all constraints and is a good proxy for the final solution. Its wide survey covers 14,500 square degrees. The limiting AB magnitudes ($5σ$ point-like source) achieved in its footprint are estimated to be 26.2 (visible) and 24.5 (near-infrared); for spectroscopy, the H$_α$ line flux limit is $2\times 10^{-16}$ erg cm$^{-2}$ s$^{-1}$ at 1600 nm; and for diffuse emission the surface brightness limits are 29.8 (visible) and 28.4 (near-infrared) mag arcsec$^{-2}$.

preprint2016arXiv

Comparison of absolute gain photometric calibration between Planck/HFI and Herschel/SPIRE at 545 and 857 GHz

We compare the absolute gain photometric calibration of the Planck/HFI and Herschel/SPIRE instruments on diffuse emission. The absolute calibration of HFI and SPIRE each relies on planet flux measurements and comparison with theoretical far-infrared emission models of planetary atmospheres. We measure the photometric cross calibration between the instruments at two overlapping bands, 545 GHz / 500 $μ$m and 857 GHz / 350 $μ$m. The SPIRE maps used have been processed in the Herschel Interactive Processing Environment (Version 12) and the HFI data are from the 2015 Public Data Release 2. For our study we used 15 large fields observed with SPIRE, which cover a total of about 120 deg^2. We have selected these fields carefully to provide high signal-to-noise ratio, avoid residual systematics in the SPIRE maps, and span a wide range of surface brightness. The HFI maps are bandpass-corrected to match the emission observed by the SPIRE bandpasses. The SPIRE maps are convolved to match the HFI beam and put on a common pixel grid. We measure the cross-calibration relative gain between the instruments using two methods in each field, pixel-to-pixel correlation and angular power spectrum measurements. The SPIRE / HFI relative gains are 1.047 ($\pm$ 0.0069) and 1.003 ($\pm$ 0.0080) at 545 and 857 GHz, respectively, indicating very good agreement between the instruments. These relative gains deviate from unity by much less than the uncertainty of the absolute extended emission calibration, which is about 6.4% and 9.5% for HFI and SPIRE, respectively, but the deviations are comparable to the values 1.4% and 5.5% for HFI and SPIRE if the uncertainty from models of the common calibrator can be discounted. Of the 5.5% uncertainty for SPIRE, 4% arises from the uncertainty of the effective beam solid angle, which impacts the adopted SPIRE point source to extended source unit conversion factor (Abridged)

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. 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. XLIII. The spectral energy distribution of dust in clusters of galaxies

Although infrared (IR) overall dust emission from clusters of galaxies has been statistically detected using data from the Infrared Astronomical Satellite (IRAS), it has not been possible to sample the spectral energy distribution (SED) of this emission over its peak, and thus to break the degeneracy between dust temperature and mass. By complementing the IRAS spectral coverage with Planck satellite data from 100 to 857 GHz, we provide new constraints on the IR spectrum of thermal dust emission in clusters of galaxies. We achieve this by using a stacking approach for a sample of several hundred objects from the Planck cluster sample; this procedure averages out fluctuations from the IR sky, allowing us to reach a significant detection of the faint cluster contribution. We also use the large frequency range probed by Planck, together with component-separation techniques, to remove the contamination from both cosmic microwave background anisotropies and the thermal Sunyaev-Zeldovich effect (tSZ) signal, which dominate below 353 GHz. By excluding dominant spurious signals or systematic effects, averaged detections are reported at frequencies between 353 and 5000 GHz. We confirm the presence of dust in clusters of galaxies at low and intermediate redshifts, yielding an SED with a shape similar to that of the Milky Way. Planck's beam does not allow us to investigate the detailed spatial distribution of this emission (e.g., whether it comes from intergalactic dust or simply the dust content of the cluster galaxies), but the radial distribution of the emission appears to follow that of the stacked SZ signal, and thus the extent of the clusters. The recovered SED allows us to constrain the dust mass responsible for the signal, as well as its temperature. We additionally explore the evolution of the IR emission as a function of cluster mass and redshift.

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

GOODS-$Herschel$: identification of the individual galaxies responsible for the 80-290$μ$m cosmic infrared background

We propose a new method of pushing $Herschel$ to its faintest detection limits using universal trends in the redshift evolution of the far infrared over 24$μ$m colours in the well-sampled GOODS-North field. An extension to other fields with less multi-wavelength information is presented. This method is applied here to raise the contribution of individually detected $Herschel$ sources to the cosmic infrared background (CIRB) by a factor 5 close to its peak at 250$μ$m and more than 3 in the 350$μ$m and 500$μ$m bands. We produce realistic mock $Herschel$ images of the deep PACS and SPIRE images of the GOODS-North field from the GOODS-$Herschel$ Key Program and use them to quantify the confusion noise at the position of individual sources, i.e., estimate a "local confusion noise". Two methods are used to identify sources with reliable photometric accuracy extracted using 24$μ$m prior positions. The clean index (CI), previously defined but validated here with simulations, which measures the presence of bright 24$μ$m neighbours and the photometric accuracy index (PAI) directly extracted from the mock $Herschel$ images. After correction for completeness, thanks to our mock $Herschel$ images, individually detected sources make up as much as 54% and 60% of the CIRB in the PACS bands down to 1.1 mJy at 100$μ$m and 2.2 mJy at 160$μ$m and 55, 33, and 13% of the CIRB in the SPIRE bands down to 2.5, 5, and 9 mJy at 250$μ$m, 350$μ$m, and 500$μ$m, respectively. The latter depths improve the detection limits of $Herschel$ by factors of 5 at 250$μ$m, and 3 at 350$μ$m and 500$μ$m as compared to the standard confusion limit. Interestingly, the dominant contributors to the CIRB in all $Herschel$ bands appear to be distant siblings of the Milky Way ($z$$\sim$0.96 for $λ$$<$300$μ$m) with a stellar mass of $M_{\star}$$\sim$9$\times$10$^{10}$M$_{\odot}$.

preprint2015arXiv

Planck 2013 results. XXII. Constraints on inflation

We analyse the implications of the Planck data for cosmic inflation. The Planck nominal mission temperature anisotropy measurements, combined with the WMAP large-angle polarization, constrain the scalar spectral index to $n_s = 0.9603 \pm 0.0073$, ruling out exact scale invariance at over 5 $σ$. Planck establishes an upper bound on the tensor-to-scalar ratio of r < 0.11 (95% CL). The Planck data thus shrink the space of allowed standard inflationary models, preferring potentials with V" < 0. Exponential potential models, the simplest hybrid inflationary models, and monomial potential models of degree n > 2 do not provide a good fit to the data. Planck does not find statistically significant running of the scalar spectral index, obtaining $d n_s/d ln k = -0.0134 \pm 0.0090$. Several analyses dropping the slow-roll approximation are carried out, including detailed model comparison and inflationary potential reconstruction. We also investigate whether the primordial power spectrum contains any features. We find that models with a parameterized oscillatory feature improve the fit $χ^2$ by ~ 10; however, Bayesian evidence does not prefer these models. We constrain several single-field inflation models with generalized Lagrangians by combining power spectrum data with bounds on $f_\mathrm{NL}$ measured by Planck. The fractional primordial contribution of CDM isocurvature modes in the curvaton and axion scenarios has upper bounds of 0.25% or 3.9% (95% CL), respectively. In models with arbitrarily correlated CDM or neutrino isocurvature modes, an anticorrelation can improve $χ^2$ by approximatively 4 as a result of slightly lowering the theoretical prediction for the $\ell<40$ multipoles relative to the higher multipoles. Nonetheless, the data are consistent with adiabatic initial conditions.

preprint2015arXiv

Planck 2013 results. XXIX. The Planck catalogue of Sunyaev-Zeldovich sources: Addendum

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

preprint2015arXiv

Planck 2013 results. XXXI. Consistency of the Planck data

The Planck design and scanning strategy provide many levels of redundancy that can be exploited to provide tests of internal consistency. One of the most important is the comparison of the 70GHz and 100GHz channels. Based on different instrument technologies, with feeds located differently in the focal plane, analysed independently by different teams using different software, and near the minimum of diffuse foreground emission, these channels are in effect two different experiments. The 143GHz channel has the lowest noise level on Planck, and is near the minimum of unresolved foreground emission. In this paper, we analyse the level of consistency achieved in the 2013 Planck data. We concentrate on comparisons between the 70/100/143GHz channel maps and power spectra, particularly over the angular scales of the first and second acoustic peaks, on maps masked for diffuse Galactic emission and for strong unresolved sources. Difference maps covering angular scales from 8deg-15arcmin are consistent with noise, and show no evidence of cosmic microwave background structure. Including small but important corrections for unresolved-source residuals, we demonstrate agreement between 70 and 100GHz power spectra averaged over 70<l<390 at the 0.8% level, and agreement between 143 and 100GHz power spectra of 0.4% over the same l range. These values are within and consistent with the overall uncertainties in calibration given in the Planck 2013 results. We also present results based on the 2013 likelihood analysis showing consistency at the 0.35% between the 100/143/217GHz power spectra. We analyse calibration procedures and beams to determine what fraction of these differences can be accounted for by known approximations or systematic errors that could be controlled even better in the future, reducing uncertainties still further. Several possible small improvements are described...(abridged)

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

Planck's Dusty GEMS: Gravitationally lensed high-redshift galaxies discovered with the Planck survey

We present an analysis of 11 bright far-IR/submm sources discovered through a combination of the Planck survey and follow-up Herschel-SPIRE imaging. Each source has a redshift z=2.2-3.6 obtained through a blind redshift search with EMIR at the IRAM 30-m telescope. Interferometry obtained at IRAM and the SMA, and optical/near-infrared imaging obtained at the CFHT and the VLT reveal morphologies consistent with strongly gravitationally lensed sources. Additional photometry was obtained with JCMT/SCUBA-2 and IRAM/GISMO at 850 um and 2 mm, respectively. All objects are bright, isolated point sources in the 18 arcsec beam of SPIRE at 250 um, with spectral energy distributions peaking either near the 350 um or the 500 um bands of SPIRE, and with apparent far-infrared luminosities of up to 3x10^14 L_sun. Their morphologies and sizes, CO line widths and luminosities, dust temperatures, and far-infrared luminosities provide additional empirical evidence that these are strongly gravitationally lensed high-redshift galaxies. We discuss their dust masses and temperatures, and use additional WISE 22-um photometry and template fitting to rule out a significant contribution of AGN heating to the total infrared luminosity. Six sources are detected in FIRST at 1.4 GHz. Four have flux densities brighter than expected from the local far-infrared-radio correlation, but in the range previously found for high-z submm galaxies, one has a deficit of FIR emission, and 6 are consistent with the local correlation. The global dust-to-gas ratios and star-formation efficiencies of our sources are predominantly in the range expected from massive, metal-rich, intense, high-redshift starbursts. An extensive multi-wavelength follow-up programme is being carried out to further characterize these sources and the intense star-formation within them.

preprint2015arXiv

Probing star formation in the dense environments of z~1 lensing halos aligned with dusty star-forming galaxies detected with the South Pole Telescope

We probe star formation in the environments of massive $\sim10^{13}\,M_{\odot}$ dark matter halos at redshifts of $z$$\sim$$1$. This star formation is linked to a sub-millimetre clustering signal which we detect in maps of the Planck High Frequency Instrument that are stacked at the positions of a sample of high-redshift ($z$$>$$2$) strongly-lensed dusty star-forming galaxies (DSFGs) selected from the South Pole Telescope (SPT) 2500 deg$^2$ survey. The clustering signal has sub-millimetre colours which are consistent with the mean redshift of the foreground lensing halos ($z$$\sim$$1$). We report a mean excess of star formation rate (SFR) compared to the field, of $(2700\pm700)\,M_{\odot}\,{yr}^{-1}$ from all galaxies contributing to this clustering signal within a radius of 3.5' from the SPT DSFGs. The magnitude of the Planck excess is in broad agreement with predictions of a current model of the cosmic infrared background. The model predicts that 80$\%$ of the excess emission measured by Planck originates from galaxies lying in the neighbouring halos of the lensing halo. Using Herschel maps of the same fields, we find a clear excess, relative to the field, of individual sources which contribute to the Planck excess. The mean excess SFR compared to the field is measured to be ($370\pm40)$$\,M_{\odot}\,{yr}^{-1}$ per resolved, clustered source. Our findings suggest that the environments around these massive $z$$\sim$$1$ lensing halos host intense star formation out to about $2\,$Mpc. The flux enhancement due to clustering should also be considered when measuring flux densities of galaxies in Planck data.

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. 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). Knowledge of the beam profiles is necessary for determining the transfer function to go from the observed to the actual sky anisotropy power spectrum. The main beam distortions affect the beam window function, complicating the reconstruction of the anisotropy power spectrum at high multipoles, whereas the sidelobes affect the low and intermediate multipoles. The in-flight assessment of the LFI main beams relies on the measurements performed during Jupiter observations. By stacking the data from multiple Jupiter transits, the main beam profiles are measured down to -20 dB at 30 and 44 GHz, and down to -25 dB at 70 GHz. The main beam solid angles are determined to better than 0.2% at each LFI frequency band. The Planck pre-launch optical model is conveniently tuned to characterize the main beams independently of any noise effects. This approach provides an optical model whose beams fully reproduce the measurements in the main beam region, but also allows a description of the beams at power levels lower than can be achieved by the Jupiter measurements themselves. The agreement between the simulated beams and the measured beams is better than 1% at each LFI frequency band. The simulated beams are used for the computation of the window functions for the effective beams. The error budget for the window functions is estimated from both main beam and sidelobe contributions, and accounts for the radiometer bandshapes. The total uncertainties in the effective beam window functions are: 2% and 1.2% at 30 and 44 GHz, respectively (at $\ell \approx 600$), and 0.7% at 70 GHz (at $\ell \approx 1000$).

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. VI. High Frequency Instrument data processing

We describe the processing of the 531 billion raw data samples from the High Frequency Instrument (hereafter HFI), which we performed to produce six temperature maps from the first 473 days of Planck-HFI survey data. These maps provide an accurate rendition of the sky emission at 100, 143, 217, 353, 545, and 857 GHz with an angular resolution ranging from 9.7 to 4.6 arcmin. The detector noise per (effective) beam solid angle is respectively, 10, 6, 12 and 39 microKelvin in HFI four lowest frequency channel (100--353 GHz) and 13 and 14 kJy/sr for the 545 and 857 GHz channels. Using the 143 GHz channel as a reference, these two high frequency channels are intercalibrated within 5% and the 353 GHz relative calibration is at the percent level. The 100 and 217 GHz channels, which together with the 143 GHz channel determine the high-multipole part of the CMB power spectrum (50 < l <2500), are intercalibrated at better than 0.2 %.

preprint2014arXiv

Planck 2013 results. VIII. HFI photometric calibration and mapmaking

This paper describes the processing applied to the HFI cleaned time-ordered data to produce photometrically calibrated maps. HFI observes the sky over a broad range of frequencies, from 100 to 857 GHz. To get the best accuracy on the calibration on such a large range, two different photometric calibration schemes have to be used. The 545 and 857 \GHz\ data are calibrated using Uranus and Neptune flux density measurements, compared with models of their atmospheric emissions to calibrate the data. The lower frequencies (below 353 GHz) are calibrated using the cosmological microwave background dipole.One of the components of this anisotropy results from the orbital motion of the satellite in the Solar System, and is therefore time-variable. Photometric calibration is thus tightly linked to mapmaking, which also addresses low frequency noise removal. The 2013 released HFI data show some evidence for apparent gain variations of the HFI bolometers' detection chain. These variations were identified by comparing observations taken more than one year apart in the same configuration. We developed an effective correction to limit its effect on calibration, and assess its accuracy. We present several methods used to estimate the precision of the photometric calibration. We distinguish relative (from one detector to another, or from one frequency to another) and absolute uncertainties. In both cases, we found that these uncertainties range from a few $10^{-3}$ to several per cents from 100 to 857 GHz. We describe the pipeline producing the maps from the HFI timelines, based on the photometric calibration parameters and we detail the scheme used to a posteriori set the zero level of the maps. We also briefly discuss the cross-calibration between HFI and the SPIRE instrument on board Herschel. We finally summarize the basic characteristics of the set of the HFI maps from the 2013 Planck data release.

preprint2014arXiv

Planck 2013 results. XI. All-sky model of thermal dust emission

This paper presents an all-sky model of dust emission from the Planck 857, 545 and 353 GHz, and IRAS 100 micron data. Using a modified black-body fit to the data we present all-sky maps of the dust optical depth, temperature, and spectral index over the 353-3000 GHz range. This model is a tight representation of the data at 5 arc min. It shows variations of the order of 30 % compared with the widely-used model of Finkbeiner, Davis, and Schlegel. The Planck data allow us to estimate the dust temperature uniformly over the whole sky, providing an improved estimate of the dust optical depth compared to previous all-sky dust model, especially in high-contrast molecular regions. An increase of the dust opacity at 353 GHz, tau_353/N_H, from the diffuse to the denser interstellar medium (ISM) is reported. It is associated with a decrease in the observed dust temperature, T_obs, that could be due at least in part to the increased dust opacity. We also report an excess of dust emission at HI column densities lower than 10^20 cm^-2 that could be the signature of dust in the warm ionized medium. In the diffuse ISM at high Galactic latitude, we report an anti-correlation between tau_353/N_H and T_obs while the dust specific luminosity, i.e., the total dust emission integrated over frequency (the radiance) per hydrogen atom, stays about constant. The implication is that in the diffuse high-latitude ISM tau_353 is not as reliable a tracer of dust column density as we conclude it is in molecular clouds where the correlation of tau_353 with dust extinction estimated using colour excess measurements on stars is strong. To estimate Galactic E(B-V) in extragalactic fields at high latitude we develop a new method based on the thermal dust radiance, instead of the dust optical depth, calibrated to E(B-V) using reddening measurements of quasars deduced from Sloan Digital Sky Survey data.

preprint2014arXiv

Planck 2013 results. XIV. Zodiacal emission

The Planck satellite provides a set of all-sky maps at nine frequencies from 30 GHz to 857 GHz. Planets, minor bodies, and diffuse interplanetary dust emission (IPD) are all observed. The IPD can be separated from Galactic and other emissions because Planck views a given point on the celestial sphere multiple times, through different columns of IPD. We use the Planck data to investigate the behaviour of zodiacal emission over the whole sky at sub-millimetre and millimetre wavelengths. We fit the Planck data to find the emissivities of the various components of the COBE zodiacal model -- a diffuse cloud, three asteroidal dust bands, a circumsolar ring, and an Earth-trailing feature. The emissivity of the diffuse cloud decreases with increasing wavelength, as expected from earlier analyses. The emissivities of the dust bands, however, decrease less rapidly, indicating that the properties of the grains in the bands are different from those in the diffuse cloud. We fit the small amount of Galactic emission seen through the telescope's far sidelobes, and place limits on possible contamination of the CMB results from both zodiacal and far-sidelobe emission. When necessary, the results are used in the Planck pipeline to make maps with zodiacal emission and far sidelobes removed. We show that the zodiacal correction to the CMB maps is small compared to the Planck CMB temperature power spectrum and give a list of flux densities for small Solar System bodies.

preprint2014arXiv

Planck 2013 results. XVI. Cosmological parameters

We present the first results based on Planck measurements of the CMB temperature and lensing-potential power spectra. The Planck spectra at high multipoles are extremely well described by the standard spatially-flat six-parameter LCDM cosmology. In this model Planck data determine the cosmological parameters to high precision. We find a low value of the Hubble constant, H0=67.3+/-1.2 km/s/Mpc and a high value of the matter density parameter, Omega_m=0.315+/-0.017 (+/-1 sigma errors) in excellent agreement with constraints from baryon acoustic oscillation (BAO) surveys. Including curvature, we find that the Universe is consistent with spatial flatness to percent-level precision using Planck CMB data alone. We present results from an analysis of extensions to the standard cosmology, using astrophysical data sets in addition to Planck and high-resolution CMB data. None of these models are favoured significantly over standard LCDM. The deviation of the scalar spectral index from unity is insensitive to the addition of tensor modes and to changes in the matter content of the Universe. We find a 95% upper limit of r<0.11 on the tensor-to-scalar ratio. There is no evidence for additional neutrino-like relativistic particles. Using BAO and CMB data, we find N_eff=3.30+/-0.27 for the effective number of relativistic degrees of freedom, and an upper limit of 0.23 eV for the summed neutrino mass. Our results are in excellent agreement with big bang nucleosynthesis and the standard value of N_eff=3.046. We find no evidence for dynamical dark energy. Despite the success of the standard LCDM model, this cosmology does not provide a good fit to the CMB power spectrum at low multipoles, as noted previously by the WMAP team. While not of decisive significance, this is an anomaly in an otherwise self-consistent analysis of the Planck temperature data.

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. XX. Cosmology from Sunyaev-Zeldovich cluster counts

We present constraints on cosmological parameters using number counts as a function of redshift for a sub-sample of 189 galaxy clusters from the Planck SZ (PSZ) catalogue. The PSZ is selected through the signature of the Sunyaev--Zeldovich (SZ) effect, and the sub-sample used here has a signal-to-noise threshold of seven, with each object confirmed as a cluster and all but one with a redshift estimate. We discuss the completeness of the sample and our construction of a likelihood analysis. Using a relation between mass $M$ and SZ signal $Y$ calibrated to X-ray measurements, we derive constraints on the power spectrum amplitude $σ_8$ and matter density parameter $Ω_{\mathrm{m}}$ in a flat $Λ$CDM model. We test the robustness of our estimates and find that possible biases in the $Y$--$M$ relation and the halo mass function are larger than the statistical uncertainties from the cluster sample. Assuming the X-ray determined mass to be biased low relative to the true mass by between zero and 30%, motivated by comparison of the observed mass scaling relations to those from a set of numerical simulations, we find that $σ_8=0.75\pm 0.03$, $Ω_{\mathrm{m}}=0.29\pm 0.02$, and $σ_8(Ω_{\mathrm{m}}/0.27)^{0.3} = 0.764 \pm 0.025$. The value of $σ_8$ is degenerate with the mass bias; if the latter is fixed to a value of 20% we find $σ_8(Ω_{\mathrm{m}}/0.27)^{0.3}=0.78\pm 0.01$ and a tighter one-dimensional range $σ_8=0.77\pm 0.02$. We find that the larger values of $σ_8$ and $Ω_{\mathrm{m}}$ preferred by Planck's measurements of the primary CMB anisotropies can be accommodated by a mass bias of about 40%. Alternatively, consistency with the primary CMB constraints can be achieved by inclusion of processes that suppress power on small scales relative to the $Λ$CDM model, such as a component of massive neutrinos (abridged).

preprint2014arXiv

Planck 2013 results. XXI. All-sky Compton parameter power spectrum and high-order statistics

We have constructed the first all-sky map of the thermal Sunyaev-Zeldovich (tSZ) effect by applying specifically tailored component separation algorithms to the 100 to 857 GHz frequency channel maps from the Planck survey. These maps show an obvious galaxy cluster tSZ signal that is well matched with blindly detected clusters in the Planck SZ catalogue. To characterize the signal in the tSZ map we have computed its angular power spectrum. At large angular scales ($\ell < 60$), the major foreground contaminant is the diffuse thermal dust emission. At small angular scales ($\ell > 500$) the clustered Cosmic Infrared Background (CIB) and residual point sources are the major contaminants. These foregrounds are carefully modelled and subtracted. We measure the tSZ power spectrum in angular scales, $0.17^{\circ} \lesssim θ\lesssim 3.0^{\circ}$, that were previously unexplored. The measured tSZ power spectrum is consistent with that expected from the Planck catalogue of SZ sources, with additional clear evidence of signal from unresolved clusters and, potentially, diffuse warm baryons. We use the tSZ power spectrum to obtain the following cosmological constraints: $σ_8(Ω_{\mathrm{m}}/0.28)^{3.2/8.1}=0.784 \pm 0.016 (68% C.L.). Marginalized band-powers of the Planck tSZ power spectrum and the best-fit model are given. The non-Gaussianity of the Compton parameter map is further characterized by computing its 1D probability distribution function and its bispectrum. These are used to place additional independent constraints on $σ_{8}$.

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 2013 results. XXIX. Planck catalogue of Sunyaev-Zeldovich sources

We describe the all-sky Planck catalogue of clusters and cluster candidates derived from Sunyaev--Zeldovich (SZ) effect detections using the first 15.5 months of Planck satellite observations. The catalogue contains 1227 entries, making it over six times the size of the Planck Early SZ (ESZ) sample and the largest SZ-selected catalogue to date. It contains 861 confirmed clusters, of which 178 have been confirmed as clusters, mostly through follow-up observations, and a further 683 are previously-known clusters. The remaining 366 have the status of cluster candidates, and we divide them into three classes according to the quality of evidence that they are likely to be true clusters. The Planck SZ catalogue is the deepest all-sky cluster catalogue, with redshifts up to about one, and spans the broadest cluster mass range from (0.1 to 1.6) 10^{15}Msun. Confirmation of cluster candidates through comparison with existing surveys or cluster catalogues is extensively described, as is the statistical characterization of the catalogue in terms of completeness and statistical reliability. The outputs of the validation process are provided as additional information. This gives, in particular, an ensemble of 813 cluster redshifts, and for all these Planck clusters we also include a mass estimated from a newly-proposed SZ-mass proxy. A refined measure of the SZ Compton parameter for the clusters with X-ray counter-parts is provided, as is an X-ray flux for all the Planck clusters not previously detected in X-ray surveys.

preprint2014arXiv

Planck intermediate results. XIII. Constraints on peculiar velocities

Using \Planck\ data combined with the Meta Catalogue of X-ray detected Clusters of galaxies (MCXC), we address the study of peculiar motions by searching for evidence of the kinetic Sunyaev-Zeldovich effect (kSZ). By implementing various filters designed to extract the kSZ generated at the positions of the clusters, we obtain consistent constraints on the radial peculiar velocity average, root mean square (rms), and local bulk flow amplitude at different depths. For the whole cluster sample of average redshift 0.18, the measured average radial peculiar velocity with respect to the cosmic microwave background (CMB) radiation at that redshift, i.e., the kSZ monopole, amounts to $72 \pm 60$ km s$^{-1}$. This constitutes less than 1% of the relative Hubble velocity of the cluster sample with respect to our local CMB frame. While the linear $Λ$CDM prediction for the typical cluster radial velocity rms at $z=0.15$ is close to 230km s$^{-1}$, the upper limit imposed by \Planck\ data on the cluster subsample corresponds to 800 km s$^{-1}$ at 95% confidence level, i.e., about three times higher. \Planck\ data also set strong constraints on the local bulk flow in volumes centred on the Local Group. There is no detection of bulk flow as measured in any comoving sphere extending to the maximum redshift covered by the cluster sample. A blind search for bulk flows in this sample has an upper limit of 254 km s$^{-1}$ (95% confidence level) dominated by CMB confusion and instrumental noise, indicating that the Universe is largely homogeneous on Gpc scales. In this context, in conjunction with supernova observations, \Planck\ is able to rule out a large class of inhomogeneous void models as alternatives to dark energy or modified gravity. The \Planck\ constraints on peculiar velocities and bulk flows are thus consistent with the $Λ$CDM scenario.

preprint2014arXiv

Planck intermediate results. XIV. Dust emission at millimetre wavelengths in the Galactic plane

We use Planck HFI data combined with ancillary radio data to study the emissivity index of the interstellar dust emission in the frequency range 10 - 353 GHz, or 3 - 0.8 mm, in the Galactic plane. We analyse the region l=20 degr - 44 degr and |b| \leq 4 degr where the free-free emission can be estimated from radio recombination line data. We fit the spectra at each sky pixel with a modified blackbody model and two spectral indices, beta_mm and beta_FIR, below and above 353 GHz respectively. We find that beta_mm is smaller than beta_FIR and we detect a correlation between this low frequency power-law index and the dust optical depth at 353 GHz, tau_353. The opacity spectral index beta_mm increases from about 1.54 in the more diffuse regions of the Galactic disk, |b| = 3 degr - 4 degr and tau_353 ~ 5 x 10^{-5}, to about 1.66 in the densest regions with an optical depth of more than one order of magnitude higher. We associate this correlation with an evolution of the dust emissivity related to the fraction of molecular gas along the line of sight. This translates into beta_mm ~ 1.54 for a medium that is mostly atomic and beta_mm ~ 1.66 when the medium is dominated by molecular gas. We find that both the Two-Level System model and the emission by ferromagnetic particles can explain the results. The results improve our understanding of the physics of interstellar dust and lead towards a complete model of the dust spectrum of the Milky Way from far-infrared to millimetre wavelengths.

preprint2014arXiv

Planck intermediate results. XVII. Emission of dust in the diffuse interstellar medium from the far-infrared to microwave frequencies

The dust-HI correlation is used to characterize the emission properties of dust in the diffuse interstellar medium. We cross-correlate sky maps from Planck, WMAP, and DIRBE, at 17 frequencies from 23 to 3000 GHz, with the Parkes survey of the 21-cm line emission of neutral atomic hydrogen, over a contiguous area of 7500 deg$^2$ centred on the southern Galactic pole. Our analysis yields four specific results. (1) The dust temperature is observed to be anti-correlated with the dust emissivity and opacity. We interpret this result as evidence for dust evolution within the diffuse ISM. The mean dust opacity is measured to be $(7.1 \pm 0.6) 10^{-27} cm^2/H \times (ν/353\, GHz)^{1.53\pm0.03}$ for $100 < ν<353$GHz. (2) We map the spectral index of dust emission at millimetre wavelengths, which is remarkably constant at $β_{mm} = 1.51\pm 0.13$. We compare it with the far infrared spectral index beta_FIR derived from greybody fits at higher frequencies, and find a systematic difference, $β_{mm}-β_{FIR} = -0.15$, which suggests that the dust SED flattens at $ν< 353\,$GHz. (3) We present spectral fits of the microwave emission correlated with HI from 23 to 353 GHz, which separate dust and anomalous microwave emission. The flattening of the dust SED can be accounted for with an additional component with a blackbody spectrum, which accounts for $(26 \pm 6)$% of the dust emission at 100 GHz and could represent magnetic dipole emission. Alternatively, it could account for an increasing contribution of carbon dust, or a flattening of the emissivity of amorphous silicates, at millimetre wavelengths. These interpretations make different predictions for the dust polarization SED. (4) We identify a Galactic contribution to the residuals of the dust-HI correlation, which we model with variations of the dust emissivity on angular scales smaller than that of our correlation analysis.

preprint2014arXiv

Planck intermediate results. XXIV. Constraints on variation of fundamental constants

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

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

Planck 2013 results. IX. HFI spectral response

The Planck High Frequency Instrument (HFI) spectral response was determined through a series of ground based tests conducted with the HFI focal plane in a cryogenic environment prior to launch. The main goal of the spectral transmission tests was to measure the relative spectral response (including out-of-band signal rejection) of all HFI detectors. This was determined by measuring the output of a continuously scanned Fourier transform spectrometer coupled with all HFI detectors. As there is no on-board spectrometer within HFI, the ground-based spectral response experiments provide the definitive data set for the relative spectral calibration of the HFI. The spectral response of the HFI is used in Planck data analysis and component separation, this includes extraction of CO emission observed within Planck bands, dust emission, Sunyaev-Zeldovich sources, and intensity to polarization leakage. The HFI spectral response data have also been used to provide unit conversion and colour correction analysis tools. Verifications of the HFI spectral response data are provided through comparisons with photometric HFI flight data. This validation includes use of HFI zodiacal emission observations to demonstrate out-of-band spectral signal rejection better than 10^8. The accuracy of the HFI relative spectral response data is verified through comparison with complementary flight-data based unit conversion coefficients and colour correction coefficients. These coefficients include those based upon HFI observations of CO, dust, and Sunyaev-Zeldovich emission. General agreement is observed between the ground-based spectral characterization of HFI and corresponding in-flight observations, within the quoted uncertainty of each; explanations are provided for any discrepancies.

preprint2013arXiv

Planck 2013 results. VII. HFI time response and beams

This paper characterizes the effective beams,the effective beam window functions and the associated errors for the Planck HFI detectors. The effective beam is the angular response including the effect of the optics,detectors,data processing and the scan strategy. The window function is the representation of this beam in the harmonic domain which is required to recover an unbiased measurement of the CMB angular power spectrum. The HFI is a scanning instrument and its effective beams are the convolution of: (a) the optical response of the telescope and feeds;(b)the processing of the time-ordered data and deconvolution of the bolometric and electronic time response; and (c) the merging of several surveys to produce maps. The time response functions are measured using observations of Jupiter and Saturn and by minimizing survey difference residuals. The scanning beam is the post-deconvolution angular response of the instrument, and is characterized with observations of Mars. The main beam solid angles are determined to better than 0.5% at each HFI frequency band. Observations of Jupiter and Saturn limit near sidelobes (within 5deg) to about 0.1% of the total solid angle. Time response residuals remain as long tails in the scanning beams, but contribute less than 0.1% of the total. The bias and uncertainty in the beam products are estimated using ensembles of simulated planet observations that include the impact of instrumental noise and known systematic effects.The correlation structure of these ensembles is well-described by five error eigenmodes that are sub-dominant to sample variance and instrumental noise in the harmonic domain. A suite of consistency tests provide confidence that the error model represents a sufficient description of the data. The total error in the effective beam window functions is below 1% at 100GHz up to ell~1500$,and below 0.5% at 143 and 217GHz up to ~2000.

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. XIII. Galactic CO emission

Rotational transition lines of CO play a major role in molecular radio astronomy and in particular in the study of star formation and the Galactic structure. Although a wealth of data exists in the Galactic plane and some well-known molecular clouds, there is no available CO high sensitivity all-sky survey to date. Such all-sky surveys can be constructed using the \Planck\ HFI data because the three lowest CO rotational transition lines at 115, 230 and 345 GHz significantly contribute to the signal of the 100, 217 and 353 GHz HFI channels respectively. Two different component separation methods are used to extract the CO maps from Planck HFI data. The maps obtained are then compared to one another and to existing external CO surveys. From these quality checks the best CO maps in terms of signal to noise and/or residual foreground contamination are selected. Three sets of velocity-integrated CO emission maps are produced: Type 1 maps of the CO (1-0), (2-1), and (3-2) rotational transitions with low foreground contamination but moderate signal-to-noise ratio; Type 2 maps for the (1-0) and (2-1) transitions with a better signal-to-noise ratio; and one Type 3 map, a line composite map with the best signal-to-noise ratio in order to locate the faintest molecular regions. The maps are described in detail. They are shown to be fully compatible with previous surveys of parts of the Galactic Plane and also of fainter regions out of the Galactic plane. The Planck HFI velocity-integrated CO maps for the (1-0), (2-1), and (3-2) rotational transitions provide an unprecedented all-sky CO view of the Galaxy. These maps are also of great interest to monitor potential CO contamination on CMB \Planck\ studies.

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. XVIII. Gravitational lensing-infrared background correlation

The multi-frequency capability of the Planck satellite provides information both on the integrated history of star formation (via the cosmic infrared background, or CIB) and on the distribution of dark matter (via the lensing effect on the cosmic microwave background, or CMB). The conjunction of these two unique probes allows us to measure directly the connection between dark and luminous matter in the high redshift (1 < z <3) Universe. We use a three-point statistic optimized to detect the correlation between these two tracers. Following a thorough discussion of possible contaminants and a suite of consistency tests, using lens reconstructions at 100, 143 and 217 GHz and CIB measurements at 100-857 GHz, we report the first detection of the correlation between the CIB and CMB lensing. The well matched redshift distribution of these two signals leads to a detection significance with a peak value of 42 σat 545 GHz and a correlation as high as 80% across these two tracers. Our full set of multi-frequency measurements (both CIB auto- and CIB-lensing cross-spectra) are consistent with a simple halo-based model, with a characteristic mass scale for the halos hosting CIB sources of log_{10}(M/M_sun) = 10.5 \pm 0.6. Leveraging the frequency dependence of our signal, we isolate the high redshift contribution to the CIB, and constrain the star formation rate (SFR) density at z>1. We measure directly the SFR density with around 2 sigma significance for three redshift bins between z=1 and 7, thus opening a new window into the study of the formation of stars at early times.

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.

preprint2013arXiv

Planck 2013 results. XXX. Cosmic infrared background measurements and implications for star formation

We present new measurements of CIB anisotropies using Planck. Combining HFI data with IRAS, the angular auto- and cross frequency power spectrum is measured from 143 to 3000 GHz, and the auto-bispectrum from 217 to 545 GHz. The total areas used to compute the CIB power spectrum and bispectrum are about 2240 and 4400 deg^2, respectively. After careful removal of the contaminants, and a complete study of systematics, the CIB power spectrum and bispectrum are measured with unprecedented signal to noise ratio from angular multipoles ell~150 to 2500, and ell~130 to 1100, respectively. Two approaches are developed for modelling CIB power spectrum anisotropies. The first approach takes advantage of the unique measurements by Planck at large angular scales, and models only the linear part of the power spectrum, with a mean bias of dark matter halos hosting dusty galaxies at a given redshift weighted by their contribution to the emissivities. The second approach is based on a model that associates star-forming galaxies with dark matter halos and their subhalos, using a parametrized relation between the dust-processed infrared luminosity and (sub-)halo mass. The two approaches simultaneously fit all auto- and cross- power spectra very well. We find that the star formation history is well constrained up to z~2. However, at higher redshift, the accuracy of the star formation history measurement is strongly degraded by the uncertainty in the spectral energy distribution of CIB galaxies. We also find that CIB galaxies have warmer temperatures as redshift increases. The CIB bispectrum is steeper than that expected from the power spectrum, although well fitted by a power law; this gives some information about the contribution of massive halos to the CIB bispectrum.

preprint2013arXiv

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

(abridged for arXiv) We make use of the Planck all-sky survey to derive number counts and spectral indices of extragalactic sources -- infrared and radio sources -- from the Planck Early Catalogue (ERCSC) at 100 to 857GHz. Our sample contains, after the 80% completeness cut, between 122 and 452 and sources, with flux densities above 0.3 and 1.9Jy at 100 and 857GHz, over about 31 to 40% of the sky. Using Planck HFI, all the sources have been classified as either dust-dominated or synchrotron-dominated on the basis of their spectral energy distributions (SED). Our sample is thus complete, flux-limited and color-selected to differentiate between the two populations. We find an approximately equal number of synchrotron and dusty sources between 217 and 353GHz; at 353GHz or higher (or 217GHz and lower) frequencies, the number is dominated by dusty (synchrotron) sources, as expected. For most of the sources, the spectral indices are also derived. We provide for the first time counts of bright sources from 353 to 857GHz and the contributions from dusty and synchrotron sources at all HFI frequencies in the key spectral range where these spectra are crossing. The observed counts are in the Euclidean regime. The number counts are compared to previously published data (earlier Planck, Herschel, BLAST, SCUBA, LABOCA, SPT, and ACT) and models taking into account both radio or infrared galaxies. We derive the multi-frequency Euclidean level and compare it to WMAP, Spitzer and IRAS results. The submillimetre number counts are not well reproduced by current evolution models of dusty galaxies, whereas the millimetre part appears reasonably well fitted by the most recent model for synchrotron-dominated sources. Finally we provide estimates of the local luminosity density of dusty galaxies, providing the first such measurements at 545 and 857GHz.

preprint2013arXiv

Planck intermediate results. XII: Diffuse Galactic components in the Gould Belt System

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

preprint2013arXiv

Planck intermediate results. XVI. Profile likelihoods for cosmological parameters

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

preprint2012arXiv

A Comprehensive View of a Strongly Lensed Planck-Associated Submillimeter Galaxy

We present high-resolution maps of stars, dust, and molecular gas in a strongly lensed submillimeter galaxy (SMG) at z = 3.259. HATLAS12--00 is selected from the Herschel-Astrophysical Terahertz Large Area Survey (H-ATLAS) as a strong lens candidate mainly based on its unusually high 500um flux density (~300 mJy). It is the only high-redshift Planck detection in the 130 deg^2 H-ATLAS Phase 1 area. Keck Adaptive Optics images reveal a quadruply imaged galaxy in the K-band while the Submillimeter Array and the Extended Very Large Array show doubly imaged 880um and CO(1-0) sources, indicating differentiated distributions of the various components in the galaxy. In the source plane, the stars reside in three major kpc-scale clumps extended over ~1.6 kpc, the dust in a compact (~1 kpc) region ~3 kpc north of the stars, and the cold molecular gas in an extended (~7 kpc) disk ~5 kpc northeast of the stars. The emission from the stars, dust, and gas are magnified by ~17, 8, and 7 times, respectively, by four lensing galaxies at z ~ 1. Intrinsically, the lensed galaxy is a warm (T_dust ~ 40-65 K), hyper-luminous (L_IR ~ 1.7e13 Lsun; SFR ~ 2000 Msun/yr), gas-rich (M_gas/M_baryon ~ 70%), young (M_stellar/SFR ~ 20 Myr), and short-lived (M_gas/SFR ~ 40 Myr) starburst, without a significant active galactic nucleus. With physical properties similar to unlensed z > 2 SMGs, HATLAS12--00 offers a detailed view of a typical SMG through a powerful cosmic microscope.

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

Planck Intermediate Results II: Comparison of Sunyaev-Zeldovich measurements from Planck and from the Arcminute Microkelvin Imager for 11 galaxy clusters

A comparison is presented of Sunyaev-Zeldovich measurements for 11 galaxy clusters as obtained by Planck and by the ground-based interferometer, the Arcminute Microkelvin Imager. Assuming a universal spherically-symmetric Generalised Navarro, Frenk & White (GNFW) model for the cluster gas pressure profile, we jointly constrain the integrated Compton-Y parameter (Y_500) and the scale radius (theta_500) of each cluster. Our resulting constraints in the Y_500-theta_500 2D parameter space derived from the two instruments overlap significantly for eight of the clusters, although, overall, there is a tendency for AMI to find the Sunyaev-Zeldovich signal to be smaller in angular size and fainter than Planck. Significant discrepancies exist for the three remaining clusters in the sample, namely A1413, A1914, and the newly-discovered Planck cluster PLCKESZ G139.59+24.18. The robustness of the analysis of both the Planck and AMI data is demonstrated through the use of detailed simulations, which also discount confusion from residual point (radio) sources and from diffuse astrophysical foregrounds as possible explanations for the discrepancies found. For a subset of our cluster sample, we have investigated the dependence of our results on the assumed pressure profile by repeating the analysis adopting the best-fitting GNFW profile shape which best matches X-ray observations. Adopting the best-fitting profile shape from the X-ray data does not, in general, resolve the discrepancies found in this subset of five clusters. Though based on a small sample, our results suggest that the adopted GNFW model may not be sufficiently flexible to describe clusters universally.

preprint2012arXiv

Planck Intermediate Results. I. Further validation of new Planck clusters with XMM-Newton

We present further results from the ongoing XMM-Newton validation follow-up of Planck cluster candidates, detailing X-ray observations of eleven candidates detected at a signal-to-noise ratio of 4.5<S/N<5.3 in the same 10-month survey maps used in the construction of the Early SZ sample. The sample was selected in order to test internal SZ quality flags, and the pertinence of these flags is discussed in light of the validation results. Ten of the candidates are found to be bona fide clusters lying below the RASS flux limit. Redshift estimates are available for all confirmed systems via X-ray Fe-line spectroscopy. They lie in the redshift range 0.19<z<0.94, demonstrating Planck's capability to detect clusters up to high z. The X-ray properties of the new clusters appear to be similar to previous new detections by Planck at lower z and higher SZ flux: the majority are X-ray underluminous for their mass, estimated using Y_X as mass proxy, and many have a disturbed morphology. We find tentative indication for Malmquist bias in the Y_SZ-Y_X relation, with a turnover at Y_SZ \sim 4 e-4 arcmin^2. We present additional new optical redshift determinations with ENO and ESO telescopes of candidates previously confirmed with XMM-Newton. The X-ray and optical redshifts for a total of 20 clusters are found to be in excellent agreement. We also show that useful lower limits can be put on cluster redshifts using X-ray data only via the use of the Y_X vs. Y_SZ and X-ray flux F_X vs. Y_SZ relations.

preprint2012arXiv

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

We examine the relation between the galaxy cluster mass M and Sunyaev-Zeldovich (SZ) effect signal D_A^2 Y for a sample of 19 objects for which weak lensing (WL) mass measurements obtained from Subaru Telescope data are available in the literature. Hydrostatic X-ray masses are derived from XMM-Newton archive data and the SZ effect signal is measured from Planck all-sky survey data. We find an M_WL-D_A^2 Y relation that is consistent in slope and normalisation with previous determinations using weak lensing masses; however, there is a normalisation offset with respect to previous measures based on hydrostatic X-ray mass-proxy relations. We verify that our SZ effect measurements are in excellent agreement with previous determinations from Planck data. For the present sample, the hydrostatic X-ray masses at R_500 are on average ~ 20 per cent larger than the corresponding weak lensing masses, at odds with expectations. We show that the mass discrepancy is driven by a difference in mass concentration as measured by the two methods, and, for the present sample, the mass discrepancy and difference in mass concentration is especially large for disturbed systems. The mass discrepancy is also linked to the offset in centres used by the X-ray and weak lensing analyses, which again is most important in disturbed systems. We outline several approaches that are needed to help achieve convergence in cluster mass measurement with X-ray and weak lensing observations.

preprint2012arXiv

Planck Intermediate Results. IV. The XMM-Newton validation programme for new Planck galaxy clusters

We present the final results from the XMM-Newton validation follow-up of new Planck galaxy cluster candidates. We observed 15 new candidates, detected with signal-to-noise ratios between 4.0 and 6.1 in the 15.5-month nominal Planck survey. The candidates were selected using ancillary data flags derived from the ROSAT All Sky Survey (RASS) and Digitized Sky Survey all-sky maps, with the aim of pushing into the low SZ flux, high-z regime and testing RASS flags as indicators of candidate reliability. 14 new clusters were detected by XMM, including 2 double systems. Redshifts lie in the range 0.2 to 0.9, with 6 clusters at z>0.5. Estimated M500 range from 2.5 10^14 to 8 10^14 Msun. We discuss our results in the context of the full XMM validation programme, in which 51 new clusters have been detected. This includes 4 double and 2 triple systems, some of which are chance projections on the sky of clusters at different z. We find that association with a RASS-BSC source is a robust indicator of the reliability of a candidate, whereas association with a FSC source does not guarantee that the SZ candidate is a bona fide cluster. Nevertheless, most Planck clusters appear in RASS maps, with a significance greater than 2 sigma being a good indication that the candidate is a real cluster. The full sample gives a Planck sensitivity threshold of Y500 ~ 4 10^-4 arcmin^2, with indication for Malmquist bias in the YX-Y500 relation below this level. The corresponding mass threshold depends on z. Systems with M500 > 5 10^14 Msun at z > 0.5 are easily detectable with Planck. The newly-detected clusters follow the YX-Y500 relation derived from X-ray selected samples. Compared to X-ray selected clusters, the new SZ clusters have a lower X-ray luminosity on average for their mass. There is no indication of departure from standard self-similar evolution in the X-ray versus SZ scaling properties. (abridged)

preprint2012arXiv

Planck Intermediate Results. IX. Detection of the Galactic haze with Planck

Using precise full-sky observations from Planck, and applying several methods of component separation, we identify and characterize the emission from the Galactic "haze" at microwave wavelengths. The haze is a distinct component of diffuse Galactic emission, roughly centered on the Galactic centre, and extends to |b| ~35 deg in Galactic latitude and |l| ~15 deg in longitude. By combining the Planck data with observations from the WMAP we are able to determine the spectrum of this emission to high accuracy, unhindered by the large systematic biases present in previous analyses. The derived spectrum is consistent with power-law emission with a spectral index of -2.55 +/- 0.05, thus excluding free-free emission as the source and instead favouring hard-spectrum synchrotron radiation from an electron population with a spectrum (number density per energy) dN/dE ~ E^-2.1. At Galactic latitudes |b|<30 deg, the microwave haze morphology is consistent with that of the Fermi gamma-ray "haze" or "bubbles," indicating that we have a multi-wavelength view of a distinct component of our Galaxy. Given both the very hard spectrum and the extended nature of the emission, it is highly unlikely that the haze electrons result from supernova shocks in the Galactic disk. Instead, a new mechanism for cosmic-ray acceleration in the centre of our Galaxy is implied.

preprint2012arXiv

Planck Intermediate Results. V. Pressure profiles of galaxy clusters from the Sunyaev-Zeldovich effect

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

preprint2012arXiv

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

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

preprint2012arXiv

Planck intermediate results. VIII. Filaments between interacting clusters

About half of the baryons of the Universe are expected to be in the form of filaments of hot and low density intergalactic medium. Most of these baryons remain undetected even by the most advanced X-ray observatories which are limited in sensitivity to the diffuse low density medium. The Planck satellite has provided hundreds of detections of the hot gas in clusters of galaxies via the thermal Sunyaev-Zel'dovich (tSZ) effect and is an ideal instrument for studying extended low density media through the tSZ effect. In this paper we use the Planck data to search for signatures of a fraction of these missing baryons between pairs of galaxy clusters. Cluster pairs are good candidates for searching for the hotter and denser phase of the intergalactic medium (which is more easily observed through the SZ effect). Using an X-ray catalogue of clusters and the Planck data, we select physical pairs of clusters as candidates. Using the Planck data we construct a local map of the tSZ effect centered on each pair of galaxy clusters. ROSAT data is used to construct X-ray maps of these pairs. After having modelled and subtracted the tSZ effect and X-ray emission for each cluster in the pair we study the residuals on both the SZ and X-ray maps. For the merging cluster pair A399-A401 we observe a significant tSZ effect signal in the intercluster region beyond the virial radii of the clusters. A joint X-ray SZ analysis allows us to constrain the temperature and density of this intercluster medium. We obtain a temperature of kT = 7.1 +- 0.9, keV (consistent with previous estimates) and a baryon density of (3.7 +- 0.2)x10^-4, cm^-3. The Planck satellite mission has provided the first SZ detection of the hot and diffuse intercluster gas.

preprint2012arXiv

Planck intermediate results. X. Physics of the hot gas in the Coma cluster

We present an analysis of Planck satellite data on the Coma Cluster observed via the Sunyaev-Zeldovich effect. Planck is able, for the first time, to detect SZ emission up to r ~ 3 X R_500. We test previously proposed models for the pressure distribution in clusters against the azimuthally averaged data. We find that the Arnaud et al. universal pressure profile does not fit Coma, and that their pressure profile for merging systems provides a good fit of the data only at r<R_500: by r=2XR_500 it underestimates the observed y profile by a factor of ~2. This may indicate that at these larger radii either i) the cluster SZ emission is contaminated by unresolved SZ sources along the line of sight or ii) the pressure profile of Coma is higher at r>R_500 than the mean pressure profile predicted by the simulations. The Planck image shows significant local steepening of the y profile in two regions about half a degree to the west and to the south-east of the cluster centre. These features are consistent with the presence of shock fronts at these radii, and indeed the western feature was previously noticed in the ROSAT PSPC mosaic as well as in the radio. Using Planck y profiles extracted from corresponding sectors we find pressure jumps of 4.5+0.4-0.2 and 5.0+1.3-0.1 in the west and southeast, respectively. Assuming Rankine-Hugoniot pressure jump conditions, we deduce that the shock waves should propagate with Mach number M_w=2.03+0.09-0.04 and M_se=2.05+0.25-0.02 in the West and Southeast, respectively. Finally, we find that the y and radio-synchrotron signals are quasi-linearly correlated on Mpc scales with small intrinsic scatter. This implies either that the energy density of cosmic-ray electrons is relatively constant throughout the cluster, or that the magnetic fields fall off much more slowly with radius than previously thought.

preprint2012arXiv

Planck Intermediate Results. XI: The gas content of dark matter halos: the Sunyaev-Zeldovich-stellar mass relation for locally brightest galaxies

We present the scaling relation between Sunyaev-Zeldovich (SZ) signal and stellar mass for almost 260,000 locally brightest galaxies (LBGs) selected from the Sloan Digital Sky Survey (SDSS). These are predominantly the central galaxies of their dark matter halos. We calibrate the stellar-to-halo mass conversion using realistic mock catalogues based on the Millennium Simulation. Applying a multi-frequency matched filter to the Planck data for each LBG, and averaging the results in bins of stellar mass, we measure the mean SZ signal down to $M_\ast\sim 2\times 10^{11} \Msolar$, with a clear indication of signal at even lower stellar mass. We derive the scaling relation between SZ signal and halo mass by assigning halo properties from our mock catalogues to the real LBGs and simulating the Planck observation process. This relation shows no evidence for deviation from a power law over a halo mass range extending from rich clusters down to $M_{500}\sim 2\times 10^{13} \Msolar$, and there is a clear indication of signal down to $M_{500}\sim 4\times 10^{12} \Msolar$. Planck's SZ detections in such low-mass halos imply that about a quarter of all baryons have now been seen in the form of hot halo gas, and that this gas must be less concentrated than the dark matter in such halos in order to remain consistent with X-ray observations. At the high-mass end, the measured SZ signal is 20% lower than found from observations of X-ray clusters, a difference consistent with Malmquist bias effects in the X-ray sample.

preprint2011arXiv

GOODS-Herschel: an infrared main sequence for star-forming galaxies

We present the deepest far-IR observations obtained with Herschel and examine the 3-500um SEDs of galaxies at 0<z<2.5, supplemented by a local reference sample from IRAS, ISO, Spitzer and AKARI data. We find that the ratio of total IR luminosity to rest-frame 8um luminosity, IR8 (=Lir/L8), follows a Gaussian distribution centered on IR8=4 and defines an IR main sequence (MS). A minority population (<20 %) of outliers producing a tail skewed toward higher values of IR8 consist of starbursts (SB) with compact projected star formation densities. IR8 can be used to separate galaxies with normal and extended modes of star formation from compact SBs with high-IR8, high projected IR surface brightness (>3x10^10 Lsun kpc^-2) and a high specific SFR (i.e., SBs). The rest-frame, UV-2700A size of these distant SBs is typically half that of MS galaxies, supporting the correlation between star formation density and SB activity that is measured for the local sample. Locally, (U)LIRGs are systematically in the SB mode, whereas most distant (U)LIRGs form stars in the "normal" MS mode. This confusion between two modes of star formation is the cause of the so-called "mid-IR excess" population of galaxies found at z>1.5 by previous studies. MS galaxies have strong PAH emission line features, a broad far-IR bump resulting from a combination of dust temperatures (Tdust~15-50 K), and an effective Tdust~31 K, as derived from the peak wavelength of their IR SED. Galaxies in the SB regime instead exhibit weak PAH EW and a sharper far-IR bump with an effective Tdust~40 K. Finally, we present evidence that the mid-to-far IR emission of X-ray AGNs is predominantly produced by star formation and that candidate dusty AGNs with a power-law emission in the mid-IR systematically occur in compact, dusty SBs. After correcting for the effect of SBs on IR8, we identify new candidates for extremely obscured AGNs.

preprint2011arXiv

Planck Early Results VIII: The all-sky Early Sunyaev-Zeldovich cluster sample

We present the first all-sky sample of galaxy clusters detected blindly by the Planck satellite through the Sunyaev-Zeldovich (SZ) effect from its six highest frequencies. This early SZ (ESZ) sample is comprised of 189 candidates, which have a high signal-to-noise ratio ranging from 6 to 29. Its high reliability (purity above 95%) is further ensured by an extensive validation process based on Planck internal quality assessments and by external cross-identification and follow-up observations. Planck provides the first measured SZ signal for about 80% of the 169 previously-known ESZ clusters. Planck furthermore releases 30 new cluster candidates, amongst which 20 meet the ESZ signal-to-noise selection criterion. At the submission date, twelve of the 20 ESZ candidates were confirmed as new clusters, with eleven confirmed using XMM-Newton snapshot observations, most of them with disturbed morphologies and low luminosities. The ESZ clusters are mostly at moderate redshifts (86% with z below 0.3) and span more than a decade in mass, up to the rarest and most massive clusters with masses above 10^15 Msol.

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 XXVI: Detection with Planck and confirmation by XMM-Newton of PLCK G266.6-27.3, an exceptionally X-ray luminous and massive galaxy cluster at z~1

We present first results on PLCK G266.6-27.3, a galaxy cluster candidate detected at a signal-to-noise ratio of 5 in the Planck All Sky survey. An XMM-Newton validation observation has allowed us to confirm that the candidate is a bona fide galaxy cluster. With these X-ray data we measure an accurate redshift, z = 0.94 +/- 0.02, and estimate the cluster mass to be M_500 = (7.8 +/- 0.8)e+14 solar masses. PLCK G266.6-27.3 is an exceptional system: its luminosity of L_X(0.5-2.0 keV)=(1.4 +/- 0.05)e+45 erg/s, equals that of the two most luminous known clusters in the z > 0.5 universe, and it is one of the most massive clusters at z~1. Moreover, unlike the majority of high-redshift clusters, PLCK G266.6-27.3 appears to be highly relaxed. This observation confirms Planck's capability of detecting high-redshift, high-mass clusters, and opens the way to the systematic study of population evolution in the exponential tail of the mass function.

preprint2011arXiv

Planck early results: first assessment of the High Frequency Instrument in-flight performance

The Planck High Frequency Instrument (HFI) is designed to measure the temperature and polarization anisotropies of the Cosmic Microwave Background and galactic foregrounds in six wide bands centered at 100, 143, 217, 353, 545 and 857 GHz at an angular resolution of 10' (100 GHz), 7' (143 GHz), and 5' (217 GHz and higher). HFI has been operating flawlessly since launch on 14 May 2009. The bolometers cooled to 100 mK as planned. The settings of the readout electronics, such as the bolometer bias current, that optimize HFI's noise performance on orbit are nearly the same as the ones chosen during ground testing. Observations of Mars, Jupiter, and Saturn verified both the optical system and the time response of the detection chains. The optical beams are close to predictions from physical optics modeling. The time response of the detection chains is close to pre-launch measurements. The detectors suffer from an unexpected high flux of cosmic rays related to low solar activity. Due to the redundancy of Planck's observations strategy, the removal of a few percent of data contaminated by glitches does not affect significantly the sensitivity. The cosmic rays heat up significantly the bolometer plate and the modulation on periods of days to months of the heat load creates a common drift of all bolometer signals which do not affect the scientific capabilities. Only the high energy cosmic rays showers induce inhomogeneous heating which is a probable source of low frequency noise.

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: The Planck mission

The European Space Agency's Planck satellite was launched on 14 May 2009, and has been surveying the sky stably and continuously since 13 August 2009. Its performance is well in line with expectations, and it will continue to gather scientific data until the end of its cryogenic lifetime. We give an overview of the history of Planck in its first year of operations, and describe some of the key performance aspects of the satellite. This paper is part of a package submitted in conjunction with Planck's Early Release Compact Source Catalogue, the first data product based on Planck to be released publicly. The package describes the scientific performance of the Planck payload, and presents results on a variety of astrophysical topics related to the sources included in the Catalogue, as well as selected topics on diffuse emission.

preprint2011arXiv

Planck Early Results: The Planck View of Nearby Galaxies

The all-sky coverage of the Planck Early Release Compact Source Catalogue (ERCSC) provides an unsurpassed survey of galaxies at submillimetre (submm) wavelengths, representing a major improvement in the numbers of galaxies detected, as well as the range of far-IR/submm wavelengths over which they have been observed. We here present the first results on the properties of nearby galaxies using these data. We match the ERCSC catalogue to IRAS-detected galaxies in the Imperial IRAS Faint Source Redshift Catalogue (IIFSCz), so that we can measure the spectral energy distributions (SEDs) of these objects from 60 to 850 microns. This produces a list of 1717 galaxies with reliable associations between Planck and IRAS, from which we select a subset of 468 for SED studies, namely those with strong detections in the three highest frequency Planck bands and no evidence of cirrus contamination. The SEDs are fitted using parametric dust models to determine the range of dust temperatures and emissivities. We find evidence for colder dust than has previously been found in external galaxies, with T<20K. Such cold temperatures are found using both the standard single temperature dust model with variable emissivity beta, or a two dust temperature model with beta fixed at 2. We also compare our results to studies of distant submm galaxies (SMGs) which have been claimed to contain cooler dust than their local counterparts. We find that including our sample of 468 galaxies significantly reduces the distinction between the two populations. Fits to SEDs of selected objects using more sophisticated templates derived from radiative transfer models confirm the presence of the colder dust found through parameteric fitting. We thus conclude that cold (T<20K) dust is a significant and largely unexplored component of many nearby galaxies.

preprint2011arXiv

Planck early results. VI. The High Frequency Instrument data processing

We describe the processing of the 336 billion raw data samples from the High Frequency Instrument (HFI) which we performed to produce six temperature maps from the first 295 days of Planck-HFI survey data. These maps provide an accurate rendition of the sky emission at 100, 143, 217, 353, 545 and 857 GHz with an angular resolution ranging from 9.9 to 4.4^2. The white noise level is around 1.5 μK degree or less in the 3 main CMB channels (100--217GHz). The photometric accuracy is better than 2% at frequencies between 100 and 353 GHz and around 7% at the two highest frequencies. The maps created by the HFI Data Processing Centre reach our goals in terms of sensitivity, resolution, and photometric accuracy. They are already sufficiently accurate and well-characterised to allow scientific analyses which are presented in an accompanying series of early papers. At this stage, HFI data appears to be of high quality and we expect that with further refinements of the data processing we should be able to achieve, or exceed, the science goals of the Planck project.

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.

preprint2010arXiv

Far-Infrared Properties of Spitzer-selected Luminous Starbursts

We present SHARC-2 350 micron data on 20 luminous z ~ 2 starbursts with S(1.2mm) > 2 mJy from the Spitzer-selected samples of Lonsdale et al. and Fiolet et al. All the sources were detected, with S(350um) > 25 mJy for 18 of them. With the data, we determine precise dust temperatures and luminosities for these galaxies using both single-temperature fits and models with power-law mass--temperature distributions. We derive appropriate formulae to use when optical depths are non-negligible. Our models provide an excellent fit to the 6um--2mm measurements of local starbursts. We find characteristic single-component temperatures T1 ~ 35.5+-2.2 K and integrated infrared (IR) luminosities around 10^(12.9+-0.1) Lsun for the SWIRE-selected sources. Molecular gas masses are estimated at 4 x 10^(10) Msun, assuming kappa(850um)=0.15 m^2/kg and a submillimeter-selected galaxy (SMG)-like gas-to-dust mass ratio. The best-fit models imply >~2 kpc emission scales. We also note a tight correlation between rest-frame 1.4 GHz radio and IR luminosities confirming star formation as the predominant power source. The far-IR properties of our sample are indistinguishable from the purely submillimeter-selected populations from current surveys. We therefore conclude that our original selection criteria, based on mid-IR colors and 24 um flux densities, provides an effective means for the study of SMGs at z ~ 1.5--2.5.

preprint2010arXiv

Simulations of the cosmic infrared and submillimeter background for future large surveys: II. Removing the low-redshift contribution to the anisotropies using stacking

Herschel and Planck are surveying the sky at unprecedented angular scales and sensitivities over large areas. But both experiments are limited by source confusion in the submillimeter. The high confusion noise in particular restricts the study of the clustering properties of the sources that dominate the cosmic infrared background. At these wavelengths, it is more appropriate to consider the statistics of the unresolved component. In particular, high clustering will contribute in excess of Poisson noise in the power spectra of CIB anisotropies. These power spectra contain contributions from sources at all redshift. We show how the stacking technique can be used to separate the different redshift contributions to the power spectra. We use simulations of CIB representative of realistic Spitzer, Herschel, Planck, and SCUBA-2 observations. We stack the 24um sources in longer wavelengths maps to measure mean colors per redshift and flux bins. The information retrieved on the mean spectral energy distribution obtained with the stacking technique is then used to clean the maps, in particular to remove the contribution of low-redshift undetected sources to the anisotropies. Using the stacking, we measure the mean flux of populations 4 to 6 times fainter than the total noise at 350um at redshifts z=1 and z=2, respectively, and as faint as 6 to 10 times fainter than the total noise at 850um at the same redshifts. In the deep Spitzer fields, the detected 24um sources up to z~2 contribute significantly to the submillimeter anisotropies. We show that the method provides excellent (using COSMOS 24um data) to good (using SWIRE 24um data) removal of the z<2 (COSMOS) and z<1 (SWIRE) anisotropies. Using this cleaning method, we then hope to have a set of large maps dominated by high redshift galaxies for galaxy evolution study (e.g., clustering, luminosity density).

preprint2010arXiv

Submillimeter Number Counts at 250, 350 and 500 microns in BLAST data

BLAST (Balloon-borne Large-Aperture Submillimeter Telescope) performed the first deep and wide extragalactic survey at 250, 350 and 500 um. The extragalactic number counts at these wavelengths are important constraints for modeling the infrared galaxies evolution. [...] We use three methods to identify the submillimeter sources. 1) Blind extraction. [...] The photometry is computed with a new simple and quick PSF fitting routine (FASTPHOT). [...] 2) Extraction with prior. [...] 3) A stacking analysis. [...] With the blind extraction, we reach 97, 83 and 76 mJy at resp. 250, 350 and 500 um with a 95% completeness. With the prior extraction, we reach 76 mJy (resp. 63 and 49 mJy) at 250 um (resp. 350 and 500 um). With the stacking analysis, we reach 6.2 mJy (resp. 5.2 and 3.5 mJy) at 250 um (resp. 350 and 500 um). The differential submillimeter number counts are derived, and start showing a turnover at flux densities decreasing with increasing wavelength. There is a very good agreement with the P(D) analysis of Patanchon et al. (2009). At bright fluxes (>100 mJy), the Lagache et al. (2004) and Le Borgne et al. (2009) models slightly overestimate the observed counts, but there is a very good agreement near the peak of differential number counts. [...] Counts are available at: http://www.ias.u-psud.fr/irgalaxies/downloads.php

preprint2010arXiv

The Cosmic Far-Infrared Background Buildup Since Redshift 2 at 70 and 160 microns in the COSMOS and GOODS fields

The Cosmic Far-Infrared Background (CIB) at wavelengths around 160 μm corresponds to the peak intensity of the whole Extragalactic Background Light, which is being measured with increasing accuracy. However, the build up of the CIB emission as a function of redshift, is still not well known. Our goal is to measure the CIB history at 70 μm and 160 μm at different redshifts, and provide constraints for infrared galaxy evolution models. We use complete deep Spitzer 24 μm catalogs down to about 80 μJy, with spectroscopic and photometric redshifts identifications, from the GOODS and COSMOS deep infrared surveys covering 2 square degrees total. After cleaning the Spitzer/MIPS 70 μm and 160 μm maps from detected sources, we stacked the far-IR images at the positions of the 24 μm sources in different redshift bins. We measured the contribution of each stacked source to the total 70 and 160 μm light, and compare with model predictions and recent far-IR measurements made with Herschel/PACS on smaller fields. We have detected components of the 70 and 160 μm backgrounds in different redshift bins up to z ~ 2. The contribution to the CIB is maximum at 0.3 <= z <= 0.9 at 160μm (and z <= 0.5 at 70 μm). A total of 81% (74%) of the 70 (160) μm background was emitted at z < 1. We estimate that the AGN relative contribution to the far-IR CIB is less than about 10% at z < 1.5. We provide a comprehensive view of the CIB buildup at 24, 70, 100, 160 μm. IR galaxy models predicting a major contribution to the CIB at z < 1 are in agreement with our measurements, while our results discard other models that predict a peak of the background at higher redshifts. Our results are available online http://www.ias.u-psud.fr/irgalaxies/ .

preprint2007arXiv

The infrared luminosity function of galaxies at redshifts z=1 and z~2 in the GOODS fields

We present the rest-frame 8 micron luminosity function (LF) at redshifts z=1 and ~2, computed from Spitzer 24 micron-selected galaxies in the GOODS fields over an area of 291 sq. arcmin. Using classification criteria based on X-ray data and IRAC colours, we identify the AGN in our sample. The rest-frame 8 micron LF for star-forming galaxies at redshifts z=1 and ~2 have the same shape as at z~0, but with a strong positive luminosity evolution. The number density of star-forming galaxies with log_{10}(nu L_nu(8 micron))>11 increases by a factor >250 from redshift z~0 to 1, and is basically the same at z=1 and ~2. The resulting rest-frame 8 micron luminosity densities associated with star formation at z=1 and ~2 are more than four and two times larger than at z~0, respectively. We also compute the total rest-frame 8 micron LF for star-forming galaxies and AGN at z~2 and show that AGN dominate its bright end, which is well-described by a power-law. Using a new calibration based on Spitzer star-forming galaxies at 0<z<0.6 and validated at higher redshifts through stacking analysis, we compute the bolometric infrared (IR) LF for star-forming galaxies at z=1 and ~2. We find that the respective bolometric IR luminosity densities are (1.2+/-0.2) x 10^9 and (6.6^{+1.2}_{-1.0}) x 10^8 L_sun Mpc^{-3}, in agreement with previous studies within the error bars. At z~2, around 90% of the IR luminosity density associated with star formation is produced by luminous and ultraluminous IR galaxies (LIRG and ULIRG), with the two populations contributing in roughly similar amounts. Finally, we discuss the consistency of our findings with other existing observational results on galaxy evolution.

preprint2006arXiv

Spitzer Observations of the Brightest Galaxies in X-ray-Luminous Clusters

We have studied the infrared properties of the brightest cluster galaxies (BCGs) located in the cores of X-ray-luminous clusters at 0.15 < z < 0.35. The majority of the BCGs are not particularly infrared-luminous compared with other massive early-type galaxies, suggesting that the cluster environment has little influence on the infrared luminosities of the BCGs. The exceptions, however, are the BCGs in the three X-ray-brightest clusters in the sample, A1835, Z3146, and A2390. These BCGs have a prominent far-infrared peak in their spectral energy distributions (SEDs), and two of them (those in A1835 and Z3146) can be classified as luminous infrared galaxies (LIRGs: L_{IR} > 10^{11} L_{sun}). Although radio AGNs are found to be prevalent among the BCGs, the infrared luminosities of these three BCGs, judged from the infrared SED signatures, are likely to be powered by star formation. Considering the overall trend that clusters with shorter radiative gas cooling times harbor more infrared-luminous BCGs, the enhanced star formation may be caused by the cooling cluster gas accreting onto the BCGs.

preprint2005arXiv

Spitzer 70 and 160-micron Observations of the Extragalactic First Look Survey

We present Spitzer 70um and 160um observations of the Spitzer extragalactic First Look Survey (xFLS). The data reduction techniques and the methods for producing co-added mosaics and source catalogs are discussed. Currently, 26% of the 70um sample and 49% of the 160um-selected sources have redshifts. The majority of sources with redshifts are star-forming galaxies at z<0.5, while about 5% have infrared colors consistent with AGN. The observed infrared colors agree with the spectral energy distribution (SEDs) of local galaxies previously determined from IRAS and ISO data. The average 160um/70um color temperature for the dust is Td~= 30+/-5 K, and the average 70um/24um spectral index is alpha~= 2.4+/-0.4. The observed infrared to radio correlation varies with redshift as expected out to z~1 based on the SEDs of local galaxies. The xFLS number counts at 70um and 160um are consistent within uncertainties with the models of galaxy evolution, but there are indications that the current models may require slight modifications. Deeper 70um observations are needed to constrain the models, and redshifts for the faint sources are required to measure the evolution of the infrared luminosity function.

preprint2004arXiv

PAHs contribution to the infrared output energy of the Universe at z~2

We present an updated phenomenological galaxy evolution model to fit the Spitzer 24, 70 and 160 microns number counts as well as all the previous mid and far infrared observations. Only a minor change of the co-moving luminosity density distribution in the previous model (Lagache, Dole, Puget 2003), combined with a slight modification of the starburst template spectra mainly between 12 and 30 microns, are required to fit all the data available. We show that the peak in the MIPS 24 micron counts is dominated by galaxies with redshift between 1 and 2, with a non negligible contribution from the z>2 galaxies (~30% at S=0.2 mJy). The very close agreement between the model and number counts at 15 and 24 microns strikingly implies that (1) the PAHs (Policyclic Aromatic Hydrocarbons) features remain prominent in the redshift band 0.5 to 2.5 and (2) the IR energy output has to be dominated by ~3 10^11 Lo to ~3 10^12 Lo galaxies from redshift 0.5 to 2.5. Combining Spitzer with the Infrared Space Observatory (ISO) deep cosmological surveys gives for the first time an unbiased view of the infrared Universe from z=0 to z=2.5.

preprint2002arXiv

Modelling the infrared galaxy evolution using a phenomenological approach

To characterise the cosmological evolution of the sources contributing to the infrared extragalactic background, we have developped a phenomenological model that constrains in a simple way the galaxy luminosity function evolution with the redshift, and fits all the existing source counts and redshift distributions, Cosmic Infrared Background intensity and fluctuations observations, from the mid-infrared to the submillimetre range. The model is based on template spectra of starburst and normal galaxies, and on the local infrared luminosity function. Although the Cosmic Infrared Background can be modeled with very different luminosity functions as long as the radiation production with redshift is the right one, the number counts, and the anisotropies of the unresolved background, imply that the luminosity function must change dramatically with redshift, with a rapid evolution of the high-luminosity sources (L>3 10^{11} L_{\odot}) from z=0 to z=1 which then stay rather constant up to redshift 5. The derived evolution of the IR luminosity function may be linked to a bimodal star formation process, one associated with the quiescent and passive phase of the galaxy evolution and one associated with the starburst phase, trigerred by merging and interactions. The latter dominates the infrared and submillimetre ouput energy of the Universe. The model is intended as a convenient tool to plan further observations, as illustrated through predictions for Herschel, Planck and ALMA observations. Our model predictions for given wavelengths, together with some usefull routines, are available for general use.