arxiv: v1 [astro-ph.co] 20 Mar 2013

Transcription

1 Astronomy & Astrophysics manuscript no. component separation c ESO 2013 March 22, 2013 arxiv: v1 [astro-ph.co] 20 Mar 2013 Panck 2013 resuts. XII. Component separation Panck Coaboration: P. A. R. Ade 90, N. Aghanim 63, C. Armitage-Capan 95, M. Arnaud 77, M. Ashdown 74,6, F. Atrio-Barandea 20, J. Aumont 63, C. Baccigaupi 89, A. J. Banday 98,11, R. B. Barreiro 70, J. G. Bartett 1,72, E. Battaner 100, K. Benabed 64,97, A. Benoît 61, A. Benoit-Lévy 28,64,97, J.-P. Bernard 11, M. Bersanei 37,52, P. Bieewicz 98,11,89, J. Bobin 77, J. J. Bock 72,12, A. Bonadi 73, L. Bonavera 70, J. R. Bond 9, J. Borri 15,93, F. R. Bouchet 64,97, F. Bouanger 63, M. Bridges 74,6,67, M. Bucher 1, C. Burigana 51,35, R. C. Buter 51, J.-F. Cardoso 78,1,64, A. Cataano 79,76, A. Chainor 67,74,13, A. Chambau 77,17,63, R.-R. Chary 60, X. Chen 60, L.-Y Chiang 66, H. C. Chiang 30,7, P. R. Christensen 85,40, S. Church 94, D. L. Cements 59, S. Coombi 64,97, L. P. L. Coombo 27,72, F. Couchot 75, A. Couais 76, B. P. Cri 72,86, M. Cruz 22, A. Curto 6,70, F. Cuttaia 51, L. Danese 89, R. D. Davies 73, R. J. Davis 73, P. de Bernardis 36, A. de Rosa 51, G. de Zotti 48,89, J. Deabrouie 1, J.-M. Deouis 64,97, F.-X. Désert 56, C. Dickinson 73, J. M. Diego 70, H. Doe 63,62, S. Donzei 52, O. Doré 72,12, M. Douspis 63, J. Dunkey 95, X. Dupac 43, G. Efstathiou 67, T. A. Enßin 82, H. K. Eriksen 68, E. Fagarone 76, F. Finei 51,53, O. Forni 98,11, M. Fraiis 50, A. A. Fraisse 30, E. Franceschi 51, S. Gaeotta 50, K. Ganga 1, M. Giard 98,11, G. Giardino 44, Y. Giraud-Héraud 1, J. Gonzáez-Nuevo 70,89, K. M. Górski 72,102, S. Gratton 74,67, A. Gregorio 38,50, A. Gruppuso 51, F. K. Hansen 68, D. Hanson 83,72,9, D. Harrison 67,74, G. Heou 12, S. Henrot-Versié 75, C. Hernández-Monteagudo 14,82, D. Herranz 70, S. R. Hidebrandt 12, E. Hivon 64,97, M. Hobson 6, W. A. Homes 72, A. Hornstrup 18, W. Hovest 82, G. Huey 33, K. M. Huffenberger 101, T. R. Jaffe 98,11, A. H. Jaffe 59, J. Jewe 72, W. C. Jones 30, M. Juvea 29, E. Keihänen 29, R. Keskitao 25,15, T. S. Kisner 81, R. Kneiss 42,8, J. Knoche 82, L. Knox 31, M. Kunz 19,63,3, H. Kurki-Suonio 29,46, G. Lagache 63, A. Lähteenmäki 2,46, J.-M. Lamarre 76, A. Lasenby 6,74, R. J. Laureijs 44, C. R. Lawrence 72, M. Le Jeune 1, S. Leach 89, J. P. Leahy 73, R. Leonardi 43, J. Lesgourgues 96,88, M. Liguori 34, P. B. Lije 68, M. Linden-Vørne 18, M. López-Caniego 70, P. M. Lubin 32, J. F. Macías-Pérez 79, B. Maffei 73, D. Maino 37,52, N. Mandoesi 51,5,35, A. Marcos-Cabaero 70, M. Maris 50, D. J. Marsha 77, P. G. Martin 9, E. Martínez-Gonzáez 70, S. Masi 36, S. Matarrese 34, F. Matthai 82, P. Mazzotta 39, P. R. Meinhod 32, A. Mechiorri 36,54, L. Mendes 43, A. Mennea 37,52, M. Migiaccio 67,74, K. Mikkesen 68, S. Mitra 58,72, M.-A. Mivie-Deschênes 63,9, A. Moneti 64, L. Montier 98,11, G. Morgante 51, D. Mortock 59, A. Moss 91, D. Munshi 90, P. Nasesky 85,40, F. Nati 36, P. Natoi 35,4,51, C. B. Netterfied 23, H. U. Nørgaard-Niesen 18, F. Novieo 73, D. Novikov 59, I. Novikov 85, I. J. O Dwyer 72, S. Osborne 94, C. A. Oxborrow 18, F. Paci 89, L. Pagano 36,54, F. Pajot 63, R. Paadini 60, D. Paoetti 51,53, B. Partridge 45, F. Pasian 50, G. Patanchon 1, T. J. Pearson 12,60, O. Perdereau 75, L. Perotto 79, F. Perrotta 89, V. Pettorino 19, F. Piacentini 36, M. Piat 1, E. Pierpaoi 27, D. Pietrobon 72, S. Paszczynski 75, P. Patania 71, E. Pointecouteau 98,11, G. Poenta 4,49, N. Ponthieu 63,56, L. Popa 65, T. Poutanen 46,29,2, G. W. Pratt 77, G. Prézeau 12,72, S. Prunet 64,97, J.-L. Puget 63, J. P. Rachen 24,82, W. T. Reach 99, R. Reboo 69,16,41, M. Reinecke 82, M. Remazeies 63,1, C. Renaut 79, A. Renzi 89, S. Ricciardi 51, T. Rier 82, I. Ristorcei 98,11, G. Rocha 72,12, C. Rosset 1, G. Roudier 1,76,72, M. Rowan-Robinson 59, J. A. Rubiño-Martín 69,41, B. Rushome 60, E. Saerno 10, M. Sandri 51, D. Santos 79, G. Savini 87, F. Schiavon 51, D. Scott 26, M. D. Seiffert 72,12, E. P. S. Sheard 13, L. D. Spencer 90, J.-L. Starck 77, R. Stompor 1, R. Sudiwaa 90, R. Sunyaev 82,92, F. Sureau 77, D. Sutton 67,74, A.-S. Suur-Uski 29,46, J.-F. Sygnet 64, J. A. Tauber 44, D. Tavagnacco 50,38, L. Terenzi 51, L. Toffoatti 21,70, M. Tomasi 52, M. Tristram 75, M. Tucci 19,75, J. Tuovinen 84, M. Türer 57, G. Umana 47, L. Vaenziano 51, J. Vaiviita 46,29,68, B. Van Tent 80, J. Varis 84, M. Vie 50,55, P. Vieva 70, F. Via 51, N. Vittorio 39, L. A. Wade 72, B. D. Wandet 64,97,33, I. K. Wehus 72, A. Wikinson 73, J.-Q. Xia 89, D. Yvon 17, A. Zacchei 50, and A. Zonca 32 (Affiiations can be found after the references) Preprint onine version: March 22, 2013 ABSTRACT Panck has produced detaied a-sky observations over nine frequency bands between 30 and 857 GHz. These observations aow robust reconstruction of the primordia cosmic microwave background (CMB) temperature fuctuations over neary the fu sky, as we as new constraints on Gaactic foregrounds, incuding therma dust and ine emission from moecuar carbon monoxide (CO). This paper describes the component separation framework adopted by Panck for many cosmoogica anayses, incuding CMB power spectrum determination and ikeihood construction on arge anguar scaes, studies of primordia non-gaussianity and statistica isotropy, the integrated Sachs-Wofe effect (ISW), and gravitationa ensing, and searches for topoogica defects. We test four foreground-ceaned CMB maps derived using quaitativey different component separation agorithms. The quaity of our reconstructions is evauated through detaied simuations and interna comparisons, and shown through various tests to be internay consistent and robust for CMB power spectrum and cosmoogica parameter estimation up to = The parameter constraints on ΛCDM cosmoogies derived from these maps are consistent with those presented in the cross-spectrum based Panck ikeihood anaysis. We choose two of the CMB maps for specific scientific goas. We aso present maps and frequency spectra of the Gaactic ow-frequency, CO, and therma dust emission. The component maps are found to provide a faithfu representation of the sky, as evauated by simuations, with the argest bias seen in the CO component at 3 %. For the ow-frequency component, the spectra index varies widey over the sky, ranging from about β = 4 to 2. Considering both morphoogy and prior knowedge of the ow frequency components, the index map aows us to associate a steep spectra index (β < 3.2) with strong anomaous microwave emission, corresponding to a spinning dust spectrum peaking beow 20 GHz, a fat index of β > 2.3 with strong free-free emission, and intermediate vaues with synchrotron emission. Key words. cosmoogy: observations Corresponding author: M. Ashdown, maja1@mrao.cam.ac.uk. 1

2 1. Introduction This paper, one of a set associated with the 2013 reease of data from the Panck 1 mission (Panck Coaboration I 2013), describes the component separation techniques appied to the Panck data to produce maps of the cosmic microwave background (CMB) temperature anisotropies (see Fig. 1) and of diffuse foregrounds. The sky at miimetre and sub-miimetre waveengths contains a weath of cosmoogica and astrophysica information. Accessing it is an inversion process, known as component separation, to extract the sources of emission contributing to a set of maps observed at different frequencies. Panck gives us a powerfu data set to unock new information in this manner by observing the entire sky from 30 to 857 GHz in nine frequency bands at higher anguar resoution and sensitivity than its predecessors. Accurate and detaied component separation is a centra objective of the mission. We divide the foregrounds into two distinct categories: diffuse emission from the Gaaxy and compact sources. The Gaactic foregrounds are the principa source of contamination of the CMB on arge anguar scaes, with fuctuation power decreasing roughy as a power aw towards higher mutipoes (Bennett et a. 2003). They are dominated by synchrotron, free-free and Anomaous Microwave Emission (AME, ascribed to spinning dust grains, at frequencies beow 70 GHz) and by rotationa ine emission from carbon monoxide (CO) moecues and therma dust emission at frequencies above 100 GHz. Extragaactic foregrounds, on the other hand, dominate the sma-scae contamination of the CMB. They arise from discrete, individuay detectabe compact sources and the coective emission from unresoved radio and infrared (IR) sources, and aso from the Sunyaev-Zedovich (SZ) effect in gaaxy custers (Panck Coaboration XVIII 2011; Panck Coaboration XXVIII 2013; Panck Coaboration XXIX 2013). In the Panck anayses, these foregrounds are deat with in a variety of ways. At the power spectrum and ikeihood eve, the extragaactic foregrounds are modeed with parameterized power spectra, appropriate to their statistica isotropy, over regions restricted to ow Gaactic emission (Panck Coaboration XV 2013). Component separation as described in the present paper aims at removing Gaactic emission to produce CMB maps covering the argest possibe sky area for studies of the argescae properties and higher-order statistics of the CMB. In addition, this component separation provides a reconstruction of the diffuse emission from our Gaaxy. Detaied studies of specific extragaactic foregrounds, such as the cosmic infrared background (CIB) (Panck Coaboration XVIII 2013) and the diffuse Sunyaev-Zedovich (SZ) signa (Panck Coaboration XXI 2013), empoy methods taiored to their particuar needs. Buiding on previous work (Leach et a. 2008), we approach CMB extraction with a phiosophy designed to ensure robustness by appying four distinct agorithms based on two different methodoogies. The first avoids any assumptions concerning the foregrounds and reies soey on a minimum variance criterion for the data component possessing a backbody spectrum (i.e., the CMB), whie the second methodoogy reies on parametric modeing of the foregrounds in either rea or harmonic space. 1 Panck ( is a project of the European Space Agency (ESA) with instruments provided by two scientific consortia funded by ESA member states (in particuar the ead countries France and Itay), with contributions from NASA (USA) and teescope refectors provided by a coaboration between ESA and a scientific consortium ed and funded by Denmark. We evauate the performance of these component separation agorithms through detaied simuations, and we examine the robustness of the recovered CMB maps by comparing them, their power spectra, and their resuting cosmoogica constraints. As a diagnostic, we aso briefy examine their higher-order statistics. The CMB resuts presented in this work serve a number appications. We use the rea-space modeing to produce a cean CMB map and power spectra on arge anguar scaes, where diffuse Gaactic emission is the main contaminant, to construct the ikeihood function at ow mutipoes; this is then combined with the high mutipoe ikeihood function that modes extragaactic foregrounds with power spectra (Panck Coaboration XV 2013). The high resoution CMB maps are used as a check on primary cosmoogica constraints (see beow), for ensing studies (Panck Coaboration XVII 2013), studies of the integrated Sachs-Wofe effect (Panck Coaboration XIX 2013), of the isotropy of the CMB (Panck Coaboration XXIII 2013), of non-gaussian statistics (Panck Coaboration XXIV 2013), in searches for topoogica defects (Panck Coaboration XXV 2013), and for examination of the geometry and topoogy of the Universe (Panck Coaboration XXVI 2013). In addition, we present maps of diffuse Gaactic emission divided into ow- and high-frequency components, as we as a moecuar CO component. We judge the adequacy of this reconstruction through simuations and by comparison with known properties of the diffuse Gaactic foregrounds. The paper is organized as foows. In Sect. 2 we discuss the expected sources of sky emission over the Panck frequency interva and how they are modeed. Then in Sect. 3 we detai the overa approach and introduce the four component separation methods. In Sect. 4 we present the Panck data set and preprocessing procedure, and we describe our simuations. This is foowed by a presentation of the derived CMB maps and their characterization in Sect. 5. Section 6 is dedicated to power spectra and cosmoogica parameter constraints obtained from these maps, and Sect. 7 to studies of higher-order statistics. Section 8 presents a reconstruction of the diffuse Gaactic foregrounds, and Sect. 9 concudes. We reegate detais of the agorithms to appendices. 2. The sky at Panck frequencies The properties of Gaactic emission vary significanty across the Panck frequency range from 30 to 857 GHz. At frequencies beow 70 GHz, the dominant radiation processes are: synchrotron emission from cosmic ray eectrons interacting with the Gaactic magnetic fied (e.g., Hasam et a. 1982; Reich & Reich 1988; Broadbent et a. 1989; Davies et a. 1996; Patania et a. 2003; Bennett et a. 2003; God et a. 2011); therma Bremsstrahung (or free-free emission) from eectron-eectron and eectron-ion scattering (e.g., Banday et a. 2003; Dickinson et a. 2003; Davies et a. 2006; Ghosh et a. 2012; Panck Coaboration Int. XII 2013; Panck Coaboration XX 2011); and AME from dust grains (Kogut 1996; Leitch et a. 1997; Banday et a. 2003; Lagache 2003; de Oiveira-Costa et a. 2004; Finkbeiner et a. 2004; Davies et a. 2006; Bonadi et a. 2007; Dober & Finkbeiner 2008; Mivie-Deschênes et a. 2008; Ysard et a. 2010; God et a. 2011; Panck Coaboration XX 2011), possiby due to their rotationa ine emission (Draine & Lazarian 1998; Ai-Haïmoud et a. 2009; Ysard & Verstraete 2010; Hoang & Lazarian 2012). Over the frequency range covered by Panck, both synchrotron and free-free spectra are we approximated by power aws in brightness temperature, T B ν β, with the synchrotron index, β synch, ranging from 3.2 to 2.8 (Davies et a. 2

3 C-R 300 µk 300 Fig. 1: Foreground-ceaned CMB maps derived by Commander-Ruer,, and. Note that the map has been fied in smoothy inside a 3 % Gaactic mask. 1996) and the free-free index, β ff, ying between 2.2 and 2.1. Less is known about the AME spectrum, but spinning dust modes with a spectrum peaking at frequencies beow 20 GHz (in brightness temperature units) adequatey describe current observations 2. Above the peak, the spectrum appears consistent with a power-aw (e.g., Banday et a. 2003; Davies et a. 2006; Dober & Finkbeiner 2008; Ghosh et a. 2012). In addition to these three, the existence of a fourth ow-frequency foreground component, known as the Gaactic haze, has been caimed, possiby due to a hard-spectrum synchrotron popuation near the Gaactic centre (e.g., Finkbeiner 2004; Dober & Finkbeiner 2008; Pietrobon et a. 2012; Panck Coaboration Int. IX 2013). At frequencies higher than 100 GHz, therma dust emission dominates over most of the sky and is commony described by a modified backbody spectrum with power-aw emissivity, ɛ ν ν β d, and temperature, T d. Both the temperature and spectra index, β d, vary spatiay. Prior to Panck, the best-fitting singe component dust mode had a temperature T d 18 K and spectra index β d 1.7 (Finkbeiner et a. 1999; Bennett et a. 2003; God et a. 2011), athough there is evidence of fattening of the spectra index from around 1.8 in the far-infrared to 1.55 in the 2 Note that we adopt brightness temperature for AME in this paper, whie many other pubications adopt fux density. When comparing peak frequencies, it is usefu to note that that a spectrum that has a maximum at 30 GHz in fux density peaks at 17 GHz in brightness temperature. microwave region (Panck Coaboration 2012), the interpretation of which is sti under study. In addition to these diffuse Gaactic components, extragaactic emission contributes at Panck frequencies. In particuar, a arge number of radio and far-infrared (FIR; Panck Coaboration XIII 2011) gaaxies, custers of gaaxies and the Cosmic Infrared Background (CIB; Panck Coaboration XVIII 2011) produce a statisticay isotropic foreground, with frequency spectra we approximated by modes simiar to those appicabe to the Gaactic foregrounds (modified backbody spectra, power aws, etc.). Except for a frequency-dependent absoute offset, which may be removed as part of the overa offset remova procedure, these extragaactic components are therefore typicay absorbed by either the ow-frequency or therma dust components during component separation. No specia treatment is given here to extragaactic foregrounds, beyond the masking of bright objects. Dedicated scientific anayses of these sources are described in detai in Panck Coaboration XVIII (2011), Panck Coaboration XXVIII (2013), and Panck Coaboration XXIX (2013). In the Panck ikeihood, extragaactic sources are modeed in terms of power spectrum tempates at high (Panck Coaboration XV 2013). Other reevant sources incude emission from moecuar couds, supernova remnants, and compact H ii regions inside our own Gaaxy, as we as the therma and kinetic SZ effects, due to inverse Compton scattering of CMB photons off free eectrons 3

4 Tabe 1: Overview and comparison of component separation agorithms. Characteristic Commander-Ruer Method Bayesian parameter Interna inear Interna tempate Spectra parameter estimation combination fitting estimation Domain Pixe Needet Pixe Spherica harmonic Channes [GHz] Effective beam FWHM [arcmin] max none Fig. 2: Combined Gaactic (CG) emission masks for the Panck data, corresponding to sky fractions of 20, 40, 60, 70, 75, 80, 90, 97, and 99 %. The masks are named CG20, etc. in ionized media. Panck provides new and important information on a these processes, as described both in the foowing and in the companion papers Panck Coaboration XIII (2013), Panck Coaboration XXI (2013), and Panck Coaboration (2013b). In particuar, Panck s frequency range, anguar resoution and sensitivity make it a powerfu probe of therma dust, resuting in new and tight constraints on dust temperature and emissivity. The same frequencies aso aow extraction of the first ever fu-sky maps of the emission resuting from the CO J=1 0, J=2 1 and J=3 2 rotationa transitions at 115, 230 and 345 GHz, respectivey (Panck Coaboration XIII 2013). The focus of this paper is to reconstruct the CMB anisotropies over a arge sky fraction, expoiting ony the Panck frequency bands. We aso present a detaied reconstruction of the therma dust emission at high frequencies, as we as CO emission ines. At ow frequencies and over the region used for CMB anaysis, the tota foreground contribution is we approximated by a singe power aw (see Sect. 8). We therefore mode the sum of a ow-frequency foregrounds by a power aw with spatiay varying spectra index whose numerica vaue in any pixe resuts from the infuence of the dominant foreground component at that ocation. The fu anaysis of diffuse foregrounds, using anciary data to resove the individua components at ow frequencies, wi be presented in a forthcoming pubication. 3. Approach to component separation The rich content of the Panck data encourages appication of severa component separation techniques. We consider four, as summarized in Tabe 1, which we cassify according to one of two different genera methodoogies. The first makes minima assumptions concerning the foregrounds and seeks ony to minimize the variance of the CMB, i.e., the sky component possessing a backbody spectrum. We impement this approach with a needet (waveet on the sphere) version of the interna inear combination (ILC) agorithm (; Deabrouie et a. 2009), and aso with a tempate-based method to remove foreground contamination from the CMB-dominant bands. These foreground tempates are constructed from the owest and highest frequency channes (Spectra Estimation Via Expectation Maximization, ; Fernández-Cobos et a. 2012). The second methodoogy uses parametric modeing of the foregrounds. In our rea space impementation, we expore mode parameters through Bayesian parameter estimation techniques, fitting a parametric signa mode per pixe (Commander; Eriksen et a. 2006, 2008); a simiar impementation is presented by Stompor et a. (2009). To estimate spectra indices robusty in pixe space, this procedure requires identica anguar resoution across a frequencies incuded in the anaysis, and is therefore imited in resoution by the 30 GHz LFI channe. However, this is sufficient to generate the ow-resoution CMB map and power spectrum sampes required for the ow mutipoe part of the Panck ikeihood function for cosmoogica parameters (Panck Coaboration XV 2013). To produce fu resoution maps, we use the resuting ow-resoution spectra parameter sampes to sove for the component ampitudes, in an extension to the method known as Ruer (we refer to the combined method as Commander-Ruer, or C-R). In our fourth technique, we impement a CMB-oriented parametric approach that fits the ampitude and spectra parameters of CMB and foregrounds in the harmonic domain (Spectra Matching Independent Component Anaysis, ; Cardoso et a. 2008). Detais of each agorithm are given in the appendices. We now turn to their appication to the data and evauate their performance using simuations. 4. Data, simuations and masks We use the data set from the first 15.5 months of Panck observations, corresponding to 2.6 sky surveys, from both the Low Frequency Instrument (LFI) and High Frequency Instrument (HFI). The primary inputs for component separation are the frequency channe maps, incuding haf-ring maps, bandpasses, and beam characteristics; a fu description of these products is given in Panck Coaboration II (2013) and Panck Coaboration VI (2013). No specia corrections are made for zodiaca ight emission (ZLE; Panck Coaboration VI 2013) in the anayses presented here. The ZLE is not stationary on the sky, since it depends on Panck s position and scanning strategy. Therefore the frequency maps contain a projected version of the emission averaged over the nomina mission. Despite this, a series of exporatory anayses showed that our agorithms naturay correct for this component within their existing mode space. It was aso found that arger CMB residuas were induced when appying a correction based on a ZLE mode than when appying no correction, most ikey due to uncertainties in the mode itsef. 4

5 Fig. 3: Summary of Component Separation (CS) confidence masks. Each pixe is encoded in terms of a sum in which Commander-Ruer equas 1 (ight bue), equas 2 (dark red), equas 4 (yeow), and equas 8 (ight red). The masks are named CS-CR75, CS-93, CS-76, and CS-89, respectivey, refecting their accepted sky fraction. The union mask (U73), used for evauation purposes in this paper, removes a cooured pixes. B C-R Fig. 4: Beam transfer functions of the four foreground-ceaned CMB maps. To evauate and vaidate our agorithms, we anayse a arge suite of reaistic simuations, the so-caed Fu Foca Pane (FFP) simuations, based on detaied modes of the instrument and sky. The version used for this data reease is denoted FFP6, and is described in Panck Coaboration ES (2013). The simuation procedure generates time streams for each detector, incorporating the sateite pointing, the individua detector beams, bandpasses, noise properties, and data fags, and then produces simuated frequency channe maps through the mapmaking process. For the input sky, we use the Panck Sky Mode (PSM), which incudes the CMB, diffuse Gaactic emission (synchrotron, freefree, therma dust, AME, and moecuar CO ines), and compact sources (therma and kinetic SZ effects, radio sources, infrared sources, the CIB, and utra-compact H ii regions). The preaunch version of the PSM is described by Deabrouie et a. (2012), and has been modified for the present work as described in Panck Coaboration ES (2013). Each FFP data set consists of three parts: the simuated observations, Monte Caro reaizations of the CMB, and Monte Caro reaizations of the instrumenta noise. For both the data and the simuations, we reconstruct the CMB and foregrounds from the fu frequency channe maps and the corresponding haf-ring maps, which are made from the data in the first haf or second haf of each stabe pointing period. The haf-ring maps can be used to obtain an estimate of the noise in each channe by taking haf of the difference between the two maps, thereby normaizing the noise eve to that of the fu map. This is referred to as the haf-ring haf-difference (HRHD) map. The signas fixed to the sky wi be canceed eaving ony the noise contribution. The HRHD map can be treated as a reaization of the same underying noise processes and it can be used to estimate the power spectrum, and other properties, of the noise. If there are noise correations between the haf-ring maps, then the estimates of the noise properties thus obtained can be biased. This is the case for HFI channes; the cosmic ray gitch remova (Panck Coaboration VI 2013; Panck Coaboration X 2013) induces correations that ead to the noise power spectrum being underestimated by a few percent at high when using the HRHD maps. 0 µk 10 Fig. 5: Standard deviation between the four foreground-ceaned CMB maps. A maps have been downgraded to a HEALPix resoution of N side = 128. The differences are typicay ess than 5 µk at high Gaactic atitudes, demonstrating that the maps are consistent over a arge part of the sky. Prior to processing the data through each component separation pipeine, we define masks for the point sources and bright Gaactic regions. Point source masking is based on the source cataogues obtained by fitering the input sky maps with the Mexican Hat Waveet 2 (MHW2) fiter and appying a 4 σ threshod for the LFI bands and a 5 σ threshod for the HFI bands (Panck Coaboration XVIII 2011; Panck Coaboration XXVIII 2013). The mask radius of each source is different for the LFI and HFI. Due to the arge beam size of LFI channes, we define a variabe masking radius for each source according to its signa-to-noise ratio (S/N) as r = (2 og(a/m)) 1/2 /(2 2 og 2) FWHM, where r is the radius, A is the S/N, and m is the maximum ampitude (given in units of the background noise eve) aowed for the tai of unmasked point sources; we set m = 0.1, which is a compromise between masking the source tais and minimizing the number of masked pixes. 5

6 C-R - C-R - C-R µk 30 Fig. 6: Pairwise differences between foreground-ceaned CMB maps. A maps have been downgraded to a HEALPix resoution of N side = 128 to show the arge-scae differences. The ine-ike discontinuities in the differences invoving is due to the two different regions used in this agorithm to cean the sky (see Appendix C for detais). For HFI, the mask radius around each source is 1.27 FWHM, using the average FWHM obtained from the effective beams. A basic set of Gaactic masks is defined as foows. We subtract a CMB estimate from the 30 and 353 GHz maps, mask point sources, and smooth the resuting maps by a Gaussian with FWHM of 5. We then threshod and combine them, generating a series of masks with different amounts of avaiabe sky. The resuting combined Gaactic (CG) masks, shown in Fig. 2, correspond to sky fractions of 20, 40, 60, 70, 75, 80, 90, 97, and 99 %, and are named CG20, etc. 5. CMB Maps We begin the discussion of our resuts by presenting the foreground-ceaned CMB maps. These maps are shown in Fig. 1 for each of the four component separation agorithms. Aready from this figure it is cear that the wide frequency coverage and high anguar resoution of Panck aow a faithfu reconstruction 6

7 C-R 30 µk 30 Fig. 7: CMB residua maps from the FFP6 simuation. A monopoe determined at high Gaactic atitude has been subtracted from the maps, and they have been downgraded to a HEALPix resoution of N side = 128 to show the arge-scae features. The residuas presented here provide a conservative estimate of those expected in the data (see text for detais). of the CMB fied over most of the sky. The fuctuations appear visuay consistent with the theoretica expectation of a Gaussian and isotropic signa everywhere except inside a sma band very cose to the Gaactic pane 3 Each CMB map is accompanied by its own confidence mask outside which the corresponding soution is considered statisticay robust, shown in Fig. 3; for a definition of each mask, see Appendices A D. Accepted sky fractions are 75, 93, 76, and 89 %, respectivey, for Commander-Ruer,,, and. These masks are denoted CS-CR75, CS-93, CS- 76, CS-89, respectivey. The union of the confidence masks accepts 73 % of the sky and is denoted U73. It is adopted as the defaut mask for evauation purposes in this paper. In addition to the CMB maps from the fu data set, the hafring frequency maps have been processed by each agorithm to provide haf-ring CMB maps. They are used to provide estimates of the instrumenta noise contribution to the foreground-ceaned maps in the power spectrum anaysis (see Sect. 6). The agorithms were aso used to process Monte Caro simuations: 1000 reaizations of the CMB and 1000 reaizations of noise. They 3 Note that, being defined in harmonic space, empoys a smooth fiing process inside a sma Gaactic mask to prevent foreground residuas from eaking from ow to high Gaactic atitudes, and therefore appears visuay different from the other three soutions in this respect; see Appendix D. are not used in the anayses presented in this paper, but are used by Panck Coaboration XXIII (2013) and Panck Coaboration XXIV (2013). The beam transfer functions of the foreground-ceaned CMB maps have been estimated for each agorithm, as shown in Fig. 4. The anguar resoution of the,, and maps corresponds to a Gaussian beam with FWHM of 5. The difference between and / is due to their different treatment of the HEALPix 4 pixe window function (Górski et a. 2005). The deviation of beam from a Gaussian shape at > 2800 is caused by the ast needet window (see Appendix B). Commander-Ruer has a arger beam, because it is defined expicity as a weighted average of frequency maps in pixe space. Its resoution is equivaent to a Gaussian beam with FWHM of approximatey The beam transfer functions have been computed assuming the best-fit beam transfer function for each frequency channe, and the uncertainties in the atter have not been propagated to these estimates. In Fig. 5 we show the standard deviation per pixe among the four foreground-ceaned CMB maps downgraded to N side = 128, and in Fig. 6 we show a pairwise difference maps. Typica differences at high Gaactic atitudes are smaer than 5 µk. Considering the difference maps in more detai, it is cear that 4 7

8 the Commander-Ruer map is the most different from the other three, whereas and are the most simiar. This is not competey unexpected, because whie Commander-Ruer uses ony frequencies between 30 and 353 GHz in its soution, the other three codes additionay incude the dust-dominated 545 and 857 GHz maps. This difference in data seection may expain some of the coherent structures seen in Fig. 6. In particuar, the most striking arge-scae feature in the difference maps invoving Commander-Ruer is a arge negative band roughy foowing the eciptic pane. This is where the ZLE (Panck Coaboration VI 2013) is brightest. Since the ZLE is aso stronger at high frequencies, having a spectrum cose to that of therma dust, it is possibe that this pattern may be an imprint of residua ZLE either in the Commander-Ruer map, or in a of the other three maps. Both cases are pausibe. The Commander-Ruer soution may not have enough high-frequency information to distinguish between ZLE and norma therma dust emission, and, by assuming a therma dust spectrum for the entire high-frequency signa at 353 GHz, over-subtracts the ZLE at ower frequencies. It is aso possibe that the other three CMB soutions have positive ZLE residuas from extrapoating the high-frequency signa mode from 857 GHz to the CMB frequencies. Without an accurate and detaied ZLE mode, it is difficut to distinguish between these two possibiities. It is of course aso possibe that the true expanation is in fact unreated to ZLE, and the correation with the eciptic pane is accidenta. In either case, it is cear that the residuas are sma in ampitude, with peak-to-peak vaues typicay smaer than 10 µk, of which by far the most is contained in a quadrupoe aigned with the eciptic. This provides additiona evidence that residua ZLE is not important for the CMB power spectrum and cosmoogica parameter estimation, athough some care is warranted when using these maps to study the statistics of the very argest anguar scaes (e.g., Panck Coaboration XXIII 2013); checking consistency among a four maps for a given appication aeviates much of this concern. We end this section by showing in Fig. 7 a set of residua maps derived by anaysing the FFP6 simuation with exacty the same anaysis approaches as appied to the data. It is evident that produces the map with owest eve of residuas. Considering the morphoogy in each case, we see that the main contaminant for Commander-Ruer is under-subtracted free-free emission, whie for both and it is oversubtracted therma dust emission, and for it is undersubtracted therma dust emission. However, at high atitudes and outside the confidence masks, the residuas are generay beow a few µk in ampitude. It is aso worth noting that each agorithm has been optimized (in terms of mode definition, ocaization parameters, etc.) for the data, and the same configuration was subsequenty used for the FFP6 simuations without further tuning. The simuations presented here therefore provide a conservative estimate of the residuas in the data. This is aso refected in the fact that the differences between CMB reconstructions for the FFP6 simuations are arger than those found in the data. See Appendix E for further detais. 6. Power spectrum and cosmoogica parameters In this section we evauate the foreground-ceaned maps in terms of CMB power spectra and cosmoogica parameters. Our purpose in doing this is to show that the maps are consistent with the high- ikeihood obtained from the cross-spectrum anaysis of detector set and frequency maps in Panck Coaboration XV (2013), and with the cosmoogica parameters derived from ( +1)C /2π [ µk 2 ] C-R Fig. 8: Anguar power spectra of the foreground-ceaned CMB maps and haf-ring haf-difference (HRHD) maps. The spectra have been evauated using the U73 mask apodized with a 30 cosine function. them in Panck Coaboration XVI (2013). This aso estabishes the consistency between Panck s cosmoogica constraints and studies of the arge-scae structure and higher order statistics of the CMB Power spectra Figure 8 shows the power spectra of the foreground-ceaned CMB maps and the corresponding HRHD maps, evauated using the U73 mask with a 30 cosine apodization. The spectra have been corrected for the effect of the mask and the beam transfer function of each agorithm has been deconvoved. The spectra of the HRHD maps give an estimate of the instrumenta noise contribution to the power spectrum of the ceaned map. The correations between the HFI haf-ring frequency maps are inherited by the haf-ring CMB maps that use them as input. At sma anguar scaes, the CMB soution comes amost entirey from data in the HFI channes, and therefore the spectrum of the CMB HRHD maps is aso biased ow. At sma anguar scaes, the effective noise eves of,, and are very simiar, and ower than that of Commander-Ruer. The ast has arger noise because it operates entirey in pixe space and therefore appies the same weights to a mutipoes. It cannot take advantage of the changing signato-noise ratio of the frequency channes with anguar scae. We can estimate the contribution of residua foregrounds to the foreground-ceaned CMB maps by making use of the FFP6 simuations. In addition to processing the simuated frequency maps, the maps of the individua input sky components were processed by the agorithms after fixing their parameters or weights to the vaues obtained from the observed maps. Figure 9 shows the power spectra of the simuated FFP6 components, in this case CMB, noise and the sum of the foreground components. The top pane shows the spectra computed using the union mask derived from the simuation with a 30 cosine apodization. The tota foreground contribution becomes comparabe to the CMB signa at The bottom pane shows the same computed with an apodized point source mask appied to the maps (i.e., 8

9 ( +1)C /2π [ µk 2 ] ( +1)C /2π [ µk 2 ] C-R C-R ( + 1)C /2π [µ K 2 ] C R Fig. 10: Estimates of the CMB power spectra from the foreground-ceaned maps, computed by XFaster. The soid ines show the spectra after subtracting the best-fit mode of residua foregrounds. The vertica dotted ine shows the maximum mutipoe ( = 2000) used in the ikeihood for fitting the foreground mode and cosmoogica parameters (see Sect for further detais). The dashed ines show the spectra before residua foreground subtraction. Chain Monte Caro samper. To speed up this process, we additionay use PICO (Parameters for the Impatient COsmoogist, Fendt & Wandet 2008), a too which interpoates the CMB power spectra and matter power spectra as a function of cosmoogica parameters. Fig. 9: Anguar power spectra of FFP6 simuated components evauated over the common mask (top) and the common point source mask (bottom), both apodized with a 30 cosine function. Three components are shown: the CMB (dashed ine); noise (dot-dashed ine); and the sum of a foregrounds (soid ine). A noninear scae is used on the horizonta axis to show a the features of the spectra. no diffuse masking, athough this mask does removes a arge part of the Gaactic pane). The residua foreground contribution is arger at a anguar scaes, but sti it ony becomes comparabe to the CMB signa at 1800 in the worst case. For both masks, has the smaest residua foreground contamination at arge anguar scaes, which is aso demonstrated in Fig. 7. A more detaied examination of the contribution of the individua foreground components to the power spectrum is in Appendix E Likeihood and cosmoogica parameters We estimate the binned power spectra with XFaster (Rocha et a. 2009, 2010, 2011) and determine cosmoogica parameter constraints using a correated Gaussian ikeihood. Parameter constraints are derived using a Metropois-Hastings Markov Mode and methods We compute the power spectrum for each foreground-ceaned map over the mutipoe range , whie parameter constraints are derived using ony ; as shown in Appendix E through simuations, modeing errors become non-negigibe between = 2000 and For parameter estimation, we adopt a standard six-parameter ΛCDM mode, and impose an informative Gaussian prior of τ = ± 0.014, since poarization data are not incuded in this anaysis. Whie the foreground-ceaned maps shoud have minima contamination from diffuse Gaactic emission, they do contain significant contamination from unresoved extragaactic sources. These contributions are most easiy modeed in terms of residua power spectra, therefore we marginaize over the corresponding parameters at the power spectrum eve. To the six ΛCDM parameters, describing the standard cosmoogy, we add two foreground parameters, A ps, the ampitude of a Poisson component (and hence constant, C = A ps ), and A c, the ampitude of a custered component with shape D = ( +1)C /2π 0.8. Both are expressed in terms of D at = 3000 in units of µk 2. The power spectrum cacuation is based on the haf-ring haf-sum (HRHS) and HRHD CMB maps (see Sect. 5); the atter is used to estimate the noise bias in the power spectra extracted from the HRHS maps. From these, we cacuate the pseudospectra, C and Ñ (Hivon et a. 2002), respectivey, after appying the U73 mask. These are used as inputs to XFaster together 9

10 Ω h 2 b τ H [km s 1 Mpc 1 ] Ω c h 2 n s A ps [µk 2 ] θ og(10 10 A s ) A [µk 2 ] c C R CamSpec Pik C R CamSpec Pik 20 0 C R CamSpec Pik Fig. 11: Comparison of cosmoogica and foreground parameter vaues estimated from the foreground-ceaned CMB maps for max = 2000 (in red) and those obtained with CamSpec and Pik ikeihoods (in bue). The vaues of the foreground parameters are not shown for CamSpec and Pik, since they use a different foreground mode. ( + 1)(C C C a m S p e c )/2π [µ K 2 ] C R To avoid aiasing of power from arge to sma scaes, which woud add an offset between the signa-pus-noise and noise pseudo-spectra at high, we use the apodized version of the U73 mask. The known mismatch in the noise eve between the spectra due to the correation between the haf-ring maps is not expicity corrected. It is eft to be absorbed into the two foreground parameters. Using the pseudo-spectra and XFaster, we then reconstruct an estimate of the power spectrum of each foreground-ceaned HRHS map, removing the noise bias as estimated from the corresponding HRHD map. To this end we appy an iterative scheme starting from a fat spectrum mode. The resut is a binned power spectrum and the associated Fisher matrix, which are then used to construct the ikeihood, approximated here by a correated Gaussian distribution. Fig. 12: Residuas of a map-based best-fit modes reative to CamSpec best-fit mode (assuming a prior on τ) for max = with the beam transfer functions provided by each method (see Fig. 4). To study consistency in the ow- range, we fit a twoparameter q n (ampitude-tit) mode reative to the Panck bestfit ΛCDM mode on the form, C = q(/ pivot ) n C bf, using a pixespace ikeihood for maps smoothed to 6 FWHM; see Panck Coaboration XV (2013) for further agorithmic detais. 10

11 Power spectrum tit, n Power spectrum ampitude, q C-R Fig. 13: Low- power spectrum ampitude and tit constraints measured reative to the best-fit Panck ΛCDM mode derived from foreground-ceaned CMB maps smoothed to 6 FWHM. The cross shows the best-fit mode (q, n) = (1, 0). between the frequency channe maps used to construct them, or their beam uncertainties. These effects are accounted for in the ikeihoods in Panck Coaboration XV (2013). We have aso adopted a very simpe two-parameter mode for the residua extragaactic foregrounds. Despite these imitations, the four CMB maps yied cosmoogica parameters in agreement with the cross-spectrum based ikeihoods for a basic six-parameter ΛCDM mode. Thus we can be confident that the CMB maps are consistent with the power spectrum anaysis. Before concuding this section, we show in Fig. 13 the resuts from a two-parameter fit of an ampitude-tit mode to each of the four foreground-ceaned maps, downgraded to 6 and repixeized at an N side = 32 grid. Ceary, the maps are virtuay identica on arge anguar scaes measured reative to cosmic variance, with any differences being smaer than 0.1 σ in terms of cosmoogica parameters. However, it is worth noting that the best-fit mode, (q, n) = (1, 0), is in some tension with the ow- spectrum, at about 1.7 σ in this pot. The same tension between arge and sma anguar scaes is observed in Panck Coaboration XV (2013) and Panck Coaboration XVI (2013) with higher statistica significance using the fu Panck ikeihood. Irrespective of physica interpretation, the cacuations presented here demonstrate that these ow- features are robust with respect to component separation techniques Resuts We perform the power spectrum and parameter estimation anaysis for both the data and the FFP6 simuations described in Sect. 4. The resuts for the atter are given in Appendix E. Figure 10 shows estimates of the anguar power spectrum for each foreground-ceaned map, with the uncertainties given by the Fisher matrix. The parameter summary given in Fig. 11 shows the parameter constraints derived using mutipoes between = 70 and 2000, and compares these to resuts obtained with the CamSpec and Pik ikeihoods (Panck Coaboration XV 2013). Differences in the power spectra at high are mosty absorbed by the two-parameter foreground mode, rendering consistent cosmoogica parameters. For exampe, the high- power excess seen in the Commander-Ruer map is we-fitted in terms of residua point sources, which makes intuitive sense, considering the ower anguar resoution of this map (see Sect. 5). However, the ΛCDM parameter uncertainties derived from the four codes are very consistent. This indicates that most of the cosmoogica information content above 1500 is degenerate with the extragaactic foreground mode, and a more sophisticated foreground treatment is required in order to recover significant cosmoogica parameter constraints from these scaes. Beyond this, deviations among cosmoogica parameters are sma and within 1 σ for a methods and most of the parameters. Further, the parameters derived from the four foregroundceaned CMB maps are in good agreement with those obtained by CamSpec and Pik using cross-spectra; departures are we within 1 σ for most parameters. Inspecting the differences between the best-fit modes derived from the four foreground-ceaned maps and from CamSpec potted in Fig. 12, we find that the reative residuas are within 40 µk 2 for a mutipoe ranges, and smaer than 20 µk 2 at high. This can be compared to the corresponding residuas for the FFP6 simuation shown in Appendix E. The ikeihood used for this anaysis does not take into account some systematic effects that wi affect our foregroundceaned CMB maps, such as reative caibration uncertainties 7. Higher-order statistics The foreground-ceaned CMB maps presented in this paper are used as inputs for most Panck anayses of higherorder statistics, incuding non-gaussianity studies (Panck Coaboration XXIV 2013), studies of statistica isotropy (Panck Coaboration XXIII 2013), gravitationa ensing by arge-scae structure (Panck Coaboration XVII 2013), and of the integrated Sachs-Wofe effect (Panck Coaboration XXIV 2013). In this section we provide a summary of the non- Gaussianity and gravitationa ensing resuts Non-Gaussianity Primordia non-gaussianity is typicay constrained in terms of the ampitude, fnl oca, of the quadratic corrections to the gravitationa potentia, as we as by means of the three-point correation function based on different triange configurations. The resuts from these cacuations for the foreground-ceaned CMB maps are presented in Panck Coaboration XXIV (2013). After subtraction of the ensing-isw correation contribution, the fina resut is fnl oca = 2.7±5.8, as estimated from the map using the KSW bispectrum estimator (Komatsu et a. 2005), consistent within 1 σ with resuts from other methods and foregroundceaned maps. Uncertainties are evauated by means of the FFP6 simuations, and potentia biases are studied using both Gaussian and non-gaussian CMB reaizations. In particuar, when a detectabe eve of primordia non-gaussianity ( fnl oca = ) is injected into the FFP6 simuations, each foreground-ceaned map yieded a positive detection within 2 σ of the expected vaue, recovering vaues of fnl oca = 8.8 ± 8.6, 19.0 ± 7.5, 11.1 ± 7.6 and 19.7 ± 7.4 for Commander-Ruer,,,, respectivey. We see that and demonstrate the best recovery of the injected non-gaussianity, and we favoured the atter for non-gaussian studies for its faster performance over. The foreground-ceaned CMB maps presented in this paper do not provide significant evidence of a non-zero vaue of 11

12 [L(L + 1)] 2 C φφ [L(L + 1)] 2 C φφ /2π [ ] 10 7 L /2π [ 10 ] Input C-R L Fig. 14: Lensing power spectrum estimates from FFP6 simuations using an apodized mask covering f sky, of the sky. fnl oca, and reaistic simuations show that the component separation methods do not suppress rea non-gaussian signatures within expected uncertainties. The impications of these resuts in terms of eary Universe physics are discussed in the reevant papers (Panck Coaboration XXIV 2013; Panck Coaboration XXII 2013) Gravitationa ensing by arge-scae structure Gravitationa ensing by the intervening matter imprints a non- Gaussian signature in the CMB, which aows the reconstruction of the gravitationa potentia integrated aong the ine of sight to the ast scattering surface. In Panck Coaboration XVII (2013), this effect has been detected at a high significance eve (greater than 25 σ) using the Panck temperature maps. Specificay, the ensing induced correations between the tota intensity and its gradients have been used to reconstruct a neary fu sky map of the ensing potentia φ, which has been used for further studies on Panck data, incuding the detection of a non-zero correation with the ISW (Panck Coaboration XXIV 2013; Panck Coaboration XIX 2013) and other tracers of arge-scae structure (notaby, significant correation with the CIB is reported in Panck Coaboration XVIII 2013), as we as the estimate of the power spectrum of the ensing potentia and the associated ikeihood. The atter was constructed using a simpe minimum variance combination of the 143 and 217 GHz maps on about 70 % of the sky, as we as subtracting dust contamination using the 857 GHz Panck channe as a tempate (Panck Coaboration XVII 2013). These ensing resuts have improved the cosmoogica constraints from Panck (Panck Coaboration XVI 2013). The foreground-ceaned CMB maps described in Sect. 5 were used to perform a ensing extraction on a arger sky fraction, reaching about 87 % of the sky. We found the ensing power spectrum to be in good agreement with the one obtained using the minimum variance combination, i.e., the signa agrees within 1 σ in the majority of the anguar domain bins, and is characterized by an equivaent uncertainty. The foreground-ceaned maps were further expoited on the baseine 70 % sky fraction for assessing the robustness of the main reconstruction against the foreground contamination (Panck Coaboration XVII 2013). We show that the component separation agorithms presented in this paper do not bias the ensing reconstruction in the case of the arge sky fraction considered here. We consider FFP6 simuations incuding noise and ensed CMB signa, propagated through each of the component separation agorithms described in Sect. 3. We perform a ensing potentia reconstruction in the pixe domain based on the CMB maps processed by the four component separation methods using the metis agorithm described in Panck Coaboration XVII (2013). This method uses the quadratic estimator presented in Okamoto & Hu (2003), which corrects for the mean-fied bias caused by extra sources of statistica anisotropy in addition to the CMB. For each method, we combine the masks of CO regions, nearby gaaxies and compact objects as defined in Panck Coaboration XVII (2013), with the CG90 mask described in Sect. 4. This procedure resuts in masks with sky fractions f sky = 0.836, 0.851, 0.850, for Commander-Ruer,, and, respectivey. We estimate the ensing potentia power spectrum, C φφ L, foowing the methodoogy described in Panck Coaboration XVII (2013). It consists of a pseudo-c estimate based on a highy-apodized version of the ensing potentia reconstruction, which has an effective avaiabe sky fraction f sky,2 = 0.648, 0.690, 0.686, for Commander-Ruer,, and, respectivey. The band-power reconstructions in 17 bins in the range are potted in Fig. 14, as we as the residuas reative to the theoretica ens power spectrum. A agorithms achieved an unbiased estimation of the underying ensing power spectrum, with χ 2 = 10.58, 17.34, 18.54, 15.30, for Commander-Ruer,,, and respectivey, with 17 degrees of freedom. The associated probabiity-toexceed (PTE) vaues are 83 %, 36 %, 29 %, 50 %. The power spectrum estimates are in remarkabe agreement with each other. However, the Commander-Ruer soution has significanty arger uncertainties, as expected from its ower signa-to-noise ratio to ensing due to its arger beam. These resuts on simuated foreground-ceaned CMB maps demonstrate that the component separation agorithms do not ater the ensing signa, and this provides a strategy for achieving a robust ensing reconstruction on the argest possibe sky coverage. The foreground-ceaned maps have been used in Panck Coaboration XVII (2013) to obtain ensing potentia estimates on 87 % of the sky. 8. Foreground components In this section we consider the diffuse Gaactic components, and present fu-sky maps of therma dust and CO emission, as we as a singe ow-frequency component map representing the sum of synchrotron, AME, and free-free emission. Our a-sky CO map is a type 3 product as presented in Panck Coaboration XIII (2013). To assess the accuracy of these maps, we once again take advantage of the FFP6 simuation. The Commander-Ruer method used in the foowing is described in Appendix A and consists of a standard parametric Bayesian MCMC anaysis at 12

13 Panck Coaboration: Panck 2013 resuts. XII. Component separation 0 µk (a) Low-frequency index (a) Low-frequency component ampitude at 30 GHz 0 µk km s 1 5 (b) CO ampitude at 100 GHz (b) Dust emissivity Fig. 16: Posterior mean spectra parameter maps derived from the ow-resoution anaysis. The top pane shows the power aw index of the ow-frequency component, and the bottom pane shows the emissivity index of the one-component therma dust mode. Note that the systematic error due to monopoe and dipoe uncertainties is significant for the dust emissivity in regions with a ow therma dust ampitude. 0.0 MJy sr (c) Therma dust ampitude at 353 GHz Fig. 15: Posterior mean foreground ampitude maps derived from the ow-resoution anaysis. From top to bottom are shown the ow-frequency, CO and therma dust emission maps. ow anguar resoution, foowed by a generaized east-squares soution for component ampitudes at high resoution Data seection and processing We ony use the seven owest Panck frequencies, from 30 to 353 GHz. The two highest channes have significanty different systematic properties than the ower frequency bands, for instance concerning caibration, ZLE, and noise correations, and they are more reevant to therma dust and CIB studies than to the present CMB anaysis. Studies of specific foregrounds (CO, therma dust, CIB etc.) using a Panck frequencies as we as anciary data are discussed in companion papers Panck Coaboration VI (2013), Panck Coaboration XXI (2013), Panck Coaboration (2013b), and additiona future pubications wi consider extensions to AME, synchrotron and freefree emission. In order to obtain unbiased estimates of the spectra parameters across a frequency bands, each map is downgraded from its native resoution to a common anguar resoution of 400 and repixeized at Nside = 256, a imit imposed by the LFI 30 GHz channe. Once the spectra indices have been determined, we reestimate the component ampitudes at native Panck resoution (see Appendix A). Athough the smoothing operation introduces noise correations between pixes, we mode the noise of the smoothed maps as uncorreated white noise with an effective standard deviation, σ(p), for each pixe p. This approximation does not bias the fina soution, because the anaysis is performed independenty for each pixe. However, it is important to note that correations between pixes are not taken into account in this anaysis. The ef13

14 0 14 Fig. 17: χ 2 per pixe for the joint CMB and foreground anaysis. The expected vaue for an acceptabe fit is 7, corresponding to the number of frequency bands used in this anaysis. The pixes with high vaues can be cassified into two types, due to either modeing errors (i.e., high residuas in the Gaactic pane) or to un-modeed correated noise (i.e., stripes crossing through ow dust emission regions). fective noise uncertainty, σ(p), is estimated using reaistic noise simuations downgraded in the same way as the data. The measured instrumenta bandpasses are taken into account by integrating the emission aws over the bandpass for each component at each Monte Caro step in the anaysis. The monopoe (zero-point) of each frequency map is not constrained by Panck, but is rather determined by postprocessing, and associated with a non-negigibe uncertainty (see Tabe 5 of Panck Coaboration I 2013). In addition, each frequency map incudes a significant monopoe contribution from isotropic extragaactic sources and CIB fuctuations not traced by oca Gaactic structure, ranging from ess than about µk at 70 GHz to severa hundreds of µk at 353 GHz. Finay, the effective dipoe in each map is associated with significant uncertainty due to the arge kinematic CMB dipoe. In order to prevent these effects from introducing modeing errors during component separation, they must be fit either prior to or jointy with the Gaactic parameters. Unfortunatey, when aowing free spectra parameters per pixe, there is a near-perfect degeneracy among the offsets, the foreground ampitudes and the spectra indices, and in order to break this degeneracy, it is necessary to reduce the number of spectra degrees of freedom. We adopt the method described by Wehus et a. (2013) for this purpose, which has the additiona advantage of making minima assumptions about the foreground spectra. In short, this method uses inear regression between data from CMBsubtracted maps evauated on pixes faing within each arge N side = 8 pixe to estimate the reative offsets, m 1 and m 2, between any two maps at each position on the sky. Each regression provides a constraint of the form m 1 = am 2 + b, where a and b are the sope and offset, respectivey, and where each vaue of m i consists of the sum of both a monopoe and a dipoe term evauated at that position. The individua monopoes and dipoes can then be reconstructed by measuring a and b in different regions of the sky, expoiting spatia variations in spectra indices, and soving jointy for two monopoes and dipoes, incuding constraints from a positions. To minimize degeneracies, a positivity prior is imposed on the fit, such that statisticay significant negative pixes are heaviy penaized. For 44 and 70 GHz, we retain the dipoe vaues determined during the mapmaking process, and do not attempt to fit them. The resuting compete set of monopoe and dipoe vaues is isted in Tabe 2. As a cross-check, we performed a dedicated Commander run in which we fitted for the dipoe at 353 GHz, together with the foreground ampitudes and spectra indices, and ony found sub- µk differences. This channe is by far the most probematic in our data set in terms of offset determination, because of the very bright dust emission at this frequency. As a resut, there is a arge reative uncertainty between the zero-eve of the dust ampitude map and the 353 GHz channe offset not accounted for in the foowing anayses. However, the sum of the two terms is we determined, and a potentia error in either therefore does not compromise the quaity of the other signa components (e.g., CMB and ow-frequency components). A simiar comment appies between the offset at 30 GHz and the zero-eve of the ow-frequency component, athough at a significanty ower eve Component modes and priors Our mode for the ow-resoution CMB anaysis incudes four independent physica components: CMB; ow-frequency emission; CO emission; and therma dust emission. It can be written schematicay in the form ( ) βf (p) ν s ν (p) = A CMB (p) + A f (p) + ν 0,f + A CO (p) f ν,co + A d (p) e hν0,d kt d (p) 1 e hν kt d (p) 1 ( ν ν 0,d ) βd (p)+1, (1) where A i (p) denotes the signa ampitude for component i at pixe p, ν 0,i is the reference frequency for each component, and ν refers to frequency. (Note that for readabiity, integration over bandpass, as we as unit conversions between antenna, fux density and thermodynamic units, is suppressed in this expression.) Thus, each component is modeed with a simpe frequency spectrum parameterized in terms of an ampitude and a sma set of free spectra parameters (a power-aw index for the owfrequency component, and an emissivity index and temperature for the therma dust component); no spatia priors are imposed. One goa of the present anaysis is to understand how we this simpe mode captures the sky signa in terms of effective components over the considered frequency range, and we expoit the FFP6 simuation (see Sect. 4) for this purpose. In order to take into account the effect of bandpass integration, each term in the above mode is evauated as an integra over the bandpass as described in Sect. 3 of Panck Coaboration IX (2013), and converted internay to thermodynamic units. Accordingy, the reference frequencies in Eq. 1 are computed as effective integras over the bandpass, such that the ampitude map, A i (p), corresponds to the foreground map observed by the reference detector, i.e., after taking into account the bandpass. In order to minimize degeneracies between the different signa components, the reference band for a given component is set to the frequency at which its reative signa-to-noise ratio is maximized. The foreground mode defined in Eq. 1 is motivated by prior knowedge about the foreground composition over the CMB frequencies as outined in Sect. 2, as is our choice of priors. In 14

15 Panck Coaboration: Panck 2013 resuts. XII. Component separation (a) CMB ampitude 25 0 (b) CO ampitude at 100 GHz 3 (c) Therma dust ampitude at 353 GHz Fig. 18: Comparison of the high-resoution Ruer (top) and ow-resoution Commander (bottom) ampitude maps for a particuary strong CO compex near the Fan region; the maps are centred on Gaactic coordinates (, b) = (110, 15 ), and the grid spacing is 5. Coumns show, eft to right: the CMB ampitude; the CO ampitude at 100 GHz and the therma dust ampitude at 353 GHz. Tabe 2: Estimated monopoes and dipoes in Gaactic coordinates, a measured in thermodynamic µk. Errors are estimated by bootstrapping, and do not account for correated errors across frequencies. In particuar, the 353 GHz monopoe uncertainty is dominated by systematic errors not incuded in these estimates. Note that the dipoes at 44 and 70 GHz are fixed at the vaues determined in the mapmaking. Frequency [ GHz] Monopoe [µk] 8±2 2±1 15 ± 1 15 ± 1 33 ± 1 86 ± ± 4 addition to the Jeffreys prior5 (Eriksen et a. 2008), we adopt Gaussian priors on a spectra parameters with centre vaues and widths attempting to strike a baance between prior knowedge and aowing the data to find the optima soution. Where needed, we have aso run dedicated anayses, either incuding particuar high signa-to-noise ratio subsets of the data or using a ower res5 The purpose of the Jeffreys prior is to normaize the parameter voume reative to the ikeihood, such that the ikeihood becomes socaed data-transated, i.e., invariant under re-parameterizations. X dipoe [µk] 4 ± 3 0±0 0±0 2±1 2±1 2±1 11 ± 10 Y dipoe [µk] 6 ± 2 0±0 0±0 5±1 7±1 11 ± 2 52 ± 12 Z dipoe [µk] 6±1 0±0 0±0 5 ± 1 6 ± 1 10 ± 2 37 ± 8 oution parameterization to increase the effective signa-to-noise in order to inform our prior choices. We now consider each foreground component in turn, and note in passing that the CMB component, by virtue of being a backbody signa, is given by a constant in thermodynamic temperature units. We approximate the ow-frequency component by a straight power aw in antenna temperature with a free spectra index per pixe, and adopt a prior of β = 3 ± 0.3 (this is the index in terms of brightness temperature). This choice is determined by 15

16 noting that the prior is in practice ony reevant at high Gaactic atitudes where the signa-to-noise ratio is ow and the dominant foreground component is expected to be synchrotron emission; in the signa-dominated and ow-atitude AME and freefree regions, the data are sufficienty strong to render the prior irreevant. For vaidation purposes, we have aso considered minor variations around this prior, such as β = 2.9 ± 0.3 and β = 3.05 ± 0.2, finding ony sma differences in the fina soutions. The reference band for the ow-frequency component is set to 30 GHz, where the ow-frequency foreground signa peaks. The fina ow-frequency ampitude map is provided in units of thermodynamic microkevin. The CO emission is modeed in terms of a singe ine ratio for each frequency. Specificay, the CO ampitude is normaized to the 100 GHz band, and defined in units of µk km s 1 (Panck Coaboration XIII 2013). The ampitude at other frequencies is determined by a singe mutipicative factor reative to this, with a numerica vaue of at 217 GHz and at 353 GHz; a other frequencies are set to zero. These vaues are obtained from a dedicated CO anaysis that incudes ony high signa-tonoise ratio CO regions covering a tota of 0.5% of the sky. The derived vaues are in good agreement with those presented by Panck Coaboration XIII (2013). Therma dust emission is modeed by a one-component modified backbody emission aw with a free emissivity spectra index, β d, and dust temperature, T d, per pixe. However, since we ony incude frequencies beow 353 GHz, the dust temperature is argey unconstrained in our fits, and we therefore impose a tight prior around the commony accepted mean vaue of T d = 18 ± 0.05 K. The ony reason we do not fix it competey to 18 K is to aow for modeing errors near the Gaactic centre. The dust emissivity prior is set to β d = 1.5±0.3, where the mean is determined by a dedicated run fitting for a singe best-fit vaue for the high-atitude sky, where the prior is reevant. The reference band for the therma dust component is 353 GHz, and the fina map is provided in units of megajansky per steradian. 30 µk 30 (a) Low-frequency component residua at 30 GHz 4 µk 4 (b) CO residua at 100 GHz 8.3. Resuts and vaidation The output of the Bayesian component separation agorithm is a set of sampes drawn from the joint posterior distribution of the mode parameters, as opposed to a singe we-defined vaue for each. For convenience, we summarize this distribution in terms of posterior mean and standard deviation maps, computed over the sampe set, after rejecting a short burn-in phase. The goodness-of-fit is monitored in terms of the χ 2 per pixe. Athough convenient, it is, however, important to note that this description does not provide a comprehensive statistica representation of the fu posterior distribution, which is intrinsicay non-gaussian. One shoud be carefu about making inferences in the ow signa-to-noise regime based on this simpified description. The ow-resoution Commander posterior mean ampitude maps are shown in Fig. 15 for the ow-frequency, CO, and therma dust components, and the spectra index maps in Fig. 16. The associated χ 2 map is potted in Fig. 17. Note that because we are samping from the posterior instead of searching for the maximum-ikeihood point, the expected number of degrees of freedom is equa to N band = 7 in this pot, not N band N par. Figure 18 compares the high-resoution Ruer soution to the ow-resoution Commander soution for CMB, CO and therma dust on a particuary strong CO compex near the Fan region, centred on Gaactic coordinates (, b) = (110, 15 ). 100 µk 100 (c) Therma dust residua at 353 GHz Fig. 19: Ampitude residua maps, A out A in, computed bindy from the FFP6 simuation. The panes show (from top to bottom) the ow-frequency residua at 30 GHz, the CO residua at 100 GHz and the therma dust residua at 353 GHz. A units are thermodynamic µk. The white ines indicate the boundary of the Commander ikeihood anaysis mask, removing 13% of the sky. Severa features can be seen here, foremost of which is that the Gaactic pane is strikingy obvious, with χ 2 vaues exceeding 10 4 for seven degrees-of-freedom in a few pixes. This is not surprising, given the very simpified mode at ow frequencies (i.e., a singe power aw accounting for AME, synchrotron, and free-free emission), as we as the assumption of a neary 16

17 constant dust temperature of 18 K. Second, there is an extended region of moderatey high χ 2 roughy aigned with a great circe going through the Eciptic South Poe, indicating the presence of correated noise in the scanning rings not accounted for in our white noise mode. Based on these, and other considerations, it is cear that parts of the sky must be masked before proceeding to CMB power spectrum and ikeihood anayses. This masking process is discussed at greater ength in Panck Coaboration XV (2013), and resuts in different masks for specific appications. The goa of our present discussion is to evauate the adequacy of the mask adopted for ow- ikeihood anaysis (L87), which is based on the fits presented here. This mask removes 13% of the sky, and is derived from a combination of χ 2 and component ampitude threshoding. For vaidation purposes, we anayse the simuations described in Section 4 in the same way as the rea data, incuding monopoe and dipoe determination, CO ine ratio estimation and spectra index estimation. Individua component maps at each observed frequency are avaiabe from the simuation process, and used for direct comparison with the reconstructed products. In Fig. 19 we show the differences between the recovered and input component maps at their respective reference frequencies. The boundary of the 13% Commander mask is traced by the white contours, and a best-fit monopoe and dipoe have been subtracted from each difference map. A difference maps are shown in units of thermodynamic µk. The top pane of Fig. 20 gives the error histograms outside the masked region for each component, normaized to the respective estimated standard deviation; if the recovered soution has both correct mean and standard deviation, these histograms shoud match a Gaussian distribution with zero mean and unit variance, indicated by the dashed back ine. Conversey, a significant bias woud be visibe as a horizonta shift in this pot, whie under-estimation of the errors woud resut in too wide a distribution and vice versa. The bottom pane shows the fractiona error (i.e., the error divided by the true input vaue) for a pixes with signa above 5 σ; the fractiona error is not a usefu quantity for noisier signas. The difference maps in Fig. 19 dispay significant errors in the Gaactic pane. For the ow-frequency component, the residuas are dominated by free-free emission, whie for therma dust the dominant contaminant is CO emission. However, outside the mask the residuas are sma, and, at east for the ow-frequency and CO components, the spatia characteristics appear simiar to instrumenta noise. This is more cear in the histograms shown in the top pane of Fig. 20; the mean and standard deviations are δ f = 0.01±1.12, δ CO = 0.00±0.87, and δ td = 0.00±2.01, respectivey, for the ow-frequency, CO and therma dust components. There is no evidence of bias outside the mask in any component, and the error estimates are accurate to 12 and 13 % for the owfrequency and CO components. Note, though, that the estimated error for the CO component is actuay arger than the true uncertainty, suggesting that the white noise approximation for the 100 GHz channe overestimates the true noise. This can occur if the correated instrumenta noise is important in regions where there is no significant CO emission. Locay re-scaing the white noise to account for spatiay varying correated noise woud correct this effect. For the therma dust component, on the other hand, the error is underestimated by a factor of 2. The expanation for this is most easiy seen from the ower pane of Fig. 19. This map is dominated by isotropic CIB fuctuations, rather than instrumenta noise. Because these fuctuations have a sighty differ- Panck Coaboration: Panck 2013 resuts. XII. Component separation Probabiity distribution Probabiity distribution Low-freq CO Therma dust Expected gaussian Normaized error [σ] Low-freq CO Therma dust Fractiona error Fig. 20: Error vaidation for component ampitudes, evauated from the FFP6 simuation. The upper pane shows histograms of the normaized errors δ = (A out A in )/σ out for the three foreground components and incuding a pixes outside the Commander ikeihood anaysis mask. The ower pane shows histograms of the fractiona error f (A out A in )/A in for pixes with a foreground detection eve above 5 σ. No evidence of significant bias is observed for any component, and the uncertainty estimates for the ow-frequency and CO components are accurate to about 12 %; the therma dust uncertainty is underestimated by a factor of 2 due to the presence of unmodeed fuctuations. ent spectrum than the dominant Gaactic dust emission, and the mode does not account for a separate CIB component, the error on the Gaactic component is underestimated. When using the Gaactic map presented here for detaied anaysis near the noise imit, taking into account these residua fuctuations is essentia, and the effective noise per pixe shoud be increased by a factor of 2. As ceary seen in Fig. 19, the residuas inside the mask are highy significant in a strict statistica sense. However, as seen in the bottom pane of Fig. 20, they are reativey sma in terms of fractiona errors. Specificay, the three histograms have means 17

18 Probabiity distribution GHz 70 GHz 217 GHz 353 GHz Therma 143 GHz Therma 217 GHz Expected Gaussian Normaized error [σ] MJy sr 1 Fig. 21: Vaidation of spectra parameters for ow-frequency foregrounds, therma dust, and CO emission, evauated from the FFP6 simuation. Each histogram shows the error distribution at the two eading sub-dominant frequencies in the form of the normaized errors δ = (A out (ν) A in (ν))/σ out (ν) for a pixes outside the Commander ikeihood anaysis mask, where A out (ν) is the predicted foreground ampitude at frequency ν given the estimated ampitude and spectra parameters, and σ out (ν) is the corresponding standard deviation computed over the sampe set. and standard deviations of f f = 0.00 ± 0.10, f CO = 0.03 ± 0.10, and f td = 0.00 ± 0.06, respectivey, for the ow-frequency, CO and therma dust components. The argest bias is observed for the CO component, for which the absoute ampitude is biased by 3 %. The bias in the ow-frequency and therma dust components is negigibe, and the fractiona uncertainties are 10 and 6 %, respectivey. This confirms that approximating the sum of the three ow-frequency components by a singe power-aw over the Panck frequency bands is reasonabe; if modeing errors dominated, one woud expect to see a significant bias in the resuting ampitude. In order to vaidate the spectra parameters, we show in Fig. 21 histograms of the normaized residuas for each foreground component evauated at its two eading sub-dominant frequencies (i.e., at 44 and 70 GHz for the ow-frequency component; at 217 and 353 GHz for the CO component; and at 143 and 217 GHz for the therma dust component). The means and standard deviations of these distributions are: δ f (44 GHz) = 0.41 ± 1.98 and δ f (70 GHz) = 0.34 ± 2.04 for the owfrequency component; δ CO (217 GHz) = 0.10 ± 0.84 and δ CO (353 GHz) = 0.51 ± 1.00 for the CO component; and δ td (143 GHz) = 0.02 ± 1.53 and δ td (217 GHz) = 0.13 ± 1.87 for the therma dust component. As expected, the effect of modeing errors is more significant at the sub-dominant frequencies than at the pivot frequencies, when measured in terms of statistica uncertainties, since the foreground signa is weaker and the confusion with the other components reativey arger. Nevertheess, we see that the absoute bias is at most 0.5 σ for the CO component at 353 GHz, whie the therma dust bias is negigibe even at 143 GHz. The estimated uncertainties are generay underestimated by up to a factor of two due to these modeing errors. Fig. 22: Top: Difference map between the estimated therma dust ampitude at 353 GHz derived by Panck Coaboration (2013b) and the ow-resoution dust map presented here, both smoothed to 40. A monopoe and dipoe term outside the Commander mask has been removed. The former incudes ony high-frequency observations (the Panck 353, 545 and 857 GHz channes and observations at 100 µm), whie the one derived here ony uses ow-frequency data (Panck GHz). Note that the coour scae ranges between 0.02 and 0.02 MJy sr 1 in this pot, whereas it ranges from 0 to 2.5 MJy sr 1 in the bottom pane of Fig. 15. Masked regions indicate pixes with a CO ampitude at 100 GHz arger than 1 µk km s 1. Bottom: T T pot between the same two maps. Finay, the efficiency of the adopted foreground mode for CMB anaysis is quantified in Appendix E in terms of power spectrum residuas and cosmoogica parameter estimation. To summarize, we find that the simpified mode, defined by Eq. 1, provides a good fit to the reaistic FFP6 simuation for most of the sky. Absoute residuas are sma, and the ampitude uncertainty estimates are accurate to around 12 %, except for the therma dust component for which unmodeed CIB fuctuations 18

19 are important. Further, we find that the rea Panck data behave both quaitativey and quantitativey very simiary to the FFP6 simuation, suggesting that this approach aso performs we on the rea sky. Panck Coaboration: Panck 2013 resuts. XII. Component separation 8.4. Interpretation and comparison with other resuts The maps shown in Figs. 15 and 16 provide a succinct summary of the average foreground properties over the Panck frequency range. We now consider their physica interpretation and compare them to products from aternative methods. First, the top pane of Fig. 22 shows a difference map between the dust map at 353 GHz derived in the present paper and one determined from ony the three highest Panck frequencies and the 100 µm IRIS map by Panck Coaboration (2013b), shown in fux density units, after removing a sma monopoe and dipoe difference. This map is to be compared to the corresponding dust ampitude map in the bottom pane of Figure 15; note that the coour range varies from 0.02 to 0.02 MJy sr 1 in the difference pot, and between 0 and 2.5 MJy sr 1 in the ampitude pot. Thus, despite the very different data sets and methods, we see that the two reconstructions agree to very high accuracy outside the Gaactic pane. Inside the Gaactic pane, the differences are dominated by residuas due to different CO modeing, seen as soid bue coours in Fig. 22; however, even in this region the differences are smaer than 5 % of the ampitude. The bottom pane shows a corresponding T T pot between the two maps, excuding any pixe for which the Commander CO ampitude at 100 GHz is arger than 1 µk km s 1. The two maps agree to 0.2% in terms of best-fit ampitude. From Fig. 16, we see that the dust emissivity ranges between 1.3 and 1.7 for most of the sky; considering ony the pixes with a posterior distribution width that is a third of the prior width (i.e., σ(β d ) < 0.1), we find a mean vaue of The two exceptions are a arge region of shaow indices northeast of the Gaactic centre, and steep indices near the Gaactic pane. The former region corresponds to a part of the sky with ow dust emission, where we expect the spectra index to be sensitive to both monopoe and dipoe residuas, as we as instrumenta systematics, such as correated 1/ f noise. The atter appears to be particuary pertinent here because the shaow index region at east partiay traces the Panck scanning strategy; as a resut, the systematic error on the spectra index in this region is considerabe. The main systematic uncertainty connected to the region of steep indices around the Gaactic pane is confusion with CO emission. The CO map shown in Figs. 15 and 18 is discussed in greater detai in Panck Coaboration XIII (2013). A distinct advantage of this particuar soution over avaiabe aternatives is its high signa-to-noise ratio per pixe, which is achieved by reducing a information into a singe vaue per pixe. Consequenty, this map serves as a unique too for foow-up CO observations. However, the assumption of a constant ine ratio over the fu sky may ead to a significant systematic uncertainty on CO ampitude per pixe. This possibiity is investigated in a forthcoming work (Panck Coaboration 2013a). Finay, the spectra index map for the ow-frequency component shown in Fig. 16 can be used to determine the dominant ow-frequency component (synchrotron, free-free or AME) as a function of position on the sky. To iustrate this connection, we once again take advantage of the FFP6 simuation for which we know the ampitude of each ow-frequency component per pixe. In the top pane of Fig. 23, we use this information to make a dominant component map ; dark bue indicates (a) Dominant ow-frequency component map (b) Low-frequency component power-aw index Fig. 23: Top: Dominant foreground component per pixe at 30 GHz in the FFP6 simuation. Dark bue indicates that synchrotron emission is the strongest component at 30 GHz, ight bue indicates that free-free dominates, and orange indicates that spinning dust (AME) is the strongest component. Bottom: The recovered ow-frequency power-aw index derived from the same simuation. that synchrotron emission is strongest at a given pixe, ight bue that free-free is strongest, and orange that spinning dust (AME) dominates. In the bottom pane, we show our derived power-aw index map from the same simuation. As expected, the correspondence between the power-aw index and the dominant component is very strong, impying that the spectra index map can be used to trace the individua components. In particuar, we see that an index beow about 3.3 refects the presence of a component consistent with a spinning dust mode peaking beow 20 GHz over the Panck frequency range, whie an index higher than around 2.3 signas the importance of free-free emission. Intermediate vaues typicay indicate synchrotron emission, athough it shoud be noted that the signa-to-noise ratio at very high atitudes is ow and the resuts are therefore prior-driven in these regions. Returning to the spectra index map shown in Fig. 16, we see a good correspondence between the rea data and the simuation. Features present in the simuation aso appear in the data. For instance, we see that the spectra index in the so-caed Fan Region (i.e., near Gaactic coordinates (, b) = (90, 20 )) is ow in both cases, and this aone provides strong evidence for the presence of AME. Further, the AME spectra index is consistent with the spinning dust interpretation. This power-aw index map may be used to identify particuar AME regions for foow-up observa- 19

20 tions. Finay, we note that regions known for strong free-free emission, such as the Gum Nebua or Zeta Ophiuchi, have spectra indices cose to 2.1 or 2.2, as expected. 9. Concusions We have tested four component separation agorithms on the Panck frequency maps to produce cean maps of the CMB anisotropies over a arge area of sky. These CMB maps are used for studies of statistics and isotropy (Panck Coaboration XXIII 2013), primordia non-gaussianity (Panck Coaboration XXIV 2013), gravitationa ensing (Panck Coaboration XVII 2013), the ISW effect (Panck Coaboration XIX 2013), cosmic geometry and topoogy (Panck Coaboration XXVI 2013), searches for cosmic defects from primordia phase transitions (Panck Coaboration XXV 2013), as we as an integra part of the ow- Panck ikeihood (Panck Coaboration XV 2013). Two of the methods, one using interna foreground tempates () and the other an ILC in needet space (), are non-parametric, extracting the CMB map by minimizing the variance of the tota contamination. The other two methods fit modes of the foregrounds to cean the CMB of their emission. One fits a parametric mode in rea space (C-R) and one fits a non-parametric in the harmonic domain (). A four methods have been demonstrated to work we both on rea and simuated data, and to yied consistent resuts. Nevertheess, there are differences between the methods, making them more or ess suitabe for specific appications. For instance, Commander-Ruer aows a joint parametric foreground estimation and CMB power spectrum estimation, with fu propagation of foreground uncertainties to cosmoogica parameters, but is imited to a ower anguar resoution than the other codes. This method has therefore been seected for the ow- Panck ikeihood (Panck Coaboration XV 2013) and to produce astrophysica component maps (Sect. 8), whie it is sub-optima for appications requiring fu anguar resoution, e.g., gravitationa ensing reconstruction or estimation of primordia non-gaussianity. For these purposes, we use the three higher-resoution maps. We take to be the eading method, based on its superior performance on the FFP6 simuation, where it has be shown to have the owest residua foreground contamination at arge scaes and to preserve primordia non-gaussianity. When subjecting foreground-ceaned Panck maps to scientific anaysis, we use the other two or three maps, as appropriate, to assess the uncertainties inherent in the choice of methods and the assumptions they make. Indeed, this is the main purpose for presenting four different CMB soutions to the genera community. The CMB anisotropies are robusty recovered over a arge fraction (73 %) of the sky and down to sma anguar scaes, reaching to mutipoes We characterize the CMB maps with anguar power spectra and cosmoogica parameter constraints. Parameter constraints from these maps are consistent with those from the Panck ikeihood function based on crossspectra and arge sky cuts (Panck Coaboration XV 2013). This agreement supports the robustness of both our component separation methodoogy and cosmoogica parameter constraints. The rea-space parametric fits of Commander-Ruer enabe us to characterize the diffuse Gaactic foregrounds. We parameterize them with a ow-frequency power-aw component, representing the sum of synchrotron, free-free, and AME emission, a high-frequency modified backbody spectrum describing therma dust emission, and a moecuar CO component. Using ony the Panck data from 30 to 353 GHz, we fit for the ampitude and spectra parameters of the three foregrounds and the CMB simutaneousy at each pixe of a 40-arcmin resoution map. The spectra parameters are the ow-frequency component power-aw exponent and the modified backbody emissivity power-aw exponent; the CO ine ratios are spatiay fixed. These parameters give us the source mixing matrix, which we then use in a direct inversion to deduce the component ampitudes at higher resoution. Through Gibbs samping, we obtain reaizations drawn from the fu posterior distribution of possibe foreground and CMB soutions, giving us a powerfu abiity to statisticay characterize our resuts. Our in-depth anaysis of the recovered CMB anisotropies is unprecedented for component separation studies, concerning both the accuracy of cosmoogica parameter constraints, and studies of eary Universe physics and structure formation through gravitationa ensing. On the other hand, the compex nature of the foreground emission over such a arge frequency range imits us to the use of reativey simpe methods when anaysing Panck data aone. An extensive study in combination with other probes of Gaactic foregrounds wi be presented in forthcoming papers. In particuar, the separation of individua components at ow frequencies requires the use of anciary data, for exampe, from the Wikinson Microwave Anisotropy Probe (WMAP) and radio surveys. Acknowedgements. The deveopment of Panck has been supported by: ESA; CNES and CNRS/INSU-IN2P3-INP (France); ASI, CNR, and INAF (Itay); NASA and DoE (USA); STFC and UKSA (UK); CSIC, MICINN, JA and RES (Spain); Tekes, AoF and CSC (Finand); DLR and MPG (Germany); CSA (Canada); DTU Space (Denmark); SER/SSO (Switzerand); RCN (Norway); SFI (Ireand); FCT/MCTES (Portuga); PRACE (EU). A description of the Panck Coaboration and a ist of its members, incuding the technica or scientific activities in which they have been invoved, can be found at project=panck&page=panck_coaboration. The authors acknowedge the support provided by the Advanced Computing and e-science team at IFCA. This work made use of the COSMOS supercomputer, part of the STFC DiRAC HPC Faciity. Some of the resuts in this paper have been derived using the HEALPix package. References Ai-Haïmoud, Y., Hirata, C. M., & Dickinson, C. 2009, MNRAS, 395, 1055 Banday, A. J., Dickinson, C., Davies, R. D., Davis, R. J., & Górski, K. M. 2003, MNRAS, 345, 897 Bennett, C. L., Hi, R. S., Hinshaw, G., et a. 2003, ApJS, 148, 97 Benoit-Lévy, A., Décheette, T., Benabed, K., et a. 2013, A&A, accepted Betoue, M., Pierpaoi, E., Deabrouie, J., Le Jeune, M., & Cardoso, J.-F. 2009, A&A, 503, 691 Bonadi, A., Ricciardi, S., Leach, S., et a. 2007, MNRAS, 382, 1791 Broadbent, A., Osborne, J. L., & Hasam, C. G. T. 1989, MNRAS, 237, 381 Cardoso, J.-F., Martin, M., Deabrouie, J., Betoue, M., & Patanchon, G. 2008, IEEE Journa of Seected Topics in Signa Processing, 2, 735 Casaponsa, B., Barreiro, R. B., Curto, A., Martínez-Gonzáez, E., & Vieva, P. 2011, MNRAS, 411, 2019 Davies, R. D., Dickinson, C., Banday, A. J., et a. 2006, MNRAS, 370, 1125 Davies, R. D., Watson, R. A., & Gutierrez, C. M. 1996, MNRAS, 278, 925 de Oiveira-Costa, A., Tegmark, M., Davies, R. D., et a. 2004, ApJ, 606, L89 Deabrouie, J., Betoue, M., Mein, J.-B., et a. 2012, ArXiv e-prints Deabrouie, J., Cardoso, J.-F., Le Jeune, M., et a. 2009, A&A, 493, 835 Dickinson, C., Davies, R. D., & Davis, R. J. 2003, MNRAS, 341, 369 Dober, G. & Finkbeiner, D. P. 2008, ApJ, 680, 1222 Draine, B. T. & Lazarian, A. 1998, ApJ, 508, 157 Eriksen, H. K., Dickinson, C., Lawrence, C. R., et a. 2006, ApJ, 641, 665 Eriksen, H. K., Huey, G., Saha, R., et a. 2007, ApJ, 656, 641 Eriksen, H. K., Jewe, J. B., Dickinson, C., et a. 2008, ApJ, 676, 10 Eriksen, H. K., O Dwyer, I. J., Jewe, J. B., et a. 2004, ApJS, 155, 227 Fendt, W. & Wandet, B. 2008, in APS Apri Meeting Abstracts, J8004 Fernández-Cobos, R., Vieva, P., Barreiro, R. B., & Martínez-Gonzáez, E. 2012, MNRAS, 420, 2162 Finkbeiner, D. P. 2004, ApJ, 614, 186 Finkbeiner, D. P., Davis, M., & Schege, D. J. 1999, ApJ, 524, 867 Finkbeiner, D. P., Langston, G. I., & Minter, A. H. 2004, ApJ, 617,

21 Ghosh, T., Banday, A. J., Jaffe, T., et a. 2012, MNRAS, 422, 3617 God, B., Odegard, N., Weiand, J. L., et a. 2011, ApJS, 192, 15 Górski, K. M., Hivon, E., Banday, A. J., et a. 2005, ApJ, 622, 759 Hasam, C., Stoffe, H., Sater, C. J., & Wison, W. E. 1982, Astronomy and Astrophysics Suppement Series, 47, 1 Hivon, E., Górski, K. M., Netterfied, C. B., P., C. B., & et a. 2002, Astrophys.J., 567 Hoang, T. & Lazarian, A. 2012, Advances in Astronomy, 2012 Jewe, J., Levin, S., & Anderson, C. H. 2004, ApJ, 609, 1 Kogut, A. 1996, in Buetin of the American Astronomica Society, Vo. 28, American Astronomica Society Meeting Abstracts, 1295 Komatsu, E., Sperge, D. N., & Wandet, B. D. 2005, ApJ, 634, 14 Lagache, G. 2003, A&A, 405, 813 Leach, S. M., Cardoso, J., Baccigaupi, C., et a. 2008, A&A, 491, 597 Leitch, E. M., Readhead, A. C. S., Pearson, T. J., & Myers, S. T. 1997, ApJ, 486, L23 Martınez-Gonzaez, E., Diego, J. M., Vieva, P., & Sik, J. 2003, MNRAS, 345 Mivie-Deschênes, M.-A., Ysard, N., Lavabre, A., et a. 2008, A&A, 490, 1093 Okamoto, T. & Hu, W. 2003, Phys. Rev., D67, Patanchon, G., Cardoso, J.-F., Deabrouie, J., & Vieva, P. 2005, MNRAS, 364, 1185 Pietrobon, D., Górski, K. M., Bartett, J., et a. 2012, ApJ, 755, 69 Panck Coaboration. 2012, Panck intermediate resuts: Dust emission at miimetre waveenghts in the Gaactic Pane (in preparation) Panck Coaboration. 2013a, In preparation Panck Coaboration. 2013b, Submitted to A&A Panck Coaboration ES. 2013, The Expanatory Suppement to the Panck 2013 resuts (ESA) Panck Coaboration I. 2013, Submitted to A&A Panck Coaboration II. 2013, Submitted to A&A Panck Coaboration Int. IX. 2013, Submitted to A&A Panck Coaboration Int. XII. 2013, Submitted to A&A Panck Coaboration IX. 2013, Submitted to A&A Panck Coaboration VI. 2013, Submitted to A&A Panck Coaboration VIII. 2013, Submitted to A&A Panck Coaboration X. 2013, Submitted to A&A Panck Coaboration XIII. 2011, A&A, 536, A13 Panck Coaboration XIII. 2013, Submitted to A&A Panck Coaboration XIX. 2013, Submitted to A&A Panck Coaboration XV. 2013, Submitted to A&A Panck Coaboration XVI. 2013, Submitted to A&A Panck Coaboration XVII. 2013, Submitted to A&A Panck Coaboration XVIII. 2011, A&A, 536, A18 Panck Coaboration XVIII. 2013, Submitted to A&A Panck Coaboration XX. 2011, A&A, 536, A20 Panck Coaboration XXI. 2013, Submitted to A&A Panck Coaboration XXII. 2013, Submitted to A&A Panck Coaboration XXIII. 2013, Submitted to A&A Panck Coaboration XXIV. 2013, Submitted to A&A Panck Coaboration XXIX. 2013, Submitted to A&A Panck Coaboration XXV. 2013, Submitted to A&A Panck Coaboration XXVI. 2013, Submitted to A&A Panck Coaboration XXVIII. 2013, Submitted to A&A Patania, P., Burigana, C., Maino, D., et a. 2003, A&A, 410, 847 Reich, P. & Reich, W. 1988, A&AS, 74, 7 Rocha, G., Contadi, C. R., Bond, J. R., & Gorski, K. M. 2009, ArXiv e-prints Rocha, G., Contadi, C. R., Bond, J. R., & Gorski, K. M. 2011, MNRAS, 414, 823R Rocha, G., Contadi, C. R., Coombo, L. P. L., et a. 2010, ArXiv e-prints Stompor, R., Leach, S., Stivoi, F., & Baccigaupi, C. 2009, MNRAS, 392, 216 Tegmark, M., de Oiveira-Costa, A., & Hamiton, A. 2003, Phys.Rev.D, 68, Tristram, M., Patanchon, G., Macías-Pérez, J. F., et a. 2005, A&A, 436, 785 Wandet, B. D., Larson, D. L., & Lakshminarayanan, A. 2004, Phys. Rev. D, 70, Wehus et a. 2013, in preparation Ysard, N., Mivie-Deschênes, M. A., & Verstraete, L. 2010, A&A, 509, L1 Ysard, N. & Verstraete, L. 2010, A&A, 509, A12 Appendix A: Physica parametrization The Commander-Ruer (C-R) approach impements Bayesian component separation in pixe space, fitting a parametric mode to the data by samping the posterior distribution for the mode parameters. For computationa reasons, the fit is performed in a two-step procedure: First, both foreground ampitudes and spectra parameters are found at ow-resoution using Markov Chain Monte Caro (MCMC)/Gibbs samping agorithms (Jewe et a. 2004; Wandet et a. 2004; Eriksen et a. 2004, 2007, 2008). Second, the ampitudes are recacuated at high resoution by soving the generaized east squares system (GLSS) per pixe with the spectra parameters fixed to the their vaues from the ow-resoution run. For the CMB-oriented anaysis presented in this paper, we ony use the seven owest Panck frequencies, i.e., from 30 to 353 GHz. We first downgrade each frequency map from its native anguar resoution to a common resoution of 40 arcminutes and repixeize at HEALPix N side = 256. Second, we set the monopoes and dipoes for each frequency band as described in Sect. 8.1 using a method that ocay records spectra indices (Wehus et a. 2013). We approximate the effective instrumenta noise as white with a root mean square (RMS) per pixe given by the Panck scanning pattern and an ampitude caibrated by smoothing simuations of the instrumenta noise incuding correations to the same resoution. For the high-resoution anaysis, the important pre-processing step is the upgrading of the effective ow-resoution mixing matrices to fu Panck resoution: this is done by repixeizing from N side = 256 to 2048 in harmonic space, ensuring that potentia pixeization effects from the ow-resoution map do not introduce sharp boundaries in the high-resoution map. Our mode for the data, a map d ν at frequency ν, consists of a inear combination of N c astrophysica components, pus instrumenta noise, N c d ν = F i ν(θ) A i + n ν, i=1 (A.1) where A i denotes a sky map vector containing the foreground ampitude map for component i normaized at a reference frequency, and F i ν(θ i ) is a diagona matrix describing the spectra emission aw for component i as a function of frequency and which depends on a (sma) set of spectra parameters, θ. The CMB signa is incuded in the sum and, as a specia case, it may be represented either in harmonic or pixe space, depending on whether the main goa of the anaysis is CMB power spectrum anaysis or component separation. The former representation is used for the Commander-based ow- ikeihood presented in Panck Coaboration XV (2013), whie the atter is used for the foreground fits presented in Section 8 of this paper. Bayes theorem specifies the posterior distribution for the mode parameters, P(A i, θ d) L(A i, θ)p(a i, θ), (A.2) where L(A i, θ) is a Gaussian ikeihood, and the prior P(A i, θ i ) depends on the appication. In this paper, the prior on spectra indices is a product of a Jeffreys prior and physica priors, as detaied in Section 8.2; no priors are imposed on the foreground ampitudes. In the ow-resoution Commander anaysis, we expoit a Gibbs samper to map out the posterior distribution (Eriksen et a. 2008), adopting the foowing minima two-step scheme: A P(A θ, d) θ P(θ A, d). (A.3) (A.4) The first conditiona distribution is a mutivariate Gaussian, whie the second distribution does not have an anaytic form and must be mapped out numericay. 21

22 The high-resoution Ruer anaysis maximizes the foreground ampitude conditiona in Eq. A.3 numericay by soving the generaized east squares system Fν i,t N 1 ν F i ν ) 1 F i,t µ N 1 µ D µ W(θ) i µd µ, A i (θ) = ( ν µ µ (A.5) were N ν is the noise covariance matrix (assumed to be diagona) of the νth channe, F i F i (θ) and we have negected the different anguar resoutions of the channes. The posterior margina average for the high-resoution ampitude maps is then given by A i = A i 1 (θ)p(θ)dθ N sampe θ,ν W(θ) i νd ν ν W i νd, a sum over the N sampe sampes of the spectra parameters θ. Once the channe weights, W ν, have been computed, processing a arge number of simuations requires negigibe computationa resources. This feature has been extensivey used for computation of the effective beam of Ruer maps: FFP6 CMB simuations for the 30 to 353 GHz channes are combined according to W ν and the effective beam transfer function is found as b 2 Cout w 2 ). Here, Cinp is the power spectrum of the input simuation before convoution with the instrumenta beam, w is the HEALPix pixe window function, and the average is taken over the set of simuations. Missing pixes are set to 0 when computing C in the above expression. The ow number ( 500) of missing pixes in the data renders the impact of such a choice negigibe at < A simiar procedure is used for defining the effective beam of the non-cmb components. The above agorithm produces a set of sampes drawn from the posterior distribution, as opposed to a direct estimate of individua component ampitudes or spectra parameters. Whie this sampe set provides a statisticay compete representation of the posterior, it is non-trivia to visuaize or to compare the distribution with externa data. For convenience, we therefore summarize the distribution in terms of mean and standard deviation maps for each component. We emphasize, however, that the distribution is significanty non-gaussian, and when searching for features in the maps at ow signa-to-noise eves, one must take into account the exact distribution. Finay, the Commander-Ruer confidence mask (see Sect. 5) is primariy defined by the product of the CG80 mask and the point source mask described in Sect. 4. We additionay remove any pixes excuded by the 13% Commander ikeihood mask described by Panck Coaboration XV (2013); however, this is amost entirey incuded within the CG80 mask, and this step therefore has very itte impact on the fina resut. To compete the mask, we remove any pixes for which the highresoution Ruer CMB map, smoothed to 40 arcminutes, differs by more than 3 σ from the ow-resoution Commander CMB map, which can happen due to spatia spectra variations on pixe scaes. /(C inp Appendix B: Interna inear combination is a method to extract the CMB (or any component with known spectra behaviour) by appying the ILC technique to muti-channe observations in needet space, that is, with weights that are aowed to vary over the sky and over the fu mutipoe range. The abiity to ineary combine input maps varying over the sky and over mutipoes is caed ocaization. In the needet framework, harmonic ocaization is achieved using a set of bandpass fiters defining a series of scaes, and spatia ocaization is achieved at each scae by defining zones over the sky. The harmonic ocaization adopted here uses nine spectra bands covering mutipoes up to = 3200 (see Fig. B.1). The spatia ocaization depends on the scae. At the coarsest scae, which incudes the mutipoes of owest degree, we use a singe zone (no ocaization), whie at the finest scaes (which incude the highest mutipoes) the sky is partitioned into 20 zones (again, see Fig. B.1). The method amounts to computing an ILC in each zone of each scae, aowing the ILC weights to adapt naturay to the varying strength of the other components as a function of position and mutipoe. A compete description of the basic method can be found in Deabrouie et a. (2009). In the present work, however, there is an important difference in the processing of the coarsest scae. Since the coarsest scae of the fiter is not ocaized, a pain map woud be equivaent to a pixe-based ILC for a the mutipoes of that scae. This procedure, however, is known to be quite susceptibe to the ILC bias due to chance correations between the CMB and foregrounds. In order to mitigate this effect, the (singe) covariance matrix which determines the ILC coefficients at the coarsest scae is not computed as a pixe average, but is rather estimated in the harmonic domain as an average over spherica harmonic coefficients using a spectra weight which equaizes the power of the CMB modes (based on a fiducia spectrum). This can be shown significanty to decrease the arge scae errors. In practice, our processing depends on severa impementation choices, as foows: Input channes: In this work, the agorithm is appied to a Panck channes from 44 to 857 GHz omitting ony the 30 GHz channe. Pre-processing of point sources: Identica to the preprocessing (see Appendix D). Masking and inpainting: The CMB map is actuay produced in a three-step process. In a first step, the weights are computed from covariance matrices evauated using a Gaactic mask removing about 2 % of the sky (and is apodized at 1 ). In a second step, those weights are appied to needet coefficients computed over the compete sky (except for point source masking/subtraction), yieding a CMB estimate over the fu sky (except for the point source mask). In short, the weights are computed over a masked sky but are appied to a fu sky (up to point sources). In a fina step, the pixes masked due to point source processing are repaced by the vaues of a constrained Gaussian reaization (inpainting). Spatia ocaization: The boundaries of the zones used for spatia ocaisation (shown at Fig. B.1) are obtained as isoeve curves of a ow resoution map of Gaactic emission. Beam contro and transfer function: As in the processing, the input maps are internay re-beamed to a 5 resoution, so the resuting CMB map is automaticay synthesized with an effective Gaussian beam of five arcminutes, according to the unbiased nature of the ILC. Using recaibration: In our current impementation, the soution uses the vaues determined by for the CMB spectrum, given in Eq. (D.6). Appendix C: Tempate fitting The origina agorithm produced cean CMB maps at severa frequencies through tempate fitting, foowed by an estimation of the CMB power spectrum from these cean maps us- 22

23 Tabe C.1: Linear coefficients, α j, and tempates used to cean individua frequency maps with. The first coumn ists the tempates constructed to produce cean maps. Before subtraction, the maps are smoothed to a common resoution. The ( ) GHz tempate is constructed at the resoution of the 44 and 70 GHz channes, in order to cean the 44 and 70 GHz maps, respectivey. For the rest of the tempates, the first map used to construct the tempates was fitered with the beam corresponding to the second map and vice versa. Note that for 100, 143 and 217 GHz channes, we give the coefficients used for the argest of the two regions considered in the ceaning, which covers 97 % of the sky. Tempate 44 GHz 70 GHz 100 GHz 143 GHz 217 GHz 353 GHz ing a method based on the Expectation Maximization agorithm (Martınez-Gonzaez et a. 2003; Leach et a. 2008; Fernández- Cobos et a. 2012). From this power spectrum, a mutifrequency Wiener-fitered CMB map was produced. For the present work, ony the first step of the method, producing cean CMB maps at different frequencies, is considered. In addition, two of these cean maps are optimay combined to produce a fina CMB map. The tempates used for ceaning are interna, i.e., they are constructed from Panck data, avoiding the need for externa data sets, which usuay compicate the anayses and may introduce inconsistencies. In the ceaning process, no assumptions about the foregrounds or noise eves are needed, rendering the technique very robust. The fitting can be done in rea or waveet space (using a fast waveet adapted to the HEALPix pixeization; Casaponsa et a. 2011) to propery dea with incompete sky coverage. For expediency, however, we fi in the sma number of unobserved pixes at each channe with the mean vaue of their neighbouring pixes before appying. We construct our tempates by subtracting two cose Panck frequency channe maps, after first smoothing them to a common resoution, and converting to CMB temperature units if necessary, to ensure that the CMB signa is propery removed. A inear combination of the tempates, t j, is then subtracted from (hitherto unused) map, d, to produce a cean CMB map at that frequency. This is done either in rea or waveet space (i.e., scae by scae) at each position on the sky, n t T c (x, ν) = d(x, ν) α j t j (x), j=1 (C.1) where n t is the number of tempates. If the ceaning is performed in rea space, the α j coefficients are obtained by minimizing the variance of the cean map, T c, outside a given mask. When working in waveet space, the ceaning is performed in the same way at each waveet scae independenty (i.e., the inear coefficients depend on the scae). Athough we excude very contaminated regions during the minimization, the subtraction is performed for a pixes and, therefore, the ceaned maps cover the fu-sky (athough foreground residuas are expected to be present in the excuded areas). An additiona eve of fexibiity may aso be considered: the inear coefficients can be fixed over the fu sky, or in severa regions. The regions are then combined in a smooth way, by weighting the pixes at the boundaries to reduce possibe discontinuities in the cean maps. Since the method is inear, we may easiy propagate the noise properties to the fina CMB map. Moreover, it is very fast and permits the generation of thousands of simuations to characterize the statistica properties of the outputs, a critica need for many cosmoogica appications. The fina CMB map retains the anguar resoution of the origina frequency map. There are severa configurations possibe for, depending the number of frequency maps to be ceaned or the number of tempates used in the fitting. Note that the production of cean maps at different frequencies is of great interest in order to test the robustness of the resuts. Therefore, to define the best strategy, one needs to find a compromise between the number of maps that can be ceaned independenty and the number of tempates that can be constructed. In particuar, we have ceaned the 143 GHz and 217 GHz maps using four tempates constructed as the difference of the foowing Panck channes (smoothed to a common resoution): (30 44), (44 70), ( ) and ( ). For simpicity, the two maps have been ceaned in rea space, since there was no significant improvement when using waveets, especiay at high atitude. In order to take into account the different spectra behaviour of the foregrounds at ow and high Gaactic atitudes, we considered two independent sky regions, using different sets of coefficients (see Tabe C.1 for the vaues of the inear coefficients for the main considered region). The first region corresponds to the brightest 3 % of Gaactic emission, whie the second region is defined by the remaining 97 % of the sky (see Sect. 4 for a detaied description of our Gaactic mask construction). For the first region, the coefficients are actuay estimated over the compete sky (we find that this is better than performing the minimization on ony the brightest 3 % of the sky, where the CMB is very sub-dominant), whie for the second region, we excude the bright 3 % sky fraction, point sources detected at any frequency and those pixes which have not been observed in a channes. Note that, for consistency, we have used the same configuration (four tempates, ceaning in rea space, two regions) for the anaysis of the FFP6 simuations. However, we find that this simpe configuration produces more contaminated CMB maps than for the data (athough the region outside the confidence mask sti has ow contamination), indicating some differences between the foreground eve in the data and in simuations. Therefore, concusions derived from the FFP6 resuts for shoud be taken with caution, since they are expected to provide overestimated residuas. 23

24 Window function Mutipoe moment, Our fina CMB map was constructed by combining the 143 and 217 GHz maps by weighting the maps in harmonic space, taking into account the noise eve, the resoution, and a rough estimation of the foreground residuas of each map (obtained from reaistic simuations). This fina map has a resoution corresponding to a Gaussian beam of 5 FWHM. Moreover, additiona cean CMB maps (at frequencies 44, 70, 100 and 353 GHz) were aso produced using different combinations of tempates. In particuar, to cean the 100 GHz map, we used the same tempates and regions as for 143 and 217 GHz. This aows us to produce three (amost) independent cean maps, in the sense that none of the three maps to be ceaned is used to construct the tempates. For 44, 70 and 353 GHz, different combinations of tempates are used and the inear coefficients are chosen to be the same over the fu sky. They are obtained by minimizing the variance of the map outside the same mask as that used to cean the centra frequency maps on the argest region. The tempates and the corresponding inear coefficients used for each of the considered frequencies are given in Tabe C.1. The cean frequency maps have been used in anayses of the isotropy and statistics of the CMB (Panck Coaboration XXIII 2013) and to obtain cosmoogica constraints from the integrated Sachs-Wofe effect (Panck Coaboration XIX 2013). In particuar, cean maps from 44 to 353 GHz were used for the stacking anaysis presented in Panck Coaboration XIX (2013), whie frequencies from 70 to 217 GHz were used for consistency tests in Panck Coaboration XXIII (2013). The confidence mask provided for is constructed by combining four types of seected regions. In particuar, it excudes zones with high residuas identified through the subtraction of cean maps at different frequencies, as we as the sources detected at a Panck channes using the Mexican Hat Waveet agorithm (Panck Coaboration XXVIII 2013). These point sources are masked with hoes of varying radius, according to the fux of each source. The pixes that are not observed by a channes are aso masked. Finay, to ensure that the area eft outside the mask is statisticay robust, we aso excude from the anaysis the brightest 20 % of Gaactic emission, eaving a usefu area of around 76 %. This provides a conservative mask for CMB anaysis; however, we point out that smaer masks coud aso be used in specific appications, such as the ensing potentia reconstruction described in Sect. 7). Appendix D: Spectra matching (Spectra Matching Independent Component Anaysis) reconstructs a CMB map as a inear combination in harmonic space of N chan input frequency maps with weights that depend on mutipoe. Given the N chan 1 vector x m of spherica harmonic coefficients for the input maps, it computes coefficients ŝ m for the CMB map as Fig. B.1: Spectra ocaization for using nine spectra window functions defining nine needet scaes (top pane). The scae-dependent spatia ocaization partitions the sky in one zone (for scae 1), two zones (for scae 2), four zones (for scae 3), or twenty zones (for scaes 5, 6, 7, 8, 9). The two-zone, 4- zone and 20-zone partitions are depicted in the ower panes. ŝ m = w x m, (D.1) where the N chan 1 vector w containing the mutipoedependent weights is chosen to give unit gain to the CMB with minimum variance. This is achieved with w = C 1 a a C 1 a, (D.2) where vector a is the spectrum of the CMB evauated at each channe (aowing for possibe inter-channe re-caibration factors) and C is the N chan N chan spectra covariance matrix of 24

25 ] [ w () i µk/µk RJ i i w i () 2 N i () [ µk 2 ] ] w i () [ µk/µk RJ Fig. D.1: Weights, w, given by to the input maps, after they are re-beamed to 5 and expressed in K RJ, as a function of mutipoe. Top pane: inear scae; bottom pane: the absoute vaue of the weights on a ogarithmic scae. x m. Taking C in Eq. (D.2) to be the sampe spectra covariance matrix Ĉ of the observations, Ĉ = x m x m, m (D.3) woud impement a simpe harmonic-domain ILC simiar to (Tegmark et a. 2003). At the argest scaes, we instead use a mode, C (θ), and determine the covariance matrix to be used in Eq. (D.2) by fitting C (θ) to Ĉ. This is done in the maximum ikeihood sense for stationary Gaussian fieds, that is, the best fit matrices, C (ˆθ), are obtained for ˆθ = arg min θ (2 + 1) ( Ĉ C (θ) 1 + og det C (θ) ). (D.4) Equations D.1-D.4 summarize the basic principes of ; its actua operation depends on a choice for the spectra mode C (θ), and on severa impementation-specific detais, which we briefy describe beow. Fig. D.2: Contribution of each input channe to the noise in the map. We mode the data as a superposition of CMB, noise and foregrounds. The atter are not parametricay modeed; instead, we represent the tota foreground emission by d tempates with arbitrary frequency spectra, anguar spectra and correations. In the spectra domain, this is equivaent to modeing the covariance matrices as C (θ) = aa C + AP A + N, (D.5) where C is the anguar power spectrum of the CMB, A is a N chan d matrix, P is a positive definite d d matrix, and N is a diagona matrix representing the noise power spectra of the data. The parameter vector θ contains a or part of the quantities in (D.5). The decomposition D.5 refects the fact that CMB, foregrounds and noise are independent components of the signa. Thus, is an ICA (independent component anaysis) method. It operates by matching the observations Ĉ to the spectra mode (D.5) using the criterion (D.4). The maxima fexibiity in a fit of mode (D.5) is obtained with a the parameters free, that is without any constraint on the spectrum C, on the diagona entries of N, on a, or on A and P. One woud ideay fit a those parameters (except for obvious degeneracies, ike that between a scae factor in a and the overa normaization of the CMB spectrum C ) over the whoe mutipoe range. In practice, this turns out to be too difficut given the arge dynamic range both over the sky and over mutipoes. We resort to a pragmatic three-step approach in which the criterion (D.4) is minimized by first fitting a, then A, and finay the inear parameters C and N. Each fit is conducted over the mutipoe ranges and the sky fraction most appropriate for the parameter of interest, as foows. We first estimate the CMB spectra aw a by fitting a mode parameters (that is, without constraint) over a cean fraction of sky ( f sky =40 %) in the range where the signa is CMB-dominated in most of the channes and the beam window functions are accuratey known. In this fit, which is done over a cean part of the sky, we use a foreground emission matrix, A, with ony four coumns. From this step, we ony retain the best fit vaue for vector a. In the second step, we estimate the foreground emissivity by fixing a to its vaue from the previous 25

26 step and fitting a the other parameters over a arge fraction of sky ( f sky =97 %) in the range where the signa is dominated by the Gaactic emission in a channes. From this second step, we retain the best fit vaue for the matrix A which, again, is adjusted without constraint other than having d = 6 coumns. In the ast step, we fit a power spectrum parameters: we fix a and A to their previousy found vaues and fit for C and P at each. Note that the first step (fitting a) amounts to re-caibrating the input maps on the basis of CMB anisotropies. For the maps in thermodynamics units, we find â = [0.9900, , , , , , , , ] (D.6) The vaue at 857 GHz is not accuratey recovered by, so we have set a 857 = 1. Since the norm of a is degenerate with a goba scae factor for the CMB anguar spectrum, it can ony be recovered by up to a scae factor. This degeneracy is fixed here by taking a 143 = 1. The re-caibration step coud have been omitted since â is very cose the unit vector. However, we found that using â improved the behavior of over using a = [1,..., 1]. Before describing impementation detais, we expain how deas with the varying resoution of the input channes, since the discussion thus far assumed that a input maps had the same resoution. Since works in the harmonic domain, it is a simpe matter to account for the beam transfer function, b i (), of the i-th input map. The CMB sky mutipoes s m contribute s m a i b i ()p i () to the harmonic coefficient x i m of the i-th map (where p i () is the pixe window function for the HEALPix map at N i side ). Therefore, in order to produce a fina CMB map at 5 resoution, cose to the highest resoution of Panck, we ony need to work with input spherica harmonics re-beamed to 5 ; that is, to appy on vectors x m with entries x i m = xi m b 5()/b i ()/p i (), where b 5 () is a five-arcminute Gaussian beam function. By construction, then produces an CMB map with an effective Gaussian beam of 5 (without the pixe window function). We now give further detais on the actua impementation of : Inputs: uses a nine Panck frequency channes from 30 to 857 GHz, harmonicay transformed up to = Pre-processing of point sources: is appied on input maps in which point sources are subtracted or masked. We start by fitting the PCCS point sources with SNR > 5 to a Gaussian shape where the source ampitude is estimated together with its position and a constant factor representing the background variance. If the fit is successfu (χ 2 2), the fitted point source is removed from the map; otherwise it is masked in a channes and the hoe is inpainted by a simpe diffusive fiing process. This is done at a frequencies except 545 and 857 GHz, where a point sources with SNR > 7.5 are masked and inpainted. Beams: When the harmonic coefficients of the input maps are re-beamed at 5, we do not appy exacty the expression x i m = xi m b 5()/b i ()/p i () mentioned above because the factor 1/b i () woud diverge at high for the owest resoution input channes. That may not be a probem in infinite precision arithmetic, but woud ead to matrices Ĉ() with extremey arge condition numbers. Instead, we re-beam with the factor 1/b i () repaced by min(1/b i (), 1000). The rebeaming of the CMB modes then is no onger perfect, but this is of course irreevant because the threshoding occurs in a regime where the signa is competey dominated by the noise, so that the contribution of the corresponding channe is aready highy attenuated by the weights (as shown in Fig. D.1). Masking: In practice, operates on a masked sky, the mask being appied after the point source processing. The mask is obtained by threshoding a heaviy smoothed version of the point source mask. The threshod is chosen to eave about 97 % of the sky. Because of the heavy smoothing, the mask has smooth contours and is ony sensitive to arge aggregates of point sources: the masked areas mosty ie in the Gaactic pane, but incude aso a few bright regions ike the Large Mageanic Coud. Inpainting: The map used in this paper has no rea power in the masked region described above. However, for convenience, an inpainted map has aso been produced by repacing the masked pixes with a constrained Gaussian reaization obtained by the method of Benoit-Lévy et a. (2013). That map appears in Panck Coaboration I (2013). Binning: In our impementation, we use binned spectra. Processing at fine scaes: Since there is itte point trying to mode the spectra covariance at high mutipoes, because the sampe estimate is sufficient, impements a simpe harmonic ILC at > 1500; that is, it appies the fiter (D.2) with C = Ĉ. Confidence mask: A confidence mask (Fig. 3) is provided with, constructed in the foowing way. The CMB map is bandpass fitered through a spectra window v() = exp[ (( 1700)/200) 2 /2]. The resut is squared and smoothed at two-degree resoution, yieding a map of the (bandpassed) variance of the CMB map. That variance is corrected for the noise contribution by subtracting the variance map for the noise obtained by the same procedure appied to the HRHD map. If the map contained ony CMB and noise, the variance map woud have a uniform vaue v() 2 b 5 () 2 C()(2 + 1)/4π = 31.3 µk 2 over the sky. The confidence map is obtained by threshoding the noise-corrected variance map at 70 µk 2. Viewed as a fiter, can be summarized by the weights w appied to each input map as a function of mutipoe. In this sense, is stricty equivaent to co-adding the input maps after convoution by specific axisymmetric kernes directy reated to the corresponding entry of w. The weights used here are shown in Figure D.1 (for input maps in units of K RJ ). We see, in particuar, the (expected) progressive attenuation of the owest resoution channes with increasing mutipoe. Figure D.2 shows the contribution of each input channe to the noise in the CMB map as a function of mutipoe. The spectra noise contribution from channe i is simpy obtained as w() 2 N i (), where w i () is the i-th entry of the weight vector w() and N i () is the anguar spectrum of the i-th noise map. More detais about are given in Cardoso et a. (2008), as we as in appications to the anaysis of WMAP (Patanchon et a. 2005) and Archeops data (Tristram et a. 2005). An appication to the measurement of the tensor to scaar ratio using CMB B-modes is discussed in Betoue et a. (2009). Within the Panck coaboration, is used to define the Pik high- ikeihood (Panck Coaboration XV 2013), but physica modes of foreground emission are used there instead of the nonparametric foreground mode used here. is aso used to 26

27 0 µk 10 Fig. E.1: Standard deviation of the four foreground-ceaned CMB maps from the FFP6 simuation. A maps have been downgraded to a HEALPix resoution of N side = 128. cross-check the HFI caibration (see Panck Coaboration VIII (2013)). Appendix E: FFP6 simuation resuts We study the performance of the component separation agorithms on the FFP6 simuation described in Sect. 4, providing additiona information beyond that in the body of the paper. Much of the anaysis presented here mirrors that shown for the data in Sects. 5 and 6, aowing a direct comparison between the two anayses. E.1. CMB maps First, we show in Fig. E.1 the standard deviation between the four foreground-ceaned FFP6 maps, simiar to that shown in Fig. 5 for the data. Figure E.2 shows a pairwise differences between the same maps, mirroring Fig. 6 for the data. These two pots highight an important point concerning the FFP6 anaysis aready mentioned in Sect. 5, namey that in near-gaactic regions, where the foregrounds are important, the interna differences between the four agorithms are arger in the FFP6 simuation than in the rea data. This is due to the fact that each component separation agorithm has been optimized with the rea data in terms of mode definition, ocaization, etc. Then, the same modes have been used for the FFP6 simuation without change. Ony the parameters within those modes are refitted to the new data set. This impies in fact that we expect each method to perform better on the data than the simuations in terms of absoute residuas, to the extent that the simuation matches the rea sky. In other words, the FFP6 simuation provides a conservative estimate of the residua errors in the rea data. E.2. Power spectrum residuas from individua components In Fig. E.3 we show the residua effect of some of the individua components on the foreground-ceaned CMB map. The therma dust emission, CIB fuctuations, point sources, and noise have been processed individuay with each agorithm. A other components (free-free, synchrotron, spinning dust, CO, therma SZ, and kinetic SZ) are shown as a singe, composite residua component. E.3. CMB power spectra and cosmoogica parameters We assess the performance of our component separation techniques by evauating cosmoogica constraints from the foreground-ceaned CMB maps derived from the FFP6 simuation. Figure E.4 shows the estimates of the anguar power spectra of the CMB maps. Figure E.5 compares the cosmoogica parameters derived from the four foreground-ceaned CMB maps, CamSpec 6, and Pik to the input (theoretica) parameters for different -ranges. The parameter space is defined by the same mode appied to the rea data in Sect. 6, incuding six ΛCDM and two foreground parameters. A deviations from input parameters are sma and within 1 σ up to = 2000, verifying that a methods work we in this mutipoe range. However, for max = 2500 we start to see significant shifts, e.g., for Ω b h 2 and n s. Further, the point source foreground parameter, A ps, reaches arge vaues, impying that assumptions concerning the high- foreground mode become important. For these reasons, we consider max = 2000 as the maximum recommended -range for these maps in the current data reease. Sti, the overa agreement is exceent between a codes and a -ranges. In particuar, we see that differences in the band power spectra at high between the different codes are mosty absorbed by the two-parameter foreground mode. For instance, the Commander-Ruer band power spectrum has more power at high due to noise or residua point sources, but this excess is we fitted by the two-parameter foreground mode, and mosty interpreted in terms of a residua point source component; this is expected, given the ower anguar resoution of this map. As mentioned above, n s and and A are to some extent sensitive to max. These parameters are degenerate with the foreground parameters. This may suggest that our C foreground tempates deviate more from the shape of the Poissonian and custered component in the CMB map. This is a imitation of the simpe foreground tempates used here. To propery describe the foreground residuas in the reconstructed maps, we shoud use a foreground power spectrum tempate taiored to each method. For instance, such tempates may be constructed by processing simuated foreground maps though each of the four pipeines. The tempates are then given by the pseudo-c of each of the processed foreground map. However, our anaysis shows that the current simpe mode provides accurate resuts when restricting the anaysis to max = Figure E.6 shows the best-fit power spectrum residuas for the CMB map, CamSpec and Pik reative to the input CMB ΛCDM mode estimated up to = These pots show that the residuas of the CMB map-based best-fit modes are comparabe to the CamSpec and Pik residuas, and smaer than 40 µk 2 for most of the range with arger deviations observed for CamSpec at 200. At higher s the residuas are smaer than 10 µk 2 for both approaches, a showing simiar trends. Thus, both the map- and spectrum-based ikeihoods recover input parameters reasonaby we, with the atter yieding sighty arger deviations from the best-fit mode of the input CMB reaization. 1 APC, AstroParticue et Cosmoogie, Université Paris Diderot, CNRS/IN2P3, CEA/rfu, Observatoire de Paris, Sorbonne Paris Cité, 10, rue Aice Domon et Léonie Duquet, Paris Cedex 13, France 6 For CamSpec, k pivot = 0.05 was adopted for this test, whie a others, input parameters and input CMB reaization incuded, use k pivot =

28 C-R - C-R - C-R µk 30 Fig. E.2: Pairwise differences between foreground-ceaned CMB maps from the FFP6 simuation. A maps have been downgraded to a HEALPix resoution of N side = 128 to show the arge-scae differences. 2 Aato University Metsähovi Radio Observatory, Metsähovintie 114, FIN Kymää, Finand 3 African Institute for Mathematica Sciences, 6-8 Merose Road, Muizenberg, Cape Town, South Africa 4 Agenzia Spaziae Itaiana Science Data Center, c/o ESRIN, via Gaieo Gaiei, Frascati, Itay 5 Agenzia Spaziae Itaiana, Viae Liegi 26, Roma, Itay 6 Astrophysics Group, Cavendish Laboratory, University of Cambridge, J J Thomson Avenue, Cambridge CB3 0HE, U.K. 7 Astrophysics & Cosmoogy Research Unit, Schoo of Mathematics, Statistics & Computer Science, University of KwaZuu-Nata, Westvie Campus, Private Bag X54001, Durban 4000, South Africa 8 Atacama Large Miimeter/submiimeter Array, ALMA Santiago Centra Offices, Aonso de Cordova 3107, Vitacura, Casia , Santiago, Chie 9 CITA, University of Toronto, 60 St. George St., Toronto, ON M5S 3H8, Canada 10 CNR - ISTI, Area dea Ricerca, via G. Moruzzi 1, Pisa, Itay 11 CNRS, IRAP, 9 Av. coone Roche, BP 44346, F Tououse cedex 4, France 12 Caifornia Institute of Technoogy, Pasadena, Caifornia, U.S.A. 28