A&A 536, A20 (2011) Astronomy DOI: 10.1051/0004-6361/201116470 & c ESO 2011 Astrophysics Planck early results Special feature

Planck early results. XX. New light on anomalous microwave emission from spinning dust grains

Planck Collaboration: P. A. R. Ade72, N. Aghanim46,M.Arnaud58, M. Ashdown56,4, J. Aumont46, C. Baccigalupi70,A.Balbi28, A. J. Banday77,7,63,R.B.Barreiro52, J. G. Bartlett3,54,E.Battaner79, K. Benabed47, A. Benoît45,J.-P.Bernard77,7, M. Bersanelli25,41, R. Bhatia5, J. J. Bock54,8, A. Bonaldi37,J.R.Bond6,J.Borrill62,73,F.R.Bouchet47, F. Boulanger46, M. Bucher3,C.Burigana40,P.Cabella28, B. Cappellini41, J.-F. Cardoso59,3,47,S.Casassus76, A. Catalano3,57, L. Cayón18, A. Challinor49,56,10, A. Chamballu43, R.-R. Chary44,X.Chen44,L.-Y.Chiang48, C. Chiang17,P.R.Christensen67,29,D.L.Clements43, S. Colombi47, F. Couchot61, A. Coulais57, B. P. Crill54,68, F. Cuttaia40,L.Danese70, R. D. Davies55,R.J.Davis55,P.deBernardis24,G.deGasperis28,A.deRosa40, G. de Zotti37,70, J. Delabrouille3, J.-M. Delouis47, C. Dickinson55, S. Donzelli41,50,O.Doré54,8,U.Dörl63, M. Douspis46, X. Dupac32, G. Efstathiou49,T.A.Enßlin63,H.K.Eriksen50, F. Finelli40, O. Forni77,7, M. Frailis39, E. Franceschi40,S.Galeotta39, K. Ganga3,44,R.T.Génova-Santos51,30,M.Giard77,7, G. Giardino33, Y. Giraud-Héraud3, J. González-Nuevo70,K.M.Górski54,81,S.Gratton56,49, A. Gregorio26, A. Gruppuso40,F.K.Hansen50,D.Harrison49,56,G.Helou8, S. Henrot-Versillé61, D. Herranz52,S.R.Hildebrandt8,60,51,E.Hivon47, M. Hobson4,W.A.Holmes54,W.Hovest63,R.J.Hoyland51, K. M. Huffenberger80,T.R.Jaffe77,7,A.H.Jaffe43,W.C.Jones17,M.Juvela16,E.Keihänen16, R. Keskitalo54,16,T.S.Kisner62, R. Kneissl31,5, L. Knox20, H. Kurki-Suonio16,35, G. Lagache46, A. Lähteenmäki1,35, J.-M. Lamarre57,A.Lasenby4,56,R.J.Laureijs33,C.R.Lawrence54, S. Leach70,R.Leonardi32,33,21, P. B. Lilje50,9, M. Linden-Vørnle12,M.López-Caniego52,P.M.Lubin21, J. F. Macías-Pérez60,C.J.MacTavish56, B. Maffei55,D.Maino25,41, N. Mandolesi40,R.Mann71,M.Maris39,D.J.Marshall77,7, E. Martínez-González52,S.Masi24,S.Matarrese23, F. Matthai63, P. Mazzotta28, P. McGehee44, P. R. Meinhold21, A. Melchiorri24, L. Mendes32, A. Mennella25,39, S. Mitra54, M.-A. Miville-Deschênes46,6, A. Moneti47,L.Montier77,7, G. Morgante40,D.Mortlock43, D. Munshi72,49, A. Murphy66, P. Naselsky67,29, P. Natoli27,2,40,C.B.Netterfield14, H. U. Nørgaard-Nielsen12, F. Noviello46, D. Novikov43, I. Novikov67,I.J.O’Dwyer54, S. Osborne75,F.Pajot46, R. Paladini74,8, B. Partridge34,F.Pasian39, G. Patanchon3,T.J.Pearson8,44, M. Peel55, O. Perdereau61,L.Perotto60,F.Perrotta70, F. Piacentini24, M. Piat3, S. Plaszczynski61,P.Platania53, E. Pointecouteau77,7,G.Polenta2,38, N. Ponthieu46, T. Poutanen35,16,1, G. Prézeau8,54, P. Procopio40, S. Prunet47, J.-L. Puget46, W. T. Reach78,R.Rebolo51,30, W. Reich64, M. Reinecke63, C. Renault60, S. Ricciardi40, T. Riller63, I. Ristorcelli77,7, G. Rocha54,8, C. Rosset3, M. Rowan-Robinson43 , J. A. Rubiño-Martín51,30, B. Rusholme44, M. Sandri40,D.Santos60,G.Savini69,D.Scott15, M. D. Seiffert54,8, P. Shellard10,G.F.Smoot19,62,3, J.-L. Starck58,11,F.Stivoli42, V. Stolyarov4, R. Stompor3,R.Sudiwala72, J.-F. Sygnet47, J. A. Tauber33, L. Terenzi40,L.Toffolatti13,M.Tomasi25,41, J.-P. Torre46, M. Tristram61, J. Tuovinen65,G.Umana36, L. Valenziano40,J.Varis65, L. Verstraete46,P.Vielva52, F. Villa40, N. Vittorio28,L.A.Wade54,B.D.Wandelt47,22,R.Watson55, A. Wilkinson55,N.Ysard16,D.Yvon11, A. Zacchei39, and A. Zonca21 (Affiliations can be found after the references)

Received 8 January 2011 / Accepted 11 May 2011

ABSTRACT

Anomalous microwave emission (AME) has been observed by numerous experiments in the frequency range ∼10–60 GHz. Using Planck maps and multi-frequency ancillary data, we have constructed spectra for two known AME regions: the Perseus and ρ Ophiuchi molecular clouds. The spectra are well fitted by a combination of free-free radiation, cosmic microwave background, thermal dust, and electric dipole radiation from small spinning dust grains. The spinning dust spectra are the most precisely measured to date, and show the high frequency side clearly for the first time. The spectra have a peak in the range 20–40 GHz and are detected at high significances of 17.1σ for Perseus and 8.4σ for ρ Ophiuchi. In Perseus, spinning dust in the dense molecular gas can account for most of the AME; the low density atomic gas appears to play a minor role. In ρ Ophiuchi, the ∼30 GHz peak is dominated by dense molecular gas, but there is an indication of an extended tail at frequencies 50–100 GHz, which can be accounted for by irradiated low density atomic gas. The dust parameters are consistent with those derived from other measurements. We have also searched the Planck map at 28.5 GHz for candidate AME regions, by subtracting a simple model of the synchrotron, free-free, and thermal dust. We present spectra for two of the candidates; S140 and S235 are bright H ii regions that show evidence for AME, and are well fitted by spinning dust models. Key words. ISM: general – Galaxy: general – radiation mechanisms: general – radio continuum: ISM – submillimeter: ISM

1. Introduction AME has been observed by numerous experiments over the frequency range ∼10–60 GHz and is found to be very closely Anomalous microwave emission (AME) is an additional com- ff correlated with far infrared (FIR) emission associated with ther- ponent of di use foreground emission that cannot be easily mal emission from dust grains (Kogut et al. 1996; Leitch et al. explained by synchrotron, free-free, or thermal dust emission. 1997; de Oliveira-Costa et al. 1997; Banday et al. 2003; Lagache  Corresponding author: C. Dickinson, 2003; de Oliveira-Costa et al. 2004; Finkbeiner 2004; Finkbeiner e-mail: [email protected] et al. 2004; Davies et al. 2006; Dobler & Finkbeiner 2008;

Article published by EDP Sciences A20, page 1 of 17 A&A 536, A20 (2011)

Miville-Deschênes et al. 2008; Gold et al. 2011; Ysard et al. components. The AME is clearly detected as a convex spec- 2010). Electric dipole radiation from small rapidly spinning trum that is well-fitted by theoretical models of spinning dust dust grains, or “spinning dust”, is thought to be emitted in emission. We investigate the parameters for the spinning dust to the microwave region of the spectrum1. Theoretical models see whether they are reasonable, given the conditions in these en- predict a peaked spectrum, at a range of possible frequencies vironments. The Planck data, in combination with ancillary data, (∼10–150 GHz) depending on the properties of the dust grains are also used to search for new regions of AME. We present two and their environment (Draine & Lazarian 1998; Ali-Haïmoud such regions and show that these newly discovered regions can et al. 2009; Dobler et al. 2009; Ysard & Verstraete 2010). The also be fitted with a contribution from spinning dust grains. spectrum is very different from those of the traditional dif- In this paper, Sects. 2 and 3 describe the Planck and ancil- fuse emission expected at these frequency ranges; for example, lary data that were used to measure the spectra of AME regions. free-free, synchrotron, and thermal dust emission have power- Section 4 presents the results for the Perseus molecular cloud law-like spectra at microwave frequencies. For these reasons, and Sect. 5 for the ρ Ophiuchi cloud. Section 6 describes prelim- spinning dust emission has become the best explanation for the inary modelling of the AME in terms of spinning dust. Section 7 AME, although other physical mechanisms could still be con- presents the search for new AME regions and gives two clear tributing at some level, such as hot free-free emission (Leitch examples. Conclusions are given in Sect. 8. et al. 1997), hard synchrotron radiation (Bennett et al. 2003), or magneto-dipole emission (Draine & Lazarian 1999). The most direct evidence for spinning dust grains has come 2. Planck data from dedicated observations of known dust clouds in the inter- Planck (Tauber et al. 2010; Planck Collaboration 2011a)isthe stellar medium (Finkbeiner et al. 2002; Casassus et al. 2006; third-generation space mission to measure the anisotropy of the Dickinson et al. 2009; Scaife et al. 2009; Dickinson et al. 2009, cosmic microwave background (CMB). It observes the sky in 2010). Perhaps the best example is that of the Perseus molec- nine frequency bands covering 30–857 GHz with high sensitiv- ular cloud. Watson et al. (2005) observed the region with the ity and angular resolution from 31 to 5. The Low Frequency COSMOSOMAS telescope covering 11–17 GHz. By combining Instrument (LFI; Mandolesi et al. 2010; Bersanelli et al. 2010; data from radio surveys, WMAP, and DIRBE, they were able to Mennella et al. 2011) covers the 30, 44, and 70 GHz bands show very clearly a strong excess at frequencies ∼10–60 GHz with amplifiers cooled to 20 K. The High Frequency Instrument and a shape similar to spinning dust models. Another clear ex- (HFI; Lamarre et al. 2010; Planck HFI Core Team 2011a)covers ample is that of the photodissociation region (PDR) associated the 100, 143, 217, 353, 545, and 857 GHz bands with bolome- with the ρ Ophiuchi molecular dust cloud. Casassus et al. (2008) ◦ ◦ ters cooled to 0.1 K. Polarisation is measured in all but the observed the region centred at (l, b) = (353.0, +17.0) at 31 GHz highest two bands (Leahy et al. 2010; Rosset et al. 2010). A with the Cosmic Background Imager (CBI) and found a close combination of radiative cooling and three mechanical cool- connection between the centimetre-wavelength emission and the ers produces the temperatures needed for the detectors and op- dust emission. Further analysis of the spectrum confirmed that tics (Planck Collaboration 2011b). Two data processing centres there was significant excess emission that was consistent with (DPCs) check and calibrate the data and make maps of the sky spinning dust models. AME has also been detected within a star- (Planck HFI Core Team 2011b; Zacchei et al. 2011). Planck’s forming region in the galaxy NGC 6946 (Murphy et al. 2010), sensitivity, angular resolution, and frequency coverage make it a which appears to be well-fitted by spinning dust models (Scaife powerful instrument for Galactic and extragalactic astrophysics et al. 2010). as well as cosmology. Early astrophysics results are given in Planck2 data now represent the next step in measuring and Planck Collaboration (2011h–z). understanding spinning dust emission. The additional data at In this paper, we start with the Planck (non CMB-subtracted) 28.5, 44.1, and 70.3 GHz (where spinning dust is expected to maps at nominal frequencies of 28.5, 40.1, 70.3, 100, 143, 217, emit strongly) and the high frequency data at 100–857 GHz al- 353, 545, and 857 GHz (see Table 1). For the spectral analyses, low an accurate model for the thermal emission to be removed. we smoothed the HEALPix maps to a common angular reso- ◦ This is crucial in defining a precise spectrum for the spinning lution of 1 and degraded to Nside = 512 (Górski et al. 2005). dust, particularly at the higher frequencies where the smallest We assume Gaussian beams and use the average beamwidths dust grains dominate. Accurate spectra for spinning dust could given in Planck HFI Core Team (2011b); Zacchei et al. (2011). potentially prove to be an important diagnostic tool for measur- Details of the beam, such as the exact shape and variation of the ing the properties of the dust (e.g., the size distribution, density beam across the sky, do not significantly affect the results pre- and average electric dipole moment) as well as the physical envi- sented here, due to the applied smoothing which dominates the ronment in the surrounding region (e.g., the interstellar radiation spatial response of the maps. Furthermore, extended emission is field). less sensitive to the details of the beam. We convert from CMB In this paper, we study two known AME regions, the Perseus thermodynamic units to Rayleigh-Jeans brightness temperature and ρ Ophiuchi molecular cloud regions, and provide spectra units using the standard conversion factors given in Planck HFI over the frequency range 0.4–3000 GHz. We show, for the first Core Team (2011b)andZacchei et al. (2011). Colour corrections time, an accurate residual spectrum after removal of the free- to account for the finite bandpass at each frequency are applied free, cosmic microwave background (CMB), and thermal dust during the modelling of the spectra for each object; these are typ-

1 ically a few per cent for LFI and 10–15% for HFI. The 100 GHz Erickson (1957) was the first to suggest the possibility of radio emis- data are significantly contaminated by the CO J = 1 → 0 line sion from spinning dust grains. at 115 GHz. We apply a nominal correction to the 100 GHz map 2 Planck (http://www.esa.int/Planck) is a project of the using the CO survey of Dame et al. (2001). We multiply the in- European Space Agency (ESA) with instruments provided by two sci- ◦ −1 entific consortia funded by ESA member states (in particular the lead tegrated 1 -smoothed line intensity map (in units of K km s ) countries France and Italy), with contributions from NASA (USA) and by 14.2(Planck HFI Core Team 2011b) and subtract it from telescope reflectors provided by a collaboration between ESA and a sci- the 100 GHz map. Instead of attempting to correct for other entific consortium led and funded by Denmark. lines that may contaminate the 217 and 353 GHz data, we do

A20, page 2 of 17 Planck Collaboration: Planck early results. XX.

Table 1. Summary of the data.

Frequency Telescope/ Angular resolution Coverage Notes [GHz] survey [arcmin] 0.408 JB/Eff/Parkesa ≈60 Full-sky NCSA desourced/destriped versionb 0.820 Dwingelooc 72 Northern sky (−7◦ <δ<+85◦) Available on the webd 1.420 Stockert/Villa-Elisae 36 Full-sky Courtesy of Reich 2.326 HartRAO f 20 Southern sky (−83◦ <δ<+13/32◦) Courtesy of Jonas (priv. comm.) 11–17 COSMOSOMASg ≈60 Northern sky (+24◦ <δ<+44◦) Filtered on large angular scales 22.8 WMAP 7-yrh ≈49 Full-sky 1◦ smoothed versionb 28.5 Planck LFIi 32.65 Full-sky 33.0 WMAP 7-yrh ≈40 Full-sky 1◦ smoothed versionb 40.7 WMAP 7-yrh ≈31 Full-sky 1◦ smoothed versionb 44.1 Planck LFIi 27.92 Full-sky 60.7 WMAP 7-yrh ≈21 Full-sky 1◦ smoothed versionb 70.3 Planck LFIi 13.01 Full-sky 93.5 WMAP 7-yrh ≈13 Full-sky 1◦ smoothed versionb 100 Planck HFI j 9.37 Full-sky 143 Planck HFI j 7.04 Full-sky 217 Planck HFI j 4.68 Full-sky 353 Planck HFI j 4.43 Full-sky 545 Planck HFI j 3.80 Full-sky 857 Planck HFI j 3.67 Full-sky 1249 COBE-DIRBEk ≈40 Full-sky LAMBDA websiteb 2141 COBE-DIRBEk ≈40 Full-sky LAMBDA websiteb 2997 COBE-DIRBEk ≈40 Full-sky LAMBDA websiteb

Notes. (a) Haslam et al. (1982); (b) http://lambda.gsfc.nasa.gov/; (c) Berkhuijsen (1972); (d) http://www.mpifr-bonn.mpg.de/ survey.html; (e) Reich (1982); Reich & Reich (1986); Reich et al. (2001); ( f ) Jonas et al. (1998); (g) Watson et al. (2005); (h) Jarosik et al. (2011); (i) Zacchei et al. (2011); ( j) Planck HFI Core Team (2011b); (k) Hauser et al. (1998). not include these channels in the modelling of spectra if they pixels. This procedure has been shown to conserve photome- are observed to be in excess by more than 2σ, relative to the try (Paradis et al., in prep.). We note that there are significant best-fit model when omitting these frequencies. Uncertainties baseline uncertainties (i.e., offsets) in these maps, but these will are taken to be the nominal uncertainties in absolute calibration not affect our results since we are subtracting a local background and gain (given in the aforementioned papers), which are added for each object, which removes any such offset. The data include in quadrature. If these are below 3%, we conservatively assume the full-sky map at 408 MHz (Haslam et al. 1982), the 820 MHz an overall uncertainty of 3%. Dwingeloo survey (Berkhuijsen 1972), the new full-sky map at Jack-knife tests, in which the data were split into two halves 1420 MHz (Reich 1982; Reich & Reich 1986; Reich et al. 2001), (in terms of time and also HFI detectors; see Planck HFI and the 2326 MHz HartRAO southern survey (Jonas et al. 1998). Core Team 2011b), provided a check of the consistency of the The 1420 MHz map was multiplied by a factor of 1.55 to account Planck data. The tests confirmed that the results were consistent for the main beam to full beam ratio (Reich & Reich 1988), and with each other to within a small fraction of the derived uncer- to bring the calibration in line with the other data. We assume a tainties. 10% overall calibration uncertainty for these data.

3. Ancillary data 3.2. COSMOSOMAS data Although Planck observes over a wide frequency range (30–857 GHz), it is still important to include ancillary data to The COSMOSOMAS experiments (Gallegos et al. 2001) extend the frequency range into the FIR, and in particular, to consisted of two circular scanning instruments COSMO11 lower frequencies. To supplement Planck data we therefore use (Hildebrandt et al. 2007)andCOSMO15(Fernández-Cerezo a range of publicly available surveys and data available within et al. 2006) which operated at the Teide observatory (altitude the team. Table 1 summarises the data used in the analysis and 2400 m, Tenerife). They have produced 0◦.8–1◦.1 resolution maps the origin of the maps that were used. In the following sections, of ≈10 000 deg2 in four frequency bands (10.9, 12.7, 14.7, we describe the various data sets in more detail. 16.3 GHz) that are ideal for filling the gap between low fre- quency radio surveys and the lowest frequency channels of 3.1. Low frequency radio data WMAP and Planck. The average sensitivity of each map is in the range 80–120 μK per beam, and the observing strategy is The radio data consist of a number of low frequency surveys in described in Gallegos et al. (2001). As COSMOSOMAS was the range 0.408–2.326 GHz (Table 1). We used the most up-to- a ground-based experiment, there is significant contamination date versions of the data sets, and where necessary, regridded from atmospheric and receiver 1/ f noise, which needs to be fil- the maps into the HEALPix format (Górski et al. 2005)byusing tered out. This results in a loss of response to angular scales a procedure which computes the surface intersection between larger than 5◦ and introduces low level correlations over the individual pixels of the survey with the intersecting HEALPix scan ring. This must be taken into account when considering

A20, page 3 of 17 A&A 536, A20 (2011)

28.5 GHz44.1 GHz 70.3 GHz 100 GHz

4C+32.14

AME−G160.26−18.62

143 GHz 857 GHz 1.4 GHz H−alpha

IC348

NGC1333

Fig. 1. Maps of the Perseus molecular cloud region at their original angular resolution. From left to right, top row: Planck 28.5; 44.1; 70.3 and 100 GHz. Bottom row: Planck 143 and 857 GHz; 1.4 GHz; and Hα. The maps cover 5◦ × 5◦ centred on (l, b) = (160◦.26, −18◦.62) and have linear colour scales. The graticule has 1◦ spacing in Galactic coordinates. The FWHM of the elliptical Gaussian model used to fit the flux density in the filtered maps (see text) is shown. The strong AME is evident at 30–70 GHz. diffuse emission at angular scales 1◦ and this will be described (Cernicharo et al. 1985) having an angular extent of 6◦ × 2◦.It in Sect. 4. A 10% overall calibration uncertainty is assumed. contains a chain of six dense cores (B 5, IC 348, B 1, NGC 1333, L 1455 and L 1448) with Av > 2 mag running northeast– southwest and with two main centres of medium-mass star 3.3. WMAP data formation associated with the reflection nebulae IC 348 and WMAP 7-year data are included in the analysis (Jarosik et al. NGC 1333 near each end (see Fig. 1). As a nearby example of 2011). The data span 23–94 GHz and thus complement Planck low-to-medium mass star formation it has been extensively stud- data, particularly the K-band (22.8 GHz) channel. We use the 1◦- ied, with the Spitzer c2d Legacy Program (Evans et al. 2003) smoothed maps available from the LAMBDA website. We apply and the COMPLETE survey (Ridge et al. 2003), which provide colour corrections to the central frequencies using the recipe de- a wealth of infrared and spectral line data to understand condi- scribed by Jarosik et al. (2003) and assume a conservative 3% tions throughout the cloud. overall calibration uncertainty. Anomalous emission was first detected as an excess in this region of the sky via the Tenerife experiment (Davies et al. 1987)at5◦ angular scales, and it contributed most of the 3.4. DIRBE data dust-correlated signal found by de Oliveira-Costa et al. (1999). ◦ To sample the peak of the blackbody curve for tempera- With the 1 beam and 11–17 GHz frequency bands of the tures 15 K we include the COBE-DIRBE data at 240 μm COSMOSOMAS experiment (Gallegos et al. 2001) it was pos- (1249 GHz), 140 μm (2141 GHz) and 100 μm (2997 GHz). The sible to locate and identify the Perseus molecular cloud as the DIRBE data are the Zodi-Subtracted Mission Average (ZSMA) strong source of this emission (Watson et al. 2005). Follow-up maps (Hauser et al. 1998) regridded into the HEALPix format. observations with the VSA interferometer (Tibbs et al. 2010)at  Colour corrections were applied as described in the DIRBE ex- 10 resolution detected five “hot-spots” of anomalous emission, planatory supplement version 2.3. Data at higher frequencies one of which corresponds to the dense region (Av  10) around are not considered, since these will be dominated by transiently the IC 348 core. The other four hot-spots are associated with heated grains which are not in thermal equilibrium with the in- nearby dense dust knots that lie on a warm ring of dust emission  terstellar radiation field and therefore cannot be modelled with (G159.6−18.5), about 90 to the south-east of IC 348. This ring, a single modified blackbody curve. Furthermore, at wavelengths which is seen at 100 μm, is thought to be warmed by the OB 40 μm the spectrum is dominated by many emission/absorption star HD 278942 (de Zeeuw et al. 1999), which is also responsi- lines. ble for a weak H ii region in the centre of the ring, and which can be seen in the low-frequency continuum surveys with an in- tegrated flux density of a few janskys. The AME appears to have 4. Perseus molecular cloud little or no polarisation, with an upper limit of 1.8% at 33 GHz (Lopez-Caraballo et al. 2011). 4.1. Introduction and maps Planck maps of the Perseus molecular cloud region, cover- The Perseus molecular cloud complex is a relatively nearby gi- ing 30–857 GHz, are shown in Fig. 1, along with the 1.4 GHz 4 ant molecular cloud of 1.3 × 10 M at a distance of 260 pc and Hα maps. For display purposes the Planck maps have

A20, page 4 of 17 Planck Collaboration: Planck early results. XX.

10.9 GHz 28.5 GHz 70.3 GHz 857 GHz

NGC1499

AME−G160.26−18.62

Fig. 2. 16◦ × 16◦ maps after filtering with the COSMOSOMAS data reduction pipeline. The maps are centred at (l, b) = (160◦, −15◦)andthe graticule spacing is 2◦ in Galactic coordinates. From left to right: COSMOSOMAS 10.9 GHz; Planck 28.5 GHz; 70.3 GHz; and 857 GHz. The elliptical Gaussian model used to calculate the spectrum is indicated. been CMB-subtracted as described in Zacchei et al. (2011)and elliptical Gaussian (1◦.6 × 1◦.0 with major-axis PA 51◦), centred Planck HFI Core Team (2011b) and the units are thermodynamic at (RA, Dec) = (55◦.33, +31◦.67) or (l, b) = (160◦.26, −18◦.62), kelvins. The strong dust-correlated AME at 30–70 GHz is evi- following Watson et al. (2005); this gives a solid angle Ω= dent; it has no obvious counterpart at 1.4 GHz but correlates well 5.52 × 10−4 sr. We will refer to this as AME-G160.26−18.62. with the higher (>100 GHz) Planck frequencies, which are dom- There is a flat-spectrum radio quasar, 4C+32.14 (NRAO 140), inated by thermal dust. The majority of the AME comes from within 1◦, with a flux density of ≈1–3 Jy, which is seen in the northeast end of the molecular cloud, which harbours the the Planck Early Release Compact Source Catalogue (ERCSC) dense IC 348 reflection nebula. Some weak free-free emission is at frequencies up to 143 GHz (Planck Collaboration 2011c). seen to the east and the contribution from the extragalactic radio We use the flux densities and beam widths taken from the source 4C+32.14 is clear and must be removed (Sect. 4.2). The ERCSC and subtract the two-dimensional Gaussian profile from Hα (and 1.4 GHz emission) morphology is not co-located with the maps before smoothing and filtering. Meanwhile in the the emission seen at 30–70 GHz, but is likely to be significantly COSMOSOMAS maps we must extrapolate the flux density of absorbed by the intervening dust. Nevertheless, the brightest Hα 4C+32.14 to the appropriate frequencies and simulate its con- pixels predict ≈0.25 mK at 30 GHz (Dickinson et al. 2003), a tribution to the filtered map; we assume 2.0 Jy at 11–19 GHz, factor of ten below the AME at 30 GHz. It is interesting that which is in good agreement with NVSS (Condon et al. 1998) there is considerable bright structure observed in the Planck HFI and recent 2.7 GHz measurements made with the Effelsberg tele- maps, which is not always associated with strong AME. This is scope (Reich & Reich 2009). The contribution of 4C+32.14 is likely to be a consequence of a different environment and dust ≈5% of the integrated flux density, so the final results are not grain population. very sensitive to the exact treatment of the subtraction or the variability of the source. Both the smoothing and filtering result in correlated noise between the pixels in the map. A correlation 4.2. Spectrum analysis (e.g., Davies et al. 2006) is used to calculate the flux To reduce correlated 1/ f noise from the atmosphere/receiver the density taking into account the correlated noise. The covariance COSMOSOMAS time-ordered data have been filtered with the matrix is estimated from the size of the smoothing kernel and suppression of the first seven harmonics in the FFT of the circu- the point source response to the scan strategy. A simulated map lar scans. To compare with flux densities from other instruments is produced for the source, with an integrated flux density of one must either correct for the flux lost in filtering or filter all 1 Jy, and a maximum likelihood method is used to find the am- plitude and error for the given frequency channel data within a the other data in the same way. With extended sources such as × ◦ × ◦ Perseus there is uncertainty in estimating a single correction fac- 19 19 pixel box (6.3 6.3) around the source position. tor to take into account the filtering process, which in turn di- As a control, we tested the simulation and filtering pipeline rectly affects the shape of the frequency spectrum. Filtering all by applying the technique to the extended California nebula H ii the data in the same way retains the spectral shape and ensures region (NGC 1499), approximately 6◦ to the north of Perseus that the spatial frequencies used to form the spectrum are com- at (l, b) = (160◦.60, −12◦.05), as seen in Fig. 2. The emission mon between the different data sets. at radio/microwave wavelengths is dominated by optically thin The first step is to smooth all the data to the same resolu- free-free emission. Flux densities at frequencies >100 GHz were tion as the lowest frequency channel of the COSMOSOMAS affected by dust emission lying adjacent to the north-west of data, i.e., a FWHM of 1◦.12. Each data set is then fil- NGC 1499 and thus were not included. The spectrum is shown tered using the COSMOSOMAS reduction pipeline and binned in Fig. 3. A simple power-law fit, with a best-fitting spectral into a map as was done with the original COSMOSOMAS index α = −0.05 ± 0.01, is overplotted. The fit is relatively data (Fernández-Cerezo et al. 2006; Hildebrandt et al. 2007). good (χ2/d.o.f. = 1.23) but there is some evidence for AME at Example maps at 10.9, 28.5, 70.3, and 857 GHz are shown in frequencies ≈30–100 GHz. Furthermore, the spectral index is Fig. 2. To estimate the flux density at each frequency the fil- slightly flatter than the value expected for optically thin free- tered maps were compared to a filtered model map, which was free emission over this frequency range, α ≈−0.1to−0.15 passed through the same scan strategy. For the model we use an (Dickinson et al. 2003). Indeed, AME has been observed to be

A20, page 5 of 17 A&A 536, A20 (2011)

103 Table 2. Flux densities for AME-G160.26−18.62 in the Perseus molec- Model Ancillary data ular cloud and residual fluxes when free-free, CMB and thermal dust WMAP components are removed. Planck

Frequency Flux density Flux density residual [GHz] [Jy] [Jy] 0.408 9.7 ± 3.30.2 ± 3.4 2 ± ± 10 0.82 9.4 4.70.5 4.8 1.42 8.0 ± 1.7 −0.5 ± 1.8 10.9 16.1 ± 1.89.4 ± 2.0 Flux density [Jy] 12.7 20.0 ± 2.213.3 ± 2.4 14.7 28.4 ± 3.121.8 ± 3.3 16.3 35.8 ± 4.029.3 ± 4.2 22.8 39.8 ± 4.233.4 ± 4.5 28.5 38.1 ± 4.031.5 ± 4.4 101 33.0 36.7 ± 3.929.7 ± 4.4 1 10 100 40.9 30.8 ± 3.322.6 ± 4.0 Frequency [GHz] 44.1 28.5 ± 3.219.6 ± 4.0 61.3 27.4 ± 3.410.9 ± 6.4 Fig. 3. Spectrum of the California nebula (NGC 1499), measured using 70.3 32.2 ± 3.98.9 ± 9.0 the filtered flux method (see text). A simple power-law fit is shown. 93.8 63.0 ± 7.87± 21 100 78 ± 15 10 ± 29 143 202 ± 22 −25 ± 80 Model 217 1050 ± 130 120 ± 320 5 Ancillary data 10 WMAP 353 3860 ± 470 −600 ± 1400 Planck 545 15 300 ± 2100 −800 ± 4900 ± − ± 4 857 48 700 6100 1000 14 000 10 1250 93 000 ± 13 000 −2400 ± 28 000 2143 117 000 ± 15 000 7000 ± 35 000 2997 53 600 ± 6700 −2000 ± 20 000 103 Notes. The fitted spinning dust model consists of two components. 2

Flux density [Jy] 10 frequencies (>100 GHz) are dominated by thermal dust emis- 101 sion. The excess at ≈10–70 GHz is evident and has a peaked (convex) spectrum, with a maximum at 25 GHz. The data allow us to fit a multi-component parametric 100 1 10 100 1000 model to the flux density spectrum. The model consists of Frequency [GHz] four components: free-free emission; thermal dust emission; a CMB fluctuation; and spinning dust emission. The sum is − Fig. 4. Spectrum of AME-G160.26 18.62 in the Perseus molecular = + + + cloud. The best-fitting model consisting of free-free (orange dashed S S ff S td S cmb S sp. (1) line), spinning dust, and thermal dust (light blue dashed line) is shown. The free-free flux density, S ff, is calculated from the brightness The two-component spinning dust model consists of high density Tff τff molecular gas (magenta dot-dashed line) and low density atomic gas temperature, , based on the optical depth, , using the stan- (green dotted line). dard formulae: 2 2kTffΩν S ff = , (2) c2 ii significant in some H regions (Dickinson et al. 2009; Todorovic´ Ω et al. 2010). We confirmed the spectral shape with a standard where k is the Boltzmann constant, is the solid angle, and ν is aperture photometry analysis (as described in Sect. 5) and will the frequency, investigate this in more detail in a future publication. Note that −τff Tff = Te(1 − e ), (3) the excess at 30–100 GHz disappears when the spectrum is cal- culated away from the nearby dust feature to the north-west of and the optical depth, τff ,isgivenby NGC 1499. Given the consistency of the flux densities between = × −2 −1.5 −2 the frequencies, and that we have obtained similar results using τff 3.014 10 Te ν EMgff, (4) standard aperture photometry, we are confident that the spectrum where Te is the electron temperature (in units of K), EM is the calculated using the filtered maps is robust. −6 − emission measure (in units of cm pc) and gff is the Gaunt fac- The spectrum for AME-G160.26 18.62 in Perseus is shown tor, which is approximated as in Fig. 4. The flux densities and associated errors are listed   in Table 2. The spectrum is well sampled across the radio, 4.955 × 10−2 gff = ln + 1.5ln(T ). (5) microwave and FIR regimes and is a significant improvement ν/GHz e on that presented in Watson et al. (2005), with the additional Planck data allowing a much more accurate spectrum of the It is the Gaunt factor that results in a slight steepening of the free- AME to be extracted. The low-frequency (<2 GHz) data show free spectrum with frequency, which is particularly pronounced a flat spectrum consistent with free-free emission while the high at 100 GHz.

A20, page 6 of 17 Planck Collaboration: Planck early results. XX.

The thermal dust is modelled as a single-component, modi- 102 β Spinning dust model fied blackbody curve, ν B(ν, Td), which we normalise using an Ancillary data optical depth at 250 μm(1.2THz),τ ,as WMAP 250 Planck

3 ν 1 β S td = 2 h τ250 (ν/1.2THz) Ω, (6) c2 ehν/kTd − 1 101 where h is the Planck constant and β is the emissivity index. The CMB anisotropy is evaluated as the differential of the black- body temperature function, in thermodynamic units (ΔTCMB), converted to flux density units, Residual flux density [Jy]   2kΩν2 S = ΔT . 0 cmb 2 CMB (7) 10 c 1 10 100 Frequency [GHz] The spinning dust model is based on an assumed theoretical −1 2 −1 emissivity curve, jν (in units of Jy sr cm H ), normalised by Fig. 5. Spectrum of AME-G160.26−18.62 in the Perseus molecular the average column density per H nucleon, NH, and solid an- cloud after subtracting the best-fit free-free, CMB and thermal dust gle, Ω, components. A theoretical spinning dust model, consisting of two com- ponents is shown as the black solid line; the magenta dot-dashed line = Ω is for high density molecular gas and the green dotted line is for low S sp NH jν . (8) density atomic gas, which is a very small contribution in this case (see Sect. 6). As described in Sect. 6,weusetheSPDUSTcode(Ali-Haïmoud et al. 2009) to model the spinning dust spectrum. The best-fitting model is plotted in Fig. 4. The fitted pa- 5. ρ Ophiuchi molecular cloud ff rameters are EM for free-free, di erential temperature of CMB 5.1. Introduction and maps perturbation (ΔTCMB), dust temperature (Td), dust emissivity in- dex (β) and 250 μm opacity (τ250). The spinning dust model The ρ Ophiuchi molecular cloud (e.g., Encrenaz 1974; Kulesa consists of two components, to represent dust associated with et al. 2005; Young et al. 2006) lies in the Gould Belt of the dense molecular gas and low density atomic gas, and will be dis- closest molecular complexes, at a distance D = 135 ± 15 pc cussed in Sect. 6. The electron temperature is assumed to have (Perryman et al. 1997). It is undergoing intermediate-mass star a typical value of 8000 K for gas in the solar neighbourhood. A formation. Ultra-violet radiation from its hottest young stars conservative 3% uncertainty in the CO correction is assumed at heats and dissociates exposed layers, but does not ionise bulk hy- 100 GHz and is not included in the fit. Similarly, the 217 GHz drogen. The most prominent PDR in ρ Ophiuchi is the ρ Oph W 4 5 −3 value is not included in the fit, due to the possible small con- filament, with nH = 10 –10 cm , which has been studied by tamination from the CO J = 2 → 1 line. From the initial spec- Liseau et al. (1999)andHabart et al. (2003). ρ Oph W is excited tral fit the spectral index is found at each frequency, which is by HD 147889 (a B2, 3 IV binary), the earliest star in the ρ Oph then used to calculate the appropriate colour corrections for each star formation complex. The line of sight to HD 147889 can be band in an iterative manner. After applying the colour correc- taken to characterise the physical conditions in the ambient dif- 20 −2 −3 tions (typically a few per cent correction), the spectrum is fit- fuse cloud, with NH ∼ 2 × 10 cm , nH ∼ 400 cm ; no accu- ted again to obtain the final flux densities, which are given in rate values are available for NH2 (van Dishoeck & Black 1989; Table 2. The uncertainty in the fitted model for the spectrum is Ka´zmierczak et al. 2010). The bulk of the mass in the ρ Ophiuchi computed by propagating the errors of each of the parameters cloud is situated in the Oph A molecular core (e.g., Wootten et al. and combining them with the other sources of error, i.e., back- 1978, 1979), located roughly 10 southwest of ρ Oph W. ground/noise residuals and an overall calibration error. The fit Only faint radiation from the Rayleigh-Jeans tail of is very good with χ2/d.o.f. = 0.95. The best-fitting parameters ∼10–100 K dust is expected at wavelengths longwards of −6 are: EM = 148 ± 13.5cm pc for free-free; Td = 18.5 ± 0.6K ∼3 mm. Yet CBI observations revealed that the ρ Oph W PDR and β = 1.65 ± 0.08 for the thermal dust; and a negligible CMB is surprisingly bright at centimetre wavelengths (Casassus et al. contribution of ΔTCMB = −6 ± 66 μK. 2008). The WMAP 33 GHz image confirms that the general lo- In Fig. 5 we show the residual spectrum after subtracting the cation of the signal seen by CBI is offset by 10 from the centroid free-free, CMB and thermal dust models. The additional uncer- of the WMAP 94 GHz data, located on the Oph A core. The sig- tainty from the subtraction of the models is added in quadrature nal seen at centimetre wavelengths has no 5 GHz counterpart; to the flux density errors to find the final uncertainties on the instead the optically thin free-free emission is located in the im- residual flux densities, which are listed in Table 2.Thisrepre- mediate vicinity of HD 147889. The other early-type stars in the sents the most precise spectrum of AME to date. The weighted complex are S 1 (B4 V) and SR 3 (B6 V), and each is surrounded average over the range 10 to 94 GHz gives a detection signifi- by bright IR emission peaks, which lack any radio counterpart. cance level of 17.1σ. The peak of the spinning dust component Planck maps of the ρ Ophiuchi molecular cloud region cov- is centred at ≈25 GHz and the low and high frequency sides are ering 30–857 GHz, along with the 1.4 GHz and Hα maps, are well-defined in the range 10–90 GHz. A physical model for the shown in Fig. 6. Strong AME is evident at 28.5–44.1 GHz, with spinning dust is presented in Sect. 6. It consists of two compo- the bulk of the emission falling within 0◦.5 of the central posi- nents (atomic and molecular), which are overplotted in Fig. 5, tion. For display purposes the Planck maps have been CMB- and is seen to be an excellent fit to the data. subtracted as described in Zacchei et al. (2011)andPlanck HFI

A20, page 7 of 17 A&A 536, A20 (2011)

28.5 GHz44.1 GHz 70.3 GHz 100 GHz

AME−G353.05+16.90

143 GHz 857 GHz 1.4 GHz H−alpha

Fig. 6. Maps of the ρ Ophiuchi molecular cloud region at their original angular resolution. From left to right, top row: Planck 28.5; 44.1; 70.3; and 100 GHz. Bottom row: Planck 143 and 857 GHz; 1.4 GHz; and Hα. The maps cover 5◦ × 5◦ centred on (l, b) = (353◦.05, +16◦.90) and have linear colour scales. The graticule has 1◦ spacing in Galactic coordinates. The circular aperture and background annulus, which were used to calculate the flux density (see text), are indicated. The strong AME is evident at 28.5 and 44.1 GHz.

Core Team (2011b) and the units are thermodynamic kelvins. background fluctuations. We also compared two halves of the The H ii region to the west is also clearly detected in the data as a jack-knife test. Results were found to be consistent to Planck LFI maps, but is well separated from the bulk of the within a fraction of the derived uncertainties. AME. In the location of the strongest AME, there is a valley We tested the aperture photometry by measuring the spec- in intensity seen at low frequencies, supported by little Hα emis- trum of well-known, bright objects. We chose bright H ii regions, sion. The AME is well correlated with Planck HFI maps, indi- which are expected to be dominated by thermal bremsstrahlung cating an association with dust. (free-free) radiation at radio frequencies and thermal dust in the sub-mm/FIR. Figure 7 shows the spectrum of the Orion nebula (M42) with an aperture of radius 60 and background annulus 5.2. Spectrum between radii of 80 and 120. At frequencies below a few gi- − ◦ gahertz, the free-free emission is becoming optically thick, with The ρ Ophiuchi molecular cloud (at declination 24 ) lies out- ≈ ∼ side the range of the COSMOSOMAS survey. We therefore a turnover frequency of 2 GHz. At frequencies 10–60 GHz, chose to use simple aperture photometry on the unfiltered maps the spectrum is well-fitted by a simple optically thin free-free spectrum (α ≈−0.1). The WMAP and Planck data points are to calculate the spectrum. The maps were first smoothed to a ≈ common resolution of 1◦ and converted to units of Jy pixel−1. in good agreement to within 1%; the residuals are consistent with zero and within a fraction of the assumed uncertainties. At The brightness in each pixel was then summed over a given  aperture to give the flux density. We chose a circular aperture frequencies 100 GHz, thermal dust emission dominates and is = ◦ ◦  well described by a modified blackbody with an emissivity index centred at (l, b) (353.05, 16.90) with a radius of 60 ; this cor- = ± = ± responds to a solid angle Ω=9.57×10−4 sr. We will refer to this β 1.71 0.12 and dust temperature Td 25.7 1.7K. source as AME-G353.05+16.90. An estimate of the background We note that the lack of AME in M42 is likely to be due brightness is subtracted using the median value in a circular an- to the difference in conditions (e.g., intense radiation field) that nulus between an inner radius of 80 and an outer radius of 120, may result in depletion of the smallest dust grains, which would as indicated in Fig. 6. This removes local background emis- result in a large reduction in the amount of spinning dust emis- sion and any residual offsets in the maps. Errors are estimated sion. This is a different situation to the more diffuse warm from the standard deviation in the background annulus, which ionised medium (WIM) studied by Dobler et al. (2009), where was verified with simulations of random sources of known flux they detected AME in the WIM at the level of ≈20% relative to density that were injected into the sky maps. The results were the free-free emission. not strongly sensitive to the exact choice of background annu- Figure 8 shows the spectrum of AME-G353.03+16.90 in lus; all values were within the 1σ uncertainties. This indicates the ρ Oph West region; the flux densities are listed in Table 3. that the emission in the aperture is bright relative to the local The same model used for the Perseus region (Sect. 4.2)was

A20, page 8 of 17 Planck Collaboration: Planck early results. XX.

106 Table 3. Flux densities for AME-G353.05+16.90 in the ρ Ophiuchi Model Ancillary data molecular cloud and residual flux densities when free-free, CMB and WMAP thermal dust components are removed. Planck 5 10 100 Frequency Flux density Flux density residual 50 [GHz] [Jy] [Jy] 0.4 −6.8 ± 9.4 −7 ± 11 0 104 1.4 −1.0 ± 6.8 −1.0 ± 8.3 -50 2.3 1.3 ± 6.91.2 ± 8.3 22.7 26.3 ± 5.524.8 ± 6.6 Residual flux density [Jy] Flux density [Jy] -100 10 100 28.5 29.6 ± 5.627.3 ± 6.6 3 10 Frequency [GHz] 33.0 30.7 ± 5.327.2 ± 6.3 40.7 27.7 ± 4.621.9 ± 5.8 44.1 27.0 ± 4.419.9 ± 5.7 60.6 26.3 ± 4.59.8 ± 6.5 2 10 70.3 31.4 ± 5.16.5 ± 7.8 1 10 100 1000 93.4 63.6 ± 8.96± 15 Frequency [GHz] 100 84 ± 13 13 ± 19 143 222 ± 31 5 ± 49 Fig. 7. Spectrum of the Orion nebula (M 42) H ii region using aperture 217 989 ± 16 130 ± 200 photometry. The model consists of optically thin free-free emission (or- 353 4100 ± 670 −180 ± 780 ange dashed line) and thermal dust emission (light blue dashed line). 545 17 100 ± 2900 400 ± 3200 The inner panel shows the residual spectrum after removal of the model, 857 58 300 ± 8800 500 ± 10 000 indicating the consistency between Planck and WMAP data. 1249 122 000 ± 22 000 −5000 ± 25 000 2141 218 000 ± 37 000 20 000 ± 41 000 106 2997 125 000 ± 26 000 −9000 ± 32 000 Model 5 Ancillary data 10 WMAP Planck 104 As expected, the impact was considerable. In particular, the ther- mal dust model is less well-constrained when the HFI data are 103 not included. For example, the uncertainty in the dust tempera- ture increases to ±2.5 K, while the uncertainty in the emissivity 2 10 index increases to ±0.36. This is crucial for making precise mea- surements of the AME at frequencies, particularly above 50 GHz 101 Flux density [Jy] where the subtraction of thermal dust is critical. 100 In Fig. 9 we show the residual spectrum after subtracting the free-free, CMB, and thermal dust models. The additional er- 10-1 ror from the subtraction of the models is added in quadrature to the flux density errors to obtain the uncertainty on the resid- 10-2 1 10 100 1000 ual flux densities, listed in Table 3. As in Perseus, the residual Frequency [GHz] spectrum has a clearly defined convex shape, peaking at a fre- quency of 30 GHz; the excess is significant at the 8.4σ level. Fig. 8. Spectrum of AME-G353.05+16.90 in the ρ Ophiuchi West The spectrum at 50–100 GHz is flatter than a single spinning molecular cloud. The best-fitting model consisting of free-free (or- dust model can easily account for, suggesting a different envi- ange dashed line), spinning dust, CMB (black double/triple-dotted line), ronment or distribution of dust grains. This may also be due to and thermal dust (light blue dashed line), is shown. The spinning dust variations in the dust emissivity index attributable to multiple model consists of two components: high density molecular gas (ma- dust components (e.g., Finkbeiner et al. 1999). To test this possi- genta dot-dashed line); and low density atomic gas (green dotted line). bility would require the inclusion of the 100/217 GHz data after / The 100 217 GHz data are contaminated by CO line emission and are careful subtraction of the CO line contamination. not included in the fit. A theoretical spinning dust model consisting of two com- ponents (associated with atomic and molecular gas) is over- plotted and provides an excellent fit to the data. Note that the fitted to the data, except for the details of the spinning dust denser molecular gas represents the dominant AME with peak at component, which will be discussed in Sect. 6. Significant CO ≈30 GHz, while the irradiated low density atomic gas accounts line contamination is visible at 100 and 217 GHz, so these two for the 50–100 GHz part of the spectrum. The physical model is bands were not included in the fit. There is minimal free-free presented in Sect. 6. emission within the aperture, and hence the error at low (0.4– 2.3 GHz) frequencies is dominated by fluctuations in the lo- cal background, as indicated by the large error in the 2.3 GHz 6. Modelling the anomalous microwave emission / flux density. Note that the 408 1420 MHz data were included in with spinning dust the fit but were consistent with zero. The best-fitting parameters −6 are: EM = 1±48 cm pc; Td = 20.7±1.6K;β = 1.75±0.18; and Comparisons of the AME with dust emission in localised re- CMB contribution ΔTCMB = 82 ± 61 μK. The fit of the overall gions (Casassus et al. 2006; Scaife et al. 2009; Tibbs et al. model is very good (χ2/d.o.f. = 0.29). 2010; Castellanos et al. 2011) and in all-sky surveys (Lagache As a test of the impact that the Planck data was having on 2003; Ysard et al. 2010) suggest small dust grains (poly- the spectrum, we repeated the analysis omitting the Planck data. cyclic aromatic hydrocarbons, PAH) as carriers of this emission,

A20, page 9 of 17 A&A 536, A20 (2011)

103 Table 4. Spinning dust model Parameters in the spinning dust models for Perseus and Ancillary data ρ Ophiuchi. WMAP Planck 102 Gas state Molecular Atomic Ionised Perseus 21 −2 NH [10 cm ] 11.7 1.3 0.4 −3 1 nH [cm ] 250 30 1 10 z [pc] 15.1 14.0 . . . G0 1 2 ... T [K] 40 100 8 × 103 6 0 xH [ppm] 112 410 10 Residual flux density [Jy] 10 xC [ppm] <1 100 . . . y 1 0.1 ... a0 [nm] 0.58 0.53 . . . 10-1 bC [ppm] 68 68 . . . 1 10 100 β ... 1.65 ... Frequency [GHz] Td [K] ... 18.5 ... −4 τ250 ... 9.4 × 10 ... Fig. 9. Spectrum of AME-G353.05+16.90 in the ρ Ophiuchi molecu- ρ Ophiuchi 21 −2 lar cloud after subtracting the best-fit free-free, CMB and thermal dust NH [10 cm ] 18.2 0.4 0.4 −3 4 components. A theoretical spinning dust model consisting of two com- nH [cm ]2× 10 200 0.5 ponents is shown as the black solid line; the magenta dot-dashed line is z [pc] 0.3 0.6 ... for dense molecular gas and the green dotted line is for irradiated low G0 0.4 400 . . . density atomic gas (see Sect. 6). The 100/217 GHz data are contami- T [K] 20 103 8 × 103 6 nated by CO line emission and are not included in the fit. xH [ppm] 9.2 373 10 xC [ppm] <1 100 . . . y 1 0.1 ... a [nm] 0.60 0.38 . . . supporting the spinning dust explanation (Draine & Lazarian 0 bC [ppm] 65 50 . . . 1998). The present and previous measurements of the AME β ... 1.75 ... show that this component peaks at 20–40 GHz with a brightness − Td [K] ... 20.7 ... 1 −3 of 0.01 to 0.1 MJy sr (see references in Sect. 1). Current spin- τ250 ... 3.2 × 10 ... ning dust models (Ali-Haïmoud et al. 2009; Ysard & Verstraete + + 2010; Hoang et al. 2010; Silsbee et al. 2011) indicate that low- Notes. xH and xC are the abundances of H and C ; y is the molecular −3 n /n a b density (nH ≤ 30 cm ), atomic diffuse gas peaks in this fre- fraction 2 H2 H; 0 is the centroid of the PAH size distribution; C is quency range. Therefore the AME has often been explained with the PAH abundance of carbon in PAHs; and z is the depth of the emis- spinning dust emission arising mostly from regions of low den- sion region along the line of sight. The radiation field is that of Mathis G sity gas, such as from the Cold Neutral Medium (CNM) or Warm et al. (1983) multiplied by the factor 0. The ionised gas accounts for the free-free emission but does not contribute to the spinning dust emis- Neutral Medium (WNM). The column density NH associated sion. with the AME can be derived from observations provided the emissivity of spinning dust is known: adopting an average of 10−23 MJy sr−1 cm2 H−1 from theoretical models, we find that 22 −2 −1 3 NH ∼ 10 cm for a flux density level of 0.1 MJy sr .Fora the size of PAHs . The largest PAHs have the lowest emissivity. −3 −3 gas density of 30 cm , the extent of the emitting region is thus In low temperature gas (G0 < 3andnH > 30 cm ), the gas- NH/nH ∼ 100 pc, much larger than the observed size of the grain collisions dominate the excitation and the damping of the present AME clouds (<10 pc). The situation is even worse for grain rotation. Thus the spinning dust spectrum becomes sensi- −3 + + the warm ionised medium (WIM) component (nH < 1cm ). tive to the abundance of the major ions, H and C , denoted xH The spinning dust emission of these regions must therefore in- and xC respectively. To compute the spinning dust emissivity, we clude a significant contribution from gas which is denser by use SPDUST4 in the case 2 of Silsbee et al. (2011). Rather than fit- at least a factor of 10, as suggested by Casassus et al. (2008) ting all the spinning dust parameters on the observed spectrum, and Castellanos et al. (2011). In the preliminary models of the we derive realistic values from observations, and fit the ampli- Perseus and ρ Ophiuchi regions presented below, we show that, tude to the spectrum. We derive xH from the ionisation balance assuming plausible physical conditions and PAH size distribu- of hydrogen assuming a standard cosmic-ray ionisation rate of tions, spinning dust emission from dense gas can explain most 5 × 10−17 s−1 H−1 (Williams et al. 1998; Wolfire et al. 2003). of the AME. Conversely, xC cannot be estimated so readily because C under- The rotation of PAHs in space is governed by IR emis- goes reactions with H2 and PAHs (Röllig et al. 2006; Wolfire sion and gas-grain interactions (Draine & Lazarian 1998). et al. 2008); we therefore take xC as a free parameter. The size Quantitative modelling of spinning dust emission requires the distribution of PAHs is assumed to be a log-normal of centroid knowledge of many (∼10) parameters describing the state of a0 and width σ = 0.4. The abundance of PAHs is set by bC,the the gas, the radiation field, and the PAH size distribution, as abundance of carbon locked up in PAHs. well as the electric dipole moment. As shown by Ali-Haïmoud et al. (2009)andYsard et al. (2010), the peak frequency of spin- 3 We assume here the electric dipole moment to be as in Draine & ning dust emission in neutral gas (and for radiation intensity Lazarian (1998), a prescription also shown to be compatible with the 2 G0 < 10 ,whereG0 is the multiplying factor relative to the AME extracted from WMAP data (Ysard et al. 2010). mean radiation intensity in the interstellar medium, as defined 4 http://www.tapir.caltech.edu/~yacine/spdust/spdust. by Mathis et al. 1983) depends chiefly on the gas density and on html

A20, page 10 of 17 Planck Collaboration: Planck early results. XX.

Now we describe the physical conditions within our two less abundant PAHs in stronger radiation fields, in line with re- regions, the Perseus and ρ Ophiuchi molecular clouds, as de- cent work suggesting that PAHs are photoevaporated from larger rived from observations, and how we take them into account carbonaceous grains (Berné et al. 2007; Compiègne et al. 2008; in our modelling (see Table 4). For both regions, the free- Velusamy & Langer 2008). The PAH abundances we assume for free contribution is that of the H ii region and the thermal both regions are within the range found by Habart et al. (2003). dust emission is described as a modified blackbody, Iν = The fit is shown in Fig. 9 and the parameters of our model are −β τ250 (λ/250 μm) Bν(T)whereτ250 is the dust opacity at summarised in Table 4. 250 μm, T the dust temperature (see Table 4), and Bν is the We note, however, that the determination of a0, G0 and nH blackbody brightness. We do not include the contribution of the is degenerate in the present modelling; you can create almost ionised gas to the spinning dust emission because of its low identical spinning dust curves by varying these parameters in density and also because PAHs have been shown to be destroyed different ways (see, e.g., Ali-Haïmoud et al. 2009). Lifting this efficiently in such environments (Lebouteiller et al. 2007, 2011). degeneracy involves treating the radiative transfer along with the In addition, we note that the current SPDUST code does not in- spinning dust motion. The present results are therefore sugges- clude ionising photons (hν ≥ 13.6 eV), which are important for tive only and we emphasise that the physical conditions of our the PAH charge and IR emission. models are derived from scales much smaller than the resolu- tion of our observations. The previous discussion also indicates that quantitative model spectra can only be obtained from a con- 6.1. Perseus molecular cloud sistent treatment of the gas state, the radiative transfer, and the At the degree scale, the Perseus line of sight intersects a spinning motion of dust grains. Finally, a combination of multi- molecular cloud surrounding the H ii region associated with the frequency data, including detailed IR measurements, is clearly HD 278942 B0 V star. Ridge et al. (2006) argued that this ionised needed to extract detailed information about the environment shell has a density of ∼1cm−3. The gas density and temper- and the dust grains. ature of the molecular component are taken from the analysis 5 of observed C2 lines carried out by Iglesias-Groth (2011) .The parameters of our model, which provides a good match to the 7. New regions of anomalous emission AME, are given in Table 4, and the corresponding emission is showninFig.5. We note that the depth, z, of the emissive re- We have searched the Planck maps for new regions of AME gion along the sightline is compatible with determinations at on the scale of ∼3◦ or smaller. To do this, one must suppress smaller scales (Ridge et al. 2006). The spinning dust emission is the other emission components that contribute to the frequencies dominated by the contribution of the dense gas; the low-density where AME is the strongest (∼20–60 GHz). Numerous methods atomic gas appears to play only a minor role. This is consistent can be used (e.g., Leach et al. 2008; Delabrouille & Cardoso with the results of Planck Collaboration (2011q), which show 2009). However, we have taken a simplistic approach to identify that, on large scales, there is bright AME associated with molec- bright AME regions. We subtracted a simple spatial model of ular clouds traced by CO everywhere in the Galaxy, while the known emission mechanisms (synchrotron, free-free, and ther- association with the diffuse and atomic phase correlated with H i mal dust) using extrapolations of existing templates from obser- is less clear. vations or theoretical expectations, and then examined the loca- In contrast to previous work where xH = 0 was assumed in tions of the remaining emission in more detail. We used maps −5 ◦ = molecular gas, our ionisation balance of H gives xH ∼ 3 × 10 , smoothed to 1 , with HEALPix Nside 256 (Górski et al. 2005). a change particularly important for Perseus. We take xC to be We subtracted the synchrotron emission from the CMB- 1 ppm: at this level the C+ ion has no influence on the spin- subtracted Planck maps using the 408 MHz map (Haslam et al. ning dust emission. Adopting larger sizes for the PAHs (a0 = 1982), while the 1420 MHz (Reich 1982; Reich & Reich 1986; 0.58 nm) than those used in previous spinning dust models Reich et al. 2001) and the 2326 MHz (Jonas et al. 1998)maps (but comparable to the diffuse interstellar medium (0.64 nm); were used in conjunction with 408 MHz data to determine the Compiègne et al. 2011), the molecular gas accounts well for the spectral index. Brightness temperature spectral indices were re- AME. stricted to lie in the range −2.0 >β>−4.0(T ∝ ν β). For the gap region where no 1420/2326MHz data exist, we assume β = −3.0. We used versions of the maps that have had no point 6.2. The ρ Ophiuchi molecular cloud sources subtracted (Platania et al. 2003) so as to allow direct The ρ Ophiuchi spectrum includes contributions from a photo- comparison with the Planck data. For free-free emission, we dissociated interface region (PDR) and the ionised gas associ- use the Hα map of Dickinson et al. (2003), which is based on ff ated with the star HD 147889. We note that IR observations have WHAM (Ha ner et al. 2003) and SHASSA (Gaustad et al. 2001) shown that the PAH emission arises from dense irradiated inter- data corrected for dust absorption. The map is also used to cor- faces (Abergel et al. 2002; Habart et al. 2005). As constrained rect the synchrotron maps for free-free emission, assuming an by Habart et al. (2003), the PDR has a density gradient between average electron temperature of 8000 K. The Rayleigh-Jeans tail the star and the molecular cloud. For the present preliminary of the thermal (vibrational) dust component is determined from model, we represent this interface with two effective media at the Planck maps at 143 and 545 GHz, thereby avoiding signif- low (200 cm−3) and high (2 × 104 cm−3) densities irradiated by icant CO contamination. The 545 GHz data were used for nor- strong and weak radiation fields. In the more excited ρ Ophiuchi malisation and the 143 GHz data were used for the spectral index + region we are able to match the AME by requiring: (a) an at- of the power-law. This was constrained to lie between 0 and 3 tenuated radiation field in the dense gas; and (b) smaller and to avoid extrapolation artefacts in some areas at high latitudes where the uncertainty on the spectral index is large. 5 A single high resolution measurement may not be representative of The full-sky map of the residuals from the 28.5 GHz the line-of-sight on degree scales, since the density and temperature Planck LFI data is shown in Fig. 10. Artefacts from the sub- may vary considerably across the beam. traction are present. The division between the 1420 and the

A20, page 11 of 17 A&A 536, A20 (2011)

AME−G173.6+2.8 Ophiuchus

Perseus AME−G107.1+5.2

Fig. 10. Residuals in the full sky Planck LFI 28.5 GHz 1◦ smoothed map after subtraction of synchrotron, free-free and thermal dust emission (see text). 12◦.5 × 12◦.5 Cut out maps are shown for the Perseus and ρ Ophiuchi molecular clouds, and the two new regions of AME, AME-G107.1+5.2 and AME-G173.6+2.8. The graticule spacing is 5◦ in Galactic coordinates.

2326 MHz maps is evident as well as the lack of observations The spectrum of the 1◦.3 degree (COSMOSOMAS) region around the south celestial pole. Other large-scale features are centred on AME-G173.6+2.8 is shown in Fig. 12. The radio flux also present. These are extended features and hence do not affect densities are at 820 MHz (Berkhuijsen 1972), 1.4 GHz (Reich our search for compact areas of AME. A large number of AME 1982) and 11–17 GHz (COSMOSOMAS). The spectrum and candidates are evident in the map, including the strong Perseus residuals of the region are derived using the COSMOSOMAS and ρ Ophiuchi AME regions described in Sects. 4 and 5.We filtering method described in Sect. 4.2. examined ∼50 regions in detail. The spectrum is well-fitted using a combination of thermal The next step in identifying new AME regions is to con- dust (with T = 19.6±0.6Kandanemissivityindexofβ = 1.54± struct a spectrum for each candidate, using aperture photometry. 0.07), free-free emission (EM = (2.47 ± 0.13) × 103 cm−6 pc) In each case there is strong thermal dust and free-free emission and spinning dust (Draine & Lazarian 1998 CNM model with 22 −2 indicative of an associated H ii region. We now illustrate the re- NH = (2.7±0.3)×10 cm ), with a small CMB contribution of sults of our search in two of the compact regions which show ΔTCMB = 130 ± 40 μK. The residual spectrum shows the signa- excess emission in the ∼20–60 GHz range. ture of spinning dust peaking in the frequency range 20–30 GHz, detected with a significance level of 6.4σ. The main uncertainty in the AME spectrum is the level of free-free emission to be sub- 7.1. AME-G173.6+2.8 tracted. As shown, the spinning dust emission in the 20–30 GHz AME-G173.6+2.8 (05h41m, +35◦51 J2000.0) is a strong, com- region is ∼2 times the free-free emission; the uncertainty is <10– pact region of dust emission towards the Galactic anticentre. 20%. The current spectral coverage defines the peak of the spec- The H ii region S235 (Sharpless 1959) lies within this area. It trum and provides useful constraints in modelling the AME. An is a region of massive star formation with a young stellar ob- ultracompact H ii region, producing free-free emission that is op- ject at its centre, and it has several compact radio components tically thick at frequencies <10 GHz, could possibly contribute within several arcmin (Nordh et al. 1984; Felli et al. 2006). to the flux density at >10 GHz. However, it cannot easily account Radio recombination line observations indicate a radial veloc- for the shape at high (40 GHz) frequencies. ity of −24.5 ± 0.7kms−1 (Paladini et al. 2003). S235 lies at a distance of 1.8 kpc and is therefore within the Perseus spiral arm. 7.2. AME-G107.1+5.2 Maps of the region at a range of frequencies are shown at the observing resolution of each in Fig. 11. The AME is strong AME-G107.1+5.2 (22h22m, +63◦23 J2000.0) is a dust complex at 28.5 GHz while the dust distribution at the higher Planck fre- containing the star-forming H ii region S140 (Sharpless 1959) −1 quencies shows the small-scale structure of the dust within the with VLSR = −8.5 ± 1.0kms . The region lies at the edge of the 28.5 GHz beam. It can be seen that the Hα at the 6 resolution Cepheus bubble (Patel et al. 1998) at a distance of 800 pc. of Finkbeiner (2003)isoffset from the dust, most likely due to Maps of the region are shown in Fig. 13 using the same for- absorption by the dust. mat as in Fig. 11. Dust is widespread over the region but with

A20, page 12 of 17 Planck Collaboration: Planck early results. XX.

105 28.5 GHz 70.3 GHz Model Ancillary data WMAP 104 Planck

103

102 Flux density [Jy]

101

545 GHz H−alpha 100 1 10 100 1000 Frequency [GHz] 102 Spinning dust model Ancillary data WMAP Planck

101

Fig. 11. Maps of the AME-G173.6+2.8 region: 28.5 GHz (linear scale); 70.3 GHz (linear scale); 545 GHz (logarithmic scale); and Hα (loga- Residual flux density [Jy] rithmic scale). The maps are 8◦ × 8◦ and the graticule spacing is 1◦ in Galactic coordinates.

100 1 10 100 a concentration at the position of the AME. There are multiple Frequency [GHz] dust components within the 28.5 GHz beam. There is a small ff ≈  o set of the 28.5 GHz feature of 15 to lower latitudes and Fig. 12. The spectrum (top) and residuals (bottom)forAME- higher longitudes than the brightest feature in the dust maps. G173.6+2.8 after removing free-free emission (orange dashed line), This indicates a higher radio emissivity of the dust in this off- thermal dust emission (light blue dashed line) and CMB anisotropies set position. There is also an offset in the Hα emission, again (not visible). A spinning dust model is shown as a magenta dot-dashed possibly due to absorption by the associated dust. line. The spectrum estimated for AME-G107.1+5.2 and the resid- ual spectrum are shown in Fig. 14. The least-squares fit to the free-free emission (EM = 150 ± 39 cm−6 pc), thermal dust (Td = 19.0 ± 2.6K,β = 2.04 ± 0.45), spinning dust (Draine & free-free flux density to be 0.40 ± 0.11; the corresponding values 21 −2 Lazarian 1998 WIM model with NH = (1.76±0.44)×10 cm ), for AME-G173.6+2.8 and AME-G107.1+5.2 are 0.67 and 1.4, and CMB (ΔTCMB = 75 ± 38 μK) is well-defined. The AME is respectively. In the case of the ratio of 33 GHz AME to 100 μm detected at a significance of 4.8σ with a shape that is consis- thermal dust emission, the VSA survey value was 1.1 × 10−4 tent with spinning dust. However, an ultracompact H ii region (or 2 × 10−4 at the 15 GHz peak of the spectrum), while the producing optically thick free-free at frequencies below 10 GHz AME-G173.6+2.8 and the AME-G107.1+5.2 values are some- may also be contributing. what larger at 3.3 × 10−4 and 3.8 × 10−4, respectively. It is how- ever not surprising that the 30 GHz AME is stronger in the two Planck H ii regions than in the VSA H ii regions, because 7.3. Assessment of the search for new AME regions the former were selected as the brightest AME sources in the We have shown that it is possible to detect bright compact AME Planck maps. The contribution from ultracompact H ii regions regions in the Planck data using a rather simplistic selection pro- (Wood & Churchwell 1989) may also be a significant contribu- cess. In order to verify the identifications it is necessary to add tor to the AME at low Galactic latitudes. ancillary data at radio and FIR frequencies. Now that identifica- Finally, we note that the derived dust temperature for these tions have been made, new ancillary data covering a larger fre- H ii regions, at Td ≈ 19 K, is closer to that of the diffuse inter- quency range will be sought. Further development of the search stellar medium than to the warm dust (Td ∼ 30–50 K) usually algorithms will also be worthwhile. found surrounding OB stars (e.g., Wood & Churchwell 1989). The strength of AME relative to the free-free and thermal The dust observed along the line of sight at ∼1◦ resolution is dust emission in these two regions can be compared with the likely to be a mixture of populations, including higher density values for H ii regions obtained by Todorovicetal.´ (2010)in regions (e.g., molecular clouds) that could dominate the signal the VSA survey of the l = 27◦–47◦ section of the Galactic on these scales. This may explain the lower average temperature plane. The latter survey found the ratio of 33 GHz AME to and the presence of significant AME in these sightlines.

A20, page 13 of 17 A&A 536, A20 (2011)

5 Model 28.5 GHz 70.3 GHz 10 Ancillary data WMAP Planck

104

103 Flux density [Jy] 102

101

545 GHz H−alpha 1 10 100 1000 Frequency [GHz] 102 Spinning dust model Ancillary data WMAP Planck

101

100 Fig. 13. Maps of the AME-G107.1+5.2 region: 28.5 GHz (linear scale); 70.3 GHz (linear scale); 545 GHz (logarithmic scale); and Hα (loga- Residual flux density [Jy] rithmic scale). The maps are 8◦ × 8◦ and the graticule spacing is 1◦ in Galactic coordinates. 10-1 1 10 100 As selection procedures develop, we expect many AME re- Frequency [GHz] gions to be discovered in a range of physical environments, which will provide data necessary to acquire a full understanding Fig. 14. The spectrum (top) and residuals (bottom)forAME- of the AME emission mechanism. G107.1+5.2 after removing free-free (orange dashed line), thermal dust emission (light blue dashed line) and CMB anisotropies (black double- dotted line). A spinning dust model is shown as a magenta dot-dashed 8. Conclusions line. Filling in the gap between earlier radio and FIR measurements, the frequency coverage of Planck has provided a unique oppor- Low-density gas is only a minor contributor to AME in Perseus; tunity to establish comprehensive spectra of regions of AME. however, in ρ Ophiuchi, although dense gas accounts for the Planck has revealed the high frequency side of the spectral peak. peak at ≈30 GHz, irradiated low-density atomic gas may be con- The new spectra are the basis for understanding the emission tributing in the range 50–100 GHz. The picture seems to be that mechanism and the environment in which it occurs. The evi- smaller PAHs are found in PDRs (G0 > 100), as suggested by dence from the present observations strongly favours the spin- recent Spitzer observations. However, the determination of the ning dust mechanism (electric dipole radiation). Planck provides PAH size is degenerate with that of nH and G0 and quantitative a rich data set that can be used as a basis for developing a realis- conclusions will only be obtained from consistent modelling of tic understanding of the AME mechanism in a range of Galactic the gas state, radiative transfer, and spinning dust. At this level of environments. modelling it is not possible to constrain the electric dipole mo- The two best-studied AME sources that have extensive an- ment of PAHs. Future work to carry out more detailed modelling cillary data are the Perseus and ρ Ophiuchi molecular clouds. is required. In each, the spectrum is well fitted by free-free, thermal Our preliminary search for new AME regions in the dust, and spinning dust, with a small contribution from the Planck data has revealed many candidates. They were uncov- CMB. Spinning dust provides a good fit to the microwave ered by subtracting synchrotron, free-free and thermal dust emis- (10–100 GHz) part of the spectrum, which peaks at ≈30 GHz. sion based on the usual spatial templates. Two of the new can- Theoretical spinning dust spectra are presented for a physi- didate regions that show AME at >5σ significance have spectra cal model consisting of molecular and atomic states. It is pos- consistent with a spinning dust spectral shape. Additional multi- sible to derive physical parameters that are consistent with the frequency data (e.g., such as the 5 GHz C-Band All-Sky Survey, environment and still provide a good fit to the data. Using pa- C-BASS6; King et al. 2010) as well as high-resolution observa- rameters constrained at smaller angular scales, the 20–40 GHz tions will be needed to understand the detailed structure of the AME peak in Perseus is well explained with spinning dust AME in these regions. −3 emission arising from dense, molecular gas (nH > 200 cm ) subjected to a few times the typical interstellar radiation field. 6 http://www.astro.caltech.edu/cbass/

A20, page 14 of 17 Planck Collaboration: Planck early results. XX.

Acknowledgements. We thank the referee, Doug Finkbeiner, for useful com- Hildebrandt, S. R., Rebolo, R., Rubiño-Martín, J. A., et al. 2007, MNRAS, 382, ments. We thank Justin Jonas for providing the 2326 MHz HartRAO map. We 594 acknowledge the use of the MPIfR Survey Sampler website at http://www. Hoang, T., Draine, B. T., & Lazarian, A. 2010, ApJ, 715, 1462 mpifr-bonn.mpg.de/survey.html. We acknowledge the use of the Legacy Iglesias-Groth, S. 2011, MNRAS, 411, 1857 Archive for Microwave Background Data Analysis (LAMBDA); support for Jarosik, N., Barnes, C., Bennett, C. L., et al. 2003, ApJS, 148, 29 LAMBDA is provided by the NASA Office of Space Science. This research has Jarosik, N., Bennett, C. L., Dunkley, J., et al. 2011, ApJS, 192, 14 made use of the NASA/IPAC Extragalactic Database (NED) which is operated by Jonas, J. L., Baart, E. E., & Nicolson, G. D. 1998, MNRAS, 297, 977 the Jet Propulsion Laboratory, California Institute of Technology, under contract Ka´zmierczak, M., Schmidt, M. R., Bondar, A., & Krełowski, J. 2010, MNRAS, with the National Aeronautics and Space Administration. This research makes 402, 2548 use of the SIMBAD database, operated at CDS, Strasbourg, France. King, O. G., Copley, C., Davies, R., et al. 2010, in SPIE Conf. Ser., 7741 The Planck Collaboration acknowledges the support of: ESA; CNES and Kogut, A., Banday, A. J., Bennett, C. L., et al. 1996, ApJ, 464, L5 CNRS/INSU-IN2P3-INP (France); ASI, CNR, and INAF (Italy); NASA and Kulesa, C. A., Hungerford, A. L., Walker, C. K., Zhang, X., & Lane, A. P. 2005, DoE (USA); STFC and UKSA (UK); CSIC, MICINN and JA (Spain); Tekes, ApJ, 625, 194 AoF and CSC (Finland); DLR and MPG (Germany); CSA (Canada); DTU Lagache, G. 2003, A&A, 405, 813 Space (Denmark); SER/SSO (Switzerland); RCN (Norway); SFI (Ireland); Lamarre, J., Puget, J., Ade, P. A. R., et al. 2010, A&A, 520, A9 FCT/MCTES (Portugal); and DEISA (EU). A detailed description of the Planck Leach, S. M., Cardoso, J., Baccigalupi, C., et al. 2008, A&A, 491, 597 Collaboration and a list of its members can be found at http://www.rssd. Leahy, J. P., Bersanelli, M., D’Arcangelo, O., et al. 2010, A&A, 520, A8 esa.int/index.php?project=PLANCK&page=Planck_Collaboration. Lebouteiller, V., Brandl, B., Bernard-Salas, J., Devost, D., & Houck, J. R. 2007, ApJ, 665, 390 Lebouteiller, V., Bernard-Salas, J., Whelan, D. G., et al. 2011, ApJ, 728, 45 Leitch, E. M., Readhead, A. C. S., Pearson, T. J., & Myers, S. T. 1997, ApJ, 486, References L23 Abergel, A., Bernard, J. P., Boulanger, F., et al. 2002, A&A, 389, 239 Liseau, R., White, G. J., Larsson, B., et al. 1999, A&A, 344, 342 Ali-Haïmoud, Y., Hirata, C. M., & Dickinson, C. 2009, MNRAS, 395, 1055 Lopez-Caraballo, C. H., Rubino-Martin, J. A., Rebolo, R., & Genova-Santos, R. Banday, A. J., Dickinson, C., Davies, R. D., Davis, R. J., & Górski, K. M. 2003, 2011, ApJ, 729, 25 MNRAS, 345, 897 Mandolesi, N., Bersanelli, M., Butler, R. C., et al. 2010, A&A, 520, A3 Bennett, C. L., Hill, R. S., Hinshaw, G., et al. 2003, ApJS, 148, 97 Mathis, J. S., Mezger, P. G., & Panagia, N. 1983, A&A, 128, 212 Berkhuijsen, E. M. 1972, A&AS, 5, 263 Mennella, A., Butler, R. C., Curto, A., et al. 2011, A&A, 536, A3 Berné, O., Joblin, C., Deville, Y., et al. 2007, A&A, 469, 575 Miville-Deschênes, M., Ysard, N., Lavabre, A., et al. 2008, A&A, 490, 1093 Bersanelli, M., Mandolesi, N., Butler, R. C., et al. 2010, A&A, 520, A4 Murphy, E. J., Helou, G., Condon, J. J., et al. 2010, ApJ, 709, L108 Casassus, S., Cabrera, G. F., Förster, F., et al. 2006, ApJ, 639, 951 Nordh, H. L., Fridlund, C. V. M., van Duinen, R. J., et al. 1984, A&A, 131, 221 Casassus, S., Dickinson, C., Cleary, K., et al. 2008, MNRAS, 391, 1075 Paladini, R., Burigana, C., Davies, R. D., et al. 2003, A&A, 397, 213 Castellanos, P., Casassus, S., Dickinson, C., et al. 2011, MNRAS, 411, 1137 Patel, N. A., Goldsmith, P. F., Heyer, M. H., Snell, R. L., & Pratap, P. 1998, ApJ, Cernicharo, J., Bachiller, R., & Duvert, G. 1985, A&A, 149, 273 507, 241 Compiègne, M., Abergel, A., Verstraete, L., & Habart, E. 2008, A&A, 491, 797 Perryman, M. A. C., Lindegren, L., Kovalevsky, J., et al. 1997, A&A, 323, L49 Compiègne, M., Verstraete, L., Jones, A., et al. 2011, A&A, 525, A103 Planck Collaboration 2011a, A&A, 536, A1 Condon, J. J., Cotton, W. D., Greisen, E. W., et al. 1998, AJ, 115, 1693 Planck Collaboration 2011b, A&A, 536, A2 Dame, T. M., Hartmann, D., & Thaddeus, P. 2001, ApJ, 547, 792 Planck Collaboration 2011c, A&A, 536, A7 Davies, R. D., Lasenby, A. N., Watson, R. A., et al. 1987, Nature, 326, 462 Planck Collaboration 2011d, A&A, 536, A8 Davies, R. D., Dickinson, C., Banday, A. J., et al. 2006, MNRAS, 370, 1125 Planck Collaboration 2011e, A&A, 536, A9 de Oliveira-Costa, A., Kogut, A., Devlin, M. J., et al. 1997, ApJ, 482, L17 Planck Collaboration 2011f, A&A, 536, A10 de Oliveira-Costa, A., Tegmark, M., Gutierrez, C. M., et al. 1999, ApJ, 527, L9 Planck Collaboration 2011g, A&A, 536, A11 de Oliveira-Costa, A., Tegmark, M., Davies, R. D., et al. 2004, ApJ, 606, L89 Planck Collaboration 2011h, A&A, 536, A12 de Zeeuw, P. T., Hoogerwerf, R., de Bruijne, J. H. J., Brown, A. G. A., & Blaauw, Planck Collaboration 2011i, A&A, 536, A13 A. 1999, AJ, 117, 354 Planck Collaboration 2011j, A&A, 536, A14 Delabrouille, J., & Cardoso, J. 2009, in Data Analysis in Cosmology, ed. Planck Collaboration 2011k, A&A, 536, A15 V. J. Martínez, E. Saar, E. Martínez-González, & M.-J. Pons-Bordería, Lect. Planck Collaboration 2011l, A&A, 536, A16 Notes Phys. (Berlin: Springer Verlag), 665, 159 Planck Collaboration 2011m, A&A, 536, A17 Dickinson, C., Davies, R. D., & Davis, R. J. 2003, MNRAS, 341, 369 Planck Collaboration 2011n, A&A, 536, A18 Dickinson, C., Davies, R. D., Allison, J. R., et al. 2009, ApJ, 690, 1585 Planck Collaboration 2011o, A&A, 536, A19 Dickinson, C., Casassus, S., Davies, R. D., et al. 2010, MNRAS, 407, 2223 Planck Collaboration 2011p, A&A, 536, A20 Dobler, G., & Finkbeiner, D. P. 2008, ApJ, 680, 1222 Planck Collaboration 2011q, A&A, 536, A21 Dobler, G., Draine, B., & Finkbeiner, D. P. 2009, ApJ, 699, 1374 Planck Collaboration 2011r, A&A, 536, A22 Draine, B. T., & Lazarian, A. 1998, ApJ, 508, 157 Planck Collaboration 2011s, A&A, 536, A23 Draine, B. T., & Lazarian, A. 1999, ApJ, 512, 740 Planck Collaboration 2011t, A&A, 536, A24 Encrenaz, P. J. 1974, ApJ, 189, L135 Planck Collaboration 2011u, A&A, 536, A25 Erickson, W. C. 1957, ApJ, 126, 480 Planck Collaboration 2011v, The Explanatory Supplement to the Planck Early Evans, II, N. J., Allen, L. E., Blake, G. A., et al. 2003, PASP, 115, 965 Release Compact Source Catalogue (ESA) Felli, M., Massi, F., Robberto, M., & Cesaroni, R. 2006, A&A, 453, 911 Planck Collaboration 2011w, A&A, 536, A26 Fernández-Cerezo, S., Gutiérrez, C. M., Rebolo, R., et al. 2006, MNRAS, 370, Planck HFI Core Team 2011a, A&A, 536, A4 15 Planck HFI Core Team 2011b, A&A, 536, A6 Finkbeiner, D. P. 2003, ApJS, 146, 407 Platania, P., Burigana, C., Maino, D., et al. 2003, A&A, 410, 847 Finkbeiner, D. P. 2004, ApJ, 614, 186 Reich, W. 1982, A&AS, 48, 219 Finkbeiner, D. P., Davis, M., & Schlegel, D. J. 1999, ApJ, 524, 867 Reich, P., & Reich, W. 1986, A&AS, 63, 205 Finkbeiner, D. P., Schlegel, D. J., Frank, C., & Heiles, C. 2002, ApJ, 566, 898 Reich, P., & Reich, W. 1988, A&AS, 74, 7 Finkbeiner, D. P., Langston, G. I., & Minter, A. H. 2004, ApJ, 617, 350 Reich, W., & Reich, P. 2009, in Cosmic Magnetic Fields: from Planets, to Stars Gallegos, J. E., Macías-Pérez, J. F., Gutiérrez, C. M., et al. 2001, MNRAS, 327, and Galaxies, IAU Symp., 259, 603 1178 Reich, P., Testori, J. C., & Reich, W. 2001, A&A, 376, 861 Gaustad, J. E., McCullough, P. R., Rosing, W., & Van Buren, D. 2001, PASP, Ridge, N. A., Schnee, S. L., Goodman, A. A., & COMPLETE Collaboration 113, 1326 2003, BAAS, 36, 580 Gold, B., Odegard, N., Weiland, J. L., et al. 2011, ApJS, 192, 15 Ridge, N. A., Schnee, S. L., Goodman, A. A., & Foster, J. B. 2006, ApJ, 643, Górski, K. M., Hivon, E., Banday, A. J., et al. 2005, ApJ, 622, 759 932 Habart, E., Boulanger, F., Verstraete, L., et al. 2003, A&A, 397, 623 Röllig, M., Ossenkopf, V., Jeyakumar, S., Stutzki, J., & Sternberg, A. 2006, Habart, E., Abergel, A., Walmsley, C. M., Teyssier, D., & Pety, J. 2005, A&A, A&A, 451, 917 437, 177 Rosset, C., Tristram, M., Ponthieu, N., et al. 2010, A&A, 520, A13 Haffner, L. M., Reynolds, R. J., Tufte, S. L., et al. 2003, ApJS, 149, 405 Scaife, A. M. M., Hurley-Walker, N., Green, D. A., et al. 2009, MNRAS, 400, Haslam, C. G. T., Salter, C. J., Stoffel, H., & Wilson, W. E. 1982, A&AS, 47, 1 1394 Hauser, M. G., Arendt, R. G., Kelsall, T., et al. 1998, ApJ, 508, 25 Scaife, A. M. M., Nikolic, B., Green, D. A., et al. 2010, MNRAS, 406, L45

A20, page 15 of 17 A&A 536, A20 (2011)

Sharpless, S. 1959, ApJS, 4, 257 25 Dipartimento di Fisica, Università degli Studi di Milano, via Silsbee, K., Ali-Haïmoud, Y., & Hirata, C. M. 2011, MNRAS, 411, 2750 Celoria, 16, Milano, Italy Tauber, J. A., Mandolesi, N., Puget, J., et al. 2010, A&A, 520, A1 26 Dipartimento di Fisica, Università degli Studi di Trieste, via A. Tibbs, C. T., Watson, R. A., Dickinson, C., et al. 2010, MNRAS, 402, 1969 Valerio 2, Trieste, Italy Todorovic,´ M., Davies, R. D., Dickinson, C., et al. 2010, MNRAS, 406, 1629 27 Dipartimento di Fisica, Università di Ferrara, via Saragat 1, 44122 van Dishoeck, E. F., & Black, J. H. 1989, ApJ, 340, 273 Ferrara, Italy Velusamy, T., & Langer, W. D. 2008, AJ, 136, 602 28 Watson, R. A., Rebolo, R., Rubiño-Martín, J. A., et al. 2005, ApJ, 624, L89 Dipartimento di Fisica, Università di Roma Tor Vergata, via della Williams, J. P., Bergin, E. A., Caselli, P., Myers, P. C., & Plume, R. 1998, ApJ, Ricerca Scientifica, 1, Roma, Italy 503, 689 29 Discovery Center, Niels Bohr Institute, Blegdamsvej 17, Wolfire, M. G., McKee, C. F., Hollenbach, D., & Tielens, A. G. G. M. 2003, ApJ, Copenhagen, Denmark 587, 278 30 Dpto. Astrofísica, Universidad de La Laguna (ULL), 38206 Wolfire, M. G., Tielens, A. G. G. M., Hollenbach, D., & Kaufman, M. J. 2008, La Laguna, Tenerife, Spain ApJ, 680, 384 31 European Southern Observatory, ESO Vitacura, Alonso de Cordova Wood, D. O. S., & Churchwell, E. 1989, ApJS, 69, 831 3107, Vitacura, Casilla 19001, Santiago, Chile Wootten, A., Evans, II, N. J., Snell, R., & vanden Bout, P. 1978, ApJ, 225, L143 32 European Space Agency, ESAC, Planck Science Office, Camino Wootten, A., Snell, R., & Glassgold, A. E. 1979, ApJ, 234, 876 / Young, K. E., Enoch, M. L., Evans, II, N. J., et al. 2006, ApJ, 644, 326 bajo del Castillo, s n, Urbanización Villafranca del Castillo, Ysard, N., & Verstraete, L. 2010, A&A, 509, A12 Villanueva de la Cañada, Madrid, Spain Ysard, N., Miville-Deschênes, M. A., & Verstraete, L. 2010, A&A, 509, L1 33 European Space Agency, ESTEC, Keplerlaan 1, 2201 AZ Zacchei, A., Maino, D., Baccigalupi, C., et al. 2011, A&A, 536, A5 Noordwijk, The Netherlands 34 Haverford College Astronomy Department, 370 Lancaster Avenue, 1 Aalto University Metsähovi Radio Observatory, Metsähovintie 114, Haverford, Pennsylvania, USA 35 02540 Kylmälä, Finland Helsinki Institute of Physics, Gustaf Hällströmin katu 2, University 2 Agenzia Spaziale Italiana Science Data Center, c/oESRIN,via of Helsinki, Helsinki, Finland 36 Galileo Galilei, Frascati, Italy INAF - Osservatorio Astrofisico di Catania, via S. Sofia 78, 3 Astroparticule et Cosmologie, CNRS (UMR7164), Université Catania, Italy 37 Denis Diderot Paris 7, Bâtiment Condorcet, 10 rue A. Domon INAF - Osservatorio Astronomico di Padova, vicolo et Léonie Duquet, Paris, France dell’Osservatorio 5, Padova, Italy 38 4 Astrophysics Group, Cavendish Laboratory, University of INAF - Osservatorio Astronomico di Roma, via di Frascati 33, Cambridge, J J Thomson Avenue, Cambridge CB3 0HE, UK Monte Porzio Catone, Italy 39 5 Atacama Large Millimeter/submillimeter Array, ALMA Santiago INAF - Osservatorio Astronomico di Trieste, via G.B. Tiepolo 11, Central Offices, Alonso de Cordova 3107, Vitacura, Casilla 763 Trieste, Italy 40 0355, Santiago, Chile INAF/IASF Bologna, via Gobetti 101, Bologna, Italy 41 6 CITA, University of Toronto, 60 St. George St., Toronto, ON M5S INAF/IASF Milano, via E. Bassini 15, Milano, Italy 42 3H8, Canada INRIA, Laboratoire de Recherche en Informatique, Université 7 CNRS, IRAP, 9 Av. colonel Roche, BP 44346, 31028 Toulouse Paris-Sud 11, Bâtiment 490, 91405 Orsay Cedex, France 43 Cedex 4, France Imperial College London, Astrophysics group, Blackett 8 California Institute of Technology, Pasadena, California, USA Laboratory, Prince Consort Road, London, SW7 2AZ, UK 44 9 Centre of Mathematics for Applications, University of Oslo, Infrared Processing and Analysis Center, California Institute of Blindern, Oslo, Norway Technology, Pasadena, CA 91125, USA 45 10 DAMTP, University of Cambridge, Centre for Mathematical Institut Néel, CNRS, Université Joseph Fourier Grenoble I, 25 rue Sciences, Wilberforce Road, Cambridge CB3 0WA, UK des Martyrs, Grenoble, France 46 11 DSM/Irfu/SPP, CEA-Saclay, 91191 Gif-sur-Yvette Cedex, France Institut d’Astrophysique Spatiale, CNRS (UMR8617) Université 12 DTU Space, National Space Institute, Juliane Mariesvej 30, Paris-Sud 11, Bâtiment 121, Orsay, France 47 Copenhagen, Denmark Institut d’Astrophysique de Paris, CNRS UMR7095, Université 13 Departamento de Física, Universidad de Oviedo, Avda. Calvo Pierre & Marie Curie, 98 bis boulevard Arago, Paris, France 48 Sotelo s/n, Oviedo, Spain Institute of Astronomy and Astrophysics, Academia Sinica, Taipei, 14 Department of Astronomy and Astrophysics, University of Taiwan 49 Toronto, 50 Saint George Street, Toronto, Ontario, Canada Institute of Astronomy, University of Cambridge, Madingley Road, 15 Department of Physics and Astronomy, University of British Cambridge CB3 0HA, UK 50 Columbia, 6224 Agricultural Road, Vancouver, British Columbia, Institute of Theoretical Astrophysics, University of Oslo, Blindern, Canada Oslo, Norway 51 16 Department of Physics, Gustaf Hällströmin katu 2a, University of Instituto de Astrofísica de Canarias, C/Vía Láctea s/n, La Laguna, Helsinki, Helsinki, Finland Tenerife, Spain 52 17 Department of Physics, Princeton University, Princeton, New Instituto de Física de Cantabria (CSIC-Universidad de Cantabria), Jersey, USA Avda. de los Castros s/n, Santander, Spain 53 18 Department of Physics, Purdue University, 525 Northwestern Istituto di Fisica del Plasma, CNR-ENEA-EURATOM Association, Avenue, West Lafayette, Indiana, USA via R. Cozzi 53, Milano, Italy 54 19 Department of Physics, University of California, Berkeley, Jet Propulsion Laboratory, California Institute of Technology, 4800 California, USA Oak Grove Drive, Pasadena, California, USA 55 20 Department of Physics, University of California, One Shields Jodrell Bank Centre for Astrophysics, Alan Turing Building, Avenue, Davis, California, USA School of Physics and Astronomy, The University of Manchester, 21 Department of Physics, University of California, Santa Barbara, Oxford Road, Manchester, M13 9PL, UK 56 California, USA Kavli Institute for Cosmology Cambridge, Madingley Road, 22 Department of Physics, University of Illinois at Cambridge, CB3 0HA, UK 57 Urbana-Champaign, 1110 West Green Street, Urbana, Illinois, USA LERMA, CNRS, Observatoire de Paris, 61 Avenue de 23 Dipartimento di Fisica G. Galilei, Università degli Studi di Padova, l’Observatoire, Paris, France 58 via Marzolo 8, 35131 Padova, Italy Laboratoire AIM, IRFU/Service d’Astrophysique - CEA/DSM - 24 Dipartimento di Fisica, Università La Sapienza, P. le A. Moro 2, CNRS - Université Paris Diderot, Bât. 709, CEA-Saclay, 91191 Roma, Italy Gif-sur-Yvette Cedex, France

A20, page 16 of 17 Planck Collaboration: Planck early results. XX.

59 Laboratoire Traitement et Communication de l’Information, CNRS 70 SISSA, Astrophysics Sector, via Bonomea 265, 34136, Trieste, (UMR 5141) and Télécom ParisTech, 46 rue Barrault 75634 Paris Italy Cedex 13, France 71 SUPA, Institute for Astronomy, University of Edinburgh, Royal 60 Laboratoire de Physique Subatomique et de Cosmologie, Observatory, Blackford Hill, Edinburgh EH9 3HJ, UK CNRS/IN2P3, Université Joseph Fourier Grenoble I, Institut 72 School of Physics and Astronomy, Cardiff University, Queens National Polytechnique de Grenoble, 53 rue des Martyrs, 38026 Buildings, The Parade, Cardiff, CF24 3AA, UK Grenoble Cedex, France 73 Space Sciences Laboratory, University of California, Berkeley, 61 Laboratoire de l’Accélérateur Linéaire, Université Paris-Sud 11, California, USA CNRS/IN2P3, Orsay, France 74 Spitzer Science Center, 1200 E. California Blvd., Pasadena, 62 Lawrence Berkeley National Laboratory, Berkeley, California, California, USA USA 75 Stanford University, Dept of Physics, Varian Physics Bldg, 382 via 63 Max-Planck-Institut für Astrophysik, Karl-Schwarzschild-Str. 1, Pueblo Mall, Stanford, California, USA 85741 Garching, Germany 76 Universidad de Chile, Casilla 36-D, Santiago, Chile 64 Max-Planck-Institut für Radioastronomie, Auf dem Hügel 69, 77 Université de Toulouse, UPS-OMP, IRAP, 31028 Toulouse Cedex 53121 Bonn, Germany 4, France 65 MilliLab, VTT Technical Research Centre of Finland, Tietotie 3, 78 Universities Space Research Association, Stratospheric Espoo, Finland Observatory for Infrared Astronomy, MS 211-3, Moffett Field, CA 66 National University of Ireland, Department of Experimental 94035, USA Physics, Maynooth, Co. Kildare, Ireland 79 University of Granada, Departamento de Física Teórica y del 67 Niels Bohr Institute, Blegdamsvej 17, Copenhagen, Denmark Cosmos, Facultad de Ciencias, Granada, Spain 68 Observational Cosmology, Mail Stop 367-17, California Institute 80 University of Miami, Knight Physics Building, 1320 Campo Sano of Technology, Pasadena, CA, 91125, USA Dr., Coral Gables, Florida, USA 69 Optical Science Laboratory, University College London, Gower 81 Warsaw University Observatory, Aleje Ujazdowskie 4, 00-478 Street, London, UK Warszawa, Poland

A20, page 17 of 17