The pre-launch Planck Sky Model: a model of sky emission at submillimetre to centimetre wavelengths
University of Paris1, DSM2, University of Toronto3, University of Paris-Sud4, International School for Advanced Studies5, INAF6, University of Cambridge7, University of Toulouse8, Institut d'Astrophysique de Paris9, Imperial College London10, University of Porto11, University of Manchester12, Max Planck Society13, European Space Agency14, École Centrale Paris15, École Polytechnique Fédérale de Lausanne16, CERN17, University of Savoy18, University of Padua19, Joseph Fourier University20, University of Rome Tor Vergata21, California Institute of Technology22, Haverford College23, Heidelberg University24, University of Oviedo25
TL;DR: The Planck Sky Model (PSM) as mentioned in this paper is a parametric model for the generation of all-sky, few arcminute resolution maps of sky emission at submillimetre to centimetre wavelengths, in both intensity and polarisation.
Abstract: We present the Planck Sky Model (PSM), a parametric model for the generation of all-sky, few arcminute resolution maps of sky emission at submillimetre to centimetre wavelengths, in both intensity and polarisation. Several options are implemented to model the cosmic microwave background, Galactic diffuse emission (synchrotron, free-free, thermal and spinning dust, CO lines), Galactic H-II regions, extragalactic radio sources, dusty galaxies, and thermal and kinetic Sunyaev-Zeldovich signals from clusters of galaxies. Each component is simulated by means of educated interpolations/extrapolations of data sets available at the time of the launch of the Planck mission, complemented by state-of-the-art models of the emission. Distinctive features of the simulations are: spatially varying spectral properties of synchrotron and dust; different spectral parameters for each point source; modeling of the clustering properties of extragalactic sources and of the power spectrum of fluctuations in the cosmic infrared background. The PSM enables the production of random realizations of the sky emission, constrained to match observational data within their uncertainties, and is implemented in a software package that is regularly updated with incoming information from observations. The model is expected to serve as a useful tool for optimizing planned microwave and sub-millimetre surveys and to test data processing and analysis pipelines. It is, in particular, used for the development and validation of data analysis pipelines within the planck collaboration. A version of the software that can be used for simulating the observations for a variety of experiments is made available on a dedicated website.
Citations
More filters
Cardiff University1, California Institute of Technology2, Atacama Large Millimeter Submillimeter Array3, University of Toronto4, Harvard University5, University of Minnesota6, University of British Columbia7, National Institute of Standards and Technology8, Stanford University9, University of California, San Diego10, University of Chicago11, Canadian Institute for Advanced Research12
TL;DR: An excess of B-mode power over the base lensed-ΛCDM expectation is found in the range 30 < ℓ < 150, inconsistent with the null hypothesis at a significance of >5σ, and it is shown that systematic contamination is much smaller than the observed excess.
Abstract: We report results from the BICEP2 experiment, a cosmic microwave background (CMB) polarimeter specifically designed to search for the signal of inflationary gravitational waves in the B -mode power spectrum around l∼80 . The telescope comprised a 26 cm aperture all-cold refracting optical system equipped with a focal plane of 512 antenna coupled transition edge sensor 150 GHz bolometers each with temperature sensitivity of ≈300 μK CMB s √ . BICEP2 observed from the South Pole for three seasons from 2010 to 2012. A low-foreground region of sky with an effective area of 380 square deg was observed to a depth of 87 nK deg in Stokes Q and U . In this paper we describe the observations, data reduction, maps, simulations, and results. We find an excess of B -mode power over the base lensed-ΛCDM expectation in the range 30 5σ . Through jackknife tests and simulations based on detailed calibration measurements we show that systematic contamination is much smaller than the observed excess. Cross correlating against WMAP 23 GHz maps we find that Galactic synchrotron makes a negligible contribution to the observed signal. We also examine a number of available models of polarized dust emission and find that at their default parameter values they predict power ∼(5–10)× smaller than the observed excess signal (with no significant cross-correlation with our maps). However, these models are not sufficiently constrained by external public data to exclude the possibility of dust emission bright enough to explain the entire excess signal. Cross correlating BICEP2 against 100 GHz maps from the BICEP1 experiment, the excess signal is confirmed with 3σ significance and its spectral index is found to be consistent with that of the CMB, disfavoring dust at 1.7σ . The observed B -mode power spectrum is well fit by a lensed-ΛCDM+tensor theoretical model with tensor-to-scalar ratio r=0.20 +0.07 −0.05 , with r=0 disfavored at 7.0σ . Accounting for the contribution of foreground, dust will shift this value downward by an amount which will be better constrained with upcoming data sets
1,954 citations
Cites methods from "The pre-launch Planck Sky Model: a ..."
...The polarization angles are taken from the PSM. DDM2: “Data Driven Model 2” (DDM2) constructed using all publicly available information from Planck....
[...]
...The Planck dust model map at 353 GHz is scaled to 150 GHz assuming a con- 30http://www3.imperial.ac.uk/people/c.contaldi/ fgpol 31 http://www.apc.univ-paris7.fr/~delabrou/PSM/ psm.html stant emissivity value of 1.6 and a constant temperature of 19.6 K....
[...]
...PSM: Planck Sky Model (PSM) (Delabrouille et al. 2013) version 1.7.8, run as a “Prediction” with 15% dust intrinsic polarization fraction and Galactic magnetic field pitch angle of −30◦31....
[...]
TL;DR: The European Space Agency'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 as discussed by the authors.
Abstract: The European Space Agency’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. In March 2013, ESA and the Planck Collaboration released the initial cosmology products based on the first 15.5 months of Planck data, along with a set of scientific and technical papers and a web-based explanatory supplement. 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 cosmic microwave background (CMB) and diffuse extragalactic foregrounds, a catalogue of compact Galactic and extragalactic sources, and a list of sources detected through the Sunyaev-Zeldovich 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 ΛCDM 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σ. 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 (σ8) derived from CMB data and that derived from Sunyaev-Zeldovich data. The Planck and WMAP power spectra are offset from each other by an average level of about 2% around the first acoustic peak. Analysis of Planck polarization data is not yet mature, therefore polarization results are not released, although the robust detection of E-mode polarization around CMB hot and cold spots is shown graphically.
1,719 citations
TL;DR: 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 as discussed by the authors.
Abstract: 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.
1,302 citations
Cites methods from "The pre-launch Planck Sky Model: a ..."
...For the fiducial realization, maps were made of the total observation (CMB, foregrounds, and noise) at each frequency for the nominal mission period, using the Planck Sky Model (Delabrouille et al. 2013)....
[...]
...For the fiducial realization, maps were made of the total observation (CMB, foregrounds, and noise) at each frequency for the nominal mission period, using the Planck Sky Model (Delabrouille et al. 2012)....
[...]
TL;DR: Strong evidence for dust and no statistically significant evidence for tensor modes is found and various model variations and extensions are probe, including adding a synchrotron component in combination with lower frequency data, and find that these make little difference to the r constraint.
Abstract: 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 deg2 patch of sky centered on RA 0h, Dec. −57.5deg. 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×353 cross-correlation in B-modes at high significance. We fit the single- and cross-frequency power spectra at frequencies above 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). 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 r0.05<0.12 at 95% confidence. Marginalizing over dust and r, lensing B-modes are detected at 7.0σ significance.
1,255 citations
TL;DR: In this article, the authors characterized the effective beams, the effective beam window functions and the associated errors for the Planck High Frequency Instrument (HFI) detectors, including the effect of the optics, detectors, data processing and the scan strategy.
Abstract: This paper characterizes the effective beams, the effective beam window functions and the associated errors for the Planck High Frequency Instrument (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 cosmic microwave background 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 transfer function; and c) the merging of several surveys to produce maps. The time response transfer 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 5 degrees) 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 solid angle. 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 errors 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 100 GHz up to multiple l similar to 1500, below 0.5% at 143 and 217 GHz up to l similar to 2000.
1,124 citations
Cites methods from "The pre-launch Planck Sky Model: a ..."
...The pre-launch version of the PSM is described by Delabrouille et al. (2013), and has been modified for the present work as described in Planck Collaboration (2013)....
[...]
References
More filters
TL;DR: In this article, a reprocessed composite of the COBE/DIRBE and IRAS/ISSA maps, with the zodiacal foreground and confirmed point sources removed, is presented.
Abstract: We present a full-sky 100 μm map that is a reprocessed composite of the COBE/DIRBE and IRAS/ISSA maps, with the zodiacal foreground and confirmed point sources removed. Before using the ISSA maps, we remove the remaining artifacts from the IRAS scan pattern. Using the DIRBE 100 and 240 μm data, we have constructed a map of the dust temperature so that the 100 μm map may be converted to a map proportional to dust column density. The dust temperature varies from 17 to 21 K, which is modest but does modify the estimate of the dust column by a factor of 5. The result of these manipulations is a map with DIRBE quality calibration and IRAS resolution. A wealth of filamentary detail is apparent on many different scales at all Galactic latitudes. In high-latitude regions, the dust map correlates well with maps of H I emission, but deviations are coherent in the sky and are especially conspicuous in regions of saturation of H I emission toward denser clouds and of formation of H2 in molecular clouds. In contrast, high-velocity H I clouds are deficient in dust emission, as expected. To generate the full-sky dust maps, we must first remove zodiacal light contamination, as well as a possible cosmic infrared background (CIB). This is done via a regression analysis of the 100 μm DIRBE map against the Leiden-Dwingeloo map of H I emission, with corrections for the zodiacal light via a suitable expansion of the DIRBE 25 μm flux. This procedure removes virtually all traces of the zodiacal foreground. For the 100 μm map no significant CIB is detected. At longer wavelengths, where the zodiacal contamination is weaker, we detect the CIB at surprisingly high flux levels of 32 ± 13 nW m-2 sr-1 at 140 μm and of 17 ± 4 nW m-2 sr-1 at 240 μm (95% confidence). This integrated flux ~2 times that extrapolated from optical galaxies in the Hubble Deep Field. The primary use of these maps is likely to be as a new estimator of Galactic extinction. To calibrate our maps, we assume a standard reddening law and use the colors of elliptical galaxies to measure the reddening per unit flux density of 100 μm emission. We find consistent calibration using the B-R color distribution of a sample of the 106 brightest cluster ellipticals, as well as a sample of 384 ellipticals with B-V and Mg line strength measurements. For the latter sample, we use the correlation of intrinsic B-V versus Mg2 index to tighten the power of the test greatly. We demonstrate that the new maps are twice as accurate as the older Burstein-Heiles reddening estimates in regions of low and moderate reddening. The maps are expected to be significantly more accurate in regions of high reddening. These dust maps will also be useful for estimating millimeter emission that contaminates cosmic microwave background radiation experiments and for estimating soft X-ray absorption. We describe how to access our maps readily for general use.
15,988 citations
"The pre-launch Planck Sky Model: a ..." refers methods in this paper
...The E(B − V) map of Schlegel et al. (1998)....
[...]
...The map of Dickinson et al. (2003) that uses a Hα emission template (a combination of WHAM (Reynolds et al. 1998) and SHASSA (Gaustad et al. 2001) data) corrected for extinction using the E(B − V) map of Schlegel et al. (1998); 2....
[...]
...The 100 μm map used is the one of Schlegel et al. (1998) or a modified version of this map, displayed in Fig....
[...]
TL;DR: In this paper, the authors presented a reprocessed composite of the COBE/DIRBE and IRAS/ISSA maps, with the zodiacal foreground and confirmed point sources removed.
Abstract: We present a full sky 100 micron map that is a reprocessed composite of the COBE/DIRBE and IRAS/ISSA maps, with the zodiacal foreground and confirmed point sources removed. Before using the ISSA maps, we remove the remaining artifacts from the IRAS scan pattern. Using the DIRBE 100 micron and 240 micron data, we have constructed a map of the dust temperature, so that the 100 micron map can be converted to a map proportional to dust column density. The result of these manipulations is a map with DIRBE-quality calibration and IRAS resolution.
To generate the full sky dust maps, we must first remove zodiacal light contamination as well as a possible cosmic infrared background (CIB). This is done via a regression analysis of the 100 micron DIRBE map against the Leiden- Dwingeloo map of H_I emission, with corrections for the zodiacal light via a suitable expansion of the DIRBE 25 micron flux. For the 100 micron map, no significant CIB is detected. In the 140 micron and 240 micron maps, where the zodiacal contamination is weaker, we detect the CIB at surprisingly high flux levels of 32 \pm 13 nW/m^2/sr at 140 micron, and 17 \pm 4 nW/m^2/sr at 240 micron (95% confidence). This integrated flux is ~2 times that extrapolated from optical galaxies in the Hubble Deep Field.
The primary use of these maps is likely to be as a new estimator of Galactic extinction. We demonstrate that the new maps are twice as accurate as the older Burstein-Heiles estimates in regions of low and moderate reddening. These dust maps will also be useful for estimating millimeter emission that contaminates CMBR experiments and for estimating soft X-ray absorption.
14,295 citations
TL;DR: In this article, a combination of seven-year data from WMAP and improved astrophysical data rigorously tests the standard cosmological model and places new constraints on its basic parameters and extensions.
Abstract: The combination of seven-year data from WMAP and improved astrophysical data rigorously tests the standard cosmological model and places new constraints on its basic parameters and extensions. By combining the WMAP data with the latest distance measurements from the baryon acoustic oscillations (BAO) in the distribution of galaxies and the Hubble constant (H0) measurement, we determine the parameters of the simplest six-parameter ΛCDM model. The power-law index of the primordial power spectrum is ns = 0.968 ± 0.012 (68% CL) for this data combination, a measurement that excludes the Harrison–Zel’dovich–Peebles spectrum by 99.5% CL. The other parameters, including those beyond the minimal set, are also consistent with, and improved from, the five-year results. We find no convincing deviations from the minimal model. The seven-year temperature power spectrum gives a better determination of the third acoustic peak, which results in a better determination of the redshift of the matter-radiation equality epoch. Notable examples of improved parameters are the total mass of neutrinos, � mν < 0.58 eV (95% CL), and the effective number of neutrino species, Neff = 4.34 +0.86 −0.88 (68% CL), which benefit from better determinations of the third peak and H0. The limit on a constant dark energy equation of state parameter from WMAP+BAO+H0, without high-redshift Type Ia supernovae, is w =− 1.10 ± 0.14 (68% CL). We detect the effect of primordial helium on the temperature power spectrum and provide a new test of big bang nucleosynthesis by measuring Yp = 0.326 ± 0.075 (68% CL). We detect, and show on the map for the first time, the tangential and radial polarization patterns around hot and cold spots of temperature fluctuations, an important test of physical processes at z = 1090 and the dominance of adiabatic scalar fluctuations. The seven-year polarization data have significantly improved: we now detect the temperature–E-mode polarization cross power spectrum at 21σ , compared with 13σ from the five-year data. With the seven-year temperature–B-mode cross power spectrum, the limit on a rotation of the polarization plane due to potential parity-violating effects has improved by 38% to Δα =− 1. 1 ± 1. 4(statistical) ± 1. 5(systematic) (68% CL). We report significant detections of the Sunyaev–Zel’dovich (SZ) effect at the locations of known clusters of galaxies. The measured SZ signal agrees well with the expected signal from the X-ray data on a cluster-by-cluster basis. However, it is a factor of 0.5–0.7 times the predictions from “universal profile” of Arnaud et al., analytical models, and hydrodynamical simulations. We find, for the first time in the SZ effect, a significant difference between the cooling-flow and non-cooling-flow clusters (or relaxed and non-relaxed clusters), which can explain some of the discrepancy. This lower amplitude is consistent with the lower-than-theoretically expected SZ power spectrum recently measured by the South Pole Telescope Collaboration.
11,309 citations
"The pre-launch Planck Sky Model: a ..." refers background in this paper
...The dependence of counts on σ8 is important and its precise value still uncertain, with a typical allowed range extending from about 0.81 as obtained from primary CMB constraints (Komatsu et al. 2011) to about 0.87 as obtained from weak lensing (Massey et al. 2007)....
[...]
TL;DR: In this article, high-resolution N-body simulations show that the density profiles of dark matter halos formed in the standard CDM cosmogony can be fit accurately by scaling a simple universal profile.
Abstract: High resolution N-body simulations show that the density profiles of dark matter halos formed in the standard CDM cosmogony can be fit accurately by scaling a simple “universal” profile. Regardless of their mass, halos are nearly isothermal over a large range in radius, but significantly shallower than r -2 near the center and steeper than r -2 in the outer regions. The characteristic overdensity of a halo correlates strongly with halo mass in a manner consistent with the mass dependence of the epoch of halo formation. Matching the shape of the rotation curves of disk galaxies with this halo structure requires (i) disk mass-to-light ratios to increase systematically with luminosity, (ii) halo circular velocities to be systematically lower than the disk rotation speed, and (iii) that the masses of halos surrounding bright galaxies depend only weakly on galaxy luminosity. This offers an attractive explanation for the puzzling lack of correlation between luminosity and dynamics in observed samples of binary galaxies and of satellite companions of bright spiral galaxies, suggesting that the structure of dark matter halos surrounding bright spirals is similar to that of cold dark matter halos.
7,622 citations
"The pre-launch Planck Sky Model: a ..." refers methods in this paper
...Baryons are pasted on clusters assuming an electron density profile (either a Navarro et al. 1996 profile, modified according to Nagai et al. 2007 and Arnaud et al. 2010 or a β-model) and using a mass-temperature relation (Pierpaoli et al. 2003), the normalisation of which is a free parameter, or a…...
[...]
TL;DR: This paper considers the requirements and implementation constraints on a framework that simultaneously enables an efficient discretization with associated hierarchical indexation and fast analysis/synthesis of functions defined on the sphere and demonstrates how these are explicitly satisfied by HEALPix.
Abstract: HEALPix the Hierarchical Equal Area isoLatitude Pixelization is a versatile structure for the pixelization of data on the sphere. An associated library of computational algorithms and visualization software supports fast scientific applications executable directly on discretized spherical maps generated from very large volumes of astronomical data. Originally developed to address the data processing and analysis needs of the present generation of cosmic microwave background experiments (e.g., BOOMERANG, WMAP), HEALPix can be expanded to meet many of the profound challenges that will arise in confrontation with the observational output of future missions and experiments, including, e.g., Planck, Herschel, SAFIR, and the Beyond Einstein inflation probe. In this paper we consider the requirements and implementation constraints on a framework that simultaneously enables an efficient discretization with associated hierarchical indexation and fast analysis/synthesis of functions defined on the sphere. We demonstrate how these are explicitly satisfied by HEALPix.
5,518 citations
"The pre-launch Planck Sky Model: a ..." refers methods in this paper
...Most of these data sets, in HEALPix format (Górski et al. 2005), have been downloaded from the LAMBDA web site....
[...]
...The HEALPix sky pixelisation scheme (Górski et al. 2005), in Galactic coordinates, is used in the current implementation....
[...]