Simulation data

Introduction

The 2015 Planck data release is supported by a set of simulated maps of the sky, by astrophysical component, and of that sky as seen by Planck (fiducial mission realizations), together with separate sets of Monte Carlo realizations of the CMB and the instrument noise. They contain the dominant instrumental (detector beam, bandpass, and correlated noise properties), scanning (pointing and flags), and analysis (map-making algorithm and implementation) effects. These simulations have been described in Planck-2015-A12.

In addition to the baseline maps made from the data from all detectors at a given frequency for the entire mission, there are a number of data cuts that are mapped both for systematics tests and to support cross-spectral analyses. These include

1. detector subsets (“detsets”), comprising the individual unpolarized detectors and the polarized detector quadruplets corresponding to each leading trailing horn pair. Note that HFI sometimes refers to full channels as detset0; here detset only refers to subsets of detectors.
2. mission subsets, comprising the surveys, years, and half-missions, with exact boundary definitions given in Planck-2015-A02 and Planck-2015-A07 for LFI and HFI, respectively.
3. half-ring subsets, comprising the data from either the first or the second half of each pointing-period ring

The various combinations of these data cuts then define 1134 maps, as enumerated in the top section of Table 1 from Planck-2015-A12. The different types of map are then named according to their included detectors (channel or detset), interval (mission, half-mission, year or survey), and ring-content (full or half-ring); for example the baseline maps are described as channel/mission/full, etc.

The simulation process consists of

1. modeling each astrophysical component of the sky emission for each Planck detector, using Planck data and the relevant characteristics of the Planck instruments.
2. simulating each detector's observation of each sky component following the Planck scanning strategy and using the best estimates of the detector's beam and noise properties (obtained in flight), then combining these timelines into a single one per detector, and projecting these simulated timelines onto observed maps (the fiducial sky), as is done with the on-orbit data;
3. generating Monte Carlo realizations of the CMB and of the noise, again following the Planck scanning strategy and using our best estimates of the detector beams and noise properties respectively.

The first step is performed by the Planck Sky Model (PSM), and the last two by the Planck Simulation Tools (PST), both of which are described in the sections below.

The production of a full focal plane (FFP) simulation, and including the many MC realizations of the CMB and the noise, requires both HFI and LFI data and includes large, computationally challenging, MC realizations. They are too large to be generated on either of the DPC's own cluster. Instead the PST consists of three distinct tools, each designed to run on the largest available supercomputers, that are used to generate the fiducial sky realization, the CMB MC, and the noise MC respectively. The simulations delivered here are part of the 8th generation FFP simulations, known as FFP8. They were primarily generated on the National Energy Research Scientific Computing Center (NERSC) in the USA and at CSC–IT Center for Science (CSC) in Finland.

The fiducial realizations include instrument noise, astrophysical foregrounds, and the lensed scalar, tensor, and non-Gaussian CMB components, and are primarily designed to support the validation and verification of analysis codes. To test our ability to detect tensor modes and non-Gaussianity, we generate five CMB realizations with various cosmologically interesting -- but undeclared -- values of the tensor-to-scalar ratio r and non-Gaussianity parameter f_{NL}. To investigate the impact of differences in the bandpasses of the detectors at any given frequency, the foreground sky is simulated using both the individual detector bandpasses and a common average bandpass, to include and exclude the effects of bandpass mismatch. To check that the PR2-2015 results are not sensitive to the exact cosmological parameters used in FFP8 we subsequently generated FFP8.1, exactly matching the PR2-2015 cosmology.

Since mapmaking is a linear operation, the easiest way to generate all of these different realizations is to build the full set of maps of each of six components:

1. the lensed scalar CMB ({\tt cmb\_scl});
2. the tensor CMB ({\tt cmb\_ten});
3. the non-Gaussian complement CMB ({\tt cmb\_ngc});
4. the forgreounds including bandpass mismatch ({\tt fg\_bpm});
5. the foregrounds excluding bandpass mismatch ({\tt fg\_nobpm});
6. the noise (noise).

We then sum these, weighting the tensor and non-Gaussian complement maps with \sqrt{r} and f_{\rm NL}, respectively, and including one of the two foreground maps, to produce 10 total maps of each type. The complete fiducial data set then comprises 18,144 maps.

While the full set of maps can be generated for the fiducial cases, for the 10^4-realization MC sets this would result in some 10^7 maps and require about $6$\,PB of storage. Instead, therefore, the number of realizations generated for each type of map is chosen to balance the improved statistics it supports against the computational cost of its generation and storage. The remaining noise MCs sample broadly across all data cuts, while the additional CMB MCs are focused on the channel/half-mission/full maps and the subset of the detset/mission/full maps required by the "commander" component separation code Planck-2015-A10.

The goal of FFP8 is to simulate the Planck mission as accurately as possible; however, there are a number of known systematic effects that are not included, either because they are removed in the pre-processing of the time-ordered data (TOD), or because they are insufficiently well-characterized to simulate reliably, or because their inclusion (simulation and removal) would be too computationally expensive. These systematic effects are discussed in detail in Planck-2015-A02 and Planck-2015-A07 and include:

1. cosmic ray glitches (HFI);
2. spurious spectral lines from the 4-K cooler electronics (HFI);
3. nonlinearity in the analogue-to-digital converter (HFI);
4. imperfect reconstruction of the focal plane geometry.

Note that if the residuals from the treatment of any of these effects could be mapped in isolation, then maps of such systematics could simply be added to the existing FFP8 maps to improve their correspondence to the real data.

The Planck Sky Model

Overall description

The Planck Sky Model, PSM, consists of a set of data and of code used to simulate sky emission at millimeter-wave frequencies; it is described in detail in Delabrouille et al., (2013), henceforth the PSM paper..

The Planck Sky Model is available here: http://www.apc.univ-paris7.fr/~delabrou/PSM/psm.html

The main simulations used to test and validate the Planck data analysis pipelines (and, in particular, component separation) makes use of simulations generated with version 1.9 of the PSM software. The total sky emission is built from the CMB plus ten foreground components, namely thermal dust, spinning dust, synchrotron, CO lines, free-free, thermal Sunyaev-Zel'dovich (SZ) effect (with first order relativistic corrections), kinetic SZ effect, radio and infrared sources, Cosmic Infrared Background (CIB).

The CMB is modeled using CAMB. It is based on adiabatic initial perturbations, with the following cosmological parameters:

The foreground components generated by the PSM are diffuse Galactic emission from thermal dust, spinning dust, relativistic electrons (synchrotron emission), free-free radiation, and CO (the $J$=1$\rightarrow$0, $J$=2$\rightarrow$1, and $J$=3$\rightarrow$2 lines at $115.27$, $230.54$, and $345.80$\,GHz, respectively),

The Galactic ISM emission comprises five components: thermal dust, spinning dust, synchrotron, free-free, CO lines (the J=1->0, J=2->1, and J=3->2 lines at $115.27$, $230.54$, and $345.80$\,GHz, respectively), and plus the cosmic infrared background (CIB), emission from radio sources, and the thermal and kinetic Sunyaev-Zeldovich (SZ) effects.

The thermal dust emission is modelled using single-frequency template maps of the intensity and polarization, together with a pixel-dependent emission law. For FFP8 the thermal dust emission templates are derived from the Planck 353 GHz observations.

The spinning-dust map used for FFP8 simulations is a simple realization of the spinning dust model, post-processed to remove negative values occurring in a few pixels because of the generation of small-scale fluctuations on top of the spinning dust template extracted from WMAP data.

The FFP8 synchrotron emission is modelled on the basis of the template emission map observed at 408 MHz by Haslam et al. (1982). This template synchrotron map is extrapolated in frequency using a spectral index map corresponding to a simple power law.

The free-free spectral dependence is modelled in FFP8 by assuming a constant electron temperature $T_{e} = 7000$ K, as described in \cite{2003MNRAS.341..369D}.

The radio sources are modelled in FFP8 in a different way from the pre-launch versions of the PSM. For ``strong radio sources (S_{30}>0.5 Jy), we use radio sources at 0.84, 1.4, or 4.85\,GHz. For sources observed at two of these frequencies, we extrapolate or interpolate to the third frequency assuming the spectral index estimated from two observed. For sources observed at only one frequency, we use differential source counts to obtain the ratio of steep- to flat-spectrum sources in each interval of flux density considered. From this ratio, we assign spectral indices (randomly) to each source within each flux density interval. Fiducial Gaussian spectral index distributions as a function of spectral class are obtained from the literature. These are then adjusted slightly until there is reasonable agreement between the PSM differential counts and the predicted model counts predicted.

For ``faint radio sources (S_{30}<= 0.5 Jy), the pre-launch PSM showed a deficit of sources resulting from inhomogeneities in surveys at different depths. We address this issue by constructing a simulated catalogue of sources at 1.4\,GHz. We replace the simulated sources by the observed ones, wherever possible. If, however, in any particular pixel, we have a shortfall of observed sources, we make up the deficit with the simulated sources. Every source in this new catalogue is given a model-derived spectral class. We thus assign a spectral index to each source based on the spectral class, and model the spectrum of each source using four power laws. We also assume some steepening of the spectral index with frequency, with fiducial values of the steepening obtained from the literature.

We combine the faint and strong radio source catalogues we constructed and compute the differential source counts on these sources between 0.005 Jy and 1 Jy. Finally we also model the polarization of these radio sources using the measured polarization fractions from the literature; for each simulated source we draw a polarization fraction at random from the list of real sources of the same spectral type.

The SZ clusters are simulated following the model of Delabrouille, Melin, and Bartlett (DMB) as implemented in the PSM. A catalogue of halos is drawn from a Poisson distribution of the mass function with a limiting mass of M_{500,true}>2x10^{13} M_\odot. We use the pressure profile from the literature to model the thermal SZ emission of each halo given its redshift and mass. We determine the cluster temperature and assume that the profiles are isothermal. These steps allow us to compute the first-order thermal relativistic correction and the kinetic SZ effect for each cluster, both of which are included in the simulation. Finally, we inject catalogued clusters following the same model, and remove from the simulation corresponding clusters in each redshift and mass range. Hence the SZ simulation features the majority of known X-ray and optical clusters, and is fully consistent with X-ray scaling laws and observed \Planck\ SZ counts.

The CIB model used to simulate FFP8 relies on the distribution of individual galaxies in template maps based on the distribution of dark matter at a range of relevant redshifts. We assume the CIB galaxies can be grouped into three different populations (proto-spheroid, spiral, starburst). Within each population, galaxies have the same SED, while the flux density is randomly distributed according to redshift-dependent number counts obtained from JCMT/SCUBA-2 observations and the Planck ERCSC, as well as observations from Herschel-SPIRE and AzTEC/ASTE. We use the Class software to generate dark matter maps at 17 different redshifts between 1 and 5.5. Since the galaxy distribution does not exactly follow the dark matter distribution, we modify the a_{lm} coefficients of dark matter anisotropies given by Class. Template maps generated from the a_{lm} coefficients are then exponentiated to avoid negative pixels. Galaxies are randomly distributed with a probability of presence proportional to the pixel values of the template maps. One map is generated for each population, at each redshift, and associated with a redshifted SED depending on the population. The emission of these maps (initially at a reference frequency) can be extrapolated to any frequency using the associated redshifted SED. By summing the emission of all maps, we can generate CIB maps at any frequency in the range of validity of our model.

See Planck-2015-A12 and references therein for a very detailed explanation of the procedures to simulate each of the components.

PSM Products

To build maps corresponding to the Planck channels, the models described above are convolved with the spectral response of the channel in question. The products given here are for the full frequency channels, and as such they are not used in the Planck specific simulations, which use only individual detector channels. The frequency channel spectral responses used (given in the RIMO), are averages of the responses of the detectors of each frequency channel weighted as they are in the mapmaking step. They are provided for the purpose of testing user's own software of simulations and component separation.

PSM maps of the CMB and of the ten foregrounds are given in the following map products:

HFI

• HFI_SimMap_cmb_2048_R1.10.fits
• HFI_SimMap_co_2048_R1.10.fits
• HFI_SimMap_firb_2048_R1.10.fits
• HFI_SimMap_strongps_2048_R1.10.fits
• HFI_SimMap_faintps_2048_R1.10.fits
• HFI_SimMap_freefree_2048_R1.10.fits
• HFI_SimMap_synchrotron_2048_R1.10.fits
• HFI_SimMap_thermaldust_2048_R1.10.fits
• HFI_SimMap_spindust_2048_R1.10.fits
• HFI_SimMap_kineticsz_2048_R1.10.fits
• HFI_SimMap_thermalsz_2048_R1.10.fits

LFI

• LFI_SimMap_cmb_1024_R1.10.fits
• LFI_SimMap_firb_1024_R1.10.fits
• LFI_SimMap_strongps_1024_R1.10.fits
• LFI_SimMap_faintps_1024_R1.10.fits
• LFI_SimMap_freefree_1024_R1.10.fits
• LFI_SimMap_synchrotron_1024_R1.10.fits
• LFI_SimMap_thermaldust_1024_R1.10.fits
• LFI_SimMap_spindust_1024_R1.10.fits
• LFI_SimMap_kineticsz_1024_R1.10.fits
• LFI_SimMap_thermalsz_1024_R1.10.fits

Each file contains a single BINTABLE extension with either a single map (for the CMB file) or one map for each HFI/LFI frequency (for the foreground components). In the latter case the columns are named F030, F044 ,F070,F100, F143, … , F857. Units are microK_CMB for the CMB, K_CMB at 30, 44 and 70 GHz and MJy/sr for the others. The structure is given below for multi-column files.

Note: Original PSM foreground components has been generated at NSIDE 2048 and using a gaussian beam of 4 arcmin, LFI maps where then smoothed to LFI resolution (32.0, 27.0 and 13.0 arcmin for the 30, 44 and 70 GHz) and donwgraded at NSIDE 1024. LFI CMB maps has been smoothed at 13.0 arcmin (70 GHz resolution) and downgraded at NSIDE 1024.

HFI FITS file structure

1. EXTNAME = 'SIM-MAP' : Data columns
Column Name	Data Type	Units	Description
F100	Real*4	MJy/sr	100GHz signal map
F143	Real*4	MJy/sr	143GHz signal map
F217	Real*4	MJy/sr	217GHz signal map
F353	Real*4	MJy/sr	353GHz signal map
F545	Real*4	MJy/sr	545GHz signal map
F857	Real*4	MJy/sr	857GHz signal map
Keyword	Data Type	Value	Description
PIXTYPE	string	HEALPIX
COMP	string	component	Astrophysical omponent
COORDSYS	string	GALACTIC	Coordinate system
ORDERING	string	NESTED	Healpix ordering
NSIDE	Int	2048	Healpix Nside for LFI and HFI, respectively
FIRSTPIX	Int*4	0	First pixel number
LASTPIX	Int*4	50331647	Last pixel number, for LFI and HFI, respectively
BAD_DATA	Real*4	-1.63750E+30	Healpix bad pixel value
BEAMTYPE	string	GAUSSIAN	Type of beam
BEAMSIZE	Real*4	size	Beam size in arcmin
PSM-VERS	string		PSM Versions used

LFI FITS file structure

1. EXTNAME = 'SIM-MAP' : Data columns
Column Name	Data Type	Units	Description
F030	Real*4	KCMB	30GHz signal map
F044	Real*4	KCMB	44GHz signal map
F070	Real*4	KCMB	70GHz signal map
Keyword	Data Type	Value	Description
PIXTYPE	string	HEALPIX
COMP	string	component	Astrophysical omponent
COORDSYS	string	GALACTIC	Coordinate system
ORDERING	string	NESTED	Healpix ordering
NSIDE	Int	1024	Healpix Nside for LFI and HFI, respectively
FIRSTPIX	Int*4	0	First pixel number
LASTPIX	Int*4	12582911	Last pixel number, for LFI and HFI, respectively
BAD_DATA	Real*4	-1.63750E+30	Healpix bad pixel value
BEAMTYPE	string	GAUSSIAN	Type of beam
BEAMS_30	Real*4	32.0	Beam size at 30 GHz in arcmin
BEAMS_44	Real*4	27.0	Beam size at 44 GHz in arcmin
BEAMS_70	Real*4	13.0	Beam size at 70 GHz in arcmin
PSM-VERS	string		PSM Versions used

The Fiducial Sky Simulations

For each detector, fiducial time-ordered data are generated separately for each of the ten PSM components using the LevelS software as follows:

• the detector's beam and PSM map are converted to spherical harmonics using beam2alm and anafast respectively;
• the beam-convolved map value is calculated over a 3-dimensional grid of sky locations and beam orientations using conviqt;
• the map-based timelines are calculated sample-by-sample by interpolating over this grid using multimod;
• the catalogue-based timelines are produced sample-by-sample by beam-convolving any point source laying within a given angular distance of the pointing at each sample time using multimod.

For each frequency, fiducial sky maps are generated for

• the total signal (i.e. sky + instrument noise), for both the nominal mission and the halfrings thereof (see details)
• the foreground sky alone (excluding CMB but including noise),
• the point source sky, and
• the noise alone

All maps are built using the MADAM destriping map-maker interfaced with the TOAST data abstraction layer . In order to construct the total timelines required by each map, for each detector TOAST reads the various component timelines separately and sums then, and, where necessary, simulates and adds a noise realization time-stream on the fly. HFI frequencies are mapped at HEALPix resolution Nside=2048 using ring-length destriping baselines, while LFI frequencies are mapped at Nside=1024 using 1s baselines.

Products delivered

A single simulation is delivered, which is divided into two types of products:

1. six files of the full sky signal at each HFI and LFI frequency, and their corresponding halfring maps:

HFI_SimMap_100_2048_R1.10_nominal.fits	HFI_SimMap_100_2048_R1.10_nominal_ringhalf_1.fits	HFI_SimMap_100_2048_R1.10_nominal_ringhalf_2.fits
HFI_SimMap_143_2048_R1.10_nominal.fits	HFI_SimMap_143_2048_R1.10_nominal_ringhalf_1.fits	HFI_SimMap_143_2048_R1.10_nominal_ringhalf_2.fits
HFI_SimMap_217_2048_R1.10_nominal.fits	HFI_SimMap_217_2048_R1.10_nominal_ringhalf_1.fits	HFI_SimMap_217_2048_R1.10_nominal_ringhalf_2.fits
HFI_SimMap_353_2048_R1.10_nominal.fits	HFI_SimMap_353_2048_R1.10_nominal_ringhalf_1.fits	HFI_SimMap_353_2048_R1.10_nominal_ringhalf_2.fits
HFI_SimMap_545_2048_R1.10_nominal.fits	HFI_SimMap_545_2048_R1.10_nominal_ringhalf_1.fits	HFI_SimMap_545_2048_R1.10_nominal_ringhalf_2.fits
HFI_SimMap_857_2048_R1.10_nominal.fits	HFI_SimMap_857_2048_R1.10_nominal_ringhalf_1.fits	HFI_SimMap_857_2048_R1.10_nominal_ringhalf_2.fits

LFI_SimMap_030_1024_R1.10_nominal.fits	LFI_SimMap_030_1024_R1.10_nominal_ringhalf_1.fits	LFI_SimMap_030_1024_R1.10_nominal_ringhalf_2.fits
LFI_SimMap_044_1024_R1.10_nominal.fits	LFI_SimMap_044_1024_R1.10_nominal_ringhalf_1.fits	LFI_SimMap_044_1024_R1.10_nominal_ringhalf_2.fits
LFI_SimMap_070_1024_R1.10_nominal.fits	LFI_SimMap_070_1024_R1.10_nominal_ringhalf_1.fits	LFI_SimMap_070_1024_R1.10_nominal_ringhalf_2.fits

These files have the same structure as the equivalent SkyMap products described in the Frequency Maps chapter, namely one BINTABLE extension with three columns containing 1) Signal, 2) hit-count, and 3) variance. Units are K_CMB for all channels.

2. Three files containing 1) the sum of all astrophysical foregrounds, 2) the point sources alone, and 3) the noise alone: which are subproducts of the above, and are in the form of the PSM maps described in the previous section.

• HFI_SimMap_foreground_2048_R1.10_nominal.fits
• HFI_SimMap_noise_2048_R1.10_nominal.fits
• HFI_SimMap_ps_2048_R1.10_nominal.fits

• LFI_SimMap_foreground_1024_R1.10_nominal.fits
• LFI_SimMap_noise_1024_R1.10_nominal.fits
• LFI_SimMap_ps_1024_R1.10_nominal.fits

These files have the same structure as the PSM output maps described above, namely a single BINTABLE extension with 6 columns named F100 -- F857 each containing the given map for that HFI band and with 3 columns named F030, F044, F070 each containing the given map for that LFI band. Units are alway K_CMB.

Note that the CMB alone is not delivered as a separate product, but it can be recovered by simple subtraction of the component maps for the total signal map.

Monte Carlo realizations of CMB and of noise

The CMB MC set is generated using FEBeCoP, which generates an effective beam for each pixel in a map at each frequency by accumulating the weights of all pixels within a fixed distance of that pixel, summed over all observations by all detectors at that frequency. It then applies this effective beam pixel-by-pixel to each of 1000 input CMB sky realizations.

The noise MC set is generated just as the fiducial noise maps, using MADAM/TOAST. In order to avoid spurious correlations within and between the 1000 realizations, each stationary interval for each detector for each realization is generated from a distinct sub-sequence of a single statistically robust, extremely long period, pseudo-random number sequence.

Products delivered

100 realizations of the CMB (lensed) and of the noise are made available. They are named

• HFI_SimMap_cmb-{nnnn}_2048_R1.nn_nominal.fits
• HFI_SimMap_noise-{nnnn}_2048_R1.nn_nominal.fits
• LFI_SimMap_cmb-{nnnn}_2048_R1.nn_nominal.fits
• LFI_SimMap_noise-{nnnn}_2048_R1.nn_nominal.fits

where nnnn ranges from 0000 to 0099.

The FITS file structure is the same as for the other similar products above, with a single BINTABLE extension with six columns, one for each HFI frequency, named F100, F143, … , F857 and with three columns, one for each LFI frequency, named F030, F044, F070. Units are always microK_CMB(NB: due to an error in the HFI file construction, the unit keywords in the headers indicate K_CMB, the "micro" is missing there).

Products delivered in PLA

In the Planck Legacy Archive, the 2015 simulations are delivered as follows :

• 1000 Monte-Carlo realizations from the ffp8 series (noise and CMB scalar) - full maps only

• the 1001 Monte-Carlo realizations from the ffp9 series (CMB scalar only) - full maps only

• the ffp8 fiducial sky for band-pass mismatch and no band-pass mismatch

The files can be found in the PLA in Maps - Advanced search - Simulated maps search.

All the rest of the ffp8 and ffp9 simulations can be found on NERSC at /project/projectdirs/planck/data/

Lensing Simulations

The lensing simulations package contains 100 realisations of the Planck 2014 "MV" lensing potential estimate, as well as the input lensing realizations. They can be used to determine error bars as well eas effective normalizations for cross-correlation with other tracers of lensing. These simulations are of the lensing convergence map contained in the Lensing map release file. The production and characterisation of this lensing potential map are described in detail in , which also describes the procedure used to generate the realizations given here.

The simulations are delivered as a gzipped tarball of approximately 8 GB in size. For delivery purposes, the package has been split into 2GB chunks using the unix command

split -d -b 2048m

After downloading the individual chunks, the full tarball can be reconstructed with the command

cat COM_SimMap_Lensing_R0.00.tar.* | tar xvf -

The contents of the tarball are described below:

Contents of COM_SimMap_Lensing_2048_R0.00.tgz

Filename	Format	Description
obs_klms/sim_????_klm.fits	HEALPIX FITS format alm, with [math] L_{\rm max} = 2048 [/math]	Contains the simulated convergence estimate [math] \hat{\kappa}_{LM} = \frac{1}{2} L(L+1)\hat{\phi}_{LM} [/math] for each simulation.
sky_klms/sim_????_klm.fits	HEALPIX FITS format alm, with [math] L_{\rm max} = 2048 [/math]	Contains the input lensing convergence for each simulation.
inputs/mask.fits.gz	HEALPIX FITS format map, with [math] N_{\rm side} = 2048 [/math]	Contains the lens reconstruction analysis mask.
inputs/cls/cl??.dat	ASCII text file, with columns = ([math]L[/math], [math]C_L [/math])	Contains the fiducial theory CMB power spectra for TT, EE, BB, [math] \kappa \kappa [/math] and [math] T \kappa [/math], with temperature and polarization in units of [math] \mu K [/math].